/ElectroWeakAnalysis/WENu/interface/GenPurposeSkimmerAcceptance.h

https://github.com/aivanov-cern/cmssw · C Header · 229 lines · 161 code · 32 blank · 36 comment · 0 complexity · 0f7b29829664ef7c9fcdef5a90486967 MD5 · raw file 

    

  1. #ifndef GenPurposeSkimmerAcceptance_H
    2. 

  2. #define GenPurposeSkimmerAcceptance_H
    3. 

  3. 

  4. // system include files
    5. 

  5. #include <memory>
    6. 

  6. 

  7. // user include files
    8. 

  8. #include "FWCore/Framework/interface/Frameworkfwd.h"
    9. 

  9. #include "FWCore/Framework/interface/EDAnalyzer.h"
    10. 

  10. 

  11. #include "FWCore/Framework/interface/Event.h"
    12. 

  12. #include "FWCore/Framework/interface/MakerMacros.h"
    13. 

  13. #include "FWCore/ParameterSet/interface/ParameterSet.h"
    14. 

  14. #include "FWCore/MessageLogger/interface/MessageLogger.h"
    15. 

  15. 

  16. // other files
    17. 

  17. // root + maths
    18. 

  18. #include "TFile.h"
    19. 

  19. #include "TBranch.h"
    20. 

  20. #include "TTree.h"
    21. 

  21. #include "TVector.h"
    22. 

  22. #include "TString.h"
    23. 

  23. #include "TMath.h"
    24. 

  24. //
    25. 

  25. #include "HLTrigger/HLTcore/interface/TriggerSummaryAnalyzerAOD.h"
    26. 

  26. #include "DataFormats/HLTReco/interface/TriggerEvent.h"
    27. 

  27. #include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
    28. 

  28. #include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"
    29. 

  29. #include "DataFormats/EgammaCandidates/interface/Electron.h"
    30. 

  30. #include "DataFormats/EgammaCandidates/interface/ElectronFwd.h"
    31. 

  31. //
    32. 

  32. #include "DataFormats/Common/interface/TriggerResults.h"
    33. 

  33. #include "FWCore/Common/interface/TriggerNames.h"
    34. 

  34. #include "DataFormats/HLTReco/interface/TriggerEventWithRefs.h"
    35. 

  35. #include "DataFormats/HLTReco/interface/TriggerTypeDefs.h"
    36. 

  36. #include "DataFormats/HLTReco/interface/TriggerRefsCollections.h"
    37. 

  37. #include "DataFormats/HLTReco/interface/TriggerFilterObjectWithRefs.h"
    38. 

  38. #include "DataFormats/HLTReco/interface/TriggerEvent.h"
    39. 

  39. #include "DataFormats/HLTReco/interface/TriggerObject.h"
    40. 

  40. 

  41. //
    42. 

  42. #include "DataFormats/Common/interface/ValueMap.h"
    43. 

  43. #include "DataFormats/GsfTrackReco/interface/GsfTrack.h"
    44. 

  44. #include "DataFormats/TrackReco/interface/Track.h"
    45. 

  45. #include "DataFormats/TrackReco/interface/TrackFwd.h"
    46. 

  46. #include "DataFormats/METReco/interface/CaloMET.h"
    47. 

  47. #include "DataFormats/METReco/interface/CaloMETFwd.h"
    48. 

  48. #include "DataFormats/METReco/interface/PFMET.h"
    49. 

  49. #include "DataFormats/METReco/interface/PFMETFwd.h"
    50. 

  50. #include "DataFormats/METReco/interface/GenMET.h"
    51. 

  51. #include "DataFormats/METReco/interface/GenMETFwd.h"
    52. 

  52. 

  53. #include "RecoEcal/EgammaCoreTools/interface/EcalClusterLazyTools.h"
    54. 

  54. #include "RecoEgamma/ElectronIdentification/interface/ElectronIDAlgo.h"
    55. 

  55. //
    56. 

  56. #include "DataFormats/HepMCCandidate/interface/GenParticle.h"
    57. 

  57. #include "DataFormats/RecoCandidate/interface/IsoDeposit.h"
    58. 

  58. #include "DataFormats/RecoCandidate/interface/IsoDepositFwd.h"
    59. 

  59. #include "DataFormats/Common/interface/ValueMap.h"
    60. 

  60. 

  61. #include "DataFormats/PatCandidates/interface/Electron.h"
    62. 

  62. #include "DataFormats/PatCandidates/interface/MET.h"
    63. 

  63. #include "DataFormats/Math/interface/deltaR.h"
    64. 

  64. //
    65. 

  65. // class decleration
    66. 

  66. //
    67. 

  67. 

  68. class GenPurposeSkimmerAcceptance : public edm::EDAnalyzer {
    69. 

  69.    public:
    70. 

  70.       explicit GenPurposeSkimmerAcceptance(const edm::ParameterSet&);
    71. 

  71.       ~GenPurposeSkimmerAcceptance();
    72. 

  72. 

  73. 

  74.    private:
    75. 

  75.       virtual void beginJob(const edm::EventSetup&) ;
    76. 

  76.       virtual void analyze(const edm::Event&, const edm::EventSetup&);
    77. 

  77.       virtual void endJob() ;
    78. 

  78. 

  79.       // ----------member data ---------------------------
    80. 

  80. 

  81.   
    82. 

  82.   std::string outputFile_;
    83. 

  83.   int tree_fills_;
    84. 

  84. 

  85.   edm::InputTag ElectronCollection_;
    86. 

  86.   edm::InputTag MCCollection_;
    87. 

  87.   edm::InputTag MetCollectionTag_;
    88. 

  88.   edm::InputTag tcMetCollectionTag_;
    89. 

  89.   edm::InputTag pfMetCollectionTag_;
    90. 

  90.   edm::InputTag t1MetCollectionTag_;
    91. 

  91.   edm::InputTag t1MetCollectionTagTwiki_;
    92. 

  92.   edm::InputTag genMetCollectionTag_;
    93. 

  93.   //
    94. 

  94.   edm::InputTag HLTCollectionE29_;
    95. 

  95.   edm::InputTag HLTCollectionE31_;
    96. 

  96.   edm::InputTag HLTTriggerResultsE29_;
    97. 

  97.   edm::InputTag HLTTriggerResultsE31_;
    98. 

  98.   edm::InputTag HLTFilterType_[25];
    99. 

  99.   std::string HLTPath_[25];
    100. 

  100. 

  101.   TTree * probe_tree;
    102. 

  102.   TFile * histofile;
    103. 

  103.   //
    104. 

  104. 

  105. 

  106.   //probe SC variables
    107. 

  107.   double probe_sc_eta_for_tree[4];
    108. 

  108.   double probe_sc_phi_for_tree[4];
    109. 

  109.   double probe_sc_et_for_tree[4];
    110. 

  110.   int probe_sc_pass_fiducial_cut[4]; 
    111. 

  111.   int probe_sc_pass_et_cut[4]; 
    112. 

  112.   
    113. 

  113.   //probe electron variables
    114. 

  114.   double probe_ele_eta_for_tree[4];
    115. 

  115.   double probe_ele_phi_for_tree[4];
    116. 

  116.   double probe_ele_et_for_tree[4];
    117. 

  117.   double probe_ele_Xvertex_for_tree[4];
    118. 

  118.   double probe_ele_Yvertex_for_tree[4];
    119. 

  119.   double probe_ele_Zvertex_for_tree[4];
    120. 

  120.   double probe_ele_tip[4];
    121. 

  121.   int probe_charge_for_tree[4];
    122. 

  122.   int probe_index_for_tree[4];
    123. 

  123.   
    124. 

  124.   //efficiency cuts
    125. 

  125.   int probe_ele_pass_fiducial_cut[4]; 
    126. 

  126.   int probe_ele_pass_et_cut[4]; 
    127. 

  127.   int probe_pass_recoEle_cut[4]; 
    128. 

  128.   int probe_pass_iso_cut[4]; 
    129. 

  129.   //
    130. 

  130.   double probe_isolation_value[4]; 
    131. 

  131.   double probe_iso_user[4]; 
    132. 

  132.   //
    133. 

  133.   double probe_ecal_isolation_value[4];
    134. 

  134.   double probe_ecal_iso_user[4];
    135. 

  135. 

  136.   double probe_hcal_isolation_value[4];
    137. 

  137.   double probe_hcal_iso_user[4];
    138. 

  138.   //
    139. 

  139.   int probe_classification_index_for_tree[4]; 
    140. 

  140.   int probe_pass_tip_cut[4];
    141. 

  141.   //
    142. 

  142.   int probe_pass_id_robust_loose[4];
    143. 

  143.   int probe_pass_id_robust_tight[4];
    144. 

  144.   int probe_pass_id_loose[4];
    145. 

  145.   int probe_pass_id_tight[4];
    146. 

  146.   double probe_ele_hoe[4];
    147. 

  147.   double probe_ele_shh[4];
    148. 

  148.   double probe_ele_sihih[4];
    149. 

  149.   double probe_ele_dhi[4];
    150. 

  150.   double probe_ele_dfi[4];
    151. 

  151.   double probe_ele_eop[4];
    152. 

  152.   double probe_ele_pin[4];
    153. 

  153.   double probe_ele_pout[4];
    154. 

  154.   double probe_ele_e5x5[4];
    155. 

  155.   double probe_ele_e2x5[4];
    156. 

  156.   double probe_ele_e1x5[4];
    157. 

  157. 

  158.   //
    159. 

  159.   int probe_pass_trigger_cut[4][25];   
    160. 

  160.   double probe_hlt_matched_dr[4];
    161. 

  161.   //
    162. 

  162.   double MCMatch_Deta_;
    163. 

  163.   double MCMatch_Dphi_;
    164. 

  164.   int probe_mc_matched[4];
    165. 

  165.   double probe_mc_matched_deta[4];
    166. 

  166.   double probe_mc_matched_dphi[4];
    167. 

  167.   double probe_mc_matched_denergy[4];
    168. 

  168.   int probe_mc_matched_mother[4];
    169. 

  169.   
    170. 

  170.  //event variables
    171. 

  171.   int numberOfHLTFilterObjects[25];
    172. 

  172. 

  173.   int event_HLTPath[25];
    174. 

  174.   
    175. 

  175.   double BarrelMaxEta;
    176. 

  176.   double EndcapMinEta;
    177. 

  177.   double EndcapMaxEta;
    178. 

  178. 

  179.   double ProbeSCMinEt;
    180. 

  180.   double ProbeRecoEleSCMaxDE; 
    181. 

  181. 

  182.   double ProbeHLTObjMaxDR;
    183. 

  183.   double RecoEleSeedBCMaxDE;
    184. 

  184.   double GsfTrackMinInnerPt;
    185. 

  185. 

  186.   int elec_number_in_event;      
    187. 

  187.   int elec_1_duplicate_removal;
    188. 

  188. 

  189.   double    event_MET, event_MET_sig; 
    190. 

  190.   //  double    event_MET_eta;
    191. 

  191.   double    event_MET_phi;
    192. 

  192.   double    event_tcMET, event_tcMET_sig, event_tcMET_phi;
    193. 

  193.   double    event_pfMET, event_pfMET_sig; 
    194. 

  194.   //  double    event_pfMET_eta;
    195. 

  195.   double    event_pfMET_phi;
    196. 

  196.   double    event_genMET, event_genMET_sig; 
    197. 

  197.   //  double    event_genMET_eta; 
    198. 

  198.   double    event_genMET_phi;
    199. 

  199.   //
    200. 

  200.   double event_t1MET, event_t1MET_phi, event_t1MET_sig;
    201. 

  201.   double event_twikiT1MET, event_twikiT1MET_phi, event_twikiT1MET_sig;
    202. 

  202.   //
    203. 

  203.   //
    204. 

  204.   // acceptance and t&p specific
    205. 

  205.   double   mc_ele_eta[10];
    206. 

  206.   double   mc_ele_phi[10];
    207. 

  207.   double   mc_ele_et[10];
    208. 

  208.   double   mc_ele_vertex_x[10];
    209. 

  209.   double   mc_ele_vertex_y[10];
    210. 

  210.   double   mc_ele_vertex_z[10];
    211. 

  211.   int   mc_ele_mother[10];
    212. 

  212.   int   mc_ele_charge[10];
    213. 

  213.   int   mc_ele_status[10];
    214. 

  214.   //
    215. 

  215.   // for the sc collections
    216. 

  216.   //  double sc0_eta[8], sc0_phi[8], sc0_et[8];
    217. 

  217.   double sc1_eta[8], sc1_phi[8], sc1_et[8];
    218. 

  218.   double sc2_eta[8], sc2_phi[8], sc2_et[8];
    219. 

  219.   double sc3_eta[8], sc3_phi[8], sc3_et[8];
    220. 

  220.   double sc4_eta[8], sc4_phi[8], sc4_et[8];
    221. 

  221.   double sc5_eta[8], sc5_phi[8], sc5_et[8];
    222. 

  222.   //  double sc6_eta[8], sc6_phi[8], sc6_et[8];
    223. 

  223.   //  double sc7_eta[8], sc7_phi[8], sc7_et[8];
    224. 

  224. 

  225. };
    226. 

  226. 

  227. 

  228. #endif
    229.