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