32 locMissingP4 += locProtonP4Init;
43 for(
size_t loc_i = 0; loc_i < 5; ++loc_i) {
44 locMissingP4 -= locParticles[loc_i]->lorentzMomentum();
46 if(locParticles[loc_i]->PID() ==
Gamma)
47 loc2g_P4 += locParticles[loc_i]->lorentzMomentum();
48 else if(locParticles[loc_i]->PID() ==
PiPlus)
49 locPiPlusP4 += locParticles[loc_i]->lorentzMomentum();
50 else if(locParticles[loc_i]->PID() ==
PiMinus)
51 locPiMinusP4 += locParticles[loc_i]->lorentzMomentum();
53 locProtonP4 += locParticles[loc_i]->lorentzMomentum();
56 locCandidateP4 = loc2g_P4 + locPiPlusP4 + locPiMinusP4;
58 double locDeltaPhi = (locProtonP4.Phi() - locCandidateP4.Phi())*180./TMath::Pi();
59 if(locDeltaPhi > 360.) locDeltaPhi -= 360.;
60 if(locDeltaPhi < 0.) locDeltaPhi += 360.;
63 if(locDeltaPhi < 160. || locDeltaPhi > 200.)
66 if(fabs(locMissingP4.M2()) > 0.1 || fabs(locMissingP4.E()) > 1.0)
bool Perform_Action(JEventLoop *locEventLoop, const DParticleCombo *locParticleCombo)
const DReaction * Get_Reaction(void) const
const DKinematicData * Get_InitialParticle_Measured(void) const
TLorentzVector DLorentzVector
bool Get_UseKinFitResultsFlag(void) const
vector< const DKinematicData * > Get_FinalParticles(const DReaction *locReaction, bool locIncludeMissingFlag=true, bool locIncludeDecayingFlag=true, Charge_t locCharge=d_AllCharges) const
vector< const DKinematicData * > Get_FinalParticles_Measured(const DReaction *locReaction, Charge_t locCharge=d_AllCharges) const
DLorentzVector lorentzMomentum(void) const
const DKinematicData * Get_InitialParticle(void) const
const DParticleComboStep * Get_ParticleComboStep(size_t locStepIndex) const
void Initialize(JEventLoop *locEventLoop)