ElectronAnalyzer.h

Go to the documentation of this file.

#ifndef DQMOffline_EGamma_ElectronAnalyzer_h
#define DQMOffline_EGamma_ElectronAnalyzer_h

#include "DQMOffline/EGamma/interface/ElectronDqmAnalyzerBase.h"


#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
#include "TrackingTools/TrajectoryState/interface/TrajectoryStateTransform.h"

#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/GsfTrackReco/interface/GsfTrack.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectronCore.h"
#include "DataFormats/EgammaReco/interface/SuperClusterFwd.h"
#include "DataFormats/EgammaReco/interface/BasicClusterFwd.h"
#include "DataFormats/EgammaReco/interface/ElectronSeed.h"
#include "DataFormats/EgammaReco/interface/ElectronSeedFwd.h"
#include "DataFormats/EcalDetId/interface/EcalSubdetector.h"
#include "DataFormats/HepMCCandidate/interface/GenParticle.h"
#include "DataFormats/HepMCCandidate/interface/GenParticleFwd.h"
#include "DataFormats/BeamSpot/interface/BeamSpot.h"
#include "DataFormats/JetReco/interface/GenJetCollection.h"
 
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/Common/interface/ValueMap.h"

class MagneticField ;

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"

class ElectronAnalyzer : public ElectronDqmAnalyzerBase
 {
  public:

    explicit ElectronAnalyzer(const edm::ParameterSet& conf);
    ~ElectronAnalyzer() override;

//    virtual void book() ;
    void bookHistograms( DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override ;
    void analyze( const edm::Event & e, const edm::EventSetup & c) override ;

  private:

    //=========================================
    // parameters
    //=========================================

    // general, collections
    int Selection_;
    edm::EDGetTokenT<reco::GsfElectronCollection> electronCollection_;
    edm::EDGetTokenT<reco::SuperClusterCollection> matchingObjectCollection_;
    edm::EDGetTokenT<reco::GsfTrackCollection> gsftrackCollection_;
    edm::EDGetTokenT<reco::TrackCollection> trackCollection_;
    edm::EDGetTokenT<reco::VertexCollection> vertexCollection_;
    edm::EDGetTokenT<reco::BeamSpot> beamSpotTag_;
    bool readAOD_; //NEW

    std::string inputFile_ ;
    std::string outputFile_ ;
    std::string inputInternalPath_ ;
    std::string outputInternalPath_ ;

    // matching
    std::string matchingCondition_; //NEW
    double maxPtMatchingObject_; // SURE ?
    double maxAbsEtaMatchingObject_; // SURE ?
    double deltaR_;

    // electron selection NEW
    double minEt_;
    double minPt_;
    double maxAbsEta_;
    bool isEB_;
    bool isEE_;
    bool isNotEBEEGap_;
    bool isEcalDriven_;
    bool isTrackerDriven_;
    double eOverPMinBarrel_;
    double eOverPMaxBarrel_;
    double eOverPMinEndcaps_;
    double eOverPMaxEndcaps_;
    double dEtaMinBarrel_;
    double dEtaMaxBarrel_;
    double dEtaMinEndcaps_;
    double dEtaMaxEndcaps_;
    double dPhiMinBarrel_;
    double dPhiMaxBarrel_;
    double dPhiMinEndcaps_;
    double dPhiMaxEndcaps_;
    double sigIetaIetaMinBarrel_;
    double sigIetaIetaMaxBarrel_;
    double sigIetaIetaMinEndcaps_;
    double sigIetaIetaMaxEndcaps_;
    double hadronicOverEmMaxBarrel_;
    double hadronicOverEmMaxEndcaps_;
    double mvaMin_;
    double tipMaxBarrel_;
    double tipMaxEndcaps_;
    double tkIso03Max_;
    double hcalIso03Depth1MaxBarrel_;
    double hcalIso03Depth1MaxEndcaps_;
    double hcalIso03Depth2MaxEndcaps_;
    double ecalIso03MaxBarrel_;
    double ecalIso03MaxEndcaps_;

    // for trigger NEW
    edm::InputTag triggerResults_;
//    std::vector<std::string > HLTPathsByName_;

    // histos limits and binning
    int nbineta; int nbineta2D; double etamin; double etamax;
    int nbinphi; int nbinphi2D; double phimin; double phimax;
    int nbinpt; int nbinpteff; int nbinpt2D; double ptmax;
    int nbinp; int nbinp2D; double pmax;
    int nbineop; int nbineop2D; double eopmax; double eopmaxsht;
    int nbindeta; double detamin; double detamax;
    int nbindphi; double dphimin; double dphimax;
    int nbindetamatch; int nbindetamatch2D; double detamatchmin; double detamatchmax;
    int nbindphimatch; int nbindphimatch2D; double dphimatchmin; double dphimatchmax;
    int nbinfhits; double fhitsmax;
    int nbinlhits; double lhitsmax;
    int nbinxyz; int nbinxyz2D;
    int nbinpoptrue; double poptruemin; double poptruemax; //NEW
    int nbinmee; double meemin; double meemax; //NEW
    int nbinhoe; double hoemin; double hoemax; //NEW
//    bool set_EfficiencyFlag ; bool set_StatOverflowFlag ;

    //=========================================
    // general attributes and utility methods
    //=========================================

    unsigned int nEvents_ ;

    float computeInvMass
     ( const reco::GsfElectron & e1,
       const reco::GsfElectron & e2 ) ;

    bool selected( const reco::GsfElectronCollection::const_iterator & gsfIter , double vertexTIP ) ;
    bool generalCut( const reco::GsfElectronCollection::const_iterator & gsfIter) ;
    bool etCut( const reco::GsfElectronCollection::const_iterator & gsfIter ) ;
    bool isolationCut( const reco::GsfElectronCollection::const_iterator & gsfIter, double vertexTIP ) ;
    bool idCut( const reco::GsfElectronCollection::const_iterator & gsfIter ) ;

//    bool trigger( const edm::Event & e ) ;
//    unsigned int nAfterTrigger_;
//    std::vector<unsigned int> HLTPathsByIndex_;

    edm::ESHandle<TrackerGeometry> pDD;
    edm::ESHandle<MagneticField> theMagField;

    float mcEnergy[10], mcEta[10], mcPhi[10], mcPt[10], mcQ[10];
    float superclusterEnergy[10], superclusterEta[10], superclusterPhi[10], superclusterEt[10];
    float seedMomentum[10], seedEta[10], seedPhi[10], seedPt[10], seedQ[10];

    //=========================================
    // histograms
    //=========================================

    // general
    MonitorElement * h2_beamSpotXvsY ;
    MonitorElement * py_nElectronsVsLs ;
    MonitorElement * py_nClustersVsLs ;
    MonitorElement * py_nGsfTracksVsLs ;
    MonitorElement * py_nTracksVsLs ;
    MonitorElement * py_nVerticesVsLs ;
    MonitorElement * h1_triggers ;

    // basic quantities
//    MonitorElement * h1_num_ ; // number of electrons in a single event
//    MonitorElement * h1_charge ;
//    MonitorElement * h1_vertexP ;
//    MonitorElement * h1_Et ;
//    MonitorElement * h1_vertexTIP ;
//    MonitorElement * h1_vertexPhi ;
//    MonitorElement * h1_vertexX ;
//    MonitorElement * h1_vertexY ;
    MonitorElement * h1_vertexPt_barrel ;
    MonitorElement * h1_vertexPt_endcaps ;
    MonitorElement * h1_vertexEta ;
    MonitorElement * h2_vertexEtaVsPhi ;
    MonitorElement * h2_vertexXvsY ;
    MonitorElement * h1_vertexZ ;

    // super-clusters
//    MonitorElement * h1_sclEn ;
//    MonitorElement * h1_sclEta ;
//    MonitorElement * h1_sclPhi ;
    MonitorElement * h1_sclEt ;

    // gsf tracks
//    MonitorElement * h1_ambiguousTracks ;
//    MonitorElement * h2ele_ambiguousTracksVsEta ;
//    MonitorElement * h2_ambiguousTracksVsPhi ;
//    MonitorElement * h2_ambiguousTracksVsPt ;
    MonitorElement * h1_chi2 ;
    MonitorElement * py_chi2VsEta ;
    MonitorElement * py_chi2VsPhi ;
//    MonitorElement * h2_chi2VsPt ;
    MonitorElement * h1_foundHits ;
    MonitorElement * py_foundHitsVsEta ;
    MonitorElement * py_foundHitsVsPhi ;
//    MonitorElement * h2_foundHitsVsPt ;
    MonitorElement * h1_lostHits ;
    MonitorElement * py_lostHitsVsEta ;
    MonitorElement * py_lostHitsVsPhi ;
//    MonitorElement * h2_lostHitsVsPt ;

    // electron matching and ID
    //MonitorElement * h_EopOut ;
    //MonitorElement * h_dEtaCl_propOut ;
    //MonitorElement * h_dPhiCl_propOut ;
//    MonitorElement * h1_Eop ;
//    MonitorElement * h2_EopVsEta ;
    MonitorElement * h1_Eop_barrel ;
    MonitorElement * h1_Eop_endcaps ;
    MonitorElement * py_EopVsPhi ;
//    MonitorElement * h1_EopVsPt ;
//    MonitorElement * h1_EeleOPout ;
//    MonitorElement * h2_EeleOPoutVsEta ;
    MonitorElement * h1_EeleOPout_barrel ;
    MonitorElement * h1_EeleOPout_endcaps ;
//    MonitorElement * h2_EeleOPoutVsPhi ;
//    MonitorElement * h2_EeleOPoutVsPt ;
//    MonitorElement * h1_dEtaSc_propVtx ;
//    MonitorElement * h2_dEtaSc_propVtxVsEta ;
    MonitorElement * h1_dEtaSc_propVtx_barrel ;
    MonitorElement * h1_dEtaSc_propVtx_endcapsPos ;
    MonitorElement * h1_dEtaSc_propVtx_endcapsNeg ;
    MonitorElement * py_dEtaSc_propVtxVsPhi ;
//    MonitorElement * h2_dEtaSc_propVtxVsPt ;
//    MonitorElement * h1_dEtaEleCl_propOut ;
//    MonitorElement * h2_dEtaEleCl_propOutVsEta ;
    MonitorElement * h1_dEtaEleCl_propOut_barrel ;
    MonitorElement * h1_dEtaEleCl_propOut_endcapsPos ;
    MonitorElement * h1_dEtaEleCl_propOut_endcapsNeg ;
//    MonitorElement * h2_dEtaEleCl_propOutVsPhi ;
//    MonitorElement * h2_dEtaEleCl_propOutVsPt ;
//    MonitorElement * h1_dPhiSc_propVtx ;
//    MonitorElement * h2_dPhiSc_propVtxVsEta ;
    MonitorElement * h1_dPhiSc_propVtx_barrel ;
    MonitorElement * h1_dPhiSc_propVtx_endcapsPos ;
    MonitorElement * h1_dPhiSc_propVtx_endcapsNeg ;
    MonitorElement * py_dPhiSc_propVtxVsPhi ;
//    MonitorElement * h2_dPhiSc_propVtxVsPt ;
//    MonitorElement * h1_dPhiEleCl_propOut ;
//    MonitorElement * h2_dPhiEleCl_propOutVsEta ;
    MonitorElement * h1_dPhiEleCl_propOut_barrel ;
    MonitorElement * h1_dPhiEleCl_propOut_endcapsPos ;
    MonitorElement * h1_dPhiEleCl_propOut_endcapsNeg ;
//    MonitorElement * h2_dPhiEleCl_propOutVsPhi ;
//    MonitorElement * h2_dPhiEleCl_propOutVsPt ;
//    MonitorElement * h1_Hoe ;
//    MonitorElement * h2_HoeVsEta ;
    MonitorElement * h1_Hoe_barrel ;
    MonitorElement * h1_Hoe_endcaps ;
    MonitorElement * py_HoeVsPhi ;
//    MonitorElement * h2_HoeVsPt ;
    MonitorElement * h1_sclSigEtaEta_barrel ;
    MonitorElement * h1_sclSigEtaEta_endcaps ;
    MonitorElement * h1_sigIEtaIEta5x5_barrel ;
    MonitorElement * h1_sigIEtaIEta5x5_endcaps ;

    // fbrem related variables
    //MonitorElement * h_outerP ;
    //MonitorElement * h_outerP_mode ;
//    MonitorElement * h_innerPt_mean ;
//    MonitorElement * h_outerPt_mean ;
//    MonitorElement * h_outerPt_mode ;
//    MonitorElement * h_PinMnPout ;
//    MonitorElement * h_PinMnPout_mode ;
    MonitorElement * h1_fbrem ;
    MonitorElement * py_fbremVsEta ;
    MonitorElement * py_fbremVsPhi ;
//    MonitorElement * h2_fbremVsPt ;
    MonitorElement * h1_classes ;

    // pflow
    MonitorElement * h1_mva ;
    MonitorElement * h1_provenance ;

    // isolation
    MonitorElement * h1_tkSumPt_dr03 ;
    MonitorElement * h1_ecalRecHitSumEt_dr03 ;
    MonitorElement * h1_hcalTowerSumEt_dr03 ;
    MonitorElement * h1_PFch_dr03 ;
    MonitorElement * h1_PFem_dr03 ;
    MonitorElement * h1_PFnh_dr03 ;
//    MonitorElement * h1_hcalDepth1TowerSumEt_dr03 ;
//    MonitorElement * h1_hcalDepth2TowerSumEt_dr03 ;
//    MonitorElement * h1_tkSumPt_dr04 ;
//    MonitorElement * h1_ecalRecHitSumEt_dr04 ;
//    MonitorElement * h1_hcalTowerSumEt_dr04 ;

    // di-electron mass
    MonitorElement * h1_mee ;
    MonitorElement * h1_mee_os ;
    MonitorElement * h1_mee_os_bb ;
    MonitorElement * h1_mee_os_ee ;
    MonitorElement * h1_mee_os_eb ;


    // histos for matching and matched objects

//    MonitorElement * h1_matchedEle_eta ;
//    MonitorElement * h1_matchedEle_eta_golden ;
//    MonitorElement * h1_matchedEle_eta_shower ;
//    //MonitorElement * h1_matchedEle_eta_bbrem ;
//    //MonitorElement * h1_matchedEle_eta_narrow ;

    MonitorElement * h1_matchedObject_Eta ;
//    MonitorElement * h1_matchedObject_AbsEta ;
    MonitorElement * h1_matchedObject_Pt ;
    MonitorElement * h1_matchedObject_Phi ;
//    MonitorElement * h1_matchedObject_Z ;

//    MonitorElement * h1_matchingObject_Num ;
    MonitorElement * h1_matchingObject_Eta ;
//    MonitorElement * h1_matchingObject_AbsEta ;
//    MonitorElement * h1_matchingObject_P ;
    MonitorElement * h1_matchingObject_Pt ;
    MonitorElement * h1_matchingObject_Phi ;
//    MonitorElement * h1_matchingObject_Z ;

 } ;

#endif