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