CMS 3D CMS Logo
GenPurposeSkimmerAcceptance.h
Go to the documentation of this file.
1 #ifndef GenPurposeSkimmerAcceptance_H
2 #define GenPurposeSkimmerAcceptance_H
3 
4 // system include files
5 #include <memory>
6 
7 // user include files
8 #include "FWCore/Framework/interface/Frameworkfwd.h"
9 #include "FWCore/Framework/interface/EDAnalyzer.h"
10 
11 #include "FWCore/Framework/interface/Event.h"
12 #include "FWCore/Framework/interface/MakerMacros.h"
13 #include "FWCore/ParameterSet/interface/ParameterSet.h"
14 #include "FWCore/MessageLogger/interface/MessageLogger.h"
15 
16 // other files
17 // root + maths
18 #include "TFile.h"
19 #include "TBranch.h"
20 #include "TTree.h"
21 #include "TVector.h"
22 #include "TString.h"
23 #include "TMath.h"
24 //
25 #include "HLTrigger/HLTcore/interface/TriggerSummaryAnalyzerAOD.h"
26 #include "DataFormats/HLTReco/interface/TriggerEvent.h"
27 #include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
28 #include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"
29 #include "DataFormats/EgammaCandidates/interface/Electron.h"
30 #include "DataFormats/EgammaCandidates/interface/ElectronFwd.h"
31 //
32 #include "DataFormats/Common/interface/TriggerResults.h"
33 #include "FWCore/Common/interface/TriggerNames.h"
34 #include "DataFormats/HLTReco/interface/TriggerEventWithRefs.h"
35 #include "DataFormats/HLTReco/interface/TriggerTypeDefs.h"
36 #include "DataFormats/HLTReco/interface/TriggerRefsCollections.h"
37 #include "DataFormats/HLTReco/interface/TriggerFilterObjectWithRefs.h"
38 #include "DataFormats/HLTReco/interface/TriggerEvent.h"
39 #include "DataFormats/HLTReco/interface/TriggerObject.h"
40 
41 //
42 #include "DataFormats/Common/interface/ValueMap.h"
43 #include "DataFormats/GsfTrackReco/interface/GsfTrack.h"
44 #include "DataFormats/TrackReco/interface/Track.h"
45 #include "DataFormats/TrackReco/interface/TrackFwd.h"
46 #include "DataFormats/METReco/interface/CaloMET.h"
47 #include "DataFormats/METReco/interface/CaloMETFwd.h"
48 #include "DataFormats/METReco/interface/PFMET.h"
49 #include "DataFormats/METReco/interface/PFMETFwd.h"
50 #include "DataFormats/METReco/interface/GenMET.h"
51 #include "DataFormats/METReco/interface/GenMETFwd.h"
52 
53 #include "RecoEcal/EgammaCoreTools/interface/EcalClusterLazyTools.h"
54 #include "RecoEgamma/ElectronIdentification/interface/ElectronIDAlgo.h"
55 //
56 #include "DataFormats/HepMCCandidate/interface/GenParticle.h"
57 #include "DataFormats/RecoCandidate/interface/IsoDeposit.h"
58 #include "DataFormats/RecoCandidate/interface/IsoDepositFwd.h"
59 #include "DataFormats/Common/interface/ValueMap.h"
60 
61 #include "DataFormats/PatCandidates/interface/Electron.h"
62 #include "DataFormats/PatCandidates/interface/MET.h"
63 #include "DataFormats/Math/interface/deltaR.h"
64 //
65 // class decleration
66 //
67 
68 class GenPurposeSkimmerAcceptance : public edm::EDAnalyzer {
69  public:
70  explicit GenPurposeSkimmerAcceptance(const edm::ParameterSet&);
71  ~GenPurposeSkimmerAcceptance();
72 
73 
74  private:
75  virtual void beginJob(const edm::EventSetup&) ;
76  virtual void analyze(const edm::Event&, const edm::EventSetup&);
77  virtual void endJob() ;
78 
79  // ----------member data ---------------------------
80 
81 
82  std::string outputFile_;
83  int tree_fills_;
84 
85  edm::InputTag ElectronCollection_;
86  edm::InputTag MCCollection_;
87  edm::InputTag MetCollectionTag_;
88  edm::InputTag tcMetCollectionTag_;
89  edm::InputTag pfMetCollectionTag_;
90  edm::InputTag t1MetCollectionTag_;
91  edm::InputTag t1MetCollectionTagTwiki_;
92  edm::InputTag genMetCollectionTag_;
93  //
94  edm::InputTag HLTCollectionE29_;
95  edm::InputTag HLTCollectionE31_;
96  edm::InputTag HLTTriggerResultsE29_;
97  edm::InputTag HLTTriggerResultsE31_;
98  edm::InputTag HLTFilterType_[25];
99  std::string HLTPath_[25];
100 
101  TTree * probe_tree;
102  TFile * histofile;
103  //
104 
105 
106  //probe SC variables
107  double probe_sc_eta_for_tree[4];
108  double probe_sc_phi_for_tree[4];
109  double probe_sc_et_for_tree[4];
110  int probe_sc_pass_fiducial_cut[4];
111  int probe_sc_pass_et_cut[4];
112 
113  //probe electron variables
114  double probe_ele_eta_for_tree[4];
115  double probe_ele_phi_for_tree[4];
116  double probe_ele_et_for_tree[4];
117  double probe_ele_Xvertex_for_tree[4];
118  double probe_ele_Yvertex_for_tree[4];
119  double probe_ele_Zvertex_for_tree[4];
120  double probe_ele_tip[4];
121  int probe_charge_for_tree[4];
122  int probe_index_for_tree[4];
123 
124  //efficiency cuts
125  int probe_ele_pass_fiducial_cut[4];
126  int probe_ele_pass_et_cut[4];
127  int probe_pass_recoEle_cut[4];
128  int probe_pass_iso_cut[4];
129  //
130  double probe_isolation_value[4];
131  double probe_iso_user[4];
132  //
133  double probe_ecal_isolation_value[4];
134  double probe_ecal_iso_user[4];
135 
136  double probe_hcal_isolation_value[4];
137  double probe_hcal_iso_user[4];
138  //
139  int probe_classification_index_for_tree[4];
140  int probe_pass_tip_cut[4];
141  //
142  int probe_pass_id_robust_loose[4];
143  int probe_pass_id_robust_tight[4];
144  int probe_pass_id_loose[4];
145  int probe_pass_id_tight[4];
146  double probe_ele_hoe[4];
147  double probe_ele_shh[4];
148  double probe_ele_sihih[4];
149  double probe_ele_dhi[4];
150  double probe_ele_dfi[4];
151  double probe_ele_eop[4];
152  double probe_ele_pin[4];
153  double probe_ele_pout[4];
154  double probe_ele_e5x5[4];
155  double probe_ele_e2x5[4];
156  double probe_ele_e1x5[4];
157 
158  //
159  int probe_pass_trigger_cut[4][25];
160  double probe_hlt_matched_dr[4];
161  //
162  double MCMatch_Deta_;
163  double MCMatch_Dphi_;
164  int probe_mc_matched[4];
165  double probe_mc_matched_deta[4];
166  double probe_mc_matched_dphi[4];
167  double probe_mc_matched_denergy[4];
168  int probe_mc_matched_mother[4];
169 
170  //event variables
171  int numberOfHLTFilterObjects[25];
172 
173  int event_HLTPath[25];
174 
175  double BarrelMaxEta;
176  double EndcapMinEta;
177  double EndcapMaxEta;
178 
179  double ProbeSCMinEt;
180  double ProbeRecoEleSCMaxDE;
181 
182  double ProbeHLTObjMaxDR;
183  double RecoEleSeedBCMaxDE;
184  double GsfTrackMinInnerPt;
185 
186  int elec_number_in_event;
187  int elec_1_duplicate_removal;
188 
189  double event_MET, event_MET_sig;
190  // double event_MET_eta;
191  double event_MET_phi;
192  double event_tcMET, event_tcMET_sig, event_tcMET_phi;
193  double event_pfMET, event_pfMET_sig;
194  // double event_pfMET_eta;
195  double event_pfMET_phi;
196  double event_genMET, event_genMET_sig;
197  // double event_genMET_eta;
198  double event_genMET_phi;
199  //
200  double event_t1MET, event_t1MET_phi, event_t1MET_sig;
201  double event_twikiT1MET, event_twikiT1MET_phi, event_twikiT1MET_sig;
202  //
203  //
204  // acceptance and t&p specific
205  double mc_ele_eta[10];
206  double mc_ele_phi[10];
207  double mc_ele_et[10];
208  double mc_ele_vertex_x[10];
209  double mc_ele_vertex_y[10];
210  double mc_ele_vertex_z[10];
211  int mc_ele_mother[10];
212  int mc_ele_charge[10];
213  int mc_ele_status[10];
214  //
215  // for the sc collections
216  // double sc0_eta[8], sc0_phi[8], sc0_et[8];
217  double sc1_eta[8], sc1_phi[8], sc1_et[8];
218  double sc2_eta[8], sc2_phi[8], sc2_et[8];
219  double sc3_eta[8], sc3_phi[8], sc3_et[8];
220  double sc4_eta[8], sc4_phi[8], sc4_et[8];
221  double sc5_eta[8], sc5_phi[8], sc5_et[8];
222  // double sc6_eta[8], sc6_phi[8], sc6_et[8];
223  // double sc7_eta[8], sc7_phi[8], sc7_et[8];
224 
225 };
226 
227 
228 #endif
229 
edm::EDAnalyzer::beginJob
virtual void beginJob()
Definition: EDAnalyzer.h:85
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
GenPurposeSkimmerAcceptance::analyze
virtual void analyze(const edm::Event &, const edm::EventSetup &)
GenPurposeSkimmerAcceptance::GenPurposeSkimmerAcceptance
GenPurposeSkimmerAcceptance(const edm::ParameterSet &)
Generated for CMSSW Reference Manual by doxygen 1.8.11