CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
GsfElectronDataAnalyzer.h
Go to the documentation of this file.
1 #ifndef GsfElectronDataAnalyzer_h
2 #define GsfElectronDataAnalyzer_h
3 
4 //
5 // Package: RecoEgamma/Examples
6 // Class: GsfElectronDataAnalyzer
7 //
8 
9 //
10 // Original Author: Ursula Berthon, Claude Charlot
11 // Created: Mon Mar 27 13:22:06 CEST 2006
12 // $Id: GsfElectronDataAnalyzer.h,v 1.19 2010/02/20 21:00:54 wmtan Exp $
13 //
14 //
15 
16 
17 #include "FWCore/Framework/interface/EDAnalyzer.h"
18 #include "FWCore/Framework/interface/Event.h"
19 #include "DataFormats/Common/interface/Handle.h"
20 #include "FWCore/Framework/interface/ESHandle.h"
21 #include "FWCore/Framework/interface/EventSetup.h"
22 #include "FWCore/ParameterSet/interface/ParameterSet.h"
23 
24 #include "DataFormats/Common/interface/EDProduct.h"
25 
26 #include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
27 
28 #include "TrackingTools/TrajectoryState/interface/TrajectoryStateTransform.h"
29 
30 class MagneticField;
31 class TFile;
32 class TH1F;
33 class TH2F;
34 class TH1I;
35 class TProfile;
36 class TTree;
37 
38 class GsfElectronDataAnalyzer : public edm::EDAnalyzer
39 {
40  public:
41 
42  explicit GsfElectronDataAnalyzer(const edm::ParameterSet& conf);
43 
44  virtual ~GsfElectronDataAnalyzer();
45 
46  virtual void beginJob();
47  virtual void endJob();
48  virtual void analyze(const edm::Event& e, const edm::EventSetup& c);
49 
50  private:
51 
52  bool trigger(const edm::Event & e);
53 
54  unsigned int nEvents_;
55  unsigned int nAfterTrigger_;
56 
57  edm::InputTag triggerResults_;
58 
59  std::vector<std::string > HLTPathsByName_;
60  std::vector<unsigned int> HLTPathsByIndex_;
61 
62  std::string outputFile_;
63  edm::InputTag electronCollection_;
64  edm::InputTag matchingObjectCollection_;
65  std::string matchingCondition_;
66  std::string type_;
67  bool readAOD_;
68  // matching
69  double maxPtMatchingObject_;
70  double maxAbsEtaMatchingObject_;
71  double deltaR_;
72 
73  TrajectoryStateTransform transformer_;
74  edm::ESHandle<TrackerGeometry> pDD;
75  edm::ESHandle<MagneticField> theMagField;
76 
77  TFile *histfile_;
78  TTree *tree_;
79  float mcEnergy[10], mcEta[10], mcPhi[10], mcPt[10], mcQ[10];
80  float superclusterEnergy[10], superclusterEta[10], superclusterPhi[10], superclusterEt[10];
81  float seedMomentum[10], seedEta[10], seedPhi[10], seedPt[10], seedQ[10];
82 
83  TH1F *h_matchingObjectNum;
84 
85  TH1F *h_matchingObjectEta;
86  TH1F *h_matchingObjectAbsEta;
87  TH1F *h_matchingObjectP;
88  TH1F *h_matchingObjectPt;
89  TH1F *h_matchingObjectPhi;
90  TH1F *h_matchingObjectZ;
91 
92  TH1F *h_ele_matchingObjectEta_matched;
93  TH1F *h_ele_matchingObjectAbsEta_matched;
94  TH1F *h_ele_matchingObjectPt_matched;
95  TH1F *h_ele_matchingObjectPhi_matched;
96  TH1F *h_ele_matchingObjectZ_matched;
97 
98  TH1F *h_ele_mee_all;
99  TH1F *h_ele_mee_os;
100  TH1F *h_ele_mee_os_ebeb;
101  TH1F *h_ele_mee_os_ebee;
102  TH1F *h_ele_mee_os_eeee;
103  TH1F *h_ele_mee_os_gg;
104  TH1F *h_ele_mee_os_gb;
105  TH1F *h_ele_mee_os_bb;
106 
107  TH2F *h_ele_E2mnE1vsMee_all;
108  TH2F *h_ele_E2mnE1vsMee_egeg_all;
109 
110  TH1F *h_ele_charge;
111  TH2F *h_ele_chargeVsEta;
112  TH2F *h_ele_chargeVsPhi;
113  TH2F *h_ele_chargeVsPt;
114  TH1F *h_ele_vertexP;
115  TH1F *h_ele_vertexPt;
116  TH1F *h_ele_Et;
117  TH2F *h_ele_vertexPtVsEta;
118  TH2F *h_ele_vertexPtVsPhi;
119  TH1F *h_ele_vertexPt_5100;
120  TH1F *h_ele_vertexEta;
121  TH2F *h_ele_vertexEtaVsPhi;
122  TH1F *h_ele_vertexAbsEta;
123  TH1F *h_ele_vertexPhi;
124  TH1F *h_ele_vertexX;
125  TH1F *h_ele_vertexY;
126  TH1F *h_ele_vertexZ;
127  TH1F *h_ele_vertexTIP;
128  TH2F *h_ele_vertexTIPVsEta;
129  TH2F *h_ele_vertexTIPVsPhi;
130  TH2F *h_ele_vertexTIPVsPt;
131 
132  TH1F *histNum_;
133 
134  TH1F *histSclEn_ ;
135  TH1F *histSclEoEmatchingObject_barrel_matched;
136  TH1F *histSclEoEmatchingObject_endcaps_matched;
137  TH1F *histSclEoEmatchingObject_barrel_new_matched;
138  TH1F *histSclEoEmatchingObject_endcaps_new_matched;
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_PtoPtmatchingObject_matched;
185  TH1F *h_ele_PoPmatchingObject_matched;
186  TH2F *h_ele_PoPmatchingObjectVsEta_matched;
187  TH2F *h_ele_PoPmatchingObjectVsPhi_matched;
188  TH2F *h_ele_PoPmatchingObjectVsPt_matched;
189  TH1F *h_ele_PoPmatchingObject_barrel_matched;
190  TH1F *h_ele_PoPmatchingObject_endcaps_matched;
191  TH1F *h_ele_PtoPtmatchingObject_barrel_matched;
192  TH1F *h_ele_PtoPtmatchingObject_endcaps_matched;
193  TH1F *h_ele_EtaMnEtamatchingObject_matched;
194  TH2F *h_ele_EtaMnEtamatchingObjectVsEta_matched;
195  TH2F *h_ele_EtaMnEtamatchingObjectVsPhi_matched;
196  TH2F *h_ele_EtaMnEtamatchingObjectVsPt_matched;
197  TH1F *h_ele_PhiMnPhimatchingObject_matched;
198  TH1F *h_ele_PhiMnPhimatchingObject2_matched;
199  TH2F *h_ele_PhiMnPhimatchingObjectVsEta_matched;
200  TH2F *h_ele_PhiMnPhimatchingObjectVsPhi_matched;
201  TH2F *h_ele_PhiMnPhimatchingObjectVsPt_matched;
202  TH1F *h_ele_PinMnPout;
203  TH1F *h_ele_PinMnPout_mode;
204  TH2F *h_ele_PinMnPoutVsEta_mode;
205  TH2F *h_ele_PinMnPoutVsPhi_mode;
206  TH2F *h_ele_PinMnPoutVsPt_mode;
207  TH2F *h_ele_PinMnPoutVsE_mode;
208  TH2F *h_ele_PinMnPoutVsChi2_mode;
209 
210  TH1F *h_ele_outerP;
211  TH1F *h_ele_outerP_mode;
212  TH2F *h_ele_outerPVsEta_mode;
213  TH1F *h_ele_outerPt;
214  TH1F *h_ele_outerPt_mode;
215  TH2F *h_ele_outerPtVsEta_mode;
216  TH2F *h_ele_outerPtVsPhi_mode;
217  TH2F *h_ele_outerPtVsPt_mode;
218  TH1F *h_ele_EoP;
219  TH1F *h_ele_EoP_barrel;
220  TH1F *h_ele_EoP_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  TH2F *h_ele_EseedOPVsEta;
228  TH2F *h_ele_EseedOPVsPhi;
229  TH2F *h_ele_EseedOPVsE;
230  TH1F *h_ele_EoPout;
231  TH1F *h_ele_EoPout_barrel;
232  TH1F *h_ele_EoPout_endcaps;
233  TH2F *h_ele_EoPoutVsEta;
234  TH2F *h_ele_EoPoutVsPhi;
235  TH2F *h_ele_EoPoutVsE;
236  TH1F *h_ele_EeleOPout;
237  TH1F *h_ele_EeleOPout_barrel;
238  TH1F *h_ele_EeleOPout_endcaps;
239  TH2F *h_ele_EeleOPoutVsEta;
240  TH2F *h_ele_EeleOPoutVsPhi;
241  TH2F *h_ele_EeleOPoutVsE;
242 
243  TH1F *h_ele_dEtaSc_propVtx;
244  TH1F *h_ele_dEtaSc_propVtx_barrel;
245  TH1F *h_ele_dEtaSc_propVtx_endcaps;
246  TH2F *h_ele_dEtaScVsEta_propVtx;
247  TH2F *h_ele_dEtaScVsPhi_propVtx;
248  TH2F *h_ele_dEtaScVsPt_propVtx;
249  TH1F *h_ele_dPhiSc_propVtx;
250  TH1F *h_ele_dPhiSc_propVtx_barrel;
251  TH1F *h_ele_dPhiSc_propVtx_endcaps;
252  TH2F *h_ele_dPhiScVsEta_propVtx;
253  TH2F *h_ele_dPhiScVsPhi_propVtx;
254  TH2F *h_ele_dPhiScVsPt_propVtx;
255  TH1F *h_ele_dEtaCl_propOut;
256  TH1F *h_ele_dEtaCl_propOut_barrel;
257  TH1F *h_ele_dEtaCl_propOut_endcaps;
258  TH2F *h_ele_dEtaClVsEta_propOut;
259  TH2F *h_ele_dEtaClVsPhi_propOut;
260  TH2F *h_ele_dEtaClVsPt_propOut;
261  TH1F *h_ele_dPhiCl_propOut;
262  TH1F *h_ele_dPhiCl_propOut_barrel;
263  TH1F *h_ele_dPhiCl_propOut_endcaps;
264  TH2F *h_ele_dPhiClVsEta_propOut;
265  TH2F *h_ele_dPhiClVsPhi_propOut;
266  TH2F *h_ele_dPhiClVsPt_propOut;
267  TH1F *h_ele_dEtaEleCl_propOut;
268  TH1F *h_ele_dEtaEleCl_propOut_barrel;
269  TH1F *h_ele_dEtaEleCl_propOut_endcaps;
270  TH2F *h_ele_dEtaEleClVsEta_propOut;
271  TH2F *h_ele_dEtaEleClVsPhi_propOut;
272  TH2F *h_ele_dEtaEleClVsPt_propOut;
273  TH1F *h_ele_dPhiEleCl_propOut;
274  TH1F *h_ele_dPhiEleCl_propOut_barrel;
275  TH1F *h_ele_dPhiEleCl_propOut_endcaps;
276  TH2F *h_ele_dPhiEleClVsEta_propOut;
277  TH2F *h_ele_dPhiEleClVsPhi_propOut;
278  TH2F *h_ele_dPhiEleClVsPt_propOut;
279 
280  TH1F *h_ele_seed_dphi2_;
281  TH2F *h_ele_seed_dphi2VsEta_;
282  TH2F *h_ele_seed_dphi2VsPt_ ;
283  TH1F *h_ele_seed_drz2_;
284  TH2F *h_ele_seed_drz2VsEta_;
285  TH2F *h_ele_seed_drz2VsPt_;
286  TH1F *h_ele_seed_subdet2_;
287 
288  TH1F *h_ele_classes;
289  TH1F *h_ele_eta;
290  TH1F *h_ele_eta_golden;
291  TH1F *h_ele_eta_bbrem;
292  TH1F *h_ele_eta_narrow;
293  TH1F *h_ele_eta_shower;
294 
295  TH1F *h_ele_HoE;
296  TH1F *h_ele_HoE_barrel;
297  TH1F *h_ele_HoE_endcaps;
298  TH1F *h_ele_HoE_fiducial;
299  TH2F *h_ele_HoEVsEta;
300  TH2F *h_ele_HoEVsPhi;
301  TH2F *h_ele_HoEVsE;
302 
303  TH1F *h_ele_fbrem;
304  TProfile *h_ele_fbremVsEta_mode;
305  TProfile *h_ele_fbremVsEta_mean;
306 
307  TH2F *h_ele_PinVsPoutGolden_mode;
308  TH2F *h_ele_PinVsPoutShowering_mode;
309  TH2F *h_ele_PinVsPoutGolden_mean;
310  TH2F *h_ele_PinVsPoutShowering_mean;
311  TH2F *h_ele_PtinVsPtoutGolden_mode;
312  TH2F *h_ele_PtinVsPtoutShowering_mode;
313  TH2F *h_ele_PtinVsPtoutGolden_mean;
314  TH2F *h_ele_PtinVsPtoutShowering_mean;
315  TH1F *histSclEoEmatchingObjectGolden_barrel;
316  TH1F *histSclEoEmatchingObjectGolden_endcaps;
317  TH1F *histSclEoEmatchingObjectShowering_barrel;
318  TH1F *histSclEoEmatchingObjectShowering_endcaps;
319 
320  TH1F *h_ele_mva;
321  TH1F *h_ele_provenance;
322 
323  TH1F *h_ele_tkSumPt_dr03;
324  TH1F *h_ele_ecalRecHitSumEt_dr03;
325  TH1F *h_ele_hcalDepth1TowerSumEt_dr03;
326  TH1F *h_ele_hcalDepth2TowerSumEt_dr03;
327  TH1F *h_ele_tkSumPt_dr04;
328  TH1F *h_ele_ecalRecHitSumEt_dr04;
329  TH1F *h_ele_hcalDepth1TowerSumEt_dr04;
330  TH1F *h_ele_hcalDepth2TowerSumEt_dr04;
331 
332  // electron selection
333  double minEt_;
334  double minPt_;
335  double maxAbsEta_;
336  bool isEB_;
337  bool isEE_;
338  bool isNotEBEEGap_;
339  bool isEcalDriven_;
340  bool isTrackerDriven_;
341  double eOverPMinBarrel_;
342  double eOverPMaxBarrel_;
343  double eOverPMinEndcaps_;
344  double eOverPMaxEndcaps_;
345  double dEtaMinBarrel_;
346  double dEtaMaxBarrel_;
347  double dEtaMinEndcaps_;
348  double dEtaMaxEndcaps_;
349  double dPhiMinBarrel_;
350  double dPhiMaxBarrel_;
351  double dPhiMinEndcaps_;
352  double dPhiMaxEndcaps_;
353  double sigIetaIetaMinBarrel_;
354  double sigIetaIetaMaxBarrel_;
355  double sigIetaIetaMinEndcaps_;
356  double sigIetaIetaMaxEndcaps_;
357  double hadronicOverEmMaxBarrel_;
358  double hadronicOverEmMaxEndcaps_;
359  double mvaMin_;
360  double tipMaxBarrel_;
361  double tipMaxEndcaps_;
362  double tkIso03Max_;
363  double hcalIso03Depth1MaxBarrel_;
364  double hcalIso03Depth1MaxEndcaps_;
365  double hcalIso03Depth2MaxEndcaps_;
366  double ecalIso03MaxBarrel_;
367  double ecalIso03MaxEndcaps_;
368 
369  // histos limits and binning
370  double etamin;
371  double etamax;
372  double phimin;
373  double phimax;
374  double ptmax;
375  double pmax;
376  double eopmax;
377  double eopmaxsht;
378  double detamin;
379  double detamax;
380  double dphimin;
381  double dphimax;
382  double detamatchmin;
383  double detamatchmax;
384  double dphimatchmin;
385  double dphimatchmax;
386  double fhitsmax;
387  double lhitsmax;
388  double poptruemin;
389  double poptruemax;
390  double meemin;
391  double meemax;
392  double hoemin;
393  double hoemax;
394  int nbineta;
395  int nbinp;
396  int nbinpt;
397  int nbinpteff;
398  int nbinphi;
399  int nbinp2D;
400  int nbinpt2D;
401  int nbineta2D;
402  int nbinphi2D;
403  int nbineop;
404  int nbineop2D;
405  int nbinfhits;
406  int nbinlhits;
407  int nbinxyz;
408  int nbindeta;
409  int nbindphi;
410  int nbindetamatch;
411  int nbindphimatch;
412  int nbindetamatch2D;
413  int nbindphimatch2D;
414  int nbinpoptrue;
415  int nbinmee;
416  int nbinhoe;
417 
418  };
419 
420 #endif
421 
422 
423 
GsfElectronDataAnalyzer::HLTPathsByIndex_
std::vector< unsigned int > HLTPathsByIndex_
Definition: GsfElectronDataAnalyzer.h:60
GsfElectronDataAnalyzer::pDD
edm::ESHandle< TrackerGeometry > pDD
Definition: GsfElectronDataAnalyzer.h:74
GsfElectronDataAnalyzer::transformer_
TrajectoryStateTransform transformer_
Definition: GsfElectronDataAnalyzer.h:73
GsfElectronDataAnalyzer::GsfElectronDataAnalyzer
GsfElectronDataAnalyzer(const edm::ParameterSet &conf)
Definition: GsfElectronDataAnalyzer.cc:58
GsfElectronDataAnalyzer::trigger
bool trigger(const edm::Event &e)
Definition: GsfElectronDataAnalyzer.cc:1472
GsfElectronDataAnalyzer::theMagField
edm::ESHandle< MagneticField > theMagField
Definition: GsfElectronDataAnalyzer.h:75
dbtoconf.conf
tuple conf
Definition: dbtoconf.py:185
GsfElectronDataAnalyzer::analyze
virtual void analyze(const edm::Event &e, const edm::EventSetup &c)
Definition: GsfElectronDataAnalyzer.cc:1054
GsfElectronDataAnalyzer::HLTPathsByName_
std::vector< std::string > HLTPathsByName_
Definition: GsfElectronDataAnalyzer.h:59
Generated for CMSSW Reference Manual by doxygen 1.8.5