db/d1b/ElectronAnalyzer_8h_source.html

#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