GenPurposeSkimmerAcceptance.h

Go to the documentation of this file.

#ifndef GenPurposeSkimmerAcceptance_H
#define GenPurposeSkimmerAcceptance_H
 
// system include files
#include <memory>
 
// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
 
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
 
// other files
// root + maths
#include "TFile.h"
#include "TBranch.h"
#include "TTree.h"
#include "TVector.h"
#include "TString.h"
#include "TMath.h"
//
#include "DataFormats/HLTReco/interface/TriggerEvent.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"
#include "DataFormats/EgammaCandidates/interface/Electron.h"
#include "DataFormats/EgammaCandidates/interface/ElectronFwd.h"
//
#include "DataFormats/Common/interface/TriggerResults.h"
#include "FWCore/Common/interface/TriggerNames.h"
#include "DataFormats/HLTReco/interface/TriggerEventWithRefs.h"
#include "DataFormats/HLTReco/interface/TriggerTypeDefs.h"
#include "DataFormats/HLTReco/interface/TriggerRefsCollections.h"
#include "DataFormats/HLTReco/interface/TriggerFilterObjectWithRefs.h"
#include "DataFormats/HLTReco/interface/TriggerEvent.h"
#include "DataFormats/HLTReco/interface/TriggerObject.h"
 
//
#include "DataFormats/Common/interface/ValueMap.h"
#include "DataFormats/GsfTrackReco/interface/GsfTrack.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackReco/interface/TrackFwd.h"
#include "DataFormats/METReco/interface/CaloMET.h"
#include "DataFormats/METReco/interface/CaloMETFwd.h"
#include "DataFormats/METReco/interface/PFMET.h"
#include "DataFormats/METReco/interface/PFMETFwd.h"
#include "DataFormats/METReco/interface/GenMET.h"
#include "DataFormats/METReco/interface/GenMETFwd.h"
 
#include "RecoEcal/EgammaCoreTools/interface/EcalClusterLazyTools.h"
//
#include "DataFormats/HepMCCandidate/interface/GenParticle.h"
#include "DataFormats/RecoCandidate/interface/IsoDeposit.h"
#include "DataFormats/RecoCandidate/interface/IsoDepositFwd.h"
#include "DataFormats/Common/interface/ValueMap.h"
 
#include "DataFormats/PatCandidates/interface/Electron.h"
#include "DataFormats/PatCandidates/interface/MET.h"
#include "DataFormats/Math/interface/deltaR.h"
//
// class decleration
//
 
class GenPurposeSkimmerAcceptance : public edm::EDAnalyzer {
public:
  explicit GenPurposeSkimmerAcceptance(const edm::ParameterSet&);
  ~GenPurposeSkimmerAcceptance();
 
private:
  virtual void beginJob(const edm::EventSetup&);
  virtual void analyze(const edm::Event&, const edm::EventSetup&);
  virtual void endJob();
 
  // ----------member data ---------------------------
 
  std::string outputFile_;
  int tree_fills_;
 
  edm::InputTag ElectronCollection_;
  edm::InputTag MCCollection_;
  edm::InputTag MetCollectionTag_;
  edm::InputTag tcMetCollectionTag_;
  edm::InputTag pfMetCollectionTag_;
  edm::InputTag t1MetCollectionTag_;
  edm::InputTag t1MetCollectionTagTwiki_;
  edm::InputTag genMetCollectionTag_;
  //
  edm::InputTag HLTCollectionE29_;
  edm::InputTag HLTCollectionE31_;
  edm::InputTag HLTTriggerResultsE29_;
  edm::InputTag HLTTriggerResultsE31_;
  edm::InputTag HLTFilterType_[25];
  std::string HLTPath_[25];
 
  TTree* probe_tree;
  TFile* histofile;
  //
 
  //probe SC variables
  double probe_sc_eta_for_tree[4];
  double probe_sc_phi_for_tree[4];
  double probe_sc_et_for_tree[4];
  int probe_sc_pass_fiducial_cut[4];
  int probe_sc_pass_et_cut[4];
 
  //probe electron variables
  double probe_ele_eta_for_tree[4];
  double probe_ele_phi_for_tree[4];
  double probe_ele_et_for_tree[4];
  double probe_ele_Xvertex_for_tree[4];
  double probe_ele_Yvertex_for_tree[4];
  double probe_ele_Zvertex_for_tree[4];
  double probe_ele_tip[4];
  int probe_charge_for_tree[4];
  int probe_index_for_tree[4];
 
  //efficiency cuts
  int probe_ele_pass_fiducial_cut[4];
  int probe_ele_pass_et_cut[4];
  int probe_pass_recoEle_cut[4];
  int probe_pass_iso_cut[4];
  //
  double probe_isolation_value[4];
  double probe_iso_user[4];
  //
  double probe_ecal_isolation_value[4];
  double probe_ecal_iso_user[4];
 
  double probe_hcal_isolation_value[4];
  double probe_hcal_iso_user[4];
  //
  int probe_classification_index_for_tree[4];
  int probe_pass_tip_cut[4];
  //
  int probe_pass_id_robust_loose[4];
  int probe_pass_id_robust_tight[4];
  int probe_pass_id_loose[4];
  int probe_pass_id_tight[4];
  double probe_ele_hoe[4];
  double probe_ele_shh[4];
  double probe_ele_sihih[4];
  double probe_ele_dhi[4];
  double probe_ele_dfi[4];
  double probe_ele_eop[4];
  double probe_ele_pin[4];
  double probe_ele_pout[4];
  double probe_ele_e5x5[4];
  double probe_ele_e2x5[4];
  double probe_ele_e1x5[4];
 
  //
  int probe_pass_trigger_cut[4][25];
  double probe_hlt_matched_dr[4];
  //
  double MCMatch_Deta_;
  double MCMatch_Dphi_;
  int probe_mc_matched[4];
  double probe_mc_matched_deta[4];
  double probe_mc_matched_dphi[4];
  double probe_mc_matched_denergy[4];
  int probe_mc_matched_mother[4];
 
  //event variables
  int numberOfHLTFilterObjects[25];
 
  int event_HLTPath[25];
 
  double BarrelMaxEta;
  double EndcapMinEta;
  double EndcapMaxEta;
 
  double ProbeSCMinEt;
  double ProbeRecoEleSCMaxDE;
 
  double ProbeHLTObjMaxDR;
  double RecoEleSeedBCMaxDE;
  double GsfTrackMinInnerPt;
 
  int elec_number_in_event;
  int elec_1_duplicate_removal;
 
  double event_MET, event_MET_sig;
  //  double    event_MET_eta;
  double event_MET_phi;
  double event_tcMET, event_tcMET_sig, event_tcMET_phi;
  double event_pfMET, event_pfMET_sig;
  //  double    event_pfMET_eta;
  double event_pfMET_phi;
  double event_genMET, event_genMET_sig;
  //  double    event_genMET_eta;
  double event_genMET_phi;
  //
  double event_t1MET, event_t1MET_phi, event_t1MET_sig;
  double event_twikiT1MET, event_twikiT1MET_phi, event_twikiT1MET_sig;
  //
  //
  // acceptance and t&p specific
  double mc_ele_eta[10];
  double mc_ele_phi[10];
  double mc_ele_et[10];
  double mc_ele_vertex_x[10];
  double mc_ele_vertex_y[10];
  double mc_ele_vertex_z[10];
  int mc_ele_mother[10];
  int mc_ele_charge[10];
  int mc_ele_status[10];
  //
  // for the sc collections
  //  double sc0_eta[8], sc0_phi[8], sc0_et[8];
  double sc1_eta[8], sc1_phi[8], sc1_et[8];
  double sc2_eta[8], sc2_phi[8], sc2_et[8];
  double sc3_eta[8], sc3_phi[8], sc3_et[8];
  double sc4_eta[8], sc4_phi[8], sc4_et[8];
  double sc5_eta[8], sc5_phi[8], sc5_et[8];
  //  double sc6_eta[8], sc6_phi[8], sc6_et[8];
  //  double sc7_eta[8], sc7_phi[8], sc7_et[8];
};
 
#endif