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::one::EDProducerBase
typedef EDProducerBase	ModuleType
Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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
virtual void	dqmBeginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
virtual void	dqmEndLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final
void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &, edm::EventSetup const &) final
void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, edm::Accumulator >
	EDProducer ()=default
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDProducerBase () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
	~ProducerBase () noexcept(false) override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

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
void	dqmEndRun (edm::Run const &run, edm::EventSetup const &eSetup) override
void	endJob () 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 edm::ProducerBase
ProducesCollector	producesCollector ()
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

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::PFTauDiscriminator >	AntiEleInputTag_
edm::EDGetTokenT< reco::PFTauDiscriminator >	AntiMuInputTag_
edm::EDGetTokenT< reco::BeamSpot >	BsInputTag_
edm::EDGetTokenT< reco::PFTauDiscriminator >	comb3TInputTag_
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_
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 edm::one::EDProducerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 77 of file L1TTauOffline.h.

Member Enumeration Documentation

enum L1TTauOffline::PlotConfig

Enumerator
nVertex
ETvsET
PHIvsPHI

Definition at line 82 of file L1TTauOffline.h.

{ nVertex, ETvsET, PHIvsPHI };

Constructor & Destructor Documentation

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::PFTauDiscriminator>(ps.getUntrackedParameter<edm::InputTag>("antiMuInputTag"))),
      AntiEleInputTag_(consumes<reco::PFTauDiscriminator>(ps.getUntrackedParameter<edm::InputTag>("antiEleInputTag"))),
      DecayModeFindingInputTag_(
          consumes<reco::PFTauDiscriminator>(ps.getUntrackedParameter<edm::InputTag>("decayModeFindingInputTag"))),
      comb3TInputTag_(consumes<reco::PFTauDiscriminator>(ps.getUntrackedParameter<edm::InputTag>("comb3TInputTag"))),
      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 ( )

override

Definition at line 113 of file L1TTauOffline.cc.

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

Member Function Documentation

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 164 of file L1TTauOffline.cc.

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

                                                                          {
  m_MaxTauEta = 2.1;
  m_MaxL1tTauDR = 0.5;
  m_MaxHltTauDR = 0.5;

  edm::Handle<reco::PFTauCollection> taus;
  e.getByToken(theTauCollection_, taus);

  if (!taus.isValid()) {
    edm::LogWarning("L1TTauOffline") << "invalid collection: reco::PFTauCollection " << std::endl;
    return;
  }

  edm::Handle<reco::MuonCollection> muons;
  e.getByToken(MuonInputTag_, muons);

  if (!muons.isValid()) {
    edm::LogWarning("L1TTauOffline") << "invalid collection: reco::MuonCollection " << std::endl;
    return;
  }

  edm::Handle<reco::BeamSpot> beamSpot;
  e.getByToken(BsInputTag_, beamSpot);

  if (!beamSpot.isValid()) {
    edm::LogWarning("L1TTauOffline") << "invalid collection: reco::BeamSpot " << std::endl;
    return;
  }

  edm::Handle<reco::VertexCollection> vertex;
  e.getByToken(VtxInputTag_, vertex);

  if (!vertex.isValid()) {
    edm::LogWarning("L1TTauOffline") << "invalid collection: reco::VertexCollection " << std::endl;
    return;
  }

  edm::Handle<l1t::TauBxCollection> l1tCands;
  e.getByToken(stage2CaloLayer2TauToken_, l1tCands);

  if (!l1tCands.isValid()) {
    edm::LogWarning("L1TTauOffline") << "invalid collection: l1t::TauBxCollection " << std::endl;
    return;
  }

  edm::Handle<edm::TriggerResults> trigResults;
  e.getByToken(triggerResults_, trigResults);

  if (!trigResults.isValid()) {
    edm::LogWarning("L1TTauOffline") << "invalid collection: edm::TriggerResults " << std::endl;
    return;
  }

  edm::Handle<trigger::TriggerEvent> trigEvent;
  e.getByToken(triggerEvent_, trigEvent);

  if (!trigEvent.isValid()) {
    edm::LogWarning("L1TTauOffline") << "invalid collection: trigger::TriggerEvent " << std::endl;
    return;
  }

  edm::Handle<reco::PFMETCollection> mets;
  e.getByToken(MetInputTag_, mets);

  if (!mets.isValid()) {
    edm::LogWarning("L1TTauOffline") << "invalid collection: reco::PFMETCollection " << std::endl;
    return;
  }

  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
}

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 132 of file L1TTauOffline.cc.

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

                                                                                                {
  edm::LogInfo("L1TTauOffline") << "L1TTauOffline::bookHistograms" << std::endl;

  // book at beginRun
  bookTauHistos(ibooker);

  for (auto trigNamesIt = triggerPath_.begin(); trigNamesIt != triggerPath_.end(); trigNamesIt++) {
    std::string tNameTmp = (*trigNamesIt);
    std::string tNamePattern = "";
    std::size_t found0 = tNameTmp.find("*");
    if (found0 != std::string::npos)
      tNamePattern = tNameTmp.substr(0, tNameTmp.size() - 1);
    else
      tNamePattern = tNameTmp;

    int tIndex = -1;

    for (unsigned ipath = 0; ipath < m_hltConfig.size(); ++ipath) {
      std::string tmpName = m_hltConfig.triggerName(ipath);

      std::size_t found = tmpName.find(tNamePattern);
      if (found != std::string::npos) {
        tIndex = int(ipath);
        m_trigIndices.push_back(tIndex);
      }
    }
  }
}

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

private

Definition at line 330 of file L1TTauOffline.cc.

References dqmoffline::l1t::HistDefinition::binsX, dqmoffline::l1t::HistDefinition::binsY, dqm::dqmstoreimpl::DQMStore::IBooker::book1D(), dqm::dqmstoreimpl::DQMStore::IBooker::book2D(), dqm::dqmstoreimpl::DQMStore::IBooker::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::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, tauEfficiencyBins_, tauEfficiencyThresholds_, MessageLogger_cff::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().

                                                          {
  ibooker.cd();
  ibooker.setCurrentFolder(histFolder_);
  dqmoffline::l1t::HistDefinition nVertexDef = histDefinitions_[PlotConfig::nVertex];
  h_nVertex_ = ibooker.book1D(nVertexDef.name, nVertexDef.title, nVertexDef.nbinsX, nVertexDef.xmin, nVertexDef.xmax);
  h_tagAndProbeMass_ = ibooker.book1D("tagAndProbeMass", "Invariant mass of tag & probe pair", 100, 40, 140);

  dqmoffline::l1t::HistDefinition templateETvsET = histDefinitions_[PlotConfig::ETvsET];
  h_L1TauETvsTauET_EB_ = ibooker.book2D("L1TauETvsTauET_EB",
                                        "L1 Tau E_{T} vs PFTau E_{T} (EB); PFTau E_{T} (GeV); L1 Tau E_{T} (GeV)",
                                        templateETvsET.nbinsX,
                                        &templateETvsET.binsX[0],
                                        templateETvsET.nbinsY,
                                        &templateETvsET.binsY[0]);
  h_L1TauETvsTauET_EE_ = ibooker.book2D("L1TauETvsTauET_EE",
                                        "L1 Tau E_{T} vs PFTau E_{T} (EE); PFTau E_{T} (GeV); L1 Tau E_{T} (GeV)",
                                        templateETvsET.nbinsX,
                                        &templateETvsET.binsX[0],
                                        templateETvsET.nbinsY,
                                        &templateETvsET.binsY[0]);
  h_L1TauETvsTauET_EB_EE_ = ibooker.book2D("L1TauETvsTauET_EB_EE",
                                           "L1 Tau E_{T} vs PFTau E_{T} (EB+EE); PFTau E_{T} (GeV); L1 Tau E_{T} (GeV)",
                                           templateETvsET.nbinsX,
                                           &templateETvsET.binsX[0],
                                           templateETvsET.nbinsY,
                                           &templateETvsET.binsY[0]);

  dqmoffline::l1t::HistDefinition templatePHIvsPHI = histDefinitions_[PlotConfig::PHIvsPHI];
  h_L1TauPhivsTauPhi_EB_ =
      ibooker.book2D("L1TauPhivsTauPhi_EB",
                     "#phi_{tau}^{L1} vs #phi_{tau}^{offline} (EB); #phi_{tau}^{offline}; #phi_{tau}^{L1}",
                     templatePHIvsPHI.nbinsX,
                     templatePHIvsPHI.xmin,
                     templatePHIvsPHI.xmax,
                     templatePHIvsPHI.nbinsY,
                     templatePHIvsPHI.ymin,
                     templatePHIvsPHI.ymax);
  h_L1TauPhivsTauPhi_EE_ =
      ibooker.book2D("L1TauPhivsTauPhi_EE",
                     "#phi_{tau}^{L1} vs #phi_{tau}^{offline} (EE); #phi_{tau}^{offline}; #phi_{tau}^{L1}",
                     templatePHIvsPHI.nbinsX,
                     templatePHIvsPHI.xmin,
                     templatePHIvsPHI.xmax,
                     templatePHIvsPHI.nbinsY,
                     templatePHIvsPHI.ymin,
                     templatePHIvsPHI.ymax);
  h_L1TauPhivsTauPhi_EB_EE_ =
      ibooker.book2D("L1TauPhivsTauPhi_EB_EE",
                     "#phi_{tau}^{L1} vs #phi_{tau}^{offline} (EB+EE); #phi_{tau}^{offline}; #phi_{tau}^{L1}",
                     templatePHIvsPHI.nbinsX,
                     templatePHIvsPHI.xmin,
                     templatePHIvsPHI.xmax,
                     templatePHIvsPHI.nbinsY,
                     templatePHIvsPHI.ymin,
                     templatePHIvsPHI.ymax);

  h_L1TauEtavsTauEta_ =
      ibooker.book2D("L1TauEtavsTauEta", "L1 Tau #eta vs PFTau #eta; PFTau #eta; L1 Tau #eta", 100, -3, 3, 100, -3, 3);

  // tau resolutions
  h_resolutionTauET_EB_ = ibooker.book1D(
      "resolutionTauET_EB", "tau ET resolution (EB); (L1 Tau E_{T} - PFTau E_{T})/PFTau E_{T}; events", 50, -1, 1.5);
  h_resolutionTauET_EE_ = ibooker.book1D(
      "resolutionTauET_EE", "tau ET resolution (EE); (L1 Tau E_{T} - PFTau E_{T})/PFTau E_{T}; events", 50, -1, 1.5);
  h_resolutionTauET_EB_EE_ =
      ibooker.book1D("resolutionTauET_EB_EE",
                     "tau ET resolution (EB+EE); (L1 Tau E_{T} - PFTau E_{T})/PFTau E_{T}; events",
                     50,
                     -1,
                     1.5);

  h_resolutionTauPhi_EB_ = ibooker.book1D("resolutionTauPhi_EB",
                                          "#phi_{tau} resolution (EB); #phi_{tau}^{L1} - #phi_{tau}^{offline}; events",
                                          120,
                                          -0.3,
                                          0.3);
  h_resolutionTauPhi_EE_ = ibooker.book1D(
      "resolutionTauPhi_EE", "tau #phi resolution (EE); #phi_{tau}^{L1} - #phi_{tau}^{offline}; events", 120, -0.3, 0.3);
  h_resolutionTauPhi_EB_EE_ =
      ibooker.book1D("resolutionTauPhi_EB_EE",
                     "tau #phi resolution (EB+EE); #phi_{tau}^{L1} - #phi_{tau}^{offline}; events",
                     120,
                     -0.3,
                     0.3);

  h_resolutionTauEta_ =
      ibooker.book1D("resolutionTauEta", "tau #eta resolution  (EB); L1 Tau #eta - PFTau #eta; events", 120, -0.3, 0.3);

  // tau turn-ons
  ibooker.setCurrentFolder(efficiencyFolder_);
  std::vector<float> tauBins(tauEfficiencyBins_.begin(), tauEfficiencyBins_.end());
  int nBins = tauBins.size() - 1;
  float* tauBinArray = &(tauBins[0]);

  for (auto threshold : tauEfficiencyThresholds_) {
    std::string str_threshold = std::to_string(int(threshold));
    h_efficiencyIsoTauET_EB_pass_[threshold] =
        ibooker.book1D("efficiencyIsoTauET_EB_threshold_" + str_threshold + "_Num",
                       "iso tau efficiency (EB); PFTau E_{T} (GeV); events",
                       nBins,
                       tauBinArray);
    h_efficiencyIsoTauET_EE_pass_[threshold] =
        ibooker.book1D("efficiencyIsoTauET_EE_threshold_" + str_threshold + "_Num",
                       "iso tau efficiency (EE); PFTau E_{T} (GeV); events",
                       nBins,
                       tauBinArray);
    h_efficiencyIsoTauET_EB_EE_pass_[threshold] =
        ibooker.book1D("efficiencyIsoTauET_EB_EE_threshold_" + str_threshold + "_Num",
                       "iso tau efficiency (EB+EE); PFTau E_{T} (GeV); events",
                       nBins,
                       tauBinArray);

    h_efficiencyIsoTauET_EB_total_[threshold] =
        ibooker.book1D("efficiencyIsoTauET_EB_threshold_" + str_threshold + "_Den",
                       "iso tau efficiency (EB); PFTau E_{T} (GeV); events",
                       nBins,
                       tauBinArray);
    h_efficiencyIsoTauET_EE_total_[threshold] =
        ibooker.book1D("efficiencyIsoTauET_EE_threshold_" + str_threshold + "_Den",
                       "iso tau efficiency (EE); PFTau E_{T} (GeV); events",
                       nBins,
                       tauBinArray);
    h_efficiencyIsoTauET_EB_EE_total_[threshold] =
        ibooker.book1D("efficiencyIsoTauET_EB_EE_threshold_" + str_threshold + "_Den",
                       "iso tau efficiency (EB+EE); PFTau E_{T} (GeV); events",
                       nBins,
                       tauBinArray);

    // non iso
    h_efficiencyNonIsoTauET_EB_pass_[threshold] =
        ibooker.book1D("efficiencyNonIsoTauET_EB_threshold_" + str_threshold + "_Num",
                       "inclusive tau efficiency (EB); PFTau E_{T} (GeV); events",
                       nBins,
                       tauBinArray);
    h_efficiencyNonIsoTauET_EE_pass_[threshold] =
        ibooker.book1D("efficiencyNonIsoTauET_EE_threshold_" + str_threshold + "_Num",
                       "inclusive tau efficiency (EE); PFTau E_{T} (GeV); events",
                       nBins,
                       tauBinArray);
    h_efficiencyNonIsoTauET_EB_EE_pass_[threshold] =
        ibooker.book1D("efficiencyNonIsoTauET_EB_EE_threshold_" + str_threshold + "_Num",
                       "inclusive tau efficiency (EB+EE); PFTau E_{T} (GeV); events",
                       nBins,
                       tauBinArray);

    h_efficiencyNonIsoTauET_EB_total_[threshold] =
        ibooker.book1D("efficiencyNonIsoTauET_EB_threshold_" + str_threshold + "_Den",
                       "inclusive tau efficiency (EB); PFTau E_{T} (GeV); events",
                       nBins,
                       tauBinArray);
    h_efficiencyNonIsoTauET_EE_total_[threshold] =
        ibooker.book1D("efficiencyNonIsoTauET_EE_threshold_" + str_threshold + "_Den",
                       "inclusive tau efficiency (EE); PFTau E_{T} (GeV); events",
                       nBins,
                       tauBinArray);
    h_efficiencyNonIsoTauET_EB_EE_total_[threshold] =
        ibooker.book1D("efficiencyNonIsoTauET_EB_EE_threshold_" + str_threshold + "_Den",
                       "inclusive tau efficiency (EB+EE); PFTau E_{T} (GeV); events",
                       nBins,
                       tauBinArray);
  }

  ibooker.cd();

  return;
}

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

private

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

private

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

private

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

overrideprotected

Definition at line 120 of file L1TTauOffline.cc.

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

{
  bool changed = true;
  m_hltConfig.init(run, iSetup, trigProcess_, changed);

  edm::LogInfo("L1TTauOffline") << "L1TTauOffline::beginRun" << std::endl;
}

void L1TTauOffline::dqmEndRun ( edm::Run const &  run, edm::EventSetup const &  eSetup  )

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 323 of file L1TTauOffline.cc.

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

void L1TTauOffline::endJob ( void  )

overrideprotectedvirtual

Reimplemented from edm::one::EDProducerBase.

Definition at line 677 of file L1TTauOffline.cc.

References normalise2DHistogramsToBinArea().

{ normalise2DHistogramsToBinArea(); }

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

protected

Definition at line 497 of file L1TTauOffline.cc.

References reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthY(), edm::HandleBase::isValid(), reco::BeamSpot::position(), reco::Vertex::position(), BeamMonitor_cff::primaryVertex, and reco::BeamSpot::sigmaZ().

Referenced by analyze().

                                                                                            {
  reco::Vertex::Point posVtx;
  reco::Vertex::Error errVtx;

  bool hasPrimaryVertex = false;

  if (vertex.isValid()) {
    for (auto vertexIt = vertex->begin(); vertexIt != vertex->end(); ++vertexIt) {
      if (vertexIt->isValid() && !vertexIt->isFake()) {
        posVtx = vertexIt->position();
        errVtx = vertexIt->error();
        hasPrimaryVertex = true;
        break;
      }
    }
  }

  if (!hasPrimaryVertex) {
    posVtx = beamSpot->position();
    errVtx(0, 0) = beamSpot->BeamWidthX();
    errVtx(1, 1) = beamSpot->BeamWidthY();
    errVtx(2, 2) = beamSpot->sigmaZ();
  }

  const reco::Vertex primaryVertex(posVtx, errVtx);

  return primaryVertex;
}

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

protected

Definition at line 624 of file L1TTauOffline.cc.

References AntiEleInputTag_, AntiMuInputTag_, comb3TInputTag_, DecayModeFindingInputTag_, HCALHighEnergyHPDFilter_cfi::energy, PVValHelper::eta, edm::Event::getByToken(), edm::HandleBase::isValid(), m_ProbeTaus, m_TightMuons, phi, and DiDispStaMuonMonitor_cfi::pt.

Referenced by analyze().

                                                           {
  m_ProbeTaus.clear();

  edm::Handle<reco::PFTauDiscriminator> 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::PFTauDiscriminator> antiele;
  iEvent.getByToken(AntiEleInputTag_, antiele);
  if (!antiele.isValid()) {
    edm::LogWarning("L1TTauOffline") << "invalid collection: reco::PFTauDiscriminator " << std::endl;
    return;
  }

  edm::Handle<reco::PFTauDiscriminator> 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;
    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] > 0.5 &&
          (*antiele)[tauCandidate] > 0.5 && (*dmf)[tauCandidate] > 0.5 && (*comb3T)[tauCandidate] > 0.5) {
        if (mymu.DeltaR(mytau) > 0.5 && (mymu + mytau).M() > 40 && (mymu + mytau).M() < 80 &&
            m_TightMuons[0]->charge() * tauIt->charge() < 0) {
          m_ProbeTaus.push_back(&(*tauIt));
        }
      }
    }
  }
}

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

protected

Definition at line 556 of file L1TTauOffline.cc.

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

Referenced by analyze().

                                                                                {
  m_TauL1tPairs.clear();

  vector<l1t::Tau> l1tContainer;
  l1tContainer.reserve(l1tCands->size() + 1);

  for (auto tau = l1tCands->begin(0); tau != l1tCands->end(0); ++tau) {
    l1tContainer.push_back(*tau);
  }

  for (auto probeTauIt = m_ProbeTaus.begin(); probeTauIt != m_ProbeTaus.end(); ++probeTauIt) {
    TauL1TPair pairBestCand((*probeTauIt), nullptr);

    for (auto l1tIt = l1tContainer.begin(); l1tIt != l1tContainer.end(); ++l1tIt) {
      TauL1TPair pairTmpCand((*probeTauIt), &(*l1tIt));

      if (pairTmpCand.dR() < m_MaxL1tTauDR && pairTmpCand.l1tPt() > pairBestCand.l1tPt())
        pairBestCand = pairTmpCand;
    }

    m_TauL1tPairs.push_back(pairBestCand);
  }
}

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

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

Referenced by analyze().

                                                                                   {
  m_TightMuons.clear();

  const reco::PFMET* pfmet = nullptr;
  pfmet = &(mets->front());

  int nb_mu = 0;

  for (auto muonIt2 = muons->begin(); muonIt2 != muons->end(); ++muonIt2) {
    if (fabs(muonIt2->eta()) < 2.4 && muonIt2->pt() > 10 && muon::isLooseMuon((*muonIt2)) &&
        (muonIt2->pfIsolationR04().sumChargedHadronPt +
         max(muonIt2->pfIsolationR04().sumNeutralHadronEt + muonIt2->pfIsolationR04().sumPhotonEt -
                 0.5 * muonIt2->pfIsolationR04().sumPUPt,
             0.0)) /
                muonIt2->pt() <
            0.3) {
      ++nb_mu;
    }
  }
  bool foundTightMu = false;
  for (auto muonIt = muons->begin(); muonIt != muons->end(); ++muonIt) {
    if (!matchHlt(trigEvent, &(*muonIt)))
      continue;
    float muiso = (muonIt->pfIsolationR04().sumChargedHadronPt +
                   max(muonIt->pfIsolationR04().sumNeutralHadronEt + muonIt->pfIsolationR04().sumPhotonEt -
                           0.5 * muonIt->pfIsolationR04().sumPUPt,
                       0.0)) /
                  muonIt->pt();

    if (muiso < 0.1 && nb_mu < 2 && !foundTightMu && fabs(muonIt->eta()) < 2.1 && muonIt->pt() > 24 &&
        muon::isLooseMuon((*muonIt))) {
      float mt = sqrt(pow(muonIt->pt() + pfmet->pt(), 2) - pow(muonIt->px() + pfmet->px(), 2) -
                      pow(muonIt->py() + pfmet->py(), 2));
      if (mt < 30) {
        m_TightMuons.push_back(&(*muonIt));
        foundTightMu = true;
      }
    }
  }
}

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

protected

Definition at line 527 of file L1TTauOffline.cc.

References PbPb_ZMuSkimMuonDPG_cff::deltaR, trigger::TriggerEvent::filterIds(), trigger::TriggerEvent::filterIndex(), trigger::TriggerEvent::filterKeys(), trigger::TriggerEvent::getObjects(), HLT_2018_cff::InputTag, m_hltConfig, m_MaxHltTauDR, m_trigIndices, HLTConfigProvider::moduleLabels(), HLTConfigProvider::size(), TriggerAnalyzer::trigObjs, and trigProcess_.

Referenced by getTightMuons().

                                                                                                     {
  double matchDeltaR = 9999;

  trigger::TriggerObjectCollection trigObjs = triggerEvent->getObjects();

  for (auto trigIndexIt = m_trigIndices.begin(); trigIndexIt != m_trigIndices.end(); ++trigIndexIt) {
    const vector<string> moduleLabels(m_hltConfig.moduleLabels(*trigIndexIt));
    const unsigned moduleIndex = m_hltConfig.size((*trigIndexIt)) - 2;

    const unsigned hltFilterIndex = triggerEvent->filterIndex(InputTag(moduleLabels[moduleIndex], "", trigProcess_));

    if (hltFilterIndex < triggerEvent->sizeFilters()) {
      const Keys triggerKeys(triggerEvent->filterKeys(hltFilterIndex));
      const Vids triggerVids(triggerEvent->filterIds(hltFilterIndex));

      const unsigned nTriggers = triggerVids.size();
      for (size_t iTrig = 0; iTrig < nTriggers; ++iTrig) {
        const TriggerObject trigObject = trigObjs[triggerKeys[iTrig]];

        double dRtmp = deltaR((*muon), trigObject);
        if (dRtmp < matchDeltaR)
          matchDeltaR = dRtmp;
      }
    }
  }

  return (matchDeltaR < m_MaxHltTauDR);
}

void L1TTauOffline::normalise2DHistogramsToBinArea ( )

private

Definition at line 679 of file L1TTauOffline.cc.

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

Referenced by endJob().

                                                   {
  std::vector<MonitorElement*> monElementstoNormalize = {h_L1TauETvsTauET_EB_,
                                                         h_L1TauETvsTauET_EE_,
                                                         h_L1TauETvsTauET_EB_EE_,
                                                         h_L1TauPhivsTauPhi_EB_,
                                                         h_L1TauPhivsTauPhi_EE_,
                                                         h_L1TauPhivsTauPhi_EB_EE_,
                                                         h_L1TauEtavsTauEta_};

  for (auto mon : monElementstoNormalize) {
    if (mon != nullptr) {
      auto h = mon->getTH2F();
      if (h != nullptr) {
        h->Scale(1, "width");
      }
    }
  }
}

Member Data Documentation

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

private

Definition at line 130 of file L1TTauOffline.h.

Referenced by getProbeTaus().

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

private

Definition at line 129 of file L1TTauOffline.h.

Referenced by getProbeTaus().

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

private

Definition at line 136 of file L1TTauOffline.h.

Referenced by analyze().

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

private

Definition at line 132 of file L1TTauOffline.h.

Referenced by getProbeTaus().

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

private

Definition at line 131 of file L1TTauOffline.h.

Referenced by getProbeTaus().

std::string L1TTauOffline::efficiencyFolder_

private

Definition at line 142 of file L1TTauOffline.h.

Referenced by bookTauHistos().

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

private

Definition at line 195 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

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

private

Definition at line 205 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

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

private

Definition at line 193 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

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

private

Definition at line 203 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

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

private

Definition at line 194 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

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

private

Definition at line 204 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

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

private

Definition at line 199 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

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

private

Definition at line 209 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

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

private

Definition at line 197 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

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

private

Definition at line 207 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

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

private

Definition at line 198 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

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

private

Definition at line 208 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

MonitorElement* L1TTauOffline::h_L1TauEtavsTauEta_

private

Definition at line 179 of file L1TTauOffline.h.

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

MonitorElement* L1TTauOffline::h_L1TauETvsTauET_EB_

private

Definition at line 171 of file L1TTauOffline.h.

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

MonitorElement* L1TTauOffline::h_L1TauETvsTauET_EB_EE_

private

Definition at line 173 of file L1TTauOffline.h.

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

MonitorElement* L1TTauOffline::h_L1TauETvsTauET_EE_

private

Definition at line 172 of file L1TTauOffline.h.

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

MonitorElement* L1TTauOffline::h_L1TauPhivsTauPhi_EB_

private

Definition at line 175 of file L1TTauOffline.h.

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

MonitorElement* L1TTauOffline::h_L1TauPhivsTauPhi_EB_EE_

private

Definition at line 177 of file L1TTauOffline.h.

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

MonitorElement* L1TTauOffline::h_L1TauPhivsTauPhi_EE_

private

Definition at line 176 of file L1TTauOffline.h.

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

MonitorElement* L1TTauOffline::h_nVertex_

private

Definition at line 167 of file L1TTauOffline.h.

Referenced by bookTauHistos().

MonitorElement* L1TTauOffline::h_resolutionTauET_EB_

private

Definition at line 182 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

MonitorElement* L1TTauOffline::h_resolutionTauET_EB_EE_

private

Definition at line 184 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

MonitorElement* L1TTauOffline::h_resolutionTauET_EE_

private

Definition at line 183 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

MonitorElement* L1TTauOffline::h_resolutionTauEta_

private

Definition at line 190 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

MonitorElement* L1TTauOffline::h_resolutionTauPhi_EB_

private

Definition at line 186 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

MonitorElement* L1TTauOffline::h_resolutionTauPhi_EB_EE_

private

Definition at line 188 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

MonitorElement* L1TTauOffline::h_resolutionTauPhi_EE_

private

Definition at line 187 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

MonitorElement* L1TTauOffline::h_tagAndProbeMass_

private

Definition at line 168 of file L1TTauOffline.h.

Referenced by bookTauHistos().

dqmoffline::l1t::HistDefinitions L1TTauOffline::histDefinitions_

private

Definition at line 146 of file L1TTauOffline.h.

Referenced by bookTauHistos().

std::string L1TTauOffline::histFolder_

private

Definition at line 141 of file L1TTauOffline.h.

Referenced by bookTauHistos().

edm::ESHandle<MagneticField> L1TTauOffline::m_BField

private

Definition at line 123 of file L1TTauOffline.h.

Referenced by analyze().

HLTConfigProvider L1TTauOffline::m_hltConfig

private

Definition at line 121 of file L1TTauOffline.h.

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

BXVector<l1t::Tau> L1TTauOffline::m_L1tL1tTaus

private

Definition at line 155 of file L1TTauOffline.h.

std::vector<int> L1TTauOffline::m_L1tPtCuts

private

Definition at line 158 of file L1TTauOffline.h.

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

private

Definition at line 153 of file L1TTauOffline.h.

float L1TTauOffline::m_MaxHltTauDR

private

Definition at line 162 of file L1TTauOffline.h.

Referenced by analyze(), and matchHlt().

float L1TTauOffline::m_MaxL1tTauDR

private

Definition at line 161 of file L1TTauOffline.h.

Referenced by analyze(), and getTauL1tPairs().

float L1TTauOffline::m_MaxTauEta

private

Definition at line 160 of file L1TTauOffline.h.

Referenced by analyze().

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

private

Definition at line 149 of file L1TTauOffline.h.

Referenced by getProbeTaus(), and getTauL1tPairs().

edm::ESHandle<Propagator> L1TTauOffline::m_propagatorAlong

private

Definition at line 124 of file L1TTauOffline.h.

edm::ESHandle<Propagator> L1TTauOffline::m_propagatorOpposite

private

Definition at line 125 of file L1TTauOffline.h.

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

private

Definition at line 154 of file L1TTauOffline.h.

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

private

Definition at line 152 of file L1TTauOffline.h.

std::vector<TauL1TPair> L1TTauOffline::m_TauL1tPairs

private

Definition at line 150 of file L1TTauOffline.h.

Referenced by analyze(), and getTauL1tPairs().

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

private

Definition at line 148 of file L1TTauOffline.h.

Referenced by getProbeTaus(), and getTightMuons().

std::vector<int> L1TTauOffline::m_trigIndices

private

Definition at line 164 of file L1TTauOffline.h.

Referenced by bookHistograms(), and matchHlt().

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

private

Definition at line 134 of file L1TTauOffline.h.

Referenced by analyze().

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

private

Definition at line 133 of file L1TTauOffline.h.

Referenced by analyze().

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

static

Initial value:

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

Definition at line 84 of file L1TTauOffline.h.

Referenced by TauL1TPair::dR().

math::XYZPoint L1TTauOffline::PVPoint_

private

Definition at line 119 of file L1TTauOffline.h.

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

private

Definition at line 143 of file L1TTauOffline.h.

Referenced by analyze().

std::vector<double> L1TTauOffline::tauEfficiencyBins_

private

Definition at line 145 of file L1TTauOffline.h.

Referenced by bookTauHistos().

std::vector<int> L1TTauOffline::tauEfficiencyThresholds_

private

Definition at line 144 of file L1TTauOffline.h.

Referenced by analyze(), and bookTauHistos().

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

private

Definition at line 128 of file L1TTauOffline.h.

Referenced by analyze().

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

private

Definition at line 137 of file L1TTauOffline.h.

Referenced by analyze().

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

private

Definition at line 140 of file L1TTauOffline.h.

Referenced by bookHistograms().

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

private

Definition at line 139 of file L1TTauOffline.h.

Referenced by analyze().

std::string L1TTauOffline::trigProcess_

private

Definition at line 138 of file L1TTauOffline.h.

Referenced by dqmBeginRun(), and matchHlt().

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

private

Definition at line 135 of file L1TTauOffline.h.

Referenced by analyze().