LepHTMonitor Class Reference

Inheritance diagram for LepHTMonitor:

Public Types
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache
Public Types inherited from TriggerDQMBase
typedef dqm::legacy::DQMStore	DQMStore
typedef dqm::legacy::MonitorElement	MonitorElement

Public Member Functions
	LepHTMonitor (const edm::ParameterSet &ps)
	~LepHTMonitor () override throw ()
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final
Public Member Functions inherited from TriggerDQMBase
void	bookME (DQMStore::IBooker &, ObjME &me, const std::string &histname, const std::string &histtitle, const std::vector< double > &binningX, const bool bookDen=true)
void	bookME (DQMStore::IBooker &, ObjME &me, const std::string &histname, const std::string &histtitle, const std::vector< double > &binningX, const std::vector< double > &binningY, const bool bookDen=true)
void	bookME (DQMStore::IBooker &, ObjME &me, const std::string &histname, const std::string &histtitle, const uint nbins, const double xmin, const double xmax, const bool bookDen=true)
void	bookME (DQMStore::IBooker &, ObjME &me, const std::string &histname, const std::string &histtitle, const uint nbinsX, const double xmin, const double xmax, const double ymin, const double ymax, const bool bookDen=true)
void	bookME (DQMStore::IBooker &, ObjME &me, const std::string &histname, const std::string &histtitle, const uint nbinsX, const double xmin, const double xmax, const uint nbinsY, const double ymin, const double ymax, const bool bookDen=true)
void	setMETitle (ObjME &me, const std::string &titleX, const std::string &titleY)
	TriggerDQMBase ()=default
virtual	~TriggerDQMBase ()=default

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &eSetup) override
void	bookHistograms (DQMStore::IBooker &ibooker, const edm::Run &, const edm::EventSetup &) override
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Attributes
std::unique_ptr< GenericTriggerEventFlag >	den_HT_genTriggerEventFlag_
std::unique_ptr< GenericTriggerEventFlag >	den_lep_genTriggerEventFlag_
float	etabins_max_
float	etabins_min_
const std::string	folderName_
MonitorElement *	h_lepEtaPhiTurnOn_den_
MonitorElement *	h_lepEtaPhiTurnOn_num_
MonitorElement *	h_lepEtaTurnOn_den_
MonitorElement *	h_lepEtaTurnOn_num_
MonitorElement *	h_lepPhiTurnOn_den_
MonitorElement *	h_lepPhiTurnOn_num_
MonitorElement *	h_lepPtTurnOn_den_
MonitorElement *	h_lepPtTurnOn_num_
MonitorElement *	h_NPVTurnOn_den_
MonitorElement *	h_NPVTurnOn_num_
MonitorElement *	h_pfHTTurnOn_den_
MonitorElement *	h_pfHTTurnOn_num_
bool	hltPathsAreValid_
std::vector< double >	htbins_
double	htCut_
double	jetEtaCut_
double	jetPtCut_
double	lep_counting_threshold_
double	lep_d0_cut_b_
double	lep_d0_cut_e_
double	lep_dz_cut_b_
double	lep_dz_cut_e_
double	lep_eta_cut_
double	lep_iso_cut_
double	lep_pt_plateau_
double	metCut_
int	muonIDlevel_
int	nbins_eta_
int	nbins_npv_
int	nbins_phi_
double	nelsCut_
double	nmusCut_
float	npvbins_max_
float	npvbins_min_
std::unique_ptr< GenericTriggerEventFlag >	num_genTriggerEventFlag_
float	phibins_max_
float	phibins_min_
std::vector< double >	ptbins_
const bool	requireValidHLTPaths_
edm::EDGetTokenT< reco::BeamSpot >	theBeamSpot_
edm::InputTag	theBeamSpotTag_
edm::EDGetTokenT< reco::ConversionCollection >	theConversionCollection_
edm::InputTag	theConversionCollectionTag_
edm::EDGetTokenT< edm::View< reco::GsfElectron > >	theElectronCollection_
edm::InputTag	theElectronTag_
edm::EDGetTokenT< edm::ValueMap< bool > >	theElectronVIDMap_
edm::InputTag	theElectronVIDTag_
edm::EDGetTokenT< reco::JetTagCollection >	theJetTagCollection_
edm::InputTag	theJetTagTag_
edm::EDGetTokenT< reco::MuonCollection >	theMuonCollection_
edm::InputTag	theMuonTag_
edm::EDGetTokenT< reco::PFJetCollection >	thePfJetCollection_
edm::InputTag	thePfJetTag_
edm::EDGetTokenT< reco::PFMETCollection >	thePfMETCollection_
edm::InputTag	thePfMETTag_
edm::EDGetTokenT< reco::VertexCollection >	theVertexCollection_
edm::InputTag	theVertexCollectionTag_

Additional Inherited Members
Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)
static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)
Static Public Member Functions inherited from TriggerDQMBase
static void	fillHistoLSPSetDescription (edm::ParameterSetDescription &pset)
static void	fillHistoPSetDescription (edm::ParameterSetDescription &pset)
static MEbinning	getHistoLSPSet (const edm::ParameterSet &pset)
static MEbinning	getHistoPSet (const edm::ParameterSet &pset)
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Definition at line 36 of file LepHTMonitor.cc.

Member Typedef Documentation

DQMStore

typedef dqm::reco::DQMStore LepHTMonitor::DQMStore

Definition at line 39 of file LepHTMonitor.cc.

MonitorElement

typedef dqm::reco::MonitorElement LepHTMonitor::MonitorElement

Definition at line 38 of file LepHTMonitor.cc.

Constructor & Destructor Documentation

LepHTMonitor()

LepHTMonitor::LepHTMonitor

(

const edm::ParameterSet &

ps

)

Definition at line 229 of file LepHTMonitor.cc.

    : theElectronTag_(ps.getParameter<edm::InputTag>("electronCollection")),
      theElectronCollection_(consumes<edm::View<reco::GsfElectron> >(theElectronTag_)),
      theElectronVIDTag_(ps.getParameter<edm::InputTag>("electronVID")),
      theElectronVIDMap_(consumes<edm::ValueMap<bool> >(theElectronVIDTag_)),
      theMuonTag_(ps.getParameter<edm::InputTag>("muonCollection")),
      theMuonCollection_(consumes<reco::MuonCollection>(theMuonTag_)),
      thePfMETTag_(ps.getParameter<edm::InputTag>("pfMetCollection")),
      thePfMETCollection_(consumes<reco::PFMETCollection>(thePfMETTag_)),
      thePfJetTag_(ps.getParameter<edm::InputTag>("pfJetCollection")),
      thePfJetCollection_(consumes<reco::PFJetCollection>(thePfJetTag_)),
      theJetTagTag_(ps.getParameter<edm::InputTag>("jetTagCollection")),
      theJetTagCollection_(consumes<reco::JetTagCollection>(theJetTagTag_)),
      theVertexCollectionTag_(ps.getParameter<edm::InputTag>("vertexCollection")),
      theVertexCollection_(consumes<reco::VertexCollection>(theVertexCollectionTag_)),
      theConversionCollectionTag_(ps.getParameter<edm::InputTag>("conversionCollection")),
      theConversionCollection_(consumes<reco::ConversionCollection>(theConversionCollectionTag_)),
      theBeamSpotTag_(ps.getParameter<edm::InputTag>("beamSpot")),
      theBeamSpot_(consumes<reco::BeamSpot>(theBeamSpotTag_)),
 
      num_genTriggerEventFlag_(new GenericTriggerEventFlag(
          ps.getParameter<edm::ParameterSet>("numGenericTriggerEventPSet"), consumesCollector(), *this)),
      den_lep_genTriggerEventFlag_(new GenericTriggerEventFlag(
          ps.getParameter<edm::ParameterSet>("den_lep_GenericTriggerEventPSet"), consumesCollector(), *this)),
      den_HT_genTriggerEventFlag_(new GenericTriggerEventFlag(
          ps.getParameter<edm::ParameterSet>("den_HT_GenericTriggerEventPSet"), consumesCollector(), *this)),
 
      folderName_(ps.getParameter<std::string>("folderName")),
      requireValidHLTPaths_(ps.getParameter<bool>("requireValidHLTPaths")),
      hltPathsAreValid_(false),
 
      muonIDlevel_(ps.getUntrackedParameter<int>("muonIDlevel")),
      jetPtCut_(ps.getUntrackedParameter<double>("jetPtCut")),
      jetEtaCut_(ps.getUntrackedParameter<double>("jetEtaCut")),
      metCut_(ps.getUntrackedParameter<double>("metCut")),
      htCut_(ps.getUntrackedParameter<double>("htCut")),
 
      nmusCut_(ps.getUntrackedParameter<double>("nmus")),
      nelsCut_(ps.getUntrackedParameter<double>("nels")),
      lep_pt_plateau_(ps.getUntrackedParameter<double>("leptonPtPlateau")),
      lep_counting_threshold_(ps.getUntrackedParameter<double>("leptonCountingThreshold")),
      lep_iso_cut_(ps.getUntrackedParameter<double>("lepIsoCut")),
      lep_eta_cut_(ps.getUntrackedParameter<double>("lepEtaCut")),
      lep_d0_cut_b_(ps.getUntrackedParameter<double>("lep_d0_cut_b")),
      lep_dz_cut_b_(ps.getUntrackedParameter<double>("lep_dz_cut_b")),
      lep_d0_cut_e_(ps.getUntrackedParameter<double>("lep_d0_cut_e")),
      lep_dz_cut_e_(ps.getUntrackedParameter<double>("lep_dz_cut_e")),
      ptbins_(ps.getParameter<std::vector<double> >("ptbins")),
      htbins_(ps.getParameter<std::vector<double> >("htbins")),
 
      nbins_eta_(ps.getUntrackedParameter<int>("nbins_eta")),
      nbins_phi_(ps.getUntrackedParameter<int>("nbins_phi")),
      nbins_npv_(ps.getUntrackedParameter<int>("nbins_npv")),
      etabins_min_(ps.getUntrackedParameter<double>("etabins_min")),
      etabins_max_(ps.getUntrackedParameter<double>("etabins_max")),
      phibins_min_(ps.getUntrackedParameter<double>("phibins_min")),
      phibins_max_(ps.getUntrackedParameter<double>("phibins_max")),
      npvbins_min_(ps.getUntrackedParameter<double>("npvbins_min")),
      npvbins_max_(ps.getUntrackedParameter<double>("npvbins_max")),
 
      h_pfHTTurnOn_num_(nullptr),
      h_pfHTTurnOn_den_(nullptr),
      h_lepPtTurnOn_num_(nullptr),
      h_lepPtTurnOn_den_(nullptr),
      h_lepEtaTurnOn_num_(nullptr),
      h_lepEtaTurnOn_den_(nullptr),
      h_lepPhiTurnOn_num_(nullptr),
      h_lepPhiTurnOn_den_(nullptr),
      h_NPVTurnOn_num_(nullptr),
      h_NPVTurnOn_den_(nullptr) {
  edm::LogInfo("LepHTMonitor") << "Constructor LepHTMonitor::LepHTMonitor\n";
}

~LepHTMonitor()

LepHTMonitor::~LepHTMonitor ( ) throw ( )

override

Definition at line 302 of file LepHTMonitor.cc.

                                    {
  if (num_genTriggerEventFlag_) {
    num_genTriggerEventFlag_.reset();
  }
  if (den_lep_genTriggerEventFlag_) {
    den_lep_genTriggerEventFlag_.reset();
  }
  if (den_HT_genTriggerEventFlag_) {
    den_HT_genTriggerEventFlag_.reset();
  }
}

References den_HT_genTriggerEventFlag_, den_lep_genTriggerEventFlag_, and num_genTriggerEventFlag_.

Member Function Documentation

analyze()

void LepHTMonitor::analyze ( const edm::Event &  e, const edm::EventSetup &  eSetup  )

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 416 of file LepHTMonitor.cc.

                                                                         {
  // if valid HLT paths are required,
  // analyze event only if all paths are valid
  if (requireValidHLTPaths_ and (not hltPathsAreValid_)) {
    return;
  }
 
  // Find whether main and auxilliary triggers fired
  bool hasFired = false;
  bool hasFiredAuxiliary = false;
  bool hasFiredLeptonAuxiliary = false;
  if (den_lep_genTriggerEventFlag_->on() && den_lep_genTriggerEventFlag_->accept(e, eSetup))
    hasFiredLeptonAuxiliary = true;
  if (den_HT_genTriggerEventFlag_->on() && den_HT_genTriggerEventFlag_->accept(e, eSetup))
    hasFiredAuxiliary = true;
  if (num_genTriggerEventFlag_->on() && num_genTriggerEventFlag_->accept(e, eSetup))
    hasFired = true;
 
  if (!(hasFiredAuxiliary || hasFiredLeptonAuxiliary))
    return;
  int npv = 0;
  //Vertex
  edm::Handle<reco::VertexCollection> VertexCollection;
  if (not theVertexCollectionTag_.label().empty()) {
    e.getByToken(theVertexCollection_, VertexCollection);
    if (!VertexCollection.isValid()) {
      edm::LogWarning("LepHTMonitor") << "Invalid VertexCollection: " << theVertexCollectionTag_.label() << '\n';
    } else
      npv = VertexCollection->size();
  }
 
  //Get electron ID map
  edm::Handle<edm::ValueMap<bool> > ele_id_decisions;
  if (not theElectronVIDTag_.label().empty()) {
    e.getByToken(theElectronVIDMap_, ele_id_decisions);
    if (!ele_id_decisions.isValid()) {
      edm::LogWarning("LepHTMonitor") << "Invalid Electron VID map: " << theElectronVIDTag_.label() << '\n';
    }
  }
 
  //Conversions
  edm::Handle<reco::ConversionCollection> ConversionCollection;
  if (not theConversionCollectionTag_.label().empty()) {
    e.getByToken(theConversionCollection_, ConversionCollection);
    if (!ConversionCollection.isValid()) {
      edm::LogWarning("LepHTMonitor") << "Invalid ConversionCollection: " << theConversionCollectionTag_.label()
                                      << '\n';
    }
  }
 
  //Beam Spot
  edm::Handle<reco::BeamSpot> BeamSpot;
  if (not theBeamSpotTag_.label().empty()) {
    e.getByToken(theBeamSpot_, BeamSpot);
    if (!BeamSpot.isValid()) {
      edm::LogWarning("LepHTMonitor") << "Invalid BeamSpot: " << theBeamSpotTag_.label() << '\n';
    }
  }
 
  //MET
  edm::Handle<reco::PFMETCollection> pfMETCollection;
  if (not thePfMETTag_.label().empty()) {
    e.getByToken(thePfMETCollection_, pfMETCollection);
    if (!pfMETCollection.isValid()) {
      edm::LogWarning("LepHTMonitor") << "Invalid PFMETCollection: " << thePfMETTag_.label() << '\n';
    }
  }
 
  //Jets
  edm::Handle<reco::PFJetCollection> pfJetCollection;
  if (not thePfJetTag_.label().empty()) {
    e.getByToken(thePfJetCollection_, pfJetCollection);
    if (!pfJetCollection.isValid()) {
      edm::LogWarning("LepHTMonitor") << "Invalid PFJetCollection: " << thePfJetTag_.label() << '\n';
    }
  }
 
  //Electron
  edm::Handle<edm::View<reco::GsfElectron> > ElectronCollection;
  if (not theElectronTag_.label().empty()) {
    e.getByToken(theElectronCollection_, ElectronCollection);
    if (!ElectronCollection.isValid()) {
      edm::LogWarning("LepHTMonitor") << "Invalid GsfElectronCollection: " << theElectronTag_.label() << '\n';
    }
  }
 
  //Muon
  edm::Handle<reco::MuonCollection> MuonCollection;
  if (not theMuonTag_.label().empty()) {
    e.getByToken(theMuonCollection_, MuonCollection);
    if (!MuonCollection.isValid()) {
      edm::LogWarning("LepHTMonitor") << "Invalid MuonCollection: " << theMuonTag_.label() << '\n';
    }
  }
 
  //Get offline HT
  double pfHT = -1.0;
  if (pfJetCollection.isValid()) {
    pfHT = 0.0;
    for (auto const &pfjet : *pfJetCollection) {
      if (pfjet.pt() < jetPtCut_)
        continue;
      if (std::abs(pfjet.eta()) > jetEtaCut_)
        continue;
      pfHT += pfjet.pt();
    }
  }
 
  //Get offline MET
  double pfMET = -1.0;
  if (pfMETCollection.isValid() && !pfMETCollection->empty()) {
    pfMET = pfMETCollection->front().et();
  }
 
  //Find offline leptons and keep track of pt,eta of leading and trailing leptons
  double lep_max_pt = -1.0;
  double lep_eta = 0;
  double lep_phi = 0;
  double trailing_ele_eta = 0;
  double trailing_ele_phi = 0;
  double trailing_mu_eta = 0;
  double trailing_mu_phi = 0;
  double min_ele_pt = -1.0;
  double min_mu_pt = -1.0;
  int nels = 0;
  int nmus = 0;
  if (VertexCollection.isValid() && !VertexCollection->empty()) {  //for quality checks
    //Try to find a reco electron
    if (ElectronCollection.isValid() && ConversionCollection.isValid() && BeamSpot.isValid() &&
        ele_id_decisions.isValid()) {
      size_t index = 0;
      for (auto const &electron : *ElectronCollection) {
        const auto el = ElectronCollection->ptrAt(index);
        bool pass_id = (*ele_id_decisions)[el];
        if (isGood(electron,
                   VertexCollection->front().position(),
                   BeamSpot->position(),
                   ConversionCollection,
                   pass_id,
                   lep_counting_threshold_,
                   lep_iso_cut_,
                   lep_eta_cut_,
                   lep_d0_cut_b_,
                   lep_dz_cut_b_,
                   lep_d0_cut_e_,
                   lep_dz_cut_e_)) {
          if (electron.pt() > lep_max_pt) {
            lep_max_pt = electron.pt();
            lep_eta = electron.eta();
            lep_phi = electron.phi();
          }
          if (electron.pt() < min_ele_pt || min_ele_pt < 0) {
            min_ele_pt = electron.pt();
            trailing_ele_eta = electron.eta();
            trailing_ele_phi = electron.phi();
          }
          nels++;
        }
        index++;
      }
    }
 
    //Try to find a reco muon
    if (MuonCollection.isValid()) {
      for (auto const &muon : *MuonCollection) {
        if (isGood(muon,
                   VertexCollection->front(),
                   lep_counting_threshold_,
                   lep_iso_cut_,
                   lep_eta_cut_,
                   lep_d0_cut_b_,
                   lep_dz_cut_b_,
                   muonIDlevel_)) {
          if (muon.pt() > lep_max_pt) {
            lep_max_pt = muon.pt();
            lep_eta = muon.eta();
            lep_phi = muon.phi();
          }
          if (muon.pt() < min_mu_pt || min_mu_pt < 0) {
            min_mu_pt = muon.pt();
            trailing_mu_eta = muon.eta();
            trailing_mu_phi = muon.phi();
          }
          nmus++;
        }
      }
    }
  }
 
  //Fill single lepton triggers with leading lepton pT
  float lep_pt = lep_max_pt;
 
  //For dilepton triggers, use trailing rather than leading lepton
  if (nmusCut_ >= 2) {
    lep_pt = min_mu_pt;
    lep_eta = trailing_mu_eta;
    lep_phi = trailing_mu_phi;
  }
  if (nelsCut_ >= 2) {
    lep_pt = min_ele_pt;
    lep_eta = trailing_ele_eta;
    lep_phi = trailing_ele_phi;
  }
  if (nelsCut_ >= 1 && nmusCut_ >= 1) {
    if (min_ele_pt < min_mu_pt) {
      lep_pt = min_ele_pt;
      lep_eta = trailing_ele_eta;
      lep_phi = trailing_ele_phi;
    } else {
      lep_pt = min_mu_pt;
      lep_eta = trailing_mu_eta;
      lep_phi = trailing_mu_phi;
    }
  }
 
  const bool nleps_cut = nels >= nelsCut_ && nmus >= nmusCut_;
  bool lep_plateau = lep_pt > lep_pt_plateau_ || lep_pt_plateau_ < 0.0;
 
  //Fill lepton pT and eta histograms
  if (hasFiredLeptonAuxiliary || !e.isRealData()) {
    if (nleps_cut && (pfMET > metCut_ || metCut_ < 0.0) && (pfHT > htCut_ || htCut_ < 0.0)) {
      if (h_lepPtTurnOn_den_) {
        if (lep_pt > ptbins_.back())
          lep_pt = ptbins_.back() - 1;  //Overflow protection
        h_lepPtTurnOn_den_->Fill(lep_pt);
      }
      if (h_lepPtTurnOn_num_ && hasFired)
        h_lepPtTurnOn_num_->Fill(lep_pt);
 
      if (lep_plateau) {
        //Fill Eta and Phi histograms for leptons above pT threshold
        if (h_lepEtaTurnOn_den_)
          h_lepEtaTurnOn_den_->Fill(lep_eta);
        if (h_lepEtaTurnOn_num_ && hasFired)
          h_lepEtaTurnOn_num_->Fill(lep_eta);
        if (h_lepPhiTurnOn_den_)
          h_lepPhiTurnOn_den_->Fill(lep_phi);
        if (h_lepPhiTurnOn_num_ && hasFired)
          h_lepPhiTurnOn_num_->Fill(lep_phi);
        if (h_lepEtaPhiTurnOn_den_)
          h_lepEtaPhiTurnOn_den_->Fill(lep_eta, lep_phi);
        if (h_lepEtaPhiTurnOn_num_ && hasFired)
          h_lepEtaPhiTurnOn_num_->Fill(lep_eta, lep_phi);
 
        //Fill NPV histograms
        if (h_NPVTurnOn_den_)
          h_NPVTurnOn_den_->Fill(npv);
        if (h_NPVTurnOn_num_ && hasFired)
          h_NPVTurnOn_num_->Fill(npv);
      }
    }
  }
 
  //Fill HT turn-on histograms
  if (hasFiredAuxiliary || !e.isRealData()) {
    if (nleps_cut && lep_plateau) {
      if (h_pfHTTurnOn_den_) {
        if (pfHT > htbins_.back())
          pfHT = htbins_.back() - 1;  //Overflow protection
        h_pfHTTurnOn_den_->Fill(pfHT);
      }
      if (h_pfHTTurnOn_num_ && hasFired)
        h_pfHTTurnOn_num_->Fill(pfHT);
    }
  }
}

References funct::abs(), align::BeamSpot, den_HT_genTriggerEventFlag_, den_lep_genTriggerEventFlag_, MillePedeFileConverter_cfg::e, metsig::electron, dqm::impl::MonitorElement::Fill(), h_lepEtaPhiTurnOn_den_, h_lepEtaPhiTurnOn_num_, h_lepEtaTurnOn_den_, h_lepEtaTurnOn_num_, h_lepPhiTurnOn_den_, h_lepPhiTurnOn_num_, h_lepPtTurnOn_den_, h_lepPtTurnOn_num_, h_NPVTurnOn_den_, h_NPVTurnOn_num_, h_pfHTTurnOn_den_, h_pfHTTurnOn_num_, hltPathsAreValid_, htbins_, htCut_, edm::HandleBase::isValid(), jetEtaCut_, jetPtCut_, edm::InputTag::label(), lep_counting_threshold_, lep_d0_cut_b_, lep_d0_cut_e_, lep_dz_cut_b_, lep_dz_cut_e_, lep_eta_cut_, lep_iso_cut_, lep_pt_plateau_, metCut_, muonIDlevel_, LepHTMonitor_cff::nels, nelsCut_, LepHTMonitor_cff::nmus, nmusCut_, num_genTriggerEventFlag_, ExoticaDQM_cfi::pfJetCollection, pfMETBenchmark_cfi::pfMET, B2GDQM_cfi::pfMETCollection, ptbins_, requireValidHLTPaths_, theBeamSpot_, theBeamSpotTag_, theConversionCollection_, theConversionCollectionTag_, theElectronCollection_, theElectronTag_, theElectronVIDMap_, theElectronVIDTag_, theMuonCollection_, theMuonTag_, thePfJetCollection_, thePfJetTag_, thePfMETCollection_, thePfMETTag_, theVertexCollection_, and theVertexCollectionTag_.

bookHistograms()

void LepHTMonitor::bookHistograms ( DQMStore::IBooker &  ibooker, const edm::Run &  iRun, const edm::EventSetup &  iSetup  )

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 314 of file LepHTMonitor.cc.

                                                                                                           {
  // Initialize trigger flags
  if (num_genTriggerEventFlag_ && num_genTriggerEventFlag_->on()) {
    num_genTriggerEventFlag_->initRun(iRun, iSetup);
  }
  if (den_lep_genTriggerEventFlag_ && den_lep_genTriggerEventFlag_->on()) {
    den_lep_genTriggerEventFlag_->initRun(iRun, iSetup);
  }
  if (den_HT_genTriggerEventFlag_ && den_HT_genTriggerEventFlag_->on()) {
    den_HT_genTriggerEventFlag_->initRun(iRun, iSetup);
  }
 
  // check if every HLT path specified in numerator and denominator has a valid match in the HLT Menu
  hltPathsAreValid_ = ((num_genTriggerEventFlag_ && num_genTriggerEventFlag_->on() &&
                        num_genTriggerEventFlag_->allHLTPathsAreValid()) &&
                       ((den_lep_genTriggerEventFlag_ && den_lep_genTriggerEventFlag_->on() &&
                         den_lep_genTriggerEventFlag_->allHLTPathsAreValid()) ||
                        (den_HT_genTriggerEventFlag_ && den_HT_genTriggerEventFlag_->on() &&
                         den_HT_genTriggerEventFlag_->allHLTPathsAreValid())));
 
  // if valid HLT paths are required,
  // create DQM outputs only if all paths are valid
  if (requireValidHLTPaths_ and (not hltPathsAreValid_)) {
    return;
  }
 
  // book at beginRun
  ibooker.cd();
  ibooker.setCurrentFolder("HLT/SUSY/LepHT/" + folderName_);
 
  bool is_mu = false;
  bool is_ele = false;
  if (theElectronTag_.label().empty() and not theMuonTag_.label().empty()) {
    is_mu = true;
  } else if (not theElectronTag_.label().empty() and theMuonTag_.label().empty()) {
    is_ele = true;
  }
 
  // Cosmetic axis names
  std::string lepton = "lepton", Lepton = "Lepton";
  if (is_mu && !is_ele) {
    lepton = "muon";
    Lepton = "Muon";
  } else if (is_ele && !is_mu) {
    lepton = "electron";
    Lepton = "Electron";
  }
 
  //Convert to vfloat for picky TH1F constructor
  vector<float> f_ptbins;
  for (double ptbin : ptbins_)
    f_ptbins.push_back(static_cast<float>(ptbin));
  vector<float> f_htbins;
  for (double htbin : htbins_)
    f_htbins.push_back(static_cast<float>(htbin));
 
  //num and den hists to be divided in harvesting step to make turn on curves
  h_pfHTTurnOn_num_ =
      ibooker.book1D("pfHTTurnOn_num", "Numerator;Offline H_{T} [GeV];", f_htbins.size() - 1, f_htbins.data());
  h_pfHTTurnOn_den_ =
      ibooker.book1D("pfHTTurnOn_den", "Denominator;Offline H_{T} [GeV];", f_htbins.size() - 1, f_htbins.data());
 
  h_lepPtTurnOn_num_ = ibooker.book1D("lepPtTurnOn_num",
                                      ("Numerator;Offline " + lepton + " p_{T} [GeV];").c_str(),
                                      f_ptbins.size() - 1,
                                      f_ptbins.data());
  h_lepPtTurnOn_den_ = ibooker.book1D("lepPtTurnOn_den",
                                      ("Denominator;Offline " + lepton + " p_{T} [GeV];").c_str(),
                                      f_ptbins.size() - 1,
                                      f_ptbins.data());
  h_lepEtaTurnOn_num_ =
      ibooker.book1D("lepEtaTurnOn_num", "Numerator;Offline lepton #eta;", nbins_eta_, etabins_min_, etabins_max_);
  h_lepEtaTurnOn_den_ =
      ibooker.book1D("lepEtaTurnOn_den", "Denominator;Offline lepton #eta;", nbins_eta_, etabins_min_, etabins_max_);
  h_lepPhiTurnOn_num_ =
      ibooker.book1D("lepPhiTurnOn_num", "Numerator;Offline lepton #phi;", nbins_phi_, phibins_min_, phibins_max_);
  h_lepPhiTurnOn_den_ =
      ibooker.book1D("lepPhiTurnOn_den", "Denominator;Offline lepton #phi;", nbins_phi_, phibins_min_, phibins_max_);
 
  h_lepEtaPhiTurnOn_num_ = ibooker.book2D("lepEtaPhiTurnOn_num",
                                          "Numerator;Offline lepton #eta;Offline lepton #phi;",
                                          nbins_eta_ / 2,
                                          etabins_min_,
                                          etabins_max_,
                                          nbins_phi_ / 2,
                                          phibins_min_,
                                          phibins_max_);
  h_lepEtaPhiTurnOn_den_ = ibooker.book2D("lepEtaPhiTurnOn_den",
                                          "Denominator;Offline lepton #eta;Offline lepton #phi;",
                                          nbins_eta_ / 2,
                                          etabins_min_,
                                          etabins_max_,
                                          nbins_phi_ / 2,
                                          phibins_min_,
                                          phibins_max_);
 
  h_NPVTurnOn_num_ = ibooker.book1D("NPVTurnOn_num", "Numerator;N_{PV};", nbins_npv_, npvbins_min_, npvbins_max_);
  h_NPVTurnOn_den_ = ibooker.book1D("NPVTurnOn_den", "Denominator;N_{PV};", nbins_npv_, npvbins_min_, npvbins_max_);
 
  ibooker.cd();
}

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::NavigatorBase::cd(), den_HT_genTriggerEventFlag_, den_lep_genTriggerEventFlag_, etabins_max_, etabins_min_, folderName_, h_lepEtaPhiTurnOn_den_, h_lepEtaPhiTurnOn_num_, h_lepEtaTurnOn_den_, h_lepEtaTurnOn_num_, h_lepPhiTurnOn_den_, h_lepPhiTurnOn_num_, h_lepPtTurnOn_den_, h_lepPtTurnOn_num_, h_NPVTurnOn_den_, h_NPVTurnOn_num_, h_pfHTTurnOn_den_, h_pfHTTurnOn_num_, hltPathsAreValid_, htbins_, edm::InputTag::label(), nbins_eta_, nbins_npv_, nbins_phi_, npvbins_max_, npvbins_min_, num_genTriggerEventFlag_, phibins_max_, phibins_min_, ptbins_, requireValidHLTPaths_, dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, theElectronTag_, and theMuonTag_.

Member Data Documentation

den_HT_genTriggerEventFlag_

std::unique_ptr<GenericTriggerEventFlag> LepHTMonitor::den_HT_genTriggerEventFlag_

private

Definition at line 70 of file LepHTMonitor.cc.

Referenced by analyze(), bookHistograms(), and ~LepHTMonitor().

den_lep_genTriggerEventFlag_

std::unique_ptr<GenericTriggerEventFlag> LepHTMonitor::den_lep_genTriggerEventFlag_

private

Definition at line 69 of file LepHTMonitor.cc.

Referenced by analyze(), bookHistograms(), and ~LepHTMonitor().

etabins_max_

float LepHTMonitor::etabins_max_

private

Definition at line 100 of file LepHTMonitor.cc.

Referenced by bookHistograms().

etabins_min_

float LepHTMonitor::etabins_min_

private

Definition at line 99 of file LepHTMonitor.cc.

Referenced by bookHistograms().

folderName_

const std::string LepHTMonitor::folderName_

private

Definition at line 72 of file LepHTMonitor.cc.

Referenced by bookHistograms().

h_lepEtaPhiTurnOn_den_

MonitorElement* LepHTMonitor::h_lepEtaPhiTurnOn_den_

private

Definition at line 116 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

h_lepEtaPhiTurnOn_num_

MonitorElement* LepHTMonitor::h_lepEtaPhiTurnOn_num_

private

Definition at line 115 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

h_lepEtaTurnOn_den_

MonitorElement* LepHTMonitor::h_lepEtaTurnOn_den_

private

Definition at line 112 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

h_lepEtaTurnOn_num_

MonitorElement* LepHTMonitor::h_lepEtaTurnOn_num_

private

Definition at line 111 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

h_lepPhiTurnOn_den_

MonitorElement* LepHTMonitor::h_lepPhiTurnOn_den_

private

Definition at line 114 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

h_lepPhiTurnOn_num_

MonitorElement* LepHTMonitor::h_lepPhiTurnOn_num_

private

Definition at line 113 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

h_lepPtTurnOn_den_

MonitorElement* LepHTMonitor::h_lepPtTurnOn_den_

private

Definition at line 110 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

h_lepPtTurnOn_num_

MonitorElement* LepHTMonitor::h_lepPtTurnOn_num_

private

Definition at line 109 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

h_NPVTurnOn_den_

MonitorElement* LepHTMonitor::h_NPVTurnOn_den_

private

Definition at line 118 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

h_NPVTurnOn_num_

MonitorElement* LepHTMonitor::h_NPVTurnOn_num_

private

Definition at line 117 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

h_pfHTTurnOn_den_

MonitorElement* LepHTMonitor::h_pfHTTurnOn_den_

private

Definition at line 108 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

h_pfHTTurnOn_num_

MonitorElement* LepHTMonitor::h_pfHTTurnOn_num_

private

Definition at line 107 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

hltPathsAreValid_

bool LepHTMonitor::hltPathsAreValid_

private

Definition at line 75 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

htbins_

std::vector<double> LepHTMonitor::htbins_

private

Definition at line 95 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

htCut_

double LepHTMonitor::htCut_

private

Definition at line 82 of file LepHTMonitor.cc.

Referenced by analyze().

jetEtaCut_

double LepHTMonitor::jetEtaCut_

private

Definition at line 80 of file LepHTMonitor.cc.

Referenced by analyze().

jetPtCut_

double LepHTMonitor::jetPtCut_

private

Definition at line 79 of file LepHTMonitor.cc.

Referenced by analyze().

lep_counting_threshold_

double LepHTMonitor::lep_counting_threshold_

private

Definition at line 86 of file LepHTMonitor.cc.

Referenced by analyze().

lep_d0_cut_b_

double LepHTMonitor::lep_d0_cut_b_

private

Definition at line 89 of file LepHTMonitor.cc.

Referenced by analyze().

lep_d0_cut_e_

double LepHTMonitor::lep_d0_cut_e_

private

Definition at line 91 of file LepHTMonitor.cc.

Referenced by analyze().

lep_dz_cut_b_

double LepHTMonitor::lep_dz_cut_b_

private

Definition at line 90 of file LepHTMonitor.cc.

Referenced by analyze().

lep_dz_cut_e_

double LepHTMonitor::lep_dz_cut_e_

private

Definition at line 92 of file LepHTMonitor.cc.

Referenced by analyze().

lep_eta_cut_

double LepHTMonitor::lep_eta_cut_

private

Definition at line 88 of file LepHTMonitor.cc.

Referenced by analyze().

lep_iso_cut_

double LepHTMonitor::lep_iso_cut_

private

Definition at line 87 of file LepHTMonitor.cc.

Referenced by analyze().

lep_pt_plateau_

double LepHTMonitor::lep_pt_plateau_

private

Definition at line 85 of file LepHTMonitor.cc.

Referenced by analyze().

metCut_

double LepHTMonitor::metCut_

private

Definition at line 81 of file LepHTMonitor.cc.

Referenced by analyze().

muonIDlevel_

int LepHTMonitor::muonIDlevel_

private

Definition at line 77 of file LepHTMonitor.cc.

Referenced by analyze().

nbins_eta_

int LepHTMonitor::nbins_eta_

private

Definition at line 96 of file LepHTMonitor.cc.

Referenced by bookHistograms().

nbins_npv_

int LepHTMonitor::nbins_npv_

private

Definition at line 98 of file LepHTMonitor.cc.

Referenced by bookHistograms().

nbins_phi_

int LepHTMonitor::nbins_phi_

private

Definition at line 97 of file LepHTMonitor.cc.

Referenced by bookHistograms().

nelsCut_

double LepHTMonitor::nelsCut_

private

Definition at line 84 of file LepHTMonitor.cc.

Referenced by analyze().

nmusCut_

double LepHTMonitor::nmusCut_

private

Definition at line 83 of file LepHTMonitor.cc.

Referenced by analyze().

npvbins_max_

float LepHTMonitor::npvbins_max_

private

Definition at line 104 of file LepHTMonitor.cc.

Referenced by bookHistograms().

npvbins_min_

float LepHTMonitor::npvbins_min_

private

Definition at line 103 of file LepHTMonitor.cc.

Referenced by bookHistograms().

num_genTriggerEventFlag_

std::unique_ptr<GenericTriggerEventFlag> LepHTMonitor::num_genTriggerEventFlag_

private

Definition at line 68 of file LepHTMonitor.cc.

Referenced by analyze(), bookHistograms(), and ~LepHTMonitor().

phibins_max_

float LepHTMonitor::phibins_max_

private

Definition at line 102 of file LepHTMonitor.cc.

Referenced by bookHistograms().

phibins_min_

float LepHTMonitor::phibins_min_

private

Definition at line 101 of file LepHTMonitor.cc.

Referenced by bookHistograms().

ptbins_

std::vector<double> LepHTMonitor::ptbins_

private

Definition at line 94 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

requireValidHLTPaths_

const bool LepHTMonitor::requireValidHLTPaths_

private

Definition at line 74 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

theBeamSpot_

edm::EDGetTokenT<reco::BeamSpot> LepHTMonitor::theBeamSpot_

private

Definition at line 66 of file LepHTMonitor.cc.

Referenced by analyze().

theBeamSpotTag_

edm::InputTag LepHTMonitor::theBeamSpotTag_

private

Definition at line 65 of file LepHTMonitor.cc.

Referenced by analyze().

theConversionCollection_

edm::EDGetTokenT<reco::ConversionCollection> LepHTMonitor::theConversionCollection_

private

Definition at line 64 of file LepHTMonitor.cc.

Referenced by analyze().

theConversionCollectionTag_

edm::InputTag LepHTMonitor::theConversionCollectionTag_

private

Definition at line 63 of file LepHTMonitor.cc.

Referenced by analyze().

theElectronCollection_

edm::EDGetTokenT<edm::View<reco::GsfElectron> > LepHTMonitor::theElectronCollection_

private

Definition at line 50 of file LepHTMonitor.cc.

Referenced by analyze().

theElectronTag_

edm::InputTag LepHTMonitor::theElectronTag_

private

Definition at line 49 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

theElectronVIDMap_

edm::EDGetTokenT<edm::ValueMap<bool> > LepHTMonitor::theElectronVIDMap_

private

Definition at line 52 of file LepHTMonitor.cc.

Referenced by analyze().

theElectronVIDTag_

edm::InputTag LepHTMonitor::theElectronVIDTag_

private

Definition at line 51 of file LepHTMonitor.cc.

Referenced by analyze().

theJetTagCollection_

edm::EDGetTokenT<reco::JetTagCollection> LepHTMonitor::theJetTagCollection_

private

Definition at line 60 of file LepHTMonitor.cc.

theJetTagTag_

edm::InputTag LepHTMonitor::theJetTagTag_

private

Definition at line 59 of file LepHTMonitor.cc.

theMuonCollection_

edm::EDGetTokenT<reco::MuonCollection> LepHTMonitor::theMuonCollection_

private

Definition at line 54 of file LepHTMonitor.cc.

Referenced by analyze().

theMuonTag_

edm::InputTag LepHTMonitor::theMuonTag_

private

Definition at line 53 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

thePfJetCollection_

edm::EDGetTokenT<reco::PFJetCollection> LepHTMonitor::thePfJetCollection_

private

Definition at line 58 of file LepHTMonitor.cc.

Referenced by analyze().

thePfJetTag_

edm::InputTag LepHTMonitor::thePfJetTag_

private

Definition at line 57 of file LepHTMonitor.cc.

Referenced by analyze().

thePfMETCollection_

edm::EDGetTokenT<reco::PFMETCollection> LepHTMonitor::thePfMETCollection_

private

Definition at line 56 of file LepHTMonitor.cc.

Referenced by analyze().

thePfMETTag_

edm::InputTag LepHTMonitor::thePfMETTag_

private

Definition at line 55 of file LepHTMonitor.cc.

Referenced by analyze().

theVertexCollection_

edm::EDGetTokenT<reco::VertexCollection> LepHTMonitor::theVertexCollection_

private

Definition at line 62 of file LepHTMonitor.cc.

Referenced by analyze().

theVertexCollectionTag_

edm::InputTag LepHTMonitor::theVertexCollectionTag_

private

Definition at line 61 of file LepHTMonitor.cc.

Referenced by analyze().