CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
ElectronMcMiniAODSignalValidator.cc
Go to the documentation of this file.
1 // user include files
2 #include "Validation/RecoEgamma/plugins/ElectronMcMiniAODSignalValidator.h"
3 #include "CLHEP/Units/GlobalPhysicalConstants.h"
4 #include "FWCore/MessageLogger/interface/MessageLogger.h"
5 
6 // user include files
7 
8 using namespace reco;
9 using namespace pat;
10 
11 ElectronMcSignalValidatorMiniAOD::ElectronMcSignalValidatorMiniAOD(const edm::ParameterSet& iConfig)
12  : ElectronDqmAnalyzerBase(iConfig) {
13  mcTruthCollection_ = consumes<edm::View<reco::GenParticle> >(
14  iConfig.getParameter<edm::InputTag>("mcTruthCollection")); // prunedGenParticles
15  electronToken_ =
16  consumes<pat::ElectronCollection>(iConfig.getParameter<edm::InputTag>("electrons")); // slimmedElectrons
17  electronTokenEndcaps_ =
18  consumes<pat::ElectronCollection>(iConfig.getParameter<edm::InputTag>("electrons_endcaps")); // slimmedElectrons
19 
20  edm::ParameterSet histosSet = iConfig.getParameter<edm::ParameterSet>("histosCfg");
21  edm::ParameterSet isolationSet = iConfig.getParameter<edm::ParameterSet>("isolationCfg");
22 
23  maxPt_ = iConfig.getParameter<double>("MaxPt");
24  maxAbsEta_ = iConfig.getParameter<double>("MaxAbsEta");
25  deltaR_ = iConfig.getParameter<double>("DeltaR");
26  deltaR2_ = deltaR_ * deltaR_;
27  matchingIDs_ = iConfig.getParameter<std::vector<int> >("MatchingID");
28  matchingMotherIDs_ = iConfig.getParameter<std::vector<int> >("MatchingMotherID");
29  outputInternalPath_ = iConfig.getParameter<std::string>("OutputFolderName");
30 
31  // histos bining and limits
32 
33  xyz_nbin = histosSet.getParameter<int>("Nbinxyz");
34 
35  pt_nbin = histosSet.getParameter<int>("Nbinpt");
36  pt2D_nbin = histosSet.getParameter<int>("Nbinpt2D");
37  pteff_nbin = histosSet.getParameter<int>("Nbinpteff");
38  pt_max = histosSet.getParameter<double>("Ptmax");
39 
40  fhits_nbin = histosSet.getParameter<int>("Nbinfhits");
41  fhits_max = histosSet.getParameter<double>("Fhitsmax");
42 
43  eta_nbin = histosSet.getParameter<int>("Nbineta");
44  eta2D_nbin = histosSet.getParameter<int>("Nbineta2D");
45  eta_min = histosSet.getParameter<double>("Etamin");
46  eta_max = histosSet.getParameter<double>("Etamax");
47 
48  detamatch_nbin = histosSet.getParameter<int>("Nbindetamatch");
49  detamatch2D_nbin = histosSet.getParameter<int>("Nbindetamatch2D");
50  detamatch_min = histosSet.getParameter<double>("Detamatchmin");
51  detamatch_max = histosSet.getParameter<double>("Detamatchmax");
52 
53  dphi_nbin = histosSet.getParameter<int>("Nbindphi");
54  dphi_min = histosSet.getParameter<double>("Dphimin");
55  dphi_max = histosSet.getParameter<double>("Dphimax");
56 
57  dphimatch_nbin = histosSet.getParameter<int>("Nbindphimatch");
58  dphimatch2D_nbin = histosSet.getParameter<int>("Nbindphimatch2D");
59  dphimatch_min = histosSet.getParameter<double>("Dphimatchmin");
60  dphimatch_max = histosSet.getParameter<double>("Dphimatchmax");
61 
62  hoe_nbin = histosSet.getParameter<int>("Nbinhoe");
63  hoe_min = histosSet.getParameter<double>("Hoemin");
64  hoe_max = histosSet.getParameter<double>("Hoemax");
65 
66  mee_nbin = histosSet.getParameter<int>("Nbinmee");
67  mee_min = histosSet.getParameter<double>("Meemin");
68  mee_max = histosSet.getParameter<double>("Meemax");
69 
70  poptrue_nbin = histosSet.getParameter<int>("Nbinpoptrue");
71  poptrue_min = histosSet.getParameter<double>("Poptruemin");
72  poptrue_max = histosSet.getParameter<double>("Poptruemax");
73 
74  set_EfficiencyFlag = histosSet.getParameter<bool>("EfficiencyFlag");
75  set_StatOverflowFlag = histosSet.getParameter<bool>("StatOverflowFlag");
76 
77  ele_nbin = histosSet.getParameter<int>("NbinELE");
78  ele_min = histosSet.getParameter<double>("ELE_min");
79  ele_max = histosSet.getParameter<double>("ELE_max");
80 
81  // so to please coverity...
82 
83  h1_recEleNum = nullptr;
84 
85  h1_ele_vertexPt = nullptr;
86  h1_ele_vertexEta = nullptr;
87  h1_ele_vertexPt_nocut = nullptr;
88 
89  h1_scl_SigIEtaIEta_mAOD = nullptr;
90  h1_scl_SigIEtaIEta_mAOD_barrel = nullptr;
91  h1_scl_SigIEtaIEta_mAOD_endcaps = nullptr;
92 
93  h2_ele_foundHitsVsEta = nullptr;
94  h2_ele_foundHitsVsEta_mAOD = nullptr;
95 
96  h2_ele_PoPtrueVsEta = nullptr;
97  h2_ele_sigmaIetaIetaVsPt = nullptr;
98 
99  h1_ele_HoE_mAOD = nullptr;
100  h1_ele_HoE_mAOD_barrel = nullptr;
101  h1_ele_HoE_mAOD_endcaps = nullptr;
102  h1_ele_mee_all = nullptr;
103  h1_ele_mee_os = nullptr;
104 
105  h1_ele_fbrem_mAOD = nullptr;
106  h1_ele_fbrem_mAOD_barrel = nullptr;
107  h1_ele_fbrem_mAOD_endcaps = nullptr;
108 
109  h1_ele_dEtaSc_propVtx_mAOD = nullptr;
110  h1_ele_dEtaSc_propVtx_mAOD_barrel = nullptr;
111  h1_ele_dEtaSc_propVtx_mAOD_endcaps = nullptr;
112  h1_ele_dPhiCl_propOut_mAOD = nullptr;
113  h1_ele_dPhiCl_propOut_mAOD_barrel = nullptr;
114  h1_ele_dPhiCl_propOut_mAOD_endcaps = nullptr;
115 
116  h1_ele_chargedHadronRelativeIso_mAOD = nullptr;
117  h1_ele_chargedHadronRelativeIso_mAOD_barrel = nullptr;
118  h1_ele_chargedHadronRelativeIso_mAOD_endcaps = nullptr;
119  h1_ele_neutralHadronRelativeIso_mAOD = nullptr;
120  h1_ele_neutralHadronRelativeIso_mAOD_barrel = nullptr;
121  h1_ele_neutralHadronRelativeIso_mAOD_endcaps = nullptr;
122  h1_ele_photonRelativeIso_mAOD = nullptr;
123  h1_ele_photonRelativeIso_mAOD_barrel = nullptr;
124  h1_ele_photonRelativeIso_mAOD_endcaps = nullptr;
125 }
126 
127 ElectronMcSignalValidatorMiniAOD::~ElectronMcSignalValidatorMiniAOD() {}
128 
129 void ElectronMcSignalValidatorMiniAOD::bookHistograms(DQMStore::IBooker& iBooker,
130  edm::Run const&,
131  edm::EventSetup const&) {
132  iBooker.setCurrentFolder(outputInternalPath_);
133 
134  setBookIndex(-1);
135  setBookPrefix("h");
136  setBookEfficiencyFlag(set_EfficiencyFlag);
137  setBookStatOverflowFlag(set_StatOverflowFlag);
138 
139  // rec event collections sizes
140  h1_recEleNum = bookH1(iBooker, "recEleNum", "# rec electrons", ele_nbin, ele_min, ele_max, "N_{ele}");
141  // matched electrons
142  setBookPrefix("h_mc");
143  setBookPrefix("h_ele");
144  h1_ele_vertexPt =
145  bookH1withSumw2(iBooker, "vertexPt", "ele transverse momentum", pt_nbin, 0., pt_max, "p_{T vertex} (GeV/c)");
146  h1_ele_vertexEta = bookH1withSumw2(iBooker, "vertexEta", "ele momentum eta", eta_nbin, eta_min, eta_max, "#eta");
147  h1_ele_vertexPt_nocut =
148  bookH1withSumw2(iBooker, "vertexPt_nocut", "pT of prunned electrons", pt_nbin, 0., 20., "p_{T vertex} (GeV/c)");
149  h2_ele_PoPtrueVsEta = bookH2withSumw2(iBooker,
150  "PoPtrueVsEta",
151  "ele momentum / gen momentum vs eta",
152  eta2D_nbin,
153  eta_min,
154  eta_max,
155  50,
156  poptrue_min,
157  poptrue_max);
158  h2_ele_sigmaIetaIetaVsPt =
159  bookH2(iBooker, "sigmaIetaIetaVsPt", "SigmaIetaIeta vs pt", 100, 0., pt_max, 100, 0., 0.05);
160 
161  // all electrons
162  setBookPrefix("h_ele");
163  h1_ele_mee_all = bookH1withSumw2(iBooker,
164  "mee_all",
165  "ele pairs invariant mass, all reco electrons",
166  mee_nbin,
167  mee_min,
168  mee_max,
169  "m_{ee} (GeV/c^{2})",
170  "Events",
171  "ELE_LOGY E1 P");
172  h1_ele_mee_os = bookH1withSumw2(iBooker,
173  "mee_os",
174  "ele pairs invariant mass, opp. sign",
175  mee_nbin,
176  mee_min,
177  mee_max,
178  "m_{e^{+}e^{-}} (GeV/c^{2})",
179  "Events",
180  "ELE_LOGY E1 P");
181 
182  // matched electron, superclusters
183  setBookPrefix("h_scl");
184  h1_scl_SigIEtaIEta_mAOD = bookH1withSumw2(iBooker,
185  "SigIEtaIEta_mAOD",
186  "ele supercluster sigma ieta ieta",
187  100,
188  0.,
189  0.05,
190  "#sigma_{i#eta i#eta}",
191  "Events",
192  "ELE_LOGY E1 P");
193  h1_scl_SigIEtaIEta_mAOD_barrel = bookH1withSumw2(iBooker,
194  "SigIEtaIEta_mAOD_barrel",
195  "ele supercluster sigma ieta ieta, barrel",
196  100,
197  0.,
198  0.05,
199  "#sigma_{i#eta i#eta}",
200  "Events",
201  "ELE_LOGY E1 P");
202  h1_scl_SigIEtaIEta_mAOD_endcaps = bookH1withSumw2(iBooker,
203  "SigIEtaIEta_mAOD_endcaps",
204  "ele supercluster sigma ieta ieta, endcaps",
205  100,
206  0.,
207  0.05,
208  "#sigma_{i#eta i#eta}",
209  "Events",
210  "ELE_LOGY E1 P");
211 
212  // matched electron, gsf tracks
213  setBookPrefix("h_ele");
214  h2_ele_foundHitsVsEta = bookH2(iBooker,
215  "foundHitsVsEta",
216  "ele track # found hits vs eta",
217  eta2D_nbin,
218  eta_min,
219  eta_max,
220  fhits_nbin,
221  0.,
222  fhits_max);
223  h2_ele_foundHitsVsEta_mAOD = bookH2(iBooker,
224  "foundHitsVsEta_mAOD",
225  "ele track # found hits vs eta",
226  eta2D_nbin,
227  eta_min,
228  eta_max,
229  fhits_nbin,
230  0.,
231  fhits_max);
232 
233  // matched electrons, matching
234  setBookPrefix("h_ele");
235  h1_ele_HoE_mAOD = bookH1withSumw2(iBooker,
236  "HoE_mAOD",
237  "ele hadronic energy / em energy",
238  hoe_nbin,
239  hoe_min,
240  hoe_max,
241  "H/E",
242  "Events",
243  "ELE_LOGY E1 P");
244  h1_ele_HoE_mAOD_barrel = bookH1withSumw2(iBooker,
245  "HoE_mAOD_barrel",
246  "ele hadronic energy / em energy, barrel",
247  hoe_nbin,
248  hoe_min,
249  hoe_max,
250  "H/E",
251  "Events",
252  "ELE_LOGY E1 P");
253  h1_ele_HoE_mAOD_endcaps = bookH1withSumw2(iBooker,
254  "HoE_mAOD_endcaps",
255  "ele hadronic energy / em energy, endcaps",
256  hoe_nbin,
257  hoe_min,
258  hoe_max,
259  "H/E",
260  "Events",
261  "ELE_LOGY E1 P");
262  h1_ele_dEtaSc_propVtx_mAOD = bookH1withSumw2(iBooker,
263  "dEtaSc_propVtx_mAOD",
264  "ele #eta_{sc} - #eta_{tr}, prop from vertex",
265  detamatch_nbin,
266  detamatch_min,
267  detamatch_max,
268  "#eta_{sc} - #eta_{tr}",
269  "Events",
270  "ELE_LOGY E1 P");
271  h1_ele_dEtaSc_propVtx_mAOD_barrel = bookH1withSumw2(iBooker,
272  "dEtaSc_propVtx_mAOD_barrel",
273  "ele #eta_{sc} - #eta_{tr}, prop from vertex, barrel",
274  detamatch_nbin,
275  detamatch_min,
276  detamatch_max,
277  "#eta_{sc} - #eta_{tr}",
278  "Events",
279  "ELE_LOGY E1 P");
280  h1_ele_dEtaSc_propVtx_mAOD_endcaps = bookH1withSumw2(iBooker,
281  "dEtaSc_propVtx_mAOD_endcaps",
282  "ele #eta_{sc} - #eta_{tr}, prop from vertex, endcaps",
283  detamatch_nbin,
284  detamatch_min,
285  detamatch_max,
286  "#eta_{sc} - #eta_{tr}",
287  "Events",
288  "ELE_LOGY E1 P");
289  h1_ele_dPhiCl_propOut_mAOD = bookH1withSumw2(iBooker,
290  "dPhiCl_propOut_mAOD",
291  "ele #phi_{cl} - #phi_{tr}, prop from outermost",
292  dphimatch_nbin,
293  dphimatch_min,
294  dphimatch_max,
295  "#phi_{seedcl} - #phi_{tr} (rad)",
296  "Events",
297  "ELE_LOGY E1 P");
298  h1_ele_dPhiCl_propOut_mAOD_barrel = bookH1withSumw2(iBooker,
299  "dPhiCl_propOut_mAOD_barrel",
300  "ele #phi_{cl} - #phi_{tr}, prop from outermost, barrel",
301  dphimatch_nbin,
302  dphimatch_min,
303  dphimatch_max,
304  "#phi_{seedcl} - #phi_{tr} (rad)",
305  "Events",
306  "ELE_LOGY E1 P");
307  h1_ele_dPhiCl_propOut_mAOD_endcaps = bookH1withSumw2(iBooker,
308  "dPhiCl_propOut_mAOD_endcaps",
309  "ele #phi_{cl} - #phi_{tr}, prop from outermost, endcaps",
310  dphimatch_nbin,
311  dphimatch_min,
312  dphimatch_max,
313  "#phi_{seedcl} - #phi_{tr} (rad)",
314  "Events",
315  "ELE_LOGY E1 P");
316 
317  // fbrem
318  h1_ele_fbrem_mAOD = bookH1withSumw2(
319  iBooker, "fbrem_mAOD", "ele brem fraction, mode of GSF components", 100, 0., 1., "P_{in} - P_{out} / P_{in}");
320  h1_ele_fbrem_mAOD_barrel = bookH1withSumw2(iBooker,
321  "fbrem_mAOD_barrel",
322  "ele brem fraction for barrel, mode of GSF components",
323  100,
324  0.,
325  1.,
326  "P_{in} - P_{out} / P_{in}");
327  h1_ele_fbrem_mAOD_endcaps = bookH1withSumw2(iBooker,
328  "fbrem_mAOD_endcaps",
329  "ele brem franction for endcaps, mode of GSF components",
330  100,
331  0.,
332  1.,
333  "P_{in} - P_{out} / P_{in}");
334 
335  // -- pflow over pT
336  h1_ele_chargedHadronRelativeIso_mAOD = bookH1withSumw2(iBooker,
337  "chargedHadronRelativeIso_mAOD",
338  "chargedHadronRelativeIso",
339  100,
340  0.0,
341  2.,
342  "chargedHadronRelativeIso",
343  "Events",
344  "ELE_LOGY E1 P");
345  h1_ele_chargedHadronRelativeIso_mAOD_barrel = bookH1withSumw2(iBooker,
346  "chargedHadronRelativeIso_mAOD_barrel",
347  "chargedHadronRelativeIso for barrel",
348  100,
349  0.0,
350  2.,
351  "chargedHadronRelativeIso_barrel",
352  "Events",
353  "ELE_LOGY E1 P");
354  h1_ele_chargedHadronRelativeIso_mAOD_endcaps = bookH1withSumw2(iBooker,
355  "chargedHadronRelativeIso_mAOD_endcaps",
356  "chargedHadronRelativeIso for endcaps",
357  100,
358  0.0,
359  2.,
360  "chargedHadronRelativeIso_endcaps",
361  "Events",
362  "ELE_LOGY E1 P");
363  h1_ele_neutralHadronRelativeIso_mAOD = bookH1withSumw2(iBooker,
364  "neutralHadronRelativeIso_mAOD",
365  "neutralHadronRelativeIso",
366  100,
367  0.0,
368  2.,
369  "neutralHadronRelativeIso",
370  "Events",
371  "ELE_LOGY E1 P");
372  h1_ele_neutralHadronRelativeIso_mAOD_barrel = bookH1withSumw2(iBooker,
373  "neutralHadronRelativeIso_mAOD_barrel",
374  "neutralHadronRelativeIso for barrel",
375  100,
376  0.0,
377  2.,
378  "neutralHadronRelativeIso_barrel",
379  "Events",
380  "ELE_LOGY E1 P");
381  h1_ele_neutralHadronRelativeIso_mAOD_endcaps = bookH1withSumw2(iBooker,
382  "neutralHadronRelativeIso_mAOD_endcaps",
383  "neutralHadronRelativeIso for endcaps",
384  100,
385  0.0,
386  2.,
387  "neutralHadronRelativeIso_endcaps",
388  "Events",
389  "ELE_LOGY E1 P");
390  h1_ele_photonRelativeIso_mAOD = bookH1withSumw2(iBooker,
391  "photonRelativeIso_mAOD",
392  "photonRelativeIso",
393  100,
394  0.0,
395  2.,
396  "photonRelativeIso",
397  "Events",
398  "ELE_LOGY E1 P");
399  h1_ele_photonRelativeIso_mAOD_barrel = bookH1withSumw2(iBooker,
400  "photonRelativeIso_mAOD_barrel",
401  "photonRelativeIso for barrel",
402  100,
403  0.0,
404  2.,
405  "photonRelativeIso_barrel",
406  "Events",
407  "ELE_LOGY E1 P");
408  h1_ele_photonRelativeIso_mAOD_endcaps = bookH1withSumw2(iBooker,
409  "photonRelativeIso_mAOD_endcaps",
410  "photonRelativeIso for endcaps",
411  100,
412  0.0,
413  2.,
414  "photonRelativeIso_endcaps",
415  "Events",
416  "ELE_LOGY E1 P");
417 }
418 
419 void ElectronMcSignalValidatorMiniAOD::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
420  // get collections
421  auto electrons = iEvent.getHandle(electronToken_);
422  auto electrons_endcaps = iEvent.getHandle(electronTokenEndcaps_);
423  auto genParticles = iEvent.getHandle(mcTruthCollection_);
424 
425  edm::Handle<pat::ElectronCollection> mergedElectrons;
426 
427  edm::LogInfo("ElectronMcSignalValidatorMiniAOD::analyze")
428  << "Treating event " << iEvent.id() << " with " << electrons.product()->size() << " electrons";
429  edm::LogInfo("ElectronMcSignalValidatorMiniAOD::analyze")
430  << "Treating event " << iEvent.id() << " with " << electrons_endcaps.product()->size()
431  << " multi slimmed electrons";
432 
433  h1_recEleNum->Fill((*electrons).size());
434 
435  //===============================================
436  // all rec electrons
437  //===============================================
438 
439  pat::Electron gsfElectron;
440  pat::ElectronCollection::const_iterator el1;
441  std::vector<pat::Electron>::const_iterator el3;
442  std::vector<pat::Electron>::const_iterator el4;
443 
444  //===============================================
445  // get a vector with EB & EE
446  //===============================================
447  std::vector<pat::Electron> localCollection;
448  int iBarrels = 0;
449  int iEndcaps = 0;
450 
451  // looking for EB
452  for (el1 = electrons->begin(); el1 != electrons->end(); el1++) {
453  if (el1->isEB()) {
454  localCollection.push_back(*el1);
455  iBarrels += 1;
456  }
457  }
458 
459  // looking for EE
460  for (el1 = electrons_endcaps->begin(); el1 != electrons_endcaps->end(); el1++) {
461  if (el1->isEE()) {
462  localCollection.push_back(*el1);
463  iEndcaps += 1;
464  }
465  }
466 
467  for (el3 = localCollection.begin(); el3 != localCollection.end(); el3++) {
468  for (el4 = el3 + 1; el4 != localCollection.end(); el4++) {
469  math::XYZTLorentzVector p12 = el3->p4() + el4->p4();
470  float mee2 = p12.Dot(p12);
471  h1_ele_mee_all->Fill(sqrt(mee2));
472  if (el3->charge() * el4->charge() < 0.) {
473  h1_ele_mee_os->Fill(sqrt(mee2));
474  }
475  }
476  }
477 
478  //===============================================
479  // charge mis-ID
480  //===============================================
481 
482  int mcNum = 0, gamNum = 0, eleNum = 0;
483  bool matchingMotherID;
484 
485  //===============================================
486  // association mc-reco
487  //===============================================
488 
489  for (size_t i = 0; i < genParticles->size(); i++) {
490  // number of mc particles
491  mcNum++;
492 
493  // counts photons
494  if ((*genParticles)[i].pdgId() == 22) {
495  gamNum++;
496  }
497 
498  // select requested mother matching gen particle
499  // always include single particle with no mother
500  const Candidate* mother = (*genParticles)[i].mother(0);
501  matchingMotherID = false;
502  for (unsigned int ii = 0; ii < matchingMotherIDs_.size(); ii++) {
503  if (mother == nullptr) {
504  matchingMotherID = true;
505  } else if (mother->pdgId() == matchingMotherIDs_[ii]) {
506  if (mother->numberOfDaughters() <= 2) {
507  matchingMotherID = true;
508  }
509  } // end of mother if test
510 
511  } // end of for loop
512  if (!matchingMotherID) {
513  continue;
514  }
515 
516  // electron preselection
517  if ((*genParticles)[i].pt() > maxPt_ || std::abs((*genParticles)[i].eta()) > maxAbsEta_) {
518  continue;
519  }
520  eleNum++;
521 
522  // find best matched electron
523  bool okGsfFound = false;
524  bool passMiniAODSelection = true;
525  double gsfOkRatio = 999999.;
526  bool isEBflag = false;
527  bool isEEflag = false;
528  pat::Electron bestGsfElectron;
529 
530  for (el3 = localCollection.begin(); el3 != localCollection.end(); el3++) {
531  double dphi = el3->phi() - (*genParticles)[i].phi();
532  if (std::abs(dphi) > CLHEP::pi) {
533  dphi = dphi < 0 ? (CLHEP::twopi) + dphi : dphi - CLHEP::twopi;
534  }
535  double deltaR2 = (el3->eta() - (*genParticles)[i].eta()) * (el3->eta() - (*genParticles)[i].eta()) + dphi * dphi;
536  if (deltaR2 < deltaR2_) {
537  if ((((*genParticles)[i].pdgId() == 11) && (el3->charge() < 0.)) ||
538  (((*genParticles)[i].pdgId() == -11) && (el3->charge() > 0.))) {
539  double tmpGsfRatio = el3->p() / (*genParticles)[i].p();
540  if (std::abs(tmpGsfRatio - 1) < std::abs(gsfOkRatio - 1)) {
541  gsfOkRatio = tmpGsfRatio;
542  bestGsfElectron = *el3;
543  okGsfFound = true;
544  }
545  }
546  }
547  } // end *electrons loop
548 
549  if (okGsfFound) {
550  //------------------------------------
551  // analysis when the mc track is found
552  //------------------------------------
553  passMiniAODSelection = bestGsfElectron.pt() >= 5.;
554  double one_over_pt = 1. / bestGsfElectron.pt();
555  isEBflag = bestGsfElectron.isEB();
556  isEEflag = bestGsfElectron.isEE();
557 
558  // electron related distributions
559  h1_ele_vertexPt->Fill(bestGsfElectron.pt());
560  h1_ele_vertexEta->Fill(bestGsfElectron.eta());
561  if ((bestGsfElectron.scSigmaIEtaIEta() == 0.) && (bestGsfElectron.fbrem() == 0.))
562  h1_ele_vertexPt_nocut->Fill(bestGsfElectron.pt());
563  // track related distributions
564  h2_ele_foundHitsVsEta->Fill(bestGsfElectron.eta(), bestGsfElectron.gsfTrack()->numberOfValidHits());
565 
566  // generated distributions for matched electrons
567  h2_ele_PoPtrueVsEta->Fill(bestGsfElectron.eta(), bestGsfElectron.p() / (*genParticles)[i].p());
568 
569  if (passMiniAODSelection) { // Pt > 5.
570  h2_ele_sigmaIetaIetaVsPt->Fill(bestGsfElectron.pt(), bestGsfElectron.scSigmaIEtaIEta());
571 
572  // supercluster related distributions
573  h1_scl_SigIEtaIEta_mAOD->Fill(bestGsfElectron.scSigmaIEtaIEta());
574  h1_ele_dEtaSc_propVtx_mAOD->Fill(bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
575  h1_ele_dPhiCl_propOut_mAOD->Fill(bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
576 
577  // track related distributions
578  h2_ele_foundHitsVsEta_mAOD->Fill(bestGsfElectron.eta(), bestGsfElectron.gsfTrack()->numberOfValidHits());
579 
580  // match distributions
581  h1_ele_HoE_mAOD->Fill(bestGsfElectron.hcalOverEcal());
582 
583  // fbrem
584  h1_ele_fbrem_mAOD->Fill(bestGsfElectron.fbrem());
585 
586  // -- pflow over pT
587 
588  h1_ele_chargedHadronRelativeIso_mAOD->Fill(bestGsfElectron.pfIsolationVariables().sumChargedHadronPt *
589  one_over_pt);
590  h1_ele_neutralHadronRelativeIso_mAOD->Fill(bestGsfElectron.pfIsolationVariables().sumNeutralHadronEt *
591  one_over_pt);
592  h1_ele_photonRelativeIso_mAOD->Fill(bestGsfElectron.pfIsolationVariables().sumPhotonEt * one_over_pt);
593 
594  if (isEBflag) {
595  // supercluster related distributions
596  h1_scl_SigIEtaIEta_mAOD_barrel->Fill(bestGsfElectron.scSigmaIEtaIEta());
597  h1_ele_dEtaSc_propVtx_mAOD_barrel->Fill(bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
598  h1_ele_dPhiCl_propOut_mAOD_barrel->Fill(bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
599  // match distributions
600  h1_ele_HoE_mAOD_barrel->Fill(bestGsfElectron.hcalOverEcal());
601  // fbrem
602  h1_ele_fbrem_mAOD_barrel->Fill(bestGsfElectron.fbrem());
603 
604  h1_ele_chargedHadronRelativeIso_mAOD_barrel->Fill(bestGsfElectron.pfIsolationVariables().sumChargedHadronPt *
605  one_over_pt);
606  h1_ele_neutralHadronRelativeIso_mAOD_barrel->Fill(bestGsfElectron.pfIsolationVariables().sumNeutralHadronEt *
607  one_over_pt);
608  h1_ele_photonRelativeIso_mAOD_barrel->Fill(bestGsfElectron.pfIsolationVariables().sumPhotonEt * one_over_pt);
609  }
610 
611  // supercluster related distributions
612  else if (isEEflag) {
613  h1_scl_SigIEtaIEta_mAOD_endcaps->Fill(bestGsfElectron.scSigmaIEtaIEta());
614  h1_ele_dEtaSc_propVtx_mAOD_endcaps->Fill(bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
615  h1_ele_dPhiCl_propOut_mAOD_endcaps->Fill(bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
616  // match distributions
617  h1_ele_HoE_mAOD_endcaps->Fill(bestGsfElectron.hcalOverEcal());
618  // fbrem
619  h1_ele_fbrem_mAOD_endcaps->Fill(bestGsfElectron.fbrem());
620  h1_ele_chargedHadronRelativeIso_mAOD_endcaps->Fill(bestGsfElectron.pfIsolationVariables().sumChargedHadronPt *
621  one_over_pt);
622  h1_ele_neutralHadronRelativeIso_mAOD_endcaps->Fill(bestGsfElectron.pfIsolationVariables().sumNeutralHadronEt *
623  one_over_pt);
624  h1_ele_photonRelativeIso_mAOD_endcaps->Fill(bestGsfElectron.pfIsolationVariables().sumPhotonEt * one_over_pt);
625  }
626  }
627  }
628  //} // end loop i_elec
629  } // end loop size_t i
630 }
reco::GsfElectron::pfIsolationVariables
const PflowIsolationVariables & pfIsolationVariables() const
Definition: GsfElectron.h:729
ElectronMcSignalValidatorMiniAOD_cfi.electrons_endcaps
tuple electrons_endcaps
Definition: ElectronMcSignalValidatorMiniAOD_cfi.py:60
mps_fire.i
i
Definition: mps_fire.py:428
reco::LeafCandidate::pt
double pt() const final
transverse momentum
Definition: LeafCandidate.h:146
reco::GsfElectron::p4
const LorentzVector & p4(P4Kind kind) const
Definition: GsfElectron.cc:217
reco::Candidate::mother
virtual const Candidate * mother(size_type i=0) const =0
return pointer to mother
dqm::implementation::NavigatorBase::setCurrentFolder
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:32
pat::Electron::gsfTrack
reco::GsfTrackRef gsfTrack() const override
override the reco::GsfElectron::gsfTrack method, to access the internal storage of the supercluster ...
reco::GsfElectron::fbrem
float fbrem() const
Definition: GsfElectron.h:809
genParticleCandidates2GenParticles_cfi.genParticles
tuple genParticles
Definition: genParticleCandidates2GenParticles_cfi.py:3
dt_dqm_sourceclient_common_cff.reco
tuple reco
Definition: dt_dqm_sourceclient_common_cff.py:111
cuy.ii
int ii
Definition: cuy.py:589
ElectronMcSignalValidatorMiniAOD::ElectronMcSignalValidatorMiniAOD
ElectronMcSignalValidatorMiniAOD(const edm::ParameterSet &)
Definition: ElectronMcMiniAODSignalValidator.cc:11
HI_PhotonSkim_cff.electrons
tuple electrons
Definition: HI_PhotonSkim_cff.py:77
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
reco::GsfElectron::isEE
bool isEE() const
Definition: GsfElectron.h:329
reco::GsfElectron::isEB
bool isEB() const
Definition: GsfElectron.h:328
pi
const Double_t pi
Definition: trackSplitPlot.h:36
dqm::impl::MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:292
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
ElectronMcSignalValidatorMiniAOD::mcTruthCollection_
edm::EDGetTokenT< edm::View< reco::GenParticle > > mcTruthCollection_
Definition: ElectronMcMiniAODSignalValidator.h:27
ElectronDqmAnalyzerBase::bookH1
MonitorElement * bookH1(DQMStore::IBooker &, const std::string &name, const std::string &title, int nchX, double lowX, double highX, const std::string &titleX="", const std::string &titleY="Events", Option_t *option="E1 P")
Definition: ElectronDqmAnalyzerBase.cc:57
edm::Event::getHandle
Handle< PROD > getHandle(EDGetTokenT< PROD > token) const
Definition: Event.h:563
reco::Candidate::numberOfDaughters
virtual size_type numberOfDaughters() const =0
number of daughters
reco::GsfElectron::deltaEtaSuperClusterTrackAtVtx
float deltaEtaSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:225
iEvent
int iEvent
Definition: GenABIO.cc:224
reco::LeafCandidate::p
double p() const final
magnitude of momentum vector
Definition: LeafCandidate.h:123
reco::GsfElectron::deltaPhiSeedClusterTrackAtCalo
float deltaPhiSeedClusterTrackAtCalo() const
Definition: GsfElectron.h:229
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
reco::GsfElectron::PflowIsolationVariables::sumPhotonEt
float sumPhotonEt
sum pt of PF photons // old float photonIso ;
Definition: GsfElectron.h:665
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
ElectronDqmAnalyzerBase::setBookPrefix
void setBookPrefix(const std::string &)
Definition: ElectronDqmAnalyzerBase.cc:30
reco::GsfElectron::PflowIsolationVariables::sumNeutralHadronEt
float sumNeutralHadronEt
sum pt of neutral hadrons // old float neutralHadronIso ;
Definition: GsfElectron.h:664
reco::GsfElectron::scSigmaIEtaIEta
float scSigmaIEtaIEta() const
Definition: GsfElectron.h:496
ElectronMcSignalValidatorMiniAOD::electronTokenEndcaps_
edm::EDGetTokenT< pat::ElectronCollection > electronTokenEndcaps_
Definition: ElectronMcMiniAODSignalValidator.h:29
ElectronDqmAnalyzerBase::setBookEfficiencyFlag
void setBookEfficiencyFlag(const bool &)
Definition: ElectronDqmAnalyzerBase.cc:34
reco::Candidate::pdgId
virtual int pdgId() const =0
PDG identifier.
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
reco::deltaR2
constexpr auto deltaR2(const T1 &t1, const T2 &t2) -> decltype(t1.eta())
Definition: deltaR.h:16
ElectronDqmAnalyzerBase::bookH2withSumw2
MonitorElement * bookH2withSumw2(DQMStore::IBooker &, const std::string &name, const std::string &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, const std::string &titleX="", const std::string &titleY="", Option_t *option="COLZ")
Definition: ElectronDqmAnalyzerBase.cc:140
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
pat::Electron
Analysis-level electron class.
Definition: Electron.h:51
ElectronMcSignalValidatorMiniAOD::analyze
void analyze(const edm::Event &, const edm::EventSetup &) override
Definition: ElectronMcMiniAODSignalValidator.cc:419
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:59
ElectronDqmAnalyzerBase::bookH2
MonitorElement * bookH2(DQMStore::IBooker &, const std::string &name, const std::string &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, const std::string &titleX="", const std::string &titleY="", Option_t *option="COLZ")
Definition: ElectronDqmAnalyzerBase.cc:111
ElectronDqmAnalyzerBase::bookH1withSumw2
MonitorElement * bookH1withSumw2(DQMStore::IBooker &, const std::string &name, const std::string &title, int nchX, double lowX, double highX, const std::string &titleX="", const std::string &titleY="Events", Option_t *option="E1 P")
Definition: ElectronDqmAnalyzerBase.cc:83
reco::GsfElectron::hcalOverEcal
float hcalOverEcal(const ShowerShape &ss, int depth) const
Definition: GsfElectron.h:425
ElectronMcSignalValidatorMiniAOD::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: ElectronMcMiniAODSignalValidator.cc:129
ElectronMcSignalValidatorMiniAOD::electronToken_
edm::EDGetTokenT< pat::ElectronCollection > electronToken_
Definition: ElectronMcMiniAODSignalValidator.h:28
reco::LeafCandidate::phi
double phi() const final
momentum azimuthal angle
Definition: LeafCandidate.h:148
ElectronDqmAnalyzerBase::setBookStatOverflowFlag
void setBookStatOverflowFlag(const bool &)
Definition: ElectronDqmAnalyzerBase.cc:36
reco::GsfElectron::PflowIsolationVariables::sumChargedHadronPt
float sumChargedHadronPt
sum-pt of charged Hadron // old float chargedHadronIso ;
Definition: GsfElectron.h:663
edm::Run
Definition: Run.h:45
reco::LeafCandidate::eta
double eta() const final
momentum pseudorapidity
Definition: LeafCandidate.h:152
Generated for CMSSW Reference Manual by doxygen 1.8.5