CMS 3D CMS Logo
GsfElectronDataAnalyzer.cc
Go to the documentation of this file.
1 
2 // -*- C++ -*-
3 //
4 // Package: RecoEgamma/Examples
5 // Class: GsfElectronDataAnalyzer
6 //
14 //
15 // Original Author: Ursula Berthon, Claude Charlot
16 // Created: Mon Mar 27 13:22:06 CEST 2006
17 //
18 //
19 
20 #include "DataFormats/BeamSpot/interface/BeamSpot.h"
21 #include "DataFormats/Common/interface/Handle.h"
22 #include "DataFormats/Common/interface/TriggerResults.h"
23 #include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
24 #include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"
25 #include "DataFormats/EgammaReco/interface/ElectronSeed.h"
26 #include "DataFormats/EgammaReco/interface/ElectronSeedFwd.h"
27 #include "DataFormats/EgammaReco/interface/SuperClusterFwd.h"
28 #include "DataFormats/GsfTrackReco/interface/GsfTrack.h"
29 #include "FWCore/Common/interface/TriggerNames.h"
30 #include "FWCore/Framework/interface/EDAnalyzer.h"
31 #include "FWCore/Framework/interface/ESHandle.h"
32 #include "FWCore/Framework/interface/Event.h"
33 #include "FWCore/Framework/interface/EventSetup.h"
34 #include "FWCore/Framework/interface/one/EDAnalyzer.h"
35 #include "FWCore/MessageLogger/interface/MessageLogger.h"
36 #include "FWCore/ParameterSet/interface/ParameterSet.h"
37 #include "FWCore/Utilities/interface/InputTag.h"
38 #include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
39 #include "TrackingTools/TrajectoryState/interface/TrajectoryStateTransform.h"
40 
41 #include "CLHEP/Units/GlobalPhysicalConstants.h"
42 
43 #include "TMath.h"
44 #include "TFile.h"
45 #include "TH1F.h"
46 #include "TH1I.h"
47 #include "TH2F.h"
48 #include "TProfile.h"
49 #include "TTree.h"
50 
51 #include <iostream>
52 
53 class GsfElectronDataAnalyzer : public edm::one::EDAnalyzer<> {
54 public:
55  explicit GsfElectronDataAnalyzer(const edm::ParameterSet &conf);
56 
57  ~GsfElectronDataAnalyzer() override;
58 
59  void beginJob() override;
60  void endJob() override;
61  void analyze(const edm::Event &e, const edm::EventSetup &c) override;
62 
63 private:
64  bool trigger(const edm::Event &e);
65 
66  unsigned int nEvents_;
67  unsigned int nAfterTrigger_;
68 
69  edm::InputTag triggerResults_;
70  edm::InputTag beamSpot_;
71 
72  std::vector<std::string> HLTPathsByName_;
73  std::vector<unsigned int> HLTPathsByIndex_;
74 
75  std::string outputFile_;
76  edm::InputTag electronCollection_;
77  edm::InputTag matchingObjectCollection_;
78  std::string matchingCondition_;
79  std::string type_;
80  bool readAOD_;
81  // matching
82  double maxPtMatchingObject_;
83  double maxAbsEtaMatchingObject_;
84  double deltaR_;
85 
86  TrajectoryStateTransform transformer_;
87  edm::ESHandle<TrackerGeometry> pDD;
88  edm::ESHandle<MagneticField> theMagField;
89 
90  TFile *histfile_;
91  TTree *tree_;
92  float mcEnergy[10], mcEta[10], mcPhi[10], mcPt[10], mcQ[10];
93  float superclusterEnergy[10], superclusterEta[10], superclusterPhi[10], superclusterEt[10];
94  float seedMomentum[10], seedEta[10], seedPhi[10], seedPt[10], seedQ[10];
95 
96  TH1F *h_matchingObjectNum;
97 
98  TH1F *h_matchingObjectEta;
99  TH1F *h_matchingObjectAbsEta;
100  TH1F *h_matchingObjectP;
101  TH1F *h_matchingObjectPt;
102  TH1F *h_matchingObjectPhi;
103  TH1F *h_matchingObjectZ;
104 
105  TH1F *h_ele_matchingObjectEta_matched;
106  TH1F *h_ele_matchingObjectAbsEta_matched;
107  TH1F *h_ele_matchingObjectPt_matched;
108  TH1F *h_ele_matchingObjectPhi_matched;
109  TH1F *h_ele_matchingObjectZ_matched;
110 
111  TH1F *h_ele_mee_all;
112  TH1F *h_ele_mee_os;
113  TH1F *h_ele_mee_os_ebeb;
114  TH1F *h_ele_mee_os_ebee;
115  TH1F *h_ele_mee_os_eeee;
116  TH1F *h_ele_mee_os_gg;
117  TH1F *h_ele_mee_os_gb;
118  TH1F *h_ele_mee_os_bb;
119 
120  TH2F *h_ele_E2mnE1vsMee_all;
121  TH2F *h_ele_E2mnE1vsMee_egeg_all;
122 
123  TH1F *h_ele_charge;
124  TH2F *h_ele_chargeVsEta;
125  TH2F *h_ele_chargeVsPhi;
126  TH2F *h_ele_chargeVsPt;
127  TH1F *h_ele_vertexP;
128  TH1F *h_ele_vertexPt;
129  TH1F *h_ele_Et;
130  TH2F *h_ele_vertexPtVsEta;
131  TH2F *h_ele_vertexPtVsPhi;
132  TH1F *h_ele_vertexPt_5100;
133  TH1F *h_ele_vertexEta;
134  TH2F *h_ele_vertexEtaVsPhi;
135  TH1F *h_ele_vertexAbsEta;
136  TH1F *h_ele_vertexPhi;
137  TH1F *h_ele_vertexX;
138  TH1F *h_ele_vertexY;
139  TH1F *h_ele_vertexZ;
140  TH1F *h_ele_vertexTIP;
141  TH2F *h_ele_vertexTIPVsEta;
142  TH2F *h_ele_vertexTIPVsPhi;
143  TH2F *h_ele_vertexTIPVsPt;
144 
145  TH1F *histNum_;
146 
147  TH1F *histSclEn_;
148  TH1F *histSclEoEmatchingObject_barrel_matched;
149  TH1F *histSclEoEmatchingObject_endcaps_matched;
150  TH1F *histSclEoEmatchingObject_barrel_new_matched;
151  TH1F *histSclEoEmatchingObject_endcaps_new_matched;
152  TH1F *histSclEt_;
153  TH2F *histSclEtVsEta_;
154  TH2F *histSclEtVsPhi_;
155  TH2F *histSclEtaVsPhi_;
156  TH1F *histSclEta_;
157  TH1F *histSclPhi_;
158  TH1F *histSclSigEtaEta_;
159  TH1F *histSclSigEtaEta_barrel_;
160  TH1F *histSclSigEtaEta_endcaps_;
161  TH1F *histSclSigIEtaIEta_;
162  TH1F *histSclSigIEtaIEta_barrel_;
163  TH1F *histSclSigIEtaIEta_endcaps_;
164  TH1F *histSclE1x5_;
165  TH1F *histSclE1x5_barrel_;
166  TH1F *histSclE1x5_endcaps_;
167  TH1F *histSclE2x5max_;
168  TH1F *histSclE2x5max_barrel_;
169  TH1F *histSclE2x5max_endcaps_;
170  TH1F *histSclE5x5_;
171  TH1F *histSclE5x5_barrel_;
172  TH1F *histSclE5x5_endcaps_;
173 
174  TH1F *h_ele_ambiguousTracks;
175  TH2F *h_ele_ambiguousTracksVsEta;
176  TH2F *h_ele_ambiguousTracksVsPhi;
177  TH2F *h_ele_ambiguousTracksVsPt;
178  TH1F *h_ele_foundHits;
179  TH1F *h_ele_foundHits_barrel;
180  TH1F *h_ele_foundHits_endcaps;
181  TH2F *h_ele_foundHitsVsEta;
182  TH2F *h_ele_foundHitsVsPhi;
183  TH2F *h_ele_foundHitsVsPt;
184  TH1F *h_ele_lostHits;
185  TH1F *h_ele_lostHits_barrel;
186  TH1F *h_ele_lostHits_endcaps;
187  TH2F *h_ele_lostHitsVsEta;
188  TH2F *h_ele_lostHitsVsPhi;
189  TH2F *h_ele_lostHitsVsPt;
190  TH1F *h_ele_chi2;
191  TH1F *h_ele_chi2_barrel_;
192  TH1F *h_ele_chi2_endcaps_;
193  TH2F *h_ele_chi2VsEta;
194  TH2F *h_ele_chi2VsPhi;
195  TH2F *h_ele_chi2VsPt;
196 
197  TH1F *h_ele_PtoPtmatchingObject_matched;
198  TH1F *h_ele_PoPmatchingObject_matched;
199  TH2F *h_ele_PoPmatchingObjectVsEta_matched;
200  TH2F *h_ele_PoPmatchingObjectVsPhi_matched;
201  TH2F *h_ele_PoPmatchingObjectVsPt_matched;
202  TH1F *h_ele_PoPmatchingObject_barrel_matched;
203  TH1F *h_ele_PoPmatchingObject_endcaps_matched;
204  TH1F *h_ele_PtoPtmatchingObject_barrel_matched;
205  TH1F *h_ele_PtoPtmatchingObject_endcaps_matched;
206  TH1F *h_ele_EtaMnEtamatchingObject_matched;
207  TH2F *h_ele_EtaMnEtamatchingObjectVsEta_matched;
208  TH2F *h_ele_EtaMnEtamatchingObjectVsPhi_matched;
209  TH2F *h_ele_EtaMnEtamatchingObjectVsPt_matched;
210  TH1F *h_ele_PhiMnPhimatchingObject_matched;
211  TH1F *h_ele_PhiMnPhimatchingObject2_matched;
212  TH2F *h_ele_PhiMnPhimatchingObjectVsEta_matched;
213  TH2F *h_ele_PhiMnPhimatchingObjectVsPhi_matched;
214  TH2F *h_ele_PhiMnPhimatchingObjectVsPt_matched;
215  TH1F *h_ele_PinMnPout;
216  TH1F *h_ele_PinMnPout_mode;
217  TH2F *h_ele_PinMnPoutVsEta_mode;
218  TH2F *h_ele_PinMnPoutVsPhi_mode;
219  TH2F *h_ele_PinMnPoutVsPt_mode;
220  TH2F *h_ele_PinMnPoutVsE_mode;
221  TH2F *h_ele_PinMnPoutVsChi2_mode;
222 
223  TH1F *h_ele_outerP;
224  TH1F *h_ele_outerP_mode;
225  TH2F *h_ele_outerPVsEta_mode;
226  TH1F *h_ele_outerPt;
227  TH1F *h_ele_outerPt_mode;
228  TH2F *h_ele_outerPtVsEta_mode;
229  TH2F *h_ele_outerPtVsPhi_mode;
230  TH2F *h_ele_outerPtVsPt_mode;
231  TH1F *h_ele_EoP;
232  TH1F *h_ele_EoP_barrel;
233  TH1F *h_ele_EoP_endcaps;
234  TH2F *h_ele_EoPVsEta;
235  TH2F *h_ele_EoPVsPhi;
236  TH2F *h_ele_EoPVsE;
237  TH1F *h_ele_EseedOP;
238  TH1F *h_ele_EseedOP_barrel;
239  TH1F *h_ele_EseedOP_endcaps;
240  TH2F *h_ele_EseedOPVsEta;
241  TH2F *h_ele_EseedOPVsPhi;
242  TH2F *h_ele_EseedOPVsE;
243  TH1F *h_ele_EoPout;
244  TH1F *h_ele_EoPout_barrel;
245  TH1F *h_ele_EoPout_endcaps;
246  TH2F *h_ele_EoPoutVsEta;
247  TH2F *h_ele_EoPoutVsPhi;
248  TH2F *h_ele_EoPoutVsE;
249  TH1F *h_ele_EeleOPout;
250  TH1F *h_ele_EeleOPout_barrel;
251  TH1F *h_ele_EeleOPout_endcaps;
252  TH2F *h_ele_EeleOPoutVsEta;
253  TH2F *h_ele_EeleOPoutVsPhi;
254  TH2F *h_ele_EeleOPoutVsE;
255 
256  TH1F *h_ele_dEtaSc_propVtx;
257  TH1F *h_ele_dEtaSc_propVtx_barrel;
258  TH1F *h_ele_dEtaSc_propVtx_endcaps;
259  TH2F *h_ele_dEtaScVsEta_propVtx;
260  TH2F *h_ele_dEtaScVsPhi_propVtx;
261  TH2F *h_ele_dEtaScVsPt_propVtx;
262  TH1F *h_ele_dPhiSc_propVtx;
263  TH1F *h_ele_dPhiSc_propVtx_barrel;
264  TH1F *h_ele_dPhiSc_propVtx_endcaps;
265  TH2F *h_ele_dPhiScVsEta_propVtx;
266  TH2F *h_ele_dPhiScVsPhi_propVtx;
267  TH2F *h_ele_dPhiScVsPt_propVtx;
268  TH1F *h_ele_dEtaCl_propOut;
269  TH1F *h_ele_dEtaCl_propOut_barrel;
270  TH1F *h_ele_dEtaCl_propOut_endcaps;
271  TH2F *h_ele_dEtaClVsEta_propOut;
272  TH2F *h_ele_dEtaClVsPhi_propOut;
273  TH2F *h_ele_dEtaClVsPt_propOut;
274  TH1F *h_ele_dPhiCl_propOut;
275  TH1F *h_ele_dPhiCl_propOut_barrel;
276  TH1F *h_ele_dPhiCl_propOut_endcaps;
277  TH2F *h_ele_dPhiClVsEta_propOut;
278  TH2F *h_ele_dPhiClVsPhi_propOut;
279  TH2F *h_ele_dPhiClVsPt_propOut;
280  TH1F *h_ele_dEtaEleCl_propOut;
281  TH1F *h_ele_dEtaEleCl_propOut_barrel;
282  TH1F *h_ele_dEtaEleCl_propOut_endcaps;
283  TH2F *h_ele_dEtaEleClVsEta_propOut;
284  TH2F *h_ele_dEtaEleClVsPhi_propOut;
285  TH2F *h_ele_dEtaEleClVsPt_propOut;
286  TH1F *h_ele_dPhiEleCl_propOut;
287  TH1F *h_ele_dPhiEleCl_propOut_barrel;
288  TH1F *h_ele_dPhiEleCl_propOut_endcaps;
289  TH2F *h_ele_dPhiEleClVsEta_propOut;
290  TH2F *h_ele_dPhiEleClVsPhi_propOut;
291  TH2F *h_ele_dPhiEleClVsPt_propOut;
292 
293  TH1F *h_ele_seed_dphi2_;
294  TH2F *h_ele_seed_dphi2VsEta_;
295  TH2F *h_ele_seed_dphi2VsPt_;
296  TH1F *h_ele_seed_drz2_;
297  TH2F *h_ele_seed_drz2VsEta_;
298  TH2F *h_ele_seed_drz2VsPt_;
299  TH1F *h_ele_seed_subdet2_;
300 
301  TH1F *h_ele_classes;
302  TH1F *h_ele_eta;
303  TH1F *h_ele_eta_golden;
304  TH1F *h_ele_eta_bbrem;
305  TH1F *h_ele_eta_narrow;
306  TH1F *h_ele_eta_shower;
307 
308  TH1F *h_ele_HoE;
309  TH1F *h_ele_HoE_barrel;
310  TH1F *h_ele_HoE_endcaps;
311  TH1F *h_ele_HoE_fiducial;
312  TH2F *h_ele_HoEVsEta;
313  TH2F *h_ele_HoEVsPhi;
314  TH2F *h_ele_HoEVsE;
315 
316  TH1F *h_ele_fbrem;
317  TProfile *h_ele_fbremVsEta_mode;
318  TProfile *h_ele_fbremVsEta_mean;
319 
320  TH2F *h_ele_PinVsPoutGolden_mode;
321  TH2F *h_ele_PinVsPoutShowering_mode;
322  TH2F *h_ele_PinVsPoutGolden_mean;
323  TH2F *h_ele_PinVsPoutShowering_mean;
324  TH2F *h_ele_PtinVsPtoutGolden_mode;
325  TH2F *h_ele_PtinVsPtoutShowering_mode;
326  TH2F *h_ele_PtinVsPtoutGolden_mean;
327  TH2F *h_ele_PtinVsPtoutShowering_mean;
328  TH1F *histSclEoEmatchingObjectGolden_barrel;
329  TH1F *histSclEoEmatchingObjectGolden_endcaps;
330  TH1F *histSclEoEmatchingObjectShowering_barrel;
331  TH1F *histSclEoEmatchingObjectShowering_endcaps;
332 
333  TH1F *h_ele_mva;
334  TH1F *h_ele_provenance;
335 
336  TH1F *h_ele_tkSumPt_dr03;
337  TH1F *h_ele_ecalRecHitSumEt_dr03;
338  TH1F *h_ele_hcalDepth1TowerSumEt_dr03;
339  TH1F *h_ele_hcalDepth2TowerSumEt_dr03;
340  TH1F *h_ele_tkSumPt_dr04;
341  TH1F *h_ele_ecalRecHitSumEt_dr04;
342  TH1F *h_ele_hcalDepth1TowerSumEt_dr04;
343  TH1F *h_ele_hcalDepth2TowerSumEt_dr04;
344 
345  // electron selection
346  double minEt_;
347  double minPt_;
348  double maxAbsEta_;
349  bool isEB_;
350  bool isEE_;
351  bool isNotEBEEGap_;
352  bool isEcalDriven_;
353  bool isTrackerDriven_;
354  double eOverPMinBarrel_;
355  double eOverPMaxBarrel_;
356  double eOverPMinEndcaps_;
357  double eOverPMaxEndcaps_;
358  double dEtaMinBarrel_;
359  double dEtaMaxBarrel_;
360  double dEtaMinEndcaps_;
361  double dEtaMaxEndcaps_;
362  double dPhiMinBarrel_;
363  double dPhiMaxBarrel_;
364  double dPhiMinEndcaps_;
365  double dPhiMaxEndcaps_;
366  double sigIetaIetaMinBarrel_;
367  double sigIetaIetaMaxBarrel_;
368  double sigIetaIetaMinEndcaps_;
369  double sigIetaIetaMaxEndcaps_;
370  double hadronicOverEmMaxBarrel_;
371  double hadronicOverEmMaxEndcaps_;
372  double mvaMin_;
373  double tipMaxBarrel_;
374  double tipMaxEndcaps_;
375  double tkIso03Max_;
376  double hcalIso03Depth1MaxBarrel_;
377  double hcalIso03Depth1MaxEndcaps_;
378  double hcalIso03Depth2MaxEndcaps_;
379  double ecalIso03MaxBarrel_;
380  double ecalIso03MaxEndcaps_;
381 
382  // histos limits and binning
383  double etamin;
384  double etamax;
385  double phimin;
386  double phimax;
387  double ptmax;
388  double pmax;
389  double eopmax;
390  double eopmaxsht;
391  double detamin;
392  double detamax;
393  double dphimin;
394  double dphimax;
395  double detamatchmin;
396  double detamatchmax;
397  double dphimatchmin;
398  double dphimatchmax;
399  double fhitsmax;
400  double lhitsmax;
401  double poptruemin;
402  double poptruemax;
403  double meemin;
404  double meemax;
405  double hoemin;
406  double hoemax;
407  int nbineta;
408  int nbinp;
409  int nbinpt;
410  int nbinpteff;
411  int nbinphi;
412  int nbinp2D;
413  int nbinpt2D;
414  int nbineta2D;
415  int nbinphi2D;
416  int nbineop;
417  int nbineop2D;
418  int nbinfhits;
419  int nbinlhits;
420  int nbinxyz;
421  int nbindeta;
422  int nbindphi;
423  int nbindetamatch;
424  int nbindphimatch;
425  int nbindetamatch2D;
426  int nbindphimatch2D;
427  int nbinpoptrue;
428  int nbinmee;
429  int nbinhoe;
430 };
431 
432 #include "FWCore/Framework/interface/MakerMacros.h"
433 DEFINE_FWK_MODULE(GsfElectronDataAnalyzer);
434 
435 using namespace reco;
436 
437 GsfElectronDataAnalyzer::GsfElectronDataAnalyzer(const edm::ParameterSet &conf)
438  : beamSpot_(conf.getParameter<edm::InputTag>("beamSpot")) {
439  outputFile_ = conf.getParameter<std::string>("outputFile");
440  histfile_ = new TFile(outputFile_.c_str(), "RECREATE");
441  electronCollection_ = conf.getParameter<edm::InputTag>("electronCollection");
442  readAOD_ = conf.getParameter<bool>("readAOD");
443 
444  matchingObjectCollection_ = conf.getParameter<edm::InputTag>("matchingObjectCollection");
445  matchingCondition_ = conf.getParameter<std::string>("matchingCondition");
446  // currently only one possible matching conditions
447  assert(matchingCondition_ == "Cone");
448  maxPtMatchingObject_ = conf.getParameter<double>("MaxPtMatchingObject");
449  maxAbsEtaMatchingObject_ = conf.getParameter<double>("MaxAbsEtaMatchingObject");
450  deltaR_ = conf.getParameter<double>("DeltaR");
451 
452  triggerResults_ = conf.getParameter<edm::InputTag>("triggerResults");
453  HLTPathsByName_ = conf.getParameter<std::vector<std::string> >("hltPaths");
454  HLTPathsByIndex_.resize(HLTPathsByName_.size());
455 
456  minEt_ = conf.getParameter<double>("MinEt");
457  minPt_ = conf.getParameter<double>("MinPt");
458  maxAbsEta_ = conf.getParameter<double>("MaxAbsEta");
459  isEB_ = conf.getParameter<bool>("SelectEB");
460  isEE_ = conf.getParameter<bool>("SelectEE");
461  isNotEBEEGap_ = conf.getParameter<bool>("SelectNotEBEEGap");
462  isEcalDriven_ = conf.getParameter<bool>("SelectEcalDriven");
463  isTrackerDriven_ = conf.getParameter<bool>("SelectTrackerDriven");
464  eOverPMinBarrel_ = conf.getParameter<double>("MinEOverPBarrel");
465  eOverPMaxBarrel_ = conf.getParameter<double>("MaxEOverPBarrel");
466  eOverPMinEndcaps_ = conf.getParameter<double>("MinEOverPEndcaps");
467  eOverPMaxEndcaps_ = conf.getParameter<double>("MaxEOverPEndcaps");
468  dEtaMinBarrel_ = conf.getParameter<double>("MinDetaBarrel");
469  dEtaMaxBarrel_ = conf.getParameter<double>("MaxDetaBarrel");
470  dEtaMinEndcaps_ = conf.getParameter<double>("MinDetaEndcaps");
471  dEtaMaxEndcaps_ = conf.getParameter<double>("MaxDetaEndcaps");
472  dPhiMinBarrel_ = conf.getParameter<double>("MinDphiBarrel");
473  dPhiMaxBarrel_ = conf.getParameter<double>("MaxDphiBarrel");
474  dPhiMinEndcaps_ = conf.getParameter<double>("MinDphiEndcaps");
475  dPhiMaxEndcaps_ = conf.getParameter<double>("MaxDphiEndcaps");
476  sigIetaIetaMinBarrel_ = conf.getParameter<double>("MinSigIetaIetaBarrel");
477  sigIetaIetaMaxBarrel_ = conf.getParameter<double>("MaxSigIetaIetaBarrel");
478  sigIetaIetaMinEndcaps_ = conf.getParameter<double>("MinSigIetaIetaEndcaps");
479  sigIetaIetaMaxEndcaps_ = conf.getParameter<double>("MaxSigIetaIetaEndcaps");
480  hadronicOverEmMaxBarrel_ = conf.getParameter<double>("MaxHoEBarrel");
481  hadronicOverEmMaxEndcaps_ = conf.getParameter<double>("MaxHoEEndcaps");
482  mvaMin_ = conf.getParameter<double>("MinMVA");
483  tipMaxBarrel_ = conf.getParameter<double>("MaxTipBarrel");
484  tipMaxEndcaps_ = conf.getParameter<double>("MaxTipEndcaps");
485  tkIso03Max_ = conf.getParameter<double>("MaxTkIso03");
486  hcalIso03Depth1MaxBarrel_ = conf.getParameter<double>("MaxHcalIso03Depth1Barrel");
487  hcalIso03Depth1MaxEndcaps_ = conf.getParameter<double>("MaxHcalIso03Depth1Endcaps");
488  hcalIso03Depth2MaxEndcaps_ = conf.getParameter<double>("MaxHcalIso03Depth2Endcaps");
489  ecalIso03MaxBarrel_ = conf.getParameter<double>("MaxEcalIso03Barrel");
490  ecalIso03MaxEndcaps_ = conf.getParameter<double>("MaxEcalIso03Endcaps");
491 
492  edm::ParameterSet pset = conf.getParameter<edm::ParameterSet>("HistosConfigurationData");
493 
494  etamin = pset.getParameter<double>("Etamin");
495  etamax = pset.getParameter<double>("Etamax");
496  phimin = pset.getParameter<double>("Phimin");
497  phimax = pset.getParameter<double>("Phimax");
498  ptmax = pset.getParameter<double>("Ptmax");
499  pmax = pset.getParameter<double>("Pmax");
500  eopmax = pset.getParameter<double>("Eopmax");
501  eopmaxsht = pset.getParameter<double>("Eopmaxsht");
502  detamin = pset.getParameter<double>("Detamin");
503  detamax = pset.getParameter<double>("Detamax");
504  dphimin = pset.getParameter<double>("Dphimin");
505  dphimax = pset.getParameter<double>("Dphimax");
506  detamatchmin = pset.getParameter<double>("Detamatchmin");
507  detamatchmax = pset.getParameter<double>("Detamatchmax");
508  dphimatchmin = pset.getParameter<double>("Dphimatchmin");
509  dphimatchmax = pset.getParameter<double>("Dphimatchmax");
510  fhitsmax = pset.getParameter<double>("Fhitsmax");
511  lhitsmax = pset.getParameter<double>("Lhitsmax");
512  nbineta = pset.getParameter<int>("Nbineta");
513  nbineta2D = pset.getParameter<int>("Nbineta2D");
514  nbinp = pset.getParameter<int>("Nbinp");
515  nbinpt = pset.getParameter<int>("Nbinpt");
516  nbinp2D = pset.getParameter<int>("Nbinp2D");
517  nbinpt2D = pset.getParameter<int>("Nbinpt2D");
518  nbinpteff = pset.getParameter<int>("Nbinpteff");
519  nbinphi = pset.getParameter<int>("Nbinphi");
520  nbinphi2D = pset.getParameter<int>("Nbinphi2D");
521  nbineop = pset.getParameter<int>("Nbineop");
522  nbineop2D = pset.getParameter<int>("Nbineop2D");
523  nbinfhits = pset.getParameter<int>("Nbinfhits");
524  nbinlhits = pset.getParameter<int>("Nbinlhits");
525  nbinxyz = pset.getParameter<int>("Nbinxyz");
526  nbindeta = pset.getParameter<int>("Nbindeta");
527  nbindphi = pset.getParameter<int>("Nbindphi");
528  nbindetamatch = pset.getParameter<int>("Nbindetamatch");
529  nbindphimatch = pset.getParameter<int>("Nbindphimatch");
530  nbindetamatch2D = pset.getParameter<int>("Nbindetamatch2D");
531  nbindphimatch2D = pset.getParameter<int>("Nbindphimatch2D");
532  nbinpoptrue = pset.getParameter<int>("Nbinpoptrue");
533  poptruemin = pset.getParameter<double>("Poptruemin");
534  poptruemax = pset.getParameter<double>("Poptruemax");
535  nbinmee = pset.getParameter<int>("Nbinmee");
536  meemin = pset.getParameter<double>("Meemin");
537  meemax = pset.getParameter<double>("Meemax");
538  nbinhoe = pset.getParameter<int>("Nbinhoe");
539  hoemin = pset.getParameter<double>("Hoemin");
540  hoemax = pset.getParameter<double>("Hoemax");
541 }
542 
543 GsfElectronDataAnalyzer::~GsfElectronDataAnalyzer() {
544  // do anything here that needs to be done at desctruction time
545  // (e.g. close files, deallocate resources etc.)
546  histfile_->Write();
547  histfile_->Close();
548 }
549 
550 void GsfElectronDataAnalyzer::beginJob() {
551  histfile_->cd();
552 
553  nEvents_ = 0;
554  nAfterTrigger_ = 0;
555 
556  // matching object
557  std::string::size_type locSC = matchingObjectCollection_.label().find("SuperCluster", 0);
558  std::string type_;
559  if (locSC != std::string::npos) {
560  std::cout << "Matching objects are SuperClusters " << std::endl;
561  type_ = "SC";
562  } else {
563  std::cout << "Didn't recognize input matching objects!! " << std::endl;
564  }
565 
566  std::string htitle, hlabel;
567  hlabel = "h_" + type_ + "Num";
568  htitle = "# " + type_ + "s";
569  h_matchingObjectNum = new TH1F(hlabel.c_str(), htitle.c_str(), nbinfhits, 0., fhitsmax);
570 
571  // rec event
572 
573  histNum_ = new TH1F("h_recEleNum", "# rec electrons", 20, 0., 20.);
574 
575  // matching object distributions
576  hlabel = "h_" + type_ + "_eta";
577  htitle = type_ + " #eta";
578  h_matchingObjectEta = new TH1F(hlabel.c_str(), htitle.c_str(), nbineta, etamin, etamax);
579  hlabel = "h_" + type_ + "_abseta";
580  htitle = type_ + " |#eta|";
581  h_matchingObjectAbsEta = new TH1F(hlabel.c_str(), htitle.c_str(), nbineta / 2, 0., etamax);
582  hlabel = "h_" + type_ + "_P";
583  htitle = type_ + " p";
584  h_matchingObjectP = new TH1F(hlabel.c_str(), htitle.c_str(), nbinp, 0., pmax);
585  hlabel = "h_" + type_ + "_Pt";
586  htitle = type_ + " pt";
587  h_matchingObjectPt = new TH1F(hlabel.c_str(), htitle.c_str(), nbinpteff, 5., ptmax);
588  hlabel = "h_" + type_ + "_phi";
589  htitle = type_ + " phi";
590  h_matchingObjectPhi = new TH1F(hlabel.c_str(), htitle.c_str(), nbinphi, phimin, phimax);
591  hlabel = "h_" + type_ + "_z";
592  htitle = type_ + " z";
593  h_matchingObjectZ = new TH1F(hlabel.c_str(), htitle.c_str(), nbinxyz, -25, 25);
594 
595  // mee
596  h_ele_mee_all =
597  new TH1F("h_ele_mee_all", "ele pairs invariant mass, all charge combinations", nbinmee, meemin, meemax);
598  h_ele_mee_os = new TH1F("h_ele_mee_os", "ele pairs invariant mass, opposite sign", nbinmee, meemin, meemax);
599  h_ele_mee_os_ebeb =
600  new TH1F("h_ele_mee_os_ebeb", "ele pairs invariant mass, opp. sign, EB-EB", nbinmee, meemin, meemax);
601  h_ele_mee_os_ebeb->Sumw2();
602  h_ele_mee_os_ebee =
603  new TH1F("h_ele_mee_os_ebee", "ele pairs invariant mass, opp. sign, EB-EE", nbinmee, meemin, meemax);
604  h_ele_mee_os_ebee->Sumw2();
605  h_ele_mee_os_eeee =
606  new TH1F("h_ele_mee_os_eeee", "ele pairs invariant mass, opp. sign, EE-EE", nbinmee, meemin, meemax);
607  h_ele_mee_os_eeee->Sumw2();
608  h_ele_mee_os_gg =
609  new TH1F("h_ele_mee_os_gg", "ele pairs invariant mass, opp. sign, good-good", nbinmee, meemin, meemax);
610  h_ele_mee_os_gg->Sumw2();
611  h_ele_mee_os_gb =
612  new TH1F("h_ele_mee_os_gb", "ele pairs invariant mass, opp. sign, good-bad", nbinmee, meemin, meemax);
613  h_ele_mee_os_gb->Sumw2();
614  h_ele_mee_os_bb =
615  new TH1F("h_ele_mee_os_bb", "ele pairs invariant mass, opp. sign, bad-bad", nbinmee, meemin, meemax);
616  h_ele_mee_os_bb->Sumw2();
617 
618  // duplicates
619  h_ele_E2mnE1vsMee_all = new TH2F("h_ele_E2mnE1vsMee_all",
620  "E2 - E1 vs ele pairs invariant mass, all electrons",
621  nbinmee,
622  meemin,
623  meemax,
624  100,
625  -50.,
626  50.);
627  h_ele_E2mnE1vsMee_egeg_all = new TH2F("h_ele_E2mnE1vsMee_egeg_all",
628  "E2 - E1 vs ele pairs invariant mass, ecal driven pairs, all electrons",
629  nbinmee,
630  meemin,
631  meemax,
632  100,
633  -50.,
634  50.);
635 
636  // recoed and matched electrons
637  h_ele_charge = new TH1F("h_ele_charge", "ele charge", 5, -2., 2.);
638  h_ele_chargeVsEta = new TH2F("h_ele_chargeVsEta", "ele charge vs eta", nbineta2D, etamin, etamax, 5, -2., 2.);
639  h_ele_chargeVsPhi = new TH2F("h_ele_chargeVsPhi", "ele charge vs phi", nbinphi2D, phimin, phimax, 5, -2., 2.);
640  h_ele_chargeVsPt = new TH2F("h_ele_chargeVsPt", "ele charge vs pt", nbinpt, 0., 100., 5, -2., 2.);
641  h_ele_vertexP = new TH1F("h_ele_vertexP", "ele momentum", nbinp, 0., pmax);
642  h_ele_vertexPt = new TH1F("h_ele_vertexPt", "ele transverse momentum", nbinpt, 0., ptmax);
643  h_ele_Et = new TH1F("h_ele_Et", "ele SC transverse energy", nbinpt, 0., ptmax);
644  h_ele_vertexPtVsEta =
645  new TH2F("h_ele_vertexPtVsEta", "ele transverse momentum vs eta", nbineta2D, etamin, etamax, nbinpt2D, 0., ptmax);
646  h_ele_vertexPtVsPhi =
647  new TH2F("h_ele_vertexPtVsPhi", "ele transverse momentum vs phi", nbinphi2D, phimin, phimax, nbinpt2D, 0., ptmax);
648  h_ele_matchingObjectPt_matched =
649  new TH1F("h_ele_matchingObjectPt_matched", "Efficiency vs matching SC E_{T}", nbinpteff, 5., ptmax);
650  h_ele_matchingObjectPt_matched->Sumw2();
651  h_ele_vertexEta = new TH1F("h_ele_vertexEta", "ele momentum eta", nbineta, etamin, etamax);
652  h_ele_vertexEtaVsPhi =
653  new TH2F("h_ele_vertexEtaVsPhi", "ele momentum eta vs phi", nbineta2D, etamin, etamax, nbinphi2D, phimin, phimax);
654  h_ele_matchingObjectAbsEta_matched =
655  new TH1F("h_ele_matchingObjectAbsEta_matched", "Efficiency vs matching SC |#eta|", nbineta / 2, 0., 2.5);
656  h_ele_matchingObjectAbsEta_matched->Sumw2();
657  h_ele_matchingObjectEta_matched =
658  new TH1F("h_ele_matchingObjectEta_matched", "Efficiency vs matching SC #eta", nbineta, etamin, etamax);
659  h_ele_matchingObjectEta_matched->Sumw2();
660  h_ele_matchingObjectPhi_matched =
661  new TH1F("h_ele_matchingObjectPhi_matched", "Efficiency vs matching SC phi", nbinphi, phimin, phimax);
662  h_ele_matchingObjectPhi_matched->Sumw2();
663  h_ele_vertexPhi = new TH1F("h_ele_vertexPhi", "ele momentum #phi", nbinphi, phimin, phimax);
664  h_ele_vertexX = new TH1F("h_ele_vertexX", "ele vertex x", nbinxyz, -0.1, 0.1);
665  h_ele_vertexY = new TH1F("h_ele_vertexY", "ele vertex y", nbinxyz, -0.1, 0.1);
666  h_ele_vertexZ = new TH1F("h_ele_vertexZ", "ele vertex z", nbinxyz, -25, 25);
667  h_ele_matchingObjectZ_matched =
668  new TH1F("h_ele_matchingObjectZ_matched", "Efficiency vs matching SC z", nbinxyz, -25, 25);
669  h_ele_matchingObjectZ_matched->Sumw2();
670  h_ele_vertexTIP = new TH1F("h_ele_vertexTIP", "ele transverse impact parameter (wrt bs)", 90, 0., 0.15);
671  h_ele_vertexTIPVsEta = new TH2F("h_ele_vertexTIPVsEta",
672  "ele transverse impact parameter (wrt bs) vs eta",
673  nbineta2D,
674  etamin,
675  etamax,
676  45,
677  0.,
678  0.15);
679  h_ele_vertexTIPVsPhi = new TH2F("h_ele_vertexTIPVsPhi",
680  "ele transverse impact parameter (wrt bs) vs phi",
681  nbinphi2D,
682  phimin,
683  phimax,
684  45,
685  0.,
686  0.15);
687  h_ele_vertexTIPVsPt = new TH2F("h_ele_vertexTIPVsPt",
688  "ele transverse impact parameter (wrt bs) vs transverse momentum",
689  nbinpt2D,
690  0.,
691  ptmax,
692  45,
693  0.,
694  0.15);
695  h_ele_PoPmatchingObject_matched = new TH1F(
696  "h_ele_PoPmatchingObject_matched", "ele momentum / matching SC energy", nbinpoptrue, poptruemin, poptruemax);
697  h_ele_PtoPtmatchingObject_matched = new TH1F("h_ele_PtoPtmatchingObject_matched",
698  "ele trans momentum / matching SC trans energy",
699  nbinpoptrue,
700  poptruemin,
701  poptruemax);
702  h_ele_PoPmatchingObjectVsEta_matched = new TH2F("h_ele_PoPmatchingObjectVsEta_matched",
703  "ele momentum / matching SC energy vs eta",
704  nbineta2D,
705  etamin,
706  etamax,
707  50,
708  poptruemin,
709  poptruemax);
710  h_ele_PoPmatchingObjectVsPhi_matched = new TH2F("h_ele_PoPmatchingObjectVsPhi_matched",
711  "ele momentum / matching SC energy vs phi",
712  nbinphi2D,
713  phimin,
714  phimax,
715  50,
716  poptruemin,
717  poptruemax);
718  h_ele_PoPmatchingObjectVsPt_matched = new TH2F("h_ele_PoPmatchingObjectVsPt_matched",
719  "ele momentum / matching SC energy vs eta",
720  nbinpt2D,
721  0.,
722  ptmax,
723  50,
724  poptruemin,
725  poptruemax);
726  h_ele_PoPmatchingObject_barrel_matched = new TH1F("h_ele_PoPmatchingObject_barrel_matched",
727  "ele momentum / matching SC energy, barrel",
728  nbinpoptrue,
729  poptruemin,
730  poptruemax);
731  h_ele_PoPmatchingObject_endcaps_matched = new TH1F("h_ele_PoPmatchingObject_endcaps_matched",
732  "ele momentum / matching SC energy, endcaps",
733  nbinpoptrue,
734  poptruemin,
735  poptruemax);
736  h_ele_PtoPtmatchingObject_barrel_matched = new TH1F("h_ele_PtoPmatchingObject_barrel_matched",
737  "ele trans momentum / matching SC trans energy, barrel",
738  nbinpoptrue,
739  poptruemin,
740  poptruemax);
741  h_ele_PtoPtmatchingObject_endcaps_matched = new TH1F("h_ele_PtoPmatchingObject_endcaps_matched",
742  "ele trans momentum / matching SC trans energy, endcaps",
743  nbinpoptrue,
744  poptruemin,
745  poptruemax);
746  h_ele_EtaMnEtamatchingObject_matched = new TH1F(
747  "h_ele_EtaMnEtamatchingObject_matched", "ele momentum eta - matching SC eta", nbindeta, detamin, detamax);
748  h_ele_EtaMnEtamatchingObjectVsEta_matched = new TH2F("h_ele_EtaMnEtamatchingObjectVsEta_matched",
749  "ele momentum eta - matching SC eta vs eta",
750  nbineta2D,
751  etamin,
752  etamax,
753  nbindeta / 2,
754  detamin,
755  detamax);
756  h_ele_EtaMnEtamatchingObjectVsPhi_matched = new TH2F("h_ele_EtaMnEtamatchingObjectVsPhi_matched",
757  "ele momentum eta - matching SC eta vs phi",
758  nbinphi2D,
759  phimin,
760  phimax,
761  nbindeta / 2,
762  detamin,
763  detamax);
764  h_ele_EtaMnEtamatchingObjectVsPt_matched = new TH2F("h_ele_EtaMnEtamatchingObjectVsPt_matched",
765  "ele momentum eta - matching SC eta vs pt",
766  nbinpt,
767  0.,
768  ptmax,
769  nbindeta / 2,
770  detamin,
771  detamax);
772  h_ele_PhiMnPhimatchingObject_matched = new TH1F(
773  "h_ele_PhiMnPhimatchingObject_matched", "ele momentum phi - matching SC phi", nbindphi, dphimin, dphimax);
774  h_ele_PhiMnPhimatchingObject2_matched = new TH1F("h_ele_PhiMnPhimatchingObject2_matched",
775  "ele momentum phi - matching SC phi",
776  nbindphimatch2D,
777  dphimatchmin,
778  dphimatchmax);
779  h_ele_PhiMnPhimatchingObjectVsEta_matched = new TH2F("h_ele_PhiMnPhimatchingObjectVsEta_matched",
780  "ele momentum phi - matching SC phi vs eta",
781  nbineta2D,
782  etamin,
783  etamax,
784  nbindphi / 2,
785  dphimin,
786  dphimax);
787  h_ele_PhiMnPhimatchingObjectVsPhi_matched = new TH2F("h_ele_PhiMnPhimatchingObjectVsPhi_matched",
788  "ele momentum phi - matching SC phi vs phi",
789  nbinphi2D,
790  phimin,
791  phimax,
792  nbindphi / 2,
793  dphimin,
794  dphimax);
795  h_ele_PhiMnPhimatchingObjectVsPt_matched = new TH2F("h_ele_PhiMnPhimatchingObjectVsPt_matched",
796  "ele momentum phi - matching SC phi vs pt",
797  nbinpt2D,
798  0.,
799  ptmax,
800  nbindphi / 2,
801  dphimin,
802  dphimax);
803 
804  // matched electron, superclusters
805  histSclEn_ = new TH1F("h_scl_energy", "ele supercluster energy", nbinp, 0., pmax);
806  histSclEoEmatchingObject_barrel_matched = new TH1F(
807  "h_scl_EoEmatchingObject_barrel_matched", "ele supercluster energy / matching SC energy, barrel", 50, 0.2, 1.2);
808  histSclEoEmatchingObject_endcaps_matched = new TH1F(
809  "h_scl_EoEmatchingObject_endcaps_matched", "ele supercluster energy / matching SC energy, endcaps", 50, 0.2, 1.2);
810  histSclEoEmatchingObject_barrel_new_matched = new TH1F("h_scl_EoEmatchingObject_barrel_new_matched",
811  "ele supercluster energy / matching SC energy, barrel",
812  nbinpoptrue,
813  poptruemin,
814  poptruemax);
815  histSclEoEmatchingObject_endcaps_new_matched = new TH1F("h_scl_EoEmatchingObject_endcaps_new_matched",
816  "ele supercluster energy / matching SC energy, endcaps",
817  nbinpoptrue,
818  poptruemin,
819  poptruemax);
820  histSclEt_ = new TH1F("h_scl_et", "ele supercluster transverse energy", nbinpt, 0., ptmax);
821  histSclEtVsEta_ = new TH2F(
822  "h_scl_etVsEta", "ele supercluster transverse energy vs eta", nbineta2D, etamin, etamax, nbinpt, 0., ptmax);
823  histSclEtVsPhi_ = new TH2F(
824  "h_scl_etVsPhi", "ele supercluster transverse energy vs phi", nbinphi2D, phimin, phimax, nbinpt, 0., ptmax);
825  histSclEtaVsPhi_ =
826  new TH2F("h_scl_etaVsPhi", "ele supercluster eta vs phi", nbinphi2D, phimin, phimax, nbineta2D, etamin, etamax);
827  histSclEta_ = new TH1F("h_scl_eta", "ele supercluster eta", nbineta, etamin, etamax);
828  histSclPhi_ = new TH1F("h_scl_phi", "ele supercluster phi", nbinphi, phimin, phimax);
829 
830  histSclSigEtaEta_ = new TH1F("h_scl_sigetaeta", "ele supercluster sigma eta eta", 100, 0., 0.05);
831  histSclSigIEtaIEta_barrel_ =
832  new TH1F("h_scl_sigietaieta_barrel", "ele supercluster sigma ieta ieta, barrel", 100, 0., 0.05);
833  histSclSigIEtaIEta_endcaps_ =
834  new TH1F("h_scl_sigietaieta_endcaps", "ele supercluster sigma ieta ieta, endcaps", 100, 0., 0.05);
835  histSclE1x5_ = new TH1F("h_scl_E1x5", "ele supercluster energy in 1x5", nbinp, 0., pmax);
836  histSclE1x5_barrel_ = new TH1F("h_scl_E1x_barrel5", "ele supercluster energy in 1x5 barrel", nbinp, 0., pmax);
837  histSclE1x5_endcaps_ = new TH1F("h_scl_E1x5_endcaps", "ele supercluster energy in 1x5 endcaps", nbinp, 0., pmax);
838  histSclE2x5max_ = new TH1F("h_scl_E2x5max", "ele supercluster energy in 2x5 max", nbinp, 0., pmax);
839  histSclE2x5max_barrel_ =
840  new TH1F("h_scl_E2x5max_barrel", "ele supercluster energy in 2x5 max barrel", nbinp, 0., pmax);
841  histSclE2x5max_endcaps_ =
842  new TH1F("h_scl_E2x5max_endcaps", "ele supercluster energy in 2x5 max endcaps", nbinp, 0., pmax);
843  histSclE5x5_ = new TH1F("h_scl_E5x5", "ele supercluster energy in 5x5", nbinp, 0., pmax);
844  histSclE5x5_barrel_ = new TH1F("h_scl_E5x5_barrel", "ele supercluster energy in 5x5 barrel", nbinp, 0., pmax);
845  histSclE5x5_endcaps_ = new TH1F("h_scl_E5x5_endcaps", "ele supercluster energy in 5x5 endcaps", nbinp, 0., pmax);
846 
847  // matched electron, gsf tracks
848  h_ele_ambiguousTracks = new TH1F("h_ele_ambiguousTracks", "ele # ambiguous tracks", 5, 0., 5.);
849  h_ele_ambiguousTracksVsEta =
850  new TH2F("h_ele_ambiguousTracksVsEta", "ele # ambiguous tracks vs eta", nbineta2D, etamin, etamax, 5, 0., 5.);
851  h_ele_ambiguousTracksVsPhi =
852  new TH2F("h_ele_ambiguousTracksVsPhi", "ele # ambiguous tracks vs phi", nbinphi2D, phimin, phimax, 5, 0., 5.);
853  h_ele_ambiguousTracksVsPt =
854  new TH2F("h_ele_ambiguousTracksVsPt", "ele # ambiguous tracks vs pt", nbinpt2D, 0., ptmax, 5, 0., 5.);
855  h_ele_foundHits = new TH1F("h_ele_foundHits", "ele track # found hits", nbinfhits, 0., fhitsmax);
856  h_ele_foundHitsVsEta = new TH2F(
857  "h_ele_foundHitsVsEta", "ele track # found hits vs eta", nbineta2D, etamin, etamax, nbinfhits, 0., fhitsmax);
858  h_ele_foundHitsVsPhi = new TH2F(
859  "h_ele_foundHitsVsPhi", "ele track # found hits vs phi", nbinphi2D, phimin, phimax, nbinfhits, 0., fhitsmax);
860  h_ele_foundHitsVsPt =
861  new TH2F("h_ele_foundHitsVsPt", "ele track # found hits vs pt", nbinpt2D, 0., ptmax, nbinfhits, 0., fhitsmax);
862  h_ele_lostHits = new TH1F("h_ele_lostHits", "ele track # lost hits", 5, 0., 5.);
863  h_ele_lostHitsVsEta = new TH2F(
864  "h_ele_lostHitsVsEta", "ele track # lost hits vs eta", nbineta2D, etamin, etamax, nbinlhits, 0., lhitsmax);
865  h_ele_lostHitsVsPhi = new TH2F(
866  "h_ele_lostHitsVsPhi", "ele track # lost hits vs eta", nbinphi2D, phimin, phimax, nbinlhits, 0., lhitsmax);
867  h_ele_lostHitsVsPt =
868  new TH2F("h_ele_lostHitsVsPt", "ele track # lost hits vs eta", nbinpt2D, 0., ptmax, nbinlhits, 0., lhitsmax);
869  h_ele_chi2 = new TH1F("h_ele_chi2", "ele track #chi^{2}", 100, 0., 15.);
870  h_ele_chi2VsEta = new TH2F("h_ele_chi2VsEta", "ele track #chi^{2} vs eta", nbineta2D, etamin, etamax, 50, 0., 15.);
871  h_ele_chi2VsPhi = new TH2F("h_ele_chi2VsPhi", "ele track #chi^{2} vs phi", nbinphi2D, phimin, phimax, 50, 0., 15.);
872  h_ele_chi2VsPt = new TH2F("h_ele_chi2VsPt", "ele track #chi^{2} vs pt", nbinpt2D, 0., ptmax, 50, 0., 15.);
873  h_ele_PinMnPout = new TH1F("h_ele_PinMnPout", "ele track inner p - outer p, mean", nbinp, 0., 200.);
874  h_ele_PinMnPout_mode = new TH1F("h_ele_PinMnPout_mode", "ele track inner p - outer p, mode", nbinp, 0., 100.);
875  h_ele_PinMnPoutVsEta_mode = new TH2F("h_ele_PinMnPoutVsEta_mode",
876  "ele track inner p - outer p vs eta, mode",
877  nbineta2D,
878  etamin,
879  etamax,
880  nbinp2D,
881  0.,
882  100.);
883  h_ele_PinMnPoutVsPhi_mode = new TH2F("h_ele_PinMnPoutVsPhi_mode",
884  "ele track inner p - outer p vs phi, mode",
885  nbinphi2D,
886  phimin,
887  phimax,
888  nbinp2D,
889  0.,
890  100.);
891  h_ele_PinMnPoutVsPt_mode = new TH2F(
892  "h_ele_PinMnPoutVsPt_mode", "ele track inner p - outer p vs pt, mode", nbinpt2D, 0., ptmax, nbinp2D, 0., 100.);
893  h_ele_PinMnPoutVsE_mode = new TH2F(
894  "h_ele_PinMnPoutVsE_mode", "ele track inner p - outer p vs E, mode", nbinp2D, 0., 200., nbinp2D, 0., 100.);
895  h_ele_PinMnPoutVsChi2_mode = new TH2F(
896  "h_ele_PinMnPoutVsChi2_mode", "ele track inner p - outer p vs track chi2, mode", 50, 0., 20., nbinp2D, 0., 100.);
897  h_ele_outerP = new TH1F("h_ele_outerP", "ele track outer p, mean", nbinp, 0., pmax);
898  h_ele_outerP_mode = new TH1F("h_ele_outerP_mode", "ele track outer p, mode", nbinp, 0., pmax);
899  h_ele_outerPVsEta_mode =
900  new TH2F("h_ele_outerPVsEta_mode", "ele track outer p vs eta mode", nbineta2D, etamin, etamax, 50, 0., pmax);
901  h_ele_outerPt = new TH1F("h_ele_outerPt", "ele track outer p_{T}, mean", nbinpt, 0., ptmax);
902  h_ele_outerPt_mode = new TH1F("h_ele_outerPt_mode", "ele track outer p_{T}, mode", nbinpt, 0., ptmax);
903  h_ele_outerPtVsEta_mode = new TH2F(
904  "h_ele_outerPtVsEta_mode", "ele track outer p_{T} vs eta, mode", nbineta2D, etamin, etamax, nbinpt2D, 0., ptmax);
905  h_ele_outerPtVsPhi_mode = new TH2F(
906  "h_ele_outerPtVsPhi_mode", "ele track outer p_{T} vs phi, mode", nbinphi2D, phimin, phimax, nbinpt2D, 0., ptmax);
907  h_ele_outerPtVsPt_mode =
908  new TH2F("h_ele_outerPtVsPt_mode", "ele track outer p_{T} vs pt, mode", nbinpt2D, 0., 100., nbinpt2D, 0., ptmax);
909 
910  // matched electrons, matching
911  h_ele_EoP = new TH1F("h_ele_EoP", "ele E/P_{vertex}", nbineop, 0., eopmax);
912  h_ele_EoPVsEta =
913  new TH2F("h_ele_EoPVsEta", "ele E/P_{vertex} vs eta", nbineta2D, etamin, etamax, nbineop2D, 0., eopmaxsht);
914  h_ele_EoPVsPhi =
915  new TH2F("h_ele_EoPVsPhi", "ele E/P_{vertex} vs phi", nbinphi2D, phimin, phimax, nbineop2D, 0., eopmaxsht);
916  h_ele_EoPVsE = new TH2F("h_ele_EoPVsE", "ele E/P_{vertex} vs E", 50, 0., pmax, 50, 0., 5.);
917  h_ele_EseedOP = new TH1F("h_ele_EseedOP", "ele E_{seed}/P_{vertex}", nbineop, 0., eopmax);
918  h_ele_EseedOPVsEta = new TH2F(
919  "h_ele_EseedOPVsEta", "ele E_{seed}/P_{vertex} vs eta", nbineta2D, etamin, etamax, nbineop2D, 0., eopmaxsht);
920  h_ele_EseedOPVsPhi = new TH2F(
921  "h_ele_EseedOPVsPhi", "ele E_{seed}/P_{vertex} vs phi", nbinphi2D, phimin, phimax, nbineop2D, 0., eopmaxsht);
922  h_ele_EseedOPVsE = new TH2F("h_ele_EseedOPVsE", "ele E_{seed}/P_{vertex} vs E", 50, 0., pmax, 50, 0., 5.);
923  h_ele_EoPout = new TH1F("h_ele_EoPout", "ele E/P_{out}", nbineop, 0., eopmax);
924  h_ele_EoPoutVsEta =
925  new TH2F("h_ele_EoPoutVsEta", "ele E/P_{out} vs eta", nbineta2D, etamin, etamax, nbineop2D, 0., eopmaxsht);
926  h_ele_EoPoutVsPhi =
927  new TH2F("h_ele_EoPoutVsPhi", "ele E/P_{out} vs phi", nbinphi2D, phimin, phimax, nbineop2D, 0., eopmaxsht);
928  h_ele_EoPoutVsE = new TH2F("h_ele_EoPoutVsE", "ele E/P_{out} vs E", nbinp2D, 0., pmax, nbineop2D, 0., eopmaxsht);
929  h_ele_EeleOPout = new TH1F("h_ele_EeleOPout", "ele E_{ele}/P_{out}", nbineop, 0., eopmax);
930  h_ele_EeleOPoutVsEta = new TH2F(
931  "h_ele_EeleOPoutVsEta", "ele E_{ele}/P_{out} vs eta", nbineta2D, etamin, etamax, nbineop2D, 0., eopmaxsht);
932  h_ele_EeleOPoutVsPhi = new TH2F(
933  "h_ele_EeleOPoutVsPhi", "ele E_{ele}/P_{out} vs phi", nbinphi2D, phimin, phimax, nbineop2D, 0., eopmaxsht);
934  h_ele_EeleOPoutVsE =
935  new TH2F("h_ele_EeleOPoutVsE", "ele E_{ele}/P_{out} vs E", nbinp2D, 0., pmax, nbineop2D, 0., eopmaxsht);
936  h_ele_dEtaSc_propVtx = new TH1F(
937  "h_ele_dEtaSc_propVtx", "ele #eta_{sc} - #eta_{tr}, prop from vertex", nbindetamatch, detamatchmin, detamatchmax);
938  h_ele_dEtaScVsEta_propVtx = new TH2F("h_ele_dEtaScVsEta_propVtx",
939  "ele #eta_{sc} - #eta_{tr} vs eta, prop from vertex",
940  nbineta2D,
941  etamin,
942  etamax,
943  nbindetamatch2D,
944  detamatchmin,
945  detamatchmax);
946  h_ele_dEtaScVsPhi_propVtx = new TH2F("h_ele_dEtaScVsPhi_propVtx",
947  "ele #eta_{sc} - #eta_{tr} vs phi, prop from vertex",
948  nbinphi2D,
949  phimin,
950  phimax,
951  nbindetamatch2D,
952  detamatchmin,
953  detamatchmax);
954  h_ele_dEtaScVsPt_propVtx = new TH2F("h_ele_dEtaScVsPt_propVtx",
955  "ele #eta_{sc} - #eta_{tr} vs pt, prop from vertex",
956  nbinpt2D,
957  0.,
958  ptmax,
959  nbindetamatch2D,
960  detamatchmin,
961  detamatchmax);
962  h_ele_dPhiSc_propVtx = new TH1F(
963  "h_ele_dPhiSc_propVtx", "ele #phi_{sc} - #phi_{tr}, prop from vertex", nbindphimatch, dphimatchmin, dphimatchmax);
964  h_ele_dPhiScVsEta_propVtx = new TH2F("h_ele_dPhiScVsEta_propVtx",
965  "ele #phi_{sc} - #phi_{tr} vs eta, prop from vertex",
966  nbineta2D,
967  etamin,
968  etamax,
969  nbindphimatch2D,
970  dphimatchmin,
971  dphimatchmax);
972  h_ele_dPhiScVsPhi_propVtx = new TH2F("h_ele_dPhiScVsPhi_propVtx",
973  "ele #phi_{sc} - #phi_{tr} vs phi, prop from vertex",
974  nbinphi2D,
975  phimin,
976  phimax,
977  nbindphimatch2D,
978  dphimatchmin,
979  dphimatchmax);
980  h_ele_dPhiScVsPt_propVtx = new TH2F("h_ele_dPhiScVsPt_propVtx",
981  "ele #phi_{sc} - #phi_{tr} vs pt, prop from vertex",
982  nbinpt2D,
983  0.,
984  ptmax,
985  nbindphimatch2D,
986  dphimatchmin,
987  dphimatchmax);
988  h_ele_dEtaCl_propOut = new TH1F("h_ele_dEtaCl_propOut",
989  "ele #eta_{cl} - #eta_{tr}, prop from outermost",
990  nbindetamatch,
991  detamatchmin,
992  detamatchmax);
993  h_ele_dEtaClVsEta_propOut = new TH2F("h_ele_dEtaClVsEta_propOut",
994  "ele #eta_{cl} - #eta_{tr} vs eta, prop from out",
995  nbineta2D,
996  etamin,
997  etamax,
998  nbindetamatch2D,
999  detamatchmin,
1000  detamatchmax);
1001  h_ele_dEtaClVsPhi_propOut = new TH2F("h_ele_dEtaClVsPhi_propOut",
1002  "ele #eta_{cl} - #eta_{tr} vs phi, prop from out",
1003  nbinphi2D,
1004  phimin,
1005  phimax,
1006  nbindetamatch2D,
1007  detamatchmin,
1008  detamatchmax);
1009  h_ele_dEtaClVsPt_propOut = new TH2F("h_ele_dEtaScVsPt_propOut",
1010  "ele #eta_{cl} - #eta_{tr} vs pt, prop from out",
1011  nbinpt2D,
1012  0.,
1013  ptmax,
1014  nbindetamatch2D,
1015  detamatchmin,
1016  detamatchmax);
1017  h_ele_dPhiCl_propOut = new TH1F("h_ele_dPhiCl_propOut",
1018  "ele #phi_{cl} - #phi_{tr}, prop from outermost",
1019  nbindphimatch,
1020  dphimatchmin,
1021  dphimatchmax);
1022  h_ele_dPhiClVsEta_propOut = new TH2F("h_ele_dPhiClVsEta_propOut",
1023  "ele #phi_{cl} - #phi_{tr} vs eta, prop from out",
1024  nbineta2D,
1025  etamin,
1026  etamax,
1027  nbindphimatch2D,
1028  dphimatchmin,
1029  dphimatchmax);
1030  h_ele_dPhiClVsPhi_propOut = new TH2F("h_ele_dPhiClVsPhi_propOut",
1031  "ele #phi_{cl} - #phi_{tr} vs phi, prop from out",
1032  nbinphi2D,
1033  phimin,
1034  phimax,
1035  nbindphimatch2D,
1036  dphimatchmin,
1037  dphimatchmax);
1038  h_ele_dPhiClVsPt_propOut = new TH2F("h_ele_dPhiSClsPt_propOut",
1039  "ele #phi_{cl} - #phi_{tr} vs pt, prop from out",
1040  nbinpt2D,
1041  0.,
1042  ptmax,
1043  nbindphimatch2D,
1044  dphimatchmin,
1045  dphimatchmax);
1046  h_ele_dEtaEleCl_propOut = new TH1F("h_ele_dEtaEleCl_propOut",
1047  "ele #eta_{EleCl} - #eta_{tr}, prop from outermost",
1048  nbindetamatch,
1049  detamatchmin,
1050  detamatchmax);
1051  h_ele_dEtaEleClVsEta_propOut = new TH2F("h_ele_dEtaEleClVsEta_propOut",
1052  "ele #eta_{EleCl} - #eta_{tr} vs eta, prop from out",
1053  nbineta2D,
1054  etamin,
1055  etamax,
1056  nbindetamatch2D,
1057  detamatchmin,
1058  detamatchmax);
1059  h_ele_dEtaEleClVsPhi_propOut = new TH2F("h_ele_dEtaEleClVsPhi_propOut",
1060  "ele #eta_{EleCl} - #eta_{tr} vs phi, prop from out",
1061  nbinphi2D,
1062  phimin,
1063  phimax,
1064  nbindetamatch2D,
1065  detamatchmin,
1066  detamatchmax);
1067  h_ele_dEtaEleClVsPt_propOut = new TH2F("h_ele_dEtaScVsPt_propOut",
1068  "ele #eta_{EleCl} - #eta_{tr} vs pt, prop from out",
1069  nbinpt2D,
1070  0.,
1071  ptmax,
1072  nbindetamatch2D,
1073  detamatchmin,
1074  detamatchmax);
1075  h_ele_dPhiEleCl_propOut = new TH1F("h_ele_dPhiEleCl_propOut",
1076  "ele #phi_{EleCl} - #phi_{tr}, prop from outermost",
1077  nbindphimatch,
1078  dphimatchmin,
1079  dphimatchmax);
1080  h_ele_dPhiEleClVsEta_propOut = new TH2F("h_ele_dPhiEleClVsEta_propOut",
1081  "ele #phi_{EleCl} - #phi_{tr} vs eta, prop from out",
1082  nbineta2D,
1083  etamin,
1084  etamax,
1085  nbindphimatch2D,
1086  dphimatchmin,
1087  dphimatchmax);
1088  h_ele_dPhiEleClVsPhi_propOut = new TH2F("h_ele_dPhiEleClVsPhi_propOut",
1089  "ele #phi_{EleCl} - #phi_{tr} vs phi, prop from out",
1090  nbinphi2D,
1091  phimin,
1092  phimax,
1093  nbindphimatch2D,
1094  dphimatchmin,
1095  dphimatchmax);
1096  h_ele_dPhiEleClVsPt_propOut = new TH2F("h_ele_dPhiSEleClsPt_propOut",
1097  "ele #phi_{EleCl} - #phi_{tr} vs pt, prop from out",
1098  nbinpt2D,
1099  0.,
1100  ptmax,
1101  nbindphimatch2D,
1102  dphimatchmin,
1103  dphimatchmax);
1104 
1105  h_ele_HoE = new TH1F("h_ele_HoE", "ele hadronic energy / em energy", nbinhoe, hoemin, hoemax);
1106  h_ele_HoE_fiducial =
1107  new TH1F("h_ele_HoE_fiducial", "ele hadronic energy / em energy, fiducial region", nbinhoe, hoemin, hoemax);
1108  h_ele_HoEVsEta = new TH2F(
1109  "h_ele_HoEVsEta", "ele hadronic energy / em energy vs eta", nbineta, etamin, etamax, nbinhoe, hoemin, hoemax);
1110  h_ele_HoEVsPhi = new TH2F(
1111  "h_ele_HoEVsPhi", "ele hadronic energy / em energy vs phi", nbinphi2D, phimin, phimax, nbinhoe, hoemin, hoemax);
1112  h_ele_HoEVsE =
1113  new TH2F("h_ele_HoEVsE", "ele hadronic energy / em energy vs E", nbinp, 0., 300., nbinhoe, hoemin, hoemax);
1114 
1115  h_ele_seed_dphi2_ = new TH1F("h_ele_seedDphi2", "ele seed dphi 2nd layer", 50, -0.003, +0.003);
1116  h_ele_seed_dphi2VsEta_ =
1117  new TH2F("h_ele_seedDphi2VsEta", "ele seed dphi 2nd layer vs eta", nbineta2D, etamin, etamax, 50, -0.003, +0.003);
1118  h_ele_seed_dphi2VsPt_ =
1119  new TH2F("h_ele_seedDphi2VsPt", "ele seed dphi 2nd layer vs pt", nbinpt2D, 0., ptmax, 50, -0.003, +0.003);
1120  h_ele_seed_drz2_ = new TH1F("h_ele_seedDrz2", "ele seed dr/dz 2nd layer", 50, -0.03, +0.03);
1121  h_ele_seed_drz2VsEta_ =
1122  new TH2F("h_ele_seedDrz2VsEta", "ele seed dr/dz 2nd layer vs eta", nbineta2D, etamin, etamax, 50, -0.03, +0.03);
1123  h_ele_seed_drz2VsPt_ =
1124  new TH2F("h_ele_seedDrz2VsPt", "ele seed dr/dz 2nd layer vs pt", nbinpt2D, 0., ptmax, 50, -0.03, +0.03);
1125  h_ele_seed_subdet2_ = new TH1F("h_ele_seedSubdet2", "ele seed subdet 2nd layer", 10, 0., 10.);
1126 
1127  // classes
1128  h_ele_classes = new TH1F("h_ele_classes", "electron classes", 20, 0.0, 20.);
1129  h_ele_eta = new TH1F("h_ele_eta", "ele electron eta", nbineta / 2, 0.0, etamax);
1130  h_ele_eta_golden = new TH1F("h_ele_eta_golden", "ele electron eta golden", nbineta / 2, 0.0, etamax);
1131  h_ele_eta_bbrem = new TH1F("h_ele_eta_bbrem", "ele electron eta bbrem", nbineta / 2, 0.0, etamax);
1132  h_ele_eta_narrow = new TH1F("h_ele_eta_narrow", "ele electron eta narrow", nbineta / 2, 0.0, etamax);
1133  h_ele_eta_shower = new TH1F("h_ele_eta_show", "ele electron eta showering", nbineta / 2, 0.0, etamax);
1134  h_ele_PinVsPoutGolden_mode = new TH2F("h_ele_PinVsPoutGolden_mode",
1135  "ele track inner p vs outer p vs eta, golden, mode",
1136  nbinp2D,
1137  0.,
1138  pmax,
1139  50,
1140  0.,
1141  pmax);
1142  h_ele_PinVsPoutShowering_mode = new TH2F("h_ele_PinVsPoutShowering_mode",
1143  "ele track inner p vs outer p vs eta, Showering, mode",
1144  nbinp2D,
1145  0.,
1146  pmax,
1147  50,
1148  0.,
1149  pmax);
1150  h_ele_PinVsPoutGolden_mean = new TH2F("h_ele_PinVsPoutGolden_mean",
1151  "ele track inner p vs outer p vs eta, golden, mean",
1152  nbinp2D,
1153  0.,
1154  pmax,
1155  50,
1156  0.,
1157  pmax);
1158  h_ele_PinVsPoutShowering_mean = new TH2F("h_ele_PinVsPoutShowering_mean",
1159  "ele track inner p vs outer p vs eta, Showering, mean",
1160  nbinp2D,
1161  0.,
1162  pmax,
1163  50,
1164  0.,
1165  pmax);
1166  h_ele_PtinVsPtoutGolden_mode = new TH2F("h_ele_PtinVsPtoutGolden_mode",
1167  "ele track inner pt vs outer pt vs eta, golden, mode",
1168  nbinpt2D,
1169  0.,
1170  ptmax,
1171  50,
1172  0.,
1173  ptmax);
1174  h_ele_PtinVsPtoutShowering_mode = new TH2F("h_ele_PtinVsPtoutShowering_mode",
1175  "ele track inner pt vs outer pt vs eta, showering, mode",
1176  nbinpt2D,
1177  0.,
1178  ptmax,
1179  50,
1180  0.,
1181  ptmax);
1182  h_ele_PtinVsPtoutGolden_mean = new TH2F("h_ele_PtinVsPtoutGolden_mean",
1183  "ele track inner pt vs outer pt vs eta, golden, mean",
1184  nbinpt2D,
1185  0.,
1186  ptmax,
1187  50,
1188  0.,
1189  ptmax);
1190  h_ele_PtinVsPtoutShowering_mean = new TH2F("h_ele_PtinVsPtoutShowering_mean",
1191  "ele track inner pt vs outer pt vs eta, showering, mean",
1192  nbinpt2D,
1193  0.,
1194  ptmax,
1195  50,
1196  0.,
1197  ptmax);
1198  histSclEoEmatchingObjectGolden_barrel = new TH1F("h_scl_EoEmatchingObject golden, barrel",
1199  "ele supercluster energy over matchingObject energy, golden, barrel",
1200  100,
1201  0.2,
1202  1.2);
1203  histSclEoEmatchingObjectGolden_endcaps =
1204  new TH1F("h_scl_EoEmatchingObject golden, endcaps",
1205  "ele supercluster energy over matchingObject energy, golden, endcaps",
1206  100,
1207  0.2,
1208  1.2);
1209  histSclEoEmatchingObjectShowering_barrel =
1210  new TH1F("h_scl_EoEmatchingObject Showering, barrel",
1211  "ele supercluster energy over matchingObject energy, showering, barrel",
1212  100,
1213  0.2,
1214  1.2);
1215  histSclEoEmatchingObjectShowering_endcaps =
1216  new TH1F("h_scl_EoEmatchingObject Showering, endcaps",
1217  "ele supercluster energy over matchingObject energy, showering, endcaps",
1218  100,
1219  0.2,
1220  1.2);
1221 
1222  // isolation
1223  h_ele_tkSumPt_dr03 = new TH1F("h_ele_tkSumPt_dr03", "tk isolation sum, dR=0.3", 100, 0.0, 20.);
1224  h_ele_ecalRecHitSumEt_dr03 = new TH1F("h_ele_ecalRecHitSumEt_dr03", "ecal isolation sum, dR=0.3", 100, 0.0, 20.);
1225  h_ele_hcalDepth1TowerSumEt_dr03 =
1226  new TH1F("h_ele_hcalDepth1TowerSumEt_dr03", "hcal depth1 isolation sum, dR=0.3", 100, 0.0, 20.);
1227  h_ele_hcalDepth2TowerSumEt_dr03 =
1228  new TH1F("h_ele_hcalDepth2TowerSumEt_dr03", "hcal depth2 isolation sum, dR=0.3", 100, 0.0, 20.);
1229  h_ele_tkSumPt_dr04 = new TH1F("h_ele_tkSumPt_dr04", "hcal isolation sum", 100, 0.0, 20.);
1230  h_ele_ecalRecHitSumEt_dr04 = new TH1F("h_ele_ecalRecHitSumEt_dr04", "ecal isolation sum, dR=0.4", 100, 0.0, 20.);
1231  h_ele_hcalDepth1TowerSumEt_dr04 =
1232  new TH1F("h_ele_hcalDepth1TowerSumEt_dr04", "hcal depth1 isolation sum, dR=0.4", 100, 0.0, 20.);
1233  h_ele_hcalDepth2TowerSumEt_dr04 =
1234  new TH1F("h_ele_hcalDepth2TowerSumEt_dr04", "hcal depth2 isolation sum, dR=0.4", 100, 0.0, 20.);
1235 
1236  // fbrem
1237  h_ele_fbrem = new TH1F("h_ele_fbrem", "ele brem fraction, mode", 100, 0., 1.);
1238  h_ele_fbremVsEta_mode =
1239  new TProfile("h_ele_fbremvsEtamode", "mean ele brem fraction vs eta, mode", nbineta2D, etamin, etamax, 0., 1.);
1240  h_ele_fbremVsEta_mean =
1241  new TProfile("h_ele_fbremvsEtamean", "mean ele brem fraction vs eta, mean", nbineta2D, etamin, etamax, 0., 1.);
1242 
1243  // e/g et pflow electrons
1244  h_ele_mva = new TH1F("h_ele_mva", "ele identification mva", 100, -1., 1.);
1245  h_ele_provenance = new TH1F("h_ele_provenance", "ele provenance", 5, -2., 3.);
1246 
1247  // histos titles
1248  h_matchingObjectNum->GetXaxis()->SetTitle("N_{SC}");
1249  h_matchingObjectNum->GetYaxis()->SetTitle("Events");
1250  h_matchingObjectEta->GetXaxis()->SetTitle("#eta_{SC}");
1251  h_matchingObjectEta->GetYaxis()->SetTitle("Events");
1252  h_matchingObjectP->GetXaxis()->SetTitle("E_{SC} (GeV)");
1253  h_matchingObjectP->GetYaxis()->SetTitle("Events");
1254  h_ele_foundHits->GetXaxis()->SetTitle("N_{hits}");
1255  h_ele_foundHits->GetYaxis()->SetTitle("Events");
1256  h_ele_ambiguousTracks->GetXaxis()->SetTitle("N_{ambiguous tracks}");
1257  h_ele_ambiguousTracks->GetYaxis()->SetTitle("Events");
1258  h_ele_lostHits->GetXaxis()->SetTitle("N_{lost hits}");
1259  h_ele_lostHits->GetYaxis()->SetTitle("Events");
1260  h_ele_chi2->GetXaxis()->SetTitle("#Chi^{2}");
1261  h_ele_chi2->GetYaxis()->SetTitle("Events");
1262  h_ele_charge->GetXaxis()->SetTitle("charge");
1263  h_ele_charge->GetYaxis()->SetTitle("Events");
1264  h_ele_vertexP->GetXaxis()->SetTitle("p_{vertex} (GeV/c)");
1265  h_ele_vertexP->GetYaxis()->SetTitle("Events");
1266  h_ele_vertexPt->GetXaxis()->SetTitle("p_{T vertex} (GeV/c)");
1267  h_ele_vertexPt->GetYaxis()->SetTitle("Events");
1268  h_ele_Et->GetXaxis()->SetTitle("E_{T} (GeV)");
1269  h_ele_Et->GetYaxis()->SetTitle("Events");
1270  h_ele_vertexEta->GetXaxis()->SetTitle("#eta");
1271  h_ele_vertexEta->GetYaxis()->SetTitle("Events");
1272  h_ele_vertexPhi->GetXaxis()->SetTitle("#phi (rad)");
1273  h_ele_vertexPhi->GetYaxis()->SetTitle("Events");
1274  h_ele_PoPmatchingObject_matched->GetXaxis()->SetTitle("P/E_{SC}");
1275  h_ele_PoPmatchingObject_matched->GetYaxis()->SetTitle("Events");
1276  h_ele_PoPmatchingObject_barrel_matched->GetXaxis()->SetTitle("P/E_{SC}");
1277  h_ele_PoPmatchingObject_barrel_matched->GetYaxis()->SetTitle("Events");
1278  h_ele_PoPmatchingObject_endcaps_matched->GetXaxis()->SetTitle("P/E_{SC}");
1279  h_ele_PoPmatchingObject_endcaps_matched->GetYaxis()->SetTitle("Events");
1280  h_ele_PtoPtmatchingObject_matched->GetXaxis()->SetTitle("P_{T}/E_{T}^{SC}");
1281  h_ele_PtoPtmatchingObject_matched->GetYaxis()->SetTitle("Events");
1282  h_ele_PtoPtmatchingObject_barrel_matched->GetXaxis()->SetTitle("P_{T}/E_{T}^{SC}");
1283  h_ele_PtoPtmatchingObject_barrel_matched->GetYaxis()->SetTitle("Events");
1284  h_ele_PtoPtmatchingObject_endcaps_matched->GetXaxis()->SetTitle("P_{T}/E_{T}^{SC}");
1285  h_ele_PtoPtmatchingObject_endcaps_matched->GetYaxis()->SetTitle("Events");
1286  histSclSigEtaEta_->GetXaxis()->SetTitle("#sigma_{#eta #eta}");
1287  histSclSigEtaEta_->GetYaxis()->SetTitle("Events");
1288  histSclSigIEtaIEta_barrel_->GetXaxis()->SetTitle("#sigma_{i#eta i#eta}");
1289  histSclSigIEtaIEta_barrel_->GetYaxis()->SetTitle("Events");
1290  histSclSigIEtaIEta_endcaps_->GetXaxis()->SetTitle("#sigma_{i#eta i#eta}");
1291  histSclSigIEtaIEta_endcaps_->GetYaxis()->SetTitle("Events");
1292  histSclE1x5_->GetXaxis()->SetTitle("E1x5 (GeV)");
1293  histSclE1x5_->GetYaxis()->SetTitle("Events");
1294  histSclE1x5_barrel_->GetXaxis()->SetTitle("E1x5 (GeV)");
1295  histSclE1x5_barrel_->GetYaxis()->SetTitle("Events");
1296  histSclE1x5_endcaps_->GetXaxis()->SetTitle("E1x5 (GeV)");
1297  histSclE1x5_endcaps_->GetYaxis()->SetTitle("Events");
1298  histSclE2x5max_->GetXaxis()->SetTitle("E2x5 (GeV)");
1299  histSclE2x5max_->GetYaxis()->SetTitle("Events");
1300  histSclE2x5max_barrel_->GetXaxis()->SetTitle("E2x5 (GeV)");
1301  histSclE2x5max_barrel_->GetYaxis()->SetTitle("Events");
1302  histSclE2x5max_endcaps_->GetXaxis()->SetTitle("E2x5 (GeV)");
1303  histSclE2x5max_endcaps_->GetYaxis()->SetTitle("Events");
1304  histSclE5x5_->GetXaxis()->SetTitle("E5x5 (GeV)");
1305  histSclE5x5_->GetYaxis()->SetTitle("Events");
1306  histSclE5x5_barrel_->GetXaxis()->SetTitle("E5x5 (GeV)");
1307  histSclE5x5_barrel_->GetYaxis()->SetTitle("Events");
1308  histSclE5x5_endcaps_->GetXaxis()->SetTitle("E5x5 (GeV)");
1309  histSclE5x5_endcaps_->GetYaxis()->SetTitle("Events");
1310  h_ele_EtaMnEtamatchingObject_matched->GetXaxis()->SetTitle("#eta_{rec} - #eta_{SC}");
1311  h_ele_EtaMnEtamatchingObject_matched->GetYaxis()->SetTitle("Events");
1312  h_ele_PhiMnPhimatchingObject_matched->GetXaxis()->SetTitle("#phi_{rec} - #phi_{SC} (rad)");
1313  h_ele_PhiMnPhimatchingObject_matched->GetYaxis()->SetTitle("Events");
1314  h_ele_PinMnPout->GetXaxis()->SetTitle("P_{vertex} - P_{out} (GeV/c)");
1315  h_ele_PinMnPout->GetYaxis()->SetTitle("Events");
1316  h_ele_PinMnPout_mode->GetXaxis()->SetTitle("P_{vertex} - P_{out}, mode (GeV/c)");
1317  h_ele_PinMnPout_mode->GetYaxis()->SetTitle("Events");
1318  h_ele_outerP->GetXaxis()->SetTitle("P_{out} (GeV/c)");
1319  h_ele_outerP->GetYaxis()->SetTitle("Events");
1320  h_ele_outerP_mode->GetXaxis()->SetTitle("P_{out} (GeV/c)");
1321  h_ele_outerP_mode->GetYaxis()->SetTitle("Events");
1322  h_ele_outerPt->GetXaxis()->SetTitle("P_{T out} (GeV/c)");
1323  h_ele_outerPt->GetYaxis()->SetTitle("Events");
1324  h_ele_outerPt_mode->GetXaxis()->SetTitle("P_{T out} (GeV/c)");
1325  h_ele_outerPt_mode->GetYaxis()->SetTitle("Events");
1326  h_ele_EoP->GetXaxis()->SetTitle("E/P_{vertex}");
1327  h_ele_EoP->GetYaxis()->SetTitle("Events");
1328  h_ele_EseedOP->GetXaxis()->SetTitle("E_{seed}/P_{vertex}");
1329  h_ele_EseedOP->GetYaxis()->SetTitle("Events");
1330  h_ele_EoPout->GetXaxis()->SetTitle("E_{seed}/P_{out}");
1331  h_ele_EoPout->GetYaxis()->SetTitle("Events");
1332  h_ele_EeleOPout->GetXaxis()->SetTitle("E_{ele}/P_{out}");
1333  h_ele_EeleOPout->GetYaxis()->SetTitle("Events");
1334  h_ele_vertexX->GetXaxis()->SetTitle("x (cm)");
1335  h_ele_vertexX->GetYaxis()->SetTitle("Events");
1336  h_ele_vertexY->GetXaxis()->SetTitle("y (cm)");
1337  h_ele_vertexY->GetYaxis()->SetTitle("Events");
1338  h_ele_vertexZ->GetXaxis()->SetTitle("z (cm)");
1339  h_ele_vertexZ->GetYaxis()->SetTitle("Events");
1340  h_ele_vertexTIP->GetXaxis()->SetTitle("TIP (cm)");
1341  h_ele_vertexTIP->GetYaxis()->SetTitle("Events");
1342  h_ele_vertexTIPVsEta->GetYaxis()->SetTitle("TIP (cm)");
1343  h_ele_vertexTIPVsEta->GetXaxis()->SetTitle("#eta");
1344  h_ele_vertexTIPVsPhi->GetYaxis()->SetTitle("TIP (cm)");
1345  h_ele_vertexTIPVsPhi->GetXaxis()->SetTitle("#phi (rad)");
1346  h_ele_vertexTIPVsPt->GetYaxis()->SetTitle("TIP (cm)");
1347  h_ele_vertexTIPVsEta->GetXaxis()->SetTitle("p_{T} (GeV/c)");
1348  h_ele_dEtaSc_propVtx->GetXaxis()->SetTitle("#eta_{sc} - #eta_{tr}");
1349  h_ele_dEtaSc_propVtx->GetYaxis()->SetTitle("Events");
1350  h_ele_dEtaCl_propOut->GetXaxis()->SetTitle("#eta_{seedcl} - #eta_{tr}");
1351  h_ele_dEtaCl_propOut->GetYaxis()->SetTitle("Events");
1352  h_ele_dEtaEleCl_propOut->GetXaxis()->SetTitle("#eta_{elecl} - #eta_{tr}");
1353  h_ele_dEtaEleCl_propOut->GetYaxis()->SetTitle("Events");
1354  h_ele_dPhiSc_propVtx->GetXaxis()->SetTitle("#phi_{sc} - #phi_{tr} (rad)");
1355  h_ele_dPhiSc_propVtx->GetYaxis()->SetTitle("Events");
1356  h_ele_dPhiCl_propOut->GetXaxis()->SetTitle("#phi_{seedcl} - #phi_{tr} (rad)");
1357  h_ele_dPhiCl_propOut->GetYaxis()->SetTitle("Events");
1358  h_ele_dPhiEleCl_propOut->GetXaxis()->SetTitle("#phi_{elecl} - #phi_{tr} (rad)");
1359  h_ele_dPhiEleCl_propOut->GetYaxis()->SetTitle("Events");
1360  h_ele_HoE->GetXaxis()->SetTitle("H/E");
1361  h_ele_HoE->GetYaxis()->SetTitle("Events");
1362  h_ele_HoE_fiducial->GetXaxis()->SetTitle("H/E");
1363  h_ele_HoE_fiducial->GetYaxis()->SetTitle("Events");
1364  h_ele_fbrem->GetXaxis()->SetTitle("P_{in} - P_{out} / P_{in}");
1365  h_ele_fbrem->GetYaxis()->SetTitle("Events");
1366  h_ele_seed_dphi2_->GetXaxis()->SetTitle("#phi_{hit}-#phi_{pred} (rad)");
1367  h_ele_seed_dphi2_->GetYaxis()->SetTitle("Events");
1368  h_ele_seed_drz2_->GetXaxis()->SetTitle("r(z)_{hit}-r(z)_{pred} (cm)");
1369  h_ele_seed_drz2_->GetYaxis()->SetTitle("Events");
1370  h_ele_seed_subdet2_->GetXaxis()->SetTitle("2nd hit subdet Id");
1371  h_ele_seed_subdet2_->GetYaxis()->SetTitle("Events");
1372  h_ele_classes->GetXaxis()->SetTitle("class Id");
1373  h_ele_classes->GetYaxis()->SetTitle("Events");
1374  h_ele_mee_all->GetXaxis()->SetTitle("m_{ee} (GeV/c^{2})");
1375  h_ele_mee_all->GetYaxis()->SetTitle("Events");
1376  h_ele_mee_os->GetXaxis()->SetTitle("m_{e^{+}e^{-}} (GeV/c^{2})");
1377  h_ele_mee_os->GetYaxis()->SetTitle("Events");
1378  h_ele_mee_os_ebeb->GetXaxis()->SetTitle("m_{e^{+}e^{-}} (GeV/c^{2})");
1379  h_ele_mee_os_ebeb->GetYaxis()->SetTitle("Events");
1380  h_ele_mee_os_ebee->GetXaxis()->SetTitle("m_{e^{+}e^{-}} (GeV/c^{2})");
1381  h_ele_mee_os_ebee->GetYaxis()->SetTitle("Events");
1382  h_ele_mee_os_eeee->GetXaxis()->SetTitle("m_{e^{+}e^{-}} (GeV/c^{2})");
1383  h_ele_mee_os_eeee->GetYaxis()->SetTitle("Events");
1384  h_ele_mee_os_gg->GetXaxis()->SetTitle("m_{e^{+}e^{-}} (GeV/c^{2})");
1385  h_ele_mee_os_gg->GetYaxis()->SetTitle("Events");
1386  h_ele_mee_os_gb->GetXaxis()->SetTitle("m_{e^{+}e^{-}} (GeV/c^{2})");
1387  h_ele_mee_os_gb->GetYaxis()->SetTitle("Events");
1388  h_ele_mee_os_bb->GetXaxis()->SetTitle("m_{e^{+}e^{-}} (GeV/c^{2})");
1389  h_ele_mee_os_bb->GetYaxis()->SetTitle("Events");
1390  h_ele_E2mnE1vsMee_all->GetXaxis()->SetTitle("m_{e^{+}e^{-}} (GeV/c^{2})");
1391  h_ele_E2mnE1vsMee_all->GetYaxis()->SetTitle("E2 - E1 (GeV)");
1392  h_ele_E2mnE1vsMee_egeg_all->GetXaxis()->SetTitle("m_{e^{+}e^{-}} (GeV/c^{2})");
1393  h_ele_E2mnE1vsMee_egeg_all->GetYaxis()->SetTitle("E2 - E1 (GeV)");
1394  histNum_->GetXaxis()->SetTitle("N_{ele}");
1395  histNum_->GetYaxis()->SetTitle("Events");
1396  h_ele_fbremVsEta_mode->GetXaxis()->SetTitle("#eta");
1397  h_ele_fbremVsEta_mean->GetXaxis()->SetTitle("#eta");
1398 }
1399 
1400 void GsfElectronDataAnalyzer::endJob() {
1401  histfile_->cd();
1402  std::cout << "efficiency calculation " << std::endl;
1403  // efficiency vs eta
1404  TH1F *h_ele_etaEff = (TH1F *)h_ele_matchingObjectEta_matched->Clone("h_ele_etaEff");
1405  h_ele_etaEff->Reset();
1406  h_ele_etaEff->Divide(h_ele_matchingObjectEta_matched, h_matchingObjectEta, 1, 1, "b");
1407  h_ele_etaEff->Print();
1408  h_ele_etaEff->GetXaxis()->SetTitle("#eta");
1409  h_ele_etaEff->GetYaxis()->SetTitle("Efficiency");
1410 
1411  // efficiency vs z
1412  TH1F *h_ele_zEff = (TH1F *)h_ele_matchingObjectZ_matched->Clone("h_ele_zEff");
1413  h_ele_zEff->Reset();
1414  h_ele_zEff->Divide(h_ele_matchingObjectZ_matched, h_matchingObjectZ, 1, 1, "b");
1415  h_ele_zEff->Print();
1416  h_ele_zEff->GetXaxis()->SetTitle("z (cm)");
1417  h_ele_zEff->GetYaxis()->SetTitle("Efficiency");
1418 
1419  // efficiency vs |eta|
1420  TH1F *h_ele_absetaEff = (TH1F *)h_ele_matchingObjectAbsEta_matched->Clone("h_ele_absetaEff");
1421  h_ele_absetaEff->Reset();
1422  h_ele_absetaEff->Divide(h_ele_matchingObjectAbsEta_matched, h_matchingObjectAbsEta, 1, 1, "b");
1423  h_ele_absetaEff->GetXaxis()->SetTitle("|#eta|");
1424  h_ele_absetaEff->GetYaxis()->SetTitle("Efficiency");
1425 
1426  // efficiency vs pt
1427  TH1F *h_ele_ptEff = (TH1F *)h_ele_matchingObjectPt_matched->Clone("h_ele_ptEff");
1428  h_ele_ptEff->Reset();
1429  h_ele_ptEff->Divide(h_ele_matchingObjectPt_matched, h_matchingObjectPt, 1, 1, "b");
1430  h_ele_ptEff->GetXaxis()->SetTitle("p_{T} (GeV/c)");
1431  h_ele_ptEff->GetYaxis()->SetTitle("Efficiency");
1432 
1433  // efficiency vs phi
1434  TH1F *h_ele_phiEff = (TH1F *)h_ele_matchingObjectPhi_matched->Clone("h_ele_phiEff");
1435  h_ele_phiEff->Reset();
1436  h_ele_phiEff->Divide(h_ele_matchingObjectPhi_matched, h_matchingObjectPhi, 1, 1, "b");
1437  h_ele_phiEff->GetXaxis()->SetTitle("#phi (rad)");
1438  h_ele_phiEff->GetYaxis()->SetTitle("Efficiency");
1439 
1440  // classes
1441  TH1F *h_ele_eta_goldenFrac = (TH1F *)h_ele_eta_golden->Clone("h_ele_eta_goldenFrac");
1442  h_ele_eta_goldenFrac->Reset();
1443  h_ele_eta_goldenFrac->Divide(h_ele_eta_golden, h_ele_eta, 1, 1);
1444  h_ele_eta_goldenFrac->GetXaxis()->SetTitle("|#eta|");
1445  h_ele_eta_goldenFrac->GetYaxis()->SetTitle("Fraction of electrons");
1446  h_ele_eta_goldenFrac->SetTitle("fraction of golden electrons vs eta");
1447  TH1F *h_ele_eta_bbremFrac = (TH1F *)h_ele_eta_bbrem->Clone("h_ele_eta_bbremFrac");
1448  h_ele_eta_bbremFrac->Reset();
1449  h_ele_eta_bbremFrac->Divide(h_ele_eta_bbrem, h_ele_eta, 1, 1);
1450  h_ele_eta_bbremFrac->GetXaxis()->SetTitle("|#eta|");
1451  h_ele_eta_bbremFrac->GetYaxis()->SetTitle("Fraction of electrons");
1452  h_ele_eta_bbremFrac->SetTitle("fraction of big brem electrons vs eta");
1453  TH1F *h_ele_eta_narrowFrac = (TH1F *)h_ele_eta_narrow->Clone("h_ele_eta_narrowFrac");
1454  h_ele_eta_narrowFrac->Reset();
1455  h_ele_eta_narrowFrac->Divide(h_ele_eta_narrow, h_ele_eta, 1, 1);
1456  h_ele_eta_narrowFrac->GetXaxis()->SetTitle("|#eta|");
1457  h_ele_eta_narrowFrac->GetYaxis()->SetTitle("Fraction of electrons");
1458  h_ele_eta_narrowFrac->SetTitle("fraction of narrow electrons vs eta");
1459  TH1F *h_ele_eta_showerFrac = (TH1F *)h_ele_eta_shower->Clone("h_ele_eta_showerFrac");
1460  h_ele_eta_showerFrac->Reset();
1461  h_ele_eta_showerFrac->Divide(h_ele_eta_shower, h_ele_eta, 1, 1);
1462  h_ele_eta_showerFrac->GetXaxis()->SetTitle("|#eta|");
1463  h_ele_eta_showerFrac->GetYaxis()->SetTitle("Fraction of electrons");
1464  h_ele_eta_showerFrac->SetTitle("fraction of showering electrons vs eta");
1465 
1466  // fbrem
1467  TH1F *h_ele_xOverX0VsEta = new TH1F("h_ele_xOverx0VsEta", "mean X/X_0 vs eta", nbineta / 2, 0.0, 2.5);
1468  for (int ibin = 1; ibin < h_ele_fbremVsEta_mean->GetNbinsX() + 1; ibin++) {
1469  double xOverX0 = 0.;
1470  if (h_ele_fbremVsEta_mean->GetBinContent(ibin) > 0.)
1471  xOverX0 = -log(h_ele_fbremVsEta_mean->GetBinContent(ibin));
1472  h_ele_xOverX0VsEta->SetBinContent(ibin, xOverX0);
1473  }
1474 
1475  //profiles from 2D histos
1476  TProfile *p_ele_PoPmatchingObjectVsEta_matched = h_ele_PoPmatchingObjectVsEta_matched->ProfileX();
1477  p_ele_PoPmatchingObjectVsEta_matched->SetTitle("mean ele momentum / matching SC energy vs eta");
1478  p_ele_PoPmatchingObjectVsEta_matched->GetXaxis()->SetTitle("#eta");
1479  p_ele_PoPmatchingObjectVsEta_matched->GetYaxis()->SetTitle("<P/E_{matching SC}>");
1480  p_ele_PoPmatchingObjectVsEta_matched->Write();
1481  TProfile *p_ele_PoPmatchingObjectVsPhi_matched = h_ele_PoPmatchingObjectVsPhi_matched->ProfileX();
1482  p_ele_PoPmatchingObjectVsPhi_matched->SetTitle("mean ele momentum / gen momentum vs phi");
1483  p_ele_PoPmatchingObjectVsPhi_matched->GetXaxis()->SetTitle("#phi (rad)");
1484  p_ele_PoPmatchingObjectVsPhi_matched->GetYaxis()->SetTitle("<P/E_{matching SC}>");
1485  p_ele_PoPmatchingObjectVsPhi_matched->Write();
1486  TProfile *p_ele_EtaMnEtamatchingObjectVsEta_matched = h_ele_EtaMnEtamatchingObjectVsEta_matched->ProfileX();
1487  p_ele_EtaMnEtamatchingObjectVsEta_matched->GetXaxis()->SetTitle("#eta");
1488  p_ele_EtaMnEtamatchingObjectVsEta_matched->GetYaxis()->SetTitle("<#eta_{rec} - #eta_{matching SC}>");
1489  p_ele_EtaMnEtamatchingObjectVsEta_matched->Write();
1490  TProfile *p_ele_EtaMnEtamatchingObjectVsPhi_matched = h_ele_EtaMnEtamatchingObjectVsPhi_matched->ProfileX();
1491  p_ele_EtaMnEtamatchingObjectVsPhi_matched->GetXaxis()->SetTitle("#phi");
1492  p_ele_EtaMnEtamatchingObjectVsPhi_matched->GetYaxis()->SetTitle("<#eta_{rec} - #eta_{matching SC}>");
1493  p_ele_EtaMnEtamatchingObjectVsPhi_matched->Write();
1494  TProfile *p_ele_PhiMnPhimatchingObjectVsEta_matched = h_ele_PhiMnPhimatchingObjectVsEta_matched->ProfileX();
1495  p_ele_PhiMnPhimatchingObjectVsEta_matched->GetXaxis()->SetTitle("#eta");
1496  p_ele_PhiMnPhimatchingObjectVsEta_matched->GetYaxis()->SetTitle("<#phi_{rec} - #phi_{matching SC}> (rad)");
1497  p_ele_PhiMnPhimatchingObjectVsEta_matched->Write();
1498  TProfile *p_ele_PhiMnPhimatchingObjectVsPhi_matched = h_ele_PhiMnPhimatchingObjectVsPhi_matched->ProfileX();
1499  p_ele_PhiMnPhimatchingObjectVsPhi_matched->GetXaxis()->SetTitle("#phi");
1500  p_ele_PhiMnPhimatchingObjectVsPhi_matched->GetYaxis()->SetTitle("<#phi_{rec} - #phi_{matching SC}> (rad)");
1501  p_ele_PhiMnPhimatchingObjectVsPhi_matched->Write();
1502  TProfile *p_ele_vertexPtVsEta = h_ele_vertexPtVsEta->ProfileX();
1503  p_ele_vertexPtVsEta->GetXaxis()->SetTitle("#eta");
1504  p_ele_vertexPtVsEta->GetYaxis()->SetTitle("<p_{T}> (GeV/c)");
1505  p_ele_vertexPtVsEta->Write();
1506  TProfile *p_ele_vertexPtVsPhi = h_ele_vertexPtVsPhi->ProfileX();
1507  p_ele_vertexPtVsPhi->GetXaxis()->SetTitle("#phi (rad)");
1508  p_ele_vertexPtVsPhi->GetYaxis()->SetTitle("<p_{T}> (GeV/c)");
1509  p_ele_vertexPtVsPhi->Write();
1510  TProfile *p_ele_EoPVsEta = h_ele_EoPVsEta->ProfileX();
1511  p_ele_EoPVsEta->GetXaxis()->SetTitle("#eta");
1512  p_ele_EoPVsEta->GetYaxis()->SetTitle("<E/P_{vertex}>");
1513  p_ele_EoPVsEta->Write();
1514  TProfile *p_ele_EoPVsPhi = h_ele_EoPVsPhi->ProfileX();
1515  p_ele_EoPVsPhi->GetXaxis()->SetTitle("#phi (rad)");
1516  p_ele_EoPVsPhi->GetYaxis()->SetTitle("<E/P_{vertex}>");
1517  p_ele_EoPVsPhi->Write();
1518  TProfile *p_ele_EoPoutVsEta = h_ele_EoPoutVsEta->ProfileX();
1519  p_ele_EoPoutVsEta->GetXaxis()->SetTitle("#eta");
1520  p_ele_EoPoutVsEta->GetYaxis()->SetTitle("<E_{seed}/P_{out}>");
1521  p_ele_EoPoutVsEta->Write();
1522  TProfile *p_ele_EoPoutVsPhi = h_ele_EoPoutVsPhi->ProfileX();
1523  p_ele_EoPoutVsPhi->GetXaxis()->SetTitle("#phi (rad)");
1524  p_ele_EoPoutVsPhi->GetYaxis()->SetTitle("<E_{seed}/P_{out}>");
1525  p_ele_EoPoutVsPhi->Write();
1526  TProfile *p_ele_EeleOPoutVsEta = h_ele_EeleOPoutVsEta->ProfileX();
1527  p_ele_EeleOPoutVsEta->SetTitle("mean ele Eele/pout vs eta");
1528  p_ele_EeleOPoutVsEta->GetXaxis()->SetTitle("#eta");
1529  p_ele_EeleOPoutVsEta->GetYaxis()->SetTitle("<E_{ele}/P_{out}>");
1530  p_ele_EeleOPoutVsEta->Write();
1531  TProfile *p_ele_EeleOPoutVsPhi = h_ele_EeleOPoutVsPhi->ProfileX();
1532  p_ele_EeleOPoutVsPhi->SetTitle("mean ele Eele/pout vs phi");
1533  p_ele_EeleOPoutVsPhi->GetXaxis()->SetTitle("#phi (rad)");
1534  p_ele_EeleOPoutVsPhi->GetYaxis()->SetTitle("<E_{ele}/P_{out}>");
1535  p_ele_EeleOPoutVsPhi->Write();
1536  TProfile *p_ele_HoEVsEta = h_ele_HoEVsEta->ProfileX();
1537  p_ele_HoEVsEta->GetXaxis()->SetTitle("#eta");
1538  p_ele_HoEVsEta->GetYaxis()->SetTitle("<H/E>");
1539  p_ele_HoEVsEta->Write();
1540  TProfile *p_ele_HoEVsPhi = h_ele_HoEVsPhi->ProfileX();
1541  p_ele_HoEVsPhi->GetXaxis()->SetTitle("#phi (rad)");
1542  p_ele_HoEVsPhi->GetYaxis()->SetTitle("<H/E>");
1543  p_ele_HoEVsPhi->Write();
1544  TProfile *p_ele_chi2VsEta = h_ele_chi2VsEta->ProfileX();
1545  p_ele_chi2VsEta->GetXaxis()->SetTitle("#eta");
1546  p_ele_chi2VsEta->GetYaxis()->SetTitle("<#Chi^{2}>");
1547  p_ele_chi2VsEta->Write();
1548  TProfile *p_ele_chi2VsPhi = h_ele_chi2VsPhi->ProfileX();
1549  p_ele_chi2VsPhi->GetXaxis()->SetTitle("#phi (rad)");
1550  p_ele_chi2VsPhi->GetYaxis()->SetTitle("<#Chi^{2}>");
1551  p_ele_chi2VsPhi->Write();
1552  TProfile *p_ele_foundHitsVsEta = h_ele_foundHitsVsEta->ProfileX();
1553  p_ele_foundHitsVsEta->GetXaxis()->SetTitle("#eta");
1554  p_ele_foundHitsVsEta->GetYaxis()->SetTitle("<N_{hits}>");
1555  p_ele_foundHitsVsEta->Write();
1556  TProfile *p_ele_foundHitsVsPhi = h_ele_foundHitsVsPhi->ProfileX();
1557  p_ele_foundHitsVsPhi->GetXaxis()->SetTitle("#phi (rad)");
1558  p_ele_foundHitsVsPhi->GetYaxis()->SetTitle("<N_{hits}>");
1559  p_ele_foundHitsVsPhi->Write();
1560  TProfile *p_ele_lostHitsVsEta = h_ele_lostHitsVsEta->ProfileX();
1561  p_ele_lostHitsVsEta->GetXaxis()->SetTitle("#eta");
1562  p_ele_lostHitsVsEta->GetYaxis()->SetTitle("<N_{hits}>");
1563  p_ele_lostHitsVsEta->Write();
1564  TProfile *p_ele_lostHitsVsPhi = h_ele_lostHitsVsPhi->ProfileX();
1565  p_ele_lostHitsVsPhi->GetXaxis()->SetTitle("#phi (rad)");
1566  p_ele_lostHitsVsPhi->GetYaxis()->SetTitle("<N_{hits}>");
1567  p_ele_lostHitsVsPhi->Write();
1568  TProfile *p_ele_vertexTIPVsEta = h_ele_vertexTIPVsEta->ProfileX();
1569  p_ele_vertexTIPVsEta->SetTitle("mean tip (wrt gen vtx) vs eta");
1570  p_ele_vertexTIPVsEta->GetXaxis()->SetTitle("#eta");
1571  p_ele_vertexTIPVsEta->GetYaxis()->SetTitle("<TIP> (cm)");
1572  p_ele_vertexTIPVsEta->Write();
1573  TProfile *p_ele_vertexTIPVsPhi = h_ele_vertexTIPVsPhi->ProfileX();
1574  p_ele_vertexTIPVsPhi->SetTitle("mean tip (wrt gen vtx) vs phi");
1575  p_ele_vertexTIPVsPhi->GetXaxis()->SetTitle("#phi");
1576  p_ele_vertexTIPVsPhi->GetYaxis()->SetTitle("<TIP> (cm)");
1577  p_ele_vertexTIPVsPhi->Write();
1578  TProfile *p_ele_vertexTIPVsPt = h_ele_vertexTIPVsPt->ProfileX();
1579  p_ele_vertexTIPVsPt->SetTitle("mean tip (wrt gen vtx) vs phi");
1580  p_ele_vertexTIPVsPt->GetXaxis()->SetTitle("p_{T} (GeV/c)");
1581  p_ele_vertexTIPVsPt->GetYaxis()->SetTitle("<TIP> (cm)");
1582  p_ele_vertexTIPVsPt->Write();
1583 
1584  // mc truth
1585 
1586  h_matchingObjectNum->Write();
1587 
1588  // rec event
1589 
1590  histNum_->Write();
1591 
1592  // mc
1593  h_matchingObjectEta->Write();
1594  h_matchingObjectAbsEta->Write();
1595  h_matchingObjectP->Write();
1596  h_matchingObjectPt->Write();
1597  h_matchingObjectPhi->Write();
1598  h_matchingObjectZ->Write();
1599 
1600  h_ele_mee_all->Write();
1601  h_ele_mee_os->Write();
1602  h_ele_mee_os_ebeb->Write();
1603  h_ele_mee_os_ebee->Write();
1604  h_ele_mee_os_eeee->Write();
1605  h_ele_mee_os_gg->Write();
1606  h_ele_mee_os_gb->Write();
1607  h_ele_mee_os_bb->Write();
1608  h_ele_E2mnE1vsMee_all->Write();
1609  h_ele_E2mnE1vsMee_egeg_all->Write();
1610 
1611  // matched electrons
1612  h_ele_charge->Write();
1613  h_ele_chargeVsEta->Write();
1614  h_ele_chargeVsPhi->Write();
1615  h_ele_chargeVsPt->Write();
1616  h_ele_vertexP->Write();
1617  h_ele_vertexPt->Write();
1618  h_ele_Et->Write();
1619  h_ele_vertexPtVsEta->Write();
1620  h_ele_vertexPtVsPhi->Write();
1621  h_ele_matchingObjectPt_matched->Write();
1622  h_ele_vertexEta->Write();
1623  h_ele_vertexEtaVsPhi->Write();
1624  h_ele_matchingObjectAbsEta_matched->Write();
1625  h_ele_matchingObjectEta_matched->Write();
1626  h_ele_matchingObjectPhi_matched->Write();
1627  h_ele_vertexPhi->Write();
1628  h_ele_vertexX->Write();
1629  h_ele_vertexY->Write();
1630  h_ele_vertexZ->Write();
1631  h_ele_vertexTIP->Write();
1632  h_ele_matchingObjectZ_matched->Write();
1633  h_ele_vertexTIPVsEta->Write();
1634  h_ele_vertexTIPVsPhi->Write();
1635  h_ele_vertexTIPVsPt->Write();
1636  h_ele_PoPmatchingObject_matched->Write();
1637  h_ele_PtoPtmatchingObject_matched->Write();
1638  h_ele_PoPmatchingObjectVsEta_matched->Write();
1639  h_ele_PoPmatchingObjectVsPhi_matched->Write();
1640  h_ele_PoPmatchingObjectVsPt_matched->Write();
1641  h_ele_PoPmatchingObject_barrel_matched->Write();
1642  h_ele_PoPmatchingObject_endcaps_matched->Write();
1643  h_ele_PtoPtmatchingObject_barrel_matched->Write();
1644  h_ele_PtoPtmatchingObject_endcaps_matched->Write();
1645  h_ele_EtaMnEtamatchingObject_matched->Write();
1646  h_ele_EtaMnEtamatchingObjectVsEta_matched->Write();
1647  h_ele_EtaMnEtamatchingObjectVsPhi_matched->Write();
1648  h_ele_EtaMnEtamatchingObjectVsPt_matched->Write();
1649  h_ele_PhiMnPhimatchingObject_matched->Write();
1650  h_ele_PhiMnPhimatchingObject2_matched->Write();
1651  h_ele_PhiMnPhimatchingObjectVsEta_matched->Write();
1652  h_ele_PhiMnPhimatchingObjectVsPhi_matched->Write();
1653  h_ele_PhiMnPhimatchingObjectVsPt_matched->Write();
1654 
1655  // matched electron, superclusters
1656  histSclEn_->Write();
1657  histSclEoEmatchingObject_barrel_matched->Write();
1658  histSclEoEmatchingObject_endcaps_matched->Write();
1659  histSclEoEmatchingObject_barrel_new_matched->Write();
1660  histSclEoEmatchingObject_endcaps_new_matched->Write();
1661  histSclEt_->Write();
1662  histSclEtVsEta_->Write();
1663  histSclEtVsPhi_->Write();
1664  histSclEtaVsPhi_->Write();
1665  histSclEta_->Write();
1666  histSclPhi_->Write();
1667  histSclSigEtaEta_->Write();
1668  histSclSigIEtaIEta_barrel_->Write();
1669  histSclSigIEtaIEta_endcaps_->Write();
1670  histSclE1x5_->Write();
1671  histSclE1x5_barrel_->Write();
1672  histSclE1x5_endcaps_->Write();
1673  histSclE2x5max_->Write();
1674  histSclE2x5max_barrel_->Write();
1675  histSclE2x5max_endcaps_->Write();
1676  histSclE5x5_->Write();
1677  histSclE5x5_barrel_->Write();
1678  histSclE5x5_endcaps_->Write();
1679 
1680  // matched electron, gsf tracks
1681  h_ele_ambiguousTracks->Write();
1682  h_ele_ambiguousTracksVsEta->Write();
1683  h_ele_ambiguousTracksVsPhi->Write();
1684  h_ele_ambiguousTracksVsPt->Write();
1685  h_ele_foundHits->Write();
1686  h_ele_foundHitsVsEta->Write();
1687  h_ele_foundHitsVsPhi->Write();
1688  h_ele_foundHitsVsPt->Write();
1689  h_ele_lostHits->Write();
1690  h_ele_lostHitsVsEta->Write();
1691  h_ele_lostHitsVsPhi->Write();
1692  h_ele_lostHitsVsPt->Write();
1693  h_ele_chi2->Write();
1694  h_ele_chi2VsEta->Write();
1695  h_ele_chi2VsPhi->Write();
1696  h_ele_chi2VsPt->Write();
1697  h_ele_PinMnPout->Write();
1698  h_ele_PinMnPout_mode->Write();
1699  h_ele_PinMnPoutVsEta_mode->Write();
1700  h_ele_PinMnPoutVsPhi_mode->Write();
1701  h_ele_PinMnPoutVsPt_mode->Write();
1702  h_ele_PinMnPoutVsE_mode->Write();
1703  h_ele_PinMnPoutVsChi2_mode->Write();
1704  h_ele_outerP->Write();
1705  h_ele_outerP_mode->Write();
1706  h_ele_outerPVsEta_mode->Write();
1707  h_ele_outerPt->Write();
1708  h_ele_outerPt_mode->Write();
1709  h_ele_outerPtVsEta_mode->Write();
1710  h_ele_outerPtVsPhi_mode->Write();
1711  h_ele_outerPtVsPt_mode->Write();
1712 
1713  // matched electrons, matching
1714  h_ele_EoP->Write();
1715  h_ele_EoPVsEta->Write();
1716  h_ele_EoPVsPhi->Write();
1717  h_ele_EoPVsE->Write();
1718  h_ele_EseedOP->Write();
1719  h_ele_EseedOPVsEta->Write();
1720  h_ele_EseedOPVsPhi->Write();
1721  h_ele_EseedOPVsE->Write();
1722  h_ele_EoPout->Write();
1723  h_ele_EoPoutVsEta->Write();
1724  h_ele_EoPoutVsPhi->Write();
1725  h_ele_EoPoutVsE->Write();
1726  h_ele_EeleOPout->Write();
1727  h_ele_EeleOPoutVsEta->Write();
1728  h_ele_EeleOPoutVsPhi->Write();
1729  h_ele_EeleOPoutVsE->Write();
1730  h_ele_dEtaSc_propVtx->Write();
1731  h_ele_dEtaScVsEta_propVtx->Write();
1732  h_ele_dEtaScVsPhi_propVtx->Write();
1733  h_ele_dEtaScVsPt_propVtx->Write();
1734  h_ele_dPhiSc_propVtx->Write();
1735  h_ele_dPhiScVsEta_propVtx->Write();
1736  h_ele_dPhiScVsPhi_propVtx->Write();
1737  h_ele_dPhiScVsPt_propVtx->Write();
1738  h_ele_dEtaCl_propOut->Write();
1739  h_ele_dEtaClVsEta_propOut->Write();
1740  h_ele_dEtaClVsPhi_propOut->Write();
1741  h_ele_dEtaClVsPt_propOut->Write();
1742  h_ele_dPhiCl_propOut->Write();
1743  h_ele_dPhiClVsEta_propOut->Write();
1744  h_ele_dPhiClVsPhi_propOut->Write();
1745  h_ele_dPhiClVsPt_propOut->Write();
1746  h_ele_dEtaEleCl_propOut->Write();
1747  h_ele_dEtaEleClVsEta_propOut->Write();
1748  h_ele_dEtaEleClVsPhi_propOut->Write();
1749  h_ele_dEtaEleClVsPt_propOut->Write();
1750  h_ele_dPhiEleCl_propOut->Write();
1751  h_ele_dPhiEleClVsEta_propOut->Write();
1752  h_ele_dPhiEleClVsPhi_propOut->Write();
1753  h_ele_dPhiEleClVsPt_propOut->Write();
1754  h_ele_HoE->Write();
1755  h_ele_HoE_fiducial->Write();
1756  h_ele_HoEVsEta->Write();
1757  h_ele_HoEVsPhi->Write();
1758  h_ele_HoEVsE->Write();
1759 
1760  h_ele_seed_dphi2_->Write();
1761  h_ele_seed_subdet2_->Write();
1762  TProfile *p_ele_seed_dphi2VsEta_ = h_ele_seed_dphi2VsEta_->ProfileX();
1763  p_ele_seed_dphi2VsEta_->SetTitle("mean ele seed dphi 2nd layer vs eta");
1764  p_ele_seed_dphi2VsEta_->GetXaxis()->SetTitle("#eta");
1765  p_ele_seed_dphi2VsEta_->GetYaxis()->SetTitle("<#phi_{pred} - #phi_{hit}, 2nd layer> (rad)");
1766  p_ele_seed_dphi2VsEta_->SetMinimum(-0.004);
1767  p_ele_seed_dphi2VsEta_->SetMaximum(0.004);
1768  p_ele_seed_dphi2VsEta_->Write();
1769  TProfile *p_ele_seed_dphi2VsPt_ = h_ele_seed_dphi2VsPt_->ProfileX();
1770  p_ele_seed_dphi2VsPt_->SetTitle("mean ele seed dphi 2nd layer vs pt");
1771  p_ele_seed_dphi2VsPt_->GetXaxis()->SetTitle("p_{T} (GeV/c)");
1772  p_ele_seed_dphi2VsPt_->GetYaxis()->SetTitle("<#phi_{pred} - #phi_{hit}, 2nd layer> (rad)");
1773  p_ele_seed_dphi2VsPt_->Write();
1774  p_ele_seed_dphi2VsPt_->SetMinimum(-0.004);
1775  p_ele_seed_dphi2VsPt_->SetMaximum(0.004);
1776  h_ele_seed_drz2_->Write();
1777  TProfile *p_ele_seed_drz2VsEta_ = h_ele_seed_drz2VsEta_->ProfileX();
1778  p_ele_seed_drz2VsEta_->SetTitle("mean ele seed dr(dz) 2nd layer vs eta");
1779  p_ele_seed_drz2VsEta_->GetXaxis()->SetTitle("#eta");
1780  p_ele_seed_drz2VsEta_->GetYaxis()->SetTitle("<r(z)_{pred} - r(z)_{hit}, 2nd layer> (cm)");
1781  p_ele_seed_drz2VsEta_->SetMinimum(-0.15);
1782  p_ele_seed_drz2VsEta_->SetMaximum(0.15);
1783  p_ele_seed_drz2VsEta_->Write();
1784  TProfile *p_ele_seed_drz2VsPt_ = h_ele_seed_drz2VsPt_->ProfileX();
1785  p_ele_seed_drz2VsPt_->SetTitle("mean ele seed dr(dz) 2nd layer vs pt");
1786  p_ele_seed_drz2VsPt_->GetXaxis()->SetTitle("p_{T} (GeV/c)");
1787  p_ele_seed_drz2VsPt_->GetYaxis()->SetTitle("<r(z)_{pred} - r(z)_{hit}, 2nd layer> (cm)");
1788  p_ele_seed_drz2VsPt_->SetMinimum(-0.15);
1789  p_ele_seed_drz2VsPt_->SetMaximum(0.15);
1790  p_ele_seed_drz2VsPt_->Write();
1791 
1792  // classes
1793  h_ele_classes->Write();
1794  h_ele_eta->Write();
1795  h_ele_eta_golden->Write();
1796  h_ele_eta_bbrem->Write();
1797  h_ele_eta_narrow->Write();
1798  h_ele_eta_shower->Write();
1799  h_ele_PinVsPoutGolden_mode->Write();
1800  h_ele_PinVsPoutShowering_mode->Write();
1801  h_ele_PinVsPoutGolden_mean->Write();
1802  h_ele_PinVsPoutShowering_mean->Write();
1803  h_ele_PtinVsPtoutGolden_mode->Write();
1804  h_ele_PtinVsPtoutShowering_mode->Write();
1805  h_ele_PtinVsPtoutGolden_mean->Write();
1806  h_ele_PtinVsPtoutShowering_mean->Write();
1807  histSclEoEmatchingObjectGolden_barrel->Write();
1808  histSclEoEmatchingObjectGolden_endcaps->Write();
1809  histSclEoEmatchingObjectShowering_barrel->Write();
1810  histSclEoEmatchingObjectShowering_endcaps->Write();
1811 
1812  // fbrem
1813  h_ele_fbrem->Write();
1814  h_ele_fbremVsEta_mode->Write();
1815  h_ele_fbremVsEta_mean->Write();
1816  h_ele_etaEff->Write();
1817  h_ele_zEff->Write();
1818  h_ele_phiEff->Write();
1819  h_ele_absetaEff->Write();
1820  h_ele_ptEff->Write();
1821  h_ele_eta_goldenFrac->Write();
1822  h_ele_eta_bbremFrac->Write();
1823  h_ele_eta_narrowFrac->Write();
1824  h_ele_eta_showerFrac->Write();
1825  h_ele_xOverX0VsEta->Write();
1826 
1827  // e/g et pflow electrons
1828  h_ele_mva->Write();
1829  h_ele_provenance->Write();
1830 
1831  // isolation
1832  h_ele_tkSumPt_dr03->GetXaxis()->SetTitle("TkIsoSum, cone 0.3 (GeV/c)");
1833  h_ele_tkSumPt_dr03->GetYaxis()->SetTitle("Events");
1834  h_ele_tkSumPt_dr03->Write();
1835  h_ele_ecalRecHitSumEt_dr03->GetXaxis()->SetTitle("EcalIsoSum, cone 0.3 (GeV)");
1836  h_ele_ecalRecHitSumEt_dr03->GetYaxis()->SetTitle("Events");
1837  h_ele_ecalRecHitSumEt_dr03->Write();
1838  h_ele_hcalDepth1TowerSumEt_dr03->GetXaxis()->SetTitle("Hcal1IsoSum, cone 0.3 (GeV)");
1839  h_ele_hcalDepth1TowerSumEt_dr03->GetYaxis()->SetTitle("Events");
1840  h_ele_hcalDepth1TowerSumEt_dr03->Write();
1841  h_ele_hcalDepth2TowerSumEt_dr03->GetXaxis()->SetTitle("Hcal2IsoSum, cone 0.3 (GeV)");
1842  h_ele_hcalDepth2TowerSumEt_dr03->GetYaxis()->SetTitle("Events");
1843  h_ele_hcalDepth2TowerSumEt_dr03->Write();
1844  h_ele_tkSumPt_dr04->GetXaxis()->SetTitle("TkIsoSum, cone 0.4 (GeV/c)");
1845  h_ele_tkSumPt_dr04->GetYaxis()->SetTitle("Events");
1846  h_ele_tkSumPt_dr04->Write();
1847  h_ele_ecalRecHitSumEt_dr04->GetXaxis()->SetTitle("EcalIsoSum, cone 0.4 (GeV)");
1848  h_ele_ecalRecHitSumEt_dr04->GetYaxis()->SetTitle("Events");
1849  h_ele_ecalRecHitSumEt_dr04->Write();
1850  h_ele_hcalDepth1TowerSumEt_dr04->GetXaxis()->SetTitle("Hcal1IsoSum, cone 0.4 (GeV)");
1851  h_ele_hcalDepth1TowerSumEt_dr04->GetYaxis()->SetTitle("Events");
1852  h_ele_hcalDepth1TowerSumEt_dr04->Write();
1853  h_ele_hcalDepth2TowerSumEt_dr04->GetXaxis()->SetTitle("Hcal2IsoSum, cone 0.4 (GeV)");
1854  h_ele_hcalDepth2TowerSumEt_dr04->GetYaxis()->SetTitle("Events");
1855  h_ele_hcalDepth2TowerSumEt_dr04->Write();
1856 }
1857 
1858 void GsfElectronDataAnalyzer::analyze(const edm::Event &iEvent, const edm::EventSetup &iSetup) {
1859  std::cout << "analyzing new event " << std::endl;
1860  nEvents_++;
1861 
1862  // check event pass requested triggers if any
1863  if (!trigger(iEvent))
1864  return;
1865 
1866  std::cout << "new event passing trigger " << std::endl;
1867  nAfterTrigger_++;
1868 
1869  // get reco electrons
1870  edm::Handle<reco::GsfElectronCollection> gsfElectrons;
1871  iEvent.getByLabel(electronCollection_, gsfElectrons);
1872  edm::LogInfo("") << "\n\n =================> Treating event " << iEvent.id() << " Number of electrons "
1873  << gsfElectrons.product()->size();
1874 
1875  // get reco supercluster collection
1876  edm::Handle<reco::SuperClusterCollection> recoClusters;
1877  iEvent.getByLabel(matchingObjectCollection_, recoClusters);
1878 
1879  // get the beamspot from the Event:
1880  edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
1881  iEvent.getByLabel(beamSpot_, recoBeamSpotHandle);
1882  const BeamSpot bs = *recoBeamSpotHandle;
1883 
1884  histNum_->Fill((*gsfElectrons).size());
1885 
1886  // selected rec electrons
1887  for (reco::GsfElectronCollection::const_iterator gsfIter = gsfElectrons->begin(); gsfIter != gsfElectrons->end();
1888  gsfIter++) {
1889  // select electrons
1890  if (gsfIter->superCluster()->energy() / cosh(gsfIter->superCluster()->eta()) < minEt_)
1891  continue;
1892  if (std::abs(gsfIter->eta()) > maxAbsEta_)
1893  continue;
1894  if (gsfIter->pt() < minPt_)
1895  continue;
1896 
1897  if (gsfIter->isEB() && isEE_)
1898  continue;
1899  if (gsfIter->isEE() && isEB_)
1900  continue;
1901  if (gsfIter->isEBEEGap() && isNotEBEEGap_)
1902  continue;
1903 
1904  if (gsfIter->ecalDrivenSeed() && isTrackerDriven_)
1905  continue;
1906  if (gsfIter->trackerDrivenSeed() && isEcalDriven_)
1907  continue;
1908 
1909  if (gsfIter->isEB() && gsfIter->eSuperClusterOverP() < eOverPMinBarrel_)
1910  continue;
1911  if (gsfIter->isEB() && gsfIter->eSuperClusterOverP() > eOverPMaxBarrel_)
1912  continue;
1913  if (gsfIter->isEE() && gsfIter->eSuperClusterOverP() < eOverPMinEndcaps_)
1914  continue;
1915  if (gsfIter->isEE() && gsfIter->eSuperClusterOverP() > eOverPMaxEndcaps_)
1916  continue;
1917  if (gsfIter->isEB() && std::abs(gsfIter->deltaEtaSuperClusterTrackAtVtx()) < dEtaMinBarrel_)
1918  continue;
1919  if (gsfIter->isEB() && std::abs(gsfIter->deltaEtaSuperClusterTrackAtVtx()) > dEtaMaxBarrel_)
1920  continue;
1921  if (gsfIter->isEE() && std::abs(gsfIter->deltaEtaSuperClusterTrackAtVtx()) < dEtaMinEndcaps_)
1922  continue;
1923  if (gsfIter->isEE() && std::abs(gsfIter->deltaEtaSuperClusterTrackAtVtx()) > dEtaMaxEndcaps_)
1924  continue;
1925  if (gsfIter->isEB() && std::abs(gsfIter->deltaPhiSuperClusterTrackAtVtx()) < dPhiMinBarrel_)
1926  continue;
1927  if (gsfIter->isEB() && std::abs(gsfIter->deltaPhiSuperClusterTrackAtVtx()) > dPhiMaxBarrel_)
1928  continue;
1929  if (gsfIter->isEE() && std::abs(gsfIter->deltaPhiSuperClusterTrackAtVtx()) < dPhiMinEndcaps_)
1930  continue;
1931  if (gsfIter->isEE() && std::abs(gsfIter->deltaPhiSuperClusterTrackAtVtx()) > dPhiMaxEndcaps_)
1932  continue;
1933  if (gsfIter->isEB() && gsfIter->scSigmaIEtaIEta() < sigIetaIetaMinBarrel_)
1934  continue;
1935  if (gsfIter->isEB() && gsfIter->scSigmaIEtaIEta() > sigIetaIetaMaxBarrel_)
1936  continue;
1937  if (gsfIter->isEE() && gsfIter->scSigmaIEtaIEta() < sigIetaIetaMinEndcaps_)
1938  continue;
1939  if (gsfIter->isEE() && gsfIter->scSigmaIEtaIEta() > sigIetaIetaMaxEndcaps_)
1940  continue;
1941  if (gsfIter->isEB() && gsfIter->hadronicOverEm() > hadronicOverEmMaxBarrel_)
1942  continue;
1943  if (gsfIter->isEE() && gsfIter->hadronicOverEm() > hadronicOverEmMaxEndcaps_)
1944  continue;
1945  if (gsfIter->mva_e_pi() < mvaMin_)
1946  continue;
1947 
1948  double d = (gsfIter->vertex().x() - bs.position().x()) * (gsfIter->vertex().x() - bs.position().x()) +
1949  (gsfIter->vertex().y() - bs.position().y()) * (gsfIter->vertex().y() - bs.position().y());
1950  d = sqrt(d);
1951  if (gsfIter->isEB() && d > tipMaxBarrel_)
1952  continue;
1953  if (gsfIter->isEE() && d > tipMaxEndcaps_)
1954  continue;
1955 
1956  if (gsfIter->dr03TkSumPt() > tkIso03Max_)
1957  continue;
1958  if (gsfIter->isEB() && gsfIter->dr03HcalDepth1TowerSumEt() > hcalIso03Depth1MaxBarrel_)
1959  continue;
1960  if (gsfIter->isEE() && gsfIter->dr03HcalDepth1TowerSumEt() > hcalIso03Depth1MaxEndcaps_)
1961  continue;
1962  if (gsfIter->isEE() && gsfIter->dr03HcalDepth2TowerSumEt() > hcalIso03Depth2MaxEndcaps_)
1963  continue;
1964  if (gsfIter->isEB() && gsfIter->dr03EcalRecHitSumEt() > ecalIso03MaxBarrel_)
1965  continue;
1966  if (gsfIter->isEE() && gsfIter->dr03EcalRecHitSumEt() > ecalIso03MaxEndcaps_)
1967  continue;
1968 
1969  // electron related distributions
1970  h_ele_charge->Fill(gsfIter->charge());
1971  h_ele_chargeVsEta->Fill(gsfIter->eta(), gsfIter->charge());
1972  h_ele_chargeVsPhi->Fill(gsfIter->phi(), gsfIter->charge());
1973  h_ele_chargeVsPt->Fill(gsfIter->pt(), gsfIter->charge());
1974  h_ele_vertexP->Fill(gsfIter->p());
1975  h_ele_vertexPt->Fill(gsfIter->pt());
1976  h_ele_Et->Fill(gsfIter->superCluster()->energy() / cosh(gsfIter->superCluster()->eta()));
1977  h_ele_vertexPtVsEta->Fill(gsfIter->eta(), gsfIter->pt());
1978  h_ele_vertexPtVsPhi->Fill(gsfIter->phi(), gsfIter->pt());
1979  h_ele_vertexEta->Fill(gsfIter->eta());
1980 
1981  h_ele_vertexEtaVsPhi->Fill(gsfIter->phi(), gsfIter->eta());
1982  h_ele_vertexPhi->Fill(gsfIter->phi());
1983  h_ele_vertexX->Fill(gsfIter->vertex().x());
1984  h_ele_vertexY->Fill(gsfIter->vertex().y());
1985  h_ele_vertexZ->Fill(gsfIter->vertex().z());
1986  h_ele_vertexTIP->Fill(d);
1987  h_ele_vertexTIPVsEta->Fill(gsfIter->eta(), d);
1988  h_ele_vertexTIPVsPhi->Fill(gsfIter->phi(), d);
1989  h_ele_vertexTIPVsPt->Fill(gsfIter->pt(), d);
1990 
1991  // supercluster related distributions
1992  reco::SuperClusterRef sclRef = gsfIter->superCluster();
1993  if (!gsfIter->ecalDrivenSeed() && gsfIter->trackerDrivenSeed())
1994  sclRef = gsfIter->parentSuperCluster();
1995  histSclEn_->Fill(sclRef->energy());
1996  double R = TMath::Sqrt(sclRef->x() * sclRef->x() + sclRef->y() * sclRef->y() + sclRef->z() * sclRef->z());
1997  double Rt = TMath::Sqrt(sclRef->x() * sclRef->x() + sclRef->y() * sclRef->y());
1998  histSclEt_->Fill(sclRef->energy() * (Rt / R));
1999  histSclEtVsEta_->Fill(sclRef->eta(), sclRef->energy() * (Rt / R));
2000  histSclEtVsPhi_->Fill(sclRef->phi(), sclRef->energy() * (Rt / R));
2001  histSclEta_->Fill(sclRef->eta());
2002  histSclEtaVsPhi_->Fill(sclRef->phi(), sclRef->eta());
2003  histSclPhi_->Fill(sclRef->phi());
2004  histSclSigEtaEta_->Fill(gsfIter->scSigmaEtaEta());
2005  if (gsfIter->isEB())
2006  histSclSigIEtaIEta_barrel_->Fill(gsfIter->scSigmaIEtaIEta());
2007  if (gsfIter->isEE())
2008  histSclSigIEtaIEta_endcaps_->Fill(gsfIter->scSigmaIEtaIEta());
2009  histSclE1x5_->Fill(gsfIter->scE1x5());
2010  if (gsfIter->isEB())
2011  histSclE1x5_barrel_->Fill(gsfIter->scE1x5());
2012  if (gsfIter->isEE())
2013  histSclE1x5_endcaps_->Fill(gsfIter->scE1x5());
2014  histSclE2x5max_->Fill(gsfIter->scE2x5Max());
2015  if (gsfIter->isEB())
2016  histSclE2x5max_barrel_->Fill(gsfIter->scE2x5Max());
2017  if (gsfIter->isEE())
2018  histSclE2x5max_endcaps_->Fill(gsfIter->scE2x5Max());
2019  histSclE5x5_->Fill(gsfIter->scE5x5());
2020  if (gsfIter->isEB())
2021  histSclE5x5_barrel_->Fill(gsfIter->scE5x5());
2022  if (gsfIter->isEE())
2023  histSclE5x5_endcaps_->Fill(gsfIter->scE5x5());
2024 
2025  // track related distributions
2026  h_ele_ambiguousTracks->Fill(gsfIter->ambiguousGsfTracksSize());
2027  h_ele_ambiguousTracksVsEta->Fill(gsfIter->eta(), gsfIter->ambiguousGsfTracksSize());
2028  h_ele_ambiguousTracksVsPhi->Fill(gsfIter->phi(), gsfIter->ambiguousGsfTracksSize());
2029  h_ele_ambiguousTracksVsPt->Fill(gsfIter->pt(), gsfIter->ambiguousGsfTracksSize());
2030  if (!readAOD_) { // track extra does not exist in AOD
2031  h_ele_foundHits->Fill(gsfIter->gsfTrack()->numberOfValidHits());
2032  h_ele_foundHitsVsEta->Fill(gsfIter->eta(), gsfIter->gsfTrack()->numberOfValidHits());
2033  h_ele_foundHitsVsPhi->Fill(gsfIter->phi(), gsfIter->gsfTrack()->numberOfValidHits());
2034  h_ele_foundHitsVsPt->Fill(gsfIter->pt(), gsfIter->gsfTrack()->numberOfValidHits());
2035  h_ele_lostHits->Fill(gsfIter->gsfTrack()->numberOfLostHits());
2036  h_ele_lostHitsVsEta->Fill(gsfIter->eta(), gsfIter->gsfTrack()->numberOfLostHits());
2037  h_ele_lostHitsVsPhi->Fill(gsfIter->phi(), gsfIter->gsfTrack()->numberOfLostHits());
2038  h_ele_lostHitsVsPt->Fill(gsfIter->pt(), gsfIter->gsfTrack()->numberOfLostHits());
2039  h_ele_chi2->Fill(gsfIter->gsfTrack()->normalizedChi2());
2040  h_ele_chi2VsEta->Fill(gsfIter->eta(), gsfIter->gsfTrack()->normalizedChi2());
2041  h_ele_chi2VsPhi->Fill(gsfIter->phi(), gsfIter->gsfTrack()->normalizedChi2());
2042  h_ele_chi2VsPt->Fill(gsfIter->pt(), gsfIter->gsfTrack()->normalizedChi2());
2043  }
2044  // from gsf track interface, hence using mean
2045  if (!readAOD_) { // track extra does not exist in AOD
2046  h_ele_PinMnPout->Fill(gsfIter->gsfTrack()->innerMomentum().R() - gsfIter->gsfTrack()->outerMomentum().R());
2047  h_ele_outerP->Fill(gsfIter->gsfTrack()->outerMomentum().R());
2048  h_ele_outerPt->Fill(gsfIter->gsfTrack()->outerMomentum().Rho());
2049  }
2050  // from electron interface, hence using mode
2051  h_ele_PinMnPout_mode->Fill(gsfIter->trackMomentumAtVtx().R() - gsfIter->trackMomentumOut().R());
2052  h_ele_PinMnPoutVsEta_mode->Fill(gsfIter->eta(),
2053  gsfIter->trackMomentumAtVtx().R() - gsfIter->trackMomentumOut().R());
2054  h_ele_PinMnPoutVsPhi_mode->Fill(gsfIter->phi(),
2055  gsfIter->trackMomentumAtVtx().R() - gsfIter->trackMomentumOut().R());
2056  h_ele_PinMnPoutVsPt_mode->Fill(gsfIter->pt(), gsfIter->trackMomentumAtVtx().R() - gsfIter->trackMomentumOut().R());
2057  h_ele_PinMnPoutVsE_mode->Fill(gsfIter->caloEnergy(),
2058  gsfIter->trackMomentumAtVtx().R() - gsfIter->trackMomentumOut().R());
2059  if (!readAOD_) // track extra does not exist in AOD
2060  h_ele_PinMnPoutVsChi2_mode->Fill(gsfIter->gsfTrack()->normalizedChi2(),
2061  gsfIter->trackMomentumAtVtx().R() - gsfIter->trackMomentumOut().R());
2062  h_ele_outerP_mode->Fill(gsfIter->trackMomentumOut().R());
2063  h_ele_outerPVsEta_mode->Fill(gsfIter->eta(), gsfIter->trackMomentumOut().R());
2064  h_ele_outerPt_mode->Fill(gsfIter->trackMomentumOut().Rho());
2065  h_ele_outerPtVsEta_mode->Fill(gsfIter->eta(), gsfIter->trackMomentumOut().Rho());
2066  h_ele_outerPtVsPhi_mode->Fill(gsfIter->phi(), gsfIter->trackMomentumOut().Rho());
2067  h_ele_outerPtVsPt_mode->Fill(gsfIter->pt(), gsfIter->trackMomentumOut().Rho());
2068 
2069  if (!readAOD_) { // track extra does not exist in AOD
2070  edm::RefToBase<TrajectorySeed> seed = gsfIter->gsfTrack()->extra()->seedRef();
2071  ElectronSeedRef elseed = seed.castTo<ElectronSeedRef>();
2072  h_ele_seed_dphi2_->Fill(elseed->dPhiNeg(1));
2073  h_ele_seed_dphi2VsEta_->Fill(gsfIter->eta(), elseed->dPhiNeg(1));
2074  h_ele_seed_dphi2VsPt_->Fill(gsfIter->pt(), elseed->dPhiNeg(1));
2075  h_ele_seed_drz2_->Fill(elseed->dRZNeg(1));
2076  h_ele_seed_drz2VsEta_->Fill(gsfIter->eta(), elseed->dRZNeg(1));
2077  h_ele_seed_drz2VsPt_->Fill(gsfIter->pt(), elseed->dRZNeg(1));
2078  h_ele_seed_subdet2_->Fill(elseed->subDet(1));
2079  }
2080  // match distributions
2081  h_ele_EoP->Fill(gsfIter->eSuperClusterOverP());
2082  h_ele_EoPVsEta->Fill(gsfIter->eta(), gsfIter->eSuperClusterOverP());
2083  h_ele_EoPVsPhi->Fill(gsfIter->phi(), gsfIter->eSuperClusterOverP());
2084  h_ele_EoPVsE->Fill(gsfIter->caloEnergy(), gsfIter->eSuperClusterOverP());
2085  h_ele_EseedOP->Fill(gsfIter->eSeedClusterOverP());
2086  h_ele_EseedOPVsEta->Fill(gsfIter->eta(), gsfIter->eSeedClusterOverP());
2087  h_ele_EseedOPVsPhi->Fill(gsfIter->phi(), gsfIter->eSeedClusterOverP());
2088  h_ele_EseedOPVsE->Fill(gsfIter->caloEnergy(), gsfIter->eSeedClusterOverP());
2089  h_ele_EoPout->Fill(gsfIter->eSeedClusterOverPout());
2090  h_ele_EoPoutVsEta->Fill(gsfIter->eta(), gsfIter->eSeedClusterOverPout());
2091  h_ele_EoPoutVsPhi->Fill(gsfIter->phi(), gsfIter->eSeedClusterOverPout());
2092  h_ele_EoPoutVsE->Fill(gsfIter->caloEnergy(), gsfIter->eSeedClusterOverPout());
2093  h_ele_EeleOPout->Fill(gsfIter->eEleClusterOverPout());
2094  h_ele_EeleOPoutVsEta->Fill(gsfIter->eta(), gsfIter->eEleClusterOverPout());
2095  h_ele_EeleOPoutVsPhi->Fill(gsfIter->phi(), gsfIter->eEleClusterOverPout());
2096  h_ele_EeleOPoutVsE->Fill(gsfIter->caloEnergy(), gsfIter->eEleClusterOverPout());
2097  h_ele_dEtaSc_propVtx->Fill(gsfIter->deltaEtaSuperClusterTrackAtVtx());
2098  h_ele_dEtaScVsEta_propVtx->Fill(gsfIter->eta(), gsfIter->deltaEtaSuperClusterTrackAtVtx());
2099  h_ele_dEtaScVsPhi_propVtx->Fill(gsfIter->phi(), gsfIter->deltaEtaSuperClusterTrackAtVtx());
2100  h_ele_dEtaScVsPt_propVtx->Fill(gsfIter->pt(), gsfIter->deltaEtaSuperClusterTrackAtVtx());
2101  h_ele_dPhiSc_propVtx->Fill(gsfIter->deltaPhiSuperClusterTrackAtVtx());
2102  h_ele_dPhiScVsEta_propVtx->Fill(gsfIter->eta(), gsfIter->deltaPhiSuperClusterTrackAtVtx());
2103  h_ele_dPhiScVsPhi_propVtx->Fill(gsfIter->phi(), gsfIter->deltaPhiSuperClusterTrackAtVtx());
2104  h_ele_dPhiScVsPt_propVtx->Fill(gsfIter->pt(), gsfIter->deltaPhiSuperClusterTrackAtVtx());
2105  h_ele_dEtaCl_propOut->Fill(gsfIter->deltaEtaSeedClusterTrackAtCalo());
2106  h_ele_dEtaClVsEta_propOut->Fill(gsfIter->eta(), gsfIter->deltaEtaSeedClusterTrackAtCalo());
2107  h_ele_dEtaClVsPhi_propOut->Fill(gsfIter->phi(), gsfIter->deltaEtaSeedClusterTrackAtCalo());
2108  h_ele_dEtaClVsPt_propOut->Fill(gsfIter->pt(), gsfIter->deltaEtaSeedClusterTrackAtCalo());
2109  h_ele_dPhiCl_propOut->Fill(gsfIter->deltaPhiSeedClusterTrackAtCalo());
2110  h_ele_dPhiClVsEta_propOut->Fill(gsfIter->eta(), gsfIter->deltaPhiSeedClusterTrackAtCalo());
2111  h_ele_dPhiClVsPhi_propOut->Fill(gsfIter->phi(), gsfIter->deltaPhiSeedClusterTrackAtCalo());
2112  h_ele_dPhiClVsPt_propOut->Fill(gsfIter->pt(), gsfIter->deltaPhiSeedClusterTrackAtCalo());
2113  h_ele_dEtaEleCl_propOut->Fill(gsfIter->deltaEtaEleClusterTrackAtCalo());
2114  h_ele_dEtaEleClVsEta_propOut->Fill(gsfIter->eta(), gsfIter->deltaEtaEleClusterTrackAtCalo());
2115  h_ele_dEtaEleClVsPhi_propOut->Fill(gsfIter->phi(), gsfIter->deltaEtaEleClusterTrackAtCalo());
2116  h_ele_dEtaEleClVsPt_propOut->Fill(gsfIter->pt(), gsfIter->deltaEtaEleClusterTrackAtCalo());
2117  h_ele_dPhiEleCl_propOut->Fill(gsfIter->deltaPhiEleClusterTrackAtCalo());
2118  h_ele_dPhiEleClVsEta_propOut->Fill(gsfIter->eta(), gsfIter->deltaPhiEleClusterTrackAtCalo());
2119  h_ele_dPhiEleClVsPhi_propOut->Fill(gsfIter->phi(), gsfIter->deltaPhiEleClusterTrackAtCalo());
2120  h_ele_dPhiEleClVsPt_propOut->Fill(gsfIter->pt(), gsfIter->deltaPhiEleClusterTrackAtCalo());
2121  h_ele_HoE->Fill(gsfIter->hadronicOverEm());
2122  if (!gsfIter->isEBEtaGap() && !gsfIter->isEBPhiGap() && !gsfIter->isEBEEGap() && !gsfIter->isEERingGap() &&
2123  !gsfIter->isEEDeeGap())
2124  h_ele_HoE_fiducial->Fill(gsfIter->hadronicOverEm());
2125  h_ele_HoEVsEta->Fill(gsfIter->eta(), gsfIter->hadronicOverEm());
2126  h_ele_HoEVsPhi->Fill(gsfIter->phi(), gsfIter->hadronicOverEm());
2127  h_ele_HoEVsE->Fill(gsfIter->caloEnergy(), gsfIter->hadronicOverEm());
2128 
2129  //classes
2130  int eleClass = gsfIter->classification();
2131  if (gsfIter->isEE())
2132  eleClass += 10;
2133  h_ele_classes->Fill(eleClass);
2134 
2135  h_ele_eta->Fill(std::abs(gsfIter->eta()));
2136  if (gsfIter->classification() == GsfElectron::GOLDEN)
2137  h_ele_eta_golden->Fill(std::abs(gsfIter->eta()));
2138  if (gsfIter->classification() == GsfElectron::BIGBREM)
2139  h_ele_eta_bbrem->Fill(std::abs(gsfIter->eta()));
2140  //if (gsfIter->classification() == GsfElectron::OLDNARROW) h_ele_eta_narrow ->Fill(std::abs(gsfIter->eta()));
2141  if (gsfIter->classification() == GsfElectron::SHOWERING)
2142  h_ele_eta_shower->Fill(std::abs(gsfIter->eta()));
2143 
2144  //fbrem
2145  double fbrem_mean = 0.;
2146  if (!readAOD_) // track extra does not exist in AOD
2147  fbrem_mean = 1. - gsfIter->gsfTrack()->outerMomentum().R() / gsfIter->gsfTrack()->innerMomentum().R();
2148  double fbrem_mode = gsfIter->fbrem();
2149  h_ele_fbrem->Fill(fbrem_mode);
2150  h_ele_fbremVsEta_mode->Fill(gsfIter->eta(), fbrem_mode);
2151  if (!readAOD_) // track extra does not exist in AOD
2152  h_ele_fbremVsEta_mean->Fill(gsfIter->eta(), fbrem_mean);
2153 
2154  if (gsfIter->classification() == GsfElectron::GOLDEN)
2155  h_ele_PinVsPoutGolden_mode->Fill(gsfIter->trackMomentumOut().R(), gsfIter->trackMomentumAtVtx().R());
2156  if (gsfIter->classification() == GsfElectron::SHOWERING)
2157  h_ele_PinVsPoutShowering_mode->Fill(gsfIter->trackMomentumOut().R(), gsfIter->trackMomentumAtVtx().R());
2158  if (!readAOD_) // track extra does not exist in AOD
2159  if (gsfIter->classification() == GsfElectron::GOLDEN)
2160  h_ele_PinVsPoutGolden_mean->Fill(gsfIter->gsfTrack()->outerMomentum().R(),
2161  gsfIter->gsfTrack()->innerMomentum().R());
2162  if (!readAOD_) // track extra does not exist in AOD
2163  if (gsfIter->classification() == GsfElectron::SHOWERING)
2164  h_ele_PinVsPoutShowering_mean->Fill(gsfIter->gsfTrack()->outerMomentum().R(),
2165  gsfIter->gsfTrack()->innerMomentum().R());
2166  if (gsfIter->classification() == GsfElectron::GOLDEN)
2167  h_ele_PtinVsPtoutGolden_mode->Fill(gsfIter->trackMomentumOut().Rho(), gsfIter->trackMomentumAtVtx().Rho());
2168  if (gsfIter->classification() == GsfElectron::SHOWERING)
2169  h_ele_PtinVsPtoutShowering_mode->Fill(gsfIter->trackMomentumOut().Rho(), gsfIter->trackMomentumAtVtx().Rho());
2170  if (!readAOD_) // track extra does not exist in AOD
2171  if (gsfIter->classification() == GsfElectron::GOLDEN)
2172  h_ele_PtinVsPtoutGolden_mean->Fill(gsfIter->gsfTrack()->outerMomentum().Rho(),
2173  gsfIter->gsfTrack()->innerMomentum().Rho());
2174  if (!readAOD_) // track extra does not exist in AOD
2175  if (gsfIter->classification() == GsfElectron::SHOWERING)
2176  h_ele_PtinVsPtoutShowering_mean->Fill(gsfIter->gsfTrack()->outerMomentum().Rho(),
2177  gsfIter->gsfTrack()->innerMomentum().Rho());
2178 
2179  h_ele_mva->Fill(gsfIter->mva_e_pi());
2180  if (gsfIter->ecalDrivenSeed())
2181  h_ele_provenance->Fill(1.);
2182  if (gsfIter->trackerDrivenSeed())
2183  h_ele_provenance->Fill(-1.);
2184  if (gsfIter->trackerDrivenSeed() || gsfIter->ecalDrivenSeed())
2185  h_ele_provenance->Fill(0.);
2186  if (gsfIter->trackerDrivenSeed() && !gsfIter->ecalDrivenSeed())
2187  h_ele_provenance->Fill(-2.);
2188  if (!gsfIter->trackerDrivenSeed() && gsfIter->ecalDrivenSeed())
2189  h_ele_provenance->Fill(2.);
2190 
2191  h_ele_tkSumPt_dr03->Fill(gsfIter->dr03TkSumPt());
2192  h_ele_ecalRecHitSumEt_dr03->Fill(gsfIter->dr03EcalRecHitSumEt());
2193  h_ele_hcalDepth1TowerSumEt_dr03->Fill(gsfIter->dr03HcalDepth1TowerSumEt());
2194  h_ele_hcalDepth2TowerSumEt_dr03->Fill(gsfIter->dr03HcalDepth2TowerSumEt());
2195  h_ele_tkSumPt_dr04->Fill(gsfIter->dr04TkSumPt());
2196  h_ele_ecalRecHitSumEt_dr04->Fill(gsfIter->dr04EcalRecHitSumEt());
2197  h_ele_hcalDepth1TowerSumEt_dr04->Fill(gsfIter->dr04HcalDepth1TowerSumEt());
2198  h_ele_hcalDepth2TowerSumEt_dr04->Fill(gsfIter->dr04HcalDepth2TowerSumEt());
2199 
2200  float enrj1 = gsfIter->superCluster()->energy();
2201  // mee
2202  for (reco::GsfElectronCollection::const_iterator gsfIter2 = gsfIter + 1; gsfIter2 != gsfElectrons->end();
2203  gsfIter2++) {
2204  math::XYZTLorentzVector p12 = (*gsfIter).p4() + (*gsfIter2).p4();
2205  float mee2 = p12.Dot(p12);
2206  float enrj2 = gsfIter2->superCluster()->energy();
2207  h_ele_mee_all->Fill(sqrt(mee2));
2208  h_ele_E2mnE1vsMee_all->Fill(sqrt(mee2), enrj2 - enrj1);
2209  if (gsfIter->ecalDrivenSeed() && gsfIter2->ecalDrivenSeed())
2210  h_ele_E2mnE1vsMee_egeg_all->Fill(sqrt(mee2), enrj2 - enrj1);
2211  if (gsfIter->charge() * gsfIter2->charge() < 0.) {
2212  h_ele_mee_os->Fill(sqrt(mee2));
2213  if (gsfIter->isEB() && gsfIter2->isEB())
2214  h_ele_mee_os_ebeb->Fill(sqrt(mee2));
2215  if ((gsfIter->isEB() && gsfIter2->isEE()) || (gsfIter->isEE() && gsfIter2->isEB()))
2216  h_ele_mee_os_ebee->Fill(sqrt(mee2));
2217  if (gsfIter->isEE() && gsfIter2->isEE())
2218  h_ele_mee_os_eeee->Fill(sqrt(mee2));
2219  if ((gsfIter->classification()==GsfElectron::GOLDEN && gsfIter2->classification()==GsfElectron::GOLDEN) ||
2220  (gsfIter->classification()==GsfElectron::GOLDEN && gsfIter2->classification()==GsfElectron::BIGBREM) ||
2221  //(gsfIter->classification()==GsfElectron::GOLDEN && gsfIter2->classification()==GsfElectron::OLDNARROW) ||
2222  (gsfIter->classification()==GsfElectron::BIGBREM && gsfIter2->classification()==GsfElectron::GOLDEN) ||
2223  (gsfIter->classification()==GsfElectron::BIGBREM && gsfIter2->classification()==GsfElectron::BIGBREM)/* ||
2224  (gsfIter->classification()==GsfElectron::BIGBREM && gsfIter2->classification()==GsfElectron::OLDNARROW) ||
2225  (gsfIter->classification()==GsfElectron::OLDNARROW && gsfIter2->classification()==GsfElectron::GOLDEN) ||
2226  (gsfIter->classification()==GsfElectron::OLDNARROW && gsfIter2->classification()==GsfElectron::BIGBREM) ||
2227  (gsfIter->classification()==GsfElectron::OLDNARROW && gsfIter2->classification()==GsfElectron::OLDNARROW)*/)
2228  {
2229  h_ele_mee_os_gg->Fill(sqrt(mee2));
2230  } else if ((gsfIter->classification() == GsfElectron::SHOWERING &&
2231  gsfIter2->classification() == GsfElectron::SHOWERING) ||
2232  (gsfIter->classification() == GsfElectron::SHOWERING && gsfIter2->isGap()) ||
2233  (gsfIter->isGap() && gsfIter2->classification() == GsfElectron::SHOWERING) ||
2234  (gsfIter->isGap() && gsfIter2->isGap())) {
2235  h_ele_mee_os_bb->Fill(sqrt(mee2));
2236  } else {
2237  h_ele_mee_os_gb->Fill(sqrt(mee2));
2238  }
2239  }
2240  }
2241  }
2242 
2243  // association matching object-reco electrons
2244  int matchingObjectNum = 0;
2245 
2246  for (reco::SuperClusterCollection::const_iterator moIter = recoClusters->begin(); moIter != recoClusters->end();
2247  moIter++) {
2248  // number of matching objects
2249  matchingObjectNum++;
2250 
2251  if (moIter->energy() / cosh(moIter->eta()) > maxPtMatchingObject_ ||
2252  std::abs(moIter->eta()) > maxAbsEtaMatchingObject_)
2253  continue;
2254 
2255  // suppress the endcaps
2256  //if (std::abs(moIter->eta()) > 1.5) continue;
2257  // select central z
2258  //if ( std::abs((*mcIter)->production_vertex()->position().z())>50.) continue;
2259 
2260  h_matchingObjectEta->Fill(moIter->eta());
2261  h_matchingObjectAbsEta->Fill(std::abs(moIter->eta()));
2262  h_matchingObjectP->Fill(moIter->energy());
2263  h_matchingObjectPt->Fill(moIter->energy() / cosh(moIter->eta()));
2264  h_matchingObjectPhi->Fill(moIter->phi());
2265  h_matchingObjectZ->Fill(moIter->z());
2266 
2267  // looking for the best matching gsf electron
2268  bool okGsfFound = false;
2269  double gsfOkRatio = 999999.;
2270 
2271  // find matching electron
2272  reco::GsfElectron bestGsfElectron;
2273  for (reco::GsfElectronCollection::const_iterator gsfIter = gsfElectrons->begin(); gsfIter != gsfElectrons->end();
2274  gsfIter++) {
2275  // matching with a cone in eta phi
2276  if (matchingCondition_ == "Cone") {
2277  double dphi = gsfIter->phi() - moIter->phi();
2278  if (std::abs(dphi) > CLHEP::pi)
2279  dphi = dphi < 0 ? (CLHEP::twopi) + dphi : dphi - CLHEP::twopi;
2280  double deltaR = sqrt(std::pow((moIter->eta() - gsfIter->eta()), 2) + std::pow(dphi, 2));
2281  if (deltaR < deltaR_) {
2282  //if ( (genPc->pdg_id() == 11) && (gsfIter->charge() < 0.) || (genPc->pdg_id() == -11) &&
2283  //(gsfIter->charge() > 0.) ){
2284  double tmpGsfRatio = gsfIter->p() / moIter->energy();
2285  if (std::abs(tmpGsfRatio - 1) < std::abs(gsfOkRatio - 1)) {
2286  gsfOkRatio = tmpGsfRatio;
2287  bestGsfElectron = *gsfIter;
2288  okGsfFound = true;
2289  }
2290  //}
2291  }
2292  }
2293  } // loop over rec ele to look for the best one
2294 
2295  // analysis when the matching object is matched by a rec electron
2296  if (okGsfFound) {
2297  // generated distributions for matched electrons
2298  h_ele_matchingObjectPt_matched->Fill(moIter->energy() / cosh(moIter->eta()));
2299  h_ele_matchingObjectPhi_matched->Fill(moIter->phi());
2300  h_ele_matchingObjectAbsEta_matched->Fill(std::abs(moIter->eta()));
2301  h_ele_matchingObjectEta_matched->Fill(moIter->eta());
2302  h_ele_matchingObjectZ_matched->Fill(moIter->z());
2303 
2304  // comparison electron vs matching object
2305  h_ele_EtaMnEtamatchingObject_matched->Fill(bestGsfElectron.eta() - moIter->eta());
2306  h_ele_EtaMnEtamatchingObjectVsEta_matched->Fill(bestGsfElectron.eta(), bestGsfElectron.eta() - moIter->eta());
2307  h_ele_EtaMnEtamatchingObjectVsPhi_matched->Fill(bestGsfElectron.phi(), bestGsfElectron.eta() - moIter->eta());
2308  h_ele_EtaMnEtamatchingObjectVsPt_matched->Fill(bestGsfElectron.pt(), bestGsfElectron.eta() - moIter->eta());
2309  h_ele_PhiMnPhimatchingObject_matched->Fill(bestGsfElectron.phi() - moIter->phi());
2310  h_ele_PhiMnPhimatchingObject2_matched->Fill(bestGsfElectron.phi() - moIter->phi());
2311  h_ele_PhiMnPhimatchingObjectVsEta_matched->Fill(bestGsfElectron.eta(), bestGsfElectron.phi() - moIter->phi());
2312  h_ele_PhiMnPhimatchingObjectVsPhi_matched->Fill(bestGsfElectron.phi(), bestGsfElectron.phi() - moIter->phi());
2313  h_ele_PhiMnPhimatchingObjectVsPt_matched->Fill(bestGsfElectron.pt(), bestGsfElectron.phi() - moIter->phi());
2314  h_ele_PoPmatchingObject_matched->Fill(bestGsfElectron.p() / moIter->energy());
2315  h_ele_PtoPtmatchingObject_matched->Fill(bestGsfElectron.pt() / moIter->energy() / cosh(moIter->eta()));
2316  h_ele_PoPmatchingObjectVsEta_matched->Fill(bestGsfElectron.eta(), bestGsfElectron.p() / moIter->energy());
2317  h_ele_PoPmatchingObjectVsPhi_matched->Fill(bestGsfElectron.phi(), bestGsfElectron.p() / moIter->energy());
2318  h_ele_PoPmatchingObjectVsPt_matched->Fill(bestGsfElectron.py(), bestGsfElectron.p() / moIter->energy());
2319  if (bestGsfElectron.isEB())
2320  h_ele_PoPmatchingObject_barrel_matched->Fill(bestGsfElectron.p() / moIter->energy());
2321  if (bestGsfElectron.isEE())
2322  h_ele_PoPmatchingObject_endcaps_matched->Fill(bestGsfElectron.p() / moIter->energy());
2323  if (bestGsfElectron.isEB())
2324  h_ele_PtoPtmatchingObject_barrel_matched->Fill(bestGsfElectron.pt() / moIter->energy() / cosh(moIter->eta()));
2325  if (bestGsfElectron.isEE())
2326  h_ele_PtoPtmatchingObject_endcaps_matched->Fill(bestGsfElectron.pt() / moIter->energy() / cosh(moIter->eta()));
2327 
2328  reco::SuperClusterRef sclRef = bestGsfElectron.superCluster();
2329  if (bestGsfElectron.isEB())
2330  histSclEoEmatchingObject_barrel_matched->Fill(sclRef->energy() / moIter->energy());
2331  if (bestGsfElectron.isEE())
2332  histSclEoEmatchingObject_endcaps_matched->Fill(sclRef->energy() / moIter->energy());
2333  if (bestGsfElectron.isEB())
2334  histSclEoEmatchingObject_barrel_new_matched->Fill(sclRef->energy() / moIter->energy());
2335  if (bestGsfElectron.isEE())
2336  histSclEoEmatchingObject_endcaps_new_matched->Fill(sclRef->energy() / moIter->energy());
2337 
2338  // add here distributions for matched electrons as for all electrons
2339  //..
2340 
2341  } // gsf electron found
2342 
2343  } // loop overmatching object
2344 
2345  h_matchingObjectNum->Fill(matchingObjectNum);
2346 }
2347 
2348 bool GsfElectronDataAnalyzer::trigger(const edm::Event &e) {
2349  // retreive TriggerResults from the event
2350  edm::Handle<edm::TriggerResults> triggerResults;
2351  e.getByLabel(triggerResults_, triggerResults);
2352 
2353  bool accept = false;
2354 
2355  if (triggerResults.isValid()) {
2356  //std::cout << "TriggerResults found, number of HLT paths: " << triggerResults->size() << std::endl;
2357 
2358  // get trigger names
2359  const edm::TriggerNames &triggerNames = e.triggerNames(*triggerResults);
2360  if (nEvents_ == 1) {
2361  for (unsigned int i = 0; i < triggerNames.size(); i++) {
2362  std::cout << "trigger path= " << triggerNames.triggerName(i) << std::endl;
2363  }
2364  }
2365 
2366  unsigned int n = HLTPathsByName_.size();
2367  for (unsigned int i = 0; i != n; i++) {
2368  HLTPathsByIndex_[i] = triggerNames.triggerIndex(HLTPathsByName_[i]);
2369  }
2370 
2371  // empty input vectors (n==0) means any trigger paths
2372  if (n == 0) {
2373  n = triggerResults->size();
2374  HLTPathsByName_.resize(n);
2375  HLTPathsByIndex_.resize(n);
2376  for (unsigned int i = 0; i != n; i++) {
2377  HLTPathsByName_[i] = triggerNames.triggerName(i);
2378  HLTPathsByIndex_[i] = i;
2379  }
2380  }
2381 
2382  if (nEvents_ == 1) {
2383  if (n > 0) {
2384  std::cout << "HLT trigger paths requested: index, name and valididty:" << std::endl;
2385  for (unsigned int i = 0; i != n; i++) {
2386  bool validity = HLTPathsByIndex_[i] < triggerResults->size();
2387  std::cout << " " << HLTPathsByIndex_[i] << " " << HLTPathsByName_[i] << " " << validity << std::endl;
2388  }
2389  }
2390  }
2391 
2392  // count number of requested HLT paths which have fired
2393  unsigned int fired = 0;
2394  for (unsigned int i = 0; i != n; i++) {
2395  if (HLTPathsByIndex_[i] < triggerResults->size()) {
2396  if (triggerResults->accept(HLTPathsByIndex_[i])) {
2397  fired++;
2398  std::cout << "Fired HLT path= " << HLTPathsByName_[i] << std::endl;
2399  accept = true;
2400  }
2401  }
2402  }
2403  }
2404 
2405  return accept;
2406 }
GsfElectronDataAnalyzer::h_ele_HoEVsPhi
TH2F * h_ele_HoEVsPhi
Definition: GsfElectronDataAnalyzer.cc:313
GsfElectronDataAnalyzer::h_ele_E2mnE1vsMee_egeg_all
TH2F * h_ele_E2mnE1vsMee_egeg_all
Definition: GsfElectronDataAnalyzer.cc:121
GsfElectronDataAnalyzer::h_ele_PoPmatchingObject_endcaps_matched
TH1F * h_ele_PoPmatchingObject_endcaps_matched
Definition: GsfElectronDataAnalyzer.cc:203
GsfElectronDataAnalyzer::h_ele_EeleOPoutVsEta
TH2F * h_ele_EeleOPoutVsEta
Definition: GsfElectronDataAnalyzer.cc:252
GsfElectronDataAnalyzer::h_matchingObjectAbsEta
TH1F * h_matchingObjectAbsEta
Definition: GsfElectronDataAnalyzer.cc:99
GsfElectronDataAnalyzer::h_ele_PinMnPoutVsChi2_mode
TH2F * h_ele_PinMnPoutVsChi2_mode
Definition: GsfElectronDataAnalyzer.cc:221
GsfElectronDataAnalyzer::h_ele_seed_drz2VsEta_
TH2F * h_ele_seed_drz2VsEta_
Definition: GsfElectronDataAnalyzer.cc:297
reco::GsfElectron::isEE
bool isEE() const
Definition: GsfElectron.h:336
GsfElectronDataAnalyzer::h_ele_EtaMnEtamatchingObject_matched
TH1F * h_ele_EtaMnEtamatchingObject_matched
Definition: GsfElectronDataAnalyzer.cc:206
Handle.h
GsfElectronDataAnalyzer::h_ele_PinMnPout
TH1F * h_ele_PinMnPout
Definition: GsfElectronDataAnalyzer.cc:215
GsfElectronDataAnalyzer::outputFile_
std::string outputFile_
Definition: GsfElectronDataAnalyzer.cc:75
GsfElectronDataAnalyzer::h_ele_fbrem
TH1F * h_ele_fbrem
Definition: GsfElectronDataAnalyzer.cc:316
EDAnalyzer.h
mps_fire.i
i
Definition: mps_fire.py:428
GsfElectronDataAnalyzer::h_ele_EeleOPout_barrel
TH1F * h_ele_EeleOPout_barrel
Definition: GsfElectronDataAnalyzer.cc:250
GsfElectronDataAnalyzer::histSclE5x5_endcaps_
TH1F * histSclE5x5_endcaps_
Definition: GsfElectronDataAnalyzer.cc:172
reco::GsfElectron::isEB
bool isEB() const
Definition: GsfElectron.h:335
GsfElectronDataAnalyzer::etamax
double etamax
Definition: GsfElectronDataAnalyzer.cc:384
GsfElectronDataAnalyzer::h_ele_chi2VsPhi
TH2F * h_ele_chi2VsPhi
Definition: GsfElectronDataAnalyzer.cc:194
GsfElectronDataAnalyzer::readAOD_
bool readAOD_
Definition: GsfElectronDataAnalyzer.cc:80
GsfElectronDataAnalyzer::h_matchingObjectP
TH1F * h_matchingObjectP
Definition: GsfElectronDataAnalyzer.cc:100
GsfElectronDataAnalyzer::h_ele_chargeVsPt
TH2F * h_ele_chargeVsPt
Definition: GsfElectronDataAnalyzer.cc:126
GsfElectronDataAnalyzer::seedEta
float seedEta[10]
Definition: GsfElectronDataAnalyzer.cc:94
GsfElectronDataAnalyzer::h_ele_foundHits
TH1F * h_ele_foundHits
Definition: GsfElectronDataAnalyzer.cc:178
MessageLogger.h
GsfElectronDataAnalyzer::nbindphimatch2D
int nbindphimatch2D
Definition: GsfElectronDataAnalyzer.cc:426
GsfElectronDataAnalyzer::h_ele_vertexPt
TH1F * h_ele_vertexPt
Definition: GsfElectronDataAnalyzer.cc:128
dqmiodumpmetadata.n
n
Definition: dqmiodumpmetadata.py:28
GsfElectronDataAnalyzer::tipMaxBarrel_
double tipMaxBarrel_
Definition: GsfElectronDataAnalyzer.cc:373
TrackerGeometry.h
GsfElectronDataAnalyzer::h_ele_ambiguousTracksVsEta
TH2F * h_ele_ambiguousTracksVsEta
Definition: GsfElectronDataAnalyzer.cc:175
GsfElectronDataAnalyzer::h_ele_PinMnPoutVsEta_mode
TH2F * h_ele_PinMnPoutVsEta_mode
Definition: GsfElectronDataAnalyzer.cc:217
GsfElectronDataAnalyzer::h_ele_dEtaSc_propVtx_barrel
TH1F * h_ele_dEtaSc_propVtx_barrel
Definition: GsfElectronDataAnalyzer.cc:257
GsfElectronDataAnalyzer::h_ele_PhiMnPhimatchingObject_matched
TH1F * h_ele_PhiMnPhimatchingObject_matched
Definition: GsfElectronDataAnalyzer.cc:210
ESHandle.h
TriggerResults.h
GsfElectronDataAnalyzer::h_ele_hcalDepth2TowerSumEt_dr03
TH1F * h_ele_hcalDepth2TowerSumEt_dr03
Definition: GsfElectronDataAnalyzer.cc:339
GsfElectronDataAnalyzer::h_ele_lostHitsVsPhi
TH2F * h_ele_lostHitsVsPhi
Definition: GsfElectronDataAnalyzer.cc:188
GsfElectronDataAnalyzer::phimax
double phimax
Definition: GsfElectronDataAnalyzer.cc:386
GsfElectronDataAnalyzer::h_ele_PhiMnPhimatchingObjectVsPt_matched
TH2F * h_ele_PhiMnPhimatchingObjectVsPt_matched
Definition: GsfElectronDataAnalyzer.cc:214
GsfElectronDataAnalyzer::h_ele_foundHits_endcaps
TH1F * h_ele_foundHits_endcaps
Definition: GsfElectronDataAnalyzer.cc:180
GsfElectronDataAnalyzer::detamax
double detamax
Definition: GsfElectronDataAnalyzer.cc:392
GsfElectronDataAnalyzer::poptruemin
double poptruemin
Definition: GsfElectronDataAnalyzer.cc:401
GsfElectronDataAnalyzer::h_ele_seed_subdet2_
TH1F * h_ele_seed_subdet2_
Definition: GsfElectronDataAnalyzer.cc:299
GsfElectronDataAnalyzer::h_ele_PoPmatchingObjectVsPhi_matched
TH2F * h_ele_PoPmatchingObjectVsPhi_matched
Definition: GsfElectronDataAnalyzer.cc:200
GsfElectronDataAnalyzer::histSclE5x5_
TH1F * histSclE5x5_
Definition: GsfElectronDataAnalyzer.cc:170
GsfElectronDataAnalyzer::h_ele_Et
TH1F * h_ele_Et
Definition: GsfElectronDataAnalyzer.cc:129
GsfElectronDataAnalyzer::h_ele_matchingObjectPt_matched
TH1F * h_ele_matchingObjectPt_matched
Definition: GsfElectronDataAnalyzer.cc:107
GsfElectronDataAnalyzer::nbinfhits
int nbinfhits
Definition: GsfElectronDataAnalyzer.cc:418
GsfElectronDataAnalyzer::eOverPMaxEndcaps_
double eOverPMaxEndcaps_
Definition: GsfElectronDataAnalyzer.cc:357
GsfElectronDataAnalyzer::eOverPMaxBarrel_
double eOverPMaxBarrel_
Definition: GsfElectronDataAnalyzer.cc:355
edm
HLT enums.
Definition: AlignableModifier.h:19
GsfElectronDataAnalyzer::nbinmee
int nbinmee
Definition: GsfElectronDataAnalyzer.cc:428
GsfElectronDataAnalyzer::h_ele_matchingObjectEta_matched
TH1F * h_ele_matchingObjectEta_matched
Definition: GsfElectronDataAnalyzer.cc:105
GsfElectronDataAnalyzer::h_ele_fbremVsEta_mode
TProfile * h_ele_fbremVsEta_mode
Definition: GsfElectronDataAnalyzer.cc:317
GsfElectronDataAnalyzer::meemax
double meemax
Definition: GsfElectronDataAnalyzer.cc:404
GsfElectronDataAnalyzer::h_ele_HoEVsE
TH2F * h_ele_HoEVsE
Definition: GsfElectronDataAnalyzer.cc:314
GsfElectronDataAnalyzer::h_ele_dEtaEleCl_propOut_barrel
TH1F * h_ele_dEtaEleCl_propOut_barrel
Definition: GsfElectronDataAnalyzer.cc:281
GsfElectronDataAnalyzer::h_ele_HoEVsEta
TH2F * h_ele_HoEVsEta
Definition: GsfElectronDataAnalyzer.cc:312
GsfElectronDataAnalyzer::h_ele_dPhiEleClVsEta_propOut
TH2F * h_ele_dPhiEleClVsEta_propOut
Definition: GsfElectronDataAnalyzer.cc:289
gather_cfg.cout
cout
Definition: gather_cfg.py:144
GsfElectronDataAnalyzer::eOverPMinBarrel_
double eOverPMinBarrel_
Definition: GsfElectronDataAnalyzer.cc:354
GsfElectronDataAnalyzer::maxAbsEtaMatchingObject_
double maxAbsEtaMatchingObject_
Definition: GsfElectronDataAnalyzer.cc:83
GsfElectronDataAnalyzer::h_ele_dPhiEleCl_propOut_barrel
TH1F * h_ele_dPhiEleCl_propOut_barrel
Definition: GsfElectronDataAnalyzer.cc:287
HLT_FULL_cff.InputTag
InputTag
Definition: HLT_FULL_cff.py:89285
GsfElectronDataAnalyzer::matchingObjectCollection_
edm::InputTag matchingObjectCollection_
Definition: GsfElectronDataAnalyzer.cc:77
GsfElectronDataAnalyzer::h_ele_EseedOPVsEta
TH2F * h_ele_EseedOPVsEta
Definition: GsfElectronDataAnalyzer.cc:240
triggerResults
static const std::string triggerResults
Definition: EdmProvDump.cc:45
GsfElectronDataAnalyzer::h_matchingObjectNum
TH1F * h_matchingObjectNum
Definition: GsfElectronDataAnalyzer.cc:96
GsfElectronDataAnalyzer
Definition: GsfElectronDataAnalyzer.cc:53
GsfElectronDataAnalyzer::nbineop
int nbineop
Definition: GsfElectronDataAnalyzer.cc:416
GsfElectronDataAnalyzer::h_ele_PtoPtmatchingObject_endcaps_matched
TH1F * h_ele_PtoPtmatchingObject_endcaps_matched
Definition: GsfElectronDataAnalyzer.cc:205
GsfElectronDataAnalyzer::h_ele_EoPVsE
TH2F * h_ele_EoPVsE
Definition: GsfElectronDataAnalyzer.cc:236
GsfElectronDataAnalyzer::h_ele_mee_os_ebeb
TH1F * h_ele_mee_os_ebeb
Definition: GsfElectronDataAnalyzer.cc:113
GsfElectronDataAnalyzer::h_ele_vertexPhi
TH1F * h_ele_vertexPhi
Definition: GsfElectronDataAnalyzer.cc:136
cms::cuda::assert
assert(be >=bs)
GsfElectronDataAnalyzer::h_ele_dEtaClVsPt_propOut
TH2F * h_ele_dEtaClVsPt_propOut
Definition: GsfElectronDataAnalyzer.cc:273
GsfElectronDataAnalyzer::h_ele_dPhiSc_propVtx
TH1F * h_ele_dPhiSc_propVtx
Definition: GsfElectronDataAnalyzer.cc:262
GsfElectronDataAnalyzer::nbineta
int nbineta
Definition: GsfElectronDataAnalyzer.cc:407
GsfElectronDataAnalyzer::beginJob
void beginJob() override
Definition: GsfElectronDataAnalyzer.cc:550
GsfElectronDataAnalyzer::h_ele_vertexTIPVsEta
TH2F * h_ele_vertexTIPVsEta
Definition: GsfElectronDataAnalyzer.cc:141
GsfElectronDataAnalyzer::h_ele_PinVsPoutGolden_mode
TH2F * h_ele_PinVsPoutGolden_mode
Definition: GsfElectronDataAnalyzer.cc:320
GsfElectronDataAnalyzer::h_ele_vertexAbsEta
TH1F * h_ele_vertexAbsEta
Definition: GsfElectronDataAnalyzer.cc:135
GsfElectronDataAnalyzer::h_ele_ambiguousTracksVsPt
TH2F * h_ele_ambiguousTracksVsPt
Definition: GsfElectronDataAnalyzer.cc:177
GsfElectronDataAnalyzer::histSclSigIEtaIEta_barrel_
TH1F * histSclSigIEtaIEta_barrel_
Definition: GsfElectronDataAnalyzer.cc:162
GsfElectronDataAnalyzer::h_ele_dEtaCl_propOut
TH1F * h_ele_dEtaCl_propOut
Definition: GsfElectronDataAnalyzer.cc:268
GsfElectronDataAnalyzer::nbindeta
int nbindeta
Definition: GsfElectronDataAnalyzer.cc:421
EDAnalyzer.h
GsfElectronDataAnalyzer::h_matchingObjectEta
TH1F * h_matchingObjectEta
Definition: GsfElectronDataAnalyzer.cc:98
GsfElectronDataAnalyzer::h_ele_PoPmatchingObject_matched
TH1F * h_ele_PoPmatchingObject_matched
Definition: GsfElectronDataAnalyzer.cc:198
GsfElectronDataAnalyzer::dEtaMaxEndcaps_
double dEtaMaxEndcaps_
Definition: GsfElectronDataAnalyzer.cc:361
reco::LeafCandidate::pt
double pt() const final
transverse momentum
Definition: LeafCandidate.h:146
reco
fixed size matrix
Definition: AlignmentAlgorithmBase.h:45
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
edm::one::EDAnalyzer
Definition: EDAnalyzer.h:30
GsfElectronDataAnalyzer::mcEta
float mcEta[10]
Definition: GsfElectronDataAnalyzer.cc:92
GsfElectronDataAnalyzer::h_ele_E2mnE1vsMee_all
TH2F * h_ele_E2mnE1vsMee_all
Definition: GsfElectronDataAnalyzer.cc:120
electronIsolatorFromEffectiveArea_cfi.gsfElectrons
gsfElectrons
Definition: electronIsolatorFromEffectiveArea_cfi.py:4
edm::Handle< reco::GsfElectronCollection >
GsfElectronDataAnalyzer::dphimax
double dphimax
Definition: GsfElectronDataAnalyzer.cc:394
GsfElectronDataAnalyzer::histSclSigEtaEta_barrel_
TH1F * histSclSigEtaEta_barrel_
Definition: GsfElectronDataAnalyzer.cc:159
GsfElectronDataAnalyzer::h_ele_EeleOPout_endcaps
TH1F * h_ele_EeleOPout_endcaps
Definition: GsfElectronDataAnalyzer.cc:251
GsfElectronDataAnalyzer::superclusterEta
float superclusterEta[10]
Definition: GsfElectronDataAnalyzer.cc:93
GsfElectronDataAnalyzer::histSclEoEmatchingObjectShowering_endcaps
TH1F * histSclEoEmatchingObjectShowering_endcaps
Definition: GsfElectronDataAnalyzer.cc:331
GsfElectronDataAnalyzer::h_ele_lostHits_endcaps
TH1F * h_ele_lostHits_endcaps
Definition: GsfElectronDataAnalyzer.cc:186
GsfElectronDataAnalyzer::ecalIso03MaxBarrel_
double ecalIso03MaxBarrel_
Definition: GsfElectronDataAnalyzer.cc:379
GsfElectronDataAnalyzer::h_ele_dPhiScVsPhi_propVtx
TH2F * h_ele_dPhiScVsPhi_propVtx
Definition: GsfElectronDataAnalyzer.cc:266
GsfElectronDataAnalyzer::h_ele_EoPoutVsE
TH2F * h_ele_EoPoutVsE
Definition: GsfElectronDataAnalyzer.cc:248
ElectronSeedFwd.h
edm::Ref< SuperClusterCollection >
GsfElectronDataAnalyzer::h_ele_EtaMnEtamatchingObjectVsEta_matched
TH2F * h_ele_EtaMnEtamatchingObjectVsEta_matched
Definition: GsfElectronDataAnalyzer.cc:207
GsfElectronDataAnalyzer::sigIetaIetaMinBarrel_
double sigIetaIetaMinBarrel_
Definition: GsfElectronDataAnalyzer.cc:366
GsfElectronDataAnalyzer::h_ele_mee_os_gb
TH1F * h_ele_mee_os_gb
Definition: GsfElectronDataAnalyzer.cc:117
GsfElectronDataAnalyzer::GsfElectronDataAnalyzer
GsfElectronDataAnalyzer(const edm::ParameterSet &conf)
Definition: GsfElectronDataAnalyzer.cc:437
GsfElectronDataAnalyzer::minPt_
double minPt_
Definition: GsfElectronDataAnalyzer.cc:347
GsfElectronDataAnalyzer::h_ele_chi2_endcaps_
TH1F * h_ele_chi2_endcaps_
Definition: GsfElectronDataAnalyzer.cc:192
GsfElectronDataAnalyzer::electronCollection_
edm::InputTag electronCollection_
Definition: GsfElectronDataAnalyzer.cc:76
GsfElectronDataAnalyzer::h_ele_hcalDepth2TowerSumEt_dr04
TH1F * h_ele_hcalDepth2TowerSumEt_dr04
Definition: GsfElectronDataAnalyzer.cc:343
GsfElectronDataAnalyzer::h_ele_eta_bbrem
TH1F * h_ele_eta_bbrem
Definition: GsfElectronDataAnalyzer.cc:304
GsfElectronDataAnalyzer::deltaR_
double deltaR_
Definition: GsfElectronDataAnalyzer.cc:84
fileCollector.seed
seed
Definition: fileCollector.py:127
GsfElectronDataAnalyzer::h_ele_classes
TH1F * h_ele_classes
Definition: GsfElectronDataAnalyzer.cc:301
GsfElectronDataAnalyzer::h_ele_EoP
TH1F * h_ele_EoP
Definition: GsfElectronDataAnalyzer.cc:231
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:36
GsfElectronDataAnalyzer::h_ele_PtinVsPtoutShowering_mode
TH2F * h_ele_PtinVsPtoutShowering_mode
Definition: GsfElectronDataAnalyzer.cc:325
accept
bool accept(const edm::Event &event, const edm::TriggerResults &triggerTable, const std::string &triggerPath)
Definition: TopDQMHelpers.h:31
GsfElectronDataAnalyzer::eOverPMinEndcaps_
double eOverPMinEndcaps_
Definition: GsfElectronDataAnalyzer.cc:356
GsfElectronDataAnalyzer::tree_
TTree * tree_
Definition: GsfElectronDataAnalyzer.cc:91
GsfElectronDataAnalyzer::h_ele_vertexZ
TH1F * h_ele_vertexZ
Definition: GsfElectronDataAnalyzer.cc:139
MakerMacros.h
GsfElectronDataAnalyzer::dphimin
double dphimin
Definition: GsfElectronDataAnalyzer.cc:393
GsfElectronDataAnalyzer::h_ele_dEtaEleClVsPhi_propOut
TH2F * h_ele_dEtaEleClVsPhi_propOut
Definition: GsfElectronDataAnalyzer.cc:284
GsfElectronDataAnalyzer::HLTPathsByIndex_
std::vector< unsigned int > HLTPathsByIndex_
Definition: GsfElectronDataAnalyzer.cc:73
cms::cuda::bs
bs
Definition: HistoContainer.h:127
GsfElectronDataAnalyzer::h_ele_seed_drz2_
TH1F * h_ele_seed_drz2_
Definition: GsfElectronDataAnalyzer.cc:296
GsfElectronDataAnalyzer::theMagField
edm::ESHandle< MagneticField > theMagField
Definition: GsfElectronDataAnalyzer.cc:88
GsfElectronDataAnalyzer::h_ele_vertexTIPVsPhi
TH2F * h_ele_vertexTIPVsPhi
Definition: GsfElectronDataAnalyzer.cc:142
GsfElectronDataAnalyzer::histSclEtVsEta_
TH2F * histSclEtVsEta_
Definition: GsfElectronDataAnalyzer.cc:153
GsfElectronDataAnalyzer::h_ele_chargeVsEta
TH2F * h_ele_chargeVsEta
Definition: GsfElectronDataAnalyzer.cc:124
GsfElectronDataAnalyzer::h_ele_dPhiSc_propVtx_endcaps
TH1F * h_ele_dPhiSc_propVtx_endcaps
Definition: GsfElectronDataAnalyzer.cc:264
GsfElectronDataAnalyzer::h_ele_dPhiScVsEta_propVtx
TH2F * h_ele_dPhiScVsEta_propVtx
Definition: GsfElectronDataAnalyzer.cc:265
GsfElectronDataAnalyzer::nbineop2D
int nbineop2D
Definition: GsfElectronDataAnalyzer.cc:417
trigger::size_type
uint16_t size_type
Definition: TriggerTypeDefs.h:18
GsfElectronDataAnalyzer::histSclE5x5_barrel_
TH1F * histSclE5x5_barrel_
Definition: GsfElectronDataAnalyzer.cc:171
GsfElectronDataAnalyzer::dPhiMaxEndcaps_
double dPhiMaxEndcaps_
Definition: GsfElectronDataAnalyzer.cc:365
GsfElectronDataAnalyzer::hcalIso03Depth1MaxEndcaps_
double hcalIso03Depth1MaxEndcaps_
Definition: GsfElectronDataAnalyzer.cc:377
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
GsfElectronDataAnalyzer::h_ele_dEtaEleCl_propOut_endcaps
TH1F * h_ele_dEtaEleCl_propOut_endcaps
Definition: GsfElectronDataAnalyzer.cc:282
GsfElectronDataAnalyzer::detamin
double detamin
Definition: GsfElectronDataAnalyzer.cc:391
BeamSpot.h
GsfElectronDataAnalyzer::detamatchmax
double detamatchmax
Definition: GsfElectronDataAnalyzer.cc:396
GsfElectronDataAnalyzer::h_ele_EoPVsPhi
TH2F * h_ele_EoPVsPhi
Definition: GsfElectronDataAnalyzer.cc:235
GsfElectronDataAnalyzer::h_ele_matchingObjectAbsEta_matched
TH1F * h_ele_matchingObjectAbsEta_matched
Definition: GsfElectronDataAnalyzer.cc:106
GsfElectronDataAnalyzer::h_ele_dPhiScVsPt_propVtx
TH2F * h_ele_dPhiScVsPt_propVtx
Definition: GsfElectronDataAnalyzer.cc:267
GsfElectronDataAnalyzer::histfile_
TFile * histfile_
Definition: GsfElectronDataAnalyzer.cc:90
GsfElectronDataAnalyzer::nbindetamatch2D
int nbindetamatch2D
Definition: GsfElectronDataAnalyzer.cc:425
reco::LeafCandidate::py
double py() const final
y coordinate of momentum vector
Definition: LeafCandidate.h:142
L1TEGammaOffline_cfi.triggerNames
triggerNames
Definition: L1TEGammaOffline_cfi.py:40
GsfElectronDataAnalyzer::h_ele_outerP
TH1F * h_ele_outerP
Definition: GsfElectronDataAnalyzer.cc:223
GsfElectronDataAnalyzer::hcalIso03Depth2MaxEndcaps_
double hcalIso03Depth2MaxEndcaps_
Definition: GsfElectronDataAnalyzer.cc:378
GsfElectronDataAnalyzer::h_ele_EeleOPoutVsE
TH2F * h_ele_EeleOPoutVsE
Definition: GsfElectronDataAnalyzer.cc:254
GsfElectronDataAnalyzer::h_ele_hcalDepth1TowerSumEt_dr03
TH1F * h_ele_hcalDepth1TowerSumEt_dr03
Definition: GsfElectronDataAnalyzer.cc:338
GsfElectronDataAnalyzer::isEcalDriven_
bool isEcalDriven_
Definition: GsfElectronDataAnalyzer.cc:352
GsfElectronDataAnalyzer::phimin
double phimin
Definition: GsfElectronDataAnalyzer.cc:385
GsfElectronDataAnalyzer::minEt_
double minEt_
Definition: GsfElectronDataAnalyzer.cc:346
GsfElectronDataAnalyzer::h_ele_dPhiEleCl_propOut
TH1F * h_ele_dPhiEleCl_propOut
Definition: GsfElectronDataAnalyzer.cc:286
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
GsfElectronDataAnalyzer::h_ele_dPhiCl_propOut
TH1F * h_ele_dPhiCl_propOut
Definition: GsfElectronDataAnalyzer.cc:274
reco::BeamSpot
Definition: BeamSpot.h:21
GsfElectronDataAnalyzer::histSclEoEmatchingObjectShowering_barrel
TH1F * histSclEoEmatchingObjectShowering_barrel
Definition: GsfElectronDataAnalyzer.cc:330
GsfElectronDataAnalyzer::h_ele_vertexX
TH1F * h_ele_vertexX
Definition: GsfElectronDataAnalyzer.cc:137
edm::ESHandle< TrackerGeometry >
GsfElectronDataAnalyzer::h_ele_dPhiSc_propVtx_barrel
TH1F * h_ele_dPhiSc_propVtx_barrel
Definition: GsfElectronDataAnalyzer.cc:263
GsfElectronDataAnalyzer::mvaMin_
double mvaMin_
Definition: GsfElectronDataAnalyzer.cc:372
GsfElectronDataAnalyzer::histSclE1x5_endcaps_
TH1F * histSclE1x5_endcaps_
Definition: GsfElectronDataAnalyzer.cc:166
GsfElectronDataAnalyzer::h_ele_lostHits
TH1F * h_ele_lostHits
Definition: GsfElectronDataAnalyzer.cc:184
GsfElectronDataAnalyzer::histSclSigEtaEta_
TH1F * histSclSigEtaEta_
Definition: GsfElectronDataAnalyzer.cc:158
GsfElectronDataAnalyzer::maxAbsEta_
double maxAbsEta_
Definition: GsfElectronDataAnalyzer.cc:348
GsfElectronDataAnalyzer::h_ele_outerPVsEta_mode
TH2F * h_ele_outerPVsEta_mode
Definition: GsfElectronDataAnalyzer.cc:225
GsfElectronDataAnalyzer::h_ele_HoE_barrel
TH1F * h_ele_HoE_barrel
Definition: GsfElectronDataAnalyzer.cc:309
GsfElectronDataAnalyzer::seedPhi
float seedPhi[10]
Definition: GsfElectronDataAnalyzer.cc:94
reco::GsfElectron
Definition: GsfElectron.h:34
GsfElectronDataAnalyzer::pmax
double pmax
Definition: GsfElectronDataAnalyzer.cc:388
GsfElectronDataAnalyzer::h_ele_EoPoutVsEta
TH2F * h_ele_EoPoutVsEta
Definition: GsfElectronDataAnalyzer.cc:246
GsfElectron.h
GsfElectronDataAnalyzer::isEE_
bool isEE_
Definition: GsfElectronDataAnalyzer.cc:350
GsfElectronDataAnalyzer::histSclEoEmatchingObjectGolden_endcaps
TH1F * histSclEoEmatchingObjectGolden_endcaps
Definition: GsfElectronDataAnalyzer.cc:329
GsfElectronDataAnalyzer::h_ele_outerPtVsPt_mode
TH2F * h_ele_outerPtVsPt_mode
Definition: GsfElectronDataAnalyzer.cc:230
GsfElectronDataAnalyzer::h_ele_mee_os_ebee
TH1F * h_ele_mee_os_ebee
Definition: GsfElectronDataAnalyzer.cc:114
GsfElectronDataAnalyzer::histSclSigIEtaIEta_
TH1F * histSclSigIEtaIEta_
Definition: GsfElectronDataAnalyzer.cc:161
GsfElectronDataAnalyzer::histSclEoEmatchingObject_barrel_matched
TH1F * histSclEoEmatchingObject_barrel_matched
Definition: GsfElectronDataAnalyzer.cc:148
GsfElectronDataAnalyzer::histSclEoEmatchingObject_endcaps_new_matched
TH1F * histSclEoEmatchingObject_endcaps_new_matched
Definition: GsfElectronDataAnalyzer.cc:151
GsfElectronDataAnalyzer::h_ele_dPhiClVsPt_propOut
TH2F * h_ele_dPhiClVsPt_propOut
Definition: GsfElectronDataAnalyzer.cc:279
PbPb_ZMuSkimMuonDPG_cff.deltaR
deltaR
Definition: PbPb_ZMuSkimMuonDPG_cff.py:63
GsfElectronDataAnalyzer::dEtaMaxBarrel_
double dEtaMaxBarrel_
Definition: GsfElectronDataAnalyzer.cc:359
GsfElectronDataAnalyzer::h_ele_dEtaScVsPhi_propVtx
TH2F * h_ele_dEtaScVsPhi_propVtx
Definition: GsfElectronDataAnalyzer.cc:260
GsfElectronDataAnalyzer::h_ele_dPhiCl_propOut_endcaps
TH1F * h_ele_dPhiCl_propOut_endcaps
Definition: GsfElectronDataAnalyzer.cc:276
GsfElectronDataAnalyzer::mcQ
float mcQ[10]
Definition: GsfElectronDataAnalyzer.cc:92
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
GsfElectronDataAnalyzer::h_ele_dPhiEleClVsPhi_propOut
TH2F * h_ele_dPhiEleClVsPhi_propOut
Definition: GsfElectronDataAnalyzer.cc:290
GsfElectronDataAnalyzer::h_ele_PinVsPoutShowering_mean
TH2F * h_ele_PinVsPoutShowering_mean
Definition: GsfElectronDataAnalyzer.cc:323
GsfElectronDataAnalyzer::h_ele_mee_os_eeee
TH1F * h_ele_mee_os_eeee
Definition: GsfElectronDataAnalyzer.cc:115
GsfElectronDataAnalyzer::histSclEtaVsPhi_
TH2F * histSclEtaVsPhi_
Definition: GsfElectronDataAnalyzer.cc:155
GsfElectronDataAnalyzer::h_ele_outerPt_mode
TH1F * h_ele_outerPt_mode
Definition: GsfElectronDataAnalyzer.cc:227
GsfElectronDataAnalyzer::HLTPathsByName_
std::vector< std::string > HLTPathsByName_
Definition: GsfElectronDataAnalyzer.cc:72
GsfElectronDataAnalyzer::h_ele_EseedOP_endcaps
TH1F * h_ele_EseedOP_endcaps
Definition: GsfElectronDataAnalyzer.cc:239
GsfElectronDataAnalyzer::h_ele_chi2VsEta
TH2F * h_ele_chi2VsEta
Definition: GsfElectronDataAnalyzer.cc:193
GsfElectronDataAnalyzer::h_ele_outerPtVsPhi_mode
TH2F * h_ele_outerPtVsPhi_mode
Definition: GsfElectronDataAnalyzer.cc:229
GsfElectronDataAnalyzer::h_ele_vertexPtVsEta
TH2F * h_ele_vertexPtVsEta
Definition: GsfElectronDataAnalyzer.cc:130
GsfElectronDataAnalyzer::endJob
void endJob() override
Definition: GsfElectronDataAnalyzer.cc:1400
GsfElectronDataAnalyzer::nbinlhits
int nbinlhits
Definition: GsfElectronDataAnalyzer.cc:419
GsfElectronDataAnalyzer::h_ele_EeleOPoutVsPhi
TH2F * h_ele_EeleOPoutVsPhi
Definition: GsfElectronDataAnalyzer.cc:253
GsfElectronDataAnalyzer::h_ele_PtinVsPtoutShowering_mean
TH2F * h_ele_PtinVsPtoutShowering_mean
Definition: GsfElectronDataAnalyzer.cc:327
GsfElectronDataAnalyzer::hoemin
double hoemin
Definition: GsfElectronDataAnalyzer.cc:405
GsfElectronDataAnalyzer::h_ele_vertexY
TH1F * h_ele_vertexY
Definition: GsfElectronDataAnalyzer.cc:138
GsfElectronDataAnalyzer::isEB_
bool isEB_
Definition: GsfElectronDataAnalyzer.cc:349
GsfElectronDataAnalyzer::h_ele_vertexTIPVsPt
TH2F * h_ele_vertexTIPVsPt
Definition: GsfElectronDataAnalyzer.cc:143
GsfElectronDataAnalyzer::h_ele_EoPout_barrel
TH1F * h_ele_EoPout_barrel
Definition: GsfElectronDataAnalyzer.cc:244
GsfElectronDataAnalyzer::mcEnergy
float mcEnergy[10]
Definition: GsfElectronDataAnalyzer.cc:92
GsfElectronDataAnalyzer::h_ele_dEtaSc_propVtx_endcaps
TH1F * h_ele_dEtaSc_propVtx_endcaps
Definition: GsfElectronDataAnalyzer.cc:258
GsfElectronFwd.h
GsfElectronDataAnalyzer::nbinp2D
int nbinp2D
Definition: GsfElectronDataAnalyzer.cc:412
edm::ParameterSet
Definition: ParameterSet.h:47
GsfElectronDataAnalyzer::dPhiMinBarrel_
double dPhiMinBarrel_
Definition: GsfElectronDataAnalyzer.cc:362
GsfElectronDataAnalyzer::h_ele_vertexPtVsPhi
TH2F * h_ele_vertexPtVsPhi
Definition: GsfElectronDataAnalyzer.cc:131
GsfElectronDataAnalyzer::h_ele_HoE
TH1F * h_ele_HoE
Definition: GsfElectronDataAnalyzer.cc:308
GsfElectronDataAnalyzer::h_ele_PhiMnPhimatchingObjectVsEta_matched
TH2F * h_ele_PhiMnPhimatchingObjectVsEta_matched
Definition: GsfElectronDataAnalyzer.cc:212
GsfElectronDataAnalyzer::nbinpt
int nbinpt
Definition: GsfElectronDataAnalyzer.cc:409
GsfElectronDataAnalyzer::h_ele_dEtaScVsEta_propVtx
TH2F * h_ele_dEtaScVsEta_propVtx
Definition: GsfElectronDataAnalyzer.cc:259
Event.h
reco::LeafCandidate::eta
double eta() const final
momentum pseudorapidity
Definition: LeafCandidate.h:152
GsfElectronDataAnalyzer::h_ele_dEtaClVsPhi_propOut
TH2F * h_ele_dEtaClVsPhi_propOut
Definition: GsfElectronDataAnalyzer.cc:272
GsfElectronDataAnalyzer::nbinphi
int nbinphi
Definition: GsfElectronDataAnalyzer.cc:411
GsfElectronDataAnalyzer::h_ele_seed_dphi2VsPt_
TH2F * h_ele_seed_dphi2VsPt_
Definition: GsfElectronDataAnalyzer.cc:295
GsfElectronDataAnalyzer::h_ele_outerPtVsEta_mode
TH2F * h_ele_outerPtVsEta_mode
Definition: GsfElectronDataAnalyzer.cc:228
GsfElectronDataAnalyzer::h_ele_vertexEtaVsPhi
TH2F * h_ele_vertexEtaVsPhi
Definition: GsfElectronDataAnalyzer.cc:134
GsfElectronDataAnalyzer::h_ele_vertexTIP
TH1F * h_ele_vertexTIP
Definition: GsfElectronDataAnalyzer.cc:140
GsfElectronDataAnalyzer::h_ele_dEtaCl_propOut_barrel
TH1F * h_ele_dEtaCl_propOut_barrel
Definition: GsfElectronDataAnalyzer.cc:269
GsfElectronDataAnalyzer::h_ele_PinMnPoutVsPhi_mode
TH2F * h_ele_PinMnPoutVsPhi_mode
Definition: GsfElectronDataAnalyzer.cc:218
GsfElectronDataAnalyzer::tkIso03Max_
double tkIso03Max_
Definition: GsfElectronDataAnalyzer.cc:375
GsfElectronDataAnalyzer::nAfterTrigger_
unsigned int nAfterTrigger_
Definition: GsfElectronDataAnalyzer.cc:67
GsfElectronDataAnalyzer::h_ele_foundHitsVsPhi
TH2F * h_ele_foundHitsVsPhi
Definition: GsfElectronDataAnalyzer.cc:182
GsfElectronDataAnalyzer::h_ele_EoP_barrel
TH1F * h_ele_EoP_barrel
Definition: GsfElectronDataAnalyzer.cc:232
TrajectoryStateTransform
Definition: TrajectoryStateTransform.h:47
GsfElectronDataAnalyzer::histSclPhi_
TH1F * histSclPhi_
Definition: GsfElectronDataAnalyzer.cc:157
GsfElectronDataAnalyzer::hadronicOverEmMaxEndcaps_
double hadronicOverEmMaxEndcaps_
Definition: GsfElectronDataAnalyzer.cc:371
GsfElectronDataAnalyzer::h_ele_PtinVsPtoutGolden_mode
TH2F * h_ele_PtinVsPtoutGolden_mode
Definition: GsfElectronDataAnalyzer.cc:324
GsfElectronDataAnalyzer::pDD
edm::ESHandle< TrackerGeometry > pDD
Definition: GsfElectronDataAnalyzer.cc:87
GsfElectronDataAnalyzer::h_ele_foundHitsVsPt
TH2F * h_ele_foundHitsVsPt
Definition: GsfElectronDataAnalyzer.cc:183
GsfElectronDataAnalyzer::h_ele_EseedOP
TH1F * h_ele_EseedOP
Definition: GsfElectronDataAnalyzer.cc:237
GsfElectronDataAnalyzer::h_ele_PinVsPoutShowering_mode
TH2F * h_ele_PinVsPoutShowering_mode
Definition: GsfElectronDataAnalyzer.cc:321
GsfElectronDataAnalyzer::nbinphi2D
int nbinphi2D
Definition: GsfElectronDataAnalyzer.cc:415
GsfElectronDataAnalyzer::histSclE1x5_
TH1F * histSclE1x5_
Definition: GsfElectronDataAnalyzer.cc:164
GsfElectronDataAnalyzer::h_ele_EseedOP_barrel
TH1F * h_ele_EseedOP_barrel
Definition: GsfElectronDataAnalyzer.cc:238
TriggerNames.h
GsfElectronDataAnalyzer::superclusterEnergy
float superclusterEnergy[10]
Definition: GsfElectronDataAnalyzer.cc:93
iEvent
int iEvent
Definition: GenABIO.cc:224
GsfElectronDataAnalyzer::nEvents_
unsigned int nEvents_
Definition: GsfElectronDataAnalyzer.cc:66
GsfElectronDataAnalyzer::type_
std::string type_
Definition: GsfElectronDataAnalyzer.cc:79
GsfElectronDataAnalyzer::transformer_
TrajectoryStateTransform transformer_
Definition: GsfElectronDataAnalyzer.cc:86
GsfElectronDataAnalyzer::h_ele_ecalRecHitSumEt_dr04
TH1F * h_ele_ecalRecHitSumEt_dr04
Definition: GsfElectronDataAnalyzer.cc:341
GsfElectronDataAnalyzer::h_ele_eta
TH1F * h_ele_eta
Definition: GsfElectronDataAnalyzer.cc:302
GsfTrack.h
GsfElectronDataAnalyzer::h_ele_dPhiEleCl_propOut_endcaps
TH1F * h_ele_dPhiEleCl_propOut_endcaps
Definition: GsfElectronDataAnalyzer.cc:288
GsfElectronDataAnalyzer::histSclE2x5max_
TH1F * histSclE2x5max_
Definition: GsfElectronDataAnalyzer.cc:167
GsfElectronDataAnalyzer::superclusterPhi
float superclusterPhi[10]
Definition: GsfElectronDataAnalyzer.cc:93
GsfElectronDataAnalyzer::h_ele_chargeVsPhi
TH2F * h_ele_chargeVsPhi
Definition: GsfElectronDataAnalyzer.cc:125
GsfElectronDataAnalyzer::histSclEtVsPhi_
TH2F * histSclEtVsPhi_
Definition: GsfElectronDataAnalyzer.cc:154
GsfElectronDataAnalyzer::ptmax
double ptmax
Definition: GsfElectronDataAnalyzer.cc:387
GsfElectronDataAnalyzer::h_ele_dEtaScVsPt_propVtx
TH2F * h_ele_dEtaScVsPt_propVtx
Definition: GsfElectronDataAnalyzer.cc:261
edm::EventSetup
Definition: EventSetup.h:58
GsfElectronDataAnalyzer::h_ele_seed_dphi2VsEta_
TH2F * h_ele_seed_dphi2VsEta_
Definition: GsfElectronDataAnalyzer.cc:294
GsfElectronDataAnalyzer::h_ele_hcalDepth1TowerSumEt_dr04
TH1F * h_ele_hcalDepth1TowerSumEt_dr04
Definition: GsfElectronDataAnalyzer.cc:342
GsfElectronDataAnalyzer::ecalIso03MaxEndcaps_
double ecalIso03MaxEndcaps_
Definition: GsfElectronDataAnalyzer.cc:380
GsfElectronDataAnalyzer::h_ele_outerPt
TH1F * h_ele_outerPt
Definition: GsfElectronDataAnalyzer.cc:226
GsfElectronDataAnalyzer::h_ele_lostHitsVsPt
TH2F * h_ele_lostHitsVsPt
Definition: GsfElectronDataAnalyzer.cc:189
GsfElectronDataAnalyzer::beamSpot_
edm::InputTag beamSpot_
Definition: GsfElectronDataAnalyzer.cc:70
GsfElectronDataAnalyzer::h_ele_vertexP
TH1F * h_ele_vertexP
Definition: GsfElectronDataAnalyzer.cc:127
GsfElectronDataAnalyzer::sigIetaIetaMinEndcaps_
double sigIetaIetaMinEndcaps_
Definition: GsfElectronDataAnalyzer.cc:368
GsfElectronDataAnalyzer::histSclEoEmatchingObject_barrel_new_matched
TH1F * histSclEoEmatchingObject_barrel_new_matched
Definition: GsfElectronDataAnalyzer.cc:150
GsfElectronDataAnalyzer::nbinpteff
int nbinpteff
Definition: GsfElectronDataAnalyzer.cc:410
GsfElectronDataAnalyzer::h_ele_eta_narrow
TH1F * h_ele_eta_narrow
Definition: GsfElectronDataAnalyzer.cc:305
GsfElectronDataAnalyzer::histSclEta_
TH1F * histSclEta_
Definition: GsfElectronDataAnalyzer.cc:156
GsfElectronDataAnalyzer::h_matchingObjectPhi
TH1F * h_matchingObjectPhi
Definition: GsfElectronDataAnalyzer.cc:102
GsfElectronDataAnalyzer::histSclEoEmatchingObject_endcaps_matched
TH1F * histSclEoEmatchingObject_endcaps_matched
Definition: GsfElectronDataAnalyzer.cc:149
GsfElectronDataAnalyzer::histSclSigIEtaIEta_endcaps_
TH1F * histSclSigIEtaIEta_endcaps_
Definition: GsfElectronDataAnalyzer.cc:163
GsfElectronDataAnalyzer::h_ele_EoPout
TH1F * h_ele_EoPout
Definition: GsfElectronDataAnalyzer.cc:243
GsfElectronDataAnalyzer::h_ele_EeleOPout
TH1F * h_ele_EeleOPout
Definition: GsfElectronDataAnalyzer.cc:249
InputTag.h
GsfElectronDataAnalyzer::h_ele_lostHits_barrel
TH1F * h_ele_lostHits_barrel
Definition: GsfElectronDataAnalyzer.cc:185
GsfElectronDataAnalyzer::h_ele_ambiguousTracksVsPhi
TH2F * h_ele_ambiguousTracksVsPhi
Definition: GsfElectronDataAnalyzer.cc:176
GsfElectronDataAnalyzer::nbindphi
int nbindphi
Definition: GsfElectronDataAnalyzer.cc:422
GsfElectronDataAnalyzer::h_ele_EseedOPVsPhi
TH2F * h_ele_EseedOPVsPhi
Definition: GsfElectronDataAnalyzer.cc:241
GsfElectronDataAnalyzer::h_ele_PhiMnPhimatchingObject2_matched
TH1F * h_ele_PhiMnPhimatchingObject2_matched
Definition: GsfElectronDataAnalyzer.cc:211
GsfElectronDataAnalyzer::isTrackerDriven_
bool isTrackerDriven_
Definition: GsfElectronDataAnalyzer.cc:353
GsfElectronDataAnalyzer::h_ele_EoP_endcaps
TH1F * h_ele_EoP_endcaps
Definition: GsfElectronDataAnalyzer.cc:233
GsfElectronDataAnalyzer::h_matchingObjectPt
TH1F * h_matchingObjectPt
Definition: GsfElectronDataAnalyzer.cc:101
GsfElectronDataAnalyzer::histSclEn_
TH1F * histSclEn_
Definition: GsfElectronDataAnalyzer.cc:147
GsfElectronDataAnalyzer::sigIetaIetaMaxEndcaps_
double sigIetaIetaMaxEndcaps_
Definition: GsfElectronDataAnalyzer.cc:369
GsfElectronDataAnalyzer::h_ele_EoPVsEta
TH2F * h_ele_EoPVsEta
Definition: GsfElectronDataAnalyzer.cc:234
GsfElectronDataAnalyzer::maxPtMatchingObject_
double maxPtMatchingObject_
Definition: GsfElectronDataAnalyzer.cc:82
GsfElectronDataAnalyzer::h_ele_dEtaEleClVsPt_propOut
TH2F * h_ele_dEtaEleClVsPt_propOut
Definition: GsfElectronDataAnalyzer.cc:285
GsfElectronDataAnalyzer::dphimatchmin
double dphimatchmin
Definition: GsfElectronDataAnalyzer.cc:397
GsfElectronDataAnalyzer::seedMomentum
float seedMomentum[10]
Definition: GsfElectronDataAnalyzer.cc:94
GsfElectronDataAnalyzer::h_ele_PinMnPoutVsE_mode
TH2F * h_ele_PinMnPoutVsE_mode
Definition: GsfElectronDataAnalyzer.cc:220
GsfElectronDataAnalyzer::histSclE1x5_barrel_
TH1F * histSclE1x5_barrel_
Definition: GsfElectronDataAnalyzer.cc:165
GsfElectronDataAnalyzer::triggerResults_
edm::InputTag triggerResults_
Definition: GsfElectronDataAnalyzer.cc:69
GsfElectronDataAnalyzer::histSclEt_
TH1F * histSclEt_
Definition: GsfElectronDataAnalyzer.cc:152
GsfElectronDataAnalyzer::trigger
bool trigger(const edm::Event &e)
Definition: GsfElectronDataAnalyzer.cc:2348
GsfElectronDataAnalyzer::dPhiMinEndcaps_
double dPhiMinEndcaps_
Definition: GsfElectronDataAnalyzer.cc:364
GsfElectronDataAnalyzer::mcPt
float mcPt[10]
Definition: GsfElectronDataAnalyzer.cc:92
GsfElectronDataAnalyzer::histSclE2x5max_endcaps_
TH1F * histSclE2x5max_endcaps_
Definition: GsfElectronDataAnalyzer.cc:169
GsfElectronDataAnalyzer::meemin
double meemin
Definition: GsfElectronDataAnalyzer.cc:403
GsfElectronDataAnalyzer::dEtaMinBarrel_
double dEtaMinBarrel_
Definition: GsfElectronDataAnalyzer.cc:358
GsfElectronDataAnalyzer::nbinp
int nbinp
Definition: GsfElectronDataAnalyzer.cc:408
GsfElectronDataAnalyzer::h_ele_tkSumPt_dr03
TH1F * h_ele_tkSumPt_dr03
Definition: GsfElectronDataAnalyzer.cc:336
GsfElectronDataAnalyzer::histSclEoEmatchingObjectGolden_barrel
TH1F * histSclEoEmatchingObjectGolden_barrel
Definition: GsfElectronDataAnalyzer.cc:328
GsfElectronDataAnalyzer::nbinpoptrue
int nbinpoptrue
Definition: GsfElectronDataAnalyzer.cc:427
GsfElectronDataAnalyzer::nbineta2D
int nbineta2D
Definition: GsfElectronDataAnalyzer.cc:414
GsfElectronDataAnalyzer::h_ele_eta_shower
TH1F * h_ele_eta_shower
Definition: GsfElectronDataAnalyzer.cc:306
reco::LeafCandidate::phi
double phi() const final
momentum azimuthal angle
Definition: LeafCandidate.h:148
SuperClusterFwd.h
GsfElectronDataAnalyzer::h_ele_chi2VsPt
TH2F * h_ele_chi2VsPt
Definition: GsfElectronDataAnalyzer.cc:195
GsfElectronDataAnalyzer::h_ele_seed_dphi2_
TH1F * h_ele_seed_dphi2_
Definition: GsfElectronDataAnalyzer.cc:293
GsfElectronDataAnalyzer::h_ele_fbremVsEta_mean
TProfile * h_ele_fbremVsEta_mean
Definition: GsfElectronDataAnalyzer.cc:318
GsfElectronDataAnalyzer::h_ele_EseedOPVsE
TH2F * h_ele_EseedOPVsE
Definition: GsfElectronDataAnalyzer.cc:242
GsfElectronDataAnalyzer::h_ele_outerP_mode
TH1F * h_ele_outerP_mode
Definition: GsfElectronDataAnalyzer.cc:224
GsfElectronDataAnalyzer::histNum_
TH1F * histNum_
Definition: GsfElectronDataAnalyzer.cc:145
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
GsfElectronDataAnalyzer::h_ele_dPhiClVsEta_propOut
TH2F * h_ele_dPhiClVsEta_propOut
Definition: GsfElectronDataAnalyzer.cc:277
GsfElectronDataAnalyzer::eopmaxsht
double eopmaxsht
Definition: GsfElectronDataAnalyzer.cc:390
GsfElectronDataAnalyzer::h_ele_PinMnPout_mode
TH1F * h_ele_PinMnPout_mode
Definition: GsfElectronDataAnalyzer.cc:216
GsfElectronDataAnalyzer::dEtaMinEndcaps_
double dEtaMinEndcaps_
Definition: GsfElectronDataAnalyzer.cc:360
GsfElectronDataAnalyzer::h_ele_foundHits_barrel
TH1F * h_ele_foundHits_barrel
Definition: GsfElectronDataAnalyzer.cc:179
reco::LeafCandidate::p
double p() const final
magnitude of momentum vector
Definition: LeafCandidate.h:123
GsfElectronDataAnalyzer::h_ele_chi2
TH1F * h_ele_chi2
Definition: GsfElectronDataAnalyzer.cc:190
GsfElectronDataAnalyzer::h_ele_foundHitsVsEta
TH2F * h_ele_foundHitsVsEta
Definition: GsfElectronDataAnalyzer.cc:181
edm::TriggerNames
Definition: TriggerNames.h:55
GsfElectronDataAnalyzer::h_ele_lostHitsVsEta
TH2F * h_ele_lostHitsVsEta
Definition: GsfElectronDataAnalyzer.cc:187
GsfElectronDataAnalyzer::seedPt
float seedPt[10]
Definition: GsfElectronDataAnalyzer.cc:94
GsfElectronDataAnalyzer::h_ele_PinMnPoutVsPt_mode
TH2F * h_ele_PinMnPoutVsPt_mode
Definition: GsfElectronDataAnalyzer.cc:219
GsfElectronDataAnalyzer::nbinxyz
int nbinxyz
Definition: GsfElectronDataAnalyzer.cc:420
GsfElectronDataAnalyzer::h_ele_mee_os
TH1F * h_ele_mee_os
Definition: GsfElectronDataAnalyzer.cc:112
GsfElectronDataAnalyzer::h_ele_matchingObjectZ_matched
TH1F * h_ele_matchingObjectZ_matched
Definition: GsfElectronDataAnalyzer.cc:109
EventSetup.h
GsfElectronDataAnalyzer::h_ele_PtinVsPtoutGolden_mean
TH2F * h_ele_PtinVsPtoutGolden_mean
Definition: GsfElectronDataAnalyzer.cc:326
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
edm::RefToBase< TrajectorySeed >
GsfElectronDataAnalyzer::nbindphimatch
int nbindphimatch
Definition: GsfElectronDataAnalyzer.cc:424
GsfElectronDataAnalyzer::h_ele_ambiguousTracks
TH1F * h_ele_ambiguousTracks
Definition: GsfElectronDataAnalyzer.cc:174
TrajectoryStateTransform.h
reco::GsfElectron::superCluster
SuperClusterRef superCluster() const override
reference to a SuperCluster
Definition: GsfElectron.h:162
dqm-mbProfile.log
log
Definition: dqm-mbProfile.py:17
GsfElectronDataAnalyzer::h_ele_dPhiEleClVsPt_propOut
TH2F * h_ele_dPhiEleClVsPt_propOut
Definition: GsfElectronDataAnalyzer.cc:291
GsfElectronDataAnalyzer::h_ele_provenance
TH1F * h_ele_provenance
Definition: GsfElectronDataAnalyzer.cc:334
GsfElectronDataAnalyzer::h_ele_PinVsPoutGolden_mean
TH2F * h_ele_PinVsPoutGolden_mean
Definition: GsfElectronDataAnalyzer.cc:322
GsfElectronDataAnalyzer::h_ele_charge
TH1F * h_ele_charge
Definition: GsfElectronDataAnalyzer.cc:123
GsfElectronDataAnalyzer::h_ele_dEtaEleCl_propOut
TH1F * h_ele_dEtaEleCl_propOut
Definition: GsfElectronDataAnalyzer.cc:280
GsfElectronDataAnalyzer::h_ele_vertexEta
TH1F * h_ele_vertexEta
Definition: GsfElectronDataAnalyzer.cc:133
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:29
ztail.d
d
Definition: ztail.py:151
GsfElectronDataAnalyzer::h_ele_PtoPtmatchingObject_barrel_matched
TH1F * h_ele_PtoPtmatchingObject_barrel_matched
Definition: GsfElectronDataAnalyzer.cc:204
pi
const Double_t pi
Definition: trackSplitPlot.h:36
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
GsfElectronDataAnalyzer::hadronicOverEmMaxBarrel_
double hadronicOverEmMaxBarrel_
Definition: GsfElectronDataAnalyzer.cc:370
GsfElectronDataAnalyzer::h_ele_mva
TH1F * h_ele_mva
Definition: GsfElectronDataAnalyzer.cc:333
GsfElectronDataAnalyzer::h_ele_dEtaClVsEta_propOut
TH2F * h_ele_dEtaClVsEta_propOut
Definition: GsfElectronDataAnalyzer.cc:271
ParameterSet.h
GsfElectronDataAnalyzer::histSclE2x5max_barrel_
TH1F * histSclE2x5max_barrel_
Definition: GsfElectronDataAnalyzer.cc:168
GsfElectronDataAnalyzer::h_ele_dPhiCl_propOut_barrel
TH1F * h_ele_dPhiCl_propOut_barrel
Definition: GsfElectronDataAnalyzer.cc:275
GsfElectronDataAnalyzer::nbinhoe
int nbinhoe
Definition: GsfElectronDataAnalyzer.cc:429
GsfElectronDataAnalyzer::detamatchmin
double detamatchmin
Definition: GsfElectronDataAnalyzer.cc:395
GsfElectronDataAnalyzer::eopmax
double eopmax
Definition: GsfElectronDataAnalyzer.cc:389
GsfElectronDataAnalyzer::h_ele_HoE_endcaps
TH1F * h_ele_HoE_endcaps
Definition: GsfElectronDataAnalyzer.cc:310
c
auto & c
Definition: CAHitNtupletGeneratorKernelsImpl.h:46
GsfElectronDataAnalyzer::h_ele_eta_golden
TH1F * h_ele_eta_golden
Definition: GsfElectronDataAnalyzer.cc:303
GsfElectronDataAnalyzer::~GsfElectronDataAnalyzer
~GsfElectronDataAnalyzer() override
Definition: GsfElectronDataAnalyzer.cc:543
GsfElectronDataAnalyzer::etamin
double etamin
Definition: GsfElectronDataAnalyzer.cc:383
GsfElectronDataAnalyzer::h_ele_mee_all
TH1F * h_ele_mee_all
Definition: GsfElectronDataAnalyzer.cc:111
GsfElectronDataAnalyzer::nbinpt2D
int nbinpt2D
Definition: GsfElectronDataAnalyzer.cc:413
edm::Event
Definition: Event.h:73
GsfElectronDataAnalyzer::h_ele_mee_os_bb
TH1F * h_ele_mee_os_bb
Definition: GsfElectronDataAnalyzer.cc:118
GsfElectronDataAnalyzer::hcalIso03Depth1MaxBarrel_
double hcalIso03Depth1MaxBarrel_
Definition: GsfElectronDataAnalyzer.cc:376
GsfElectronDataAnalyzer::h_ele_tkSumPt_dr04
TH1F * h_ele_tkSumPt_dr04
Definition: GsfElectronDataAnalyzer.cc:340
GsfElectronDataAnalyzer::poptruemax
double poptruemax
Definition: GsfElectronDataAnalyzer.cc:402
GsfElectronDataAnalyzer::h_ele_PoPmatchingObjectVsEta_matched
TH2F * h_ele_PoPmatchingObjectVsEta_matched
Definition: GsfElectronDataAnalyzer.cc:199
GsfElectronDataAnalyzer::isNotEBEEGap_
bool isNotEBEEGap_
Definition: GsfElectronDataAnalyzer.cc:351
GsfElectronDataAnalyzer::dPhiMaxBarrel_
double dPhiMaxBarrel_
Definition: GsfElectronDataAnalyzer.cc:363
GsfElectronDataAnalyzer::h_ele_PoPmatchingObject_barrel_matched
TH1F * h_ele_PoPmatchingObject_barrel_matched
Definition: GsfElectronDataAnalyzer.cc:202
dttmaxenums::R
Definition: DTTMax.h:29
GsfElectronDataAnalyzer::tipMaxEndcaps_
double tipMaxEndcaps_
Definition: GsfElectronDataAnalyzer.cc:374
GsfElectronDataAnalyzer::h_ele_matchingObjectPhi_matched
TH1F * h_ele_matchingObjectPhi_matched
Definition: GsfElectronDataAnalyzer.cc:108
GsfElectronDataAnalyzer::seedQ
float seedQ[10]
Definition: GsfElectronDataAnalyzer.cc:94
GsfElectronDataAnalyzer::lhitsmax
double lhitsmax
Definition: GsfElectronDataAnalyzer.cc:400
GsfElectronDataAnalyzer::h_ele_dPhiClVsPhi_propOut
TH2F * h_ele_dPhiClVsPhi_propOut
Definition: GsfElectronDataAnalyzer.cc:278
edm::InputTag
Definition: InputTag.h:15
GsfElectronDataAnalyzer::h_ele_HoE_fiducial
TH1F * h_ele_HoE_fiducial
Definition: GsfElectronDataAnalyzer.cc:311
GsfElectronDataAnalyzer::h_ele_EtaMnEtamatchingObjectVsPhi_matched
TH2F * h_ele_EtaMnEtamatchingObjectVsPhi_matched
Definition: GsfElectronDataAnalyzer.cc:208
GsfElectronDataAnalyzer::h_ele_seed_drz2VsPt_
TH2F * h_ele_seed_drz2VsPt_
Definition: GsfElectronDataAnalyzer.cc:298
GsfElectronDataAnalyzer::h_ele_EoPoutVsPhi
TH2F * h_ele_EoPoutVsPhi
Definition: GsfElectronDataAnalyzer.cc:247
GsfElectronDataAnalyzer::h_ele_ecalRecHitSumEt_dr03
TH1F * h_ele_ecalRecHitSumEt_dr03
Definition: GsfElectronDataAnalyzer.cc:337
GsfElectronDataAnalyzer::sigIetaIetaMaxBarrel_
double sigIetaIetaMaxBarrel_
Definition: GsfElectronDataAnalyzer.cc:367
GsfElectronDataAnalyzer::h_ele_EtaMnEtamatchingObjectVsPt_matched
TH2F * h_ele_EtaMnEtamatchingObjectVsPt_matched
Definition: GsfElectronDataAnalyzer.cc:209
GsfElectronDataAnalyzer::fhitsmax
double fhitsmax
Definition: GsfElectronDataAnalyzer.cc:399
GsfElectronDataAnalyzer::h_ele_PtoPtmatchingObject_matched
TH1F * h_ele_PtoPtmatchingObject_matched
Definition: GsfElectronDataAnalyzer.cc:197
GsfElectronDataAnalyzer::mcPhi
float mcPhi[10]
Definition: GsfElectronDataAnalyzer.cc:92
GsfElectronDataAnalyzer::nbindetamatch
int nbindetamatch
Definition: GsfElectronDataAnalyzer.cc:423
GsfElectronDataAnalyzer::h_ele_EoPout_endcaps
TH1F * h_ele_EoPout_endcaps
Definition: GsfElectronDataAnalyzer.cc:245
muonDTDigis_cfi.pset
pset
Definition: muonDTDigis_cfi.py:27
GsfElectronDataAnalyzer::hoemax
double hoemax
Definition: GsfElectronDataAnalyzer.cc:406
ElectronSeed.h
GsfElectronDataAnalyzer::histSclSigEtaEta_endcaps_
TH1F * histSclSigEtaEta_endcaps_
Definition: GsfElectronDataAnalyzer.cc:160
GsfElectronDataAnalyzer::h_ele_vertexPt_5100
TH1F * h_ele_vertexPt_5100
Definition: GsfElectronDataAnalyzer.cc:132
GsfElectronDataAnalyzer::analyze
void analyze(const edm::Event &e, const edm::EventSetup &c) override
Definition: GsfElectronDataAnalyzer.cc:1858
GsfElectronDataAnalyzer::h_ele_dEtaCl_propOut_endcaps
TH1F * h_ele_dEtaCl_propOut_endcaps
Definition: GsfElectronDataAnalyzer.cc:270
GsfElectronDataAnalyzer::h_ele_dEtaEleClVsEta_propOut
TH2F * h_ele_dEtaEleClVsEta_propOut
Definition: GsfElectronDataAnalyzer.cc:283
GsfElectronDataAnalyzer::h_ele_mee_os_gg
TH1F * h_ele_mee_os_gg
Definition: GsfElectronDataAnalyzer.cc:116
GsfElectronDataAnalyzer::dphimatchmax
double dphimatchmax
Definition: GsfElectronDataAnalyzer.cc:398
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
GsfElectronDataAnalyzer::h_ele_PhiMnPhimatchingObjectVsPhi_matched
TH2F * h_ele_PhiMnPhimatchingObjectVsPhi_matched
Definition: GsfElectronDataAnalyzer.cc:213
GsfElectronDataAnalyzer::matchingCondition_
std::string matchingCondition_
Definition: GsfElectronDataAnalyzer.cc:78
GsfElectronDataAnalyzer::h_matchingObjectZ
TH1F * h_matchingObjectZ
Definition: GsfElectronDataAnalyzer.cc:103
GsfElectronDataAnalyzer::h_ele_dEtaSc_propVtx
TH1F * h_ele_dEtaSc_propVtx
Definition: GsfElectronDataAnalyzer.cc:256
GsfElectronDataAnalyzer::superclusterEt
float superclusterEt[10]
Definition: GsfElectronDataAnalyzer.cc:93
GsfElectronDataAnalyzer::h_ele_chi2_barrel_
TH1F * h_ele_chi2_barrel_
Definition: GsfElectronDataAnalyzer.cc:191
GsfElectronDataAnalyzer::h_ele_PoPmatchingObjectVsPt_matched
TH2F * h_ele_PoPmatchingObjectVsPt_matched
Definition: GsfElectronDataAnalyzer.cc:201
Generated for CMSSW Reference Manual by doxygen 1.8.16