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