53 h1_gen_ZMass_ =
new TH1F(
"gen_ZMass",
"Generated Z mass", 200, 0., 150.);
57 h1_gen_ZEta_ =
new TH1F(
"gen_ZEta",
"Eta of gen Z", 200, -6., 6.);
61 h1_gen_ZPhi_ =
new TH1F(
"gen_ZPhi",
"Phi of gen Z", 200, -4., 4.);
69 h1_gen_ZPt_ =
new TH1F(
"gen_ZPt",
"Pt of gen Z", 200, 0., 100.);
77 h1_reco_ZEta_ =
new TH1F(
"reco_ZEta",
"Eta of reco Z", 200, -6., 6.);
81 h1_reco_ZTheta_ =
new TH1F(
"reco_ZTheta",
"Theta of reco Z", 200, 0., 4.);
89 h1_reco_ZPhi_ =
new TH1F(
"reco_ZPhi",
"Phi of reco Z", 100, -4., 4.);
93 h1_reco_ZPt_ =
new TH1F(
"reco_ZPt",
"Pt of reco Z", 200, 0., 100.);
135 for (HepMC::GenEvent::particle_const_iterator
p = myGenEvent->particles_begin();
p != myGenEvent->particles_end();
138 if ((*p)->pdg_id() == 23 && (*p)->status() == 2) {
143 0.5 *
log(((*p)->momentum().e() + (*p)->momentum().pz()) / ((*p)->momentum().e() - (*p)->momentum().pz()));
163 h1_mcElePt_ =
new TH1F(
"mcElePt",
"p_{T} of MC electrons", 300, 0., 300.);
167 h1_mcEleEta_ =
new TH1F(
"mcEleEta",
"Eta of MC electrons", 100, -4., 4.);
171 h1_mcElePhi_ =
new TH1F(
"mcElePhi",
"Phi of MC electrons", 100, -4., 4.);
181 for (HepMC::GenEvent::particle_const_iterator
p = myGenEvent->particles_begin();
p != myGenEvent->particles_end();
183 if (
abs((*p)->pdg_id()) == 11) {
198 h1_nEleReco_ =
new TH1F(
"h1_nEleReco",
"h1_nEleReco", 20, 0, 20);
201 h1_recoEleEnergy_ =
new TH1F(
"recoEleEnergy",
"EleEnergy from SC", 300, 0., 300.);
205 h1_recoElePt_ =
new TH1F(
"recoElePt",
"p_{T} of reco electrons", 300, 0., 300.);
209 h1_recoEleEta_ =
new TH1F(
"recoEleEta",
"Eta of reco electrons", 100, -4., 4.);
213 h1_recoElePhi_ =
new TH1F(
"recoElePhi",
"Phi of reco electrons", 100, -4., 4.);
243 std::cout <<
"Start with ZeePlots::writeEleHistograms(), done file_->cd(); " << std::endl;
250 std::cout <<
"Done with ZeePlots::writeEleHistograms() " << std::endl;
258 std::cout <<
"Start with ZeePlots::writeMCEleHistograms(), done file_->cd(); " << std::endl;
265 std::cout <<
"Done with ZeePlots::writeMCEleHistograms() " << std::endl;
280 new TH1F(
"h1_HLT2ElectronRelaxed_FiredEvents",
"h1_HLT2ElectronRelaxed_FiredEvents", 5, 0, 5);
283 new TH1F(
"h1_HLT1Electron_HLT2Electron_FiredEvents",
"h1_HLT1Electron_HLT2Electron_FiredEvents", 5, 0, 5);
285 "h1_HLT1Electron_HLT2ElectronRelaxed_FiredEvents",
"h1_HLT1Electron_HLT2ElectronRelaxed_FiredEvents", 5, 0, 5);
287 "h1_HLT2Electron_HLT2ElectronRelaxed_FiredEvents",
"h1_HLT2Electron_HLT2ElectronRelaxed_FiredEvents", 5, 0, 5);
290 new TH1F(
"h1_HLT1Electron_HLT2Electron_HLT2ElectronRelaxed_FiredEvents",
291 "h1_HLT1Electron_HLT2Electron_HLT2ElectronRelaxed_FiredEvents",
302 int hltCount = hltTriggerResultHandle->
size();
304 bool aHLTResults[200] = {
false};
306 for (
int i = 0;
i < hltCount;
i++) {
307 aHLTResults[
i] = hltTriggerResultHandle->
accept(
i);
318 if (aHLTResults[32] && !aHLTResults[34] && !aHLTResults[35])
321 if (aHLTResults[34] && !aHLTResults[32] && !aHLTResults[35])
324 if (aHLTResults[35] && !aHLTResults[32] && !aHLTResults[34])
327 if (aHLTResults[32] && aHLTResults[34] && !aHLTResults[35])
330 if (aHLTResults[32] && aHLTResults[35] && !aHLTResults[34])
333 if (aHLTResults[34] && aHLTResults[35] && !aHLTResults[32])
336 if (aHLTResults[32] && aHLTResults[34] && aHLTResults[35])
345 if (myClass == 0 || myClass == 100)
348 std::cout <<
"[ZeePlots::fillEleClassesPlots]Done gold" << std::endl;
350 if (myClass == 40 || myClass == 140)
353 std::cout <<
"[ZeePlots::fillEleClassesPlots]Done crack" << std::endl;
355 if ((myClass >= 30 && myClass <= 34) || (myClass >= 130 && myClass <= 134))
358 std::cout <<
"[ZeePlots::fillEleClassesPlots]Done shower" << std::endl;
360 if (myClass == 10 || myClass == 20 || myClass == 110 || myClass == 120)
363 std::cout <<
"[ZeePlots::fillEleClassesPlots]Done" << std::endl;