LepHTMonitor Class Reference

#include <LepHTMonitor.h>

Inheritance diagram for LepHTMonitor:

Public Member Functions
	LepHTMonitor (const edm::ParameterSet &ps)
	~LepHTMonitor () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
virtual void	analyze (edm::Event const &, edm::EventSetup const &)
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &)=0
virtual void	dqmBeginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
virtual void	dqmEndLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
virtual void	dqmEndRun (edm::Run const &, edm::EventSetup const &)
void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final
void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &, edm::EventSetup const &) final
void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, edm::Accumulator >
	EDProducer ()=default
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDProducerBase () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
	~ProducerBase () noexcept(false) override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

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
void	dqmBeginRun (const edm::Run &run, const edm::EventSetup &e) override
void	dqmEndRun (const edm::Run &run, const edm::EventSetup &eSetup) override
Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Private Attributes
std::unique_ptr< GenericTriggerEventFlag >	den_HT_genTriggerEventFlag_
std::unique_ptr< GenericTriggerEventFlag >	den_lep_genTriggerEventFlag_
float	etabins_max_
float	etabins_min_
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_
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_
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
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType
Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 53 of file LepHTMonitor.h.

Constructor & Destructor Documentation

LepHTMonitor::LepHTMonitor

(

const edm::ParameterSet &

ps

)

Definition at line 132 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")),

      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 ( )

override

Definition at line 204 of file LepHTMonitor.cc.

{ edm::LogInfo("LepHTMonitor") << "Destructor LepHTMonitor::~LepHTMonitor\n"; }

Member Function Documentation

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

overrideprotected

Definition at line 294 of file LepHTMonitor.cc.

References funct::abs(), align::BeamSpot, den_HT_genTriggerEventFlag_, den_lep_genTriggerEventFlag_, metsig::electron, dqm::impl::MonitorElement::Fill(), edm::Event::getByToken(), 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_, htbins_, htCut_, edm::EventBase::isRealData(), 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, reco::BeamSpot::position(), ptbins_, theBeamSpot_, theBeamSpotTag_, theConversionCollection_, theConversionCollectionTag_, theElectronCollection_, theElectronTag_, theElectronVIDMap_, theElectronVIDTag_, theMuonCollection_, theMuonTag_, thePfJetCollection_, thePfJetTag_, thePfMETCollection_, thePfMETTag_, theVertexCollection_, and theVertexCollectionTag_.

                                                                         {
  edm::LogInfo("LepHTMonitor") << "LepHTMonitor::analyze\n";

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

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

overrideprotected

Definition at line 210 of file LepHTMonitor.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book1D(), dqm::dqmstoreimpl::DQMStore::IBooker::book2D(), dqm::dqmstoreimpl::DQMStore::IBooker::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_, htbins_, edm::InputTag::label(), nbins_eta_, nbins_npv_, nbins_phi_, npvbins_max_, npvbins_min_, num_genTriggerEventFlag_, phibins_max_, phibins_min_, ptbins_, dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, theElectronTag_, and theMuonTag_.

                                                                                                           {
  edm::LogInfo("LepHTMonitor") << "LepHTMonitor::bookHistograms\n";
  //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";
  }
  //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);

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

void LepHTMonitor::dqmBeginRun ( const edm::Run &  run, const edm::EventSetup &  e  )

overrideprotected

Definition at line 206 of file LepHTMonitor.cc.

                                                                        {
  edm::LogInfo("LepHTMonitor") << "LepHTMonitor::beginRun\n";
}

void LepHTMonitor::dqmEndRun ( const edm::Run &  run, const edm::EventSetup &  eSetup  )

overrideprotected

Definition at line 557 of file LepHTMonitor.cc.

References DEFINE_FWK_MODULE.

                                                                           {
  edm::LogInfo("LepHTMonitor") << "LepHTMonitor::endRun\n";
}

Member Data Documentation

std::unique_ptr<GenericTriggerEventFlag> LepHTMonitor::den_HT_genTriggerEventFlag_

private

Definition at line 87 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

std::unique_ptr<GenericTriggerEventFlag> LepHTMonitor::den_lep_genTriggerEventFlag_

private

Definition at line 86 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

float LepHTMonitor::etabins_max_

private

Definition at line 114 of file LepHTMonitor.h.

Referenced by bookHistograms().

float LepHTMonitor::etabins_min_

private

Definition at line 113 of file LepHTMonitor.h.

Referenced by bookHistograms().

std::string LepHTMonitor::folderName_

private

Definition at line 89 of file LepHTMonitor.h.

Referenced by bookHistograms().

MonitorElement* LepHTMonitor::h_lepEtaPhiTurnOn_den_

private

Definition at line 130 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* LepHTMonitor::h_lepEtaPhiTurnOn_num_

private

Definition at line 129 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* LepHTMonitor::h_lepEtaTurnOn_den_

private

Definition at line 126 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* LepHTMonitor::h_lepEtaTurnOn_num_

private

Definition at line 125 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* LepHTMonitor::h_lepPhiTurnOn_den_

private

Definition at line 128 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* LepHTMonitor::h_lepPhiTurnOn_num_

private

Definition at line 127 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* LepHTMonitor::h_lepPtTurnOn_den_

private

Definition at line 124 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* LepHTMonitor::h_lepPtTurnOn_num_

private

Definition at line 123 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* LepHTMonitor::h_NPVTurnOn_den_

private

Definition at line 132 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* LepHTMonitor::h_NPVTurnOn_num_

private

Definition at line 131 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* LepHTMonitor::h_pfHTTurnOn_den_

private

Definition at line 122 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* LepHTMonitor::h_pfHTTurnOn_num_

private

Definition at line 121 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

std::vector<double> LepHTMonitor::htbins_

private

Definition at line 109 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

double LepHTMonitor::htCut_

private

Definition at line 96 of file LepHTMonitor.h.

Referenced by analyze().

double LepHTMonitor::jetEtaCut_

private

Definition at line 94 of file LepHTMonitor.h.

Referenced by analyze().

double LepHTMonitor::jetPtCut_

private

Definition at line 93 of file LepHTMonitor.h.

Referenced by analyze().

double LepHTMonitor::lep_counting_threshold_

private

Definition at line 100 of file LepHTMonitor.h.

Referenced by analyze().

double LepHTMonitor::lep_d0_cut_b_

private

Definition at line 103 of file LepHTMonitor.h.

Referenced by analyze().

double LepHTMonitor::lep_d0_cut_e_

private

Definition at line 105 of file LepHTMonitor.h.

Referenced by analyze().

double LepHTMonitor::lep_dz_cut_b_

private

Definition at line 104 of file LepHTMonitor.h.

Referenced by analyze().

double LepHTMonitor::lep_dz_cut_e_

private

Definition at line 106 of file LepHTMonitor.h.

Referenced by analyze().

double LepHTMonitor::lep_eta_cut_

private

Definition at line 102 of file LepHTMonitor.h.

Referenced by analyze().

double LepHTMonitor::lep_iso_cut_

private

Definition at line 101 of file LepHTMonitor.h.

Referenced by analyze().

double LepHTMonitor::lep_pt_plateau_

private

Definition at line 99 of file LepHTMonitor.h.

Referenced by analyze().

double LepHTMonitor::metCut_

private

Definition at line 95 of file LepHTMonitor.h.

Referenced by analyze().

int LepHTMonitor::muonIDlevel_

private

Definition at line 91 of file LepHTMonitor.h.

Referenced by analyze().

int LepHTMonitor::nbins_eta_

private

Definition at line 110 of file LepHTMonitor.h.

Referenced by bookHistograms().

int LepHTMonitor::nbins_npv_

private

Definition at line 112 of file LepHTMonitor.h.

Referenced by bookHistograms().

int LepHTMonitor::nbins_phi_

private

Definition at line 111 of file LepHTMonitor.h.

Referenced by bookHistograms().

double LepHTMonitor::nelsCut_

private

Definition at line 98 of file LepHTMonitor.h.

Referenced by analyze().

double LepHTMonitor::nmusCut_

private

Definition at line 97 of file LepHTMonitor.h.

Referenced by analyze().

float LepHTMonitor::npvbins_max_

private

Definition at line 118 of file LepHTMonitor.h.

Referenced by bookHistograms().

float LepHTMonitor::npvbins_min_

private

Definition at line 117 of file LepHTMonitor.h.

Referenced by bookHistograms().

std::unique_ptr<GenericTriggerEventFlag> LepHTMonitor::num_genTriggerEventFlag_

private

Definition at line 85 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

float LepHTMonitor::phibins_max_

private

Definition at line 116 of file LepHTMonitor.h.

Referenced by bookHistograms().

float LepHTMonitor::phibins_min_

private

Definition at line 115 of file LepHTMonitor.h.

Referenced by bookHistograms().

std::vector<double> LepHTMonitor::ptbins_

private

Definition at line 108 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 84 of file LepHTMonitor.h.

Referenced by analyze().

edm::InputTag LepHTMonitor::theBeamSpotTag_

private

Definition at line 83 of file LepHTMonitor.h.

Referenced by analyze().

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

private

Definition at line 82 of file LepHTMonitor.h.

Referenced by analyze().

edm::InputTag LepHTMonitor::theConversionCollectionTag_

private

Definition at line 81 of file LepHTMonitor.h.

Referenced by analyze().

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

private

Definition at line 67 of file LepHTMonitor.h.

Referenced by analyze().

edm::InputTag LepHTMonitor::theElectronTag_

private

Definition at line 66 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 69 of file LepHTMonitor.h.

Referenced by analyze().

edm::InputTag LepHTMonitor::theElectronVIDTag_

private

Definition at line 68 of file LepHTMonitor.h.

Referenced by analyze().

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

private

Definition at line 77 of file LepHTMonitor.h.

edm::InputTag LepHTMonitor::theJetTagTag_

private

Definition at line 76 of file LepHTMonitor.h.

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

private

Definition at line 71 of file LepHTMonitor.h.

Referenced by analyze().

edm::InputTag LepHTMonitor::theMuonTag_

private

Definition at line 70 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 75 of file LepHTMonitor.h.

Referenced by analyze().

edm::InputTag LepHTMonitor::thePfJetTag_

private

Definition at line 74 of file LepHTMonitor.h.

Referenced by analyze().

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

private

Definition at line 73 of file LepHTMonitor.h.

Referenced by analyze().

edm::InputTag LepHTMonitor::thePfMETTag_

private

Definition at line 72 of file LepHTMonitor.h.

Referenced by analyze().

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

private

Definition at line 80 of file LepHTMonitor.h.

Referenced by analyze().

edm::InputTag LepHTMonitor::theVertexCollectionTag_

private

Definition at line 79 of file LepHTMonitor.h.

Referenced by analyze().