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	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () 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->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 (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, edm::Provenance::stable(), 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_FULL_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 718 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().