#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_

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

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

override

Definition at line 120 of file L1TTauOffline.cc.

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

Member Function Documentation

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 171 of file L1TTauOffline.cc.

References beam_dqm_sourceclient-live_cfg::beamSpot, BsInputTag_, counter, PVValHelper::eta, dqm::impl::MonitorElement::Fill(), edm::Event::getByToken(), 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_, mFieldToken_, MuonInputTag_, patZpeak::muons, phi, beam_dqm_sourceclient-live_cfg::primaryVertex, DiDispStaMuonMonitor_cfi::pt, stage2CaloLayer2TauToken_, AlCaHLTBitMon_QueryRunRegistry::string, metsig::tau, tauEfficiencyThresholds_, runTauDisplay::taus, theTauCollection_, dtDQMClient_cfg::threshold, triggerEvent_, triggerResults_, 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
 }

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

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

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

private

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

private

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

private

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, and trigProcess_.

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

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 edm::HandleBase::isValid(), reco::Vertex::position(), and beam_dqm_sourceclient-live_cfg::primaryVertex.

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

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 628 of file L1TTauOffline.cc.

References AntiEleInputTag_, AntiEleWP_, AntiEleWPIndex_, AntiMuInputTag_, AntiMuWP_, AntiMuWPIndex_, comb3TInputTag_, comb3TWP_, comb3TWPIndex_, DecayModeFindingInputTag_, relval_parameters_module::energy, PVValHelper::eta, edm::Event::getByToken(), mps_fire::i, edm::HandleBase::isValid(), m_ProbeTaus, m_TightMuons, edm::parameterSet(), phi, phID_, edm::Event::processHistory(), edm::Event::processHistoryID(), edm::HandleBase::provenance(), DiDispStaMuonMonitor_cfi::pt, edm::Provenance::stable(), parallelization::uint(), and MuTauSkim_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 " << std::endl;
     return;
   }
 
   edm::Handle<reco::PFTauDiscriminator> dmf;
   iEvent.getByToken(DecayModeFindingInputTag_, dmf);
   if (!dmf.isValid()) {
     edm::LogWarning("L1TTauOffline") << "invalid collection: reco::PFTauDiscriminator " << std::endl;
     return;
   }
 
   edm::Handle<reco::TauDiscriminatorContainer> antiele;
   iEvent.getByToken(AntiEleInputTag_, antiele);
   if (!antiele.isValid()) {
     edm::LogWarning("L1TTauOffline") << "invalid collection: reco::PFTauDiscriminator " << std::endl;
     return;
   }
 
   edm::Handle<reco::TauDiscriminatorContainer> comb3T;
   iEvent.getByToken(comb3TInputTag_, comb3T);
   if (!comb3T.isValid()) {
     edm::LogWarning("L1TTauOffline") << "invalid collection: reco::PFTauDiscriminator " << std::endl;
     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;
         }
       }
       {
         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 (uint 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" << std::endl;
         continue;
       }
       if ((*antiele)[tauCandidate].workingPoints.empty()) {
         edm::LogWarning("L1TTauOffline") << "This offline tau has no antiele discriminator, skipping" << std::endl;
         continue;
       }
       if ((*comb3T)[tauCandidate].workingPoints.empty()) {
         edm::LogWarning("L1TTauOffline") << "This offline tau has no comb3T discriminator, skipping" << std::endl;
         continue;
       }
 
       if (fabs(tauIt->charge()) == 1 && fabs(tauIt->eta()) < 2.1 && tauIt->pt() > 20 &&
           (*antimu)[tauCandidate].workingPoints[AntiMuWPIndex_] &&
           (*antiele)[tauCandidate].workingPoints[AntiEleWPIndex_] && (*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));
         }
       }
     }
   }
 }

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

protected

Definition at line 560 of file L1TTauOffline.cc.

References TauL1TPair::dR(), TauL1TPair::l1tPt(), m_MaxL1tTauDR, m_ProbeTaus, m_TauL1tPairs, 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);
   }
 }

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 584 of file L1TTauOffline.cc.

References muon::isLooseMuon(), m_TightMuons, matchHlt(), SiStripPI::max, funct::pow(), reco::LeafCandidate::pt(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), and mathSSE::sqrt().

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

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

protected

Definition at line 531 of file L1TTauOffline.cc.

References HLT_FULL_cff::deltaR, HLT_FULL_cff::InputTag, m_hltConfig, m_MaxHltTauDR, m_trigIndices, HLTConfigProvider::moduleLabels(), HLTConfigProvider::size(), TriggerAnalyzer::trigObjs, and trigProcess_.

Referenced by getTightMuons().

                                                                                                      {
   double matchDeltaR = 9999;
 
   trigger::TriggerObjectCollection trigObjs = triggerEvent->getObjects();
 
   for (auto trigIndexIt = m_trigIndices.begin(); trigIndexIt != m_trigIndices.end(); ++trigIndexIt) {
     const vector<string> moduleLabels(m_hltConfig.moduleLabels(*trigIndexIt));
     const unsigned moduleIndex = m_hltConfig.size((*trigIndexIt)) - 2;
 
     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);
 }

void L1TTauOffline::normalise2DHistogramsToBinArea ( )

private

Definition at line 733 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");
       }
     }
   }
 }