14 japp->RootWriteLock();
23 dHist_Pi0InvariantMass = GetOrCreate_Histogram<TH1I>(
"InvariantMass_Pi0",
";#gamma#gamma Invariant Mass (GeV/c^{2})", 600, 0.0, 0.3);
33 vector<const DNeutralParticle*> locUnusedNeutralParticles;
36 bool locPassesCutFlag =
true;
37 vector<double> locInvariantMasses;
38 for(
size_t loc_i = 0; loc_i < locUnusedNeutralParticles.size(); ++loc_i)
41 for(
size_t loc_j = loc_i + 1; loc_j < locUnusedNeutralParticles.size(); ++loc_j)
46 double locInvariantMass = locP4.M();
48 locPassesCutFlag =
false;
51 set<const DNeutralParticleHypothesis*> locPhotonSet;
52 locPhotonSet.insert(locPhoton1);
53 locPhotonSet.insert(locPhoton2);
59 locInvariantMasses.push_back(locInvariantMass);
68 for(
size_t loc_i = 0; loc_i < locInvariantMasses.size(); ++loc_i)
73 return locPassesCutFlag;
TDirectoryFile * ChangeTo_BaseDirectory(void)
TDirectoryFile * CreateAndChangeTo_ActionDirectory(void)
TLorentzVector DLorentzVector
DLorentzVector lorentzMomentum(void) const
void Get_UnusedNeutralParticles(JEventLoop *locEventLoop, const DParticleCombo *locParticleCombo, vector< const DNeutralParticle * > &locUnusedNeutralParticles) const