L1TTauOffline Class Reference

#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 >
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public 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
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

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::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 81 of file L1TTauOffline.h.

Member Enumeration Documentation

PlotConfig

enum L1TTauOffline::PlotConfig

Enumerator
nVertex
ETvsET
PHIvsPHI

Definition at line 86 of file L1TTauOffline.h.

{ nVertex, ETvsET, PHIvsPHI };

Constructor & Destructor Documentation

L1TTauOffline()

L1TTauOffline::L1TTauOffline

(

const edm::ParameterSet &

ps

)

Definition at line 45 of file L1TTauOffline.cc.

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

~L1TTauOffline()

L1TTauOffline::~L1TTauOffline ( )

override

Definition at line 119 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 170 of file L1TTauOffline.cc.

                                                                          {
  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;
  }
 
  eSetup.get<IdealMagneticFieldRecord>().get(m_BField);
  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
}

References pwdgSkimBPark_cfi::beamSpot, BXVector< T >::begin(), BsInputTag_, counter, MillePedeFileConverter_cfg::e, BXVector< T >::end(), PVValHelper::eta, dqm::impl::MonitorElement::Fill(), edm::EventSetup::get(), get, 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, 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_LED_IterMethod_cfg::threshold, B2GTnPMonitor_cfi::trigEvent, triggerEvent_, triggerResults_, bphysicsOniaDQM_cfi::vertex, and VtxInputTag_.

bookHistograms()

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 138 of file L1TTauOffline.cc.

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

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

bookTauHistos()

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

private

Definition at line 332 of file L1TTauOffline.cc.

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

References dqmoffline::l1t::HistDefinition::binsX, dqmoffline::l1t::HistDefinition::binsY, dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::NavigatorBase::cd(), efficiencyFolder_, L1TEGammaOffline_cfi::ETvsET, 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_, dqmoffline::l1t::HistDefinition::name, seedmultiplicitymonitor_newtracking_cfi::nBins, dqmoffline::l1t::HistDefinition::nbinsX, dqmoffline::l1t::HistDefinition::nbinsY, L1TEGammaOffline_cfi::nVertex, L1TEGammaOffline_cfi::PHIvsPHI, dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, tauEfficiencyBins_, tauEfficiencyThresholds_, remoteMonitoring_LED_IterMethod_cfg::threshold, dqmoffline::l1t::HistDefinition::title, dqmoffline::l1t::HistDefinition::xmax, dqmoffline::l1t::HistDefinition::xmin, dqmoffline::l1t::HistDefinition::ymax, and dqmoffline::l1t::HistDefinition::ymin.

Referenced by bookHistograms().

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

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

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

getPrimaryVertex()

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

protected

Definition at line 499 of file L1TTauOffline.cc.

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

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

Referenced by analyze().

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

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

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

Referenced by analyze().

getTauL1tPairs()

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

protected

Definition at line 558 of file L1TTauOffline.cc.

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

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

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

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

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

matchHlt()

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

protected

Definition at line 529 of file L1TTauOffline.cc.

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

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

Referenced by getTightMuons().

normalise2DHistogramsToBinArea()

void L1TTauOffline::normalise2DHistogramsToBinArea ( )

private

Definition at line 717 of file L1TTauOffline.cc.

                                                   {
  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");
      }
    }
  }
}

References 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.

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 132 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 204 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 214 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 202 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 212 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 203 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 213 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 208 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 218 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 206 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 216 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 207 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 217 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

h_L1TauEtavsTauEta_

MonitorElement* L1TTauOffline::h_L1TauEtavsTauEta_

private

Definition at line 188 of file L1TTauOffline.h.

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

h_L1TauETvsTauET_EB_

MonitorElement* L1TTauOffline::h_L1TauETvsTauET_EB_

private

Definition at line 180 of file L1TTauOffline.h.

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

h_L1TauETvsTauET_EB_EE_

MonitorElement* L1TTauOffline::h_L1TauETvsTauET_EB_EE_

private

Definition at line 182 of file L1TTauOffline.h.

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

h_L1TauETvsTauET_EE_

MonitorElement* L1TTauOffline::h_L1TauETvsTauET_EE_

private

Definition at line 181 of file L1TTauOffline.h.

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

h_L1TauPhivsTauPhi_EB_

MonitorElement* L1TTauOffline::h_L1TauPhivsTauPhi_EB_

private

Definition at line 184 of file L1TTauOffline.h.

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

h_L1TauPhivsTauPhi_EB_EE_

MonitorElement* L1TTauOffline::h_L1TauPhivsTauPhi_EB_EE_

private

Definition at line 186 of file L1TTauOffline.h.

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

h_L1TauPhivsTauPhi_EE_

MonitorElement* L1TTauOffline::h_L1TauPhivsTauPhi_EE_

private

Definition at line 185 of file L1TTauOffline.h.

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

h_nVertex_

MonitorElement* L1TTauOffline::h_nVertex_

private

Definition at line 176 of file L1TTauOffline.h.

Referenced by bookTauHistos().

h_resolutionTauET_EB_

MonitorElement* L1TTauOffline::h_resolutionTauET_EB_

private

Definition at line 191 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

h_resolutionTauET_EB_EE_

MonitorElement* L1TTauOffline::h_resolutionTauET_EB_EE_

private

Definition at line 193 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

h_resolutionTauET_EE_

MonitorElement* L1TTauOffline::h_resolutionTauET_EE_

private

Definition at line 192 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

h_resolutionTauEta_

MonitorElement* L1TTauOffline::h_resolutionTauEta_

private

Definition at line 199 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

h_resolutionTauPhi_EB_

MonitorElement* L1TTauOffline::h_resolutionTauPhi_EB_

private

Definition at line 195 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

h_resolutionTauPhi_EB_EE_

MonitorElement* L1TTauOffline::h_resolutionTauPhi_EB_EE_

private

Definition at line 197 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

h_resolutionTauPhi_EE_

MonitorElement* L1TTauOffline::h_resolutionTauPhi_EE_

private

Definition at line 196 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

h_tagAndProbeMass_

MonitorElement* L1TTauOffline::h_tagAndProbeMass_

private

Definition at line 177 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_BField

edm::ESHandle<MagneticField> L1TTauOffline::m_BField

private

Definition at line 124 of file L1TTauOffline.h.

Referenced by analyze().

m_hltConfig

HLTConfigProvider L1TTauOffline::m_hltConfig

private

Definition at line 122 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 125 of file L1TTauOffline.h.

m_propagatorOpposite

edm::ESHandle<Propagator> L1TTauOffline::m_propagatorOpposite

private

Definition at line 126 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().

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 128 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 88 of file L1TTauOffline.h.

PVPoint_

math::XYZPoint L1TTauOffline::PVPoint_

private

Definition at line 120 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 131 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().