Brief Kinematic Fitting Example
From GlueXWiki
Revision as of 14:05, 15 March 2012 by Wilevine (Talk | contribs) (Created page with "This is a {{{test}}}. z <pre> vector<DKinematicData> initial; vector<DKinematicData> final; vector<DKinematicData> initial_post; vector<DKinematicData> final_post; //Add the b...")
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();