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 >
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::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
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete
Public Member Functions inherited from TriggerDQMBase
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 std::vector< double > &binningX, 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	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	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 35 of file LepHTMonitor.cc.

Member Typedef Documentation

DQMStore

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

Definition at line 38 of file LepHTMonitor.cc.

MonitorElement

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

Definition at line 37 of file LepHTMonitor.cc.

Constructor & Destructor Documentation

LepHTMonitor()

LepHTMonitor::LepHTMonitor

(

const edm::ParameterSet &

ps

)

Definition at line 228 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 301 of file LepHTMonitor.cc.

References den_HT_genTriggerEventFlag_, den_lep_genTriggerEventFlag_, and num_genTriggerEventFlag_.

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

Member Function Documentation

analyze()

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 417 of file LepHTMonitor.cc.

References funct::abs(), align::BeamSpot, den_HT_genTriggerEventFlag_, den_lep_genTriggerEventFlag_, MillePedeFileConverter_cfg::e, HPSPFTauProducerPuppi_cfi::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, pfMET_cfi::pfMET, B2GDQM_cfi::pfMETCollection, ptbins_, requireValidHLTPaths_, theBeamSpot_, theBeamSpotTag_, theConversionCollection_, theConversionCollectionTag_, theElectronCollection_, theElectronTag_, theElectronVIDMap_, theElectronVIDTag_, theMuonCollection_, theMuonTag_, thePfJetCollection_, thePfJetTag_, thePfMETCollection_, thePfMETTag_, theVertexCollection_, and theVertexCollectionTag_.

                                                                         {
  // 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);
    }
  }
}

bookHistograms()

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 313 of file LepHTMonitor.cc.

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_.

                                                                                                           {
  // 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;
  f_ptbins.reserve(ptbins_.size());
  for (double ptbin : ptbins_)
    f_ptbins.push_back(static_cast<float>(ptbin));
  vector<float> f_htbins;
  f_htbins.reserve(htbins_.size());
  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();
}

Member Data Documentation

den_HT_genTriggerEventFlag_

std::unique_ptr<GenericTriggerEventFlag> LepHTMonitor::den_HT_genTriggerEventFlag_

private

Definition at line 69 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 68 of file LepHTMonitor.cc.

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

etabins_max_

float LepHTMonitor::etabins_max_

private

Definition at line 99 of file LepHTMonitor.cc.

Referenced by bookHistograms().

etabins_min_

float LepHTMonitor::etabins_min_

private

Definition at line 98 of file LepHTMonitor.cc.

Referenced by bookHistograms().

folderName_

const std::string LepHTMonitor::folderName_

private

Definition at line 71 of file LepHTMonitor.cc.

Referenced by bookHistograms().

h_lepEtaPhiTurnOn_den_

MonitorElement* LepHTMonitor::h_lepEtaPhiTurnOn_den_

private

Definition at line 115 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

h_lepEtaPhiTurnOn_num_

MonitorElement* LepHTMonitor::h_lepEtaPhiTurnOn_num_

private

Definition at line 114 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

h_lepEtaTurnOn_den_

MonitorElement* LepHTMonitor::h_lepEtaTurnOn_den_

private

Definition at line 111 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

h_lepEtaTurnOn_num_

MonitorElement* LepHTMonitor::h_lepEtaTurnOn_num_

private

Definition at line 110 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

h_lepPhiTurnOn_den_

MonitorElement* LepHTMonitor::h_lepPhiTurnOn_den_

private

Definition at line 113 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

h_lepPhiTurnOn_num_

MonitorElement* LepHTMonitor::h_lepPhiTurnOn_num_

private

Definition at line 112 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

h_lepPtTurnOn_den_

MonitorElement* LepHTMonitor::h_lepPtTurnOn_den_

private

Definition at line 109 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

h_lepPtTurnOn_num_

MonitorElement* LepHTMonitor::h_lepPtTurnOn_num_

private

Definition at line 108 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

h_NPVTurnOn_den_

MonitorElement* LepHTMonitor::h_NPVTurnOn_den_

private

Definition at line 117 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

h_NPVTurnOn_num_

MonitorElement* LepHTMonitor::h_NPVTurnOn_num_

private

Definition at line 116 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

h_pfHTTurnOn_den_

MonitorElement* LepHTMonitor::h_pfHTTurnOn_den_

private

Definition at line 107 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

h_pfHTTurnOn_num_

MonitorElement* LepHTMonitor::h_pfHTTurnOn_num_

private

Definition at line 106 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

hltPathsAreValid_

bool LepHTMonitor::hltPathsAreValid_

private

Definition at line 74 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

htbins_

std::vector<double> LepHTMonitor::htbins_

private

Definition at line 94 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

htCut_

double LepHTMonitor::htCut_

private

Definition at line 81 of file LepHTMonitor.cc.

Referenced by analyze().

jetEtaCut_

double LepHTMonitor::jetEtaCut_

private

Definition at line 79 of file LepHTMonitor.cc.

Referenced by analyze().

jetPtCut_

double LepHTMonitor::jetPtCut_

private

Definition at line 78 of file LepHTMonitor.cc.

Referenced by analyze().

lep_counting_threshold_

double LepHTMonitor::lep_counting_threshold_

private

Definition at line 85 of file LepHTMonitor.cc.

Referenced by analyze().

lep_d0_cut_b_

double LepHTMonitor::lep_d0_cut_b_

private

Definition at line 88 of file LepHTMonitor.cc.

Referenced by analyze().

lep_d0_cut_e_

double LepHTMonitor::lep_d0_cut_e_

private

Definition at line 90 of file LepHTMonitor.cc.

Referenced by analyze().

lep_dz_cut_b_

double LepHTMonitor::lep_dz_cut_b_

private

Definition at line 89 of file LepHTMonitor.cc.

Referenced by analyze().

lep_dz_cut_e_

double LepHTMonitor::lep_dz_cut_e_

private

Definition at line 91 of file LepHTMonitor.cc.

Referenced by analyze().

lep_eta_cut_

double LepHTMonitor::lep_eta_cut_

private

Definition at line 87 of file LepHTMonitor.cc.

Referenced by analyze().

lep_iso_cut_

double LepHTMonitor::lep_iso_cut_

private

Definition at line 86 of file LepHTMonitor.cc.

Referenced by analyze().

lep_pt_plateau_

double LepHTMonitor::lep_pt_plateau_

private

Definition at line 84 of file LepHTMonitor.cc.

Referenced by analyze().

metCut_

double LepHTMonitor::metCut_

private

Definition at line 80 of file LepHTMonitor.cc.

Referenced by analyze().

muonIDlevel_

int LepHTMonitor::muonIDlevel_

private

Definition at line 76 of file LepHTMonitor.cc.

Referenced by analyze().

nbins_eta_

int LepHTMonitor::nbins_eta_

private

Definition at line 95 of file LepHTMonitor.cc.

Referenced by bookHistograms().

nbins_npv_

int LepHTMonitor::nbins_npv_

private

Definition at line 97 of file LepHTMonitor.cc.

Referenced by bookHistograms().

nbins_phi_

int LepHTMonitor::nbins_phi_

private

Definition at line 96 of file LepHTMonitor.cc.

Referenced by bookHistograms().

nelsCut_

double LepHTMonitor::nelsCut_

private

Definition at line 83 of file LepHTMonitor.cc.

Referenced by analyze().

nmusCut_

double LepHTMonitor::nmusCut_

private

Definition at line 82 of file LepHTMonitor.cc.

Referenced by analyze().

npvbins_max_

float LepHTMonitor::npvbins_max_

private

Definition at line 103 of file LepHTMonitor.cc.

Referenced by bookHistograms().

npvbins_min_

float LepHTMonitor::npvbins_min_

private

Definition at line 102 of file LepHTMonitor.cc.

Referenced by bookHistograms().

num_genTriggerEventFlag_

std::unique_ptr<GenericTriggerEventFlag> LepHTMonitor::num_genTriggerEventFlag_

private

Definition at line 67 of file LepHTMonitor.cc.

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

phibins_max_

float LepHTMonitor::phibins_max_

private

Definition at line 101 of file LepHTMonitor.cc.

Referenced by bookHistograms().

phibins_min_

float LepHTMonitor::phibins_min_

private

Definition at line 100 of file LepHTMonitor.cc.

Referenced by bookHistograms().

ptbins_

std::vector<double> LepHTMonitor::ptbins_

private

Definition at line 93 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

requireValidHLTPaths_

const bool LepHTMonitor::requireValidHLTPaths_

private

Definition at line 73 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

theBeamSpot_

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

private

Definition at line 65 of file LepHTMonitor.cc.

Referenced by analyze().

theBeamSpotTag_

edm::InputTag LepHTMonitor::theBeamSpotTag_

private

Definition at line 64 of file LepHTMonitor.cc.

Referenced by analyze().

theConversionCollection_

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

private

Definition at line 63 of file LepHTMonitor.cc.

Referenced by analyze().

theConversionCollectionTag_

edm::InputTag LepHTMonitor::theConversionCollectionTag_

private

Definition at line 62 of file LepHTMonitor.cc.

Referenced by analyze().

theElectronCollection_

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

private

Definition at line 49 of file LepHTMonitor.cc.

Referenced by analyze().

theElectronTag_

edm::InputTag LepHTMonitor::theElectronTag_

private

Definition at line 48 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

theElectronVIDMap_

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

private

Definition at line 51 of file LepHTMonitor.cc.

Referenced by analyze().

theElectronVIDTag_

edm::InputTag LepHTMonitor::theElectronVIDTag_

private

Definition at line 50 of file LepHTMonitor.cc.

Referenced by analyze().

theJetTagCollection_

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

private

Definition at line 59 of file LepHTMonitor.cc.

theJetTagTag_

edm::InputTag LepHTMonitor::theJetTagTag_

private

Definition at line 58 of file LepHTMonitor.cc.

theMuonCollection_

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

private

Definition at line 53 of file LepHTMonitor.cc.

Referenced by analyze().

theMuonTag_

edm::InputTag LepHTMonitor::theMuonTag_

private

Definition at line 52 of file LepHTMonitor.cc.

Referenced by analyze(), and bookHistograms().

thePfJetCollection_

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

private

Definition at line 57 of file LepHTMonitor.cc.

Referenced by analyze().

thePfJetTag_

edm::InputTag LepHTMonitor::thePfJetTag_

private

Definition at line 56 of file LepHTMonitor.cc.

Referenced by analyze().

thePfMETCollection_

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

private

Definition at line 55 of file LepHTMonitor.cc.

Referenced by analyze().

thePfMETTag_

edm::InputTag LepHTMonitor::thePfMETTag_

private

Definition at line 54 of file LepHTMonitor.cc.

Referenced by analyze().

theVertexCollection_

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

private

Definition at line 61 of file LepHTMonitor.cc.

Referenced by analyze().

theVertexCollectionTag_

edm::InputTag LepHTMonitor::theVertexCollectionTag_

private

Definition at line 60 of file LepHTMonitor.cc.

Referenced by analyze().