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