14 japp->RootWriteLock();
20 dEgamma = GetOrCreate_Histogram<TH1I>(
"Egamma",
"TAGGER photon energy; E_{#gamma}", 400, 0., 12.);
22 dMM2_M3pi = GetOrCreate_Histogram<TH2I>(
"MM2_M3pi",
"MM^{2} off #pi^{+}#pi^{-}#pi^{0} vs M_{#pi^{+}#pi^{-}#pi^{0}}; M_{#pi^{+}#pi^{-}#pi^{0}}; MM^{2}", 200, 0.0, 2.0, 200, -1., 1.);
23 dProton_dEdx_P = GetOrCreate_Histogram<TH2I>(
"Proton_dEdx_P",
"dE/dx vs p; p; dE/dx",200,0,2,500,0,5);
24 dProton_P_Theta = GetOrCreate_Histogram<TH2I>(
"Proton_P_Theta",
"p vs #theta; #theta; p (GeV/c)",180,0,180,120,0,12);
26 dDeltaE_M3pi = GetOrCreate_Histogram<TH2I>(
"dDeltaE_M3pi",
"#Delta E vs M_{#pi^{+}#pi^{-}#pi^{0}}; M_{#pi^{+}#pi^{-}#pi^{0}}; #Delta E (tagger - tracks)", 200, 0.0, 2.0, 200, -5., 5.);
27 dPhi3pi_PhiP = GetOrCreate_Histogram<TH2I>(
"Phi3pi_PhiP",
"#phi_{#pi^{+}#pi^{-}#pi^{0}} vs. #phi_{p}; #phi_{p}; #phi_{#pi^{+}#pi^{-}#pi^{0}}", 360, -180.0, 180.0, 360, -180., 180.);
28 dDeltaPhi_M3pi = GetOrCreate_Histogram<TH2I>(
"DeltaPhi_M3pi",
"#Delta#phi vs M_{#pi^{+}#pi^{-}#pi^{0}}; M_{#pi^{+}#pi^{-}#pi^{0}}; #Delta#phi", 200, 0.0, 2.0, 360, 0.0, 360.0);
30 dEgamma_M3pi_ProtonTag = GetOrCreate_Histogram<TH2I>(
"dEgamma_M3pi_ProtonTag",
"E_{#gamma} vs M_{#pi^{+}#pi^{-}#pi^{0}}; M_{#pi^{+}#pi^{-}#pi^{0}}; E_{#gamma}", 200, 0.0, 2.0, 240, 0., 6.);
31 dMppizero_M2pi = GetOrCreate_Histogram<TH2I>(
"Mppizero_M2pi",
"Psuedo-Dalitz p#pi^{+}#pi^{-}#pi^{0}; M_{#pi^{+}#pi^{-}}; M_{p#pi^{0}}", 200, 0.0, 1.0, 200, 1.0, 4.0);
32 dMppiplus_M2pi = GetOrCreate_Histogram<TH2I>(
"Mppiplus_M2pi",
"Psuedo-Dalitz p#pi^{+}#pi^{-}#pi^{0}; M_{#pi^{0}#pi^{-}}; M_{p#pi^{+}}", 200, 0.0, 1.0, 200, 1.0, 4.0);
33 dMppiminus_M2pi = GetOrCreate_Histogram<TH2I>(
"Mppiminus_M2pi",
"Psuedo-Dalitz p#pi^{+}#pi^{-}#pi^{0}; M_{#pi^{0}#pi^{+}}; M_{p#pi^{-}}", 200, 0.0, 1.0, 200, 1.0, 4.0);
35 dMM2_M3pi_CoplanarTag = GetOrCreate_Histogram<TH2I>(
"MM2_M3pi_CoplanarTag",
"MM^{2} off p #pi^{+}#pi^{-}#pi^{0} vs M_{#pi^{+}#pi^{-}#pi^{0}}: Coplanar Tag; M_{#pi^{+}#pi^{-}#pi^{0}}; MM^{2}", 200, 0.0, 2.0, 200, -1., 1.);
36 dDeltaE_M3pi_CoplanarTag = GetOrCreate_Histogram<TH2I>(
"dDeltaE_M3pi_CoplanarTag",
"#Delta E vs M_{#pi^{+}#pi^{-}#pi^{0}}: Coplanar Tag; M_{#pi^{+}#pi^{-}#pi^{0}}; #Delta E (tagger - p#pi^{+}#pi^{-}#pi^{0})", 200, 0.0, 2.0, 200, -5., 5.);
37 dMM2_DeltaE_CoplanarTag = GetOrCreate_Histogram<TH2I>(
"dMM2_DeltaE_CoplanarTag",
"MM^{2} vs #Delta E: Coplanar Tag; #Delta E (tagger - p#pi^{+}#pi^{-}#pi^{0}); MM^{2}", 200, -5., 5., 200, -1., 1.);
39 dMM2_M3pi_ProtonTag = GetOrCreate_Histogram<TH2I>(
"MM2_M3pi_ProtonTag",
"MM^{2} off p #pi^{+}#pi^{-}#pi^{0} vs M_{#pi^{+}#pi^{-}#pi^{0}}: Proton Tag; M_{#pi^{+}#pi^{-}#pi^{0}}; MM^{2}", 200, 0.0, 2.0, 200, -1., 1.);
40 dDeltaE_M3pi_ProtonTag = GetOrCreate_Histogram<TH2I>(
"dDeltaE_M3pi_ProtonTag",
"#Delta E vs M_{#pi^{+}#pi^{-}#pi^{0}}: Proton Tag; M_{#pi^{+}#pi^{-}#pi^{0}}; #Delta E (tagger - p#pi^{+}#pi^{-}#pi^{0})", 200, 0.0, 2.0, 200, -5., 5.);
41 dMM2_DeltaE_ProtonTag = GetOrCreate_Histogram<TH2I>(
"dMM2_DeltaE_ProtonTag",
"MM^{2} vs #Delta E: Proton Tag; #Delta E (tagger - p#pi^{+}#pi^{-}#pi^{0}); MM^{2}", 200, -5., 5., 200, -1., 1.);
56 double locBeamPhotonEnergy = locBeamPhoton->energy();
59 locSumInitP4.SetXYZM(0, 0, 0, 0.938);
60 locSumInitP4 += locBeamPhoton->lorentzMomentum();
87 dEgamma->Fill(locBeamPhotonEnergy);
89 dMM2_M3pi->Fill(locOmegaP4.M(), locMissingP4.M2());
92 double locDeltaPhi = (locProtonP4.Phi() - locOmegaP4.Phi())*180./TMath::Pi();
93 if(locDeltaPhi > 360.) locDeltaPhi -= 360.;
94 if(locDeltaPhi < 0.) locDeltaPhi += 360.;
95 dPhi3pi_PhiP->Fill(locProtonP4.Phi()*180./TMath::Pi(), locOmegaP4.Phi()*180./TMath::Pi());
99 if(locDeltaPhi < 175. || locDeltaPhi > 185.)
106 if(locOmegaP4.M() > 0.7 && locOmegaP4.M() < 0.9)
113 if(locOmegaP4.M() > 0.7 && locOmegaP4.M() < 0.9)
117 if(fabs(locMissingP4.M2()) < 0.05 && fabs(locMissingP4.E()) < 0.5) {
120 if(locBeamPhotonEnergy > 2.5 && locBeamPhotonEnergy < 3.0){
121 dMppizero_M2pi->Fill((locPiPlusP4+locPiMinusP4).M(), (locProtonP4+locPi0P4).M());
122 dMppiplus_M2pi->Fill((locPi0P4+locPiMinusP4).M(), (locProtonP4+locPiPlusP4).M());
123 dMppiminus_M2pi->Fill((locPiPlusP4+locPi0P4).M(), (locProtonP4+locPiMinusP4).M());
128 if(fabs(locMissingP4.M2()) < 0.05 && fabs(locMissingP4.E()) < 0.5) {
131 dProton_P_Theta->Fill(locProtonP4.Vect().Theta()*180/TMath::Pi(), locProtonP4.Vect().Mag());
vector< const DKinematicData * > Get_FinalParticles_Measured(void) const
TDirectoryFile * ChangeTo_BaseDirectory(void)
const DKinematicData * Get_FinalParticle_Measured(size_t locFinalParticleIndex) const
const DReaction * Get_Reaction(void) const
TH2I * dDeltaE_M3pi_ProtonTag
void Initialize(JEventLoop *locEventLoop)
const DTrackTimeBased * Get_TrackTimeBased(void) const
TDirectoryFile * CreateAndChangeTo_ActionDirectory(void)
TH2I * dEgamma_M3pi_ProtonTag
bool Perform_Action(JEventLoop *locEventLoop, const DParticleCombo *locParticleCombo)
const DKinematicData * Get_InitialParticle_Measured(void) const
TLorentzVector DLorentzVector
const DChargedTrackHypothesis * Get_Hypothesis(Particle_t locPID) const
bool Get_UseKinFitResultsFlag(void) const
TH2I * dDeltaE_M3pi_CoplanarTag
const DAnalysisUtilities * dAnalysisUtilities
DLorentzVector Calc_FinalStateP4(const DReaction *locReaction, const DParticleCombo *locParticleCombo, size_t locStepIndex, bool locUseKinFitDataFlag) const
TH2I * dMM2_M3pi_CoplanarTag
DLorentzVector Calc_MissingP4(const DReaction *locReaction, const DParticleCombo *locParticleCombo, bool locUseKinFitDataFlag) const
TH2I * dMM2_DeltaE_CoplanarTag
vector< const DKinematicData * > Get_FinalParticles(void) const
DLorentzVector lorentzMomentum(void) const
const JObject * Get_FinalParticle_SourceObject(size_t locFinalParticleIndex) const
const DKinematicData * Get_InitialParticle(void) const
TH2I * dMM2_DeltaE_ProtonTag
TH2I * dMM2_M3pi_ProtonTag
const DParticleComboStep * Get_ParticleComboStep(size_t locStepIndex) const