#include <TopSingleLeptonHLTValidation.h>

Inheritance diagram for TopSingleLeptonHLTValidation:

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

	~TopSingleLeptonHLTValidation () override

Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

Static Public Member Functions
static void	fillDescriptions (edm::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 *	elec_

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 *	mu_

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_

Detailed Description

HLTriggerOffline/Top/interface/TopSingleLeptonHLTValidation.h

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

Implementation: harvesting

HLTriggerOffline/Top/src/TopSingleLeptonHLTValidation.cc

Description:

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

Implementation: harvesting

Definition at line 50 of file TopSingleLeptonHLTValidation.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 129 of file TopSingleLeptonHLTValidation.h.

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

     : sDir_(iConfig.getUntrackedParameter<std::string>("sDir", "HLTValidation/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_);
 }

TopSingleLeptonHLTValidation::~TopSingleLeptonHLTValidation ( )

override

Definition at line 159 of file TopSingleLeptonHLTValidation.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 TopSingleLeptonHLTValidation::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 46 of file TopSingleLeptonHLTValidation.cc.

References edm::HLTGlobalStatus::accept(), MillePedeFileConverter_cfg::e, elec_, nano_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_, mu_, extraflags_cff::muons, dataset::name, EnergyCorrector::pt, 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("TopSingleLeptonHLTValidation") << "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)
       elec_ = &(*e);
   }
   // Muons
   Handle<edm::View<reco::Muon>> muons;
   if (!iEvent.getByToken(tokMuons_, muons))
     edm::LogWarning("TopSingleLeptonHLTValidation") << "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())
       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)
       mu_ = &(*m);
   }
   // Jets
   Handle<edm::View<reco::Jet>> jets;
   if (!iEvent.getByToken(tokJets_, jets))
     edm::LogWarning("TopSingleLeptonHLTValidation") << "Jets collection not found \n";
   unsigned int nGoodJ = 0;
   if (minJets_ == 4) {
     for (edm::View<reco::Jet>::const_iterator j = jets->begin(); j != jets->end(); ++j) {
       if (nGoodJ < 1 && j->pt() < 55)
         continue;
       if (nGoodJ < 2 && j->pt() < 45)
         continue;
       if (nGoodJ < 3 && j->pt() < 35)
         continue;
       if (nGoodJ >= 3 && j->pt() < 20)
         continue;
       if (fabs(j->eta()) > etaJets_)
         continue;
       nGoodJ++;
       if (nGoodJ == minJets_)
         jet_ = &(*j);
     }
   } else {
     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("TopSingleLeptonHLTValidation") << "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(elec_->pt());
       hDenLeptonEta->Fill(elec_->eta());
     }
     if (minMuons_ > 0) {
       hDenLeptonPt->Fill(mu_->pt());
       hDenLeptonEta->Fill(mu_->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(elec_->pt());
       hNumLeptonEta->Fill(elec_->eta());
     }
     if (minMuons_ > 0) {
       hNumLeptonPt->Fill(mu_->pt());
       hNumLeptonEta->Fill(mu_->eta());
     }
     hNumJetPt->Fill(jet_->pt());
     hNumJetEta->Fill(jet_->eta());
   }
 }

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

overrideprivate

Definition at line 180 of file TopSingleLeptonHLTValidation.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 TopSingleLeptonHLTValidation::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 201 of file TopSingleLeptonHLTValidation.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 TopSingleLeptonHLTValidation::monitorPath ( const std::string & label ) const

inlineprivate

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

Definition at line 62 of file TopSingleLeptonHLTValidation.h.

References tablePrinter::labels, and triggerBinLabels().

Referenced by triggerBinLabels().

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

label

char const * label

Definition: PFTauDecayModeTools.cc:11

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

inlineprivate

set configurable labels for trigger monitoring histograms

Definition at line 111 of file TopSingleLeptonHLTValidation.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);
   }
 }