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