#include <L1TTauOffline.h>

Inheritance diagram for L1TTauOffline:

Public Types
enum	PlotConfig { nVertex, ETvsET, PHIvsPHI }

Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Public Member Functions
	L1TTauOffline (const edm::ParameterSet &ps)

	~L1TTauOffline () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final

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

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

void	beginStream (edm::StreamID id) final

	DQMEDAnalyzer ()

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

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

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

Static Public Attributes
static const std::map< std::string, unsigned int >	PlotConfigNames

Protected Member Functions
void	analyze (edm::Event const &e, edm::EventSetup const &eSetup) override

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

void	dqmBeginRun (const edm::Run &run, const edm::EventSetup &iSetup) override

const reco::Vertex	getPrimaryVertex (edm::Handle< reco::VertexCollection > const &vertex, edm::Handle< reco::BeamSpot > const &beamSpot)

void	getProbeTaus (const edm::Event &e, edm::Handle< reco::PFTauCollection > const &taus, edm::Handle< reco::MuonCollection > const &muons, const reco::Vertex &vertex)

void	getTauL1tPairs (edm::Handle< l1t::TauBxCollection > const &l1tCands)

void	getTightMuons (edm::Handle< reco::MuonCollection > const &muons, edm::Handle< reco::PFMETCollection > const &mets, const reco::Vertex &vertex, edm::Handle< trigger::TriggerEvent > const &trigEvent)

bool	matchHlt (edm::Handle< trigger::TriggerEvent > const &triggerEvent, const reco::Muon *muon)

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Member Functions
void	bookTauHistos (DQMStore::IBooker &)

double	calcDeltaPhi (double phi1, double phi2)

double	Distance (const reco::Candidate &c1, const reco::Candidate &c2)

double	DistancePhi (const reco::Candidate &c1, const reco::Candidate &c2)

void	normalise2DHistogramsToBinArea ()

Private Attributes
edm::EDGetTokenT< reco::TauDiscriminatorContainer >	AntiEleInputTag_

std::string	AntiEleWP_

int	AntiEleWPIndex_

edm::EDGetTokenT< reco::TauDiscriminatorContainer >	AntiMuInputTag_

std::string	AntiMuWP_

int	AntiMuWPIndex_

edm::EDGetTokenT< reco::BeamSpot >	BsInputTag_

edm::EDGetTokenT< reco::TauDiscriminatorContainer >	comb3TInputTag_

std::string	comb3TWP_

int	comb3TWPIndex_

edm::EDGetTokenT< reco::PFTauDiscriminator >	DecayModeFindingInputTag_

std::string	efficiencyFolder_

std::map< double, MonitorElement * >	h_efficiencyIsoTauET_EB_EE_pass_

std::map< double, MonitorElement * >	h_efficiencyIsoTauET_EB_EE_total_

std::map< double, MonitorElement * >	h_efficiencyIsoTauET_EB_pass_

std::map< double, MonitorElement * >	h_efficiencyIsoTauET_EB_total_

std::map< double, MonitorElement * >	h_efficiencyIsoTauET_EE_pass_

std::map< double, MonitorElement * >	h_efficiencyIsoTauET_EE_total_

std::map< double, MonitorElement * >	h_efficiencyNonIsoTauET_EB_EE_pass_

std::map< double, MonitorElement * >	h_efficiencyNonIsoTauET_EB_EE_total_

std::map< double, MonitorElement * >	h_efficiencyNonIsoTauET_EB_pass_

std::map< double, MonitorElement * >	h_efficiencyNonIsoTauET_EB_total_

std::map< double, MonitorElement * >	h_efficiencyNonIsoTauET_EE_pass_

std::map< double, MonitorElement * >	h_efficiencyNonIsoTauET_EE_total_

MonitorElement *	h_L1TauEtavsTauEta_

MonitorElement *	h_L1TauETvsTauET_EB_

MonitorElement *	h_L1TauETvsTauET_EB_EE_

MonitorElement *	h_L1TauETvsTauET_EE_

MonitorElement *	h_L1TauPhivsTauPhi_EB_

MonitorElement *	h_L1TauPhivsTauPhi_EB_EE_

MonitorElement *	h_L1TauPhivsTauPhi_EE_

MonitorElement *	h_nVertex_

MonitorElement *	h_resolutionTauET_EB_

MonitorElement *	h_resolutionTauET_EB_EE_

MonitorElement *	h_resolutionTauET_EE_

MonitorElement *	h_resolutionTauEta_

MonitorElement *	h_resolutionTauPhi_EB_

MonitorElement *	h_resolutionTauPhi_EB_EE_

MonitorElement *	h_resolutionTauPhi_EE_

MonitorElement *	h_tagAndProbeMass_

dqmoffline::l1t::HistDefinitions	histDefinitions_

std::string	histFolder_

bool	m_AntiEleExists

edm::ESHandle< MagneticField >	m_BField

HLTConfigProvider	m_hltConfig

BXVector< l1t::Tau >	m_L1tL1tTaus

std::vector< int >	m_L1tPtCuts

std::vector< l1t::TauBxCollection >	m_L1tTaus

float	m_MaxHltTauDR

float	m_MaxL1tTauDR

float	m_MaxTauEta

std::vector< const reco::PFTau * >	m_ProbeTaus

edm::ESHandle< Propagator >	m_propagatorAlong

edm::ESHandle< Propagator >	m_propagatorOpposite

std::vector< reco::PFTau >	m_RecoRecoTaus

std::vector< reco::PFTauCollection >	m_RecoTaus

std::vector< TauL1TPair >	m_TauL1tPairs

std::vector< const reco::Muon * >	m_TightMuons

std::vector< int >	m_trigIndices

edm::EDGetTokenT< reco::PFMETCollection >	MetInputTag_

edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	mFieldToken_

edm::EDGetTokenT< reco::MuonCollection >	MuonInputTag_

edm::ProcessHistoryID	phID_

math::XYZPoint	PVPoint_

edm::EDGetTokenT< l1t::TauBxCollection >	stage2CaloLayer2TauToken_

std::vector< double >	tauEfficiencyBins_

std::vector< int >	tauEfficiencyThresholds_

edm::EDGetTokenT< reco::PFTauCollection >	theTauCollection_

edm::EDGetTokenT< trigger::TriggerEvent >	triggerEvent_

std::vector< std::string >	triggerPath_

edm::EDGetTokenT< edm::TriggerResults >	triggerResults_

std::string	trigProcess_

edm::EDGetTokenT< reco::VertexCollection >	VtxInputTag_

Additional Inherited Members
Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Definition at line 80 of file L1TTauOffline.h.

Member Enumeration Documentation

◆ PlotConfig

enum L1TTauOffline::PlotConfig

Enumerator
nVertex
ETvsET
PHIvsPHI

Definition at line 85 of file L1TTauOffline.h.

85 { nVertex, ETvsET, PHIvsPHI };

L1TTauOffline::ETvsET

Definition: L1TTauOffline.h:85

L1TTauOffline::nVertex

Definition: L1TTauOffline.h:85

L1TTauOffline::PHIvsPHI

Definition: L1TTauOffline.h:85

Constructor & Destructor Documentation

◆ L1TTauOffline()

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

Definition at line 45 of file L1TTauOffline.cc.

References deDxTools::esConsumes(), and mFieldToken_.

     : theTauCollection_(consumes<reco::PFTauCollection>(ps.getUntrackedParameter<edm::InputTag>("tauInputTag"))),
       AntiMuInputTag_(
           consumes<reco::TauDiscriminatorContainer>(ps.getUntrackedParameter<edm::InputTag>("antiMuInputTag"))),
       AntiMuWP_(ps.getUntrackedParameter<std::string>("antiMuWP")),
       AntiEleInputTag_(
           consumes<reco::TauDiscriminatorContainer>(ps.getUntrackedParameter<edm::InputTag>("antiEleInputTag"))),
       AntiEleWP_(ps.getUntrackedParameter<std::string>("antiEleWP")),
       DecayModeFindingInputTag_(
           consumes<reco::PFTauDiscriminator>(ps.getUntrackedParameter<edm::InputTag>("decayModeFindingInputTag"))),
       comb3TInputTag_(
           consumes<reco::TauDiscriminatorContainer>(ps.getUntrackedParameter<edm::InputTag>("comb3TInputTag"))),
       comb3TWP_(ps.getUntrackedParameter<std::string>("comb3TWP")),
       MuonInputTag_(consumes<reco::MuonCollection>(ps.getUntrackedParameter<edm::InputTag>("muonInputTag"))),
       MetInputTag_(consumes<reco::PFMETCollection>(ps.getUntrackedParameter<edm::InputTag>("metInputTag"))),
       VtxInputTag_(consumes<reco::VertexCollection>(ps.getUntrackedParameter<edm::InputTag>("vtxInputTag"))),
       BsInputTag_(consumes<reco::BeamSpot>(ps.getUntrackedParameter<edm::InputTag>("bsInputTag"))),
       triggerEvent_(consumes<trigger::TriggerEvent>(ps.getUntrackedParameter<edm::InputTag>("trigInputTag"))),
       trigProcess_(ps.getUntrackedParameter<string>("trigProcess")),
       triggerResults_(consumes<edm::TriggerResults>(ps.getUntrackedParameter<edm::InputTag>("trigProcess_token"))),
       triggerPath_(ps.getUntrackedParameter<vector<std::string>>("triggerNames")),
       histFolder_(ps.getParameter<std::string>("histFolder")),
       efficiencyFolder_(histFolder_ + "/efficiency_raw"),
       stage2CaloLayer2TauToken_(consumes<l1t::TauBxCollection>(ps.getUntrackedParameter<edm::InputTag>("l1tInputTag"))),
       tauEfficiencyThresholds_(ps.getParameter<std::vector<int>>("tauEfficiencyThresholds")),
       tauEfficiencyBins_(ps.getParameter<std::vector<double>>("tauEfficiencyBins")),
       histDefinitions_(dqmoffline::l1t::readHistDefinitions(ps.getParameterSet("histDefinitions"), PlotConfigNames)),
       m_TightMuons(),
       m_ProbeTaus(),
       m_TauL1tPairs(),
       m_RecoTaus(),
       m_L1tTaus(),
       m_RecoRecoTaus(),
       m_L1tL1tTaus(),
       m_L1tPtCuts(),
       m_MaxTauEta(99999),
       m_MaxL1tTauDR(99999),
       m_MaxHltTauDR(99999),
       m_trigIndices(),
       h_nVertex_(),
       h_tagAndProbeMass_(),
       h_L1TauETvsTauET_EB_(),
       h_L1TauETvsTauET_EE_(),
       h_L1TauETvsTauET_EB_EE_(),
       h_L1TauPhivsTauPhi_EB_(),
       h_L1TauPhivsTauPhi_EE_(),
       h_L1TauPhivsTauPhi_EB_EE_(),
       h_L1TauEtavsTauEta_(),
       h_resolutionTauET_EB_(),
       h_resolutionTauET_EE_(),
       h_resolutionTauET_EB_EE_(),
       h_resolutionTauPhi_EB_(),
       h_resolutionTauPhi_EE_(),
       h_resolutionTauPhi_EB_EE_(),
       h_resolutionTauEta_(),
       h_efficiencyIsoTauET_EB_pass_(),
       h_efficiencyIsoTauET_EE_pass_(),
       h_efficiencyIsoTauET_EB_EE_pass_(),
       h_efficiencyNonIsoTauET_EB_pass_(),
       h_efficiencyNonIsoTauET_EE_pass_(),
       h_efficiencyNonIsoTauET_EB_EE_pass_(),
       h_efficiencyIsoTauET_EB_total_(),
       h_efficiencyIsoTauET_EE_total_(),
       h_efficiencyIsoTauET_EB_EE_total_(),
       h_efficiencyNonIsoTauET_EB_total_(),
       h_efficiencyNonIsoTauET_EE_total_(),
       h_efficiencyNonIsoTauET_EB_EE_total_() {
   edm::LogInfo("L1TTauOffline") << "Constructor "
                                 << "L1TTauOffline::L1TTauOffline " << std::endl;
   mFieldToken_ = esConsumes();
 }

◆ ~L1TTauOffline()

L1TTauOffline::~L1TTauOffline ( )

override

Definition at line 120 of file L1TTauOffline.cc.

                               {
   edm::LogInfo("L1TTauOffline") << "Destructor L1TTauOffline::~L1TTauOffline " << std::endl;
 }

Member Function Documentation

◆ analyze()

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 171 of file L1TTauOffline.cc.

References pwdgSkimBPark_cfi::beamSpot, BXVector< T >::begin(), BsInputTag_, counter, MillePedeFileConverter_cfg::e, BXVector< T >::end(), PVValHelper::eta, dqm::impl::MonitorElement::Fill(), edm::EventSetup::getHandle(), getPrimaryVertex(), getProbeTaus(), getTauL1tPairs(), getTightMuons(), h_efficiencyIsoTauET_EB_EE_pass_, h_efficiencyIsoTauET_EB_EE_total_, h_efficiencyIsoTauET_EB_pass_, h_efficiencyIsoTauET_EB_total_, h_efficiencyIsoTauET_EE_pass_, h_efficiencyIsoTauET_EE_total_, h_efficiencyNonIsoTauET_EB_EE_pass_, h_efficiencyNonIsoTauET_EB_EE_total_, h_efficiencyNonIsoTauET_EB_pass_, h_efficiencyNonIsoTauET_EB_total_, h_efficiencyNonIsoTauET_EE_pass_, h_efficiencyNonIsoTauET_EE_total_, h_L1TauEtavsTauEta_, h_L1TauETvsTauET_EB_, h_L1TauETvsTauET_EB_EE_, h_L1TauETvsTauET_EE_, h_L1TauPhivsTauPhi_EB_, h_L1TauPhivsTauPhi_EB_EE_, h_L1TauPhivsTauPhi_EE_, h_resolutionTauET_EB_, h_resolutionTauET_EB_EE_, h_resolutionTauET_EE_, h_resolutionTauEta_, h_resolutionTauPhi_EB_, h_resolutionTauPhi_EB_EE_, h_resolutionTauPhi_EE_, edm::HandleBase::isValid(), m_BField, m_MaxHltTauDR, m_MaxL1tTauDR, m_MaxTauEta, m_TauL1tPairs, MetInputTag_, singleTopDQM_cfi::mets, mFieldToken_, MuonInputTag_, PDWG_BPHSkim_cff::muons, phi, BeamMonitor_cff::primaryVertex, DiDispStaMuonMonitor_cfi::pt, BXVector< T >::size(), stage2CaloLayer2TauToken_, AlCaHLTBitMon_QueryRunRegistry::string, metsig::tau, tauEfficiencyThresholds_, Tau3MuMonitor_cff::taus, theTauCollection_, remoteMonitoring_LASER_era2018_cfg::threshold, cond::impl::to_string(), B2GTnPMonitor_cfi::trigEvent, triggerEvent_, triggerResults_, bphysicsOniaDQM_cfi::vertex, and VtxInputTag_.

                                                                           {
   m_MaxTauEta = 2.1;
   m_MaxL1tTauDR = 0.5;
   m_MaxHltTauDR = 0.5;
 
   edm::Handle<reco::PFTauCollection> taus;
   e.getByToken(theTauCollection_, taus);
 
   if (!taus.isValid()) {
     edm::LogWarning("L1TTauOffline") << "invalid collection: reco::PFTauCollection " << std::endl;
     return;
   }
 
   edm::Handle<reco::MuonCollection> muons;
   e.getByToken(MuonInputTag_, muons);
 
   if (!muons.isValid()) {
     edm::LogWarning("L1TTauOffline") << "invalid collection: reco::MuonCollection " << std::endl;
     return;
   }
 
   edm::Handle<reco::BeamSpot> beamSpot;
   e.getByToken(BsInputTag_, beamSpot);
 
   if (!beamSpot.isValid()) {
     edm::LogWarning("L1TTauOffline") << "invalid collection: reco::BeamSpot " << std::endl;
     return;
   }
 
   edm::Handle<reco::VertexCollection> vertex;
   e.getByToken(VtxInputTag_, vertex);
 
   if (!vertex.isValid()) {
     edm::LogWarning("L1TTauOffline") << "invalid collection: reco::VertexCollection " << std::endl;
     return;
   }
 
   edm::Handle<l1t::TauBxCollection> l1tCands;
   e.getByToken(stage2CaloLayer2TauToken_, l1tCands);
 
   if (!l1tCands.isValid()) {
     edm::LogWarning("L1TTauOffline") << "invalid collection: l1t::TauBxCollection " << std::endl;
     return;
   }
 
   edm::Handle<edm::TriggerResults> trigResults;
   e.getByToken(triggerResults_, trigResults);
 
   if (!trigResults.isValid()) {
     edm::LogWarning("L1TTauOffline") << "invalid collection: edm::TriggerResults " << std::endl;
     return;
   }
 
   edm::Handle<trigger::TriggerEvent> trigEvent;
   e.getByToken(triggerEvent_, trigEvent);
 
   if (!trigEvent.isValid()) {
     edm::LogWarning("L1TTauOffline") << "invalid collection: trigger::TriggerEvent " << std::endl;
     return;
   }
 
   edm::Handle<reco::PFMETCollection> mets;
   e.getByToken(MetInputTag_, mets);
 
   if (!mets.isValid()) {
     edm::LogWarning("L1TTauOffline") << "invalid collection: reco::PFMETCollection " << std::endl;
     return;
   }
 
   m_BField = eSetup.getHandle(mFieldToken_);
 
   const reco::Vertex primaryVertex = getPrimaryVertex(vertex, beamSpot);
 
   getTightMuons(muons, mets, primaryVertex, trigEvent);
   getProbeTaus(e, taus, muons, primaryVertex);
   getTauL1tPairs(l1tCands);
 
   vector<l1t::Tau> l1tContainer;
   l1tContainer.reserve(l1tCands->size() + 1);
 
   for (auto tau = l1tCands->begin(0); tau != l1tCands->end(0); ++tau) {
     l1tContainer.push_back(*tau);
   }
 
   for (auto tauL1tPairsIt = m_TauL1tPairs.begin(); tauL1tPairsIt != m_TauL1tPairs.end(); ++tauL1tPairsIt) {
     float eta = tauL1tPairsIt->eta();
     float phi = tauL1tPairsIt->phi();
     float pt = tauL1tPairsIt->pt();
 
     // unmatched gmt cands have l1tPt = -1.
     float l1tPt = tauL1tPairsIt->l1tPt();
 
     int counter = 0;
 
     for (auto threshold : tauEfficiencyThresholds_) {
       std::string str_threshold = std::to_string(threshold);
 
       int l1tPtCut = threshold;
       bool l1tAboveCut = (l1tPt >= l1tPtCut);
 
       stringstream ptCutToTag;
       ptCutToTag << l1tPtCut;
       string ptTag = ptCutToTag.str();
 
       if (fabs(eta) < m_MaxTauEta) {
         if (counter == 0) {
           if (fabs(eta) < 1.5) {
             h_L1TauETvsTauET_EB_->Fill(pt, l1tPt);
             h_L1TauPhivsTauPhi_EB_->Fill(phi, tauL1tPairsIt->l1tPhi());
             h_resolutionTauET_EB_->Fill((l1tPt - pt) / pt);
             h_resolutionTauPhi_EB_->Fill(tauL1tPairsIt->l1tPhi() - phi);
           } else {
             h_L1TauETvsTauET_EE_->Fill(pt, l1tPt);
             h_L1TauPhivsTauPhi_EE_->Fill(phi, tauL1tPairsIt->l1tPhi());
             h_resolutionTauET_EE_->Fill((l1tPt - pt) / pt);
             h_resolutionTauPhi_EE_->Fill(tauL1tPairsIt->l1tPhi() - phi);
           }
           h_L1TauETvsTauET_EB_EE_->Fill(pt, l1tPt);
           h_L1TauPhivsTauPhi_EB_EE_->Fill(phi, tauL1tPairsIt->l1tPhi());
           h_L1TauEtavsTauEta_->Fill(eta, tauL1tPairsIt->l1tEta());
           h_resolutionTauET_EB_EE_->Fill((l1tPt - pt) / pt);
           h_resolutionTauPhi_EB_EE_->Fill(tauL1tPairsIt->l1tPhi() - phi);
           h_resolutionTauEta_->Fill(tauL1tPairsIt->l1tEta() - eta);
 
           ++counter;
         }
 
         if (fabs(eta) < 1.5) {
           h_efficiencyNonIsoTauET_EB_total_[threshold]->Fill(pt);
           h_efficiencyIsoTauET_EB_total_[threshold]->Fill(pt);
         } else {
           h_efficiencyNonIsoTauET_EE_total_[threshold]->Fill(pt);
           h_efficiencyIsoTauET_EE_total_[threshold]->Fill(pt);
         }
         h_efficiencyNonIsoTauET_EB_EE_total_[threshold]->Fill(pt);
         h_efficiencyIsoTauET_EB_EE_total_[threshold]->Fill(pt);
 
         if (l1tAboveCut) {
           if (fabs(eta) < 1.5)
             h_efficiencyNonIsoTauET_EB_pass_[threshold]->Fill(pt);
           else
             h_efficiencyNonIsoTauET_EE_pass_[threshold]->Fill(pt);
           h_efficiencyNonIsoTauET_EB_EE_pass_[threshold]->Fill(pt);
 
           if (tauL1tPairsIt->l1tIso() > 0.5) {
             if (fabs(eta) < 1.5)
               h_efficiencyIsoTauET_EB_pass_[threshold]->Fill(pt);
             else
               h_efficiencyIsoTauET_EE_pass_[threshold]->Fill(pt);
             h_efficiencyIsoTauET_EB_EE_pass_[threshold]->Fill(pt);
           }
         }
       }
     }
   }  // loop over tau-L1 pairs
 }

◆ bookHistograms()

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 139 of file L1TTauOffline.cc.

References bookTauHistos(), newFWLiteAna::found, createfilelist::int, m_hltConfig, m_trigIndices, HLTConfigProvider::size(), AlCaHLTBitMon_QueryRunRegistry::string, HLTConfigProvider::triggerName(), and triggerPath_.

                                                                                                 {
   edm::LogInfo("L1TTauOffline") << "L1TTauOffline::bookHistograms" << std::endl;
 
   // book at beginRun
   bookTauHistos(ibooker);
 
   for (auto trigNamesIt = triggerPath_.begin(); trigNamesIt != triggerPath_.end(); trigNamesIt++) {
     std::string tNameTmp = (*trigNamesIt);
     std::string tNamePattern = "";
     std::size_t found0 = tNameTmp.find('*');
     if (found0 != std::string::npos)
       tNamePattern = tNameTmp.substr(0, tNameTmp.size() - 1);
     else
       tNamePattern = tNameTmp;
 
     int tIndex = -1;
 
     for (unsigned ipath = 0; ipath < m_hltConfig.size(); ++ipath) {
       std::string tmpName = m_hltConfig.triggerName(ipath);
 
       std::size_t found = tmpName.find(tNamePattern);
       if (found != std::string::npos) {
         tIndex = int(ipath);
         m_trigIndices.push_back(tIndex);
       }
     }
   }
 }

◆ bookTauHistos()

void L1TTauOffline::bookTauHistos ( DQMStore::IBooker & ibooker )

private

Definition at line 334 of file L1TTauOffline.cc.

Referenced by bookHistograms().

                                                           {
   ibooker.cd();
   ibooker.setCurrentFolder(histFolder_);
   dqmoffline::l1t::HistDefinition nVertexDef = histDefinitions_[PlotConfig::nVertex];
   h_nVertex_ = ibooker.book1D(nVertexDef.name, nVertexDef.title, nVertexDef.nbinsX, nVertexDef.xmin, nVertexDef.xmax);
   h_tagAndProbeMass_ = ibooker.book1D("tagAndProbeMass", "Invariant mass of tag & probe pair", 100, 40, 140);
 
   dqmoffline::l1t::HistDefinition templateETvsET = histDefinitions_[PlotConfig::ETvsET];
   h_L1TauETvsTauET_EB_ = ibooker.book2D("L1TauETvsTauET_EB",
                                         "L1 Tau E_{T} vs PFTau E_{T} (EB); PFTau E_{T} (GeV); L1 Tau E_{T} (GeV)",
                                         templateETvsET.nbinsX,
                                         &templateETvsET.binsX[0],
                                         templateETvsET.nbinsY,
                                         &templateETvsET.binsY[0]);
   h_L1TauETvsTauET_EE_ = ibooker.book2D("L1TauETvsTauET_EE",
                                         "L1 Tau E_{T} vs PFTau E_{T} (EE); PFTau E_{T} (GeV); L1 Tau E_{T} (GeV)",
                                         templateETvsET.nbinsX,
                                         &templateETvsET.binsX[0],
                                         templateETvsET.nbinsY,
                                         &templateETvsET.binsY[0]);
   h_L1TauETvsTauET_EB_EE_ = ibooker.book2D("L1TauETvsTauET_EB_EE",
                                            "L1 Tau E_{T} vs PFTau E_{T} (EB+EE); PFTau E_{T} (GeV); L1 Tau E_{T} (GeV)",
                                            templateETvsET.nbinsX,
                                            &templateETvsET.binsX[0],
                                            templateETvsET.nbinsY,
                                            &templateETvsET.binsY[0]);
 
   dqmoffline::l1t::HistDefinition templatePHIvsPHI = histDefinitions_[PlotConfig::PHIvsPHI];
   h_L1TauPhivsTauPhi_EB_ =
       ibooker.book2D("L1TauPhivsTauPhi_EB",
                      "#phi_{tau}^{L1} vs #phi_{tau}^{offline} (EB); #phi_{tau}^{offline}; #phi_{tau}^{L1}",
                      templatePHIvsPHI.nbinsX,
                      templatePHIvsPHI.xmin,
                      templatePHIvsPHI.xmax,
                      templatePHIvsPHI.nbinsY,
                      templatePHIvsPHI.ymin,
                      templatePHIvsPHI.ymax);
   h_L1TauPhivsTauPhi_EE_ =
       ibooker.book2D("L1TauPhivsTauPhi_EE",
                      "#phi_{tau}^{L1} vs #phi_{tau}^{offline} (EE); #phi_{tau}^{offline}; #phi_{tau}^{L1}",
                      templatePHIvsPHI.nbinsX,
                      templatePHIvsPHI.xmin,
                      templatePHIvsPHI.xmax,
                      templatePHIvsPHI.nbinsY,
                      templatePHIvsPHI.ymin,
                      templatePHIvsPHI.ymax);
   h_L1TauPhivsTauPhi_EB_EE_ =
       ibooker.book2D("L1TauPhivsTauPhi_EB_EE",
                      "#phi_{tau}^{L1} vs #phi_{tau}^{offline} (EB+EE); #phi_{tau}^{offline}; #phi_{tau}^{L1}",
                      templatePHIvsPHI.nbinsX,
                      templatePHIvsPHI.xmin,
                      templatePHIvsPHI.xmax,
                      templatePHIvsPHI.nbinsY,
                      templatePHIvsPHI.ymin,
                      templatePHIvsPHI.ymax);
 
   h_L1TauEtavsTauEta_ =
       ibooker.book2D("L1TauEtavsTauEta", "L1 Tau #eta vs PFTau #eta; PFTau #eta; L1 Tau #eta", 100, -3, 3, 100, -3, 3);
 
   // tau resolutions
   h_resolutionTauET_EB_ = ibooker.book1D(
       "resolutionTauET_EB", "tau ET resolution (EB); (L1 Tau E_{T} - PFTau E_{T})/PFTau E_{T}; events", 50, -1, 1.5);
   h_resolutionTauET_EE_ = ibooker.book1D(
       "resolutionTauET_EE", "tau ET resolution (EE); (L1 Tau E_{T} - PFTau E_{T})/PFTau E_{T}; events", 50, -1, 1.5);
   h_resolutionTauET_EB_EE_ =
       ibooker.book1D("resolutionTauET_EB_EE",
                      "tau ET resolution (EB+EE); (L1 Tau E_{T} - PFTau E_{T})/PFTau E_{T}; events",
                      50,
                      -1,
                      1.5);
 
   h_resolutionTauPhi_EB_ = ibooker.book1D("resolutionTauPhi_EB",
                                           "#phi_{tau} resolution (EB); #phi_{tau}^{L1} - #phi_{tau}^{offline}; events",
                                           120,
                                           -0.3,
                                           0.3);
   h_resolutionTauPhi_EE_ = ibooker.book1D(
       "resolutionTauPhi_EE", "tau #phi resolution (EE); #phi_{tau}^{L1} - #phi_{tau}^{offline}; events", 120, -0.3, 0.3);
   h_resolutionTauPhi_EB_EE_ =
       ibooker.book1D("resolutionTauPhi_EB_EE",
                      "tau #phi resolution (EB+EE); #phi_{tau}^{L1} - #phi_{tau}^{offline}; events",
                      120,
                      -0.3,
                      0.3);
 
   h_resolutionTauEta_ =
       ibooker.book1D("resolutionTauEta", "tau #eta resolution  (EB); L1 Tau #eta - PFTau #eta; events", 120, -0.3, 0.3);
 
   // tau turn-ons
   ibooker.setCurrentFolder(efficiencyFolder_);
   std::vector<float> tauBins(tauEfficiencyBins_.begin(), tauEfficiencyBins_.end());
   int nBins = tauBins.size() - 1;
   float* tauBinArray = &(tauBins[0]);
 
   for (auto threshold : tauEfficiencyThresholds_) {
     std::string str_threshold = std::to_string(int(threshold));
     h_efficiencyIsoTauET_EB_pass_[threshold] =
         ibooker.book1D("efficiencyIsoTauET_EB_threshold_" + str_threshold + "_Num",
                        "iso tau efficiency (EB); PFTau E_{T} (GeV); events",
                        nBins,
                        tauBinArray);
     h_efficiencyIsoTauET_EE_pass_[threshold] =
         ibooker.book1D("efficiencyIsoTauET_EE_threshold_" + str_threshold + "_Num",
                        "iso tau efficiency (EE); PFTau E_{T} (GeV); events",
                        nBins,
                        tauBinArray);
     h_efficiencyIsoTauET_EB_EE_pass_[threshold] =
         ibooker.book1D("efficiencyIsoTauET_EB_EE_threshold_" + str_threshold + "_Num",
                        "iso tau efficiency (EB+EE); PFTau E_{T} (GeV); events",
                        nBins,
                        tauBinArray);
 
     h_efficiencyIsoTauET_EB_total_[threshold] =
         ibooker.book1D("efficiencyIsoTauET_EB_threshold_" + str_threshold + "_Den",
                        "iso tau efficiency (EB); PFTau E_{T} (GeV); events",
                        nBins,
                        tauBinArray);
     h_efficiencyIsoTauET_EE_total_[threshold] =
         ibooker.book1D("efficiencyIsoTauET_EE_threshold_" + str_threshold + "_Den",
                        "iso tau efficiency (EE); PFTau E_{T} (GeV); events",
                        nBins,
                        tauBinArray);
     h_efficiencyIsoTauET_EB_EE_total_[threshold] =
         ibooker.book1D("efficiencyIsoTauET_EB_EE_threshold_" + str_threshold + "_Den",
                        "iso tau efficiency (EB+EE); PFTau E_{T} (GeV); events",
                        nBins,
                        tauBinArray);
 
     // non iso
     h_efficiencyNonIsoTauET_EB_pass_[threshold] =
         ibooker.book1D("efficiencyNonIsoTauET_EB_threshold_" + str_threshold + "_Num",
                        "inclusive tau efficiency (EB); PFTau E_{T} (GeV); events",
                        nBins,
                        tauBinArray);
     h_efficiencyNonIsoTauET_EE_pass_[threshold] =
         ibooker.book1D("efficiencyNonIsoTauET_EE_threshold_" + str_threshold + "_Num",
                        "inclusive tau efficiency (EE); PFTau E_{T} (GeV); events",
                        nBins,
                        tauBinArray);
     h_efficiencyNonIsoTauET_EB_EE_pass_[threshold] =
         ibooker.book1D("efficiencyNonIsoTauET_EB_EE_threshold_" + str_threshold + "_Num",
                        "inclusive tau efficiency (EB+EE); PFTau E_{T} (GeV); events",
                        nBins,
                        tauBinArray);
 
     h_efficiencyNonIsoTauET_EB_total_[threshold] =
         ibooker.book1D("efficiencyNonIsoTauET_EB_threshold_" + str_threshold + "_Den",
                        "inclusive tau efficiency (EB); PFTau E_{T} (GeV); events",
                        nBins,
                        tauBinArray);
     h_efficiencyNonIsoTauET_EE_total_[threshold] =
         ibooker.book1D("efficiencyNonIsoTauET_EE_threshold_" + str_threshold + "_Den",
                        "inclusive tau efficiency (EE); PFTau E_{T} (GeV); events",
                        nBins,
                        tauBinArray);
     h_efficiencyNonIsoTauET_EB_EE_total_[threshold] =
         ibooker.book1D("efficiencyNonIsoTauET_EB_EE_threshold_" + str_threshold + "_Den",
                        "inclusive tau efficiency (EB+EE); PFTau E_{T} (GeV); events",
                        nBins,
                        tauBinArray);
   }
 
   ibooker.cd();
 
   return;
 }

◆ calcDeltaPhi()

double L1TTauOffline::calcDeltaPhi	(	double	phi1,
		double	phi2
	)

private

◆ Distance()

double L1TTauOffline::Distance	(	const reco::Candidate &	c1,
		const reco::Candidate &	c2
	)

private

◆ DistancePhi()

double L1TTauOffline::DistancePhi	(	const reco::Candidate &	c1,
		const reco::Candidate &	c2
	)

private

◆ dqmBeginRun()

void L1TTauOffline::dqmBeginRun	(	const edm::Run &	run,
		const edm::EventSetup &	iSetup
	)

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 127 of file L1TTauOffline.cc.

References HLTConfigProvider::init(), m_hltConfig, writedatasetfile::run, and trigProcess_.

 {
   bool changed = true;
   m_hltConfig.init(run, iSetup, trigProcess_, changed);
 
   edm::LogInfo("L1TTauOffline") << "L1TTauOffline::beginRun" << std::endl;
 }

◆ getPrimaryVertex()

const reco::Vertex L1TTauOffline::getPrimaryVertex	(	edm::Handle< reco::VertexCollection > const &	vertex,
		edm::Handle< reco::BeamSpot > const &	beamSpot
	)

protected

Definition at line 501 of file L1TTauOffline.cc.

References pwdgSkimBPark_cfi::beamSpot, reco::Vertex::position(), BeamMonitor_cff::primaryVertex, and bphysicsOniaDQM_cfi::vertex.

Referenced by analyze().

                                                                                             {
   reco::Vertex::Point posVtx;
   reco::Vertex::Error errVtx;
 
   bool hasPrimaryVertex = false;
 
   if (vertex.isValid()) {
     for (auto vertexIt = vertex->begin(); vertexIt != vertex->end(); ++vertexIt) {
       if (vertexIt->isValid() && !vertexIt->isFake()) {
         posVtx = vertexIt->position();
         errVtx = vertexIt->error();
         hasPrimaryVertex = true;
         break;
       }
     }
   }
 
   if (!hasPrimaryVertex) {
     posVtx = beamSpot->position();
     errVtx(0, 0) = beamSpot->BeamWidthX();
     errVtx(1, 1) = beamSpot->BeamWidthY();
     errVtx(2, 2) = beamSpot->sigmaZ();
   }
 
   const reco::Vertex primaryVertex(posVtx, errVtx);
 
   return primaryVertex;
 }

◆ getProbeTaus()

void L1TTauOffline::getProbeTaus	(	const edm::Event &	e,
		edm::Handle< reco::PFTauCollection > const &	taus,
		edm::Handle< reco::MuonCollection > const &	muons,
		const reco::Vertex &	vertex
	)

protected

Definition at line 637 of file L1TTauOffline.cc.

References AntiEleInputTag_, AntiEleWP_, AntiEleWPIndex_, AntiMuInputTag_, AntiMuWP_, AntiMuWPIndex_, comb3TInputTag_, comb3TWP_, comb3TWPIndex_, DecayModeFindingInputTag_, HCALHighEnergyHPDFilter_cfi::energy, PVValHelper::eta, mps_fire::i, iEvent, edm::HandleBase::isValid(), m_AntiEleExists, m_ProbeTaus, m_TightMuons, edm::parameterSet(), phi, phID_, edm::HandleBase::provenance(), DiDispStaMuonMonitor_cfi::pt, edm::Provenance::stable(), Tau3MuMonitor_cff::taus, parallelization::uint, and PDWG_TauSkim_cff::workingPoints.

Referenced by analyze().

                                                            {
   m_ProbeTaus.clear();
 
   edm::Handle<reco::TauDiscriminatorContainer> antimu;
   iEvent.getByToken(AntiMuInputTag_, antimu);
   if (!antimu.isValid()) {
     edm::LogWarning("L1TTauOffline") << "invalid collection: reco::PFTauDiscriminator  anti-mu";
     return;
   }
 
   edm::Handle<reco::PFTauDiscriminator> dmf;
   iEvent.getByToken(DecayModeFindingInputTag_, dmf);
   if (!dmf.isValid()) {
     edm::LogWarning("L1TTauOffline") << "invalid collection: reco::PFTauDiscriminator decay mode finding input";
     return;
   }
   m_AntiEleExists = true;
 
   edm::Handle<reco::TauDiscriminatorContainer> antiele;
   iEvent.getByToken(AntiEleInputTag_, antiele);
   if (!antiele.isValid()) {
     //V.M. 16.3.2023. Bypassing the return option for now, as the anti-ele discr. is not available.
     m_AntiEleExists = false;
   }
 
   edm::Handle<reco::TauDiscriminatorContainer> comb3T;
   iEvent.getByToken(comb3TInputTag_, comb3T);
   if (!comb3T.isValid()) {
     edm::LogWarning("L1TTauOffline") << "invalid collection: reco::PFTauDiscriminator  comb3T";
     return;
   }
 
   if (!m_TightMuons.empty()) {
     TLorentzVector mymu;
     mymu.SetPtEtaPhiE(m_TightMuons[0]->pt(), m_TightMuons[0]->eta(), m_TightMuons[0]->phi(), m_TightMuons[0]->energy());
     int iTau = 0;
 
     // load indices from input provenance config if process history changed, in particular for the first event
     if (iEvent.processHistoryID() != phID_) {
       phID_ = iEvent.processHistoryID();
       {
         const edm::Provenance* prov = antimu.provenance();
         const std::vector<edm::ParameterSet> psetsFromProvenance =
             edm::parameterSet(prov->stable(), iEvent.processHistory())
                 .getParameter<std::vector<edm::ParameterSet>>("IDWPdefinitions");
         for (uint i = 0; i < psetsFromProvenance.size(); i++) {
           if (psetsFromProvenance[i].getParameter<std::string>("IDname") == AntiMuWP_)
             AntiMuWPIndex_ = i;
         }
       }
       {
         if (m_AntiEleExists) {
           const edm::Provenance* prov = antiele.provenance();
           const std::vector<std::string> psetsFromProvenance =
               edm::parameterSet(prov->stable(), iEvent.processHistory())
                   .getParameter<std::vector<std::string>>("workingPoints");
           for (unsigned int i = 0; i < psetsFromProvenance.size(); i++) {
             if (psetsFromProvenance[i] == AntiEleWP_)
               AntiEleWPIndex_ = i;
           }
         }
       }
       {
         const edm::Provenance* prov = comb3T.provenance();
         const std::vector<edm::ParameterSet> psetsFromProvenance =
             edm::parameterSet(prov->stable(), iEvent.processHistory())
                 .getParameter<std::vector<edm::ParameterSet>>("IDWPdefinitions");
         for (uint i = 0; i < psetsFromProvenance.size(); i++) {
           if (psetsFromProvenance[i].getParameter<std::string>("IDname") == comb3TWP_)
             comb3TWPIndex_ = i;
         }
       }
     }
 
     for (auto tauIt = taus->begin(); tauIt != taus->end(); ++tauIt, ++iTau) {
       reco::PFTauRef tauCandidate(taus, iTau);
       TLorentzVector mytau;
       mytau.SetPtEtaPhiE(tauIt->pt(), tauIt->eta(), tauIt->phi(), tauIt->energy());
 
       if ((*antimu)[tauCandidate].workingPoints.empty()) {
         edm::LogWarning("L1TTauOffline") << "This offline tau has no antimu discriminator, skipping";
         continue;
       }
       if (m_AntiEleExists) {
         if ((*antiele)[tauCandidate].workingPoints.empty()) {
           edm::LogWarning("L1TTauOffline") << "This offline tau has no antiele discriminator, skipping";
           continue;
         }
       }
       if ((*comb3T)[tauCandidate].workingPoints.empty()) {
         edm::LogWarning("L1TTauOffline") << "This offline tau has no comb3T discriminator, skipping";
         continue;
       }
       bool antiele_condition = true;
 
       if (m_AntiEleExists)
         antiele_condition = (*antiele)[tauCandidate].workingPoints[AntiEleWPIndex_];
 
       if (fabs(tauIt->charge()) == 1 && fabs(tauIt->eta()) < 2.1 && tauIt->pt() > 20 &&
           (*antimu)[tauCandidate].workingPoints[AntiMuWPIndex_] && antiele_condition && (*dmf)[tauCandidate] > 0.5 &&
           (*comb3T)[tauCandidate].workingPoints[comb3TWPIndex_]) {
         if (mymu.DeltaR(mytau) > 0.5 && (mymu + mytau).M() > 40 && (mymu + mytau).M() < 80 &&
             m_TightMuons[0]->charge() * tauIt->charge() < 0) {
           m_ProbeTaus.push_back(&(*tauIt));
         }
       }
     }
   }
 }

◆ getTauL1tPairs()

void L1TTauOffline::getTauL1tPairs ( edm::Handle< l1t::TauBxCollection > const & l1tCands )

protected

Definition at line 569 of file L1TTauOffline.cc.

References BXVector< T >::begin(), TauL1TPair::dR(), BXVector< T >::end(), TauL1TPair::l1tPt(), m_MaxL1tTauDR, m_ProbeTaus, m_TauL1tPairs, BXVector< T >::size(), and metsig::tau.

Referenced by analyze().

                                                                                 {
   m_TauL1tPairs.clear();
 
   vector<l1t::Tau> l1tContainer;
   l1tContainer.reserve(l1tCands->size() + 1);
 
   for (auto tau = l1tCands->begin(0); tau != l1tCands->end(0); ++tau) {
     l1tContainer.push_back(*tau);
   }
 
   for (auto probeTauIt = m_ProbeTaus.begin(); probeTauIt != m_ProbeTaus.end(); ++probeTauIt) {
     TauL1TPair pairBestCand((*probeTauIt), nullptr);
 
     for (auto l1tIt = l1tContainer.begin(); l1tIt != l1tContainer.end(); ++l1tIt) {
       TauL1TPair pairTmpCand((*probeTauIt), &(*l1tIt));
 
       if (pairTmpCand.dR() < m_MaxL1tTauDR && pairTmpCand.l1tPt() > pairBestCand.l1tPt())
         pairBestCand = pairTmpCand;
     }
 
     m_TauL1tPairs.push_back(pairBestCand);
   }
 }

◆ getTightMuons()

void L1TTauOffline::getTightMuons	(	edm::Handle< reco::MuonCollection > const &	muons,
		edm::Handle< reco::PFMETCollection > const &	mets,
		const reco::Vertex &	vertex,
		edm::Handle< trigger::TriggerEvent > const &	trigEvent
	)

protected

Definition at line 593 of file L1TTauOffline.cc.

References muon::isLooseMuon(), m_TightMuons, matchHlt(), SiStripPI::max, singleTopDQM_cfi::mets, TtSemiLepEvtBuilder_cfi::mt, PDWG_BPHSkim_cff::muons, funct::pow(), reco::LeafCandidate::pt(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), mathSSE::sqrt(), and B2GTnPMonitor_cfi::trigEvent.

Referenced by analyze().

                                                                                    {
   m_TightMuons.clear();
 
   const reco::PFMET* pfmet = nullptr;
   pfmet = &(mets->front());
 
   int nb_mu = 0;
 
   for (auto muonIt2 = muons->begin(); muonIt2 != muons->end(); ++muonIt2) {
     if (fabs(muonIt2->eta()) < 2.4 && muonIt2->pt() > 10 && muon::isLooseMuon((*muonIt2)) &&
         (muonIt2->pfIsolationR04().sumChargedHadronPt +
          max(muonIt2->pfIsolationR04().sumNeutralHadronEt + muonIt2->pfIsolationR04().sumPhotonEt -
                  0.5 * muonIt2->pfIsolationR04().sumPUPt,
              0.0)) /
                 muonIt2->pt() <
             0.3) {
       ++nb_mu;
     }
   }
   bool foundTightMu = false;
   for (auto muonIt = muons->begin(); muonIt != muons->end(); ++muonIt) {
     if (!matchHlt(trigEvent, &(*muonIt)))
       continue;
     float muiso = (muonIt->pfIsolationR04().sumChargedHadronPt +
                    max(muonIt->pfIsolationR04().sumNeutralHadronEt + muonIt->pfIsolationR04().sumPhotonEt -
                            0.5 * muonIt->pfIsolationR04().sumPUPt,
                        0.0)) /
                   muonIt->pt();
 
     if (muiso < 0.1 && nb_mu < 2 && !foundTightMu && fabs(muonIt->eta()) < 2.1 && muonIt->pt() > 24 &&
         muon::isLooseMuon((*muonIt))) {
       float mt = sqrt(pow(muonIt->pt() + pfmet->pt(), 2) - pow(muonIt->px() + pfmet->px(), 2) -
                       pow(muonIt->py() + pfmet->py(), 2));
       if (mt < 30) {
         m_TightMuons.push_back(&(*muonIt));
         foundTightMu = true;
       }
     }
   }
 }

◆ matchHlt()

bool L1TTauOffline::matchHlt	(	edm::Handle< trigger::TriggerEvent > const &	triggerEvent,
		const reco::Muon *	muon
	)

protected

Definition at line 531 of file L1TTauOffline.cc.

References PbPb_ZMuSkimMuonDPG_cff::deltaR, cuy::ii, ProducerED_cfi::InputTag, m_hltConfig, m_MaxHltTauDR, m_trigIndices, HLTConfigProvider::moduleLabels(), PDWG_DiPhoton_SD_cff::triggerEvent, TriggerAnalyzer::trigObjs, and trigProcess_.

Referenced by getTightMuons().

                                                                                                      {
   double matchDeltaR = 9999;
 
   trigger::TriggerObjectCollection trigObjs = triggerEvent->getObjects();
 
   for (auto& trigIndex : m_trigIndices) {
     const vector<string> moduleLabels(m_hltConfig.moduleLabels(trigIndex));
     // V.M. 2023.03.08 Same issue as in the L1TMuonDQMOffline.cc - some modules are behind hltBoolEnd, but we want the last one before the "hltBoolEnd"
     int moduleIndex = 999999;
     for (int ii = moduleLabels.size() - 1; ii > 0; ii--) {
       if (moduleLabels[ii] == "hltBoolEnd") {
         moduleIndex = ii - 1;
         break;
       }
     }
     if (moduleIndex == 999999) {
       edm::LogError("L1TMuonDQMOffline") << "Found no module label in trigger " << trigIndex;
       continue;
     }
     const unsigned hltFilterIndex = triggerEvent->filterIndex(InputTag(moduleLabels[moduleIndex], "", trigProcess_));
 
     if (hltFilterIndex < triggerEvent->sizeFilters()) {
       const Keys triggerKeys(triggerEvent->filterKeys(hltFilterIndex));
       const Vids triggerVids(triggerEvent->filterIds(hltFilterIndex));
 
       const unsigned nTriggers = triggerVids.size();
       for (size_t iTrig = 0; iTrig < nTriggers; ++iTrig) {
         const TriggerObject trigObject = trigObjs[triggerKeys[iTrig]];
 
         double dRtmp = deltaR((*muon), trigObject);
         if (dRtmp < matchDeltaR)
           matchDeltaR = dRtmp;
       }
     }
   }
   return (matchDeltaR < m_MaxHltTauDR);
 }

◆ normalise2DHistogramsToBinArea()

void L1TTauOffline::normalise2DHistogramsToBinArea ( )

private

Definition at line 750 of file L1TTauOffline.cc.

References h, h_L1TauEtavsTauEta_, h_L1TauETvsTauET_EB_, h_L1TauETvsTauET_EB_EE_, h_L1TauETvsTauET_EE_, h_L1TauPhivsTauPhi_EB_, h_L1TauPhivsTauPhi_EB_EE_, h_L1TauPhivsTauPhi_EE_, and LaserDQM_cfi::mon.

                                                    {
   std::vector<MonitorElement*> monElementstoNormalize = {h_L1TauETvsTauET_EB_,
                                                          h_L1TauETvsTauET_EE_,
                                                          h_L1TauETvsTauET_EB_EE_,
                                                          h_L1TauPhivsTauPhi_EB_,
                                                          h_L1TauPhivsTauPhi_EE_,
                                                          h_L1TauPhivsTauPhi_EB_EE_,
                                                          h_L1TauEtavsTauEta_};
 
   for (auto mon : monElementstoNormalize) {
     if (mon != nullptr) {
       auto h = mon->getTH2F();
       if (h != nullptr) {
         h->Scale(1, "width");
       }
     }
   }
 }

Member Data Documentation

◆ AntiEleInputTag_

edm::EDGetTokenT<reco::TauDiscriminatorContainer> L1TTauOffline::AntiEleInputTag_

private

Definition at line 135 of file L1TTauOffline.h.

Referenced by getProbeTaus().

◆ AntiEleWP_

std::string L1TTauOffline::AntiEleWP_

private

Definition at line 136 of file L1TTauOffline.h.

Referenced by getProbeTaus().

◆ AntiEleWPIndex_

int L1TTauOffline::AntiEleWPIndex_

private

Definition at line 137 of file L1TTauOffline.h.

Referenced by getProbeTaus().

◆ AntiMuInputTag_

edm::EDGetTokenT<reco::TauDiscriminatorContainer> L1TTauOffline::AntiMuInputTag_

private

Definition at line 131 of file L1TTauOffline.h.

Referenced by getProbeTaus().

◆ AntiMuWP_

std::string L1TTauOffline::AntiMuWP_

private

Definition at line 133 of file L1TTauOffline.h.

Referenced by getProbeTaus().

◆ AntiMuWPIndex_

int L1TTauOffline::AntiMuWPIndex_

private

Definition at line 134 of file L1TTauOffline.h.

Referenced by getProbeTaus().

◆ BsInputTag_

edm::EDGetTokenT<reco::BeamSpot> L1TTauOffline::BsInputTag_

private

Definition at line 145 of file L1TTauOffline.h.

Referenced by analyze().

◆ comb3TInputTag_

edm::EDGetTokenT<reco::TauDiscriminatorContainer> L1TTauOffline::comb3TInputTag_

private

Definition at line 139 of file L1TTauOffline.h.

Referenced by getProbeTaus().

◆ comb3TWP_

std::string L1TTauOffline::comb3TWP_

private

Definition at line 140 of file L1TTauOffline.h.

Referenced by getProbeTaus().

◆ comb3TWPIndex_

int L1TTauOffline::comb3TWPIndex_

private

Definition at line 141 of file L1TTauOffline.h.

Referenced by getProbeTaus().

◆ DecayModeFindingInputTag_

edm::EDGetTokenT<reco::PFTauDiscriminator> L1TTauOffline::DecayModeFindingInputTag_

private

Definition at line 138 of file L1TTauOffline.h.

Referenced by getProbeTaus().

◆ efficiencyFolder_

std::string L1TTauOffline::efficiencyFolder_

private

Definition at line 151 of file L1TTauOffline.h.

Referenced by bookTauHistos().

◆ h_efficiencyIsoTauET_EB_EE_pass_

std::map<double, MonitorElement*> L1TTauOffline::h_efficiencyIsoTauET_EB_EE_pass_

private

Definition at line 207 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

◆ h_efficiencyIsoTauET_EB_EE_total_

std::map<double, MonitorElement*> L1TTauOffline::h_efficiencyIsoTauET_EB_EE_total_

private

Definition at line 217 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

◆ h_efficiencyIsoTauET_EB_pass_

std::map<double, MonitorElement*> L1TTauOffline::h_efficiencyIsoTauET_EB_pass_

private

Definition at line 205 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

◆ h_efficiencyIsoTauET_EB_total_

std::map<double, MonitorElement*> L1TTauOffline::h_efficiencyIsoTauET_EB_total_

private

Definition at line 215 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

◆ h_efficiencyIsoTauET_EE_pass_

std::map<double, MonitorElement*> L1TTauOffline::h_efficiencyIsoTauET_EE_pass_

private

Definition at line 206 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

◆ h_efficiencyIsoTauET_EE_total_

std::map<double, MonitorElement*> L1TTauOffline::h_efficiencyIsoTauET_EE_total_

private

Definition at line 216 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

◆ h_efficiencyNonIsoTauET_EB_EE_pass_

std::map<double, MonitorElement*> L1TTauOffline::h_efficiencyNonIsoTauET_EB_EE_pass_

private

Definition at line 211 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

◆ h_efficiencyNonIsoTauET_EB_EE_total_

std::map<double, MonitorElement*> L1TTauOffline::h_efficiencyNonIsoTauET_EB_EE_total_

private

Definition at line 221 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

◆ h_efficiencyNonIsoTauET_EB_pass_

std::map<double, MonitorElement*> L1TTauOffline::h_efficiencyNonIsoTauET_EB_pass_

private

Definition at line 209 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

◆ h_efficiencyNonIsoTauET_EB_total_

std::map<double, MonitorElement*> L1TTauOffline::h_efficiencyNonIsoTauET_EB_total_

private

Definition at line 219 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

◆ h_efficiencyNonIsoTauET_EE_pass_

std::map<double, MonitorElement*> L1TTauOffline::h_efficiencyNonIsoTauET_EE_pass_

private

Definition at line 210 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

◆ h_efficiencyNonIsoTauET_EE_total_

std::map<double, MonitorElement*> L1TTauOffline::h_efficiencyNonIsoTauET_EE_total_

private

Definition at line 220 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

◆ h_L1TauEtavsTauEta_

MonitorElement* L1TTauOffline::h_L1TauEtavsTauEta_

private

Definition at line 191 of file L1TTauOffline.h.

Referenced by analyze(), bookTauHistos(), and normalise2DHistogramsToBinArea().

◆ h_L1TauETvsTauET_EB_

MonitorElement* L1TTauOffline::h_L1TauETvsTauET_EB_

private

Definition at line 183 of file L1TTauOffline.h.

Referenced by analyze(), bookTauHistos(), and normalise2DHistogramsToBinArea().

◆ h_L1TauETvsTauET_EB_EE_

MonitorElement* L1TTauOffline::h_L1TauETvsTauET_EB_EE_

private

Definition at line 185 of file L1TTauOffline.h.

Referenced by analyze(), bookTauHistos(), and normalise2DHistogramsToBinArea().

◆ h_L1TauETvsTauET_EE_

MonitorElement* L1TTauOffline::h_L1TauETvsTauET_EE_

private

Definition at line 184 of file L1TTauOffline.h.

Referenced by analyze(), bookTauHistos(), and normalise2DHistogramsToBinArea().

◆ h_L1TauPhivsTauPhi_EB_

MonitorElement* L1TTauOffline::h_L1TauPhivsTauPhi_EB_

private

Definition at line 187 of file L1TTauOffline.h.

Referenced by analyze(), bookTauHistos(), and normalise2DHistogramsToBinArea().

◆ h_L1TauPhivsTauPhi_EB_EE_

MonitorElement* L1TTauOffline::h_L1TauPhivsTauPhi_EB_EE_

private

Definition at line 189 of file L1TTauOffline.h.

Referenced by analyze(), bookTauHistos(), and normalise2DHistogramsToBinArea().

◆ h_L1TauPhivsTauPhi_EE_

MonitorElement* L1TTauOffline::h_L1TauPhivsTauPhi_EE_

private

Definition at line 188 of file L1TTauOffline.h.

Referenced by analyze(), bookTauHistos(), and normalise2DHistogramsToBinArea().

◆ h_nVertex_

MonitorElement* L1TTauOffline::h_nVertex_

private

Definition at line 179 of file L1TTauOffline.h.

Referenced by bookTauHistos().

◆ h_resolutionTauET_EB_

MonitorElement* L1TTauOffline::h_resolutionTauET_EB_

private

Definition at line 194 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

◆ h_resolutionTauET_EB_EE_

MonitorElement* L1TTauOffline::h_resolutionTauET_EB_EE_

private

Definition at line 196 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

◆ h_resolutionTauET_EE_

MonitorElement* L1TTauOffline::h_resolutionTauET_EE_

private

Definition at line 195 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

◆ h_resolutionTauEta_

MonitorElement* L1TTauOffline::h_resolutionTauEta_

private

Definition at line 202 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

◆ h_resolutionTauPhi_EB_

MonitorElement* L1TTauOffline::h_resolutionTauPhi_EB_

private

Definition at line 198 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

◆ h_resolutionTauPhi_EB_EE_

MonitorElement* L1TTauOffline::h_resolutionTauPhi_EB_EE_

private

Definition at line 200 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

◆ h_resolutionTauPhi_EE_

MonitorElement* L1TTauOffline::h_resolutionTauPhi_EE_

private

Definition at line 199 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

◆ h_tagAndProbeMass_

MonitorElement* L1TTauOffline::h_tagAndProbeMass_

private

Definition at line 180 of file L1TTauOffline.h.

Referenced by bookTauHistos().

◆ histDefinitions_

dqmoffline::l1t::HistDefinitions L1TTauOffline::histDefinitions_

private

Definition at line 155 of file L1TTauOffline.h.

Referenced by bookTauHistos().

◆ histFolder_

std::string L1TTauOffline::histFolder_

private

Definition at line 150 of file L1TTauOffline.h.

Referenced by bookTauHistos().

◆ m_AntiEleExists

bool L1TTauOffline::m_AntiEleExists

private

Definition at line 176 of file L1TTauOffline.h.

Referenced by getProbeTaus().

◆ m_BField

edm::ESHandle<MagneticField> L1TTauOffline::m_BField

private

Definition at line 123 of file L1TTauOffline.h.

Referenced by analyze().

◆ m_hltConfig

HLTConfigProvider L1TTauOffline::m_hltConfig

private

Definition at line 121 of file L1TTauOffline.h.

Referenced by bookHistograms(), dqmBeginRun(), and matchHlt().

◆ m_L1tL1tTaus

BXVector<l1t::Tau> L1TTauOffline::m_L1tL1tTaus

private

Definition at line 164 of file L1TTauOffline.h.

◆ m_L1tPtCuts

std::vector<int> L1TTauOffline::m_L1tPtCuts

private

Definition at line 167 of file L1TTauOffline.h.

◆ m_L1tTaus

std::vector<l1t::TauBxCollection> L1TTauOffline::m_L1tTaus

private

Definition at line 162 of file L1TTauOffline.h.

◆ m_MaxHltTauDR

float L1TTauOffline::m_MaxHltTauDR

private

Definition at line 171 of file L1TTauOffline.h.

Referenced by analyze(), and matchHlt().

◆ m_MaxL1tTauDR

float L1TTauOffline::m_MaxL1tTauDR

private

Definition at line 170 of file L1TTauOffline.h.

Referenced by analyze(), and getTauL1tPairs().

◆ m_MaxTauEta

float L1TTauOffline::m_MaxTauEta

private

Definition at line 169 of file L1TTauOffline.h.

Referenced by analyze().

◆ m_ProbeTaus

std::vector<const reco::PFTau*> L1TTauOffline::m_ProbeTaus

private

Definition at line 158 of file L1TTauOffline.h.

Referenced by getProbeTaus(), and getTauL1tPairs().

◆ m_propagatorAlong

edm::ESHandle<Propagator> L1TTauOffline::m_propagatorAlong

private

Definition at line 124 of file L1TTauOffline.h.

◆ m_propagatorOpposite

edm::ESHandle<Propagator> L1TTauOffline::m_propagatorOpposite

private

Definition at line 125 of file L1TTauOffline.h.

◆ m_RecoRecoTaus

std::vector<reco::PFTau> L1TTauOffline::m_RecoRecoTaus

private

Definition at line 163 of file L1TTauOffline.h.

◆ m_RecoTaus

std::vector<reco::PFTauCollection> L1TTauOffline::m_RecoTaus

private

Definition at line 161 of file L1TTauOffline.h.

◆ m_TauL1tPairs

std::vector<TauL1TPair> L1TTauOffline::m_TauL1tPairs

private

Definition at line 159 of file L1TTauOffline.h.

Referenced by analyze(), and getTauL1tPairs().

◆ m_TightMuons

std::vector<const reco::Muon*> L1TTauOffline::m_TightMuons

private

Definition at line 157 of file L1TTauOffline.h.

Referenced by getProbeTaus(), and getTightMuons().

◆ m_trigIndices

std::vector<int> L1TTauOffline::m_trigIndices

private

Definition at line 173 of file L1TTauOffline.h.

Referenced by bookHistograms(), and matchHlt().

◆ MetInputTag_

edm::EDGetTokenT<reco::PFMETCollection> L1TTauOffline::MetInputTag_

private

Definition at line 143 of file L1TTauOffline.h.

Referenced by analyze().

◆ mFieldToken_

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> L1TTauOffline::mFieldToken_

private

Definition at line 132 of file L1TTauOffline.h.

Referenced by analyze(), and L1TTauOffline().

◆ MuonInputTag_

edm::EDGetTokenT<reco::MuonCollection> L1TTauOffline::MuonInputTag_

private

Definition at line 142 of file L1TTauOffline.h.

Referenced by analyze().

◆ phID_

edm::ProcessHistoryID L1TTauOffline::phID_

private

Definition at line 127 of file L1TTauOffline.h.

Referenced by getProbeTaus().

◆ PlotConfigNames

const std::map< std::string, unsigned int > L1TTauOffline::PlotConfigNames

static

Initial value:

= {

{"nVertex", PlotConfig::nVertex}, {"ETvsET", PlotConfig::ETvsET}, {"PHIvsPHI", PlotConfig::PHIvsPHI}}

Definition at line 87 of file L1TTauOffline.h.

◆ PVPoint_

math::XYZPoint L1TTauOffline::PVPoint_

private

Definition at line 119 of file L1TTauOffline.h.

◆ stage2CaloLayer2TauToken_

edm::EDGetTokenT<l1t::TauBxCollection> L1TTauOffline::stage2CaloLayer2TauToken_

private

Definition at line 152 of file L1TTauOffline.h.

Referenced by analyze().

◆ tauEfficiencyBins_

std::vector<double> L1TTauOffline::tauEfficiencyBins_

private

Definition at line 154 of file L1TTauOffline.h.

Referenced by bookTauHistos().

◆ tauEfficiencyThresholds_

std::vector<int> L1TTauOffline::tauEfficiencyThresholds_

private

Definition at line 153 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

◆ theTauCollection_

edm::EDGetTokenT<reco::PFTauCollection> L1TTauOffline::theTauCollection_

private

Definition at line 130 of file L1TTauOffline.h.

Referenced by analyze().

◆ triggerEvent_

edm::EDGetTokenT<trigger::TriggerEvent> L1TTauOffline::triggerEvent_

private

Definition at line 146 of file L1TTauOffline.h.

Referenced by analyze().

◆ triggerPath_

std::vector<std::string> L1TTauOffline::triggerPath_

private

Definition at line 149 of file L1TTauOffline.h.

Referenced by bookHistograms().

◆ triggerResults_

edm::EDGetTokenT<edm::TriggerResults> L1TTauOffline::triggerResults_

private

Definition at line 148 of file L1TTauOffline.h.

Referenced by analyze().

◆ trigProcess_

std::string L1TTauOffline::trigProcess_

private

Definition at line 147 of file L1TTauOffline.h.

Referenced by dqmBeginRun(), and matchHlt().

◆ VtxInputTag_

edm::EDGetTokenT<reco::VertexCollection> L1TTauOffline::VtxInputTag_

private

Definition at line 144 of file L1TTauOffline.h.

Referenced by analyze().

Public Types

Public Member Functions

Static Public Attributes

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

◆ PlotConfig

Constructor & Destructor Documentation

◆ L1TTauOffline()

◆ ~L1TTauOffline()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ bookTauHistos()

◆ calcDeltaPhi()

◆ Distance()

◆ DistancePhi()

◆ dqmBeginRun()

◆ getPrimaryVertex()

◆ getProbeTaus()

◆ getTauL1tPairs()

◆ getTightMuons()

◆ matchHlt()

◆ normalise2DHistogramsToBinArea()

Member Data Documentation

◆ AntiEleInputTag_

◆ AntiEleWP_

◆ AntiEleWPIndex_

◆ AntiMuInputTag_

◆ AntiMuWP_

◆ AntiMuWPIndex_

◆ BsInputTag_

◆ comb3TInputTag_

◆ comb3TWP_

◆ comb3TWPIndex_

◆ DecayModeFindingInputTag_

◆ efficiencyFolder_

◆ h_efficiencyIsoTauET_EB_EE_pass_

◆ h_efficiencyIsoTauET_EB_EE_total_

◆ h_efficiencyIsoTauET_EB_pass_

◆ h_efficiencyIsoTauET_EB_total_

◆ h_efficiencyIsoTauET_EE_pass_

◆ h_efficiencyIsoTauET_EE_total_

◆ h_efficiencyNonIsoTauET_EB_EE_pass_

◆ h_efficiencyNonIsoTauET_EB_EE_total_

◆ h_efficiencyNonIsoTauET_EB_pass_

◆ h_efficiencyNonIsoTauET_EB_total_

◆ h_efficiencyNonIsoTauET_EE_pass_

◆ h_efficiencyNonIsoTauET_EE_total_

◆ h_L1TauEtavsTauEta_

◆ h_L1TauETvsTauET_EB_

◆ h_L1TauETvsTauET_EB_EE_

◆ h_L1TauETvsTauET_EE_

◆ h_L1TauPhivsTauPhi_EB_

◆ h_L1TauPhivsTauPhi_EB_EE_

◆ h_L1TauPhivsTauPhi_EE_

◆ h_nVertex_

◆ h_resolutionTauET_EB_

◆ h_resolutionTauET_EB_EE_

◆ h_resolutionTauET_EE_

◆ h_resolutionTauEta_

◆ h_resolutionTauPhi_EB_

◆ h_resolutionTauPhi_EB_EE_

◆ h_resolutionTauPhi_EE_

◆ h_tagAndProbeMass_

◆ histDefinitions_

◆ histFolder_

◆ m_AntiEleExists

◆ m_BField

◆ m_hltConfig

◆ m_L1tL1tTaus

◆ m_L1tPtCuts

◆ m_L1tTaus

◆ m_MaxHltTauDR

◆ m_MaxL1tTauDR

◆ m_MaxTauEta

◆ m_ProbeTaus