CMS 3D CMS Logo
GsfElectronFakeAnalyzer.h
Go to the documentation of this file.
1 #ifndef GsfElectronFakeAnalyzer_h
2 #define GsfElectronFakeAnalyzer_h
3 
4 //
5 // Package: RecoEgamma/Examples
6 // Class: GsfElectronFakeAnalyzer
7 //
8 
9 //
10 // Original Author: Ursula Berthon, Claude Charlot
11 // Created: Mon Mar 27 13:22:06 CEST 2006
12 //
13 //
14 
15 
16 #include "FWCore/Framework/interface/EDAnalyzer.h"
17 #include "FWCore/Framework/interface/Event.h"
18 #include "DataFormats/Common/interface/Handle.h"
19 #include "FWCore/Framework/interface/ESHandle.h"
20 #include "FWCore/Framework/interface/EventSetup.h"
21 #include "FWCore/ParameterSet/interface/ParameterSet.h"
22 #include "FWCore/Utilities/interface/InputTag.h"
23 
24 
25 #include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
26 
27 #include "TrackingTools/TrajectoryState/interface/TrajectoryStateTransform.h"
28 
29 class MagneticField;
30 class TFile;
31 class TH1F;
32 class TH2F;
33 class TH1I;
34 class TProfile;
35 class TTree;
36 
37 class GsfElectronFakeAnalyzer : public edm::EDAnalyzer
38 {
39  public:
40 
41  explicit GsfElectronFakeAnalyzer(const edm::ParameterSet& conf);
42 
43  ~GsfElectronFakeAnalyzer() override;
44 
45  void beginJob() override;
46  void endJob() override;
47  void analyze(const edm::Event& e, const edm::EventSetup& c) override;
48 
49  private:
50 
51  TrajectoryStateTransform transformer_;
52  edm::ESHandle<TrackerGeometry> pDD;
53  edm::ESHandle<MagneticField> theMagField;
54  TFile *histfile_;
55  TTree *tree_;
56  float mcEnergy[10], mcEta[10], mcPhi[10], mcPt[10], mcQ[10];
57  float superclusterEnergy[10], superclusterEta[10], superclusterPhi[10], superclusterEt[10];
58  float seedMomentum[10], seedEta[10], seedPhi[10], seedPt[10], seedQ[10];
59 
60  TH1F *h_matchingObjectNum;
61 
62  TH1F *h_matchingObjectEta;
63  TH1F *h_matchingObjectAbsEta;
64  TH1F *h_matchingObjectP;
65  TH1F *h_matchingObjectPt;
66  TH1F *h_matchingObjectPhi;
67  TH1F *h_matchingObjectZ;
68 
69  TH1F *h_ele_matchingObjectEta_matched;
70  TH1F *h_ele_matchingObjectAbsEta_matched;
71  TH1F *h_ele_matchingObjectPt_matched;
72  TH1F *h_ele_matchingObjectPhi_matched;
73  TH1F *h_ele_matchingObjectZ_matched;
74 
75  TH1F *h_ele_EoverP_all;
76  TH1F *h_ele_EoverP_all_barrel;
77  TH1F *h_ele_EoverP_all_endcaps;
78  TH1F *h_ele_EseedOP_all;
79  TH1F *h_ele_EseedOP_all_barrel;
80  TH1F *h_ele_EseedOP_all_endcaps;
81  TH1F *h_ele_EoPout_all;
82  TH1F *h_ele_EoPout_all_barrel;
83  TH1F *h_ele_EoPout_all_endcaps;
84  TH1F *h_ele_EeleOPout_all;
85  TH1F *h_ele_EeleOPout_all_barrel;
86  TH1F *h_ele_EeleOPout_all_endcaps;
87  TH1F *h_ele_dEtaSc_propVtx_all;
88  TH1F *h_ele_dEtaSc_propVtx_all_barrel;
89  TH1F *h_ele_dEtaSc_propVtx_all_endcaps;
90  TH1F *h_ele_dPhiSc_propVtx_all;
91  TH1F *h_ele_dPhiSc_propVtx_all_barrel;
92  TH1F *h_ele_dPhiSc_propVtx_all_endcaps;
93  TH1F *h_ele_dEtaCl_propOut_all;
94  TH1F *h_ele_dEtaCl_propOut_all_barrel;
95  TH1F *h_ele_dEtaCl_propOut_all_endcaps;
96  TH1F *h_ele_dPhiCl_propOut_all;
97  TH1F *h_ele_dPhiCl_propOut_all_barrel;
98  TH1F *h_ele_dPhiCl_propOut_all_endcaps;
99  TH1F *h_ele_TIP_all;
100  TH1F *h_ele_TIP_all_barrel;
101  TH1F *h_ele_TIP_all_endcaps;
102  TH1F *h_ele_HoE_all;
103  TH1F *h_ele_HoE_all_barrel;
104  TH1F *h_ele_HoE_all_endcaps;
105  TH1F *h_ele_vertexEta_all;
106  TH1F *h_ele_vertexPt_all;
107  TH1F *h_ele_mee_all;
108  TH1F *h_ele_mee_os;
109 
110  TH2F *h_ele_E2mnE1vsMee_all;
111  TH2F *h_ele_E2mnE1vsMee_egeg_all;
112 
113  TH1F *h_ele_charge;
114  TH2F *h_ele_chargeVsEta;
115  TH2F *h_ele_chargeVsPhi;
116  TH2F *h_ele_chargeVsPt;
117  TH1F *h_ele_vertexP;
118  TH1F *h_ele_vertexPt;
119  TH2F *h_ele_vertexPtVsEta;
120  TH2F *h_ele_vertexPtVsPhi;
121  TH1F *h_ele_vertexPt_5100;
122  TH1F *h_ele_vertexEta;
123  TH2F *h_ele_vertexEtaVsPhi;
124  TH1F *h_ele_vertexAbsEta;
125  TH1F *h_ele_vertexPhi;
126  TH1F *h_ele_vertexX;
127  TH1F *h_ele_vertexY;
128  TH1F *h_ele_vertexZ;
129  TH1F *h_ele_vertexTIP;
130  TH2F *h_ele_vertexTIPVsEta;
131  TH2F *h_ele_vertexTIPVsPhi;
132  TH2F *h_ele_vertexTIPVsPt;
133 
134  TH1F *histNum_;
135 
136  TH1F *histSclEn_ ;
137  TH1F *histSclEoEmatchingObject_barrel;
138  TH1F *histSclEoEmatchingObject_endcaps;
139  TH1F *histSclEt_ ;
140  TH2F *histSclEtVsEta_ ;
141  TH2F *histSclEtVsPhi_ ;
142  TH2F *histSclEtaVsPhi_ ;
143  TH1F *histSclEta_ ;
144  TH1F *histSclPhi_ ;
145  TH1F *histSclSigEtaEta_ ;
146  TH1F *histSclSigEtaEta_barrel_ ;
147  TH1F *histSclSigEtaEta_endcaps_ ;
148  TH1F *histSclSigIEtaIEta_ ;
149  TH1F *histSclSigIEtaIEta_barrel_ ;
150  TH1F *histSclSigIEtaIEta_endcaps_ ;
151  TH1F *histSclE1x5_ ;
152  TH1F *histSclE1x5_barrel_ ;
153  TH1F *histSclE1x5_endcaps_ ;
154  TH1F *histSclE2x5max_ ;
155  TH1F *histSclE2x5max_barrel_ ;
156  TH1F *histSclE2x5max_endcaps_ ;
157  TH1F *histSclE5x5_ ;
158  TH1F *histSclE5x5_barrel_ ;
159  TH1F *histSclE5x5_endcaps_ ;
160 
161  TH1F *h_ele_ambiguousTracks;
162  TH2F *h_ele_ambiguousTracksVsEta;
163  TH2F *h_ele_ambiguousTracksVsPhi;
164  TH2F *h_ele_ambiguousTracksVsPt;
165  TH1F *h_ele_foundHits;
166  TH1F *h_ele_foundHits_barrel;
167  TH1F *h_ele_foundHits_endcaps;
168  TH2F *h_ele_foundHitsVsEta;
169  TH2F *h_ele_foundHitsVsPhi;
170  TH2F *h_ele_foundHitsVsPt;
171  TH1F *h_ele_lostHits;
172  TH1F *h_ele_lostHits_barrel;
173  TH1F *h_ele_lostHits_endcaps;
174  TH2F *h_ele_lostHitsVsEta;
175  TH2F *h_ele_lostHitsVsPhi;
176  TH2F *h_ele_lostHitsVsPt;
177  TH1F *h_ele_chi2;
178  TH1F *h_ele_chi2_barrel;
179  TH1F *h_ele_chi2_endcaps;
180  TH2F *h_ele_chi2VsEta;
181  TH2F *h_ele_chi2VsPhi;
182  TH2F *h_ele_chi2VsPt;
183 
184  TH1F *h_ele_PoPmatchingObject;
185  TH2F *h_ele_PoPmatchingObjectVsEta;
186  TH2F *h_ele_PoPmatchingObjectVsPhi;
187  TH2F *h_ele_PoPmatchingObjectVsPt;
188  TH1F *h_ele_PoPmatchingObject_barrel;
189  TH1F *h_ele_PoPmatchingObject_endcaps;
190  TH1F *h_ele_EtaMnEtamatchingObject;
191  TH2F *h_ele_EtaMnEtamatchingObjectVsEta;
192  TH2F *h_ele_EtaMnEtamatchingObjectVsPhi;
193  TH2F *h_ele_EtaMnEtamatchingObjectVsPt;
194  TH1F *h_ele_PhiMnPhimatchingObject;
195  TH1F *h_ele_PhiMnPhimatchingObject2;
196  TH2F *h_ele_PhiMnPhimatchingObjectVsEta;
197  TH2F *h_ele_PhiMnPhimatchingObjectVsPhi;
198  TH2F *h_ele_PhiMnPhimatchingObjectVsPt;
199  TH1F *h_ele_PinMnPout;
200  TH1F *h_ele_PinMnPout_mode;
201  TH2F *h_ele_PinMnPoutVsEta_mode;
202  TH2F *h_ele_PinMnPoutVsPhi_mode;
203  TH2F *h_ele_PinMnPoutVsPt_mode;
204  TH2F *h_ele_PinMnPoutVsE_mode;
205  TH2F *h_ele_PinMnPoutVsChi2_mode;
206 
207  TH1F *h_ele_outerP;
208  TH1F *h_ele_outerP_mode;
209  TH2F *h_ele_outerPVsEta_mode;
210  TH1F *h_ele_outerPt;
211  TH1F *h_ele_outerPt_mode;
212  TH2F *h_ele_outerPtVsEta_mode;
213  TH2F *h_ele_outerPtVsPhi_mode;
214  TH2F *h_ele_outerPtVsPt_mode;
215  TH1F *h_ele_EoP;
216  TH1F *h_ele_EoP_barrel;
217  TH1F *h_ele_EoP_endcaps;
218  TH1F *h_ele_EoP_eg;
219  TH1F *h_ele_EoP_eg_barrel;
220  TH1F *h_ele_EoP_eg_endcaps;
221  TH2F *h_ele_EoPVsEta;
222  TH2F *h_ele_EoPVsPhi;
223  TH2F *h_ele_EoPVsE;
224  TH1F *h_ele_EseedOP;
225  TH1F *h_ele_EseedOP_barrel;
226  TH1F *h_ele_EseedOP_endcaps;
227  TH1F *h_ele_EseedOP_eg;
228  TH1F *h_ele_EseedOP_eg_barrel;
229  TH1F *h_ele_EseedOP_eg_endcaps;
230  TH2F *h_ele_EseedOPVsEta;
231  TH2F *h_ele_EseedOPVsPhi;
232  TH2F *h_ele_EseedOPVsE;
233  TH1F *h_ele_EoPout;
234  TH1F *h_ele_EoPout_barrel;
235  TH1F *h_ele_EoPout_endcaps;
236  TH1F *h_ele_EoPout_eg;
237  TH1F *h_ele_EoPout_eg_barrel;
238  TH1F *h_ele_EoPout_eg_endcaps;
239  TH2F *h_ele_EoPoutVsEta;
240  TH2F *h_ele_EoPoutVsPhi;
241  TH2F *h_ele_EoPoutVsE;
242  TH1F *h_ele_EeleOPout;
243  TH1F *h_ele_EeleOPout_barrel;
244  TH1F *h_ele_EeleOPout_endcaps;
245  TH1F *h_ele_EeleOPout_eg;
246  TH1F *h_ele_EeleOPout_eg_barrel;
247  TH1F *h_ele_EeleOPout_eg_endcaps;
248  TH2F *h_ele_EeleOPoutVsEta;
249  TH2F *h_ele_EeleOPoutVsPhi;
250  TH2F *h_ele_EeleOPoutVsE;
251 
252  TH1F *h_ele_dEtaSc_propVtx;
253  TH1F *h_ele_dEtaSc_propVtx_barrel;
254  TH1F *h_ele_dEtaSc_propVtx_endcaps;
255  TH1F *h_ele_dEtaSc_propVtx_eg;
256  TH1F *h_ele_dEtaSc_propVtx_eg_barrel;
257  TH1F *h_ele_dEtaSc_propVtx_eg_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  TH1F *h_ele_dPhiSc_propVtx_eg;
265  TH1F *h_ele_dPhiSc_propVtx_eg_barrel;
266  TH1F *h_ele_dPhiSc_propVtx_eg_endcaps;
267  TH2F *h_ele_dPhiScVsEta_propVtx;
268  TH2F *h_ele_dPhiScVsPhi_propVtx;
269  TH2F *h_ele_dPhiScVsPt_propVtx;
270  TH1F *h_ele_dEtaCl_propOut;
271  TH1F *h_ele_dEtaCl_propOut_barrel;
272  TH1F *h_ele_dEtaCl_propOut_endcaps;
273  TH1F *h_ele_dEtaCl_propOut_eg;
274  TH1F *h_ele_dEtaCl_propOut_eg_barrel;
275  TH1F *h_ele_dEtaCl_propOut_eg_endcaps;
276  TH2F *h_ele_dEtaClVsEta_propOut;
277  TH2F *h_ele_dEtaClVsPhi_propOut;
278  TH2F *h_ele_dEtaClVsPt_propOut;
279  TH1F *h_ele_dPhiCl_propOut;
280  TH1F *h_ele_dPhiCl_propOut_barrel;
281  TH1F *h_ele_dPhiCl_propOut_endcaps;
282  TH1F *h_ele_dPhiCl_propOut_eg;
283  TH1F *h_ele_dPhiCl_propOut_eg_barrel;
284  TH1F *h_ele_dPhiCl_propOut_eg_endcaps;
285  TH2F *h_ele_dPhiClVsEta_propOut;
286  TH2F *h_ele_dPhiClVsPhi_propOut;
287  TH2F *h_ele_dPhiClVsPt_propOut;
288  TH1F *h_ele_dEtaEleCl_propOut;
289  TH1F *h_ele_dEtaEleCl_propOut_barrel;
290  TH1F *h_ele_dEtaEleCl_propOut_endcaps;
291  TH1F *h_ele_dEtaEleCl_propOut_eg;
292  TH1F *h_ele_dEtaEleCl_propOut_eg_barrel;
293  TH1F *h_ele_dEtaEleCl_propOut_eg_endcaps;
294  TH2F *h_ele_dEtaEleClVsEta_propOut;
295  TH2F *h_ele_dEtaEleClVsPhi_propOut;
296  TH2F *h_ele_dEtaEleClVsPt_propOut;
297  TH1F *h_ele_dPhiEleCl_propOut;
298  TH1F *h_ele_dPhiEleCl_propOut_barrel;
299  TH1F *h_ele_dPhiEleCl_propOut_endcaps;
300  TH1F *h_ele_dPhiEleCl_propOut_eg;
301  TH1F *h_ele_dPhiEleCl_propOut_eg_barrel;
302  TH1F *h_ele_dPhiEleCl_propOut_eg_endcaps;
303  TH2F *h_ele_dPhiEleClVsEta_propOut;
304  TH2F *h_ele_dPhiEleClVsPhi_propOut;
305  TH2F *h_ele_dPhiEleClVsPt_propOut;
306 
307  TH1F *h_ele_seed_dphi2_;
308  TH2F *h_ele_seed_dphi2VsEta_;
309  TH2F *h_ele_seed_dphi2VsPt_ ;
310  TH1F *h_ele_seed_drz2_;
311  TH2F *h_ele_seed_drz2VsEta_;
312  TH2F *h_ele_seed_drz2VsPt_;
313  TH1F *h_ele_seed_subdet2_;
314 
315  TH1F *h_ele_classes;
316  TH1F *h_ele_eta;
317  TH1F *h_ele_eta_golden;
318  TH1F *h_ele_eta_bbrem;
319  TH1F *h_ele_eta_narrow;
320  TH1F *h_ele_eta_shower;
321 
322  TH1F *h_ele_HoE;
323  TH1F *h_ele_HoE_barrel;
324  TH1F *h_ele_HoE_endcaps;
325  TH1F *h_ele_HoE_eg;
326  TH1F *h_ele_HoE_eg_barrel;
327  TH1F *h_ele_HoE_eg_endcaps;
328  TH1F *h_ele_HoE_fiducial;
329  TH2F *h_ele_HoEVsEta;
330  TH2F *h_ele_HoEVsPhi;
331  TH2F *h_ele_HoEVsE;
332 
333  TH1F *h_ele_fbrem;
334  TProfile *h_ele_fbremVsEta_mode;
335  TProfile *h_ele_fbremVsEta_mean;
336 
337  TH2F *h_ele_PinVsPoutGolden_mode;
338  TH2F *h_ele_PinVsPoutShowering_mode;
339  TH2F *h_ele_PinVsPoutGolden_mean;
340  TH2F *h_ele_PinVsPoutShowering_mean;
341  TH2F *h_ele_PtinVsPtoutGolden_mode;
342  TH2F *h_ele_PtinVsPtoutShowering_mode;
343  TH2F *h_ele_PtinVsPtoutGolden_mean;
344  TH2F *h_ele_PtinVsPtoutShowering_mean;
345  TH1F *histSclEoEmatchingObjectGolden_barrel;
346  TH1F *histSclEoEmatchingObjectGolden_endcaps;
347  TH1F *histSclEoEmatchingObjectShowering_barrel;
348  TH1F *histSclEoEmatchingObjectShowering_endcaps;
349 
350  TH1F *h_ele_mva;
351  TH1F *h_ele_provenance;
352 
353  TH1F *h_ele_tkSumPt_dr03;
354  TH1F *h_ele_ecalRecHitSumEt_dr03;
355  TH1F *h_ele_hcalDepth1TowerSumEt_dr03;
356  TH1F *h_ele_hcalDepth2TowerSumEt_dr03;
357  TH1F *h_ele_tkSumPt_dr04;
358  TH1F *h_ele_ecalRecHitSumEt_dr04;
359  TH1F *h_ele_hcalDepth1TowerSumEt_dr04;
360  TH1F *h_ele_hcalDepth2TowerSumEt_dr04;
361 
362  std::string outputFile_;
363  edm::InputTag electronCollection_;
364  edm::InputTag matchingObjectCollection_;
365  edm::InputTag beamSpot_;
366  std::string type_;
367  bool readAOD_;
368 
369  double maxPt_;
370  double maxAbsEta_;
371  double deltaR_;
372 
373  // histos limits and binning
374  double etamin;
375  double etamax;
376  double phimin;
377  double phimax;
378  double ptmax;
379  double pmax;
380  double eopmax;
381  double eopmaxsht;
382  double detamin;
383  double detamax;
384  double dphimin;
385  double dphimax;
386  double detamatchmin;
387  double detamatchmax;
388  double dphimatchmin;
389  double dphimatchmax;
390  double fhitsmax;
391  double lhitsmax;
392  double meemin;
393  double meemax;
394  double hoemin;
395  double hoemax;
396  int nbineta;
397  int nbinp;
398  int nbinpt;
399  int nbinpteff;
400  int nbinphi;
401  int nbinp2D;
402  int nbinpt2D;
403  int nbineta2D;
404  int nbinphi2D;
405  int nbineop;
406  int nbineop2D;
407  int nbinfhits;
408  int nbinlhits;
409  int nbinxyz;
410  int nbindeta;
411  int nbindphi;
412  int nbindetamatch;
413  int nbindphimatch;
414  int nbindetamatch2D;
415  int nbindphimatch2D;
416  int nbinmee;
417  int nbinhoe;
418 
419  };
420 
421 #endif
422 
423 
424 
GsfElectronFakeAnalyzer::transformer_
TrajectoryStateTransform transformer_
Definition: GsfElectronFakeAnalyzer.h:51
GsfElectronFakeAnalyzer::theMagField
edm::ESHandle< MagneticField > theMagField
Definition: GsfElectronFakeAnalyzer.h:53
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
GsfElectronFakeAnalyzer::GsfElectronFakeAnalyzer
GsfElectronFakeAnalyzer(const edm::ParameterSet &conf)
Definition: GsfElectronFakeAnalyzer.cc:51
GsfElectronFakeAnalyzer::analyze
void analyze(const edm::Event &e, const edm::EventSetup &c) override
Definition: GsfElectronFakeAnalyzer.cc:1231
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
GsfElectronFakeAnalyzer::pDD
edm::ESHandle< TrackerGeometry > pDD
Definition: GsfElectronFakeAnalyzer.h:52
Generated for CMSSW Reference Manual by doxygen 1.8.11