CMS 3D CMS Logo
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 //
13 //
14 
15 
16 #include "FWCore/Framework/interface/one/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 GsfElectronDataAnalyzer : public edm::one::EDAnalyzer<>
37 {
38  public:
39 
40  explicit GsfElectronDataAnalyzer(const edm::ParameterSet& conf);
41 
42  ~GsfElectronDataAnalyzer() override;
43 
44  void beginJob() override;
45  void endJob() override;
46  void analyze(const edm::Event& e, const edm::EventSetup& c) override;
47 
48  private:
49 
50  bool trigger(const edm::Event & e);
51 
52  unsigned int nEvents_;
53  unsigned int nAfterTrigger_;
54 
55  edm::InputTag triggerResults_;
56  edm::InputTag beamSpot_;
57 
58  std::vector<std::string > HLTPathsByName_;
59  std::vector<unsigned int> HLTPathsByIndex_;
60 
61  std::string outputFile_;
62  edm::InputTag electronCollection_;
63  edm::InputTag matchingObjectCollection_;
64  std::string matchingCondition_;
65  std::string type_;
66  bool readAOD_;
67  // matching
68  double maxPtMatchingObject_;
69  double maxAbsEtaMatchingObject_;
70  double deltaR_;
71 
72  TrajectoryStateTransform transformer_;
73  edm::ESHandle<TrackerGeometry> pDD;
74  edm::ESHandle<MagneticField> theMagField;
75 
76  TFile *histfile_;
77  TTree *tree_;
78  float mcEnergy[10], mcEta[10], mcPhi[10], mcPt[10], mcQ[10];
79  float superclusterEnergy[10], superclusterEta[10], superclusterPhi[10], superclusterEt[10];
80  float seedMomentum[10], seedEta[10], seedPhi[10], seedPt[10], seedQ[10];
81 
82  TH1F *h_matchingObjectNum;
83 
84  TH1F *h_matchingObjectEta;
85  TH1F *h_matchingObjectAbsEta;
86  TH1F *h_matchingObjectP;
87  TH1F *h_matchingObjectPt;
88  TH1F *h_matchingObjectPhi;
89  TH1F *h_matchingObjectZ;
90 
91  TH1F *h_ele_matchingObjectEta_matched;
92  TH1F *h_ele_matchingObjectAbsEta_matched;
93  TH1F *h_ele_matchingObjectPt_matched;
94  TH1F *h_ele_matchingObjectPhi_matched;
95  TH1F *h_ele_matchingObjectZ_matched;
96 
97  TH1F *h_ele_mee_all;
98  TH1F *h_ele_mee_os;
99  TH1F *h_ele_mee_os_ebeb;
100  TH1F *h_ele_mee_os_ebee;
101  TH1F *h_ele_mee_os_eeee;
102  TH1F *h_ele_mee_os_gg;
103  TH1F *h_ele_mee_os_gb;
104  TH1F *h_ele_mee_os_bb;
105 
106  TH2F *h_ele_E2mnE1vsMee_all;
107  TH2F *h_ele_E2mnE1vsMee_egeg_all;
108 
109  TH1F *h_ele_charge;
110  TH2F *h_ele_chargeVsEta;
111  TH2F *h_ele_chargeVsPhi;
112  TH2F *h_ele_chargeVsPt;
113  TH1F *h_ele_vertexP;
114  TH1F *h_ele_vertexPt;
115  TH1F *h_ele_Et;
116  TH2F *h_ele_vertexPtVsEta;
117  TH2F *h_ele_vertexPtVsPhi;
118  TH1F *h_ele_vertexPt_5100;
119  TH1F *h_ele_vertexEta;
120  TH2F *h_ele_vertexEtaVsPhi;
121  TH1F *h_ele_vertexAbsEta;
122  TH1F *h_ele_vertexPhi;
123  TH1F *h_ele_vertexX;
124  TH1F *h_ele_vertexY;
125  TH1F *h_ele_vertexZ;
126  TH1F *h_ele_vertexTIP;
127  TH2F *h_ele_vertexTIPVsEta;
128  TH2F *h_ele_vertexTIPVsPhi;
129  TH2F *h_ele_vertexTIPVsPt;
130 
131  TH1F *histNum_;
132 
133  TH1F *histSclEn_ ;
134  TH1F *histSclEoEmatchingObject_barrel_matched;
135  TH1F *histSclEoEmatchingObject_endcaps_matched;
136  TH1F *histSclEoEmatchingObject_barrel_new_matched;
137  TH1F *histSclEoEmatchingObject_endcaps_new_matched;
138  TH1F *histSclEt_ ;
139  TH2F *histSclEtVsEta_ ;
140  TH2F *histSclEtVsPhi_ ;
141  TH2F *histSclEtaVsPhi_ ;
142  TH1F *histSclEta_ ;
143  TH1F *histSclPhi_ ;
144  TH1F *histSclSigEtaEta_ ;
145  TH1F *histSclSigEtaEta_barrel_ ;
146  TH1F *histSclSigEtaEta_endcaps_ ;
147  TH1F *histSclSigIEtaIEta_ ;
148  TH1F *histSclSigIEtaIEta_barrel_ ;
149  TH1F *histSclSigIEtaIEta_endcaps_ ;
150  TH1F *histSclE1x5_ ;
151  TH1F *histSclE1x5_barrel_ ;
152  TH1F *histSclE1x5_endcaps_ ;
153  TH1F *histSclE2x5max_ ;
154  TH1F *histSclE2x5max_barrel_ ;
155  TH1F *histSclE2x5max_endcaps_ ;
156  TH1F *histSclE5x5_ ;
157  TH1F *histSclE5x5_barrel_ ;
158  TH1F *histSclE5x5_endcaps_ ;
159 
160  TH1F *h_ele_ambiguousTracks;
161  TH2F *h_ele_ambiguousTracksVsEta;
162  TH2F *h_ele_ambiguousTracksVsPhi;
163  TH2F *h_ele_ambiguousTracksVsPt;
164  TH1F *h_ele_foundHits;
165  TH1F *h_ele_foundHits_barrel;
166  TH1F *h_ele_foundHits_endcaps;
167  TH2F *h_ele_foundHitsVsEta;
168  TH2F *h_ele_foundHitsVsPhi;
169  TH2F *h_ele_foundHitsVsPt;
170  TH1F *h_ele_lostHits;
171  TH1F *h_ele_lostHits_barrel;
172  TH1F *h_ele_lostHits_endcaps;
173  TH2F *h_ele_lostHitsVsEta;
174  TH2F *h_ele_lostHitsVsPhi;
175  TH2F *h_ele_lostHitsVsPt;
176  TH1F *h_ele_chi2;
177  TH1F *h_ele_chi2_barrel_;
178  TH1F *h_ele_chi2_endcaps_;
179  TH2F *h_ele_chi2VsEta;
180  TH2F *h_ele_chi2VsPhi;
181  TH2F *h_ele_chi2VsPt;
182 
183  TH1F *h_ele_PtoPtmatchingObject_matched;
184  TH1F *h_ele_PoPmatchingObject_matched;
185  TH2F *h_ele_PoPmatchingObjectVsEta_matched;
186  TH2F *h_ele_PoPmatchingObjectVsPhi_matched;
187  TH2F *h_ele_PoPmatchingObjectVsPt_matched;
188  TH1F *h_ele_PoPmatchingObject_barrel_matched;
189  TH1F *h_ele_PoPmatchingObject_endcaps_matched;
190  TH1F *h_ele_PtoPtmatchingObject_barrel_matched;
191  TH1F *h_ele_PtoPtmatchingObject_endcaps_matched;
192  TH1F *h_ele_EtaMnEtamatchingObject_matched;
193  TH2F *h_ele_EtaMnEtamatchingObjectVsEta_matched;
194  TH2F *h_ele_EtaMnEtamatchingObjectVsPhi_matched;
195  TH2F *h_ele_EtaMnEtamatchingObjectVsPt_matched;
196  TH1F *h_ele_PhiMnPhimatchingObject_matched;
197  TH1F *h_ele_PhiMnPhimatchingObject2_matched;
198  TH2F *h_ele_PhiMnPhimatchingObjectVsEta_matched;
199  TH2F *h_ele_PhiMnPhimatchingObjectVsPhi_matched;
200  TH2F *h_ele_PhiMnPhimatchingObjectVsPt_matched;
201  TH1F *h_ele_PinMnPout;
202  TH1F *h_ele_PinMnPout_mode;
203  TH2F *h_ele_PinMnPoutVsEta_mode;
204  TH2F *h_ele_PinMnPoutVsPhi_mode;
205  TH2F *h_ele_PinMnPoutVsPt_mode;
206  TH2F *h_ele_PinMnPoutVsE_mode;
207  TH2F *h_ele_PinMnPoutVsChi2_mode;
208 
209  TH1F *h_ele_outerP;
210  TH1F *h_ele_outerP_mode;
211  TH2F *h_ele_outerPVsEta_mode;
212  TH1F *h_ele_outerPt;
213  TH1F *h_ele_outerPt_mode;
214  TH2F *h_ele_outerPtVsEta_mode;
215  TH2F *h_ele_outerPtVsPhi_mode;
216  TH2F *h_ele_outerPtVsPt_mode;
217  TH1F *h_ele_EoP;
218  TH1F *h_ele_EoP_barrel;
219  TH1F *h_ele_EoP_endcaps;
220  TH2F *h_ele_EoPVsEta;
221  TH2F *h_ele_EoPVsPhi;
222  TH2F *h_ele_EoPVsE;
223  TH1F *h_ele_EseedOP;
224  TH1F *h_ele_EseedOP_barrel;
225  TH1F *h_ele_EseedOP_endcaps;
226  TH2F *h_ele_EseedOPVsEta;
227  TH2F *h_ele_EseedOPVsPhi;
228  TH2F *h_ele_EseedOPVsE;
229  TH1F *h_ele_EoPout;
230  TH1F *h_ele_EoPout_barrel;
231  TH1F *h_ele_EoPout_endcaps;
232  TH2F *h_ele_EoPoutVsEta;
233  TH2F *h_ele_EoPoutVsPhi;
234  TH2F *h_ele_EoPoutVsE;
235  TH1F *h_ele_EeleOPout;
236  TH1F *h_ele_EeleOPout_barrel;
237  TH1F *h_ele_EeleOPout_endcaps;
238  TH2F *h_ele_EeleOPoutVsEta;
239  TH2F *h_ele_EeleOPoutVsPhi;
240  TH2F *h_ele_EeleOPoutVsE;
241 
242  TH1F *h_ele_dEtaSc_propVtx;
243  TH1F *h_ele_dEtaSc_propVtx_barrel;
244  TH1F *h_ele_dEtaSc_propVtx_endcaps;
245  TH2F *h_ele_dEtaScVsEta_propVtx;
246  TH2F *h_ele_dEtaScVsPhi_propVtx;
247  TH2F *h_ele_dEtaScVsPt_propVtx;
248  TH1F *h_ele_dPhiSc_propVtx;
249  TH1F *h_ele_dPhiSc_propVtx_barrel;
250  TH1F *h_ele_dPhiSc_propVtx_endcaps;
251  TH2F *h_ele_dPhiScVsEta_propVtx;
252  TH2F *h_ele_dPhiScVsPhi_propVtx;
253  TH2F *h_ele_dPhiScVsPt_propVtx;
254  TH1F *h_ele_dEtaCl_propOut;
255  TH1F *h_ele_dEtaCl_propOut_barrel;
256  TH1F *h_ele_dEtaCl_propOut_endcaps;
257  TH2F *h_ele_dEtaClVsEta_propOut;
258  TH2F *h_ele_dEtaClVsPhi_propOut;
259  TH2F *h_ele_dEtaClVsPt_propOut;
260  TH1F *h_ele_dPhiCl_propOut;
261  TH1F *h_ele_dPhiCl_propOut_barrel;
262  TH1F *h_ele_dPhiCl_propOut_endcaps;
263  TH2F *h_ele_dPhiClVsEta_propOut;
264  TH2F *h_ele_dPhiClVsPhi_propOut;
265  TH2F *h_ele_dPhiClVsPt_propOut;
266  TH1F *h_ele_dEtaEleCl_propOut;
267  TH1F *h_ele_dEtaEleCl_propOut_barrel;
268  TH1F *h_ele_dEtaEleCl_propOut_endcaps;
269  TH2F *h_ele_dEtaEleClVsEta_propOut;
270  TH2F *h_ele_dEtaEleClVsPhi_propOut;
271  TH2F *h_ele_dEtaEleClVsPt_propOut;
272  TH1F *h_ele_dPhiEleCl_propOut;
273  TH1F *h_ele_dPhiEleCl_propOut_barrel;
274  TH1F *h_ele_dPhiEleCl_propOut_endcaps;
275  TH2F *h_ele_dPhiEleClVsEta_propOut;
276  TH2F *h_ele_dPhiEleClVsPhi_propOut;
277  TH2F *h_ele_dPhiEleClVsPt_propOut;
278 
279  TH1F *h_ele_seed_dphi2_;
280  TH2F *h_ele_seed_dphi2VsEta_;
281  TH2F *h_ele_seed_dphi2VsPt_ ;
282  TH1F *h_ele_seed_drz2_;
283  TH2F *h_ele_seed_drz2VsEta_;
284  TH2F *h_ele_seed_drz2VsPt_;
285  TH1F *h_ele_seed_subdet2_;
286 
287  TH1F *h_ele_classes;
288  TH1F *h_ele_eta;
289  TH1F *h_ele_eta_golden;
290  TH1F *h_ele_eta_bbrem;
291  TH1F *h_ele_eta_narrow;
292  TH1F *h_ele_eta_shower;
293 
294  TH1F *h_ele_HoE;
295  TH1F *h_ele_HoE_barrel;
296  TH1F *h_ele_HoE_endcaps;
297  TH1F *h_ele_HoE_fiducial;
298  TH2F *h_ele_HoEVsEta;
299  TH2F *h_ele_HoEVsPhi;
300  TH2F *h_ele_HoEVsE;
301 
302  TH1F *h_ele_fbrem;
303  TProfile *h_ele_fbremVsEta_mode;
304  TProfile *h_ele_fbremVsEta_mean;
305 
306  TH2F *h_ele_PinVsPoutGolden_mode;
307  TH2F *h_ele_PinVsPoutShowering_mode;
308  TH2F *h_ele_PinVsPoutGolden_mean;
309  TH2F *h_ele_PinVsPoutShowering_mean;
310  TH2F *h_ele_PtinVsPtoutGolden_mode;
311  TH2F *h_ele_PtinVsPtoutShowering_mode;
312  TH2F *h_ele_PtinVsPtoutGolden_mean;
313  TH2F *h_ele_PtinVsPtoutShowering_mean;
314  TH1F *histSclEoEmatchingObjectGolden_barrel;
315  TH1F *histSclEoEmatchingObjectGolden_endcaps;
316  TH1F *histSclEoEmatchingObjectShowering_barrel;
317  TH1F *histSclEoEmatchingObjectShowering_endcaps;
318 
319  TH1F *h_ele_mva;
320  TH1F *h_ele_provenance;
321 
322  TH1F *h_ele_tkSumPt_dr03;
323  TH1F *h_ele_ecalRecHitSumEt_dr03;
324  TH1F *h_ele_hcalDepth1TowerSumEt_dr03;
325  TH1F *h_ele_hcalDepth2TowerSumEt_dr03;
326  TH1F *h_ele_tkSumPt_dr04;
327  TH1F *h_ele_ecalRecHitSumEt_dr04;
328  TH1F *h_ele_hcalDepth1TowerSumEt_dr04;
329  TH1F *h_ele_hcalDepth2TowerSumEt_dr04;
330 
331  // electron selection
332  double minEt_;
333  double minPt_;
334  double maxAbsEta_;
335  bool isEB_;
336  bool isEE_;
337  bool isNotEBEEGap_;
338  bool isEcalDriven_;
339  bool isTrackerDriven_;
340  double eOverPMinBarrel_;
341  double eOverPMaxBarrel_;
342  double eOverPMinEndcaps_;
343  double eOverPMaxEndcaps_;
344  double dEtaMinBarrel_;
345  double dEtaMaxBarrel_;
346  double dEtaMinEndcaps_;
347  double dEtaMaxEndcaps_;
348  double dPhiMinBarrel_;
349  double dPhiMaxBarrel_;
350  double dPhiMinEndcaps_;
351  double dPhiMaxEndcaps_;
352  double sigIetaIetaMinBarrel_;
353  double sigIetaIetaMaxBarrel_;
354  double sigIetaIetaMinEndcaps_;
355  double sigIetaIetaMaxEndcaps_;
356  double hadronicOverEmMaxBarrel_;
357  double hadronicOverEmMaxEndcaps_;
358  double mvaMin_;
359  double tipMaxBarrel_;
360  double tipMaxEndcaps_;
361  double tkIso03Max_;
362  double hcalIso03Depth1MaxBarrel_;
363  double hcalIso03Depth1MaxEndcaps_;
364  double hcalIso03Depth2MaxEndcaps_;
365  double ecalIso03MaxBarrel_;
366  double ecalIso03MaxEndcaps_;
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 poptruemin;
388  double poptruemax;
389  double meemin;
390  double meemax;
391  double hoemin;
392  double hoemax;
393  int nbineta;
394  int nbinp;
395  int nbinpt;
396  int nbinpteff;
397  int nbinphi;
398  int nbinp2D;
399  int nbinpt2D;
400  int nbineta2D;
401  int nbinphi2D;
402  int nbineop;
403  int nbineop2D;
404  int nbinfhits;
405  int nbinlhits;
406  int nbinxyz;
407  int nbindeta;
408  int nbindphi;
409  int nbindetamatch;
410  int nbindphimatch;
411  int nbindetamatch2D;
412  int nbindphimatch2D;
413  int nbinpoptrue;
414  int nbinmee;
415  int nbinhoe;
416 
417  };
418 
419 #endif
420 
421 
422 
GsfElectronDataAnalyzer::HLTPathsByIndex_
std::vector< unsigned int > HLTPathsByIndex_
Definition: GsfElectronDataAnalyzer.h:59
GsfElectronDataAnalyzer::analyze
void analyze(const edm::Event &e, const edm::EventSetup &c) override
Definition: GsfElectronDataAnalyzer.cc:1052
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
GsfElectronDataAnalyzer::pDD
edm::ESHandle< TrackerGeometry > pDD
Definition: GsfElectronDataAnalyzer.h:73
GsfElectronDataAnalyzer::transformer_
TrajectoryStateTransform transformer_
Definition: GsfElectronDataAnalyzer.h:72
GsfElectronDataAnalyzer::GsfElectronDataAnalyzer
GsfElectronDataAnalyzer(const edm::ParameterSet &conf)
Definition: GsfElectronDataAnalyzer.cc:55
GsfElectronDataAnalyzer::trigger
bool trigger(const edm::Event &e)
Definition: GsfElectronDataAnalyzer.cc:1470
GsfElectronDataAnalyzer::theMagField
edm::ESHandle< MagneticField > theMagField
Definition: GsfElectronDataAnalyzer.h:74
GsfElectronDataAnalyzer::HLTPathsByName_
std::vector< std::string > HLTPathsByName_
Definition: GsfElectronDataAnalyzer.h:58
Generated for CMSSW Reference Manual by doxygen 1.8.11