Brief Kinematic Fitting Example
From GlueXWiki
This is a test
. z
vector<DKinematicData> initial; vector<DKinematicData> final; vector<DKinematicData> initial_post; vector<DKinematicData> final_post; //Add the beam particle to the initial state DKinematicData beam_photon; DVector3 mom(0.0, 0.0, 9.0); beam_photon.setMomentum(mom); beam_photon.setPosition(pos); beam_photon.setMass(0.0); beam_photon.setCharge(0.0); beam_photon.setMassFixed(); beam_photon.clearErrorMatrix(); initial.push_back(beam_photon); //Add the target to the initial state DKinematicData kd_targ; kd_targ.setMass(0.93827); kd_targ.setCharge(1.0); kd_targ.setMassFixed(); kd_targ.setMomentum( DVector3(0.0, 0.0, 0.0) ); kd_targ.clearErrorMatrix(); initial.push_back(kd_targ); //Add all our particles to the final state final.push_back(proton[k]); final.push_back(pip[l]); final.push_back(pip[m]); final.push_back(pim[n]); final.push_back(pim[o]); double cl=kfit.Prob();