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