#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 &ev, edm::EventSetup const &es) final

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual void	bookHistograms (DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &)=0

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) override

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

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

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)

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	endRun (const edm::Run &run, const edm::EventSetup &eSetup) override

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 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 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 54 of file LepHTMonitor.h.

Constructor & Destructor Documentation

LepHTMonitor::LepHTMonitor ( const edm::ParameterSet & ps )

Definition at line 112 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 183 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 249 of file LepHTMonitor.cc.

References funct::abs(), align::BeamSpot, den_HT_genTriggerEventFlag_, den_lep_genTriggerEventFlag_, metsig::electron, 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_, nelsCut_, nmusCut_, num_genTriggerEventFlag_, pfMETBenchmark_cfi::pfMET, 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 192 of file LepHTMonitor.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), 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_, 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 188 of file LepHTMonitor.cc.

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

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

overrideprotected

Definition at line 470 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 89 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

std::unique_ptr<GenericTriggerEventFlag> LepHTMonitor::den_lep_genTriggerEventFlag_

private

Definition at line 88 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

float LepHTMonitor::etabins_max_

private

Definition at line 117 of file LepHTMonitor.h.

Referenced by bookHistograms().

float LepHTMonitor::etabins_min_

private

Definition at line 116 of file LepHTMonitor.h.

Referenced by bookHistograms().

std::string LepHTMonitor::folderName_

private

Definition at line 91 of file LepHTMonitor.h.

Referenced by bookHistograms().

MonitorElement* LepHTMonitor::h_lepEtaPhiTurnOn_den_

private

Definition at line 134 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* LepHTMonitor::h_lepEtaPhiTurnOn_num_

private

Definition at line 133 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* LepHTMonitor::h_lepEtaTurnOn_den_

private

Definition at line 130 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* LepHTMonitor::h_lepEtaTurnOn_num_

private

Definition at line 129 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* LepHTMonitor::h_lepPhiTurnOn_den_

private

Definition at line 132 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* LepHTMonitor::h_lepPhiTurnOn_num_

private

Definition at line 131 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* LepHTMonitor::h_lepPtTurnOn_den_

private

Definition at line 128 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* LepHTMonitor::h_lepPtTurnOn_num_

private

Definition at line 127 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* LepHTMonitor::h_NPVTurnOn_den_

private

Definition at line 136 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* LepHTMonitor::h_NPVTurnOn_num_

private

Definition at line 135 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* LepHTMonitor::h_pfHTTurnOn_den_

private

Definition at line 126 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* LepHTMonitor::h_pfHTTurnOn_num_

private

Definition at line 125 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

std::vector<double> LepHTMonitor::htbins_

private

Definition at line 112 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

double LepHTMonitor::htCut_

private

Definition at line 98 of file LepHTMonitor.h.

Referenced by analyze().

double LepHTMonitor::jetEtaCut_

private

Definition at line 96 of file LepHTMonitor.h.

Referenced by analyze().

double LepHTMonitor::jetPtCut_

private

Definition at line 95 of file LepHTMonitor.h.

Referenced by analyze().

double LepHTMonitor::lep_counting_threshold_

private

Definition at line 102 of file LepHTMonitor.h.

Referenced by analyze().

double LepHTMonitor::lep_d0_cut_b_

private

Definition at line 105 of file LepHTMonitor.h.

Referenced by analyze().

double LepHTMonitor::lep_d0_cut_e_

private

Definition at line 107 of file LepHTMonitor.h.

Referenced by analyze().

double LepHTMonitor::lep_dz_cut_b_

private

Definition at line 106 of file LepHTMonitor.h.

Referenced by analyze().

double LepHTMonitor::lep_dz_cut_e_

private

Definition at line 108 of file LepHTMonitor.h.

Referenced by analyze().

double LepHTMonitor::lep_eta_cut_

private

Definition at line 104 of file LepHTMonitor.h.

Referenced by analyze().

double LepHTMonitor::lep_iso_cut_

private

Definition at line 103 of file LepHTMonitor.h.

Referenced by analyze().

double LepHTMonitor::lep_pt_plateau_

private

Definition at line 101 of file LepHTMonitor.h.

Referenced by analyze().

double LepHTMonitor::metCut_

private

Definition at line 97 of file LepHTMonitor.h.

Referenced by analyze().

int LepHTMonitor::muonIDlevel_

private

Definition at line 93 of file LepHTMonitor.h.

Referenced by analyze().

int LepHTMonitor::nbins_eta_

private

Definition at line 113 of file LepHTMonitor.h.

Referenced by bookHistograms().

int LepHTMonitor::nbins_npv_

private

Definition at line 115 of file LepHTMonitor.h.

Referenced by bookHistograms().

int LepHTMonitor::nbins_phi_

private

Definition at line 114 of file LepHTMonitor.h.

Referenced by bookHistograms().

double LepHTMonitor::nelsCut_

private

Definition at line 100 of file LepHTMonitor.h.

Referenced by analyze().

double LepHTMonitor::nmusCut_

private

Definition at line 99 of file LepHTMonitor.h.

Referenced by analyze().

float LepHTMonitor::npvbins_max_

private

Definition at line 121 of file LepHTMonitor.h.

Referenced by bookHistograms().

float LepHTMonitor::npvbins_min_

private

Definition at line 120 of file LepHTMonitor.h.

Referenced by bookHistograms().

std::unique_ptr<GenericTriggerEventFlag> LepHTMonitor::num_genTriggerEventFlag_

private

Definition at line 87 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

float LepHTMonitor::phibins_max_

private

Definition at line 119 of file LepHTMonitor.h.

Referenced by bookHistograms().

float LepHTMonitor::phibins_min_

private

Definition at line 118 of file LepHTMonitor.h.

Referenced by bookHistograms().

std::vector<double> LepHTMonitor::ptbins_

private

Definition at line 111 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 86 of file LepHTMonitor.h.

Referenced by analyze().

edm::InputTag LepHTMonitor::theBeamSpotTag_

private

Definition at line 85 of file LepHTMonitor.h.

Referenced by analyze().

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

private

Definition at line 84 of file LepHTMonitor.h.

Referenced by analyze().

edm::InputTag LepHTMonitor::theConversionCollectionTag_

private

Definition at line 83 of file LepHTMonitor.h.

Referenced by analyze().

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

private

Definition at line 69 of file LepHTMonitor.h.

Referenced by analyze().

edm::InputTag LepHTMonitor::theElectronTag_

private

Definition at line 68 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 71 of file LepHTMonitor.h.

Referenced by analyze().

edm::InputTag LepHTMonitor::theElectronVIDTag_

private

Definition at line 70 of file LepHTMonitor.h.

Referenced by analyze().

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

private

Definition at line 79 of file LepHTMonitor.h.

edm::InputTag LepHTMonitor::theJetTagTag_

private

Definition at line 78 of file LepHTMonitor.h.

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

private

Definition at line 73 of file LepHTMonitor.h.

Referenced by analyze().

edm::InputTag LepHTMonitor::theMuonTag_

private

Definition at line 72 of file LepHTMonitor.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 77 of file LepHTMonitor.h.

Referenced by analyze().

edm::InputTag LepHTMonitor::thePfJetTag_

private

Definition at line 76 of file LepHTMonitor.h.

Referenced by analyze().

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

private

Definition at line 75 of file LepHTMonitor.h.

Referenced by analyze().

edm::InputTag LepHTMonitor::thePfMETTag_

private

Definition at line 74 of file LepHTMonitor.h.

Referenced by analyze().

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

private

Definition at line 82 of file LepHTMonitor.h.

Referenced by analyze().

edm::InputTag LepHTMonitor::theVertexCollectionTag_

private

Definition at line 81 of file LepHTMonitor.h.

Referenced by analyze().

Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation