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