#include <HLTriggerOffline/Top/src/TopDiLeptonHLTValidation.cc>

Inheritance diagram for TopDiLeptonHLTValidation:

Public Member Functions
	TopDiLeptonHLTValidation (const edm::ParameterSet &)

	~TopDiLeptonHLTValidation () override

Public Member Functions inherited from DQMEDAnalyzer
void	beginRun (edm::Run const &, edm::EventSetup const &) final

void	beginStream (edm::StreamID id) final

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

	DQMEDAnalyzer (void)

void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

uint32_t	streamId () const

Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default

Public Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

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

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr< dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)

static std::shared_ptr< dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)

static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const , dqmDetails::NoCache )

static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const , dqmDetails::NoCache )

Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

std::string	monitorPath (const std::string &label) const

void	triggerBinLabels (const std::vector< std::string > &labels)
	set configurable labels for trigger monitoring histograms More...

Private Attributes
const reco::GsfElectron *	elec1_

const reco::GsfElectron *	elec2_

double	etaElectrons_

double	etaJets_

double	etaMuons_

MonitorElement *	hDenJetEta

MonitorElement *	hDenJetPt

MonitorElement *	hDenLeptonEta

MonitorElement *	hDenLeptonPt

MonitorElement *	hDenTriggerMon

MonitorElement *	hNumJetEta

MonitorElement *	hNumJetPt

MonitorElement *	hNumLeptonEta

MonitorElement *	hNumLeptonPt

MonitorElement *	hNumTriggerMon

bool	isAll_ = false

double	isoElectrons_

double	isoMuons_

bool	isSel_ = false

edm::InputTag	iTrigger_

const reco::Jet *	jet_

unsigned int	minElectrons_

unsigned int	minJets_

unsigned int	minMuons_

const reco::Muon *	mu1_

const reco::Muon *	mu2_

double	ptElectrons_

double	ptJets_

double	ptMuons_

std::string	sDir_

std::string	sElectrons_

std::string	sJets_

std::string	sMuons_

edm::EDGetTokenT< edm::View< reco::GsfElectron > >	tokElectrons_

edm::EDGetTokenT< edm::View< reco::Jet > >	tokJets_

edm::EDGetTokenT< edm::View< reco::Muon > >	tokMuons_

edm::EDGetTokenT< edm::TriggerResults >	tokTrigger_

std::vector< std::string >	vsPaths_

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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)

Detailed Description

Description:

Description: compute efficiencies of trigger paths on offline reco selection with respect to pt and eta

Implementation: harvesting

Definition at line 51 of file TopDiLeptonHLTValidation.h.

Constructor & Destructor Documentation

TopDiLeptonHLTValidation::TopDiLeptonHLTValidation ( const edm::ParameterSet & iConfig )

explicit

Definition at line 134 of file TopDiLeptonHLTValidation.h.

References iTrigger_, sElectrons_, sJets_, sMuons_, tokElectrons_, tokJets_, tokMuons_, and tokTrigger_.

                                                                                  :
   sDir_(iConfig.getUntrackedParameter<std::string>("sDir","Validation/Top/Efficiencies/")),
   sElectrons_(iConfig.getUntrackedParameter<std::string>("sElectrons","gsfElectrons")),
   ptElectrons_(iConfig.getUntrackedParameter<double>("ptElectrons",0.)),
   etaElectrons_(iConfig.getUntrackedParameter<double>("etaElectrons",0.)),
   isoElectrons_(iConfig.getUntrackedParameter<double>("isoElectrons",0.)),
   minElectrons_(iConfig.getUntrackedParameter<unsigned int>("minElectrons",0)),
   sMuons_(iConfig.getUntrackedParameter<std::string>("sMuons","muons")),
   ptMuons_(iConfig.getUntrackedParameter<double>("ptMuons",0.)),
   etaMuons_(iConfig.getUntrackedParameter<double>("etaMuons",0.)),
   isoMuons_(iConfig.getUntrackedParameter<double>("isoMuons",0.)),
   minMuons_(iConfig.getUntrackedParameter<unsigned int>("minMuons",0)),
   sJets_(iConfig.getUntrackedParameter<std::string>("sJets","ak5PFJets")),
   ptJets_(iConfig.getUntrackedParameter<double>("ptJets",0.)),
   etaJets_(iConfig.getUntrackedParameter<double>("etaJets",0.)),
   minJets_(iConfig.getUntrackedParameter<unsigned int>("minJets",0)),
   iTrigger_(iConfig.getUntrackedParameter<edm::InputTag>("iTrigger")),
   vsPaths_(iConfig.getUntrackedParameter< std::vector<std::string> >("vsPaths"))
 
 {
   // Electrons
   tokElectrons_ = consumes< edm::View<reco::GsfElectron> >(edm::InputTag(sElectrons_));
   // Muons
   tokMuons_ = consumes< edm::View<reco::Muon> >(edm::InputTag(sMuons_));
   // Jets
   tokJets_ = consumes< edm::View<reco::Jet> >(edm::InputTag(sJets_));
   // Trigger
   tokTrigger_ = consumes<edm::TriggerResults>(iTrigger_); 
 }

TopDiLeptonHLTValidation::~TopDiLeptonHLTValidation ( )

override

Definition at line 165 of file TopDiLeptonHLTValidation.h.

References DEFINE_FWK_MODULE.

 {
  
    // do anything here that needs to be done at desctruction time
    // (e.g. close files, deallocate resources etc.)
 
 }

Member Function Documentation

void TopDiLeptonHLTValidation::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 46 of file TopDiLeptonHLTValidation.cc.

References edm::HLTGlobalStatus::accept(), MillePedeFileConverter_cfg::e, elec1_, elec2_, electrons_cff::electrons, reco::LeafCandidate::eta(), etaElectrons_, etaJets_, etaMuons_, MonitorElement::Fill(), edm::Event::getByToken(), hDenJetEta, hDenJetPt, hDenLeptonEta, hDenLeptonPt, hDenTriggerMon, hNumJetEta, hNumJetPt, hNumLeptonEta, hNumLeptonPt, hNumTriggerMon, mps_fire::i, training_settings::idx, isAll_, isoElectrons_, isoMuons_, isSel_, jet_, fwrapper::jets, funct::m, minElectrons_, minJets_, minMuons_, mu1_, mu2_, extraflags_cff::muons, dataset::name, reco::LeafCandidate::pt(), ptElectrons_, ptJets_, ptMuons_, tokElectrons_, tokJets_, tokMuons_, tokTrigger_, edm::TriggerNames::triggerNames(), edm::Event::triggerNames(), and vsPaths_.

 {
   using namespace edm;
 
   isAll_ = false; isSel_ = false;
 
   // Electrons 
   Handle< edm::View<reco::GsfElectron> > electrons;
   if (!iEvent.getByToken(tokElectrons_,electrons))
     edm::LogWarning("TopDiLeptonHLTValidation") << "Electrons collection not found \n";
   unsigned int nGoodE = 0;
   for(edm::View<reco::GsfElectron>::const_iterator e = electrons->begin(); e != electrons->end(); ++e){
     if (e->pt() < ptElectrons_) continue;
     if (fabs(e->eta()) > etaElectrons_) continue;
     if ((e->dr03TkSumPt()+e->dr03EcalRecHitSumEt()+e->dr03HcalTowerSumEt())/e->pt() > isoElectrons_ ) continue;
     nGoodE++;
     if (nGoodE == 1) elec1_ = &(*e);
     if (nGoodE == 2) elec2_ = &(*e);
   }
   // Muons 
   Handle< edm::View<reco::Muon> > muons;
   if (!iEvent.getByToken(tokMuons_,muons))
     edm::LogWarning("TopDiLeptonHLTValidation") << "Muons collection not found \n";
   unsigned int nGoodM = 0;
   for(edm::View<reco::Muon>::const_iterator m = muons->begin(); m != muons->end(); ++m){
     if (!m->isPFMuon() || (!m->isGlobalMuon() && !m->isTrackerMuon())) continue;
     if (m->pt() < ptMuons_) continue;
     if (fabs(m->eta()) > etaMuons_) continue;
     if (((m->pfIsolationR04()).sumChargedHadronPt+(m->pfIsolationR04()).sumPhotonEt+(m->pfIsolationR04()).sumNeutralHadronEt)/m->pt() > isoMuons_ ) continue;
     nGoodM++;
     if (nGoodM == 1) mu1_ = &(*m);
     if (nGoodM == 2) mu2_ = &(*m);
   }
   // Jets 
   Handle< edm::View<reco::Jet> > jets;
   if (!iEvent.getByToken(tokJets_,jets)) 
     edm::LogWarning("TopDiLeptonHLTValidation") << "Jets collection not found \n";
   unsigned int nGoodJ = 0;
   for(edm::View<reco::Jet>::const_iterator j = jets->begin(); j != jets->end(); ++j){
     if (j->pt() < ptJets_) continue;
     if (fabs(j->eta()) > etaJets_) continue;
     nGoodJ++;
     if (nGoodJ == minJets_) jet_ = &(*j);
   }
 
   if (nGoodE >= minElectrons_ && nGoodM >= minMuons_ && nGoodJ >= minJets_) isAll_ = true;
 
   //Trigger
   Handle<edm::TriggerResults> triggerTable;
   if (!iEvent.getByToken(tokTrigger_,triggerTable))
       edm::LogWarning("TopDiLeptonHLTValidation") << "Trigger collection not found \n";
   const edm::TriggerNames& triggerNames = iEvent.triggerNames(*triggerTable);
   unsigned int isInteresting = 0;
   for (unsigned int i=0; i<triggerNames.triggerNames().size(); ++i) {
 
       TString name = triggerNames.triggerNames()[i].c_str();
       for (unsigned int j=0; j<vsPaths_.size(); j++) {
           if (name.Contains(TString(vsPaths_[j]), TString::kIgnoreCase)) {
             if (triggerTable->accept(i)){
               isInteresting++;  
               if (isAll_) hNumTriggerMon->Fill(j+0.5 ); 
             }  
           }
       }
   }
 
   if (isAll_ && isInteresting > 0) isSel_ = true;
   else isSel_ = false;
 
   //Histos
   if (isAll_) {
       if (minElectrons_ > 0) {
           hDenLeptonPt->Fill(elec1_->pt());
           hDenLeptonEta->Fill(elec1_->eta());
       }
       if (minElectrons_ > 1) {
           hDenLeptonPt->Fill(elec2_->pt());
           hDenLeptonEta->Fill(elec2_->eta());
       }
       if (minMuons_ > 0) {
           hDenLeptonPt->Fill(mu1_->pt());
           hDenLeptonEta->Fill(mu1_->eta());
       }
       if (minMuons_ > 1) {
           hDenLeptonPt->Fill(mu2_->pt());
           hDenLeptonEta->Fill(mu2_->eta());
       }
       hDenJetPt->Fill(jet_->pt());
       hDenJetEta->Fill(jet_->eta());
       for(unsigned int idx=0; idx<vsPaths_.size(); ++idx){
           hDenTriggerMon->Fill(idx+0.5);
       }
   }
   if (isSel_) {
       if (minElectrons_ > 0) {
           hNumLeptonPt->Fill(elec1_->pt());
           hNumLeptonEta->Fill(elec1_->eta());
       }
       if (minElectrons_ > 1) {
           hNumLeptonPt->Fill(elec2_->pt());
           hNumLeptonEta->Fill(elec2_->eta());
       }
       if (minMuons_ > 0) {
           hNumLeptonPt->Fill(mu1_->pt());
           hNumLeptonEta->Fill(mu1_->eta());
       }
       if (minMuons_ > 1) {
           hNumLeptonPt->Fill(mu2_->pt());
           hNumLeptonEta->Fill(mu2_->eta());
       }
       hNumJetPt->Fill(jet_->pt());
       hNumJetEta->Fill(jet_->eta());
   }
 }

void TopDiLeptonHLTValidation::bookHistograms	(	DQMStore::IBooker &	dbe,
		edm::Run const &	,
		edm::EventSetup const &
	)

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 163 of file TopDiLeptonHLTValidation.cc.

References DQMStore::IBooker::book1D(), hDenJetEta, hDenJetPt, hDenLeptonEta, hDenLeptonPt, hDenTriggerMon, hNumJetEta, hNumJetPt, hNumLeptonEta, hNumLeptonPt, hNumTriggerMon, sDir_, DQMStore::IBooker::setCurrentFolder(), triggerBinLabels(), and vsPaths_.

 {
   dbe.setCurrentFolder(sDir_);
   hDenLeptonPt  = dbe.book1D("PtLeptonAll", "PtLeptonAll", 50, 0., 250.);
   hDenLeptonEta = dbe.book1D("EtaLeptonAll", "EtaLeptonAll", 30, -3. , 3.);
   hDenJetPt     = dbe.book1D("PtLastJetAll", "PtLastJetAll", 60, 0., 300.);
   hDenJetEta    = dbe.book1D("EtaLastJetAll", "EtaLastJetAll", 30, -3., 3.);
   hNumLeptonPt  = dbe.book1D("PtLeptonSel", "PtLeptonSel", 50, 0., 250.);
   hNumLeptonEta = dbe.book1D("EtaLeptonSel", "EtaLeptonSel", 30, -3. , 3.);
   hNumJetPt     = dbe.book1D("PtLastJetSel", "PtLastJetSel", 60, 0., 300.);
   hNumJetEta    = dbe.book1D("EtaLastJetSel", "EtaLastJetSel", 30, -3., 3.);
   // determine number of bins for trigger monitoring
   unsigned int nPaths = vsPaths_.size();
   // monitored trigger occupancy for single lepton triggers
   hNumTriggerMon    = dbe.book1D("TriggerMonSel" , "TriggerMonSel", nPaths, 0.,  nPaths);
   hDenTriggerMon    = dbe.book1D("TriggerMonAll" , "TriggerMonAll", nPaths, 0.,  nPaths);
   // set bin labels for trigger monitoring
   triggerBinLabels(vsPaths_);
 }

void TopDiLeptonHLTValidation::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 185 of file TopDiLeptonHLTValidation.cc.

References edm::ConfigurationDescriptions::addDefault(), and edm::ParameterSetDescription::setUnknown().

                                                                                      {
   //The following says we do not know what parameters are allowed so do no validation
   // Please change this to state exactly what you do use, even if it is no parameters
   edm::ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }

std::string TopDiLeptonHLTValidation::monitorPath ( const std::string & label ) const

inlineprivate

deduce monitorPath from label, the label is expected to be of type 'selectionPath:monitorPath'

Definition at line 64 of file TopDiLeptonHLTValidation.h.

References tablePrinter::labels, and triggerBinLabels().

Referenced by triggerBinLabels().

64 { return label.substr(label.find(':')+1); };

diffTwoXMLs.label

label

Definition: diffTwoXMLs.py:42

void TopDiLeptonHLTValidation::triggerBinLabels ( const std::vector< std::string > & labels )

inlineprivate

set configurable labels for trigger monitoring histograms

Definition at line 115 of file TopDiLeptonHLTValidation.h.

References hDenTriggerMon, hNumTriggerMon, training_settings::idx, monitorPath(), and MonitorElement::setBinLabel().

Referenced by bookHistograms(), and monitorPath().

 {
   for(unsigned int idx=0; idx<labels.size(); ++idx){
     hNumTriggerMon->setBinLabel( idx+1, "["+monitorPath(labels[idx])+"]", 1);
     hDenTriggerMon->setBinLabel( idx+1, "["+monitorPath(labels[idx])+"]", 1);
   }
 }