CMS 3D CMS Logo
ElectronMcFakeValidator.cc
Go to the documentation of this file.
1 
2 // user include files
3 #include "Validation/RecoEgamma/plugins/ElectronMcFakeValidator.h"
4 
5 #include "RecoEgamma/EgammaElectronAlgos/interface/ElectronUtilities.h"
6 
7 #include "DQMServices/Core/interface/DQMStore.h"
8 
9 #include "DataFormats/GsfTrackReco/interface/GsfTrack.h"
10 #include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
11 #include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"
12 #include "DataFormats/EgammaReco/interface/BasicClusterFwd.h"
13 #include "DataFormats/EgammaReco/interface/SuperClusterFwd.h"
14 #include "DataFormats/EgammaReco/interface/ElectronSeed.h"
15 #include "DataFormats/EgammaReco/interface/ElectronSeedFwd.h"
16 #include "DataFormats/EcalDetId/interface/EcalSubdetector.h"
17 #include "DataFormats/BeamSpot/interface/BeamSpot.h"
18 #include "DataFormats/JetReco/interface/GenJetCollection.h"
19 
20 #include "DataFormats/Common/interface/Handle.h"
21 #include "DataFormats/Common/interface/ValueMap.h"
22 
23 #include "DataFormats/VertexReco/interface/Vertex.h"
24 #include "DataFormats/VertexReco/interface/VertexFwd.h"
25 
26 #include "FWCore/ParameterSet/interface/ParameterSet.h"
27 #include "FWCore/Framework/interface/EDAnalyzer.h"
28 #include "FWCore/Framework/interface/Event.h"
29 #include "FWCore/Framework/interface/MakerMacros.h"
30 #include "FWCore/MessageLogger/interface/MessageLogger.h"
31 
32 #include "CLHEP/Units/GlobalPhysicalConstants.h"
33 #include "TMath.h"
34 #include "TFile.h"
35 #include "TH1F.h"
36 #include "TH1I.h"
37 #include "TH2F.h"
38 #include "TProfile.h"
39 #include "TTree.h"
40 #include <vector>
41 #include <iostream>
42 
43 using namespace reco;
44 
45 ElectronMcFakeValidator::ElectronMcFakeValidator(const edm::ParameterSet &conf) : ElectronDqmAnalyzerBase(conf) {
46  electronCollection_ = consumes<reco::GsfElectronCollection>(conf.getParameter<edm::InputTag>("electronCollection"));
47  electronCoreCollection_ =
48  consumes<reco::GsfElectronCoreCollection>(conf.getParameter<edm::InputTag>("electronCoreCollection"));
49  electronTrackCollection_ =
50  consumes<reco::GsfTrackCollection>(conf.getParameter<edm::InputTag>("electronTrackCollection"));
51  electronSeedCollection_ =
52  consumes<reco::ElectronSeedCollection>(conf.getParameter<edm::InputTag>("electronSeedCollection"));
53  matchingObjectCollection_ =
54  consumes<reco::GenJetCollection>(conf.getParameter<edm::InputTag>("matchingObjectCollection"));
55  offlineVerticesCollection_ =
56  consumes<reco::VertexCollection>(conf.getParameter<edm::InputTag>("offlinePrimaryVertices"));
57 
58  beamSpotTag_ = consumes<reco::BeamSpot>(conf.getParameter<edm::InputTag>("beamSpot"));
59  readAOD_ = conf.getParameter<bool>("readAOD");
60 
61  isoFromDepsTk03Tag_ = consumes<edm::ValueMap<double>>(conf.getParameter<edm::InputTag>("isoFromDepsTk03"));
62  isoFromDepsTk04Tag_ = consumes<edm::ValueMap<double>>(conf.getParameter<edm::InputTag>("isoFromDepsTk04"));
63  isoFromDepsEcalFull03Tag_ =
64  consumes<edm::ValueMap<double>>(conf.getParameter<edm::InputTag>("isoFromDepsEcalFull03"));
65  isoFromDepsEcalFull04Tag_ =
66  consumes<edm::ValueMap<double>>(conf.getParameter<edm::InputTag>("isoFromDepsEcalFull04"));
67  isoFromDepsEcalReduced03Tag_ =
68  consumes<edm::ValueMap<double>>(conf.getParameter<edm::InputTag>("isoFromDepsEcalReduced03"));
69  isoFromDepsEcalReduced04Tag_ =
70  consumes<edm::ValueMap<double>>(conf.getParameter<edm::InputTag>("isoFromDepsEcalReduced04"));
71  isoFromDepsHcal03Tag_ = consumes<edm::ValueMap<double>>(conf.getParameter<edm::InputTag>("isoFromDepsHcal03"));
72  isoFromDepsHcal04Tag_ = consumes<edm::ValueMap<double>>(conf.getParameter<edm::InputTag>("isoFromDepsHcal04"));
73 
74  maxPt_ = conf.getParameter<double>("MaxPt");
75  maxAbsEta_ = conf.getParameter<double>("MaxAbsEta");
76  deltaR_ = conf.getParameter<double>("DeltaR");
77  inputFile_ = conf.getParameter<std::string>("InputFile");
78  outputFile_ = conf.getParameter<std::string>("OutputFile");
79  inputInternalPath_ = conf.getParameter<std::string>("InputFolderName");
80  outputInternalPath_ = conf.getParameter<std::string>("OutputFolderName");
81 
82  // histos bining and limits
83 
84  edm::ParameterSet histosSet = conf.getParameter<edm::ParameterSet>("histosCfg");
85 
86  xyz_nbin = histosSet.getParameter<int>("Nbinxyz");
87 
88  p_nbin = histosSet.getParameter<int>("Nbinp");
89  p2D_nbin = histosSet.getParameter<int>("Nbinp2D");
90  p_max = histosSet.getParameter<double>("Pmax");
91 
92  pt_nbin = histosSet.getParameter<int>("Nbinpt");
93  pt2D_nbin = histosSet.getParameter<int>("Nbinpt2D");
94  pteff_nbin = histosSet.getParameter<int>("Nbinpteff");
95  pt_max = histosSet.getParameter<double>("Ptmax");
96 
97  fhits_nbin = histosSet.getParameter<int>("Nbinfhits");
98  fhits_max = histosSet.getParameter<double>("Fhitsmax");
99 
100  lhits_nbin = histosSet.getParameter<int>("Nbinlhits");
101  lhits_max = histosSet.getParameter<double>("Lhitsmax");
102 
103  eop_nbin = histosSet.getParameter<int>("Nbineop");
104  eop2D_nbin = histosSet.getParameter<int>("Nbineop2D");
105  eop_max = histosSet.getParameter<double>("Eopmax");
106  eopmaxsht = histosSet.getParameter<double>("Eopmaxsht");
107 
108  eta_nbin = histosSet.getParameter<int>("Nbineta");
109  eta2D_nbin = histosSet.getParameter<int>("Nbineta2D");
110  eta_min = histosSet.getParameter<double>("Etamin");
111  eta_max = histosSet.getParameter<double>("Etamax");
112 
113  deta_nbin = histosSet.getParameter<int>("Nbindeta");
114  deta_min = histosSet.getParameter<double>("Detamin");
115  deta_max = histosSet.getParameter<double>("Detamax");
116 
117  detamatch_nbin = histosSet.getParameter<int>("Nbindetamatch");
118  detamatch2D_nbin = histosSet.getParameter<int>("Nbindetamatch2D");
119  detamatch_min = histosSet.getParameter<double>("Detamatchmin");
120  detamatch_max = histosSet.getParameter<double>("Detamatchmax");
121 
122  phi_nbin = histosSet.getParameter<int>("Nbinphi");
123  phi2D_nbin = histosSet.getParameter<int>("Nbinphi2D");
124  phi_min = histosSet.getParameter<double>("Phimin");
125  phi_max = histosSet.getParameter<double>("Phimax");
126 
127  dphi_nbin = histosSet.getParameter<int>("Nbindphi");
128  dphi_min = histosSet.getParameter<double>("Dphimin");
129  dphi_max = histosSet.getParameter<double>("Dphimax");
130 
131  dphimatch_nbin = histosSet.getParameter<int>("Nbindphimatch");
132  dphimatch2D_nbin = histosSet.getParameter<int>("Nbindphimatch2D");
133  dphimatch_min = histosSet.getParameter<double>("Dphimatchmin");
134  dphimatch_max = histosSet.getParameter<double>("Dphimatchmax");
135 
136  mee_nbin = histosSet.getParameter<int>("Nbinmee");
137  mee_min = histosSet.getParameter<double>("Meemin");
138  mee_max = histosSet.getParameter<double>("Meemax");
139 
140  hoe_nbin = histosSet.getParameter<int>("Nbinhoe");
141  hoe_min = histosSet.getParameter<double>("Hoemin");
142  hoe_max = histosSet.getParameter<double>("Hoemax");
143 
144  popmatching_nbin = histosSet.getParameter<int>("Nbinpopmatching");
145  popmatching_min = histosSet.getParameter<double>("Popmatchingmin");
146  popmatching_max = histosSet.getParameter<double>("Popmatchingmax");
147 
148  set_EfficiencyFlag = histosSet.getParameter<bool>("EfficiencyFlag");
149  set_StatOverflowFlag = histosSet.getParameter<bool>("StatOverflowFlag");
150 
151  opv_nbin = histosSet.getParameter<int>("NbinOPV");
152  opv_min = histosSet.getParameter<double>("OPV_min");
153  opv_max = histosSet.getParameter<double>("OPV_max");
154 
155  ele_nbin = histosSet.getParameter<int>("NbinELE");
156  ele_min = histosSet.getParameter<double>("ELE_min");
157  ele_max = histosSet.getParameter<double>("ELE_max");
158 
159  core_nbin = histosSet.getParameter<int>("NbinCORE");
160  core_min = histosSet.getParameter<double>("CORE_min");
161  core_max = histosSet.getParameter<double>("CORE_max");
162 
163  track_nbin = histosSet.getParameter<int>("NbinTRACK");
164  track_min = histosSet.getParameter<double>("TRACK_min");
165  track_max = histosSet.getParameter<double>("TRACK_max");
166 
167  seed_nbin = histosSet.getParameter<int>("NbinSEED");
168  seed_min = histosSet.getParameter<double>("SEED_min");
169  seed_max = histosSet.getParameter<double>("SEED_max");
170 
171  // so to please coverity
172  h1_matchingObjectNum = nullptr;
173  h1_recEleNum_ = nullptr;
174  h1_recCoreNum_ = nullptr;
175  h1_recTrackNum_ = nullptr;
176  h1_recSeedNum_ = nullptr;
177  h1_recOfflineVertices_ = nullptr;
178 
179  h1_matchingObjectEta = nullptr;
180  h1_matchingObjectAbsEta = nullptr;
181  h1_matchingObjectP = nullptr;
182  h1_matchingObjectPt = nullptr;
183  h1_matchingObjectPhi = nullptr;
184  h1_matchingObjectZ = nullptr;
185 
186  h1_ele_EoverP_all = nullptr;
187  h1_ele_EseedOP_all = nullptr;
188  h1_ele_EoPout_all = nullptr;
189  h1_ele_EeleOPout_all = nullptr;
190  h1_ele_dEtaSc_propVtx_all = nullptr;
191  h1_ele_dPhiSc_propVtx_all = nullptr;
192  h1_ele_dEtaCl_propOut_all = nullptr;
193  h1_ele_dPhiCl_propOut_all = nullptr;
194  h1_ele_TIP_all = nullptr;
195  h1_ele_HoE_all = nullptr;
196  h1_ele_vertexEta_all = nullptr;
197  h1_ele_vertexPt_all = nullptr;
198  h1_ele_mee_all = nullptr;
199  h1_ele_mee_os = nullptr;
200 
201  h2_ele_E2mnE1vsMee_all = nullptr;
202  h2_ele_E2mnE1vsMee_egeg_all = nullptr;
203 
204  h1_ele_matchingObjectEta_matched = nullptr;
205  h1_ele_matchingObjectAbsEta_matched = nullptr;
206  h1_ele_matchingObjectPt_matched = nullptr;
207  h1_ele_matchingObjectPhi_matched = nullptr;
208  h1_ele_matchingObjectZ_matched = nullptr;
209 
210  h1_ele_charge = nullptr;
211  h2_ele_chargeVsEta = nullptr;
212  h2_ele_chargeVsPhi = nullptr;
213  h2_ele_chargeVsPt = nullptr;
214  h1_ele_vertexP = nullptr;
215  h1_ele_vertexPt = nullptr;
216  h2_ele_vertexPtVsEta = nullptr;
217  h2_ele_vertexPtVsPhi = nullptr;
218  h1_ele_vertexEta = nullptr;
219  h2_ele_vertexEtaVsPhi = nullptr;
220  h1_ele_vertexAbsEta = nullptr;
221  h1_ele_vertexPhi = nullptr;
222  h1_ele_vertexX = nullptr;
223  h1_ele_vertexY = nullptr;
224  h1_ele_vertexZ = nullptr;
225  h1_ele_vertexTIP = nullptr;
226  h2_ele_vertexTIPVsEta = nullptr;
227  h2_ele_vertexTIPVsPhi = nullptr;
228  h2_ele_vertexTIPVsPt = nullptr;
229 
230  h1_ele_PoPmatchingObject = nullptr;
231  h2_ele_PoPmatchingObjectVsEta = nullptr;
232  h2_ele_PoPmatchingObjectVsPhi = nullptr;
233  h2_ele_PoPmatchingObjectVsPt = nullptr;
234  h1_ele_PoPmatchingObject_barrel = nullptr;
235  h1_ele_PoPmatchingObject_endcaps = nullptr;
236 
237  h1_ele_EtaMnEtamatchingObject = nullptr;
238  h2_ele_EtaMnEtamatchingObjectVsEta = nullptr;
239  h2_ele_EtaMnEtamatchingObjectVsPhi = nullptr;
240  h2_ele_EtaMnEtamatchingObjectVsPt = nullptr;
241  h1_ele_PhiMnPhimatchingObject = nullptr;
242  h1_ele_PhiMnPhimatchingObject2 = nullptr;
243  h2_ele_PhiMnPhimatchingObjectVsEta = nullptr;
244  h2_ele_PhiMnPhimatchingObjectVsPhi = nullptr;
245  h2_ele_PhiMnPhimatchingObjectVsPt = nullptr;
246 
247  h1_scl_En_ = nullptr;
248  h1_scl_EoEmatchingObject_barrel = nullptr;
249  h1_scl_EoEmatchingObject_endcaps = nullptr;
250  h1_scl_Et_ = nullptr;
251  h2_scl_EtVsEta_ = nullptr;
252  h2_scl_EtVsPhi_ = nullptr;
253  h2_scl_EtaVsPhi_ = nullptr;
254  h1_scl_Eta_ = nullptr;
255  h1_scl_Phi_ = nullptr;
256 
257  h1_scl_SigIEtaIEta_ = nullptr;
258  h1_scl_SigIEtaIEta_barrel_ = nullptr;
259  h1_scl_SigIEtaIEta_endcaps_ = nullptr;
260  h1_scl_full5x5_sigmaIetaIeta_ = nullptr; // new 2014.01.12
261  h1_scl_full5x5_sigmaIetaIeta_barrel_ = nullptr; // new 2014.01.12
262  h1_scl_full5x5_sigmaIetaIeta_endcaps_ = nullptr; // new 2014.01.12
263  h1_scl_E1x5_ = nullptr;
264  h1_scl_E1x5_barrel_ = nullptr;
265  h1_scl_E1x5_endcaps_ = nullptr;
266  h1_scl_E2x5max_ = nullptr;
267  h1_scl_E2x5max_barrel_ = nullptr;
268  h1_scl_E2x5max_endcaps_ = nullptr;
269  h1_scl_E5x5_ = nullptr;
270  h1_scl_E5x5_barrel_ = nullptr;
271  h1_scl_E5x5_endcaps_ = nullptr;
272 
273  h1_ele_ambiguousTracks = nullptr;
274  h2_ele_ambiguousTracksVsEta = nullptr;
275  h2_ele_ambiguousTracksVsPhi = nullptr;
276  h2_ele_ambiguousTracksVsPt = nullptr;
277  h1_ele_foundHits = nullptr;
278  h1_ele_foundHits_barrel = nullptr;
279  h1_ele_foundHits_endcaps = nullptr;
280  h2_ele_foundHitsVsEta = nullptr;
281  h2_ele_foundHitsVsPhi = nullptr;
282  h2_ele_foundHitsVsPt = nullptr;
283  h1_ele_lostHits = nullptr;
284  h1_ele_lostHits_barrel = nullptr;
285  h1_ele_lostHits_endcaps = nullptr;
286  h2_ele_lostHitsVsEta = nullptr;
287  h2_ele_lostHitsVsPhi = nullptr;
288  h2_ele_lostHitsVsPt = nullptr;
289  h1_ele_chi2 = nullptr;
290  h1_ele_chi2_barrel = nullptr;
291  h1_ele_chi2_endcaps = nullptr;
292  h2_ele_chi2VsEta = nullptr;
293  h2_ele_chi2VsPhi = nullptr;
294  h2_ele_chi2VsPt = nullptr;
295 
296  h1_ele_PinMnPout = nullptr;
297  h1_ele_PinMnPout_mode = nullptr;
298  h2_ele_PinMnPoutVsEta_mode = nullptr;
299  h2_ele_PinMnPoutVsPhi_mode = nullptr;
300  h2_ele_PinMnPoutVsPt_mode = nullptr;
301  h2_ele_PinMnPoutVsE_mode = nullptr;
302  h2_ele_PinMnPoutVsChi2_mode = nullptr;
303 
304  h1_ele_outerP = nullptr;
305  h1_ele_outerP_mode = nullptr;
306  h2_ele_outerPVsEta_mode = nullptr;
307  h1_ele_outerPt = nullptr;
308  h1_ele_outerPt_mode = nullptr;
309  h2_ele_outerPtVsEta_mode = nullptr;
310  h2_ele_outerPtVsPhi_mode = nullptr;
311  h2_ele_outerPtVsPt_mode = nullptr;
312  h1_ele_EoP = nullptr;
313  h1_ele_EoP_barrel = nullptr;
314  h1_ele_EoP_endcaps = nullptr;
315  h2_ele_EoPVsEta = nullptr;
316  h2_ele_EoPVsPhi = nullptr;
317  h2_ele_EoPVsE = nullptr;
318  h1_ele_EseedOP = nullptr;
319  h1_ele_EseedOP_barrel = nullptr;
320  h1_ele_EseedOP_endcaps = nullptr;
321  h2_ele_EseedOPVsEta = nullptr;
322  h2_ele_EseedOPVsPhi = nullptr;
323  h2_ele_EseedOPVsE = nullptr;
324  h1_ele_EoPout = nullptr;
325  h1_ele_EoPout_barrel = nullptr;
326  h1_ele_EoPout_endcaps = nullptr;
327  h2_ele_EoPoutVsEta = nullptr;
328  h2_ele_EoPoutVsPhi = nullptr;
329  h2_ele_EoPoutVsE = nullptr;
330  h1_ele_EeleOPout = nullptr;
331  h1_ele_EeleOPout_barrel = nullptr;
332  h1_ele_EeleOPout_endcaps = nullptr;
333  h2_ele_EeleOPoutVsEta = nullptr;
334  h2_ele_EeleOPoutVsPhi = nullptr;
335  h2_ele_EeleOPoutVsE = nullptr;
336 
337  h1_ele_dEtaSc_propVtx = nullptr;
338  h1_ele_dEtaSc_propVtx_barrel = nullptr;
339  h1_ele_dEtaSc_propVtx_endcaps = nullptr;
340  h2_ele_dEtaScVsEta_propVtx = nullptr;
341  h2_ele_dEtaScVsPhi_propVtx = nullptr;
342  h2_ele_dEtaScVsPt_propVtx = nullptr;
343  h1_ele_dPhiSc_propVtx = nullptr;
344  h1_ele_dPhiSc_propVtx_barrel = nullptr;
345  h1_ele_dPhiSc_propVtx_endcaps = nullptr;
346  h2_ele_dPhiScVsEta_propVtx = nullptr;
347  h2_ele_dPhiScVsPhi_propVtx = nullptr;
348  h2_ele_dPhiScVsPt_propVtx = nullptr;
349  h1_ele_dEtaCl_propOut = nullptr;
350  h1_ele_dEtaCl_propOut_barrel = nullptr;
351  h1_ele_dEtaCl_propOut_endcaps = nullptr;
352  h2_ele_dEtaClVsEta_propOut = nullptr;
353  h2_ele_dEtaClVsPhi_propOut = nullptr;
354  h2_ele_dEtaClVsPt_propOut = nullptr;
355  h1_ele_dPhiCl_propOut = nullptr;
356  h1_ele_dPhiCl_propOut_barrel = nullptr;
357  h1_ele_dPhiCl_propOut_endcaps = nullptr;
358  h2_ele_dPhiClVsEta_propOut = nullptr;
359  h2_ele_dPhiClVsPhi_propOut = nullptr;
360  h2_ele_dPhiClVsPt_propOut = nullptr;
361  h1_ele_dEtaEleCl_propOut = nullptr;
362  h1_ele_dEtaEleCl_propOut_barrel = nullptr;
363  h1_ele_dEtaEleCl_propOut_endcaps = nullptr;
364  h2_ele_dEtaEleClVsEta_propOut = nullptr;
365  h2_ele_dEtaEleClVsPhi_propOut = nullptr;
366  h2_ele_dEtaEleClVsPt_propOut = nullptr;
367  h1_ele_dPhiEleCl_propOut = nullptr;
368  h1_ele_dPhiEleCl_propOut_barrel = nullptr;
369  h1_ele_dPhiEleCl_propOut_endcaps = nullptr;
370  h2_ele_dPhiEleClVsEta_propOut = nullptr;
371  h2_ele_dPhiEleClVsPhi_propOut = nullptr;
372  h2_ele_dPhiEleClVsPt_propOut = nullptr;
373 
374  h1_ele_seed_subdet2_ = nullptr;
375  h1_ele_seed_mask_ = nullptr;
376  h1_ele_seed_mask_bpix_ = nullptr;
377  h1_ele_seed_mask_fpix_ = nullptr;
378  h1_ele_seed_mask_tec_ = nullptr;
379  h1_ele_seed_dphi2_ = nullptr;
380  h2_ele_seed_dphi2VsEta_ = nullptr;
381  h2_ele_seed_dphi2VsPt_ = nullptr;
382  h1_ele_seed_dphi2pos_ = nullptr;
383  h2_ele_seed_dphi2posVsEta_ = nullptr;
384  h2_ele_seed_dphi2posVsPt_ = nullptr;
385  h1_ele_seed_drz2_ = nullptr;
386  h2_ele_seed_drz2VsEta_ = nullptr;
387  h2_ele_seed_drz2VsPt_ = nullptr;
388  h1_ele_seed_drz2pos_ = nullptr;
389  h2_ele_seed_drz2posVsEta_ = nullptr;
390  h2_ele_seed_drz2posVsPt_ = nullptr;
391 
392  h1_ele_classes = nullptr;
393  h1_ele_eta = nullptr;
394  h1_ele_eta_golden = nullptr;
395  h1_ele_eta_bbrem = nullptr;
396  h1_ele_eta_narrow = nullptr;
397  h1_ele_eta_shower = nullptr;
398 
399  h1_ele_HoE = nullptr;
400  h1_ele_HoE_barrel = nullptr;
401  h1_ele_HoE_endcaps = nullptr;
402  h1_ele_HoE_fiducial = nullptr;
403  h2_ele_HoEVsEta = nullptr;
404  h2_ele_HoEVsPhi = nullptr;
405  h2_ele_HoEVsE = nullptr;
406  // h1_scl_ESFrac = 0 ;
407  h1_scl_ESFrac_endcaps = nullptr;
408 
409  h1_ele_fbrem = nullptr;
410  p1_ele_fbremVsEta_mode = nullptr;
411  p1_ele_fbremVsEta_mean = nullptr;
412  h1_ele_superclusterfbrem = nullptr;
413  h1_ele_superclusterfbrem_barrel = nullptr;
414  h1_ele_superclusterfbrem_endcaps = nullptr;
415  h2_ele_PinVsPoutGolden_mode = nullptr;
416  h2_ele_PinVsPoutShowering_mode = nullptr;
417  h2_ele_PinVsPoutGolden_mean = nullptr;
418  h2_ele_PinVsPoutShowering_mean = nullptr;
419  h2_ele_PtinVsPtoutGolden_mode = nullptr;
420  h2_ele_PtinVsPtoutShowering_mode = nullptr;
421  h2_ele_PtinVsPtoutGolden_mean = nullptr;
422  h2_ele_PtinVsPtoutShowering_mean = nullptr;
423  h1_scl_EoEmatchingObjectGolden_barrel = nullptr;
424  h1_scl_EoEmatchingObjectGolden_endcaps = nullptr;
425  h1_scl_EoEmatchingObjectShowering_barrel = nullptr;
426  h1_scl_EoEmatchingObjectShowering_endcaps = nullptr;
427 
428  h1_ele_mva = nullptr;
429  h1_ele_mva_isolated = nullptr;
430  h1_ele_provenance = nullptr;
431 
432  h1_ele_tkSumPt_dr03 = nullptr;
433  h1_ele_tkSumPt_dr03_barrel = nullptr;
434  h1_ele_tkSumPt_dr03_endcaps = nullptr;
435  h1_ele_ecalRecHitSumEt_dr03 = nullptr;
436  h1_ele_ecalRecHitSumEt_dr03_barrel = nullptr;
437  h1_ele_ecalRecHitSumEt_dr03_endcaps = nullptr;
438  h1_ele_hcalTowerSumEt_dr03_depth1 = nullptr;
439  h1_ele_hcalTowerSumEt_dr03_depth1_barrel = nullptr;
440  h1_ele_hcalTowerSumEt_dr03_depth1_endcaps = nullptr;
441  h1_ele_hcalTowerSumEt_dr03_depth2 = nullptr;
442  h1_ele_tkSumPt_dr04 = nullptr;
443  h1_ele_tkSumPt_dr04_barrel = nullptr;
444  h1_ele_tkSumPt_dr04_endcaps = nullptr;
445  h1_ele_ecalRecHitSumEt_dr04 = nullptr;
446  h1_ele_ecalRecHitSumEt_dr04_barrel = nullptr;
447  h1_ele_ecalRecHitSumEt_dr04_endcaps = nullptr;
448  h1_ele_hcalTowerSumEt_dr04_depth1 = nullptr;
449  h1_ele_hcalTowerSumEt_dr04_depth1_barrel = nullptr;
450  h1_ele_hcalTowerSumEt_dr04_depth1_endcaps = nullptr;
451  h1_ele_hcalTowerSumEt_dr04_depth2 = nullptr;
452 
453  h1_ele_convFlags = nullptr;
454  h1_ele_convFlags_all = nullptr;
455  h1_ele_convDist = nullptr;
456  h1_ele_convDist_all = nullptr;
457  h1_ele_convDcot = nullptr;
458  h1_ele_convDcot_all = nullptr;
459  h1_ele_convRadius = nullptr;
460  h1_ele_convRadius_all = nullptr;
461 }
462 
463 void ElectronMcFakeValidator::bookHistograms(DQMStore::IBooker &iBooker, edm::Run const &, edm::EventSetup const &) {
464  iBooker.setCurrentFolder(outputInternalPath_);
465 
466  setBookIndex(-1);
467  setBookPrefix("h");
468  setBookEfficiencyFlag(set_EfficiencyFlag);
469  setBookStatOverflowFlag(set_StatOverflowFlag);
470 
471  // matching object type
472  std::string matchingObjectType;
473  // Emilia
474  matchingObjectType = "GenJet";
475 
476  std::string htitle = "# " + matchingObjectType + "s", xtitle = "N_{" + matchingObjectType + "}";
477  h1_matchingObjectNum = bookH1withSumw2(iBooker, "matchingObjectNum", htitle, fhits_nbin, 0., fhits_max, xtitle);
478 
479  // rec event collections sizes
480  h1_recEleNum_ = bookH1(iBooker, "recEleNum", "# rec electrons", ele_nbin, ele_min, ele_max, "N_{ele}");
481  h1_recCoreNum_ = bookH1(iBooker, "recCoreNum", "# rec electron cores", core_nbin, core_min, core_max, "N_{core}");
482  h1_recTrackNum_ = bookH1(iBooker, "recTrackNum", "# rec gsf tracks", track_nbin, track_min, track_max, "N_{track}");
483  h1_recSeedNum_ = bookH1(iBooker, "recSeedNum", "# rec electron seeds", seed_nbin, seed_min, seed_max, "N_{seed}");
484  h1_recOfflineVertices_ = bookH1(
485  iBooker, "recOfflineVertices", "# rec Offline Primary Vertices", opv_nbin, opv_min, opv_max, "N_{Vertices}");
486 
487  // mc
488  h1_matchingObjectEta =
489  bookH1withSumw2(iBooker, "matchingObject_eta", matchingObjectType + " #eta", eta_nbin, eta_min, eta_max, "#eta");
490  h1_matchingObjectAbsEta =
491  bookH1withSumw2(iBooker, "matchingObject_abseta", matchingObjectType + " |#eta|", eta_nbin / 2, 0., eta_max);
492  h1_matchingObjectP =
493  bookH1withSumw2(iBooker, "matchingObject_P", matchingObjectType + " p", p_nbin, 0., p_max, "p (GeV/c)");
494  h1_matchingObjectPt =
495  bookH1withSumw2(iBooker, "matchingObject_Pt", matchingObjectType + " pt", pteff_nbin, 5., pt_max);
496  h1_matchingObjectPhi =
497  bookH1withSumw2(iBooker, "matchingObject_phi", matchingObjectType + " phi", phi_nbin, phi_min, phi_max);
498  h1_matchingObjectZ = bookH1withSumw2(iBooker, "matchingObject_z", matchingObjectType + " z", xyz_nbin, -25, 25);
499 
500  setBookPrefix("h_ele");
501 
502  // all electrons
503  h1_ele_EoverP_all = bookH1withSumw2(iBooker,
504  "EoverP_all",
505  "ele E/P_{vertex}, all reco electrons",
506  eop_nbin,
507  0.,
508  eop_max,
509  "E/P_{vertex}",
510  "Events",
511  "ELE_LOGY E1 P");
512  h1_ele_EseedOP_all = bookH1withSumw2(iBooker,
513  "EseedOP_all",
514  "ele E_{seed}/P_{vertex}, all reco electrons",
515  eop_nbin,
516  0.,
517  eop_max,
518  "E_{seed}/P_{vertex}",
519  "Events",
520  "ELE_LOGY E1 P");
521  h1_ele_EoPout_all = bookH1withSumw2(iBooker,
522  "EoPout_all",
523  "ele E_{seed}/P_{out}, all reco electrons",
524  eop_nbin,
525  0.,
526  eop_max,
527  "E_{seed}/P_{out}",
528  "Events",
529  "ELE_LOGY E1 P");
530  h1_ele_EeleOPout_all = bookH1withSumw2(iBooker,
531  "EeleOPout_all",
532  "ele E_{ele}/P_{out}, all reco electrons",
533  eop_nbin,
534  0.,
535  eop_max,
536  "E_{ele}/P_{out}",
537  "Events",
538  "ELE_LOGY E1 P");
539  h1_ele_dEtaSc_propVtx_all = bookH1withSumw2(iBooker,
540  "dEtaSc_propVtx_all",
541  "ele #eta_{sc} - #eta_{tr}, prop from vertex, all reco electrons",
542  detamatch_nbin,
543  detamatch_min,
544  detamatch_max,
545  "#eta_{sc} - #eta_{tr}",
546  "Events",
547  "ELE_LOGY E1 P");
548  h1_ele_dPhiSc_propVtx_all = bookH1withSumw2(iBooker,
549  "dPhiSc_propVtx_all",
550  "ele #phi_{sc} - #phi_{tr}, prop from vertex, all reco electrons",
551  dphimatch_nbin,
552  dphimatch_min,
553  dphimatch_max,
554  "#phi_{sc} - #phi_{tr} (rad)",
555  "Events",
556  "ELE_LOGY E1 P");
557  h1_ele_dEtaCl_propOut_all = bookH1withSumw2(iBooker,
558  "dEtaCl_propOut_all",
559  "ele #eta_{cl} - #eta_{tr}, prop from outermost, all reco electrons",
560  detamatch_nbin,
561  detamatch_min,
562  detamatch_max,
563  "#eta_{sc} - #eta_{tr}",
564  "Events",
565  "ELE_LOGY E1 P");
566  h1_ele_dPhiCl_propOut_all = bookH1withSumw2(iBooker,
567  "dPhiCl_propOut_all",
568  "ele #phi_{cl} - #phi_{tr}, prop from outermost, all reco electrons",
569  dphimatch_nbin,
570  dphimatch_min,
571  dphimatch_max,
572  "#phi_{sc} - #phi_{tr} (rad)",
573  "Events",
574  "ELE_LOGY E1 P");
575  h1_ele_TIP_all = bookH1withSumw2(iBooker,
576  "TIP_all",
577  "ele vertex transverse radius, all reco electrons",
578  100,
579  0.,
580  0.2,
581  "r_{T} (cm)",
582  "Events",
583  "ELE_LOGY E1 P");
584  h1_ele_HoE_all = bookH1withSumw2(iBooker,
585  "HoE_all",
586  "ele hadronic energy / em energy, all reco electrons",
587  hoe_nbin,
588  hoe_min,
589  hoe_max,
590  "H/E",
591  "Events",
592  "ELE_LOGY E1 P");
593  h1_ele_HoE_bc_all = bookH1withSumw2(iBooker,
594  "HoE_bc_all",
595  "ele hadronic energy / em energy, all reco electrons, behind cluster",
596  hoe_nbin,
597  hoe_min,
598  hoe_max,
599  "H/E",
600  "Events",
601  "ELE_LOGY E1 P");
602  h1_ele_vertexEta_all = bookH1withSumw2(iBooker,
603  "vertexEta_all",
604  "ele eta, all reco electrons",
605  eta_nbin,
606  eta_min,
607  eta_max,
608  "",
609  "Events",
610  "ELE_LOGY E1 P");
611  h1_ele_vertexPt_all = bookH1withSumw2(
612  iBooker, "vertexPt_all", "ele p_{T}, all reco electrons", pteff_nbin, 5., pt_max, "", "Events", "ELE_LOGY E1 P");
613  h1_ele_mee_all = bookH1withSumw2(iBooker,
614  "mee_all",
615  "ele pairs invariant mass, all reco electrons",
616  mee_nbin,
617  mee_min,
618  mee_max,
619  "m_{ee} (GeV/c^{2})");
620  h1_ele_mee_os = bookH1withSumw2(iBooker,
621  "mee_os",
622  "ele pairs invariant mass, opp. sign",
623  mee_nbin,
624  mee_min,
625  mee_max,
626  "m_{e^{+}e^{-}} (GeV/c^{2})");
627 
628  // duplicates
629  h2_ele_E2mnE1vsMee_all = bookH2(iBooker,
630  "E2mnE1vsMee_all",
631  "E2 - E1 vs ele pairs invariant mass, all electrons",
632  mee_nbin,
633  mee_min,
634  mee_max,
635  100,
636  -50.,
637  50.,
638  "m_{e^{+}e^{-}} (GeV/c^{2})",
639  "E2 - E1 (GeV)");
640  h2_ele_E2mnE1vsMee_egeg_all = bookH2(iBooker,
641  "E2mnE1vsMee_egeg_all",
642  "E2 - E1 vs ele pairs invariant mass, ecal driven pairs, all electrons",
643  mee_nbin,
644  mee_min,
645  mee_max,
646  100,
647  -50.,
648  50.,
649  "m_{e^{+}e^{-}} (GeV/c^{2})",
650  "E2 - E1 (GeV)");
651 
652  // matched electrons
653 
654  htitle = "Efficiency vs matching " + matchingObjectType + " ";
655  h1_ele_matchingObjectEta_matched =
656  bookH1withSumw2(iBooker, "matchingObjectEta_matched", htitle + "#eta", eta_nbin, eta_min, eta_max);
657  h1_ele_matchingObjectAbsEta_matched =
658  bookH1withSumw2(iBooker, "matchingObjectAbsEta_matched", htitle + "|#eta|", eta_nbin / 2, 0., eta_max);
659  h1_ele_matchingObjectPt_matched =
660  bookH1(iBooker, "matchingObjectPt_matched", htitle + "p_{T}", pteff_nbin, 5., pt_max);
661  h1_ele_matchingObjectPhi_matched =
662  bookH1withSumw2(iBooker, "matchingObjectPhi_matched", htitle + "phi", phi_nbin, phi_min, phi_max);
663  h1_ele_matchingObjectZ_matched = bookH1withSumw2(iBooker, "matchingObjectZ_matched", htitle + "z", xyz_nbin, -25, 25);
664 
665  h1_ele_charge = bookH1withSumw2(iBooker, "charge", "ele charge", 5, -2.5, 2.5, "charge");
666  h2_ele_chargeVsEta = bookH2(iBooker, "chargeVsEta", "ele charge vs eta", eta2D_nbin, eta_min, eta_max, 5, -2., 2.);
667  h2_ele_chargeVsPhi = bookH2(iBooker, "chargeVsPhi", "ele charge vs phi", phi2D_nbin, phi_min, phi_max, 5, -2., 2.);
668  h2_ele_chargeVsPt = bookH2(iBooker, "chargeVsPt", "ele charge vs pt", pt_nbin, 0., 100., 5, -2., 2.);
669  h1_ele_vertexP = bookH1withSumw2(iBooker, "vertexP", "ele momentum", p_nbin, 0., p_max, "p_{vertex} (GeV/c)");
670  h1_ele_vertexPt =
671  bookH1withSumw2(iBooker, "vertexPt", "ele transverse momentum", pt_nbin, 0., pt_max, "p_{T vertex} (GeV/c)");
672  h2_ele_vertexPtVsEta = bookH2(
673  iBooker, "vertexPtVsEta", "ele transverse momentum vs eta", eta2D_nbin, eta_min, eta_max, pt2D_nbin, 0., pt_max);
674  h2_ele_vertexPtVsPhi = bookH2(
675  iBooker, "vertexPtVsPhi", "ele transverse momentum vs phi", phi2D_nbin, phi_min, phi_max, pt2D_nbin, 0., pt_max);
676  h1_ele_vertexEta = bookH1withSumw2(iBooker, "vertexEta", "ele momentum eta", eta_nbin, eta_min, eta_max, "#eta");
677  h2_ele_vertexEtaVsPhi = bookH2(
678  iBooker, "vertexEtaVsPhi", "ele momentum eta vs phi", eta2D_nbin, eta_min, eta_max, phi2D_nbin, phi_min, phi_max);
679  h1_ele_vertexPhi =
680  bookH1withSumw2(iBooker, "vertexPhi", "ele momentum #phi", phi_nbin, phi_min, phi_max, "#phi (rad)");
681  h1_ele_vertexX = bookH1withSumw2(iBooker, "vertexX", "ele vertex x", xyz_nbin, -0.6, 0.6, "x (cm)");
682  h1_ele_vertexY = bookH1withSumw2(iBooker, "vertexY", "ele vertex y", xyz_nbin, -0.6, 0.6, "y (cm)");
683  h1_ele_vertexZ = bookH1withSumw2(iBooker, "vertexZ", "ele vertex z", xyz_nbin, -25, 25, "z (cm)");
684  h1_ele_vertexTIP = bookH1withSumw2(iBooker,
685  "vertexTIP",
686  "ele transverse impact parameter (wrt gen vtx)",
687  90,
688  0.,
689  0.15,
690  "TIP (cm)",
691  "Events",
692  "ELE_LOGY E1 P");
693  h2_ele_vertexTIPVsEta = bookH2(iBooker,
694  "vertexTIPVsEta",
695  "ele transverse impact parameter (wrt gen vtx) vs eta",
696  eta2D_nbin,
697  eta_min,
698  eta_max,
699  45,
700  0.,
701  0.15,
702  "#eta",
703  "TIP (cm)");
704  h2_ele_vertexTIPVsPhi = bookH2(iBooker,
705  "vertexTIPVsPhi",
706  "ele transverse impact parameter (wrt gen vtx) vs phi",
707  phi2D_nbin,
708  phi_min,
709  phi_max,
710  45,
711  0.,
712  0.15,
713  "#phi (rad)",
714  "TIP (cm)");
715  h2_ele_vertexTIPVsPt = bookH2(iBooker,
716  "vertexTIPVsPt",
717  "ele transverse impact parameter (wrt gen vtx) vs transverse momentum",
718  pt2D_nbin,
719  0.,
720  pt_max,
721  45,
722  0.,
723  0.15,
724  "p_{T} (GeV/c)",
725  "TIP (cm)");
726 
727  htitle = "Electron / Matching " + matchingObjectType + ", momemtum";
728  xtitle = "P / P_{" + matchingObjectType + "}";
729  h1_ele_PoPmatchingObject =
730  bookH1withSumw2(iBooker, "PoPmatchingObject", htitle, popmatching_nbin, popmatching_min, popmatching_max, xtitle);
731  h2_ele_PoPmatchingObjectVsEta = bookH2(iBooker,
732  "PoPmatchingObjectVsEta",
733  htitle + ",vs eta",
734  eta2D_nbin,
735  eta_min,
736  eta_max,
737  50,
738  popmatching_min,
739  popmatching_max);
740  h2_ele_PoPmatchingObjectVsPhi = bookH2(iBooker,
741  "PoPmatchingObjectVsPhi",
742  htitle + ",vs phi",
743  phi2D_nbin,
744  phi_min,
745  phi_max,
746  50,
747  popmatching_min,
748  popmatching_max);
749  h2_ele_PoPmatchingObjectVsPt = bookH2(
750  iBooker, "PoPmatchingObjectVsPt", htitle + ",vs eta", pt2D_nbin, 0., pt_max, 50, popmatching_min, popmatching_max);
751  h1_ele_PoPmatchingObject_barrel = bookH1withSumw2(iBooker,
752  "PoPmatchingObject_barrel",
753  htitle + ", barrel",
754  popmatching_nbin,
755  popmatching_min,
756  popmatching_max,
757  xtitle);
758  h1_ele_PoPmatchingObject_endcaps = bookH1withSumw2(iBooker,
759  "PoPmatchingObject_endcaps",
760  htitle + ", endcaps",
761  popmatching_nbin,
762  popmatching_min,
763  popmatching_max,
764  xtitle);
765  htitle = "Ele - " + matchingObjectType + ", ";
766  xtitle = "#eta - #eta_{" + matchingObjectType + "}";
767  h1_ele_EtaMnEtamatchingObject =
768  bookH1withSumw2(iBooker, "EtamatchingObjectEtaTrue", htitle + "eta", deta_nbin, deta_min, deta_max, xtitle);
769  h2_ele_EtaMnEtamatchingObjectVsEta = bookH2(iBooker,
770  "EtaMnEtamatchingObjectVsEta",
771  htitle + "eta, vs eta",
772  eta2D_nbin,
773  eta_min,
774  eta_max,
775  deta_nbin / 2,
776  deta_min,
777  deta_max);
778  h2_ele_EtaMnEtamatchingObjectVsPhi = bookH2(iBooker,
779  "EtaMnEtamatchingObjectVsPhi",
780  htitle + "eta, vs phi",
781  phi2D_nbin,
782  phi_min,
783  phi_max,
784  deta_nbin / 2,
785  deta_min,
786  deta_max);
787  h2_ele_EtaMnEtamatchingObjectVsPt = bookH2(iBooker,
788  "EtaMnEtamatchingObjectVsPt",
789  htitle + "eta,, vs pt",
790  pt_nbin,
791  0.,
792  pt_max,
793  deta_nbin / 2,
794  deta_min,
795  deta_max);
796  xtitle = "#phi - #phi_{" + matchingObjectType + "} (rad)";
797  h1_ele_PhiMnPhimatchingObject =
798  bookH1withSumw2(iBooker, "PhiMnPhimatchingObject", htitle + "phi", dphi_nbin, dphi_min, dphi_max, xtitle);
799  h1_ele_PhiMnPhimatchingObject2 =
800  bookH1(iBooker, "PhiMnPhimatchingObject2", htitle + "phi", dphimatch2D_nbin, dphimatch_min, dphimatch_max);
801  h2_ele_PhiMnPhimatchingObjectVsEta = bookH2(iBooker,
802  "PhiMnPhimatchingObjectVsEta",
803  htitle + "phi, vs eta",
804  eta2D_nbin,
805  eta_min,
806  eta_max,
807  dphi_nbin / 2,
808  dphi_min,
809  dphi_max);
810  h2_ele_PhiMnPhimatchingObjectVsPhi = bookH2(iBooker,
811  "PhiMnPhimatchingObjectVsPhi",
812  htitle + "phi, vs phi",
813  phi2D_nbin,
814  phi_min,
815  phi_max,
816  dphi_nbin / 2,
817  dphi_min,
818  dphi_max);
819  h2_ele_PhiMnPhimatchingObjectVsPt = bookH2(iBooker,
820  "PhiMnPhimatchingObjectVsPt",
821  htitle + "phi, vs pt",
822  pt2D_nbin,
823  0.,
824  pt_max,
825  dphi_nbin / 2,
826  dphi_min,
827  dphi_max);
828 
829  // matched electron, superclusters
830 
831  setBookPrefix("h_scl");
832 
833  h1_scl_En_ = bookH1withSumw2(iBooker, "energy", "ele supercluster energy", p_nbin, 0., p_max);
834  htitle = "Ele supercluster / " + matchingObjectType + ", energy";
835  xtitle = "E/E_{" + matchingObjectType + "}";
836  h1_scl_EoEmatchingObject_barrel =
837  bookH1withSumw2(iBooker, "EoEmatchingObject_barrel", htitle + ", barrel", 50, 0.2, 1.2, xtitle);
838  h1_scl_EoEmatchingObject_endcaps =
839  bookH1withSumw2(iBooker, "EoEmatchingObject_endcaps", htitle + ", endcaps", 50, 0.2, 1.2, xtitle);
840  h1_scl_Et_ = bookH1withSumw2(iBooker, "et", "ele supercluster transverse energy", pt_nbin, 0., pt_max);
841  h2_scl_EtVsEta_ = bookH2(iBooker,
842  "etVsEta",
843  "ele supercluster transverse energy vs eta",
844  eta2D_nbin,
845  eta_min,
846  eta_max,
847  pt_nbin,
848  0.,
849  pt_max);
850  h2_scl_EtVsPhi_ = bookH2(iBooker,
851  "etVsPhi",
852  "ele supercluster transverse energy vs phi",
853  phi2D_nbin,
854  phi_min,
855  phi_max,
856  pt_nbin,
857  0.,
858  pt_max);
859  h2_scl_EtaVsPhi_ = bookH2(
860  iBooker, "etaVsPhi", "ele supercluster eta vs phi", phi2D_nbin, phi_min, phi_max, eta2D_nbin, eta_min, eta_max);
861  h1_scl_Eta_ = bookH1withSumw2(iBooker, "eta", "ele supercluster eta", eta_nbin, eta_min, eta_max);
862  h1_scl_Phi_ = bookH1withSumw2(iBooker, "phi", "ele supercluster phi", phi_nbin, phi_min, phi_max);
863  h1_scl_SigIEtaIEta_ = bookH1withSumw2(iBooker,
864  "sigietaieta",
865  "ele supercluster sigma ieta ieta",
866  100,
867  0.,
868  0.05,
869  "#sigma_{i#eta i#eta}",
870  "Events",
871  "ELE_LOGY E1 P");
872  h1_scl_SigIEtaIEta_barrel_ = bookH1withSumw2(iBooker,
873  "sigietaieta_barrel",
874  "ele supercluster sigma ieta ieta, barrel",
875  100,
876  0.,
877  0.05,
878  "#sigma_{i#eta i#eta}",
879  "Events",
880  "ELE_LOGY E1 P");
881  h1_scl_SigIEtaIEta_endcaps_ = bookH1withSumw2(iBooker,
882  "sigietaieta_endcaps",
883  "ele supercluster sigma ieta ieta, endcaps",
884  100,
885  0.,
886  0.05,
887  "#sigma_{i#eta i#eta}",
888  "Events",
889  "ELE_LOGY E1 P");
890  // new 2014.01.12
891  h1_scl_full5x5_sigmaIetaIeta_ = bookH1withSumw2(iBooker,
892  "full5x5_sigietaieta",
893  "ele supercluster full5x5 sigma ieta ieta",
894  100,
895  0.,
896  0.05,
897  "#sigma_{i#eta i#eta}",
898  "Events",
899  "ELE_LOGY E1 P");
900  h1_scl_full5x5_sigmaIetaIeta_barrel_ = bookH1withSumw2(iBooker,
901  "full5x5_sigietaieta_barrel",
902  "ele supercluster full5x5 sigma ieta ieta, barrel",
903  100,
904  0.,
905  0.05,
906  "#sigma_{i#eta i#eta}",
907  "Events",
908  "ELE_LOGY E1 P");
909  h1_scl_full5x5_sigmaIetaIeta_endcaps_ = bookH1withSumw2(iBooker,
910  "full5x5_sigietaieta_endcaps",
911  "ele supercluster full5x5 sigma ieta ieta, endcaps",
912  100,
913  0.,
914  0.05,
915  "#sigma_{i#eta i#eta}",
916  "Events",
917  "ELE_LOGY E1 P");
918  // new 2014.01.12
919  h1_scl_E1x5_ = bookH1withSumw2(
920  iBooker, "E1x5", "ele supercluster energy in 1x5", p_nbin, 0., p_max, "E1x5 (GeV)", "Events", "ELE_LOGY E1 P");
921  h1_scl_E1x5_barrel_ = bookH1withSumw2(iBooker,
922  "E1x5_barrel",
923  "ele supercluster energy in 1x5 barrel",
924  p_nbin,
925  0.,
926  p_max,
927  "E1x5 (GeV)",
928  "Events",
929  "ELE_LOGY E1 P");
930  h1_scl_E1x5_endcaps_ = bookH1withSumw2(iBooker,
931  "E1x5_endcaps",
932  "ele supercluster energy in 1x5 endcaps",
933  p_nbin,
934  0.,
935  p_max,
936  "E1x5 (GeV)",
937  "Events",
938  "ELE_LOGY E1 P");
939  h1_scl_E2x5max_ = bookH1withSumw2(iBooker,
940  "E2x5max",
941  "ele supercluster energy in 2x5 max",
942  p_nbin,
943  0.,
944  p_max,
945  "E2x5 (GeV)",
946  "Events",
947  "ELE_LOGY E1 P");
948  h1_scl_E2x5max_barrel_ = bookH1withSumw2(iBooker,
949  "E2x5max_barrel",
950  "ele supercluster energy in 2x5 _max barrel",
951  p_nbin,
952  0.,
953  p_max,
954  "E2x5 (GeV)",
955  "Events",
956  "ELE_LOGY E1 P");
957  h1_scl_E2x5max_endcaps_ = bookH1withSumw2(iBooker,
958  "E2x5max_endcaps",
959  "ele supercluster energy in 2x5 _max endcaps",
960  p_nbin,
961  0.,
962  p_max,
963  "E2x5 (GeV)",
964  "Events",
965  "ELE_LOGY E1 P");
966  h1_scl_E5x5_ = bookH1withSumw2(
967  iBooker, "E5x5", "ele supercluster energy in 5x5", p_nbin, 0., p_max, "E5x5 (GeV)", "Events", "ELE_LOGY E1 P");
968  h1_scl_E5x5_barrel_ = bookH1withSumw2(iBooker,
969  "E5x5_barrel",
970  "ele supercluster energy in 5x5 barrel",
971  p_nbin,
972  0.,
973  p_max,
974  "E5x5 (GeV)",
975  "Events",
976  "ELE_LOGY E1 P");
977  h1_scl_E5x5_endcaps_ = bookH1withSumw2(iBooker,
978  "E5x5_endcaps",
979  "ele supercluster energy in 5x5 endcaps",
980  p_nbin,
981  0.,
982  p_max,
983  "E5x5 (GeV)",
984  "Events",
985  "ELE_LOGY E1 P");
986 
987  // matched electron, gsf tracks
988 
989  setBookPrefix("h_ele");
990 
991  h1_ele_ambiguousTracks = bookH1withSumw2(iBooker,
992  "ambiguousTracks",
993  "ele # ambiguous tracks",
994  5,
995  0.,
996  5.,
997  "N_{ambiguous tracks}",
998  "Events",
999  "ELE_LOGY E1 P");
1000  h2_ele_ambiguousTracksVsEta = bookH2(
1001  iBooker, "ambiguousTracksVsEta", "ele # ambiguous tracks vs eta", eta2D_nbin, eta_min, eta_max, 5, 0., 5.);
1002  h2_ele_ambiguousTracksVsPhi = bookH2(
1003  iBooker, "ambiguousTracksVsPhi", "ele # ambiguous tracks vs phi", phi2D_nbin, phi_min, phi_max, 5, 0., 5.);
1004  h2_ele_ambiguousTracksVsPt =
1005  bookH2(iBooker, "ambiguousTracksVsPt", "ele # ambiguous tracks vs pt", pt2D_nbin, 0., pt_max, 5, 0., 5.);
1006  h1_ele_foundHits =
1007  bookH1withSumw2(iBooker, "foundHits", "ele track # found hits", fhits_nbin, 0., fhits_max, "N_{hits}");
1008  h2_ele_foundHitsVsEta = bookH2(iBooker,
1009  "foundHitsVsEta",
1010  "ele track # found hits vs eta",
1011  eta2D_nbin,
1012  eta_min,
1013  eta_max,
1014  fhits_nbin,
1015  0.,
1016  fhits_max);
1017  h2_ele_foundHitsVsPhi = bookH2(iBooker,
1018  "foundHitsVsPhi",
1019  "ele track # found hits vs phi",
1020  phi2D_nbin,
1021  phi_min,
1022  phi_max,
1023  fhits_nbin,
1024  0.,
1025  fhits_max);
1026  h2_ele_foundHitsVsPt = bookH2(
1027  iBooker, "foundHitsVsPt", "ele track # found hits vs pt", pt2D_nbin, 0., pt_max, fhits_nbin, 0., fhits_max);
1028  h1_ele_lostHits = bookH1withSumw2(iBooker, "lostHits", "ele track # lost hits", 5, 0., 5., "N_{lost hits}");
1029  h2_ele_lostHitsVsEta = bookH2(
1030  iBooker, "lostHitsVsEta", "ele track # lost hits vs eta", eta2D_nbin, eta_min, eta_max, lhits_nbin, 0., lhits_max);
1031  h2_ele_lostHitsVsPhi = bookH2(
1032  iBooker, "lostHitsVsPhi", "ele track # lost hits vs eta", phi2D_nbin, phi_min, phi_max, lhits_nbin, 0., lhits_max);
1033  h2_ele_lostHitsVsPt =
1034  bookH2(iBooker, "lostHitsVsPt", "ele track # lost hits vs eta", pt2D_nbin, 0., pt_max, lhits_nbin, 0., lhits_max);
1035  h1_ele_chi2 =
1036  bookH1withSumw2(iBooker, "chi2", "ele track #chi^{2}", 100, 0., 15., "#Chi^{2}", "Events", "ELE_LOGY E1 P");
1037  h2_ele_chi2VsEta =
1038  bookH2(iBooker, "chi2VsEta", "ele track #chi^{2} vs eta", eta2D_nbin, eta_min, eta_max, 50, 0., 15.);
1039  h2_ele_chi2VsPhi =
1040  bookH2(iBooker, "chi2VsPhi", "ele track #chi^{2} vs phi", phi2D_nbin, phi_min, phi_max, 50, 0., 15.);
1041  h2_ele_chi2VsPt = bookH2(iBooker, "chi2VsPt", "ele track #chi^{2} vs pt", pt2D_nbin, 0., pt_max, 50, 0., 15.);
1042  h1_ele_PinMnPout = bookH1withSumw2(iBooker,
1043  "PinMnPout",
1044  "ele track inner p - outer p, mean of GSF components",
1045  p_nbin,
1046  0.,
1047  200.,
1048  "P_{vertex} - P_{out} (GeV/c)");
1049  h1_ele_PinMnPout_mode = bookH1withSumw2(iBooker,
1050  "PinMnPout_mode",
1051  "ele track inner p - outer p, mode of GSF components",
1052  p_nbin,
1053  0.,
1054  100.,
1055  "P_{vertex} - P_{out}, mode of GSF components (GeV/c)");
1056  h2_ele_PinMnPoutVsEta_mode = bookH2(iBooker,
1057  "PinMnPoutVsEta_mode",
1058  "ele track inner p - outer p vs eta, mode of GSF components",
1059  eta2D_nbin,
1060  eta_min,
1061  eta_max,
1062  p2D_nbin,
1063  0.,
1064  100.);
1065  h2_ele_PinMnPoutVsPhi_mode = bookH2(iBooker,
1066  "PinMnPoutVsPhi_mode",
1067  "ele track inner p - outer p vs phi, mode of GSF components",
1068  phi2D_nbin,
1069  phi_min,
1070  phi_max,
1071  p2D_nbin,
1072  0.,
1073  100.);
1074  h2_ele_PinMnPoutVsPt_mode = bookH2(iBooker,
1075  "PinMnPoutVsPt_mode",
1076  "ele track inner p - outer p vs pt, mode of GSF components",
1077  pt2D_nbin,
1078  0.,
1079  pt_max,
1080  p2D_nbin,
1081  0.,
1082  100.);
1083  h2_ele_PinMnPoutVsE_mode = bookH2(iBooker,
1084  "PinMnPoutVsE_mode",
1085  "ele track inner p - outer p vs E, mode of GSF components",
1086  p2D_nbin,
1087  0.,
1088  200.,
1089  p2D_nbin,
1090  0.,
1091  100.);
1092  h2_ele_PinMnPoutVsChi2_mode = bookH2(iBooker,
1093  "PinMnPoutVsChi2_mode",
1094  "ele track inner p - outer p vs track chi2, mode of GSF components",
1095  50,
1096  0.,
1097  20.,
1098  p2D_nbin,
1099  0.,
1100  100.);
1101  h1_ele_outerP = bookH1withSumw2(
1102  iBooker, "outerP", "ele track outer p, mean of GSF components", p_nbin, 0., p_max, "P_{out} (GeV/c)");
1103  h1_ele_outerP_mode = bookH1withSumw2(
1104  iBooker, "outerP_mode", "ele track outer p, mode of GSF components", p_nbin, 0., p_max, "P_{out} (GeV/c)");
1105  h2_ele_outerPVsEta_mode =
1106  bookH2(iBooker, "outerPVsEta_mode", "ele track outer p vs eta mode", eta2D_nbin, eta_min, eta_max, 50, 0., p_max);
1107  h1_ele_outerPt = bookH1withSumw2(
1108  iBooker, "outerPt", "ele track outer p_{T}, mean of GSF components", pt_nbin, 0., pt_max, "P_{T out} (GeV/c)");
1109  h1_ele_outerPt_mode = bookH1withSumw2(iBooker,
1110  "outerPt_mode",
1111  "ele track outer p_{T}, mode of GSF components",
1112  pt_nbin,
1113  0.,
1114  pt_max,
1115  "P_{T out} (GeV/c)");
1116  h2_ele_outerPtVsEta_mode = bookH2(iBooker,
1117  "outerPtVsEta_mode",
1118  "ele track outer p_{T} vs eta, mode of GSF components",
1119  eta2D_nbin,
1120  eta_min,
1121  eta_max,
1122  pt2D_nbin,
1123  0.,
1124  pt_max);
1125  h2_ele_outerPtVsPhi_mode = bookH2(iBooker,
1126  "outerPtVsPhi_mode",
1127  "ele track outer p_{T} vs phi, mode of GSF components",
1128  phi2D_nbin,
1129  phi_min,
1130  phi_max,
1131  pt2D_nbin,
1132  0.,
1133  pt_max);
1134  h2_ele_outerPtVsPt_mode = bookH2(iBooker,
1135  "outerPtVsPt_mode",
1136  "ele track outer p_{T} vs pt, mode of GSF components",
1137  pt2D_nbin,
1138  0.,
1139  100.,
1140  pt2D_nbin,
1141  0.,
1142  pt_max);
1143 
1144  // matched electrons, matching
1145  h1_ele_EoP = bookH1withSumw2(
1146  iBooker, "EoP", "ele E/P_{vertex}", eop_nbin, 0., eop_max, "E/P_{vertex}", "Events", "ELE_LOGY E1 P");
1147  h1_ele_EoP_barrel = bookH1withSumw2(iBooker,
1148  "EoP_barrel",
1149  "ele E/P_{vertex} barrel",
1150  eop_nbin,
1151  0.,
1152  eop_max,
1153  "E/P_{vertex}",
1154  "Events",
1155  "ELE_LOGY E1 P");
1156  h1_ele_EoP_endcaps = bookH1withSumw2(iBooker,
1157  "EoP_endcaps",
1158  "ele E/P_{vertex} endcaps",
1159  eop_nbin,
1160  0.,
1161  eop_max,
1162  "E/P_{vertex}",
1163  "Events",
1164  "ELE_LOGY E1 P");
1165  h2_ele_EoPVsEta =
1166  bookH2(iBooker, "EoPVsEta", "ele E/P_{vertex} vs eta", eta2D_nbin, eta_min, eta_max, eop2D_nbin, 0., eopmaxsht);
1167  h2_ele_EoPVsPhi =
1168  bookH2(iBooker, "EoPVsPhi", "ele E/P_{vertex} vs phi", phi2D_nbin, phi_min, phi_max, eop2D_nbin, 0., eopmaxsht);
1169  h2_ele_EoPVsE = bookH2(iBooker, "EoPVsE", "ele E/P_{vertex} vs E", 50, 0., p_max, 50, 0., 5.);
1170  h1_ele_EseedOP = bookH1withSumw2(iBooker,
1171  "EseedOP",
1172  "ele E_{seed}/P_{vertex}",
1173  eop_nbin,
1174  0.,
1175  eop_max,
1176  "E_{seed}/P_{vertex}",
1177  "Events",
1178  "ELE_LOGY E1 P");
1179  h1_ele_EseedOP_barrel = bookH1withSumw2(iBooker,
1180  "EseedOP_barrel",
1181  "ele E_{seed}/P_{vertex} barrel",
1182  eop_nbin,
1183  0.,
1184  eop_max,
1185  "E_{seed}/P_{vertex}",
1186  "Events",
1187  "ELE_LOGY E1 P");
1188  h1_ele_EseedOP_endcaps = bookH1withSumw2(iBooker,
1189  "EseedOP_endcaps",
1190  "ele E_{seed}/P_{vertex} endcaps",
1191  eop_nbin,
1192  0.,
1193  eop_max,
1194  "E_{seed}/P_{vertex}",
1195  "Events",
1196  "ELE_LOGY E1 P");
1197  h2_ele_EseedOPVsEta = bookH2(iBooker,
1198  "EseedOPVsEta",
1199  "ele E_{seed}/P_{vertex} vs eta",
1200  eta2D_nbin,
1201  eta_min,
1202  eta_max,
1203  eop2D_nbin,
1204  0.,
1205  eopmaxsht);
1206  h2_ele_EseedOPVsPhi = bookH2(iBooker,
1207  "EseedOPVsPhi",
1208  "ele E_{seed}/P_{vertex} vs phi",
1209  phi2D_nbin,
1210  phi_min,
1211  phi_max,
1212  eop2D_nbin,
1213  0.,
1214  eopmaxsht);
1215  h2_ele_EseedOPVsE = bookH2(iBooker, "EseedOPVsE", "ele E_{seed}/P_{vertex} vs E", 50, 0., p_max, 50, 0., 5.);
1216  h1_ele_EoPout = bookH1withSumw2(
1217  iBooker, "EoPout", "ele E_{seed}/P_{out}", eop_nbin, 0., eop_max, "E_{seed}/P_{out}", "Events", "ELE_LOGY E1 P");
1218  h1_ele_EoPout_barrel = bookH1withSumw2(iBooker,
1219  "EoPout_barrel",
1220  "ele E_{seed}/P_{out} barrel",
1221  eop_nbin,
1222  0.,
1223  eop_max,
1224  "E_{seed}/P_{out}",
1225  "Events",
1226  "ELE_LOGY E1 P");
1227  h1_ele_EoPout_endcaps = bookH1withSumw2(iBooker,
1228  "EoPout_endcaps",
1229  "ele E_{seed}/P_{out} endcaps",
1230  eop_nbin,
1231  0.,
1232  eop_max,
1233  "E_{seed}/P_{out}",
1234  "Events",
1235  "ELE_LOGY E1 P");
1236  h2_ele_EoPoutVsEta = bookH2(
1237  iBooker, "EoPoutVsEta", "ele E_{seed}/P_{out} vs eta", eta2D_nbin, eta_min, eta_max, eop2D_nbin, 0., eopmaxsht);
1238  h2_ele_EoPoutVsPhi = bookH2(
1239  iBooker, "EoPoutVsPhi", "ele E_{seed}/P_{out} vs phi", phi2D_nbin, phi_min, phi_max, eop2D_nbin, 0., eopmaxsht);
1240  h2_ele_EoPoutVsE =
1241  bookH2(iBooker, "EoPoutVsE", "ele E_{seed}/P_{out} vs E", p2D_nbin, 0., p_max, eop2D_nbin, 0., eopmaxsht);
1242  h1_ele_EeleOPout = bookH1withSumw2(
1243  iBooker, "EeleOPout", "ele E_{ele}/P_{out}", eop_nbin, 0., eop_max, "E_{ele}/P_{out}", "Events", "ELE_LOGY E1 P");
1244  h1_ele_EeleOPout_barrel = bookH1withSumw2(iBooker,
1245  "EeleOPout_barrel",
1246  "ele E_{ele}/P_{out} barrel",
1247  eop_nbin,
1248  0.,
1249  eop_max,
1250  "E_{ele}/P_{out}",
1251  "Events",
1252  "ELE_LOGY E1 P");
1253  h1_ele_EeleOPout_endcaps = bookH1withSumw2(iBooker,
1254  "EeleOPout_endcaps",
1255  "ele E_{ele}/P_{out} endcaps",
1256  eop_nbin,
1257  0.,
1258  eop_max,
1259  "E_{ele}/P_{out}",
1260  "Events",
1261  "ELE_LOGY E1 P");
1262  h2_ele_EeleOPoutVsEta = bookH2(
1263  iBooker, "EeleOPoutVsEta", "ele E_{ele}/P_{out} vs eta", eta2D_nbin, eta_min, eta_max, eop2D_nbin, 0., eopmaxsht);
1264  h2_ele_EeleOPoutVsPhi = bookH2(
1265  iBooker, "EeleOPoutVsPhi", "ele E_{ele}/P_{out} vs phi", phi2D_nbin, phi_min, phi_max, eop2D_nbin, 0., eopmaxsht);
1266  h2_ele_EeleOPoutVsE =
1267  bookH2(iBooker, "EeleOPoutVsE", "ele E_{ele}/P_{out} vs E", p2D_nbin, 0., p_max, eop2D_nbin, 0., eopmaxsht);
1268  h1_ele_dEtaSc_propVtx = bookH1withSumw2(iBooker,
1269  "dEtaSc_propVtx",
1270  "ele #eta_{sc} - #eta_{tr}, prop from vertex",
1271  detamatch_nbin,
1272  detamatch_min,
1273  detamatch_max,
1274  "#eta_{sc} - #eta_{tr}",
1275  "Events",
1276  "ELE_LOGY E1 P");
1277  h1_ele_dEtaSc_propVtx_barrel = bookH1withSumw2(iBooker,
1278  "dEtaSc_propVtx_barrel",
1279  "ele #eta_{sc} - #eta_{tr}, prop from vertex, barrel",
1280  detamatch_nbin,
1281  detamatch_min,
1282  detamatch_max,
1283  "#eta_{sc} - #eta_{tr}",
1284  "Events",
1285  "ELE_LOGY E1 P");
1286  h1_ele_dEtaSc_propVtx_endcaps = bookH1withSumw2(iBooker,
1287  "dEtaSc_propVtx_endcaps",
1288  "ele #eta_{sc} - #eta_{tr}, prop from vertex, endcaps",
1289  detamatch_nbin,
1290  detamatch_min,
1291  detamatch_max,
1292  "#eta_{sc} - #eta_{tr}",
1293  "Events",
1294  "ELE_LOGY E1 P");
1295  h2_ele_dEtaScVsEta_propVtx = bookH2(iBooker,
1296  "dEtaScVsEta_propVtx",
1297  "ele #eta_{sc} - #eta_{tr} vs eta, prop from vertex",
1298  eta2D_nbin,
1299  eta_min,
1300  eta_max,
1301  detamatch2D_nbin,
1302  detamatch_min,
1303  detamatch_max);
1304  h2_ele_dEtaScVsPhi_propVtx = bookH2(iBooker,
1305  "dEtaScVsPhi_propVtx",
1306  "ele #eta_{sc} - #eta_{tr} vs phi, prop from vertex",
1307  phi2D_nbin,
1308  phi_min,
1309  phi_max,
1310  detamatch2D_nbin,
1311  detamatch_min,
1312  detamatch_max);
1313  h2_ele_dEtaScVsPt_propVtx = bookH2(iBooker,
1314  "dEtaScVsPt_propVtx",
1315  "ele #eta_{sc} - #eta_{tr} vs pt, prop from vertex",
1316  pt2D_nbin,
1317  0.,
1318  pt_max,
1319  detamatch2D_nbin,
1320  detamatch_min,
1321  detamatch_max);
1322  h1_ele_dPhiSc_propVtx = bookH1withSumw2(iBooker,
1323  "dPhiSc_propVtx",
1324  "ele #phi_{sc} - #phi_{tr}, prop from vertex",
1325  dphimatch_nbin,
1326  dphimatch_min,
1327  dphimatch_max,
1328  "#phi_{sc} - #phi_{tr} (rad)",
1329  "Events",
1330  "ELE_LOGY E1 P");
1331  h1_ele_dPhiSc_propVtx_barrel = bookH1withSumw2(iBooker,
1332  "dPhiSc_propVtx_barrel",
1333  "ele #phi_{sc} - #phi_{tr}, prop from vertex, barrel",
1334  dphimatch_nbin,
1335  dphimatch_min,
1336  dphimatch_max,
1337  "#phi_{sc} - #phi_{tr} (rad)",
1338  "Events",
1339  "ELE_LOGY E1 P");
1340  h1_ele_dPhiSc_propVtx_endcaps = bookH1withSumw2(iBooker,
1341  "dPhiSc_propVtx_endcaps",
1342  "ele #phi_{sc} - #phi_{tr}, prop from vertex, endcaps",
1343  dphimatch_nbin,
1344  dphimatch_min,
1345  dphimatch_max,
1346  "#phi_{sc} - #phi_{tr} (rad)",
1347  "Events",
1348  "ELE_LOGY E1 P");
1349  h2_ele_dPhiScVsEta_propVtx = bookH2(iBooker,
1350  "dPhiScVsEta_propVtx",
1351  "ele #phi_{sc} - #phi_{tr} vs eta, prop from vertex",
1352  eta2D_nbin,
1353  eta_min,
1354  eta_max,
1355  dphimatch2D_nbin,
1356  dphimatch_min,
1357  dphimatch_max);
1358  h2_ele_dPhiScVsPhi_propVtx = bookH2(iBooker,
1359  "dPhiScVsPhi_propVtx",
1360  "ele #phi_{sc} - #phi_{tr} vs phi, prop from vertex",
1361  phi2D_nbin,
1362  phi_min,
1363  phi_max,
1364  dphimatch2D_nbin,
1365  dphimatch_min,
1366  dphimatch_max);
1367  h2_ele_dPhiScVsPt_propVtx = bookH2(iBooker,
1368  "dPhiScVsPt_propVtx",
1369  "ele #phi_{sc} - #phi_{tr} vs pt, prop from vertex",
1370  pt2D_nbin,
1371  0.,
1372  pt_max,
1373  dphimatch2D_nbin,
1374  dphimatch_min,
1375  dphimatch_max);
1376  h1_ele_dEtaCl_propOut = bookH1withSumw2(iBooker,
1377  "dEtaCl_propOut",
1378  "ele #eta_{cl} - #eta_{tr}, prop from outermost",
1379  detamatch_nbin,
1380  detamatch_min,
1381  detamatch_max,
1382  "#eta_{seedcl} - #eta_{tr}",
1383  "Events",
1384  "ELE_LOGY E1 P");
1385  h1_ele_dEtaCl_propOut_barrel = bookH1withSumw2(iBooker,
1386  "dEtaCl_propOut_barrel",
1387  "ele #eta_{cl} - #eta_{tr}, prop from outermost, barrel",
1388  detamatch_nbin,
1389  detamatch_min,
1390  detamatch_max,
1391  "#eta_{seedcl} - #eta_{tr}",
1392  "Events",
1393  "ELE_LOGY E1 P");
1394  h1_ele_dEtaCl_propOut_endcaps = bookH1withSumw2(iBooker,
1395  "dEtaCl_propOut_endcaps",
1396  "ele #eta_{cl} - #eta_{tr}, prop from outermost, endcaps",
1397  detamatch_nbin,
1398  detamatch_min,
1399  detamatch_max,
1400  "#eta_{seedcl} - #eta_{tr}",
1401  "Events",
1402  "ELE_LOGY E1 P");
1403  h2_ele_dEtaClVsEta_propOut = bookH2(iBooker,
1404  "dEtaClVsEta_propOut",
1405  "ele #eta_{cl} - #eta_{tr} vs eta, prop from out",
1406  eta2D_nbin,
1407  eta_min,
1408  eta_max,
1409  detamatch2D_nbin,
1410  detamatch_min,
1411  detamatch_max);
1412  h2_ele_dEtaClVsPhi_propOut = bookH2(iBooker,
1413  "dEtaClVsPhi_propOut",
1414  "ele #eta_{cl} - #eta_{tr} vs phi, prop from out",
1415  phi2D_nbin,
1416  phi_min,
1417  phi_max,
1418  detamatch2D_nbin,
1419  detamatch_min,
1420  detamatch_max);
1421  h2_ele_dEtaClVsPt_propOut = bookH2(iBooker,
1422  "dEtaScVsPt_propOut",
1423  "ele #eta_{cl} - #eta_{tr} vs pt, prop from out",
1424  pt2D_nbin,
1425  0.,
1426  pt_max,
1427  detamatch2D_nbin,
1428  detamatch_min,
1429  detamatch_max);
1430  h1_ele_dPhiCl_propOut = bookH1withSumw2(iBooker,
1431  "dPhiCl_propOut",
1432  "ele #phi_{cl} - #phi_{tr}, prop from outermost",
1433  dphimatch_nbin,
1434  dphimatch_min,
1435  dphimatch_max,
1436  "#phi_{seedcl} - #phi_{tr} (rad)",
1437  "Events",
1438  "ELE_LOGY E1 P");
1439  h1_ele_dPhiCl_propOut_barrel = bookH1withSumw2(iBooker,
1440  "dPhiCl_propOut_barrel",
1441  "ele #phi_{cl} - #phi_{tr}, prop from outermost, barrel",
1442  dphimatch_nbin,
1443  dphimatch_min,
1444  dphimatch_max,
1445  "#phi_{seedcl} - #phi_{tr} (rad)",
1446  "Events",
1447  "ELE_LOGY E1 P");
1448  h1_ele_dPhiCl_propOut_endcaps = bookH1withSumw2(iBooker,
1449  "dPhiCl_propOut_endcaps",
1450  "ele #phi_{cl} - #phi_{tr}, prop from outermost, endcaps",
1451  dphimatch_nbin,
1452  dphimatch_min,
1453  dphimatch_max,
1454  "#phi_{seedcl} - #phi_{tr} (rad)",
1455  "Events",
1456  "ELE_LOGY E1 P");
1457  h2_ele_dPhiClVsEta_propOut = bookH2(iBooker,
1458  "dPhiClVsEta_propOut",
1459  "ele #phi_{cl} - #phi_{tr} vs eta, prop from out",
1460  eta2D_nbin,
1461  eta_min,
1462  eta_max,
1463  dphimatch2D_nbin,
1464  dphimatch_min,
1465  dphimatch_max);
1466  h2_ele_dPhiClVsPhi_propOut = bookH2(iBooker,
1467  "dPhiClVsPhi_propOut",
1468  "ele #phi_{cl} - #phi_{tr} vs phi, prop from out",
1469  phi2D_nbin,
1470  phi_min,
1471  phi_max,
1472  dphimatch2D_nbin,
1473  dphimatch_min,
1474  dphimatch_max);
1475  h2_ele_dPhiClVsPt_propOut = bookH2(iBooker,
1476  "dPhiSClsPt_propOut",
1477  "ele #phi_{cl} - #phi_{tr} vs pt, prop from out",
1478  pt2D_nbin,
1479  0.,
1480  pt_max,
1481  dphimatch2D_nbin,
1482  dphimatch_min,
1483  dphimatch_max);
1484  h1_ele_dEtaEleCl_propOut = bookH1withSumw2(iBooker,
1485  "dEtaEleCl_propOut",
1486  "ele #eta_{EleCl} - #eta_{tr}, prop from outermost",
1487  detamatch_nbin,
1488  detamatch_min,
1489  detamatch_max,
1490  "#eta_{elecl} - #eta_{tr}",
1491  "Events",
1492  "ELE_LOGY E1 P");
1493  h1_ele_dEtaEleCl_propOut_barrel = bookH1withSumw2(iBooker,
1494  "dEtaEleCl_propOut_barrel",
1495  "ele #eta_{EleCl} - #eta_{tr}, prop from outermost, barrel",
1496  detamatch_nbin,
1497  detamatch_min,
1498  detamatch_max,
1499  "#eta_{elecl} - #eta_{tr}",
1500  "Events",
1501  "ELE_LOGY E1 P");
1502  h1_ele_dEtaEleCl_propOut_endcaps = bookH1withSumw2(iBooker,
1503  "dEtaEleCl_propOut_endcaps",
1504  "ele #eta_{EleCl} - #eta_{tr}, prop from outermost, endcaps",
1505  detamatch_nbin,
1506  detamatch_min,
1507  detamatch_max,
1508  "#eta_{elecl} - #eta_{tr}",
1509  "Events",
1510  "ELE_LOGY E1 P");
1511  h2_ele_dEtaEleClVsEta_propOut = bookH2(iBooker,
1512  "dEtaEleClVsEta_propOut",
1513  "ele #eta_{EleCl} - #eta_{tr} vs eta, prop from out",
1514  eta2D_nbin,
1515  eta_min,
1516  eta_max,
1517  detamatch2D_nbin,
1518  detamatch_min,
1519  detamatch_max);
1520  h2_ele_dEtaEleClVsPhi_propOut = bookH2(iBooker,
1521  "dEtaEleClVsPhi_propOut",
1522  "ele #eta_{EleCl} - #eta_{tr} vs phi, prop from out",
1523  phi2D_nbin,
1524  phi_min,
1525  phi_max,
1526  detamatch2D_nbin,
1527  detamatch_min,
1528  detamatch_max);
1529  h2_ele_dEtaEleClVsPt_propOut = bookH2(iBooker,
1530  "dEtaScVsPt_propOut",
1531  "ele #eta_{EleCl} - #eta_{tr} vs pt, prop from out",
1532  pt2D_nbin,
1533  0.,
1534  pt_max,
1535  detamatch2D_nbin,
1536  detamatch_min,
1537  detamatch_max);
1538  h1_ele_dPhiEleCl_propOut = bookH1withSumw2(iBooker,
1539  "dPhiEleCl_propOut",
1540  "ele #phi_{EleCl} - #phi_{tr}, prop from outermost",
1541  dphimatch_nbin,
1542  dphimatch_min,
1543  dphimatch_max,
1544  "#phi_{elecl} - #phi_{tr} (rad)",
1545  "Events",
1546  "ELE_LOGY E1 P");
1547  h1_ele_dPhiEleCl_propOut_barrel = bookH1withSumw2(iBooker,
1548  "dPhiEleCl_propOut_barrel",
1549  "ele #phi_{EleCl} - #phi_{tr}, prop from outermost, barrel",
1550  dphimatch_nbin,
1551  dphimatch_min,
1552  dphimatch_max,
1553  "#phi_{elecl} - #phi_{tr} (rad)",
1554  "Events",
1555  "ELE_LOGY E1 P");
1556  h1_ele_dPhiEleCl_propOut_endcaps = bookH1withSumw2(iBooker,
1557  "dPhiEleCl_propOut_endcaps",
1558  "ele #phi_{EleCl} - #phi_{tr}, prop from outermost, endcaps",
1559  dphimatch_nbin,
1560  dphimatch_min,
1561  dphimatch_max,
1562  "#phi_{elecl} - #phi_{tr} (rad)",
1563  "Events",
1564  "ELE_LOGY E1 P");
1565  h2_ele_dPhiEleClVsEta_propOut = bookH2(iBooker,
1566  "dPhiEleClVsEta_propOut",
1567  "ele #phi_{EleCl} - #phi_{tr} vs eta, prop from out",
1568  eta2D_nbin,
1569  eta_min,
1570  eta_max,
1571  dphimatch2D_nbin,
1572  dphimatch_min,
1573  dphimatch_max);
1574  h2_ele_dPhiEleClVsPhi_propOut = bookH2(iBooker,
1575  "dPhiEleClVsPhi_propOut",
1576  "ele #phi_{EleCl} - #phi_{tr} vs phi, prop from out",
1577  phi2D_nbin,
1578  phi_min,
1579  phi_max,
1580  dphimatch2D_nbin,
1581  dphimatch_min,
1582  dphimatch_max);
1583  h2_ele_dPhiEleClVsPt_propOut = bookH2(iBooker,
1584  "dPhiSEleClsPt_propOut",
1585  "ele #phi_{EleCl} - #phi_{tr} vs pt, prop from out",
1586  pt2D_nbin,
1587  0.,
1588  pt_max,
1589  dphimatch2D_nbin,
1590  dphimatch_min,
1591  dphimatch_max);
1592  h1_ele_HoE = bookH1withSumw2(
1593  iBooker, "HoE", "ele hadronic energy / em energy", hoe_nbin, hoe_min, hoe_max, "H/E", "Events", "ELE_LOGY E1 P");
1594  h1_ele_HoE_barrel = bookH1withSumw2(iBooker,
1595  "HoE_barrel",
1596  "ele hadronic energy / em energy, barrel",
1597  hoe_nbin,
1598  hoe_min,
1599  hoe_max,
1600  "H/E",
1601  "Events",
1602  "ELE_LOGY E1 P");
1603  h1_ele_HoE_endcaps = bookH1withSumw2(iBooker,
1604  "HoE_endcaps",
1605  "ele hadronic energy / em energy, endcaps",
1606  hoe_nbin,
1607  hoe_min,
1608  hoe_max,
1609  "H/E",
1610  "Events",
1611  "ELE_LOGY E1 P");
1612  h1_ele_HoE_bc = bookH1withSumw2(iBooker,
1613  "HoE_bc",
1614  "ele hadronic energy / em energy behind cluster",
1615  hoe_nbin,
1616  hoe_min,
1617  hoe_max,
1618  "H/E",
1619  "Events",
1620  "ELE_LOGY E1 P");
1621  h1_ele_HoE_bc_barrel = bookH1withSumw2(iBooker,
1622  "HoE_bc_barrel",
1623  "ele hadronic energy / em energy, behind cluster barrel",
1624  hoe_nbin,
1625  hoe_min,
1626  hoe_max,
1627  "H/E",
1628  "Events",
1629  "ELE_LOGY E1 P");
1630  h1_ele_HoE_bc_endcaps = bookH1withSumw2(iBooker,
1631  "HoE_bc_endcaps",
1632  "ele hadronic energy / em energy, behind cluster, endcaps",
1633  hoe_nbin,
1634  hoe_min,
1635  hoe_max,
1636  "H/E",
1637  "Events",
1638  "ELE_LOGY E1 P");
1639  h1_ele_hcalDepth1OverEcalBc = bookH1withSumw2(iBooker,
1640  "hcalDepth1OverEcalBc",
1641  "hcalDepth1OverEcalBc",
1642  hoe_nbin,
1643  hoe_min,
1644  hoe_max,
1645  "H/E",
1646  "Events",
1647  "ELE_LOGY E1 P");
1648  h1_ele_hcalDepth1OverEcalBc_barrel = bookH1withSumw2(iBooker,
1649  "hcalDepth1OverEcalBc_barrel",
1650  "hcalDepth1OverEcalBc_barrel",
1651  hoe_nbin,
1652  hoe_min,
1653  hoe_max,
1654  "H/E",
1655  "Events",
1656  "ELE_LOGY E1 P");
1657  h1_ele_hcalDepth1OverEcalBc_endcaps = bookH1withSumw2(iBooker,
1658  "hcalDepth1OverEcalBc_endcaps",
1659  "hcalDepth1OverEcalBc_endcaps",
1660  hoe_nbin,
1661  hoe_min,
1662  hoe_max,
1663  "H/E",
1664  "Events",
1665  "ELE_LOGY E1 P");
1666  h1_ele_hcalDepth2OverEcalBc = bookH1withSumw2(iBooker,
1667  "hcalDepth2OverEcalBc",
1668  "hcalDepth2OverEcalBc",
1669  hoe_nbin,
1670  hoe_min,
1671  hoe_max,
1672  "H/E",
1673  "Events",
1674  "ELE_LOGY E1 P");
1675  h1_ele_hcalDepth2OverEcalBc_barrel = bookH1withSumw2(iBooker,
1676  "hcalDepth2OverEcalBc_barrel",
1677  "hcalDepth2OverEcalBc_barrel",
1678  hoe_nbin,
1679  hoe_min,
1680  hoe_max,
1681  "H/E",
1682  "Events",
1683  "ELE_LOGY E1 P");
1684  h1_ele_hcalDepth2OverEcalBc_endcaps = bookH1withSumw2(iBooker,
1685  "hcalDepth2OverEcalBc_endcaps",
1686  "hcalDepth2OverEcalBc_endcaps",
1687  hoe_nbin,
1688  hoe_min,
1689  hoe_max,
1690  "H/E",
1691  "Events",
1692  "ELE_LOGY E1 P");
1693  h1_ele_HoE_fiducial = bookH1withSumw2(iBooker,
1694  "HoE_fiducial",
1695  "ele hadronic energy / em energy, fiducial region",
1696  hoe_nbin,
1697  hoe_min,
1698  hoe_max,
1699  "H/E",
1700  "Events",
1701  "ELE_LOGY E1 P");
1702  h2_ele_HoEVsEta = bookH2(iBooker,
1703  "HoEVsEta",
1704  "ele hadronic energy / em energy vs eta",
1705  eta_nbin,
1706  eta_min,
1707  eta_max,
1708  hoe_nbin,
1709  hoe_min,
1710  hoe_max);
1711  h2_ele_HoEVsPhi = bookH2(iBooker,
1712  "HoEVsPhi",
1713  "ele hadronic energy / em energy vs phi",
1714  phi2D_nbin,
1715  phi_min,
1716  phi_max,
1717  hoe_nbin,
1718  hoe_min,
1719  hoe_max);
1720  h2_ele_HoEVsE =
1721  bookH2(iBooker, "HoEVsE", "ele hadronic energy / em energy vs E", p_nbin, 0., 300., hoe_nbin, hoe_min, hoe_max);
1722  setBookPrefix("h_scl");
1723  // h1_scl_ESFrac = bookH1withSumw2(iBooker, "ESFrac","Preshower over SC raw energy",100,0.,0.8,"E_{PS} / E^{raw}_{SC}","Events","ELE_LOGY E1 P");
1724  h1_scl_ESFrac_endcaps = bookH1withSumw2(iBooker,
1725  "ESFrac_endcaps",
1726  "Preshower over SC raw energy , endcaps",
1727  100,
1728  0.,
1729  0.8,
1730  "E_{PS} / E^{raw}_{SC}",
1731  "Events",
1732  "ELE_LOGY E1 P");
1733 
1734  // seeds
1735  setBookPrefix("h_ele");
1736  h1_ele_seed_subdet2_ =
1737  bookH1withSumw2(iBooker, "seedSubdet2", "ele seed subdet 2nd layer", 11, -0.5, 10.5, "2nd hit subdet Id");
1738  h1_ele_seed_mask_ = bookH1withSumw2(iBooker, "seedMask", "ele seed hits mask", 13, -0.5, 12.5);
1739  h1_ele_seed_mask_bpix_ =
1740  bookH1withSumw2(iBooker, "seedMask_Bpix", "ele seed hits mask when subdet2 is bpix", 13, -0.5, 12.5);
1741  h1_ele_seed_mask_fpix_ =
1742  bookH1withSumw2(iBooker, "seedMask_Fpix", "ele seed hits mask when subdet2 is bpix", 13, -0.5, 12.5);
1743  h1_ele_seed_mask_tec_ =
1744  bookH1withSumw2(iBooker, "seedMask_Tec", "ele seed hits mask when subdet2 is bpix", 13, -0.5, 12.5);
1745  h1_ele_seed_dphi2_ = bookH1withSumw2(
1746  iBooker, "seedDphi2", "ele seed dphi 2nd layer", 50, -0.003, +0.003, "#phi_{hit}-#phi_{pred} (rad)");
1747  h2_ele_seed_dphi2VsEta_ = bookH2(
1748  iBooker, "seedDphi2_VsEta", "ele seed dphi 2nd layer vs eta", eta2D_nbin, eta_min, eta_max, 50, -0.003, +0.003);
1749  h2_ele_seed_dphi2VsPt_ =
1750  bookH2(iBooker, "seedDphi2_VsPt", "ele seed dphi 2nd layer vs pt", pt2D_nbin, 0., pt_max, 50, -0.003, +0.003);
1751  h1_ele_seed_dphi2pos_ = bookH1withSumw2(
1752  iBooker, "seedDphi2Pos", "ele seed dphi 2nd layer positron", 50, -0.003, +0.003, "#phi_{hit}-#phi_{pred} (rad)");
1753  h2_ele_seed_dphi2posVsEta_ = bookH2(iBooker,
1754  "seedDphi2Pos_VsEta",
1755  "ele seed dphi 2nd layer positron vs eta",
1756  eta2D_nbin,
1757  eta_min,
1758  eta_max,
1759  50,
1760  -0.003,
1761  +0.003);
1762  h2_ele_seed_dphi2posVsPt_ = bookH2(
1763  iBooker, "seedDphi2Pos_VsPt", "ele seed dphi 2nd layer positron vs pt", pt2D_nbin, 0., pt_max, 50, -0.003, +0.003);
1764  h1_ele_seed_drz2_ = bookH1withSumw2(
1765  iBooker, "seedDrz2", "ele seed dr (dz) 2nd layer", 50, -0.03, +0.03, "r(z)_{hit}-r(z)_{pred} (cm)");
1766  h2_ele_seed_drz2VsEta_ = bookH2(
1767  iBooker, "seedDrz2_VsEta", "ele seed dr/dz 2nd layer vs eta", eta2D_nbin, eta_min, eta_max, 50, -0.03, +0.03);
1768  h2_ele_seed_drz2VsPt_ =
1769  bookH2(iBooker, "seedDrz2_VsPt", "ele seed dr/dz 2nd layer vs pt", pt2D_nbin, 0., pt_max, 50, -0.03, +0.03);
1770  h1_ele_seed_drz2pos_ = bookH1withSumw2(
1771  iBooker, "seedDrz2Pos", "ele seed dr (dz) 2nd layer positron", 50, -0.03, +0.03, "r(z)_{hit}-r(z)_{pred} (cm)");
1772  h2_ele_seed_drz2posVsEta_ = bookH2(iBooker,
1773  "seedDrz2Pos_VsEta",
1774  "ele seed dr/dz 2nd layer positron vs eta",
1775  eta2D_nbin,
1776  eta_min,
1777  eta_max,
1778  50,
1779  -0.03,
1780  +0.03);
1781  h2_ele_seed_drz2posVsPt_ = bookH2(
1782  iBooker, "seedDrz2Pos_VsPt", "ele seed dr/dz 2nd layer positron vs pt", pt2D_nbin, 0., pt_max, 50, -0.03, +0.03);
1783 
1784  // classes
1785  h1_ele_classes = bookH1withSumw2(iBooker, "classes", "ele classes", 20, 0.0, 20., "class Id");
1786  h1_ele_eta = bookH1withSumw2(iBooker, "eta", "ele electron eta", eta_nbin / 2, 0.0, eta_max);
1787  h1_ele_eta_golden = bookH1withSumw2(iBooker, "eta_golden", "ele electron eta golden", eta_nbin / 2, 0.0, eta_max);
1788  h1_ele_eta_bbrem = bookH1withSumw2(iBooker, "eta_bbrem", "ele electron eta bbrem", eta_nbin / 2, 0.0, eta_max);
1789  h1_ele_eta_shower = bookH1withSumw2(iBooker, "eta_shower", "ele electron eta showering", eta_nbin / 2, 0.0, eta_max);
1790  h2_ele_PinVsPoutGolden_mode = bookH2(iBooker,
1791  "PinVsPoutGolden_mode",
1792  "ele track inner p vs outer p vs eta, golden, mode of GSF components",
1793  p2D_nbin,
1794  0.,
1795  p_max,
1796  50,
1797  0.,
1798  p_max);
1799  h2_ele_PinVsPoutShowering_mode = bookH2(iBooker,
1800  "PinVsPoutShowering_mode",
1801  "ele track inner p vs outer p vs eta, showering, mode of GSF components",
1802  p2D_nbin,
1803  0.,
1804  p_max,
1805  50,
1806  0.,
1807  p_max);
1808  h2_ele_PinVsPoutGolden_mean = bookH2(iBooker,
1809  "PinVsPoutGolden_mean",
1810  "ele track inner p vs outer p vs eta, golden, mean of GSF components",
1811  p2D_nbin,
1812  0.,
1813  p_max,
1814  50,
1815  0.,
1816  p_max);
1817  h2_ele_PinVsPoutShowering_mean = bookH2(iBooker,
1818  "PinVsPoutShowering_mean",
1819  "ele track inner p vs outer p vs eta, showering, mean of GSF components",
1820  p2D_nbin,
1821  0.,
1822  p_max,
1823  50,
1824  0.,
1825  p_max);
1826  h2_ele_PtinVsPtoutGolden_mode = bookH2(iBooker,
1827  "PtinVsPtoutGolden_mode",
1828  "ele track inner pt vs outer pt vs eta, golden, mode of GSF components",
1829  pt2D_nbin,
1830  0.,
1831  pt_max,
1832  50,
1833  0.,
1834  pt_max);
1835  h2_ele_PtinVsPtoutShowering_mode = bookH2(iBooker,
1836  "PtinVsPtoutShowering_mode",
1837  "ele track inner pt vs outer pt vs eta, showering, mode of GSF components",
1838  pt2D_nbin,
1839  0.,
1840  pt_max,
1841  50,
1842  0.,
1843  pt_max);
1844  h2_ele_PtinVsPtoutGolden_mean = bookH2(iBooker,
1845  "PtinVsPtoutGolden_mean",
1846  "ele track inner pt vs outer pt vs eta, golden, mean of GSF components",
1847  pt2D_nbin,
1848  0.,
1849  pt_max,
1850  50,
1851  0.,
1852  pt_max);
1853  h2_ele_PtinVsPtoutShowering_mean = bookH2(iBooker,
1854  "PtinVsPtoutShowering_mean",
1855  "ele track inner pt vs outer pt vs eta, showering, mean of GSF components",
1856  pt2D_nbin,
1857  0.,
1858  pt_max,
1859  50,
1860  0.,
1861  pt_max);
1862  setBookPrefix("h_scl");
1863  h1_scl_EoEmatchingObjectGolden_barrel = bookH1withSumw2(iBooker,
1864  "EoEmatchingObject_golden_barrel",
1865  "ele supercluster energy / gen energy, golden, barrel",
1866  popmatching_nbin,
1867  popmatching_min,
1868  popmatching_max);
1869  h1_scl_EoEmatchingObjectGolden_endcaps = bookH1withSumw2(iBooker,
1870  "EoEmatchingObject_golden_endcaps",
1871  "ele supercluster energy / gen energy, golden, endcaps",
1872  popmatching_nbin,
1873  popmatching_min,
1874  popmatching_max);
1875  h1_scl_EoEmatchingObjectShowering_barrel = bookH1withSumw2(iBooker,
1876  "EoEmatchingObject_showering_barrel",
1877  "ele supercluster energy / gen energy, showering, barrel",
1878  popmatching_nbin,
1879  popmatching_min,
1880  popmatching_max);
1881  h1_scl_EoEmatchingObjectShowering_endcaps =
1882  bookH1withSumw2(iBooker,
1883  "EoEmatchingObject_showering_endcaps",
1884  "ele supercluster energy / gen energy, showering, endcaps",
1885  popmatching_nbin,
1886  popmatching_min,
1887  popmatching_max);
1888  setBookPrefix("h_ele");
1889 
1890  // isolation
1891  h1_ele_tkSumPt_dr03 = bookH1withSumw2(iBooker,
1892  "tkSumPt_dr03",
1893  "tk isolation sum, dR=0.3",
1894  100,
1895  0.0,
1896  20.,
1897  "TkIsoSum, cone 0.3 (GeV/c)",
1898  "Events",
1899  "ELE_LOGY E1 P");
1900  h1_ele_tkSumPt_dr03_barrel = bookH1withSumw2(iBooker,
1901  "tkSumPt_dr03_barrel",
1902  "tk isolation sum, dR=0.3, barrel",
1903  100,
1904  0.0,
1905  20.,
1906  "TkIsoSum, cone 0.3 (GeV/c)",
1907  "Events",
1908  "ELE_LOGY E1 P");
1909  h1_ele_tkSumPt_dr03_endcaps = bookH1withSumw2(iBooker,
1910  "tkSumPt_dr03_endcaps",
1911  "tk isolation sum, dR=0.3, endcaps",
1912  100,
1913  0.0,
1914  20.,
1915  "TkIsoSum, cone 0.3 (GeV/c)",
1916  "Events",
1917  "ELE_LOGY E1 P");
1918  h1_ele_ecalRecHitSumEt_dr03 = bookH1withSumw2(iBooker,
1919  "ecalRecHitSumEt_dr03",
1920  "ecal isolation sum, dR=0.3",
1921  100,
1922  0.0,
1923  20.,
1924  "EcalIsoSum, cone 0.3 (GeV)",
1925  "Events",
1926  "ELE_LOGY E1 P");
1927  h1_ele_ecalRecHitSumEt_dr03_barrel = bookH1withSumw2(iBooker,
1928  "ecalRecHitSumEt_dr03_barrel",
1929  "ecal isolation sum, dR=0.3, barrel",
1930  100,
1931  0.0,
1932  20.,
1933  "EcalIsoSum, cone 0.3 (GeV)",
1934  "Events",
1935  "ELE_LOGY E1 P");
1936  h1_ele_ecalRecHitSumEt_dr03_endcaps = bookH1withSumw2(iBooker,
1937  "ecalRecHitSumEt_dr03_endcaps",
1938  "ecal isolation sum, dR=0.3, endcaps",
1939  100,
1940  0.0,
1941  20.,
1942  "EcalIsoSum, cone 0.3 (GeV)",
1943  "Events",
1944  "ELE_LOGY E1 P");
1945  h1_ele_hcalTowerSumEt_dr03_depth1 = bookH1withSumw2(iBooker,
1946  "hcalTowerSumEt_dr03_depth1",
1947  "hcal depth1 isolation sum, dR=0.3",
1948  100,
1949  0.0,
1950  20.,
1951  "Hcal1IsoSum, cone 0.3 (GeV)",
1952  "Events",
1953  "ELE_LOGY E1 P");
1954  h1_ele_hcalTowerSumEt_dr03_depth1_barrel = bookH1withSumw2(iBooker,
1955  "hcalTowerSumEt_dr03_depth1_barrel",
1956  "hcal depth1 isolation sum, dR=0.3, barrel",
1957  100,
1958  0.0,
1959  20.,
1960  "Hcal1IsoSum, cone 0.3 (GeV)",
1961  "Events",
1962  "ELE_LOGY E1 P");
1963  h1_ele_hcalTowerSumEt_dr03_depth1_endcaps = bookH1withSumw2(iBooker,
1964  "hcalTowerSumEt_dr03_depth1_endcaps",
1965  "hcal depth1 isolation sum, dR=0.3, endcaps",
1966  100,
1967  0.0,
1968  20.,
1969  "Hcal1IsoSum, cone 0.3 (GeV)",
1970  "Events",
1971  "ELE_LOGY E1 P");
1972  h1_ele_hcalTowerSumEt_dr03_depth2 = bookH1withSumw2(iBooker,
1973  "hcalTowerSumEt_dr03_depth2",
1974  "hcal depth2 isolation sum, dR=0.3",
1975  100,
1976  0.0,
1977  20.,
1978  "Hcal2IsoSum, cone 0.3 (GeV)",
1979  "Events",
1980  "ELE_LOGY E1 P");
1981  h1_ele_tkSumPt_dr04 = bookH1withSumw2(iBooker,
1982  "tkSumPt_dr04",
1983  "tk isolation sum, dR=0.4",
1984  100,
1985  0.0,
1986  20.,
1987  "TkIsoSum, cone 0.4 (GeV/c)",
1988  "Events",
1989  "ELE_LOGY E1 P");
1990  h1_ele_tkSumPt_dr04_barrel = bookH1withSumw2(iBooker,
1991  "tkSumPt_dr04_barrel",
1992  "tk isolation sum, dR=0.4, barrel",
1993  100,
1994  0.0,
1995  20.,
1996  "TkIsoSum, cone 0.4 (GeV/c)",
1997  "Events",
1998  "ELE_LOGY E1 P");
1999  h1_ele_tkSumPt_dr04_endcaps = bookH1withSumw2(iBooker,
2000  "tkSumPt_dr04_endcaps",
2001  "tk isolation sum, dR=0.4, endcaps",
2002  100,
2003  0.0,
2004  20.,
2005  "TkIsoSum, cone 0.4 (GeV/c)",
2006  "Events",
2007  "ELE_LOGY E1 P");
2008  h1_ele_ecalRecHitSumEt_dr04 = bookH1withSumw2(iBooker,
2009  "ecalRecHitSumEt_dr04",
2010  "ecal isolation sum, dR=0.4",
2011  100,
2012  0.0,
2013  20.,
2014  "EcalIsoSum, cone 0.4 (GeV)",
2015  "Events",
2016  "ELE_LOGY E1 P");
2017  h1_ele_ecalRecHitSumEt_dr04_barrel = bookH1withSumw2(iBooker,
2018  "ecalRecHitSumEt_dr04_barrel",
2019  "ecal isolation sum, dR=0.4, barrel",
2020  100,
2021  0.0,
2022  20.,
2023  "EcalIsoSum, cone 0.4 (GeV)",
2024  "Events",
2025  "ELE_LOGY E1 P");
2026  h1_ele_ecalRecHitSumEt_dr04_endcaps = bookH1withSumw2(iBooker,
2027  "ecalRecHitSumEt_dr04_endcaps",
2028  "ecal isolation sum, dR=0.4, endcaps",
2029  100,
2030  0.0,
2031  20.,
2032  "EcalIsoSum, cone 0.4 (GeV)",
2033  "Events",
2034  "ELE_LOGY E1 P");
2035  h1_ele_hcalTowerSumEt_dr04_depth1 = bookH1withSumw2(iBooker,
2036  "hcalTowerSumEt_dr04_depth1",
2037  "hcal depth1 isolation sum, dR=0.4",
2038  100,
2039  0.0,
2040  20.,
2041  "Hcal1IsoSum, cone 0.4 (GeV)",
2042  "Events",
2043  "ELE_LOGY E1 P");
2044  h1_ele_hcalTowerSumEt_dr04_depth1_barrel = bookH1withSumw2(iBooker,
2045  "hcalTowerSumEt_dr04_depth1_barrel",
2046  "hcal depth1 isolation sum, dR=0.4, barrel",
2047  100,
2048  0.0,
2049  20.,
2050  "Hcal1IsoSum, cone 0.4 (GeV)",
2051  "Events",
2052  "ELE_LOGY E1 P");
2053  h1_ele_hcalTowerSumEt_dr04_depth1_endcaps = bookH1withSumw2(iBooker,
2054  "hcalTowerSumEt_dr04_depth1_endcaps",
2055  "hcal depth1 isolation sum, dR=0.4, endcaps",
2056  100,
2057  0.0,
2058  20.,
2059  "Hcal1IsoSum, cone 0.4 (GeV)",
2060  "Events",
2061  "ELE_LOGY E1 P");
2062  h1_ele_hcalTowerSumEt_dr04_depth2 = bookH1withSumw2(iBooker,
2063  "hcalTowerSumEt_dr04_depth2",
2064  "hcal depth2 isolation sum, dR=0.4",
2065  100,
2066  0.0,
2067  20.,
2068  "Hcal2IsoSum, cone 0.4 (GeV)",
2069  "Events",
2070  "ELE_LOGY E1 P");
2071 
2072  // newHCAL
2073  // isolation new hcal
2074  h1_ele_hcalTowerSumEtBc_dr03_depth1 = bookH1withSumw2(iBooker,
2075  "hcalTowerSumEtBc_dr03_depth1",
2076  "hcal depth1 isolation sum behind cluster, dR=0.3",
2077  100,
2078  0.0,
2079  20.,
2080  "Hcal1IsoSum, cone 0.3 (GeV)",
2081  "Events",
2082  "ELE_LOGY E1 P");
2083  h1_ele_hcalTowerSumEtBc_dr03_depth1_barrel =
2084  bookH1withSumw2(iBooker,
2085  "hcalTowerSumEtBc_dr03_depth1_barrel",
2086  "hcal depth1 isolation sum behind cluster, dR=0.3, barrel",
2087  100,
2088  0.0,
2089  20.,
2090  "Hcal1IsoSum, cone 0.3 (GeV)",
2091  "Events",
2092  "ELE_LOGY E1 P");
2093  h1_ele_hcalTowerSumEtBc_dr03_depth1_endcaps =
2094  bookH1withSumw2(iBooker,
2095  "hcalTowerSumEtBc_dr03_depth1_endcaps",
2096  "hcal depth1 isolation sum behind cluster, dR=0.3, endcaps",
2097  100,
2098  0.0,
2099  20.,
2100  "Hcal1IsoSum, cone 0.3 (GeV)",
2101  "Events",
2102  "ELE_LOGY E1 P");
2103 
2104  h1_ele_hcalTowerSumEtBc_dr04_depth1 = bookH1withSumw2(iBooker,
2105  "hcalTowerSumEtBc_dr04_depth1",
2106  "hcal depth1 isolation sum behind cluster, dR=0.4",
2107  100,
2108  0.0,
2109  20.,
2110  "Hcal1IsoSum, cone 0.4 (GeV)",
2111  "Events",
2112  "ELE_LOGY E1 P");
2113  h1_ele_hcalTowerSumEtBc_dr04_depth1_barrel =
2114  bookH1withSumw2(iBooker,
2115  "hcalTowerSumEtBc_dr04_depth1_barrel",
2116  "hcal depth1 isolation sum behind cluster, dR=0.4, barrel",
2117  100,
2118  0.0,
2119  20.,
2120  "Hcal1IsoSum, cone 0.4 (GeV)",
2121  "Events",
2122  "ELE_LOGY E1 P");
2123  h1_ele_hcalTowerSumEtBc_dr04_depth1_endcaps =
2124  bookH1withSumw2(iBooker,
2125  "hcalTowerSumEtBc_dr04_depth1_endcaps",
2126  "hcal depth1 isolation sum behind cluster, dR=0.4, endcaps",
2127  100,
2128  0.0,
2129  20.,
2130  "Hcal1IsoSum, cone 0.4 (GeV)",
2131  "Events",
2132  "ELE_LOGY E1 P");
2133 
2134  h1_ele_hcalTowerSumEtBc_dr03_depth2 = bookH1withSumw2(iBooker,
2135  "hcalTowerSumEtBc_dr03_depth2",
2136  "hcal depth2 isolation sum behind cluster, dR=0.3",
2137  100,
2138  0.0,
2139  20.,
2140  "Hcal1IsoSum, cone 0.3 (GeV)",
2141  "Events",
2142  "ELE_LOGY E1 P");
2143  h1_ele_hcalTowerSumEtBc_dr03_depth2_barrel =
2144  bookH1withSumw2(iBooker,
2145  "hcalTowerSumEtBc_dr03_depth2_barrel",
2146  "hcal depth2 isolation sum behind cluster, dR=0.3, barrel",
2147  100,
2148  0.0,
2149  20.,
2150  "Hcal1IsoSum, cone 0.3 (GeV)",
2151  "Events",
2152  "ELE_LOGY E1 P");
2153  h1_ele_hcalTowerSumEtBc_dr03_depth2_endcaps =
2154  bookH1withSumw2(iBooker,
2155  "hcalTowerSumEtBc_dr03_depth2_endcaps",
2156  "hcal depth2 isolation sum behind cluster, dR=0.3, endcaps",
2157  100,
2158  0.0,
2159  20.,
2160  "Hcal1IsoSum, cone 0.3 (GeV)",
2161  "Events",
2162  "ELE_LOGY E1 P");
2163 
2164  h1_ele_hcalTowerSumEtBc_dr04_depth2 = bookH1withSumw2(iBooker,
2165  "hcalTowerSumEtBc_dr04_depth2",
2166  "hcal depth2 isolation sum behind cluster, dR=0.4",
2167  100,
2168  0.0,
2169  20.,
2170  "Hcal1IsoSum, cone 0.4 (GeV)",
2171  "Events",
2172  "ELE_LOGY E1 P");
2173  h1_ele_hcalTowerSumEtBc_dr04_depth2_barrel =
2174  bookH1withSumw2(iBooker,
2175  "hcalTowerSumEtBc_dr04_depth2_barrel",
2176  "hcal depth2 isolation sum behind cluster, dR=0.4, barrel",
2177  100,
2178  0.0,
2179  20.,
2180  "Hcal1IsoSum, cone 0.4 (GeV)",
2181  "Events",
2182  "ELE_LOGY E1 P");
2183  h1_ele_hcalTowerSumEtBc_dr04_depth2_endcaps =
2184  bookH1withSumw2(iBooker,
2185  "hcalTowerSumEtBc_dr04_depth2_endcaps",
2186  "hcal depth2 isolation sum behind cluster, dR=0.4, endcaps",
2187  100,
2188  0.0,
2189  20.,
2190  "Hcal1IsoSum, cone 0.4 (GeV)",
2191  "Events",
2192  "ELE_LOGY E1 P");
2193 
2194  // fbrem
2195  h1_ele_fbrem = bookH1withSumw2(
2196  iBooker, "fbrem", "ele brem fraction, mode of GSF components", 100, 0., 1., "P_{in} - P_{out} / P_{in}");
2197  h1_ele_fbrem_barrel = bookH1withSumw2(iBooker,
2198  "fbrem_barrel",
2199  "ele brem fraction for barrel, mode of GSF components",
2200  100,
2201  0.,
2202  1.,
2203  "P_{in} - P_{out} / P_{in}");
2204  h1_ele_fbrem_endcaps = bookH1withSumw2(iBooker,
2205  "fbrem_endcaps",
2206  "ele brem franction for endcaps, mode of GSF components",
2207  100,
2208  0.,
2209  1.,
2210  "P_{in} - P_{out} / P_{in}");
2211  p1_ele_fbremVsEta_mode = bookP1(iBooker,
2212  "fbremvsEtamode",
2213  "mean ele brem fraction vs eta, mode of GSF components",
2214  eta2D_nbin,
2215  eta_min,
2216  eta_max,
2217  0.,
2218  1.,
2219  "#eta",
2220  "<P_{in} - P_{out} / P_{in}>");
2221  p1_ele_fbremVsEta_mean = bookP1(iBooker,
2222  "fbremvsEtamean",
2223  "mean ele brem fraction vs eta, mean of GSF components",
2224  eta2D_nbin,
2225  eta_min,
2226  eta_max,
2227  0.,
2228  1.,
2229  "#eta",
2230  "<P_{in} - P_{out} / P_{in}>");
2231  h1_ele_superclusterfbrem =
2232  bookH1withSumw2(iBooker, "superclusterfbrem", "supercluster brem fraction", 100, 0., 1., "1 - E_{ele} / E_{SC}");
2233  h1_ele_superclusterfbrem_barrel = bookH1withSumw2(
2234  iBooker, "superclusterfbrem_barrel", "supercluster brem fraction for barrel", 100, 0., 1., "1 - E_{ele} / E_{SC}");
2235  h1_ele_superclusterfbrem_endcaps = bookH1withSumw2(iBooker,
2236  "superclusterfbrem_endcaps",
2237  "supercluster brem franction for endcaps",
2238  100,
2239  0.,
2240  1.,
2241  "1 - E_{ele} / E_{SC}");
2242  // e/g et pflow electrons
2243  h1_ele_mva = bookH1withSumw2(iBooker, "mva", "ele identification mva", 100, -1., 1.);
2244  h1_ele_mva_barrel = bookH1withSumw2(iBooker, "mva_barrel", "ele identification mva barrel", 100, -1., 1.);
2245  h1_ele_mva_endcaps = bookH1withSumw2(iBooker, "mva_endcaps", "ele identification mva endcaps", 100, -1., 1.);
2246  h1_ele_mva_isolated = bookH1withSumw2(iBooker, "mva_isolated", "ele identification mva isolated", 100, -1., 1.);
2247  h1_ele_mva_barrel_isolated =
2248  bookH1withSumw2(iBooker, "mva_isolated_barrel", "ele identification mva isolated barrel", 100, -1., 1.);
2249  h1_ele_mva_endcaps_isolated =
2250  bookH1withSumw2(iBooker, "mva_isolated_endcaps", "ele identification mva isolated endcaps", 100, -1., 1.);
2251  h1_ele_provenance = bookH1withSumw2(iBooker, "provenance", "ele provenance", 5, -2., 3.);
2252  h1_ele_provenance_barrel = bookH1withSumw2(iBooker, "provenance_barrel", "ele provenance barrel", 5, -2., 3.);
2253  h1_ele_provenance_endcaps = bookH1withSumw2(iBooker, "provenance_endcaps", "ele provenance endcaps", 5, -2., 3.);
2254 
2255  // pflow isolation variables
2256  h1_ele_chargedHadronIso = bookH1withSumw2(
2257  iBooker, "chargedHadronIso", "chargedHadronIso", 100, 0.0, 20., "chargedHadronIso", "Events", "ELE_LOGY E1 P");
2258  h1_ele_chargedHadronIso_barrel = bookH1withSumw2(iBooker,
2259  "chargedHadronIso_barrel",
2260  "chargedHadronIso for barrel",
2261  100,
2262  0.0,
2263  20.,
2264  "chargedHadronIso_barrel",
2265  "Events",
2266  "ELE_LOGY E1 P");
2267  h1_ele_chargedHadronIso_endcaps = bookH1withSumw2(iBooker,
2268  "chargedHadronIso_endcaps",
2269  "chargedHadronIso for endcaps",
2270  100,
2271  0.0,
2272  20.,
2273  "chargedHadronIso_endcaps",
2274  "Events",
2275  "ELE_LOGY E1 P");
2276  h1_ele_neutralHadronIso = bookH1withSumw2(
2277  iBooker, "neutralHadronIso", "neutralHadronIso", 21, 0.0, 20., "neutralHadronIso", "Events", "ELE_LOGY E1 P");
2278  h1_ele_neutralHadronIso_barrel = bookH1withSumw2(iBooker,
2279  "neutralHadronIso_barrel",
2280  "neutralHadronIso for barrel",
2281  21,
2282  0.0,
2283  20.,
2284  "neutralHadronIso_barrel",
2285  "Events",
2286  "ELE_LOGY E1 P");
2287  h1_ele_neutralHadronIso_endcaps = bookH1withSumw2(iBooker,
2288  "neutralHadronIso_endcaps",
2289  "neutralHadronIso for endcaps",
2290  21,
2291  0.0,
2292  20.,
2293  "neutralHadronIso_endcaps",
2294  "Events",
2295  "ELE_LOGY E1 P");
2296  h1_ele_photonIso =
2297  bookH1withSumw2(iBooker, "photonIso", "photonIso", 100, 0.0, 20., "photonIso", "Events", "ELE_LOGY E1 P");
2298  h1_ele_photonIso_barrel = bookH1withSumw2(
2299  iBooker, "photonIso_barrel", "photonIso for barrel", 100, 0.0, 20., "photonIso_barrel", "Events", "ELE_LOGY E1 P");
2300  h1_ele_photonIso_endcaps = bookH1withSumw2(iBooker,
2301  "photonIso_endcaps",
2302  "photonIso for endcaps",
2303  100,
2304  0.0,
2305  20.,
2306  "photonIso_endcaps",
2307  "Events",
2308  "ELE_LOGY E1 P");
2309  // -- pflow over pT
2310  h1_ele_chargedHadronRelativeIso = bookH1withSumw2(iBooker,
2311  "chargedHadronRelativeIso",
2312  "chargedHadronRelativeIso",
2313  100,
2314  0.0,
2315  2.,
2316  "chargedHadronRelativeIso",
2317  "Events",
2318  "ELE_LOGY E1 P");
2319  h1_ele_chargedHadronRelativeIso_barrel = bookH1withSumw2(iBooker,
2320  "chargedHadronRelativeIso_barrel",
2321  "chargedHadronRelativeIso for barrel",
2322  100,
2323  0.0,
2324  2.,
2325  "chargedHadronRelativeIso_barrel",
2326  "Events",
2327  "ELE_LOGY E1 P");
2328  h1_ele_chargedHadronRelativeIso_endcaps = bookH1withSumw2(iBooker,
2329  "chargedHadronRelativeIso_endcaps",
2330  "chargedHadronRelativeIso for endcaps",
2331  100,
2332  0.0,
2333  2.,
2334  "chargedHadronRelativeIso_endcaps",
2335  "Events",
2336  "ELE_LOGY E1 P");
2337  h1_ele_neutralHadronRelativeIso = bookH1withSumw2(iBooker,
2338  "neutralHadronRelativeIso",
2339  "neutralHadronRelativeIso",
2340  100,
2341  0.0,
2342  2.,
2343  "neutralHadronRelativeIso",
2344  "Events",
2345  "ELE_LOGY E1 P");
2346  h1_ele_neutralHadronRelativeIso_barrel = bookH1withSumw2(iBooker,
2347  "neutralHadronRelativeIso_barrel",
2348  "neutralHadronRelativeIso for barrel",
2349  100,
2350  0.0,
2351  2.,
2352  "neutralHadronRelativeIso_barrel",
2353  "Events",
2354  "ELE_LOGY E1 P");
2355  h1_ele_neutralHadronRelativeIso_endcaps = bookH1withSumw2(iBooker,
2356  "neutralHadronRelativeIso_endcaps",
2357  "neutralHadronRelativeIso for endcaps",
2358  100,
2359  0.0,
2360  2.,
2361  "neutralHadronRelativeIso_endcaps",
2362  "Events",
2363  "ELE_LOGY E1 P");
2364  h1_ele_photonRelativeIso = bookH1withSumw2(
2365  iBooker, "photonRelativeIso", "photonRelativeIso", 100, 0.0, 2., "photonRelativeIso", "Events", "ELE_LOGY E1 P");
2366  h1_ele_photonRelativeIso_barrel = bookH1withSumw2(iBooker,
2367  "photonRelativeIso_barrel",
2368  "photonRelativeIso for barrel",
2369  100,
2370  0.0,
2371  2.,
2372  "photonRelativeIso_barrel",
2373  "Events",
2374  "ELE_LOGY E1 P");
2375  h1_ele_photonRelativeIso_endcaps = bookH1withSumw2(iBooker,
2376  "photonRelativeIso_endcaps",
2377  "photonRelativeIso for endcaps",
2378  100,
2379  0.0,
2380  2.,
2381  "photonRelativeIso_endcaps",
2382  "Events",
2383  "ELE_LOGY E1 P");
2384 
2385  // conversion rejection information
2386  h1_ele_convFlags = bookH1withSumw2(iBooker, "convFlags", "conversion rejection flag", 5, -2.5, 2.5);
2387  h1_ele_convFlags_all =
2388  bookH1withSumw2(iBooker, "convFlags_all", "conversion rejection flag, all electrons", 5, -2.5, 2.5);
2389  h1_ele_convDist = bookH1withSumw2(iBooker, "convDist", "distance to the conversion partner", 100, -15., 15.);
2390  h1_ele_convDist_all =
2391  bookH1withSumw2(iBooker, "convDist_all", "distance to the conversion partner, all electrons", 100, -15., 15.);
2392  h1_ele_convDcot = bookH1withSumw2(
2393  iBooker, "convDcot", "difference of cot(angle) with the conversion partner", 100, -CLHEP::pi / 2., CLHEP::pi / 2.);
2394  h1_ele_convDcot_all = bookH1withSumw2(iBooker,
2395  "convDcot_all",
2396  "difference of cot(angle) with the conversion partner, all electrons",
2397  100,
2398  -CLHEP::pi / 2.,
2399  CLHEP::pi / 2.);
2400  h1_ele_convRadius = bookH1withSumw2(iBooker, "convRadius", "signed conversion radius", 100, 0., 130.);
2401  h1_ele_convRadius_all =
2402  bookH1withSumw2(iBooker, "convRadius_all", "signed conversion radius, all electrons", 100, 0., 130.);
2403 }
2404 
2405 ElectronMcFakeValidator::~ElectronMcFakeValidator() {}
2406 
2407 //=========================================================================
2408 // Main methods
2409 //=========================================================================
2410 
2411 void ElectronMcFakeValidator::analyze(const edm::Event &iEvent, const edm::EventSetup &iSetup) {
2412  // get reco electrons
2413  auto gsfElectrons = iEvent.getHandle(electronCollection_);
2414  auto gsfElectronCores = iEvent.getHandle(electronCoreCollection_);
2415  auto gsfElectronTracks = iEvent.getHandle(electronTrackCollection_);
2416  auto gsfElectronSeeds = iEvent.getHandle(electronSeedCollection_);
2417 
2418  auto isoFromDepsTk03Handle = iEvent.getHandle(isoFromDepsTk03Tag_);
2419  auto isoFromDepsTk04Handle = iEvent.getHandle(isoFromDepsTk04Tag_);
2420  auto isoFromDepsEcalFull03Handle = iEvent.getHandle(isoFromDepsEcalFull03Tag_);
2421  auto isoFromDepsEcalFull04Handle = iEvent.getHandle(isoFromDepsEcalFull04Tag_);
2422  auto isoFromDepsEcalReduced03Handle = iEvent.getHandle(isoFromDepsEcalReduced03Tag_);
2423  auto isoFromDepsEcalReduced04Handle = iEvent.getHandle(isoFromDepsEcalReduced04Tag_);
2424  auto isoFromDepsHcal03Handle = iEvent.getHandle(isoFromDepsHcal03Tag_);
2425  auto isoFromDepsHcal04Handle = iEvent.getHandle(isoFromDepsHcal04Tag_);
2426 
2427  /*edm::Handle<reco::VertexCollection> vertexCollectionHandle;
2428  iEvent.getByToken(offlineVerticesCollection_, vertexCollectionHandle);*/
2429  auto vertexCollectionHandle = iEvent.getHandle(offlineVerticesCollection_);
2430  if (!vertexCollectionHandle.isValid()) {
2431  edm::LogInfo("ElectronMcFakeValidator::analyze") << "vertexCollectionHandle KO";
2432  } else {
2433  edm::LogInfo("ElectronMcFakeValidator::analyze") << "vertexCollectionHandle OK";
2434  }
2435 
2436  // get gen jets
2437  auto genJets = iEvent.getHandle(matchingObjectCollection_);
2438 
2439  // get the beamspot from the Event:
2440  auto recoBeamSpotHandle = iEvent.getHandle(beamSpotTag_);
2441  const BeamSpot bs = *recoBeamSpotHandle;
2442 
2443  edm::LogInfo("ElectronMcFakeValidator::analyze")
2444  << "Treating event " << iEvent.id() << " with " << gsfElectrons.product()->size() << " electrons";
2445  h1_recEleNum_->Fill((*gsfElectrons).size());
2446  h1_recCoreNum_->Fill((*gsfElectronCores).size());
2447  h1_recTrackNum_->Fill((*gsfElectronTracks).size());
2448  h1_recSeedNum_->Fill((*gsfElectronSeeds).size());
2449  h1_recOfflineVertices_->Fill((*vertexCollectionHandle).size());
2450 
2451  // all rec electrons
2452  reco::GsfElectronCollection::const_iterator gsfIter;
2453  for (gsfIter = gsfElectrons->begin(); gsfIter != gsfElectrons->end(); gsfIter++) {
2454  // preselect electrons
2455  if (gsfIter->pt() > maxPt_ || std::abs(gsfIter->eta()) > maxAbsEta_) {
2456  continue;
2457  }
2458 
2459  h1_ele_EoverP_all->Fill(gsfIter->eSuperClusterOverP());
2460  h1_ele_EseedOP_all->Fill(gsfIter->eSeedClusterOverP());
2461  h1_ele_EoPout_all->Fill(gsfIter->eSeedClusterOverPout());
2462  h1_ele_EeleOPout_all->Fill(gsfIter->eEleClusterOverPout());
2463  h1_ele_dEtaSc_propVtx_all->Fill(gsfIter->deltaEtaSuperClusterTrackAtVtx());
2464  h1_ele_dPhiSc_propVtx_all->Fill(gsfIter->deltaPhiSuperClusterTrackAtVtx());
2465  h1_ele_dEtaCl_propOut_all->Fill(gsfIter->deltaEtaSeedClusterTrackAtCalo());
2466  h1_ele_dPhiCl_propOut_all->Fill(gsfIter->deltaPhiSeedClusterTrackAtCalo());
2467  h1_ele_HoE_all->Fill(gsfIter->hadronicOverEm());
2468  h1_ele_HoE_bc_all->Fill(gsfIter->hcalOverEcalBc());
2469  double d = gsfIter->vertex().x() * gsfIter->vertex().x() + gsfIter->vertex().y() * gsfIter->vertex().y();
2470  h1_ele_TIP_all->Fill(sqrt(d));
2471  h1_ele_vertexEta_all->Fill(gsfIter->eta());
2472  h1_ele_vertexPt_all->Fill(gsfIter->pt());
2473  float enrj1 = gsfIter->ecalEnergy();
2474 
2475  // mee
2476  reco::GsfElectronCollection::const_iterator gsfIter2;
2477  for (gsfIter2 = gsfIter + 1; gsfIter2 != gsfElectrons->end(); gsfIter2++) {
2478  math::XYZTLorentzVector p12 = (*gsfIter).p4() + (*gsfIter2).p4();
2479  float mee2 = p12.Dot(p12);
2480  h1_ele_mee_all->Fill(sqrt(mee2));
2481  float enrj2 = gsfIter2->ecalEnergy();
2482  h2_ele_E2mnE1vsMee_all->Fill(sqrt(mee2), enrj2 - enrj1);
2483  if (gsfIter->ecalDrivenSeed() && gsfIter2->ecalDrivenSeed()) {
2484  h2_ele_E2mnE1vsMee_egeg_all->Fill(sqrt(mee2), enrj2 - enrj1);
2485  }
2486  if (gsfIter->charge() * gsfIter2->charge() < 0.) {
2487  h1_ele_mee_os->Fill(sqrt(mee2));
2488  }
2489  }
2490 
2491  // conversion rejection
2492  int flags = gsfIter->convFlags();
2493  if (flags == -9999) {
2494  flags = -1;
2495  }
2496  h1_ele_convFlags_all->Fill(flags);
2497  if (flags >= 0.) {
2498  h1_ele_convDist_all->Fill(gsfIter->convDist());
2499  h1_ele_convDcot_all->Fill(gsfIter->convDcot());
2500  h1_ele_convRadius_all->Fill(gsfIter->convRadius());
2501  }
2502  }
2503 
2504  // association matching object-reco electrons
2505  int matchingObjectNum = 0;
2506  reco::GenJetCollection::const_iterator moIter;
2507  for (moIter = genJets->begin(); moIter != genJets->end(); ++moIter) {
2508  // number of matching objects
2509  matchingObjectNum++;
2510 
2511  if (moIter->energy() / cosh(moIter->eta()) > maxPt_ || std::abs(moIter->eta()) > maxAbsEta_) {
2512  continue;
2513  }
2514 
2515  // suppress the endcaps
2516  //if (std::abs(moIter->eta()) > 1.5) continue;
2517  // select central z
2518  //if ( std::abs((*mcIter)->production_vertex()->position().z())>50.) continue;
2519 
2520  h1_matchingObjectEta->Fill(moIter->eta());
2521  h1_matchingObjectAbsEta->Fill(std::abs(moIter->eta()));
2522  h1_matchingObjectP->Fill(moIter->energy());
2523  h1_matchingObjectPt->Fill(moIter->energy() / cosh(moIter->eta()));
2524  h1_matchingObjectPhi->Fill(moIter->phi());
2525  h1_matchingObjectZ->Fill(moIter->vz());
2526 
2527  // looking for the best matching gsf electron
2528  bool okGsfFound = false;
2529  double gsfOkRatio = 999999.;
2530 
2531  // find best matched electron
2532  reco::GsfElectron bestGsfElectron;
2533  reco::GsfElectronRef bestGsfElectronRef;
2534  reco::GsfElectronCollection::const_iterator gsfIter;
2535  reco::GsfElectronCollection::size_type iElectron;
2536  for (gsfIter = gsfElectrons->begin(), iElectron = 0; gsfIter != gsfElectrons->end(); gsfIter++, iElectron++) {
2537  double dphi = gsfIter->phi() - moIter->phi();
2538  if (std::abs(dphi) > CLHEP::pi) {
2539  dphi = dphi < 0 ? (CLHEP::twopi) + dphi : dphi - CLHEP::twopi;
2540  }
2541  double deltaR = sqrt(pow((gsfIter->eta() - moIter->eta()), 2) + pow(dphi, 2));
2542  if (deltaR < deltaR_) {
2543  //if ( (genPc->pdg_id() == 11) && (gsfIter->charge() < 0.) || (genPc->pdg_id() == -11) &&
2544  //(gsfIter->charge() > 0.) ){
2545  double tmpGsfRatio = gsfIter->p() / moIter->energy();
2546  if (std::abs(tmpGsfRatio - 1) < std::abs(gsfOkRatio - 1)) {
2547  gsfOkRatio = tmpGsfRatio;
2548  bestGsfElectronRef = reco::GsfElectronRef(gsfElectrons, iElectron);
2549  bestGsfElectron = *gsfIter;
2550  okGsfFound = true;
2551  }
2552  //}
2553  }
2554  } // loop over rec ele to look for the best one
2555 
2556  // analysis when the matching object is matched by a rec electron
2557  if (okGsfFound) {
2558  // electron related distributions
2559  h1_ele_charge->Fill(bestGsfElectron.charge());
2560  h2_ele_chargeVsEta->Fill(bestGsfElectron.eta(), bestGsfElectron.charge());
2561  h2_ele_chargeVsPhi->Fill(bestGsfElectron.phi(), bestGsfElectron.charge());
2562  h2_ele_chargeVsPt->Fill(bestGsfElectron.pt(), bestGsfElectron.charge());
2563  h1_ele_vertexP->Fill(bestGsfElectron.p());
2564  h1_ele_vertexPt->Fill(bestGsfElectron.pt());
2565  h2_ele_vertexPtVsEta->Fill(bestGsfElectron.eta(), bestGsfElectron.pt());
2566  h2_ele_vertexPtVsPhi->Fill(bestGsfElectron.phi(), bestGsfElectron.pt());
2567  h1_ele_vertexEta->Fill(bestGsfElectron.eta());
2568  // generated distributions for matched electrons
2569  h1_ele_matchingObjectPt_matched->Fill(moIter->energy() / cosh(moIter->eta()));
2570  h1_ele_matchingObjectPhi_matched->Fill(moIter->phi());
2571  h1_ele_matchingObjectAbsEta_matched->Fill(std::abs(moIter->eta()));
2572  h1_ele_matchingObjectEta_matched->Fill(moIter->eta());
2573  h2_ele_vertexEtaVsPhi->Fill(bestGsfElectron.phi(), bestGsfElectron.eta());
2574  h1_ele_vertexPhi->Fill(bestGsfElectron.phi());
2575  h1_ele_vertexX->Fill(bestGsfElectron.vertex().x());
2576  h1_ele_vertexY->Fill(bestGsfElectron.vertex().y());
2577  h1_ele_vertexZ->Fill(bestGsfElectron.vertex().z());
2578  h1_ele_matchingObjectZ_matched->Fill(moIter->vz());
2579  double d =
2580  (bestGsfElectron.vertex().x() - bs.position().x()) * (bestGsfElectron.vertex().x() - bs.position().x()) +
2581  (bestGsfElectron.vertex().y() - bs.position().y()) * (bestGsfElectron.vertex().y() - bs.position().y());
2582  d = sqrt(d);
2583  h1_ele_vertexTIP->Fill(d);
2584  h2_ele_vertexTIPVsEta->Fill(bestGsfElectron.eta(), d);
2585  h2_ele_vertexTIPVsPhi->Fill(bestGsfElectron.phi(), d);
2586  h2_ele_vertexTIPVsPt->Fill(bestGsfElectron.pt(), d);
2587  h1_ele_EtaMnEtamatchingObject->Fill(bestGsfElectron.eta() - moIter->eta());
2588  h2_ele_EtaMnEtamatchingObjectVsEta->Fill(bestGsfElectron.eta(), bestGsfElectron.eta() - moIter->eta());
2589  h2_ele_EtaMnEtamatchingObjectVsPhi->Fill(bestGsfElectron.phi(), bestGsfElectron.eta() - moIter->eta());
2590  h2_ele_EtaMnEtamatchingObjectVsPt->Fill(bestGsfElectron.pt(), bestGsfElectron.eta() - moIter->eta());
2591  h1_ele_PhiMnPhimatchingObject->Fill(bestGsfElectron.phi() - moIter->phi());
2592  h1_ele_PhiMnPhimatchingObject2->Fill(bestGsfElectron.phi() - moIter->phi());
2593  h2_ele_PhiMnPhimatchingObjectVsEta->Fill(bestGsfElectron.eta(), bestGsfElectron.phi() - moIter->phi());
2594  h2_ele_PhiMnPhimatchingObjectVsPhi->Fill(bestGsfElectron.phi(), bestGsfElectron.phi() - moIter->phi());
2595  h2_ele_PhiMnPhimatchingObjectVsPt->Fill(bestGsfElectron.pt(), bestGsfElectron.phi() - moIter->phi());
2596  h1_ele_PoPmatchingObject->Fill(bestGsfElectron.p() / moIter->energy());
2597  h2_ele_PoPmatchingObjectVsEta->Fill(bestGsfElectron.eta(), bestGsfElectron.p() / moIter->energy());
2598  h2_ele_PoPmatchingObjectVsPhi->Fill(bestGsfElectron.phi(), bestGsfElectron.p() / moIter->energy());
2599  h2_ele_PoPmatchingObjectVsPt->Fill(bestGsfElectron.py(), bestGsfElectron.p() / moIter->energy());
2600  if (bestGsfElectron.isEB())
2601  h1_ele_PoPmatchingObject_barrel->Fill(bestGsfElectron.p() / moIter->energy());
2602  if (bestGsfElectron.isEE())
2603  h1_ele_PoPmatchingObject_endcaps->Fill(bestGsfElectron.p() / moIter->energy());
2604 
2605  // supercluster related distributions
2606  reco::SuperClusterRef sclRef = bestGsfElectron.superCluster();
2607  if (!bestGsfElectron.ecalDrivenSeed() && bestGsfElectron.trackerDrivenSeed())
2608  sclRef = bestGsfElectron.parentSuperCluster();
2609  if (sclRef.isNonnull()) {
2610  h1_scl_En_->Fill(sclRef->energy());
2611  double R = TMath::Sqrt(sclRef->x() * sclRef->x() + sclRef->y() * sclRef->y() + sclRef->z() * sclRef->z());
2612  double Rt = TMath::Sqrt(sclRef->x() * sclRef->x() + sclRef->y() * sclRef->y());
2613  h1_scl_Et_->Fill(sclRef->energy() * (Rt / R));
2614  h2_scl_EtVsEta_->Fill(sclRef->eta(), sclRef->energy() * (Rt / R));
2615  h2_scl_EtVsPhi_->Fill(sclRef->phi(), sclRef->energy() * (Rt / R));
2616  if (bestGsfElectron.isEB())
2617  h1_scl_EoEmatchingObject_barrel->Fill(sclRef->energy() / moIter->energy());
2618  if (bestGsfElectron.isEE())
2619  h1_scl_EoEmatchingObject_endcaps->Fill(sclRef->energy() / moIter->energy());
2620  h1_scl_Eta_->Fill(sclRef->eta());
2621  h2_scl_EtaVsPhi_->Fill(sclRef->phi(), sclRef->eta());
2622  h1_scl_Phi_->Fill(sclRef->phi());
2623  /*New from 06 05 2016*/
2624  // h1_scl_ESFrac->Fill( sclRef->preshowerEnergy() / sclRef->rawEnergy() );
2625  if (bestGsfElectron.isEE())
2626  h1_scl_ESFrac_endcaps->Fill(sclRef->preshowerEnergy() / sclRef->rawEnergy());
2627  }
2628  h1_scl_SigIEtaIEta_->Fill(bestGsfElectron.scSigmaIEtaIEta());
2629  if (bestGsfElectron.isEB())
2630  h1_scl_SigIEtaIEta_barrel_->Fill(bestGsfElectron.scSigmaIEtaIEta());
2631  if (bestGsfElectron.isEE())
2632  h1_scl_SigIEtaIEta_endcaps_->Fill(bestGsfElectron.scSigmaIEtaIEta());
2633  h1_scl_full5x5_sigmaIetaIeta_->Fill(bestGsfElectron.full5x5_sigmaIetaIeta());
2634  if (bestGsfElectron.isEB())
2635  h1_scl_full5x5_sigmaIetaIeta_barrel_->Fill(bestGsfElectron.full5x5_sigmaIetaIeta());
2636  if (bestGsfElectron.isEE())
2637  h1_scl_full5x5_sigmaIetaIeta_endcaps_->Fill(bestGsfElectron.full5x5_sigmaIetaIeta());
2638  h1_scl_E1x5_->Fill(bestGsfElectron.scE1x5());
2639  if (bestGsfElectron.isEB())
2640  h1_scl_E1x5_barrel_->Fill(bestGsfElectron.scE1x5());
2641  if (bestGsfElectron.isEE())
2642  h1_scl_E1x5_endcaps_->Fill(bestGsfElectron.scE1x5());
2643  h1_scl_E2x5max_->Fill(bestGsfElectron.scE2x5Max());
2644  if (bestGsfElectron.isEB())
2645  h1_scl_E2x5max_barrel_->Fill(bestGsfElectron.scE2x5Max());
2646  if (bestGsfElectron.isEE())
2647  h1_scl_E2x5max_endcaps_->Fill(bestGsfElectron.scE2x5Max());
2648  h1_scl_E5x5_->Fill(bestGsfElectron.scE5x5());
2649  if (bestGsfElectron.isEB())
2650  h1_scl_E5x5_barrel_->Fill(bestGsfElectron.scE5x5());
2651  if (bestGsfElectron.isEE())
2652  h1_scl_E5x5_endcaps_->Fill(bestGsfElectron.scE5x5());
2653 
2654  // track related distributions
2655  h1_ele_ambiguousTracks->Fill(bestGsfElectron.ambiguousGsfTracksSize());
2656  h2_ele_ambiguousTracksVsEta->Fill(bestGsfElectron.eta(), bestGsfElectron.ambiguousGsfTracksSize());
2657  h2_ele_ambiguousTracksVsPhi->Fill(bestGsfElectron.phi(), bestGsfElectron.ambiguousGsfTracksSize());
2658  h2_ele_ambiguousTracksVsPt->Fill(bestGsfElectron.pt(), bestGsfElectron.ambiguousGsfTracksSize());
2659  if (!readAOD_) { // track extra does not exist in AOD
2660  h1_ele_foundHits->Fill(bestGsfElectron.gsfTrack()->numberOfValidHits());
2661  h2_ele_foundHitsVsEta->Fill(bestGsfElectron.eta(), bestGsfElectron.gsfTrack()->numberOfValidHits());
2662  h2_ele_foundHitsVsPhi->Fill(bestGsfElectron.phi(), bestGsfElectron.gsfTrack()->numberOfValidHits());
2663  h2_ele_foundHitsVsPt->Fill(bestGsfElectron.pt(), bestGsfElectron.gsfTrack()->numberOfValidHits());
2664  h1_ele_lostHits->Fill(bestGsfElectron.gsfTrack()->numberOfLostHits());
2665  h2_ele_lostHitsVsEta->Fill(bestGsfElectron.eta(), bestGsfElectron.gsfTrack()->numberOfLostHits());
2666  h2_ele_lostHitsVsPhi->Fill(bestGsfElectron.phi(), bestGsfElectron.gsfTrack()->numberOfLostHits());
2667  h2_ele_lostHitsVsPt->Fill(bestGsfElectron.pt(), bestGsfElectron.gsfTrack()->numberOfLostHits());
2668  h1_ele_chi2->Fill(bestGsfElectron.gsfTrack()->normalizedChi2());
2669  h2_ele_chi2VsEta->Fill(bestGsfElectron.eta(), bestGsfElectron.gsfTrack()->normalizedChi2());
2670  h2_ele_chi2VsPhi->Fill(bestGsfElectron.phi(), bestGsfElectron.gsfTrack()->normalizedChi2());
2671  h2_ele_chi2VsPt->Fill(bestGsfElectron.pt(), bestGsfElectron.gsfTrack()->normalizedChi2());
2672  }
2673  // from gsf track interface, hence using mean
2674  if (!readAOD_) { // track extra does not exist in AOD
2675  h1_ele_PinMnPout->Fill(bestGsfElectron.gsfTrack()->innerMomentum().R() -
2676  bestGsfElectron.gsfTrack()->outerMomentum().R());
2677  h1_ele_outerP->Fill(bestGsfElectron.gsfTrack()->outerMomentum().R());
2678  h1_ele_outerPt->Fill(bestGsfElectron.gsfTrack()->outerMomentum().Rho());
2679  }
2680  // from electron interface, hence using mode
2681  h1_ele_PinMnPout_mode->Fill(bestGsfElectron.trackMomentumAtVtx().R() - bestGsfElectron.trackMomentumOut().R());
2682  h2_ele_PinMnPoutVsEta_mode->Fill(
2683  bestGsfElectron.eta(), bestGsfElectron.trackMomentumAtVtx().R() - bestGsfElectron.trackMomentumOut().R());
2684  h2_ele_PinMnPoutVsPhi_mode->Fill(
2685  bestGsfElectron.phi(), bestGsfElectron.trackMomentumAtVtx().R() - bestGsfElectron.trackMomentumOut().R());
2686  h2_ele_PinMnPoutVsPt_mode->Fill(
2687  bestGsfElectron.pt(), bestGsfElectron.trackMomentumAtVtx().R() - bestGsfElectron.trackMomentumOut().R());
2688  h2_ele_PinMnPoutVsE_mode->Fill(bestGsfElectron.caloEnergy(),
2689  bestGsfElectron.trackMomentumAtVtx().R() - bestGsfElectron.trackMomentumOut().R());
2690  if (!readAOD_) // track extra does not exist in AOD
2691  h2_ele_PinMnPoutVsChi2_mode->Fill(
2692  bestGsfElectron.gsfTrack()->normalizedChi2(),
2693  bestGsfElectron.trackMomentumAtVtx().R() - bestGsfElectron.trackMomentumOut().R());
2694  h1_ele_outerP_mode->Fill(bestGsfElectron.trackMomentumOut().R());
2695  h2_ele_outerPVsEta_mode->Fill(bestGsfElectron.eta(), bestGsfElectron.trackMomentumOut().R());
2696  h1_ele_outerPt_mode->Fill(bestGsfElectron.trackMomentumOut().Rho());
2697  h2_ele_outerPtVsEta_mode->Fill(bestGsfElectron.eta(), bestGsfElectron.trackMomentumOut().Rho());
2698  h2_ele_outerPtVsPhi_mode->Fill(bestGsfElectron.phi(), bestGsfElectron.trackMomentumOut().Rho());
2699  h2_ele_outerPtVsPt_mode->Fill(bestGsfElectron.pt(), bestGsfElectron.trackMomentumOut().Rho());
2700 
2701  if (!readAOD_) { // track extra does not exist in AOD
2702  edm::RefToBase<TrajectorySeed> seed = bestGsfElectron.gsfTrack()->extra()->seedRef();
2703  ElectronSeedRef elseed = seed.castTo<ElectronSeedRef>();
2704  h1_ele_seed_subdet2_->Fill(elseed->subDet(1));
2705  h1_ele_seed_mask_->Fill(elseed->hitsMask());
2706  if (elseed->subDet(1) == 1) {
2707  h1_ele_seed_mask_bpix_->Fill(elseed->hitsMask());
2708  } else if (elseed->subDet(1) == 2) {
2709  h1_ele_seed_mask_fpix_->Fill(elseed->hitsMask());
2710  } else if (elseed->subDet(1) == 6) {
2711  h1_ele_seed_mask_tec_->Fill(elseed->hitsMask());
2712  }
2713 
2714  if (elseed->dPhiNeg(1) != std::numeric_limits<float>::infinity()) {
2715  h1_ele_seed_dphi2_->Fill(elseed->dPhiNeg(1));
2716  h2_ele_seed_dphi2VsEta_->Fill(bestGsfElectron.eta(), elseed->dPhiNeg(1));
2717  h2_ele_seed_dphi2VsPt_->Fill(bestGsfElectron.pt(), elseed->dPhiNeg(1));
2718  }
2719  if (elseed->dPhiPos(1) != std::numeric_limits<float>::infinity()) {
2720  h1_ele_seed_dphi2pos_->Fill(elseed->dPhiPos(1));
2721  h2_ele_seed_dphi2posVsEta_->Fill(bestGsfElectron.eta(), elseed->dPhiPos(1));
2722  h2_ele_seed_dphi2posVsPt_->Fill(bestGsfElectron.pt(), elseed->dPhiPos(1));
2723  }
2724  if (elseed->dRZNeg(1) != std::numeric_limits<float>::infinity()) {
2725  h1_ele_seed_drz2_->Fill(elseed->dRZNeg(1));
2726  h2_ele_seed_drz2VsEta_->Fill(bestGsfElectron.eta(), elseed->dRZNeg(1));
2727  h2_ele_seed_drz2VsPt_->Fill(bestGsfElectron.pt(), elseed->dRZNeg(1));
2728  }
2729  if (elseed->dRZPos(1) != std::numeric_limits<float>::infinity()) {
2730  h1_ele_seed_drz2pos_->Fill(elseed->dRZPos(1));
2731  h2_ele_seed_drz2posVsEta_->Fill(bestGsfElectron.eta(), elseed->dRZPos(1));
2732  h2_ele_seed_drz2posVsPt_->Fill(bestGsfElectron.pt(), elseed->dRZPos(1));
2733  }
2734  }
2735  // match distributions
2736  h1_ele_EoP->Fill(bestGsfElectron.eSuperClusterOverP());
2737  if (bestGsfElectron.isEB())
2738  h1_ele_EoP_barrel->Fill(bestGsfElectron.eSuperClusterOverP());
2739  if (bestGsfElectron.isEE())
2740  h1_ele_EoP_endcaps->Fill(bestGsfElectron.eSuperClusterOverP());
2741  h2_ele_EoPVsEta->Fill(bestGsfElectron.eta(), bestGsfElectron.eSuperClusterOverP());
2742  h2_ele_EoPVsPhi->Fill(bestGsfElectron.phi(), bestGsfElectron.eSuperClusterOverP());
2743  h2_ele_EoPVsE->Fill(bestGsfElectron.caloEnergy(), bestGsfElectron.eSuperClusterOverP());
2744  h1_ele_EseedOP->Fill(bestGsfElectron.eSeedClusterOverP());
2745  if (bestGsfElectron.isEB())
2746  h1_ele_EseedOP_barrel->Fill(bestGsfElectron.eSeedClusterOverP());
2747  if (bestGsfElectron.isEE())
2748  h1_ele_EseedOP_endcaps->Fill(bestGsfElectron.eSeedClusterOverP());
2749  h2_ele_EseedOPVsEta->Fill(bestGsfElectron.eta(), bestGsfElectron.eSeedClusterOverP());
2750  h2_ele_EseedOPVsPhi->Fill(bestGsfElectron.phi(), bestGsfElectron.eSeedClusterOverP());
2751  h2_ele_EseedOPVsE->Fill(bestGsfElectron.caloEnergy(), bestGsfElectron.eSeedClusterOverP());
2752  h1_ele_EoPout->Fill(bestGsfElectron.eSeedClusterOverPout());
2753  if (bestGsfElectron.isEB())
2754  h1_ele_EoPout_barrel->Fill(bestGsfElectron.eSeedClusterOverPout());
2755  if (bestGsfElectron.isEE())
2756  h1_ele_EoPout_endcaps->Fill(bestGsfElectron.eSeedClusterOverPout());
2757  h2_ele_EoPoutVsEta->Fill(bestGsfElectron.eta(), bestGsfElectron.eSeedClusterOverPout());
2758  h2_ele_EoPoutVsPhi->Fill(bestGsfElectron.phi(), bestGsfElectron.eSeedClusterOverPout());
2759  h2_ele_EoPoutVsE->Fill(bestGsfElectron.caloEnergy(), bestGsfElectron.eSeedClusterOverPout());
2760  h1_ele_EeleOPout->Fill(bestGsfElectron.eEleClusterOverPout());
2761  if (bestGsfElectron.isEB())
2762  h1_ele_EeleOPout_barrel->Fill(bestGsfElectron.eEleClusterOverPout());
2763  if (bestGsfElectron.isEE())
2764  h1_ele_EeleOPout_endcaps->Fill(bestGsfElectron.eEleClusterOverPout());
2765  h2_ele_EeleOPoutVsEta->Fill(bestGsfElectron.eta(), bestGsfElectron.eEleClusterOverPout());
2766  h2_ele_EeleOPoutVsPhi->Fill(bestGsfElectron.phi(), bestGsfElectron.eEleClusterOverPout());
2767  h2_ele_EeleOPoutVsE->Fill(bestGsfElectron.caloEnergy(), bestGsfElectron.eEleClusterOverPout());
2768  h1_ele_dEtaSc_propVtx->Fill(bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
2769  if (bestGsfElectron.isEB())
2770  h1_ele_dEtaSc_propVtx_barrel->Fill(bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
2771  if (bestGsfElectron.isEE())
2772  h1_ele_dEtaSc_propVtx_endcaps->Fill(bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
2773  h2_ele_dEtaScVsEta_propVtx->Fill(bestGsfElectron.eta(), bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
2774  h2_ele_dEtaScVsPhi_propVtx->Fill(bestGsfElectron.phi(), bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
2775  h2_ele_dEtaScVsPt_propVtx->Fill(bestGsfElectron.pt(), bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
2776  h1_ele_dPhiSc_propVtx->Fill(bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
2777  if (bestGsfElectron.isEB())
2778  h1_ele_dPhiSc_propVtx_barrel->Fill(bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
2779  if (bestGsfElectron.isEE())
2780  h1_ele_dPhiSc_propVtx_endcaps->Fill(bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
2781  h2_ele_dPhiScVsEta_propVtx->Fill(bestGsfElectron.eta(), bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
2782  h2_ele_dPhiScVsPhi_propVtx->Fill(bestGsfElectron.phi(), bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
2783  h2_ele_dPhiScVsPt_propVtx->Fill(bestGsfElectron.pt(), bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
2784  h1_ele_dEtaCl_propOut->Fill(bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
2785  if (bestGsfElectron.isEB())
2786  h1_ele_dEtaCl_propOut_barrel->Fill(bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
2787  if (bestGsfElectron.isEE())
2788  h1_ele_dEtaCl_propOut_endcaps->Fill(bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
2789  h2_ele_dEtaClVsEta_propOut->Fill(bestGsfElectron.eta(), bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
2790  h2_ele_dEtaClVsPhi_propOut->Fill(bestGsfElectron.phi(), bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
2791  h2_ele_dEtaClVsPt_propOut->Fill(bestGsfElectron.pt(), bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
2792  h1_ele_dPhiCl_propOut->Fill(bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
2793  if (bestGsfElectron.isEB())
2794  h1_ele_dPhiCl_propOut_barrel->Fill(bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
2795  if (bestGsfElectron.isEE())
2796  h1_ele_dPhiCl_propOut_endcaps->Fill(bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
2797  h2_ele_dPhiClVsEta_propOut->Fill(bestGsfElectron.eta(), bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
2798  h2_ele_dPhiClVsPhi_propOut->Fill(bestGsfElectron.phi(), bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
2799  h2_ele_dPhiClVsPt_propOut->Fill(bestGsfElectron.pt(), bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
2800  h1_ele_dEtaEleCl_propOut->Fill(bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
2801  if (bestGsfElectron.isEB())
2802  h1_ele_dEtaEleCl_propOut_barrel->Fill(bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
2803  if (bestGsfElectron.isEE())
2804  h1_ele_dEtaEleCl_propOut_endcaps->Fill(bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
2805  h2_ele_dEtaEleClVsEta_propOut->Fill(bestGsfElectron.eta(), bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
2806  h2_ele_dEtaEleClVsPhi_propOut->Fill(bestGsfElectron.phi(), bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
2807  h2_ele_dEtaEleClVsPt_propOut->Fill(bestGsfElectron.pt(), bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
2808  h1_ele_dPhiEleCl_propOut->Fill(bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
2809  if (bestGsfElectron.isEB())
2810  h1_ele_dPhiEleCl_propOut_barrel->Fill(bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
2811  if (bestGsfElectron.isEE())
2812  h1_ele_dPhiEleCl_propOut_endcaps->Fill(bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
2813  h2_ele_dPhiEleClVsEta_propOut->Fill(bestGsfElectron.eta(), bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
2814  h2_ele_dPhiEleClVsPhi_propOut->Fill(bestGsfElectron.phi(), bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
2815  h2_ele_dPhiEleClVsPt_propOut->Fill(bestGsfElectron.pt(), bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
2816  h1_ele_HoE->Fill(bestGsfElectron.hadronicOverEm());
2817  h1_ele_HoE_bc->Fill(bestGsfElectron.hcalOverEcalBc());
2818  if (bestGsfElectron.isEB())
2819  h1_ele_HoE_bc_barrel->Fill(bestGsfElectron.hcalOverEcalBc());
2820  if (bestGsfElectron.isEE())
2821  h1_ele_HoE_bc_endcaps->Fill(bestGsfElectron.hcalOverEcalBc());
2822  if (bestGsfElectron.isEB())
2823  h1_ele_HoE_barrel->Fill(bestGsfElectron.hadronicOverEm());
2824  if (bestGsfElectron.isEE())
2825  h1_ele_HoE_endcaps->Fill(bestGsfElectron.hadronicOverEm());
2826  if (!bestGsfElectron.isEBEtaGap() && !bestGsfElectron.isEBPhiGap() && !bestGsfElectron.isEBEEGap() &&
2827  !bestGsfElectron.isEERingGap() && !bestGsfElectron.isEEDeeGap())
2828  h1_ele_HoE_fiducial->Fill(bestGsfElectron.hadronicOverEm());
2829  h2_ele_HoEVsEta->Fill(bestGsfElectron.eta(), bestGsfElectron.hadronicOverEm());
2830  h2_ele_HoEVsPhi->Fill(bestGsfElectron.phi(), bestGsfElectron.hadronicOverEm());
2831  h2_ele_HoEVsE->Fill(bestGsfElectron.caloEnergy(), bestGsfElectron.hadronicOverEm());
2832 
2833  //classes
2834  int eleClass = bestGsfElectron.classification();
2835  if (bestGsfElectron.isEE())
2836  eleClass += 10;
2837  h1_ele_classes->Fill(eleClass);
2838 
2839  h1_ele_eta->Fill(std::abs(bestGsfElectron.eta()));
2840  if (bestGsfElectron.classification() == GsfElectron::GOLDEN)
2841  h1_ele_eta_golden->Fill(std::abs(bestGsfElectron.eta()));
2842  if (bestGsfElectron.classification() == GsfElectron::BIGBREM)
2843  h1_ele_eta_bbrem->Fill(std::abs(bestGsfElectron.eta()));
2844  //if (bestGsfElectron.classification() == GsfElectron::OLDNARROW) h1_ele_eta_narrow->Fill(std::abs(bestGsfElectron.eta()));
2845  if (bestGsfElectron.classification() == GsfElectron::SHOWERING)
2846  h1_ele_eta_shower->Fill(std::abs(bestGsfElectron.eta()));
2847 
2848  // fbrem
2849 
2850  double fbrem_mode = bestGsfElectron.fbrem();
2851  h1_ele_fbrem->Fill(fbrem_mode);
2852  p1_ele_fbremVsEta_mode->Fill(bestGsfElectron.eta(), fbrem_mode);
2853 
2854  if (bestGsfElectron.isEB()) {
2855  double fbrem_mode_barrel = bestGsfElectron.fbrem();
2856  h1_ele_fbrem_barrel->Fill(fbrem_mode_barrel);
2857  }
2858 
2859  if (bestGsfElectron.isEE()) {
2860  double fbrem_mode_endcaps = bestGsfElectron.fbrem();
2861  h1_ele_fbrem_endcaps->Fill(fbrem_mode_endcaps);
2862  }
2863 
2864  // new 2014/02/12
2865  double superclusterfbrem_mode = bestGsfElectron.superClusterFbrem();
2866  h1_ele_superclusterfbrem->Fill(superclusterfbrem_mode);
2867 
2868  if (bestGsfElectron.isEB()) {
2869  double superclusterfbrem_mode_barrel = bestGsfElectron.superClusterFbrem();
2870  h1_ele_superclusterfbrem_barrel->Fill(superclusterfbrem_mode_barrel);
2871  }
2872 
2873  if (bestGsfElectron.isEE()) {
2874  double superclusterfbrem_mode_endcaps = bestGsfElectron.superClusterFbrem();
2875  h1_ele_superclusterfbrem_endcaps->Fill(superclusterfbrem_mode_endcaps);
2876  }
2877 
2878  if (!readAOD_) // track extra does not exist in AOD
2879  {
2880  double fbrem_mean =
2881  1. - bestGsfElectron.gsfTrack()->outerMomentum().R() / bestGsfElectron.gsfTrack()->innerMomentum().R();
2882  p1_ele_fbremVsEta_mean->Fill(bestGsfElectron.eta(), fbrem_mean);
2883  }
2884 
2885  //
2886 
2887  if (bestGsfElectron.classification() == GsfElectron::GOLDEN)
2888  h2_ele_PinVsPoutGolden_mode->Fill(bestGsfElectron.trackMomentumOut().R(),
2889  bestGsfElectron.trackMomentumAtVtx().R());
2890  if (bestGsfElectron.classification() == GsfElectron::SHOWERING)
2891  h2_ele_PinVsPoutShowering_mode->Fill(bestGsfElectron.trackMomentumOut().R(),
2892  bestGsfElectron.trackMomentumAtVtx().R());
2893  if (!readAOD_) // track extra does not exist in AOD
2894  if (bestGsfElectron.classification() == GsfElectron::GOLDEN)
2895  h2_ele_PinVsPoutGolden_mean->Fill(bestGsfElectron.gsfTrack()->outerMomentum().R(),
2896  bestGsfElectron.gsfTrack()->innerMomentum().R());
2897  if (!readAOD_) // track extra does not exist in AOD
2898  if (bestGsfElectron.classification() == GsfElectron::SHOWERING)
2899  h2_ele_PinVsPoutShowering_mean->Fill(bestGsfElectron.gsfTrack()->outerMomentum().R(),
2900  bestGsfElectron.gsfTrack()->innerMomentum().R());
2901  if (bestGsfElectron.classification() == GsfElectron::GOLDEN)
2902  h2_ele_PtinVsPtoutGolden_mode->Fill(bestGsfElectron.trackMomentumOut().Rho(),
2903  bestGsfElectron.trackMomentumAtVtx().Rho());
2904  if (bestGsfElectron.classification() == GsfElectron::SHOWERING)
2905  h2_ele_PtinVsPtoutShowering_mode->Fill(bestGsfElectron.trackMomentumOut().Rho(),
2906  bestGsfElectron.trackMomentumAtVtx().Rho());
2907  if (!readAOD_) // track extra does not exist in AOD
2908  if (bestGsfElectron.classification() == GsfElectron::GOLDEN)
2909  h2_ele_PtinVsPtoutGolden_mean->Fill(bestGsfElectron.gsfTrack()->outerMomentum().Rho(),
2910  bestGsfElectron.gsfTrack()->innerMomentum().Rho());
2911  if (!readAOD_) // track extra does not exist in AOD
2912  if (bestGsfElectron.classification() == GsfElectron::SHOWERING)
2913  h2_ele_PtinVsPtoutShowering_mean->Fill(bestGsfElectron.gsfTrack()->outerMomentum().Rho(),
2914  bestGsfElectron.gsfTrack()->innerMomentum().Rho());
2915 
2916  h1_ele_mva->Fill(bestGsfElectron.mva_e_pi());
2917  if (bestGsfElectron.isEB())
2918  h1_ele_mva_barrel->Fill(bestGsfElectron.mva_e_pi());
2919  if (bestGsfElectron.isEE())
2920  h1_ele_mva_endcaps->Fill(bestGsfElectron.mva_e_pi());
2921  h1_ele_mva_isolated->Fill(bestGsfElectron.mva_Isolated());
2922  if (bestGsfElectron.isEB())
2923  h1_ele_mva_barrel_isolated->Fill(bestGsfElectron.mva_Isolated());
2924  if (bestGsfElectron.isEE())
2925  h1_ele_mva_endcaps_isolated->Fill(bestGsfElectron.mva_Isolated());
2926  if (bestGsfElectron.ecalDrivenSeed())
2927  h1_ele_provenance->Fill(1.);
2928  if (bestGsfElectron.trackerDrivenSeed())
2929  h1_ele_provenance->Fill(-1.);
2930  if (bestGsfElectron.trackerDrivenSeed() || bestGsfElectron.ecalDrivenSeed())
2931  h1_ele_provenance->Fill(0.);
2932  if (bestGsfElectron.trackerDrivenSeed() && !bestGsfElectron.ecalDrivenSeed())
2933  h1_ele_provenance->Fill(-2.);
2934  if (!bestGsfElectron.trackerDrivenSeed() && bestGsfElectron.ecalDrivenSeed())
2935  h1_ele_provenance->Fill(2.);
2936  if (bestGsfElectron.ecalDrivenSeed() && bestGsfElectron.isEB())
2937  h1_ele_provenance_barrel->Fill(1.);
2938  if (bestGsfElectron.trackerDrivenSeed() && bestGsfElectron.isEB())
2939  h1_ele_provenance_barrel->Fill(-1.);
2940  if ((bestGsfElectron.trackerDrivenSeed() || bestGsfElectron.ecalDrivenSeed()) && bestGsfElectron.isEB())
2941  h1_ele_provenance_barrel->Fill(0.);
2942  if (bestGsfElectron.trackerDrivenSeed() && !bestGsfElectron.ecalDrivenSeed() && bestGsfElectron.isEB())
2943  h1_ele_provenance_barrel->Fill(-2.);
2944  if (!bestGsfElectron.trackerDrivenSeed() && bestGsfElectron.ecalDrivenSeed() && bestGsfElectron.isEB())
2945  h1_ele_provenance_barrel->Fill(2.);
2946  if (bestGsfElectron.ecalDrivenSeed() && bestGsfElectron.isEE())
2947  h1_ele_provenance_endcaps->Fill(1.);
2948  if (bestGsfElectron.trackerDrivenSeed() && bestGsfElectron.isEE())
2949  h1_ele_provenance_endcaps->Fill(-1.);
2950  if ((bestGsfElectron.trackerDrivenSeed() || bestGsfElectron.ecalDrivenSeed()) && bestGsfElectron.isEE())
2951  h1_ele_provenance_endcaps->Fill(0.);
2952  if (bestGsfElectron.trackerDrivenSeed() && !bestGsfElectron.ecalDrivenSeed() && bestGsfElectron.isEE())
2953  h1_ele_provenance_endcaps->Fill(-2.);
2954  if (!bestGsfElectron.trackerDrivenSeed() && bestGsfElectron.ecalDrivenSeed() && bestGsfElectron.isEE())
2955  h1_ele_provenance_endcaps->Fill(2.);
2956 
2957  // Pflow isolation
2958  h1_ele_chargedHadronIso->Fill(bestGsfElectron.pfIsolationVariables().sumChargedHadronPt);
2959  if (bestGsfElectron.isEB())
2960  h1_ele_chargedHadronIso_barrel->Fill(bestGsfElectron.pfIsolationVariables().sumChargedHadronPt);
2961  if (bestGsfElectron.isEE())
2962  h1_ele_chargedHadronIso_endcaps->Fill(bestGsfElectron.pfIsolationVariables().sumChargedHadronPt);
2963 
2964  h1_ele_neutralHadronIso->Fill(bestGsfElectron.pfIsolationVariables().sumNeutralHadronEt);
2965  if (bestGsfElectron.isEB())
2966  h1_ele_neutralHadronIso_barrel->Fill(bestGsfElectron.pfIsolationVariables().sumNeutralHadronEt);
2967  if (bestGsfElectron.isEE())
2968  h1_ele_neutralHadronIso_endcaps->Fill(bestGsfElectron.pfIsolationVariables().sumNeutralHadronEt);
2969 
2970  h1_ele_photonIso->Fill(bestGsfElectron.pfIsolationVariables().sumPhotonEt);
2971  if (bestGsfElectron.isEB())
2972  h1_ele_photonIso_barrel->Fill(bestGsfElectron.pfIsolationVariables().sumPhotonEt);
2973  if (bestGsfElectron.isEE())
2974  h1_ele_photonIso_endcaps->Fill(bestGsfElectron.pfIsolationVariables().sumPhotonEt);
2975 
2976  // -- pflow over pT
2977  h1_ele_chargedHadronRelativeIso->Fill(bestGsfElectron.pfIsolationVariables().sumChargedHadronPt /
2978  bestGsfElectron.pt());
2979  if (bestGsfElectron.isEB())
2980  h1_ele_chargedHadronRelativeIso_barrel->Fill(bestGsfElectron.pfIsolationVariables().sumChargedHadronPt /
2981  bestGsfElectron.pt());
2982  if (bestGsfElectron.isEE())
2983  h1_ele_chargedHadronRelativeIso_endcaps->Fill(bestGsfElectron.pfIsolationVariables().sumChargedHadronPt /
2984  bestGsfElectron.pt());
2985 
2986  h1_ele_neutralHadronRelativeIso->Fill(bestGsfElectron.pfIsolationVariables().sumNeutralHadronEt /
2987  bestGsfElectron.pt());
2988  if (bestGsfElectron.isEB())
2989  h1_ele_neutralHadronRelativeIso_barrel->Fill(bestGsfElectron.pfIsolationVariables().sumNeutralHadronEt /
2990  bestGsfElectron.pt());
2991  if (bestGsfElectron.isEE())
2992  h1_ele_neutralHadronRelativeIso_endcaps->Fill(bestGsfElectron.pfIsolationVariables().sumNeutralHadronEt /
2993  bestGsfElectron.pt());
2994 
2995  h1_ele_photonRelativeIso->Fill(bestGsfElectron.pfIsolationVariables().sumPhotonEt / bestGsfElectron.pt());
2996  if (bestGsfElectron.isEB())
2997  h1_ele_photonRelativeIso_barrel->Fill(bestGsfElectron.pfIsolationVariables().sumPhotonEt /
2998  bestGsfElectron.pt());
2999  if (bestGsfElectron.isEE())
3000  h1_ele_photonRelativeIso_endcaps->Fill(bestGsfElectron.pfIsolationVariables().sumPhotonEt /
3001  bestGsfElectron.pt());
3002 
3003  // isolation
3004  h1_ele_tkSumPt_dr03->Fill(bestGsfElectron.dr03TkSumPt());
3005  if (bestGsfElectron.isEB())
3006  h1_ele_tkSumPt_dr03_barrel->Fill(bestGsfElectron.dr03TkSumPt());
3007  if (bestGsfElectron.isEE())
3008  h1_ele_tkSumPt_dr03_endcaps->Fill(bestGsfElectron.dr03TkSumPt());
3009  h1_ele_ecalRecHitSumEt_dr03->Fill(bestGsfElectron.dr03EcalRecHitSumEt());
3010  if (bestGsfElectron.isEB())
3011  h1_ele_ecalRecHitSumEt_dr03_barrel->Fill(bestGsfElectron.dr03EcalRecHitSumEt());
3012  if (bestGsfElectron.isEE())
3013  h1_ele_ecalRecHitSumEt_dr03_endcaps->Fill(bestGsfElectron.dr03EcalRecHitSumEt());
3014  h1_ele_hcalTowerSumEt_dr03_depth1->Fill(bestGsfElectron.dr03HcalTowerSumEt(1));
3015  if (bestGsfElectron.isEB())
3016  h1_ele_hcalTowerSumEt_dr03_depth1_barrel->Fill(bestGsfElectron.dr03HcalTowerSumEt(1));
3017  if (bestGsfElectron.isEE())
3018  h1_ele_hcalTowerSumEt_dr03_depth1_endcaps->Fill(bestGsfElectron.dr03HcalTowerSumEt(1));
3019  h1_ele_hcalTowerSumEt_dr03_depth2->Fill(bestGsfElectron.dr03HcalTowerSumEt(2));
3020  h1_ele_hcalTowerSumEtBc_dr03_depth1->Fill(bestGsfElectron.dr03HcalTowerSumEtBc(1));
3021  if (bestGsfElectron.isEB())
3022  h1_ele_hcalTowerSumEtBc_dr03_depth1_barrel->Fill(bestGsfElectron.dr03HcalTowerSumEtBc(1));
3023  if (bestGsfElectron.isEE())
3024  h1_ele_hcalTowerSumEtBc_dr03_depth1_endcaps->Fill(bestGsfElectron.dr03HcalTowerSumEtBc(1));
3025  h1_ele_hcalTowerSumEtBc_dr03_depth2->Fill(bestGsfElectron.dr03HcalTowerSumEtBc(2));
3026  if (bestGsfElectron.isEB())
3027  h1_ele_hcalTowerSumEtBc_dr03_depth2_barrel->Fill(bestGsfElectron.dr03HcalTowerSumEtBc(2));
3028  if (bestGsfElectron.isEE())
3029  h1_ele_hcalTowerSumEtBc_dr03_depth2_endcaps->Fill(bestGsfElectron.dr03HcalTowerSumEtBc(2));
3030  h1_ele_tkSumPt_dr04->Fill(bestGsfElectron.dr04TkSumPt());
3031  if (bestGsfElectron.isEB())
3032  h1_ele_tkSumPt_dr04_barrel->Fill(bestGsfElectron.dr04TkSumPt());
3033  if (bestGsfElectron.isEE())
3034  h1_ele_tkSumPt_dr04_endcaps->Fill(bestGsfElectron.dr04TkSumPt());
3035  h1_ele_ecalRecHitSumEt_dr04->Fill(bestGsfElectron.dr04EcalRecHitSumEt());
3036  if (bestGsfElectron.isEB())
3037  h1_ele_ecalRecHitSumEt_dr04_barrel->Fill(bestGsfElectron.dr04EcalRecHitSumEt());
3038  if (bestGsfElectron.isEE())
3039  h1_ele_ecalRecHitSumEt_dr04_endcaps->Fill(bestGsfElectron.dr04EcalRecHitSumEt());
3040  h1_ele_hcalTowerSumEt_dr04_depth1->Fill(bestGsfElectron.dr04HcalTowerSumEt(1));
3041  if (bestGsfElectron.isEB())
3042  h1_ele_hcalTowerSumEt_dr04_depth1_barrel->Fill(bestGsfElectron.dr04HcalTowerSumEt(1));
3043  if (bestGsfElectron.isEE())
3044  h1_ele_hcalTowerSumEt_dr04_depth1_endcaps->Fill(bestGsfElectron.dr04HcalTowerSumEt(1));
3045  h1_ele_hcalTowerSumEt_dr04_depth2->Fill(bestGsfElectron.dr04HcalTowerSumEt(2));
3046  h1_ele_hcalTowerSumEtBc_dr04_depth1->Fill(bestGsfElectron.dr04HcalTowerSumEtBc(1));
3047  if (bestGsfElectron.isEB())
3048  h1_ele_hcalTowerSumEtBc_dr04_depth1_barrel->Fill(bestGsfElectron.dr04HcalTowerSumEtBc(1));
3049  if (bestGsfElectron.isEE())
3050  h1_ele_hcalTowerSumEtBc_dr04_depth1_endcaps->Fill(bestGsfElectron.dr04HcalTowerSumEtBc(1));
3051  h1_ele_hcalTowerSumEtBc_dr04_depth2->Fill(bestGsfElectron.dr04HcalTowerSumEtBc(2));
3052  if (bestGsfElectron.isEB())
3053  h1_ele_hcalTowerSumEtBc_dr04_depth2_barrel->Fill(bestGsfElectron.dr04HcalTowerSumEtBc(2));
3054  if (bestGsfElectron.isEE())
3055  h1_ele_hcalTowerSumEtBc_dr04_depth2_endcaps->Fill(bestGsfElectron.dr04HcalTowerSumEtBc(2));
3056 
3057  h1_ele_hcalDepth1OverEcalBc->Fill(bestGsfElectron.hcalOverEcalBc(1));
3058  if (bestGsfElectron.isEB())
3059  h1_ele_hcalDepth1OverEcalBc_barrel->Fill(bestGsfElectron.hcalOverEcalBc(1));
3060  if (bestGsfElectron.isEE())
3061  h1_ele_hcalDepth1OverEcalBc_endcaps->Fill(bestGsfElectron.hcalOverEcalBc(1));
3062  h1_ele_hcalDepth2OverEcalBc->Fill(bestGsfElectron.hcalOverEcalBc(2));
3063  if (bestGsfElectron.isEB())
3064  h1_ele_hcalDepth2OverEcalBc_barrel->Fill(bestGsfElectron.hcalOverEcalBc(2));
3065  if (bestGsfElectron.isEE())
3066  h1_ele_hcalDepth2OverEcalBc_endcaps->Fill(bestGsfElectron.hcalOverEcalBc(2));
3067 
3068  // conversion rejection
3069  int flags = bestGsfElectron.convFlags();
3070  if (flags == -9999) {
3071  flags = -1;
3072  }
3073  h1_ele_convFlags->Fill(flags);
3074  if (flags >= 0.) {
3075  h1_ele_convDist->Fill(bestGsfElectron.convDist());
3076  h1_ele_convDcot->Fill(bestGsfElectron.convDcot());
3077  h1_ele_convRadius->Fill(bestGsfElectron.convRadius());
3078  }
3079 
3080  } // gsf electron found
3081 
3082  } // loop overmatching object
3083 
3084  h1_matchingObjectNum->Fill(matchingObjectNum);
3085 }
reco::GsfElectron::deltaEtaSeedClusterTrackAtCalo
float deltaEtaSeedClusterTrackAtCalo() const
Definition: GsfElectron.h:234
ElectronMcFakeValidator::phi2D_nbin
int phi2D_nbin
Definition: ElectronMcFakeValidator.h:102
reco::GsfElectron::convFlags
int convFlags() const
Definition: GsfElectron.h:651
ElectronMcFakeValidator::h1_ele_convDcot_all
MonitorElement * h1_ele_convDcot_all
Definition: ElectronMcFakeValidator.h:486
ElectronMcFakeValidator::h1_ele_matchingObjectPt_matched
MonitorElement * h1_ele_matchingObjectPt_matched
Definition: ElectronMcFakeValidator.h:181
ElectronMcFakeValidator::h1_ele_eta
MonitorElement * h1_ele_eta
Definition: ElectronMcFakeValidator.h:368
ElectronMcFakeValidator::dphimatch_min
double dphimatch_min
Definition: ElectronMcFakeValidator.h:110
reco::GsfElectron::convDcot
float convDcot() const
Definition: GsfElectron.h:654
ElectronMcFakeValidator::h1_scl_EoEmatchingObjectShowering_barrel
MonitorElement * h1_scl_EoEmatchingObjectShowering_barrel
Definition: ElectronMcFakeValidator.h:406
ElectronMcFakeValidator::h1_ele_chargedHadronIso
MonitorElement * h1_ele_chargedHadronIso
Definition: ElectronMcFakeValidator.h:420
reco::GsfElectron::isEE
bool isEE() const
Definition: GsfElectron.h:337
ElectronMcFakeValidator::ele_min
double ele_min
Definition: ElectronMcFakeValidator.h:131
reco::GsfElectron::dr04EcalRecHitSumEt
float dr04EcalRecHitSumEt() const
Definition: GsfElectron.h:607
ElectronMcFakeValidator::h1_ele_matchingObjectZ_matched
MonitorElement * h1_ele_matchingObjectZ_matched
Definition: ElectronMcFakeValidator.h:183
ElectronMcFakeValidator::h1_scl_Eta_
MonitorElement * h1_scl_Eta_
Definition: ElectronMcFakeValidator.h:229
GenJetCollection.h
L1CaloJetHTTProducer_cfi.genJets
genJets
Definition: L1CaloJetHTTProducer_cfi.py:7
ElectronMcFakeValidator::h1_ele_provenance
MonitorElement * h1_ele_provenance
Definition: ElectronMcFakeValidator.h:415
ElectronMcFakeValidator::h1_ele_TIP_all
MonitorElement * h1_ele_TIP_all
Definition: ElectronMcFakeValidator.h:168
ElectronMcFakeValidator::isoFromDepsHcal03Tag_
edm::EDGetTokenT< edm::ValueMap< double > > isoFromDepsHcal03Tag_
Definition: ElectronMcFakeValidator.h:57
ElectronMcFakeValidator::h2_ele_EtaMnEtamatchingObjectVsPt
MonitorElement * h2_ele_EtaMnEtamatchingObjectVsPt
Definition: ElectronMcFakeValidator.h:215
Handle.h
ElectronMcFakeValidator::h1_ele_hcalTowerSumEtBc_dr03_depth2_endcaps
MonitorElement * h1_ele_hcalTowerSumEtBc_dr03_depth2_endcaps
Definition: ElectronMcFakeValidator.h:473
ElectronMcFakeValidator::readAOD_
bool readAOD_
Definition: ElectronMcFakeValidator.h:49
ElectronMcFakeValidator::electronCoreCollection_
edm::EDGetTokenT< reco::GsfElectronCoreCollection > electronCoreCollection_
Definition: ElectronMcFakeValidator.h:43
ElectronDqmAnalyzerBase::setBookEfficiencyFlag
void setBookEfficiencyFlag(const bool &)
Definition: ElectronDqmAnalyzerBase.cc:34
ElectronMcFakeValidator::h1_ele_chargedHadronRelativeIso_barrel
MonitorElement * h1_ele_chargedHadronRelativeIso_barrel
Definition: ElectronMcFakeValidator.h:431
reco::GsfElectron::dr03TkSumPt
float dr03TkSumPt() const
Definition: GsfElectron.h:565
ElectronMcFakeValidator::eop2D_nbin
int eop2D_nbin
Definition: ElectronMcFakeValidator.h:113
ElectronMcFakeValidator::h1_scl_E2x5max_
MonitorElement * h1_scl_E2x5max_
Definition: ElectronMcFakeValidator.h:241
ElectronMcFakeValidator::h2_ele_seed_dphi2VsPt_
MonitorElement * h2_ele_seed_dphi2VsPt_
Definition: ElectronMcFakeValidator.h:356
reco::GsfElectron::isEB
bool isEB() const
Definition: GsfElectron.h:336
reco::GsfElectron::scE1x5
float scE1x5() const
Definition: GsfElectron.h:505
ElectronMcFakeValidator::h2_ele_dEtaClVsPhi_propOut
MonitorElement * h2_ele_dEtaClVsPhi_propOut
Definition: ElectronMcFakeValidator.h:328
reco::GsfElectron::gsfTrack
GsfTrackRef gsfTrack() const override
reference to a GsfTrack
Definition: GsfElectron.h:164
MessageLogger.h
reco::GsfElectron::trackMomentumOut
math::XYZVectorF trackMomentumOut() const
Definition: GsfElectron.h:278
ElectronMcFakeValidator::h1_ele_seed_drz2_
MonitorElement * h1_ele_seed_drz2_
Definition: ElectronMcFakeValidator.h:360
ElectronMcFakeValidator::h2_ele_dPhiEleClVsEta_propOut
MonitorElement * h2_ele_dPhiEleClVsEta_propOut
Definition: ElectronMcFakeValidator.h:345
ElectronMcFakeValidator::h1_ele_dEtaSc_propVtx_endcaps
MonitorElement * h1_ele_dEtaSc_propVtx_endcaps
Definition: ElectronMcFakeValidator.h:314
ElectronMcFakeValidator::h2_ele_ambiguousTracksVsPt
MonitorElement * h2_ele_ambiguousTracksVsPt
Definition: ElectronMcFakeValidator.h:251
ElectronMcFakeValidator::h1_ele_seed_mask_bpix_
MonitorElement * h1_ele_seed_mask_bpix_
Definition: ElectronMcFakeValidator.h:351
ElectronMcFakeValidator::h2_ele_vertexTIPVsPhi
MonitorElement * h2_ele_vertexTIPVsPhi
Definition: ElectronMcFakeValidator.h:202
ElectronMcFakeValidator::mee_nbin
int mee_nbin
Definition: ElectronMcFakeValidator.h:116
ElectronMcFakeValidator::h2_ele_outerPVsEta_mode
MonitorElement * h2_ele_outerPVsEta_mode
Definition: ElectronMcFakeValidator.h:281
ElectronMcFakeValidator::h1_ele_hcalTowerSumEt_dr04_depth1
MonitorElement * h1_ele_hcalTowerSumEt_dr04_depth1
Definition: ElectronMcFakeValidator.h:456
ElectronMcFakeValidator::h1_ele_chi2_endcaps
MonitorElement * h1_ele_chi2_endcaps
Definition: ElectronMcFakeValidator.h:266
ElectronMcFakeValidator::h2_ele_EoPoutVsE
MonitorElement * h2_ele_EoPoutVsE
Definition: ElectronMcFakeValidator.h:304
ElectronMcFakeValidator::h1_scl_Et_
MonitorElement * h1_scl_Et_
Definition: ElectronMcFakeValidator.h:225
ElectronMcFakeValidator::h1_ele_hcalDepth2OverEcalBc_endcaps
MonitorElement * h1_ele_hcalDepth2OverEcalBc_endcaps
Definition: ElectronMcFakeValidator.h:466
ElectronMcFakeValidator::h2_ele_HoEVsE
MonitorElement * h2_ele_HoEVsE
Definition: ElectronMcFakeValidator.h:383
ElectronMcFakeValidator::h1_matchingObjectAbsEta
MonitorElement * h1_matchingObjectAbsEta
Definition: ElectronMcFakeValidator.h:154
reco::GsfElectron::dr03HcalTowerSumEtBc
float dr03HcalTowerSumEtBc(int depth=0) const
Definition: GsfElectron.h:601
reco::GsfElectron::deltaPhiSuperClusterTrackAtVtx
float deltaPhiSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:236
ElectronMcFakeValidator::h1_ele_HoE_bc_endcaps
MonitorElement * h1_ele_HoE_bc_endcaps
Definition: ElectronMcFakeValidator.h:379
ElectronMcFakeValidator::h1_ele_hcalTowerSumEtBc_dr03_depth2
MonitorElement * h1_ele_hcalTowerSumEtBc_dr03_depth2
Definition: ElectronMcFakeValidator.h:471
ElectronMcFakeValidator::h1_ele_mva_endcaps_isolated
MonitorElement * h1_ele_mva_endcaps_isolated
Definition: ElectronMcFakeValidator.h:414
ElectronMcFakeValidator::h1_ele_matchingObjectAbsEta_matched
MonitorElement * h1_ele_matchingObjectAbsEta_matched
Definition: ElectronMcFakeValidator.h:180
reco::GsfElectron::parentSuperCluster
virtual SuperClusterRef parentSuperCluster() const
Definition: GsfElectron.h:168
ElectronMcFakeValidator::h1_ele_hcalTowerSumEtBc_dr04_depth2_endcaps
MonitorElement * h1_ele_hcalTowerSumEtBc_dr04_depth2_endcaps
Definition: ElectronMcFakeValidator.h:479
ElectronMcFakeValidator::h1_ele_lostHits
MonitorElement * h1_ele_lostHits
Definition: ElectronMcFakeValidator.h:258
ElectronMcFakeValidator::h2_ele_EoPVsEta
MonitorElement * h2_ele_EoPVsEta
Definition: ElectronMcFakeValidator.h:290
ElectronMcFakeValidator::h1_ele_EseedOP_all
MonitorElement * h1_ele_EseedOP_all
Definition: ElectronMcFakeValidator.h:161
ElectronMcFakeValidator::h1_ele_hcalTowerSumEt_dr03_depth1_barrel
MonitorElement * h1_ele_hcalTowerSumEt_dr03_depth1_barrel
Definition: ElectronMcFakeValidator.h:447
ElectronMcFakeValidator::h2_ele_EtaMnEtamatchingObjectVsEta
MonitorElement * h2_ele_EtaMnEtamatchingObjectVsEta
Definition: ElectronMcFakeValidator.h:213
ElectronMcFakeValidator::set_EfficiencyFlag
bool set_EfficiencyFlag
Definition: ElectronMcFakeValidator.h:125
ElectronMcFakeValidator::dphimatch2D_nbin
int dphimatch2D_nbin
Definition: ElectronMcFakeValidator.h:109
ElectronMcFakeValidator::h1_ele_dPhiSc_propVtx
MonitorElement * h1_ele_dPhiSc_propVtx
Definition: ElectronMcFakeValidator.h:318
reco::GsfElectron::isEERingGap
bool isEERingGap() const
Definition: GsfElectron.h:345
edm::Run
Definition: Run.h:45
ElectronMcFakeValidator::h1_recEleNum_
MonitorElement * h1_recEleNum_
Definition: ElectronMcFakeValidator.h:146
ElectronDqmAnalyzerBase
Definition: ElectronDqmAnalyzerBase.h:19
ElectronMcFakeValidator::h1_ele_photonIso_barrel
MonitorElement * h1_ele_photonIso_barrel
Definition: ElectronMcFakeValidator.h:427
ElectronMcFakeValidator::h2_ele_foundHitsVsPhi
MonitorElement * h2_ele_foundHitsVsPhi
Definition: ElectronMcFakeValidator.h:256
ElectronMcFakeValidator::h1_matchingObjectP
MonitorElement * h1_matchingObjectP
Definition: ElectronMcFakeValidator.h:155
ElectronMcFakeValidator::h2_ele_lostHitsVsEta
MonitorElement * h2_ele_lostHitsVsEta
Definition: ElectronMcFakeValidator.h:261
ElectronMcFakeValidator::h1_ele_hcalTowerSumEtBc_dr04_depth1_endcaps
MonitorElement * h1_ele_hcalTowerSumEtBc_dr04_depth1_endcaps
Definition: ElectronMcFakeValidator.h:476
ElectronMcFakeValidator::phi_nbin
int phi_nbin
Definition: ElectronMcFakeValidator.h:101
ElectronMcFakeValidator::ElectronMcFakeValidator
ElectronMcFakeValidator(const edm::ParameterSet &conf)
Definition: ElectronMcFakeValidator.cc:45
ElectronMcFakeValidator::h1_matchingObjectNum
MonitorElement * h1_matchingObjectNum
Definition: ElectronMcFakeValidator.h:145
ElectronMcFakeValidator::h1_ele_provenance_barrel
MonitorElement * h1_ele_provenance_barrel
Definition: ElectronMcFakeValidator.h:416
ElectronMcFakeValidator::h2_ele_vertexPtVsEta
MonitorElement * h2_ele_vertexPtVsEta
Definition: ElectronMcFakeValidator.h:191
ElectronMcFakeValidator::h2_ele_EoPVsE
MonitorElement * h2_ele_EoPVsE
Definition: ElectronMcFakeValidator.h:292
ElectronMcFakeValidator::h2_ele_PhiMnPhimatchingObjectVsEta
MonitorElement * h2_ele_PhiMnPhimatchingObjectVsEta
Definition: ElectronMcFakeValidator.h:218
ElectronMcFakeValidator::h2_ele_vertexEtaVsPhi
MonitorElement * h2_ele_vertexEtaVsPhi
Definition: ElectronMcFakeValidator.h:194
ElectronMcFakeValidator::lhits_nbin
int lhits_nbin
Definition: ElectronMcFakeValidator.h:88
ElectronMcFakeValidator::h1_ele_HoE_bc_all
MonitorElement * h1_ele_HoE_bc_all
Definition: ElectronMcFakeValidator.h:170
ElectronMcFakeValidator::h1_ele_EoP_barrel
MonitorElement * h1_ele_EoP_barrel
Definition: ElectronMcFakeValidator.h:288
ElectronMcFakeValidator::h2_ele_chi2VsPhi
MonitorElement * h2_ele_chi2VsPhi
Definition: ElectronMcFakeValidator.h:268
ElectronMcFakeValidator::mee_min
double mee_min
Definition: ElectronMcFakeValidator.h:117
reco::GsfElectron::dr04HcalTowerSumEt
float dr04HcalTowerSumEt(int depth=0) const
Definition: GsfElectron.h:608
ElectronMcFakeValidator::h1_ele_outerPt
MonitorElement * h1_ele_outerPt
Definition: ElectronMcFakeValidator.h:282
reco::GsfElectron::scE2x5Max
float scE2x5Max() const
Definition: GsfElectron.h:506
ElectronMcFakeValidator::h1_scl_E2x5max_endcaps_
MonitorElement * h1_scl_E2x5max_endcaps_
Definition: ElectronMcFakeValidator.h:243
dqm::implementation::NavigatorBase::setCurrentFolder
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:32
ElectronMcFakeValidator::dphi_min
double dphi_min
Definition: ElectronMcFakeValidator.h:106
ElectronMcFakeValidator::maxAbsEta_
double maxAbsEta_
Definition: ElectronMcFakeValidator.h:68
ElectronMcFakeValidator::h1_ele_EoP_endcaps
MonitorElement * h1_ele_EoP_endcaps
Definition: ElectronMcFakeValidator.h:289
ElectronMcFakeValidator::h2_ele_dPhiEleClVsPt_propOut
MonitorElement * h2_ele_dPhiEleClVsPt_propOut
Definition: ElectronMcFakeValidator.h:347
ElectronMcFakeValidator::deta_min
double deta_min
Definition: ElectronMcFakeValidator.h:95
ElectronMcFakeValidator::xyz_nbin
int xyz_nbin
Definition: ElectronMcFakeValidator.h:78
ElectronMcFakeValidator::ele_max
double ele_max
Definition: ElectronMcFakeValidator.h:132
ElectronMcFakeValidator::h2_ele_dEtaScVsPhi_propVtx
MonitorElement * h2_ele_dEtaScVsPhi_propVtx
Definition: ElectronMcFakeValidator.h:316
ElectronMcFakeValidator::h2_ele_PhiMnPhimatchingObjectVsPhi
MonitorElement * h2_ele_PhiMnPhimatchingObjectVsPhi
Definition: ElectronMcFakeValidator.h:219
ElectronMcFakeValidator::h1_ele_PhiMnPhimatchingObject2
MonitorElement * h1_ele_PhiMnPhimatchingObject2
Definition: ElectronMcFakeValidator.h:217
ElectronMcFakeValidator::track_min
double track_min
Definition: ElectronMcFakeValidator.h:137
ElectronMcFakeValidator::h2_ele_PinMnPoutVsPhi_mode
MonitorElement * h2_ele_PinMnPoutVsPhi_mode
Definition: ElectronMcFakeValidator.h:274
ElectronMcFakeValidator::h1_ele_matchingObjectEta_matched
MonitorElement * h1_ele_matchingObjectEta_matched
Definition: ElectronMcFakeValidator.h:179
ElectronMcFakeValidator::h1_ele_neutralHadronIso_endcaps
MonitorElement * h1_ele_neutralHadronIso_endcaps
Definition: ElectronMcFakeValidator.h:425
ElectronMcFakeValidator::h1_ele_lostHits_endcaps
MonitorElement * h1_ele_lostHits_endcaps
Definition: ElectronMcFakeValidator.h:260
ElectronMcFakeValidator::h2_ele_E2mnE1vsMee_all
MonitorElement * h2_ele_E2mnE1vsMee_all
Definition: ElectronMcFakeValidator.h:176
ElectronMcFakeValidator::pteff_nbin
int pteff_nbin
Definition: ElectronMcFakeValidator.h:84
ElectronMcFakeValidator::h1_ele_dEtaCl_propOut_all
MonitorElement * h1_ele_dEtaCl_propOut_all
Definition: ElectronMcFakeValidator.h:166
ElectronMcFakeValidator::h1_ele_EoPout_barrel
MonitorElement * h1_ele_EoPout_barrel
Definition: ElectronMcFakeValidator.h:300
DQMStore.h
ElectronMcFakeValidator::h2_ele_dPhiEleClVsPhi_propOut
MonitorElement * h2_ele_dPhiEleClVsPhi_propOut
Definition: ElectronMcFakeValidator.h:346
ElectronMcFakeValidator::h2_ele_dEtaEleClVsPhi_propOut
MonitorElement * h2_ele_dEtaEleClVsPhi_propOut
Definition: ElectronMcFakeValidator.h:340
ElectronMcFakeValidator::h2_ele_PoPmatchingObjectVsPt
MonitorElement * h2_ele_PoPmatchingObjectVsPt
Definition: ElectronMcFakeValidator.h:208
ElectronMcFakeValidator::h1_ele_vertexEta
MonitorElement * h1_ele_vertexEta
Definition: ElectronMcFakeValidator.h:193
ElectronMcFakeValidator::h1_ele_EseedOP
MonitorElement * h1_ele_EseedOP
Definition: ElectronMcFakeValidator.h:293
ElectronMcFakeValidator::h1_ele_EeleOPout
MonitorElement * h1_ele_EeleOPout
Definition: ElectronMcFakeValidator.h:305
ElectronMcFakeValidator::eopmaxsht
double eopmaxsht
Definition: ElectronMcFakeValidator.h:115
reco::GsfElectron::deltaPhiSeedClusterTrackAtCalo
float deltaPhiSeedClusterTrackAtCalo() const
Definition: GsfElectron.h:237
ElectronMcFakeValidator::h2_ele_EseedOPVsE
MonitorElement * h2_ele_EseedOPVsE
Definition: ElectronMcFakeValidator.h:298
ElectronMcFakeValidator::isoFromDepsEcalReduced04Tag_
edm::EDGetTokenT< edm::ValueMap< double > > isoFromDepsEcalReduced04Tag_
Definition: ElectronMcFakeValidator.h:56
ElectronMcFakeValidator::h2_ele_foundHitsVsEta
MonitorElement * h2_ele_foundHitsVsEta
Definition: ElectronMcFakeValidator.h:255
ElectronMcFakeValidator::h1_ele_EoPout_all
MonitorElement * h1_ele_EoPout_all
Definition: ElectronMcFakeValidator.h:162
ElectronMcFakeValidator::h1_ele_dEtaSc_propVtx
MonitorElement * h1_ele_dEtaSc_propVtx
Definition: ElectronMcFakeValidator.h:312
EDAnalyzer.h
ElectronMcFakeValidator::h2_ele_vertexPtVsPhi
MonitorElement * h2_ele_vertexPtVsPhi
Definition: ElectronMcFakeValidator.h:192
ElectronMcFakeValidator::h2_ele_seed_drz2VsPt_
MonitorElement * h2_ele_seed_drz2VsPt_
Definition: ElectronMcFakeValidator.h:362
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
ElectronMcFakeValidator::inputInternalPath_
std::string inputInternalPath_
Definition: ElectronMcFakeValidator.h:73
reco::LeafCandidate::pt
double pt() const final
transverse momentum
Definition: LeafCandidate.h:146
reco
fixed size matrix
Definition: AlignmentAlgorithmBase.h:46
ElectronUtilities.h
ElectronMcFakeValidator::h1_ele_superclusterfbrem
MonitorElement * h1_ele_superclusterfbrem
Definition: ElectronMcFakeValidator.h:390
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
ElectronMcFakeValidator::h1_ele_convFlags
MonitorElement * h1_ele_convFlags
Definition: ElectronMcFakeValidator.h:481
ElectronMcFakeValidator::seed_max
double seed_max
Definition: ElectronMcFakeValidator.h:141
infinity
const double infinity
Definition: CSCChamberFitter.cc:10
ElectronMcFakeValidator::h1_ele_provenance_endcaps
MonitorElement * h1_ele_provenance_endcaps
Definition: ElectronMcFakeValidator.h:417
ElectronMcFakeValidator::h2_ele_seed_dphi2VsEta_
MonitorElement * h2_ele_seed_dphi2VsEta_
Definition: ElectronMcFakeValidator.h:355
ElectronMcFakeValidator::h1_ele_tkSumPt_dr03_endcaps
MonitorElement * h1_ele_tkSumPt_dr03_endcaps
Definition: ElectronMcFakeValidator.h:442
ElectronMcFakeValidator::h1_ele_hcalTowerSumEtBc_dr04_depth1_barrel
MonitorElement * h1_ele_hcalTowerSumEtBc_dr04_depth1_barrel
Definition: ElectronMcFakeValidator.h:475
reco::GsfElectron::hadronicOverEm
float hadronicOverEm() const
Definition: GsfElectron.h:508
ElectronMcFakeValidator::core_min
double core_min
Definition: ElectronMcFakeValidator.h:134
electronIsolatorFromEffectiveArea_cfi.gsfElectrons
gsfElectrons
Definition: electronIsolatorFromEffectiveArea_cfi.py:4
ElectronMcFakeValidator::h1_ele_dPhiCl_propOut
MonitorElement * h1_ele_dPhiCl_propOut
Definition: ElectronMcFakeValidator.h:330
ElectronMcFakeValidator::h1_ele_dEtaCl_propOut_barrel
MonitorElement * h1_ele_dEtaCl_propOut_barrel
Definition: ElectronMcFakeValidator.h:325
reco::GsfElectron::fbrem
float fbrem() const
Definition: GsfElectron.h:796
ElectronMcFakeValidator::h2_ele_PinMnPoutVsPt_mode
MonitorElement * h2_ele_PinMnPoutVsPt_mode
Definition: ElectronMcFakeValidator.h:275
ElectronMcFakeValidator::detamatch_min
double detamatch_min
Definition: ElectronMcFakeValidator.h:99
ElectronMcFakeValidator::p1_ele_fbremVsEta_mode
MonitorElement * p1_ele_fbremVsEta_mode
Definition: ElectronMcFakeValidator.h:393
ElectronMcFakeValidator::h2_scl_EtaVsPhi_
MonitorElement * h2_scl_EtaVsPhi_
Definition: ElectronMcFakeValidator.h:228
ElectronMcFakeValidator::h1_scl_EoEmatchingObjectGolden_endcaps
MonitorElement * h1_scl_EoEmatchingObjectGolden_endcaps
Definition: ElectronMcFakeValidator.h:405
ElectronMcFakeValidator::seed_min
double seed_min
Definition: ElectronMcFakeValidator.h:140
BasicClusterFwd.h
ElectronMcFakeValidator::deta_max
double deta_max
Definition: ElectronMcFakeValidator.h:96
ElectronSeedFwd.h
edm::Ref
Definition: AssociativeIterator.h:58
ElectronMcFakeValidator::h1_ele_fbrem_barrel
MonitorElement * h1_ele_fbrem_barrel
Definition: ElectronMcFakeValidator.h:388
reco::GsfElectron::mva_e_pi
float mva_e_pi() const
Definition: GsfElectron.h:736
reco::GsfElectron::dr04HcalTowerSumEtBc
float dr04HcalTowerSumEtBc(int depth=0) const
Definition: GsfElectron.h:609
ElectronMcFakeValidator::h1_ele_foundHits
MonitorElement * h1_ele_foundHits
Definition: ElectronMcFakeValidator.h:252
ElectronMcFakeValidator::h1_recCoreNum_
MonitorElement * h1_recCoreNum_
Definition: ElectronMcFakeValidator.h:147
ElectronMcFakeValidator::h1_scl_full5x5_sigmaIetaIeta_endcaps_
MonitorElement * h1_scl_full5x5_sigmaIetaIeta_endcaps_
Definition: ElectronMcFakeValidator.h:237
ElectronMcFakeValidator::h1_ele_PhiMnPhimatchingObject
MonitorElement * h1_ele_PhiMnPhimatchingObject
Definition: ElectronMcFakeValidator.h:216
ElectronMcFakeValidator::h1_ele_dEtaCl_propOut
MonitorElement * h1_ele_dEtaCl_propOut
Definition: ElectronMcFakeValidator.h:324
reco::GsfElectron::mva_Isolated
float mva_Isolated() const
Definition: GsfElectron.h:735
reco::GsfElectron::hcalOverEcalBc
float hcalOverEcalBc(const ShowerShape &ss, int depth) const
Definition: GsfElectron.h:450
ElectronMcFakeValidator::h1_ele_eta_golden
MonitorElement * h1_ele_eta_golden
Definition: ElectronMcFakeValidator.h:369
ElectronMcFakeValidator::h1_ele_dEtaSc_propVtx_barrel
MonitorElement * h1_ele_dEtaSc_propVtx_barrel
Definition: ElectronMcFakeValidator.h:313
ElectronMcFakeValidator::h1_ele_tkSumPt_dr03_barrel
MonitorElement * h1_ele_tkSumPt_dr03_barrel
Definition: ElectronMcFakeValidator.h:441
ElectronMcFakeValidator::h1_ele_seed_mask_fpix_
MonitorElement * h1_ele_seed_mask_fpix_
Definition: ElectronMcFakeValidator.h:352
ElectronMcFakeValidator::electronSeedCollection_
edm::EDGetTokenT< reco::ElectronSeedCollection > electronSeedCollection_
Definition: ElectronMcFakeValidator.h:45
fileCollector.seed
seed
Definition: fileCollector.py:127
ElectronMcFakeValidator::h2_ele_outerPtVsPt_mode
MonitorElement * h2_ele_outerPtVsPt_mode
Definition: ElectronMcFakeValidator.h:286
ElectronMcFakeValidator::h1_ele_dEtaSc_propVtx_all
MonitorElement * h1_ele_dEtaSc_propVtx_all
Definition: ElectronMcFakeValidator.h:164
ElectronMcFakeValidator::h1_ele_EseedOP_endcaps
MonitorElement * h1_ele_EseedOP_endcaps
Definition: ElectronMcFakeValidator.h:295
ElectronMcFakeValidator::h1_ele_seed_drz2pos_
MonitorElement * h1_ele_seed_drz2pos_
Definition: ElectronMcFakeValidator.h:363
MakerMacros.h
reco::GsfElectron::PflowIsolationVariables::sumPhotonEt
float sumPhotonEt
sum pt of PF photons // old float photonIso ;
Definition: GsfElectron.h:673
ElectronMcFakeValidator::lhits_max
double lhits_max
Definition: ElectronMcFakeValidator.h:89
cms::cuda::bs
bs
Definition: HistoContainer.h:76
ElectronMcFakeValidator::dphimatch_nbin
int dphimatch_nbin
Definition: ElectronMcFakeValidator.h:108
ElectronMcFakeValidator::h1_ele_ecalRecHitSumEt_dr03_barrel
MonitorElement * h1_ele_ecalRecHitSumEt_dr03_barrel
Definition: ElectronMcFakeValidator.h:444
reco::GsfElectron::ambiguousGsfTracksSize
GsfTrackRefVector::size_type ambiguousGsfTracksSize() const
Definition: GsfElectron.h:753
ElectronMcFakeValidator::h2_ele_PhiMnPhimatchingObjectVsPt
MonitorElement * h2_ele_PhiMnPhimatchingObjectVsPt
Definition: ElectronMcFakeValidator.h:220
ElectronMcFakeValidator::h2_ele_dPhiClVsPt_propOut
MonitorElement * h2_ele_dPhiClVsPt_propOut
Definition: ElectronMcFakeValidator.h:335
ElectronMcFakeValidator::h1_ele_HoE_barrel
MonitorElement * h1_ele_HoE_barrel
Definition: ElectronMcFakeValidator.h:376
ElectronMcFakeValidator::maxPt_
double maxPt_
Definition: ElectronMcFakeValidator.h:67
ElectronMcFakeValidator::h1_ele_mva_barrel_isolated
MonitorElement * h1_ele_mva_barrel_isolated
Definition: ElectronMcFakeValidator.h:413
trigger::size_type
uint16_t size_type
Definition: TriggerTypeDefs.h:18
ElectronMcFakeValidator::h1_ele_EoverP_all
MonitorElement * h1_ele_EoverP_all
Definition: ElectronMcFakeValidator.h:160
ElectronMcFakeValidator::h1_scl_E5x5_barrel_
MonitorElement * h1_scl_E5x5_barrel_
Definition: ElectronMcFakeValidator.h:245
ElectronMcFakeValidator::isoFromDepsTk03Tag_
edm::EDGetTokenT< edm::ValueMap< double > > isoFromDepsTk03Tag_
Definition: ElectronMcFakeValidator.h:51
ElectronMcFakeValidator::h1_ele_dPhiSc_propVtx_all
MonitorElement * h1_ele_dPhiSc_propVtx_all
Definition: ElectronMcFakeValidator.h:165
ElectronMcFakeValidator::h1_ele_fbrem_endcaps
MonitorElement * h1_ele_fbrem_endcaps
Definition: ElectronMcFakeValidator.h:389
ElectronMcFakeValidator::fhits_max
double fhits_max
Definition: ElectronMcFakeValidator.h:87
ElectronMcFakeValidator::h1_scl_Phi_
MonitorElement * h1_scl_Phi_
Definition: ElectronMcFakeValidator.h:230
BeamSpot.h
ElectronMcFakeValidator::h1_ele_ecalRecHitSumEt_dr04_endcaps
MonitorElement * h1_ele_ecalRecHitSumEt_dr04_endcaps
Definition: ElectronMcFakeValidator.h:455
ElectronMcFakeValidator::p1_ele_fbremVsEta_mean
MonitorElement * p1_ele_fbremVsEta_mean
Definition: ElectronMcFakeValidator.h:394
ElectronMcFakeValidator::h1_ele_tkSumPt_dr04_barrel
MonitorElement * h1_ele_tkSumPt_dr04_barrel
Definition: ElectronMcFakeValidator.h:451
ElectronMcFakeValidator::h2_ele_EeleOPoutVsEta
MonitorElement * h2_ele_EeleOPoutVsEta
Definition: ElectronMcFakeValidator.h:308
reco::LeafCandidate::py
double py() const final
y coordinate of momentum vector
Definition: LeafCandidate.h:142
ElectronMcFakeValidator::h1_ele_vertexPt
MonitorElement * h1_ele_vertexPt
Definition: ElectronMcFakeValidator.h:190
ElectronMcFakeValidator::h2_ele_PinVsPoutShowering_mean
MonitorElement * h2_ele_PinVsPoutShowering_mean
Definition: ElectronMcFakeValidator.h:399
ElectronMcFakeValidator::phi_min
double phi_min
Definition: ElectronMcFakeValidator.h:103
ElectronMcFakeValidator::opv_min
double opv_min
Definition: ElectronMcFakeValidator.h:128
ElectronMcFakeValidator::isoFromDepsEcalFull03Tag_
edm::EDGetTokenT< edm::ValueMap< double > > isoFromDepsEcalFull03Tag_
Definition: ElectronMcFakeValidator.h:53
ElectronMcFakeValidator::h1_matchingObjectPhi
MonitorElement * h1_matchingObjectPhi
Definition: ElectronMcFakeValidator.h:157
ElectronMcFakeValidator::h1_ele_eta_shower
MonitorElement * h1_ele_eta_shower
Definition: ElectronMcFakeValidator.h:372
ElectronMcFakeValidator::h2_ele_ambiguousTracksVsEta
MonitorElement * h2_ele_ambiguousTracksVsEta
Definition: ElectronMcFakeValidator.h:249
ElectronMcFakeValidator::h2_ele_EtaMnEtamatchingObjectVsPhi
MonitorElement * h2_ele_EtaMnEtamatchingObjectVsPhi
Definition: ElectronMcFakeValidator.h:214
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
ElectronMcFakeValidator::h1_ele_hcalTowerSumEtBc_dr03_depth1
MonitorElement * h1_ele_hcalTowerSumEtBc_dr03_depth1
Definition: ElectronMcFakeValidator.h:468
ElectronMcFakeValidator::h1_ele_tkSumPt_dr04
MonitorElement * h1_ele_tkSumPt_dr04
Definition: ElectronMcFakeValidator.h:450
ElectronMcFakeValidator::h1_scl_SigIEtaIEta_endcaps_
MonitorElement * h1_scl_SigIEtaIEta_endcaps_
Definition: ElectronMcFakeValidator.h:234
ElectronMcFakeValidator::h1_ele_EoPout
MonitorElement * h1_ele_EoPout
Definition: ElectronMcFakeValidator.h:299
dqm::impl::MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:290
reco::GsfElectron::trackMomentumAtVtx
math::XYZVectorF trackMomentumAtVtx() const
Definition: GsfElectron.h:276
reco::BeamSpot
Definition: BeamSpot.h:21
ElectronMcFakeValidator::h2_ele_PinMnPoutVsChi2_mode
MonitorElement * h2_ele_PinMnPoutVsChi2_mode
Definition: ElectronMcFakeValidator.h:277
ElectronMcFakeValidator::h1_ele_superclusterfbrem_endcaps
MonitorElement * h1_ele_superclusterfbrem_endcaps
Definition: ElectronMcFakeValidator.h:392
ElectronMcFakeValidator::h1_ele_foundHits_barrel
MonitorElement * h1_ele_foundHits_barrel
Definition: ElectronMcFakeValidator.h:253
reco::GsfElectron::ecalDrivenSeed
bool ecalDrivenSeed() const
Definition: GsfElectron.h:166
ElectronMcFakeValidator::h1_ele_photonRelativeIso_endcaps
MonitorElement * h1_ele_photonRelativeIso_endcaps
Definition: ElectronMcFakeValidator.h:438
ElectronDqmAnalyzerBase::setBookIndex
void setBookIndex(short)
Definition: ElectronDqmAnalyzerBase.cc:32
ElectronMcFakeValidator::electronTrackCollection_
edm::EDGetTokenT< reco::GsfTrackCollection > electronTrackCollection_
Definition: ElectronMcFakeValidator.h:44
ElectronMcFakeValidator::h2_ele_dEtaClVsPt_propOut
MonitorElement * h2_ele_dEtaClVsPt_propOut
Definition: ElectronMcFakeValidator.h:329
reco::GsfElectron::isEEDeeGap
bool isEEDeeGap() const
Definition: GsfElectron.h:344
ElectronMcFakeValidator::h1_ele_HoE_bc
MonitorElement * h1_ele_HoE_bc
Definition: ElectronMcFakeValidator.h:375
ElectronMcFakeValidator::h1_ele_HoE
MonitorElement * h1_ele_HoE
Definition: ElectronMcFakeValidator.h:374
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
ElectronMcFakeValidator::h2_ele_E2mnE1vsMee_egeg_all
MonitorElement * h2_ele_E2mnE1vsMee_egeg_all
Definition: ElectronMcFakeValidator.h:177
ElectronMcFakeValidator::popmatching_nbin
int popmatching_nbin
Definition: ElectronMcFakeValidator.h:122
ElectronMcFakeValidator::isoFromDepsHcal04Tag_
edm::EDGetTokenT< edm::ValueMap< double > > isoFromDepsHcal04Tag_
Definition: ElectronMcFakeValidator.h:58
ElectronMcFakeValidator::h1_scl_EoEmatchingObject_barrel
MonitorElement * h1_scl_EoEmatchingObject_barrel
Definition: ElectronMcFakeValidator.h:223
reco::GsfElectronRef
edm::Ref< GsfElectronCollection > GsfElectronRef
reference to an object in a collection of GsfElectron objects
Definition: GsfElectronFwd.h:21
ElectronMcFakeValidator::h2_ele_PtinVsPtoutShowering_mode
MonitorElement * h2_ele_PtinVsPtoutShowering_mode
Definition: ElectronMcFakeValidator.h:401
hgcalPlots.xtitle
xtitle
Definition: hgcalPlots.py:93
ElectronMcFakeValidator::h1_ele_dEtaCl_propOut_endcaps
MonitorElement * h1_ele_dEtaCl_propOut_endcaps
Definition: ElectronMcFakeValidator.h:326
ElectronMcFakeValidator::~ElectronMcFakeValidator
~ElectronMcFakeValidator() override
Definition: ElectronMcFakeValidator.cc:2405
ElectronMcFakeValidator::h1_ele_hcalTowerSumEt_dr03_depth2
MonitorElement * h1_ele_hcalTowerSumEt_dr03_depth2
Definition: ElectronMcFakeValidator.h:449
ElectronMcFakeValidator::h1_ele_dEtaEleCl_propOut_barrel
MonitorElement * h1_ele_dEtaEleCl_propOut_barrel
Definition: ElectronMcFakeValidator.h:337
reco::GsfElectron::convDist
float convDist() const
Definition: GsfElectron.h:653
ElectronMcFakeValidator::h1_ele_vertexP
MonitorElement * h1_ele_vertexP
Definition: ElectronMcFakeValidator.h:189
ElectronMcFakeValidator::h1_ele_chi2
MonitorElement * h1_ele_chi2
Definition: ElectronMcFakeValidator.h:264
ElectronMcFakeValidator::isoFromDepsEcalReduced03Tag_
edm::EDGetTokenT< edm::ValueMap< double > > isoFromDepsEcalReduced03Tag_
Definition: ElectronMcFakeValidator.h:55
reco::GsfElectron
Definition: GsfElectron.h:35
ElectronMcFakeValidator::h2_ele_lostHitsVsPhi
MonitorElement * h2_ele_lostHitsVsPhi
Definition: ElectronMcFakeValidator.h:262
ElectronMcFakeValidator::h1_ele_chargedHadronIso_endcaps
MonitorElement * h1_ele_chargedHadronIso_endcaps
Definition: ElectronMcFakeValidator.h:422
ElectronMcFakeValidator::h1_ele_vertexZ
MonitorElement * h1_ele_vertexZ
Definition: ElectronMcFakeValidator.h:199
ElectronMcFakeValidator::h1_ele_superclusterfbrem_barrel
MonitorElement * h1_ele_superclusterfbrem_barrel
Definition: ElectronMcFakeValidator.h:391
ElectronMcFakeValidator::h1_ele_charge
MonitorElement * h1_ele_charge
Definition: ElectronMcFakeValidator.h:185
GsfElectron.h
ElectronMcFakeValidator::h1_ele_mva_isolated
MonitorElement * h1_ele_mva_isolated
Definition: ElectronMcFakeValidator.h:412
ElectronMcFakeValidator::h1_ele_neutralHadronRelativeIso_endcaps
MonitorElement * h1_ele_neutralHadronRelativeIso_endcaps
Definition: ElectronMcFakeValidator.h:435
ElectronMcFakeValidator::h1_ele_mva
MonitorElement * h1_ele_mva
Definition: ElectronMcFakeValidator.h:409
PbPb_ZMuSkimMuonDPG_cff.deltaR
deltaR
Definition: PbPb_ZMuSkimMuonDPG_cff.py:63
ElectronMcFakeValidator::h2_ele_outerPtVsPhi_mode
MonitorElement * h2_ele_outerPtVsPhi_mode
Definition: ElectronMcFakeValidator.h:285
ElectronMcFakeValidator::h2_ele_dEtaScVsEta_propVtx
MonitorElement * h2_ele_dEtaScVsEta_propVtx
Definition: ElectronMcFakeValidator.h:315
ElectronMcFakeValidator::popmatching_min
double popmatching_min
Definition: ElectronMcFakeValidator.h:123
ElectronMcFakeValidator::h2_ele_dPhiScVsEta_propVtx
MonitorElement * h2_ele_dPhiScVsEta_propVtx
Definition: ElectronMcFakeValidator.h:321
reco::GsfElectron::superClusterFbrem
float superClusterFbrem() const
Definition: GsfElectron.h:790
ElectronMcFakeValidator::h2_ele_PoPmatchingObjectVsPhi
MonitorElement * h2_ele_PoPmatchingObjectVsPhi
Definition: ElectronMcFakeValidator.h:207
ElectronMcFakeValidator::h1_ele_lostHits_barrel
MonitorElement * h1_ele_lostHits_barrel
Definition: ElectronMcFakeValidator.h:259
EcalSubdetector.h
ElectronMcFakeValidator::h2_ele_EseedOPVsPhi
MonitorElement * h2_ele_EseedOPVsPhi
Definition: ElectronMcFakeValidator.h:297
ElectronMcFakeValidator::h2_ele_EoPoutVsPhi
MonitorElement * h2_ele_EoPoutVsPhi
Definition: ElectronMcFakeValidator.h:303
ElectronMcFakeValidator::h1_ele_hcalTowerSumEt_dr03_depth1
MonitorElement * h1_ele_hcalTowerSumEt_dr03_depth1
Definition: ElectronMcFakeValidator.h:446
ElectronMcFakeValidator::h1_scl_E1x5_endcaps_
MonitorElement * h1_scl_E1x5_endcaps_
Definition: ElectronMcFakeValidator.h:240
ElectronMcFakeValidator::h2_ele_EoPVsPhi
MonitorElement * h2_ele_EoPVsPhi
Definition: ElectronMcFakeValidator.h:291
ElectronMcFakeValidator::h2_ele_PtinVsPtoutShowering_mean
MonitorElement * h2_ele_PtinVsPtoutShowering_mean
Definition: ElectronMcFakeValidator.h:403
ElectronMcFakeValidator::h1_ele_hcalTowerSumEtBc_dr04_depth2
MonitorElement * h1_ele_hcalTowerSumEtBc_dr04_depth2
Definition: ElectronMcFakeValidator.h:477
ElectronMcFakeValidator::h1_ele_hcalTowerSumEtBc_dr03_depth1_endcaps
MonitorElement * h1_ele_hcalTowerSumEtBc_dr03_depth1_endcaps
Definition: ElectronMcFakeValidator.h:470
ElectronMcFakeValidator::eta_min
double eta_min
Definition: ElectronMcFakeValidator.h:92
ElectronMcFakeValidator::p2D_nbin
int p2D_nbin
Definition: ElectronMcFakeValidator.h:80
ElectronMcFakeValidator::h1_ele_hcalTowerSumEt_dr04_depth2
MonitorElement * h1_ele_hcalTowerSumEt_dr04_depth2
Definition: ElectronMcFakeValidator.h:459
ElectronMcFakeValidator::h2_ele_PinVsPoutGolden_mode
MonitorElement * h2_ele_PinVsPoutGolden_mode
Definition: ElectronMcFakeValidator.h:396
GsfElectronFwd.h
ElectronMcFakeValidator::detamatch2D_nbin
int detamatch2D_nbin
Definition: ElectronMcFakeValidator.h:98
edm::ParameterSet
Definition: ParameterSet.h:47
ElectronMcFakeValidator::h2_scl_EtVsPhi_
MonitorElement * h2_scl_EtVsPhi_
Definition: ElectronMcFakeValidator.h:227
ElectronMcFakeValidator::h1_ele_dPhiEleCl_propOut_barrel
MonitorElement * h1_ele_dPhiEleCl_propOut_barrel
Definition: ElectronMcFakeValidator.h:343
ElectronMcFakeValidator::pt_nbin
int pt_nbin
Definition: ElectronMcFakeValidator.h:82
reco::GsfElectron::deltaEtaSuperClusterTrackAtVtx
float deltaEtaSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:233
ElectronMcFakeValidator::eta_nbin
int eta_nbin
Definition: ElectronMcFakeValidator.h:90
ElectronMcFakeValidator::set_StatOverflowFlag
bool set_StatOverflowFlag
Definition: ElectronMcFakeValidator.h:126
ElectronMcFakeValidator::h1_ele_hcalDepth1OverEcalBc_barrel
MonitorElement * h1_ele_hcalDepth1OverEcalBc_barrel
Definition: ElectronMcFakeValidator.h:462
ElectronMcFakeValidator::h1_ele_eta_narrow
MonitorElement * h1_ele_eta_narrow
Definition: ElectronMcFakeValidator.h:371
ElectronMcFakeValidator::eta_max
double eta_max
Definition: ElectronMcFakeValidator.h:93
Event.h
reco::LeafCandidate::eta
double eta() const final
momentum pseudorapidity
Definition: LeafCandidate.h:152
ElectronMcFakeValidator::h1_matchingObjectEta
MonitorElement * h1_matchingObjectEta
Definition: ElectronMcFakeValidator.h:153
ElectronMcFakeValidator::h1_ele_neutralHadronIso
MonitorElement * h1_ele_neutralHadronIso
Definition: ElectronMcFakeValidator.h:423
ElectronMcFakeValidator::h2_ele_PoPmatchingObjectVsEta
MonitorElement * h2_ele_PoPmatchingObjectVsEta
Definition: ElectronMcFakeValidator.h:206
ElectronMcFakeValidator::detamatch_max
double detamatch_max
Definition: ElectronMcFakeValidator.h:100
ElectronMcFakeValidator::h1_ele_mee_os
MonitorElement * h1_ele_mee_os
Definition: ElectronMcFakeValidator.h:174
ElectronMcFakeValidator::h1_ele_fbrem
MonitorElement * h1_ele_fbrem
Definition: ElectronMcFakeValidator.h:387
ElectronMcFakeValidator::h1_ele_hcalDepth1OverEcalBc
MonitorElement * h1_ele_hcalDepth1OverEcalBc
Definition: ElectronMcFakeValidator.h:461
reco::GsfElectron::full5x5_sigmaIetaIeta
float full5x5_sigmaIetaIeta() const
Definition: GsfElectron.h:479
ElectronMcFakeValidator::h1_ele_chargedHadronIso_barrel
MonitorElement * h1_ele_chargedHadronIso_barrel
Definition: ElectronMcFakeValidator.h:421
reco::GsfElectron::scE5x5
float scE5x5() const
Definition: GsfElectron.h:507
ElectronMcFakeValidator::popmatching_max
double popmatching_max
Definition: ElectronMcFakeValidator.h:124
ElectronMcFakeValidator::h1_ele_photonIso_endcaps
MonitorElement * h1_ele_photonIso_endcaps
Definition: ElectronMcFakeValidator.h:428
edm::Ref::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:238
ElectronMcFakeValidator::h1_ele_PoPmatchingObject_barrel
MonitorElement * h1_ele_PoPmatchingObject_barrel
Definition: ElectronMcFakeValidator.h:209
ElectronMcFakeValidator::eop_nbin
int eop_nbin
Definition: ElectronMcFakeValidator.h:112
ElectronMcFakeValidator::h1_ele_hcalTowerSumEtBc_dr04_depth1
MonitorElement * h1_ele_hcalTowerSumEtBc_dr04_depth1
Definition: ElectronMcFakeValidator.h:474
ElectronMcFakeValidator::h1_ele_seed_subdet2_
MonitorElement * h1_ele_seed_subdet2_
Definition: ElectronMcFakeValidator.h:349
ElectronMcFakeValidator::h1_ele_ambiguousTracks
MonitorElement * h1_ele_ambiguousTracks
Definition: ElectronMcFakeValidator.h:248
ElectronMcFakeValidator::h1_ele_EtaMnEtamatchingObject
MonitorElement * h1_ele_EtaMnEtamatchingObject
Definition: ElectronMcFakeValidator.h:212
ElectronMcFakeValidator::matchingObjectCollection_
edm::EDGetTokenT< reco::GenJetCollection > matchingObjectCollection_
Definition: ElectronMcFakeValidator.h:46
ElectronMcFakeValidator::h1_ele_neutralHadronRelativeIso_barrel
MonitorElement * h1_ele_neutralHadronRelativeIso_barrel
Definition: ElectronMcFakeValidator.h:434
ElectronMcFakeValidator::offlineVerticesCollection_
edm::EDGetTokenT< reco::VertexCollection > offlineVerticesCollection_
Definition: ElectronMcFakeValidator.h:47
ElectronMcFakeValidator::dphi_nbin
int dphi_nbin
Definition: ElectronMcFakeValidator.h:105
ElectronMcFakeValidator::h1_ele_convRadius_all
MonitorElement * h1_ele_convRadius_all
Definition: ElectronMcFakeValidator.h:488
iEvent
int iEvent
Definition: GenABIO.cc:224
ElectronMcFakeValidator::h1_ele_eta_bbrem
MonitorElement * h1_ele_eta_bbrem
Definition: ElectronMcFakeValidator.h:370
ElectronMcFakeValidator::analyze
void analyze(const edm::Event &e, const edm::EventSetup &c) override
Definition: ElectronMcFakeValidator.cc:2411
reco::GsfElectron::eSeedClusterOverP
float eSeedClusterOverP() const
Definition: GsfElectron.h:230
reco::GsfElectron::dr04TkSumPt
float dr04TkSumPt() const
Definition: GsfElectron.h:605
ElectronMcFakeValidator::h2_ele_PinVsPoutGolden_mean
MonitorElement * h2_ele_PinVsPoutGolden_mean
Definition: ElectronMcFakeValidator.h:398
GsfTrack.h
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
ElectronMcFakeValidator::electronCollection_
edm::EDGetTokenT< reco::GsfElectronCollection > electronCollection_
Definition: ElectronMcFakeValidator.h:42
ElectronMcFakeValidator::h2_ele_outerPtVsEta_mode
MonitorElement * h2_ele_outerPtVsEta_mode
Definition: ElectronMcFakeValidator.h:284
ElectronMcFakeValidator::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: ElectronMcFakeValidator.cc:463
ElectronMcFakeValidator::eta2D_nbin
int eta2D_nbin
Definition: ElectronMcFakeValidator.h:91
ElectronMcFakeValidator::h1_ele_mee_all
MonitorElement * h1_ele_mee_all
Definition: ElectronMcFakeValidator.h:173
ElectronMcFakeValidator::h1_scl_ESFrac_endcaps
MonitorElement * h1_scl_ESFrac_endcaps
Definition: ElectronMcFakeValidator.h:385
ElectronMcFakeValidator::dphimatch_max
double dphimatch_max
Definition: ElectronMcFakeValidator.h:111
reco::LeafCandidate::charge
int charge() const final
electric charge
Definition: LeafCandidate.h:106
reco::GsfElectron::isEBEtaGap
bool isEBEtaGap() const
Definition: GsfElectron.h:341
ElectronMcFakeValidator::isoFromDepsTk04Tag_
edm::EDGetTokenT< edm::ValueMap< double > > isoFromDepsTk04Tag_
Definition: ElectronMcFakeValidator.h:52
ElectronMcFakeValidator::h2_ele_HoEVsPhi
MonitorElement * h2_ele_HoEVsPhi
Definition: ElectronMcFakeValidator.h:382
ElectronMcFakeValidator::h2_ele_seed_drz2posVsPt_
MonitorElement * h2_ele_seed_drz2posVsPt_
Definition: ElectronMcFakeValidator.h:365
edm::EventSetup
Definition: EventSetup.h:58
ElectronMcFakeValidator::h1_scl_SigIEtaIEta_
MonitorElement * h1_scl_SigIEtaIEta_
Definition: ElectronMcFakeValidator.h:232
ElectronMcFakeValidator::h1_scl_SigIEtaIEta_barrel_
MonitorElement * h1_scl_SigIEtaIEta_barrel_
Definition: ElectronMcFakeValidator.h:233
ElectronMcFakeValidator::deta_nbin
int deta_nbin
Definition: ElectronMcFakeValidator.h:94
ElectronMcFakeValidator::h2_ele_dPhiScVsPhi_propVtx
MonitorElement * h2_ele_dPhiScVsPhi_propVtx
Definition: ElectronMcFakeValidator.h:322
ElectronMcFakeValidator::h1_ele_dPhiCl_propOut_endcaps
MonitorElement * h1_ele_dPhiCl_propOut_endcaps
Definition: ElectronMcFakeValidator.h:332
ElectronDqmAnalyzerBase::setBookStatOverflowFlag
void setBookStatOverflowFlag(const bool &)
Definition: ElectronDqmAnalyzerBase.cc:36
ElectronMcFakeValidator::p_max
double p_max
Definition: ElectronMcFakeValidator.h:81
ElectronMcFakeValidator::h1_ele_hcalTowerSumEt_dr03_depth1_endcaps
MonitorElement * h1_ele_hcalTowerSumEt_dr03_depth1_endcaps
Definition: ElectronMcFakeValidator.h:448
ElectronMcFakeValidator::h2_ele_chi2VsEta
MonitorElement * h2_ele_chi2VsEta
Definition: ElectronMcFakeValidator.h:267
ElectronMcFakeValidator::track_max
double track_max
Definition: ElectronMcFakeValidator.h:138
ElectronMcFakeValidator::mee_max
double mee_max
Definition: ElectronMcFakeValidator.h:118
ElectronMcFakeValidator::h2_ele_seed_drz2posVsEta_
MonitorElement * h2_ele_seed_drz2posVsEta_
Definition: ElectronMcFakeValidator.h:364
ElectronMcFakeValidator::hoe_min
double hoe_min
Definition: ElectronMcFakeValidator.h:120
ElectronMcFakeValidator::h1_ele_classes
MonitorElement * h1_ele_classes
Definition: ElectronMcFakeValidator.h:367
ElectronMcFakeValidator::h2_ele_dPhiClVsPhi_propOut
MonitorElement * h2_ele_dPhiClVsPhi_propOut
Definition: ElectronMcFakeValidator.h:334
ElectronMcFakeValidator::h1_scl_EoEmatchingObjectShowering_endcaps
MonitorElement * h1_scl_EoEmatchingObjectShowering_endcaps
Definition: ElectronMcFakeValidator.h:407
ElectronMcFakeValidator::hoe_max
double hoe_max
Definition: ElectronMcFakeValidator.h:121
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
ElectronMcFakeValidator::h2_ele_HoEVsEta
MonitorElement * h2_ele_HoEVsEta
Definition: ElectronMcFakeValidator.h:381
ElectronMcFakeValidator::h1_ele_ecalRecHitSumEt_dr04_barrel
MonitorElement * h1_ele_ecalRecHitSumEt_dr04_barrel
Definition: ElectronMcFakeValidator.h:454
ElectronMcFakeValidator::h1_ele_dPhiSc_propVtx_endcaps
MonitorElement * h1_ele_dPhiSc_propVtx_endcaps
Definition: ElectronMcFakeValidator.h:320
reco::GsfElectron::PflowIsolationVariables::sumNeutralHadronEt
float sumNeutralHadronEt
sum pt of neutral hadrons // old float neutralHadronIso ;
Definition: GsfElectron.h:672
ElectronMcFakeValidator::h2_ele_chargeVsEta
MonitorElement * h2_ele_chargeVsEta
Definition: ElectronMcFakeValidator.h:186
ElectronMcFakeValidator::h1_matchingObjectPt
MonitorElement * h1_matchingObjectPt
Definition: ElectronMcFakeValidator.h:156
ElectronMcFakeValidator::h1_scl_E1x5_
MonitorElement * h1_scl_E1x5_
Definition: ElectronMcFakeValidator.h:238
ElectronMcFakeValidator::h1_ele_vertexEta_all
MonitorElement * h1_ele_vertexEta_all
Definition: ElectronMcFakeValidator.h:171
reco::LeafCandidate::vertex
const Point & vertex() const override
vertex position (overwritten by PF...)
Definition: LeafCandidate.h:165
ElectronMcFakeValidator::h1_scl_E2x5max_barrel_
MonitorElement * h1_scl_E2x5max_barrel_
Definition: ElectronMcFakeValidator.h:242
ElectronMcFakeValidator::h1_ele_vertexX
MonitorElement * h1_ele_vertexX
Definition: ElectronMcFakeValidator.h:197
ElectronMcFakeValidator::phi_max
double phi_max
Definition: ElectronMcFakeValidator.h:104
ElectronMcFakeValidator::h1_ele_mva_endcaps
MonitorElement * h1_ele_mva_endcaps
Definition: ElectronMcFakeValidator.h:411
ElectronMcFakeValidator::h2_ele_dEtaScVsPt_propVtx
MonitorElement * h2_ele_dEtaScVsPt_propVtx
Definition: ElectronMcFakeValidator.h:317
ElectronMcFakeValidator.h
ElectronMcFakeValidator::h1_ele_PoPmatchingObject
MonitorElement * h1_ele_PoPmatchingObject
Definition: ElectronMcFakeValidator.h:205
ValueMap.h
ElectronMcFakeValidator::h1_recOfflineVertices_
MonitorElement * h1_recOfflineVertices_
Definition: ElectronMcFakeValidator.h:151
VertexFwd.h
ElectronMcFakeValidator::h2_ele_PinMnPoutVsEta_mode
MonitorElement * h2_ele_PinMnPoutVsEta_mode
Definition: ElectronMcFakeValidator.h:273
ElectronMcFakeValidator::h2_ele_EeleOPoutVsPhi
MonitorElement * h2_ele_EeleOPoutVsPhi
Definition: ElectronMcFakeValidator.h:309
ElectronMcFakeValidator::core_nbin
int core_nbin
Definition: ElectronMcFakeValidator.h:133
ElectronDqmAnalyzerBase::bookP1
MonitorElement * bookP1(DQMStore::IBooker &, const std::string &name, const std::string &title, int nchX, double lowX, double highX, double lowY, double highY, const std::string &titleX="", const std::string &titleY="", Option_t *option="E1 P")
Definition: ElectronDqmAnalyzerBase.cc:171
ElectronMcFakeValidator::h1_scl_E1x5_barrel_
MonitorElement * h1_scl_E1x5_barrel_
Definition: ElectronMcFakeValidator.h:239
ElectronMcFakeValidator::outputFile_
std::string outputFile_
Definition: ElectronMcFakeValidator.h:72
ElectronMcFakeValidator::h1_ele_neutralHadronIso_barrel
MonitorElement * h1_ele_neutralHadronIso_barrel
Definition: ElectronMcFakeValidator.h:424
ElectronMcFakeValidator::opv_max
double opv_max
Definition: ElectronMcFakeValidator.h:129
ElectronMcFakeValidator::h1_ele_EeleOPout_endcaps
MonitorElement * h1_ele_EeleOPout_endcaps
Definition: ElectronMcFakeValidator.h:307
ElectronMcFakeValidator::h1_ele_vertexPt_all
MonitorElement * h1_ele_vertexPt_all
Definition: ElectronMcFakeValidator.h:172
ElectronMcFakeValidator::h1_ele_ecalRecHitSumEt_dr04
MonitorElement * h1_ele_ecalRecHitSumEt_dr04
Definition: ElectronMcFakeValidator.h:453
ElectronMcFakeValidator::h2_ele_dEtaEleClVsEta_propOut
MonitorElement * h2_ele_dEtaEleClVsEta_propOut
Definition: ElectronMcFakeValidator.h:339
ElectronMcFakeValidator::outputInternalPath_
std::string outputInternalPath_
Definition: ElectronMcFakeValidator.h:74
ElectronMcFakeValidator::h1_ele_dPhiEleCl_propOut_endcaps
MonitorElement * h1_ele_dPhiEleCl_propOut_endcaps
Definition: ElectronMcFakeValidator.h:344
ElectronMcFakeValidator::h1_ele_seed_mask_
MonitorElement * h1_ele_seed_mask_
Definition: ElectronMcFakeValidator.h:350
ElectronMcFakeValidator::h1_ele_hcalDepth2OverEcalBc_barrel
MonitorElement * h1_ele_hcalDepth2OverEcalBc_barrel
Definition: ElectronMcFakeValidator.h:465
ElectronMcFakeValidator::h1_ele_convFlags_all
MonitorElement * h1_ele_convFlags_all
Definition: ElectronMcFakeValidator.h:482
ElectronMcFakeValidator::pt2D_nbin
int pt2D_nbin
Definition: ElectronMcFakeValidator.h:83
ElectronMcFakeValidator::h2_ele_chi2VsPt
MonitorElement * h2_ele_chi2VsPt
Definition: ElectronMcFakeValidator.h:269
ElectronMcFakeValidator::h1_ele_outerP_mode
MonitorElement * h1_ele_outerP_mode
Definition: ElectronMcFakeValidator.h:280
ElectronMcFakeValidator::h2_ele_dEtaEleClVsPt_propOut
MonitorElement * h2_ele_dEtaEleClVsPt_propOut
Definition: ElectronMcFakeValidator.h:341
ElectronMcFakeValidator::h1_ele_photonRelativeIso
MonitorElement * h1_ele_photonRelativeIso
Definition: ElectronMcFakeValidator.h:436
ElectronMcFakeValidator::h1_matchingObjectZ
MonitorElement * h1_matchingObjectZ
Definition: ElectronMcFakeValidator.h:158
ElectronMcFakeValidator::h1_ele_convRadius
MonitorElement * h1_ele_convRadius
Definition: ElectronMcFakeValidator.h:487
ElectronMcFakeValidator::dphi_max
double dphi_max
Definition: ElectronMcFakeValidator.h:107
ElectronMcFakeValidator::h1_ele_dPhiEleCl_propOut
MonitorElement * h1_ele_dPhiEleCl_propOut
Definition: ElectronMcFakeValidator.h:342
reco::GsfElectron::convRadius
float convRadius() const
Definition: GsfElectron.h:655
ElectronMcFakeValidator::h2_ele_lostHitsVsPt
MonitorElement * h2_ele_lostHitsVsPt
Definition: ElectronMcFakeValidator.h:263
ElectronMcFakeValidator::h1_ele_neutralHadronRelativeIso
MonitorElement * h1_ele_neutralHadronRelativeIso
Definition: ElectronMcFakeValidator.h:433
reco::LeafCandidate::phi
double phi() const final
momentum azimuthal angle
Definition: LeafCandidate.h:148
ElectronMcFakeValidator::h1_ele_EeleOPout_all
MonitorElement * h1_ele_EeleOPout_all
Definition: ElectronMcFakeValidator.h:163
SuperClusterFwd.h
ElectronMcFakeValidator::hoe_nbin
int hoe_nbin
Definition: ElectronMcFakeValidator.h:119
ElectronMcFakeValidator::h1_ele_PinMnPout_mode
MonitorElement * h1_ele_PinMnPout_mode
Definition: ElectronMcFakeValidator.h:272
ElectronMcFakeValidator::h1_ele_dPhiCl_propOut_all
MonitorElement * h1_ele_dPhiCl_propOut_all
Definition: ElectronMcFakeValidator.h:167
ElectronMcFakeValidator::h2_ele_seed_drz2VsEta_
MonitorElement * h2_ele_seed_drz2VsEta_
Definition: ElectronMcFakeValidator.h:361
Vertex.h
ElectronMcFakeValidator::h1_ele_mva_barrel
MonitorElement * h1_ele_mva_barrel
Definition: ElectronMcFakeValidator.h:410
ElectronMcFakeValidator::h2_ele_seed_dphi2posVsEta_
MonitorElement * h2_ele_seed_dphi2posVsEta_
Definition: ElectronMcFakeValidator.h:358
reco::GsfElectron::isEBEEGap
bool isEBEEGap() const
Definition: GsfElectron.h:339
ElectronMcFakeValidator::h2_ele_vertexTIPVsEta
MonitorElement * h2_ele_vertexTIPVsEta
Definition: ElectronMcFakeValidator.h:201
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
ElectronMcFakeValidator::h2_ele_EoPoutVsEta
MonitorElement * h2_ele_EoPoutVsEta
Definition: ElectronMcFakeValidator.h:302
ElectronMcFakeValidator::h2_ele_foundHitsVsPt
MonitorElement * h2_ele_foundHitsVsPt
Definition: ElectronMcFakeValidator.h:257
ElectronMcFakeValidator::h1_scl_E5x5_
MonitorElement * h1_scl_E5x5_
Definition: ElectronMcFakeValidator.h:244
ElectronMcFakeValidator::h1_ele_vertexTIP
MonitorElement * h1_ele_vertexTIP
Definition: ElectronMcFakeValidator.h:200
ElectronMcFakeValidator::h2_ele_chargeVsPt
MonitorElement * h2_ele_chargeVsPt
Definition: ElectronMcFakeValidator.h:188
ElectronMcFakeValidator::h1_scl_EoEmatchingObject_endcaps
MonitorElement * h1_scl_EoEmatchingObject_endcaps
Definition: ElectronMcFakeValidator.h:224
reco::LeafCandidate::p
double p() const final
magnitude of momentum vector
Definition: LeafCandidate.h:123
ElectronMcFakeValidator::h1_ele_tkSumPt_dr03
MonitorElement * h1_ele_tkSumPt_dr03
Definition: ElectronMcFakeValidator.h:440
ElectronMcFakeValidator::deltaR_
double deltaR_
Definition: ElectronMcFakeValidator.h:69
ElectronMcFakeValidator::pt_max
double pt_max
Definition: ElectronMcFakeValidator.h:85
ElectronMcFakeValidator::h1_ele_hcalTowerSumEt_dr04_depth1_barrel
MonitorElement * h1_ele_hcalTowerSumEt_dr04_depth1_barrel
Definition: ElectronMcFakeValidator.h:457
reco::GsfElectron::deltaEtaEleClusterTrackAtCalo
float deltaEtaEleClusterTrackAtCalo() const
Definition: GsfElectron.h:235
ElectronMcFakeValidator::h2_ele_dEtaClVsEta_propOut
MonitorElement * h2_ele_dEtaClVsEta_propOut
Definition: ElectronMcFakeValidator.h:327
ElectronMcFakeValidator::h1_ele_hcalTowerSumEtBc_dr04_depth2_barrel
MonitorElement * h1_ele_hcalTowerSumEtBc_dr04_depth2_barrel
Definition: ElectronMcFakeValidator.h:478
reco::GsfElectron::eSuperClusterOverP
float eSuperClusterOverP() const
Definition: GsfElectron.h:229
reco::GsfElectron::classification
Classification classification() const
Definition: GsfElectron.h:792
ElectronMcFakeValidator::h2_ele_dPhiScVsPt_propVtx
MonitorElement * h2_ele_dPhiScVsPt_propVtx
Definition: ElectronMcFakeValidator.h:323
ElectronMcFakeValidator::h1_ele_photonRelativeIso_barrel
MonitorElement * h1_ele_photonRelativeIso_barrel
Definition: ElectronMcFakeValidator.h:437
ElectronMcFakeValidator::h1_ele_dPhiSc_propVtx_barrel
MonitorElement * h1_ele_dPhiSc_propVtx_barrel
Definition: ElectronMcFakeValidator.h:319
ElectronMcFakeValidator::core_max
double core_max
Definition: ElectronMcFakeValidator.h:135
reco::GsfElectron::eEleClusterOverPout
float eEleClusterOverPout() const
Definition: GsfElectron.h:232
ElectronMcFakeValidator::h1_scl_E5x5_endcaps_
MonitorElement * h1_scl_E5x5_endcaps_
Definition: ElectronMcFakeValidator.h:246
ElectronMcFakeValidator::h1_ele_vertexY
MonitorElement * h1_ele_vertexY
Definition: ElectronMcFakeValidator.h:198
reco::GsfElectron::trackerDrivenSeed
bool trackerDrivenSeed() const
Definition: GsfElectron.h:167
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
edm::RefToBase< TrajectorySeed >
ElectronMcFakeValidator::h2_ele_ambiguousTracksVsPhi
MonitorElement * h2_ele_ambiguousTracksVsPhi
Definition: ElectronMcFakeValidator.h:250
reco::GsfElectron::isEBPhiGap
bool isEBPhiGap() const
Definition: GsfElectron.h:342
reco::GsfElectron::dr03HcalTowerSumEt
float dr03HcalTowerSumEt(int depth=0) const
Definition: GsfElectron.h:584
reco::GsfElectron::eSeedClusterOverPout
float eSeedClusterOverPout() const
Definition: GsfElectron.h:231
reco::GsfElectron::superCluster
SuperClusterRef superCluster() const override
reference to a SuperCluster
Definition: GsfElectron.h:163
ElectronMcFakeValidator::h1_ele_EoPout_endcaps
MonitorElement * h1_ele_EoPout_endcaps
Definition: ElectronMcFakeValidator.h:301
ElectronMcFakeValidator::h1_ele_hcalTowerSumEtBc_dr03_depth2_barrel
MonitorElement * h1_ele_hcalTowerSumEtBc_dr03_depth2_barrel
Definition: ElectronMcFakeValidator.h:472
dqm::implementation::IBooker
Definition: DQMStore.h:43
ElectronMcFakeValidator::h1_ele_tkSumPt_dr04_endcaps
MonitorElement * h1_ele_tkSumPt_dr04_endcaps
Definition: ElectronMcFakeValidator.h:452
ElectronMcFakeValidator::h1_ele_hcalTowerSumEt_dr04_depth1_endcaps
MonitorElement * h1_ele_hcalTowerSumEt_dr04_depth1_endcaps
Definition: ElectronMcFakeValidator.h:458
ElectronMcFakeValidator::h1_ele_EoP
MonitorElement * h1_ele_EoP
Definition: ElectronMcFakeValidator.h:287
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:29
ElectronMcFakeValidator::h1_ele_ecalRecHitSumEt_dr03_endcaps
MonitorElement * h1_ele_ecalRecHitSumEt_dr03_endcaps
Definition: ElectronMcFakeValidator.h:445
ElectronMcFakeValidator::h2_ele_PinMnPoutVsE_mode
MonitorElement * h2_ele_PinMnPoutVsE_mode
Definition: ElectronMcFakeValidator.h:276
ztail.d
d
Definition: ztail.py:151
ElectronMcFakeValidator::h1_ele_chargedHadronRelativeIso
MonitorElement * h1_ele_chargedHadronRelativeIso
Definition: ElectronMcFakeValidator.h:430
ElectronMcFakeValidator::isoFromDepsEcalFull04Tag_
edm::EDGetTokenT< edm::ValueMap< double > > isoFromDepsEcalFull04Tag_
Definition: ElectronMcFakeValidator.h:54
ElectronMcFakeValidator::h1_ele_dEtaEleCl_propOut_endcaps
MonitorElement * h1_ele_dEtaEleCl_propOut_endcaps
Definition: ElectronMcFakeValidator.h:338
ElectronMcFakeValidator::h1_ele_seed_dphi2pos_
MonitorElement * h1_ele_seed_dphi2pos_
Definition: ElectronMcFakeValidator.h:357
ElectronMcFakeValidator::h1_ele_HoE_all
MonitorElement * h1_ele_HoE_all
Definition: ElectronMcFakeValidator.h:169
pi
const Double_t pi
Definition: trackSplitPlot.h:36
ElectronMcFakeValidator::h1_scl_En_
MonitorElement * h1_scl_En_
Definition: ElectronMcFakeValidator.h:222
ElectronMcFakeValidator::h1_recTrackNum_
MonitorElement * h1_recTrackNum_
Definition: ElectronMcFakeValidator.h:148
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
ElectronMcFakeValidator::h2_ele_chargeVsPhi
MonitorElement * h2_ele_chargeVsPhi
Definition: ElectronMcFakeValidator.h:187
ElectronMcFakeValidator::inputFile_
std::string inputFile_
Definition: ElectronMcFakeValidator.h:71
ElectronMcFakeValidator::h1_ele_PinMnPout
MonitorElement * h1_ele_PinMnPout
Definition: ElectronMcFakeValidator.h:271
reco::GsfElectron::pfIsolationVariables
const PflowIsolationVariables & pfIsolationVariables() const
Definition: GsfElectron.h:721
ElectronMcFakeValidator::h1_ele_HoE_endcaps
MonitorElement * h1_ele_HoE_endcaps
Definition: ElectronMcFakeValidator.h:377
ElectronMcFakeValidator::track_nbin
int track_nbin
Definition: ElectronMcFakeValidator.h:136
ElectronMcFakeValidator::h1_ele_vertexPhi
MonitorElement * h1_ele_vertexPhi
Definition: ElectronMcFakeValidator.h:196
ParameterSet.h
ElectronMcFakeValidator::h1_ele_hcalDepth1OverEcalBc_endcaps
MonitorElement * h1_ele_hcalDepth1OverEcalBc_endcaps
Definition: ElectronMcFakeValidator.h:463
ElectronDqmAnalyzerBase::setBookPrefix
void setBookPrefix(const std::string &)
Definition: ElectronDqmAnalyzerBase.cc:30
ElectronMcFakeValidator::h2_ele_dPhiClVsEta_propOut
MonitorElement * h2_ele_dPhiClVsEta_propOut
Definition: ElectronMcFakeValidator.h:333
ElectronMcFakeValidator::h1_ele_PoPmatchingObject_endcaps
MonitorElement * h1_ele_PoPmatchingObject_endcaps
Definition: ElectronMcFakeValidator.h:210
ElectronMcFakeValidator::ele_nbin
int ele_nbin
Definition: ElectronMcFakeValidator.h:130
ElectronMcFakeValidator::h1_ele_chargedHadronRelativeIso_endcaps
MonitorElement * h1_ele_chargedHadronRelativeIso_endcaps
Definition: ElectronMcFakeValidator.h:432
ElectronMcFakeValidator::h1_ele_matchingObjectPhi_matched
MonitorElement * h1_ele_matchingObjectPhi_matched
Definition: ElectronMcFakeValidator.h:182
ElectronMcFakeValidator::p_nbin
int p_nbin
Definition: ElectronMcFakeValidator.h:79
ElectronMcFakeValidator::h1_recSeedNum_
MonitorElement * h1_recSeedNum_
Definition: ElectronMcFakeValidator.h:149
ElectronMcFakeValidator::h2_ele_vertexTIPVsPt
MonitorElement * h2_ele_vertexTIPVsPt
Definition: ElectronMcFakeValidator.h:203
ElectronMcFakeValidator::h1_scl_EoEmatchingObjectGolden_barrel
MonitorElement * h1_scl_EoEmatchingObjectGolden_barrel
Definition: ElectronMcFakeValidator.h:404
ElectronMcFakeValidator::seed_nbin
int seed_nbin
Definition: ElectronMcFakeValidator.h:139
ElectronMcFakeValidator::detamatch_nbin
int detamatch_nbin
Definition: ElectronMcFakeValidator.h:97
ElectronMcFakeValidator::eop_max
double eop_max
Definition: ElectronMcFakeValidator.h:114
ElectronMcFakeValidator::h1_ele_dEtaEleCl_propOut
MonitorElement * h1_ele_dEtaEleCl_propOut
Definition: ElectronMcFakeValidator.h:336
ElectronMcFakeValidator::h1_scl_full5x5_sigmaIetaIeta_barrel_
MonitorElement * h1_scl_full5x5_sigmaIetaIeta_barrel_
Definition: ElectronMcFakeValidator.h:236
ElectronMcFakeValidator::h1_ele_EeleOPout_barrel
MonitorElement * h1_ele_EeleOPout_barrel
Definition: ElectronMcFakeValidator.h:306
HLT_FULL_cff.flags
flags
Definition: HLT_FULL_cff.py:13168
ElectronMcFakeValidator::h1_ele_ecalRecHitSumEt_dr03
MonitorElement * h1_ele_ecalRecHitSumEt_dr03
Definition: ElectronMcFakeValidator.h:443
edm::Event
Definition: Event.h:73
ElectronMcFakeValidator::h1_ele_outerP
MonitorElement * h1_ele_outerP
Definition: ElectronMcFakeValidator.h:279
reco::GsfElectron::caloEnergy
float caloEnergy() const
Definition: GsfElectron.h:886
ElectronMcFakeValidator::h1_ele_HoE_fiducial
MonitorElement * h1_ele_HoE_fiducial
Definition: ElectronMcFakeValidator.h:380
ElectronMcFakeValidator::h2_scl_EtVsEta_
MonitorElement * h2_scl_EtVsEta_
Definition: ElectronMcFakeValidator.h:226
reco::GsfElectron::scSigmaIEtaIEta
float scSigmaIEtaIEta() const
Definition: GsfElectron.h:504
ElectronMcFakeValidator::h1_ele_dPhiCl_propOut_barrel
MonitorElement * h1_ele_dPhiCl_propOut_barrel
Definition: ElectronMcFakeValidator.h:331
ElectronMcFakeValidator::h1_ele_convDist_all
MonitorElement * h1_ele_convDist_all
Definition: ElectronMcFakeValidator.h:484
ElectronMcFakeValidator::h2_ele_PtinVsPtoutGolden_mean
MonitorElement * h2_ele_PtinVsPtoutGolden_mean
Definition: ElectronMcFakeValidator.h:402
ElectronMcFakeValidator::h1_ele_EseedOP_barrel
MonitorElement * h1_ele_EseedOP_barrel
Definition: ElectronMcFakeValidator.h:294
ElectronMcFakeValidator::opv_nbin
int opv_nbin
Definition: ElectronMcFakeValidator.h:127
ElectronMcFakeValidator::h1_ele_seed_mask_tec_
MonitorElement * h1_ele_seed_mask_tec_
Definition: ElectronMcFakeValidator.h:353
reco::GsfElectron::deltaPhiEleClusterTrackAtCalo
float deltaPhiEleClusterTrackAtCalo() const
Definition: GsfElectron.h:238
ElectronMcFakeValidator::h2_ele_seed_dphi2posVsPt_
MonitorElement * h2_ele_seed_dphi2posVsPt_
Definition: ElectronMcFakeValidator.h:359
ElectronMcFakeValidator::h1_ele_photonIso
MonitorElement * h1_ele_photonIso
Definition: ElectronMcFakeValidator.h:426
dttmaxenums::R
Definition: DTTMax.h:29
ElectronMcFakeValidator::h2_ele_PinVsPoutShowering_mode
MonitorElement * h2_ele_PinVsPoutShowering_mode
Definition: ElectronMcFakeValidator.h:397
ElectronMcFakeValidator::h1_ele_HoE_bc_barrel
MonitorElement * h1_ele_HoE_bc_barrel
Definition: ElectronMcFakeValidator.h:378
edm::InputTag
Definition: InputTag.h:15
ElectronMcFakeValidator::h1_scl_full5x5_sigmaIetaIeta_
MonitorElement * h1_scl_full5x5_sigmaIetaIeta_
Definition: ElectronMcFakeValidator.h:235
reco::GsfElectron::dr03EcalRecHitSumEt
float dr03EcalRecHitSumEt() const
Definition: GsfElectron.h:567
ElectronMcFakeValidator::beamSpotTag_
edm::EDGetTokenT< reco::BeamSpot > beamSpotTag_
Definition: ElectronMcFakeValidator.h:48
ElectronMcFakeValidator::fhits_nbin
int fhits_nbin
Definition: ElectronMcFakeValidator.h:86
ElectronMcFakeValidator::h1_ele_foundHits_endcaps
MonitorElement * h1_ele_foundHits_endcaps
Definition: ElectronMcFakeValidator.h:254
reco::GsfElectron::PflowIsolationVariables::sumChargedHadronPt
float sumChargedHadronPt
sum-pt of charged Hadron // old float chargedHadronIso ;
Definition: GsfElectron.h:671
ElectronSeed.h
ElectronMcFakeValidator::h1_ele_hcalDepth2OverEcalBc
MonitorElement * h1_ele_hcalDepth2OverEcalBc
Definition: ElectronMcFakeValidator.h:464
ElectronMcFakeValidator::h1_ele_chi2_barrel
MonitorElement * h1_ele_chi2_barrel
Definition: ElectronMcFakeValidator.h:265
ElectronMcFakeValidator::h1_ele_seed_dphi2_
MonitorElement * h1_ele_seed_dphi2_
Definition: ElectronMcFakeValidator.h:354
ElectronMcFakeValidator::h1_ele_hcalTowerSumEtBc_dr03_depth1_barrel
MonitorElement * h1_ele_hcalTowerSumEtBc_dr03_depth1_barrel
Definition: ElectronMcFakeValidator.h:469
ElectronMcFakeValidator::h2_ele_EeleOPoutVsE
MonitorElement * h2_ele_EeleOPoutVsE
Definition: ElectronMcFakeValidator.h:310
ElectronMcFakeValidator::h1_ele_vertexAbsEta
MonitorElement * h1_ele_vertexAbsEta
Definition: ElectronMcFakeValidator.h:195
ElectronMcFakeValidator::h2_ele_PtinVsPtoutGolden_mode
MonitorElement * h2_ele_PtinVsPtoutGolden_mode
Definition: ElectronMcFakeValidator.h:400
ElectronMcFakeValidator::h1_ele_convDist
MonitorElement * h1_ele_convDist
Definition: ElectronMcFakeValidator.h:483
ElectronMcFakeValidator::h1_ele_outerPt_mode
MonitorElement * h1_ele_outerPt_mode
Definition: ElectronMcFakeValidator.h:283
ElectronMcFakeValidator::h1_ele_convDcot
MonitorElement * h1_ele_convDcot
Definition: ElectronMcFakeValidator.h:485
ElectronMcFakeValidator::h2_ele_EseedOPVsEta
MonitorElement * h2_ele_EseedOPVsEta
Definition: ElectronMcFakeValidator.h:296
Generated for CMSSW Reference Manual by doxygen 1.8.16