#include <BaseTnPEfficiencyTask.h>

Inheritance diagram for BaseTnPEfficiencyTask:

Public Member Functions
	BaseTnPEfficiencyTask (const edm::ParameterSet &config)
	Constructor. More...

	~BaseTnPEfficiencyTask () override
	Destructor. More...

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

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

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

void	beginStream (edm::StreamID id) final

	DQMEDAnalyzer ()

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

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

virtual bool	getCanSaveByLumi ()

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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

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

Protected Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &context) override
	Analyze. More...

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

void	dqmBeginRun (const edm::Run &run, const edm::EventSetup &context) override
	BeginRun. More...

bool	hasTrigger (std::vector< int > &trigIndices, const trigger::TriggerObjectCollection &trigObjs, edm::Handle< trigger::TriggerEvent > &trigEvent, const reco::Muon &muon)

virtual std::string	topFolder () const =0
	Return the top folder. More...

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Protected Attributes
const double	m_borderCut

const bool	m_detailedAnalysis

const double	m_dxCut

std::map< std::string, MonitorElement * >	m_histos

const edm::EDGetTokenT< reco::MuonCollection >	m_muToken

int	m_nEvents

std::vector< std::vector< unsigned > >	m_probeIndices

std::vector< std::vector< unsigned > >	m_tagIndices

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

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Private Attributes
const double	m_dxyCut

const double	m_dzCut

const double	m_highPairMassCut

HLTConfigProvider	m_hltConfig

const double	m_lowPairMassCut

const edm::EDGetTokenT< std::vector< reco::Vertex > >	m_primaryVerticesToken

const StringCutObjectSelector< reco::Candidate, true >	m_probeSelector

const StringCutObjectSelector< reco::Muon, true >	m_tagSelector

const edm::EDGetTokenT< trigger::TriggerEvent >	m_triggerEventToken

const edm::EDGetTokenT< edm::TriggerResults >	m_triggerResultsToken

std::vector< int >	m_trigIndices

const std::string	m_trigName

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

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

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

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

Detailed Description

Definition at line 35 of file BaseTnPEfficiencyTask.h.

Constructor & Destructor Documentation

◆ BaseTnPEfficiencyTask()

BaseTnPEfficiencyTask::BaseTnPEfficiencyTask ( const edm::ParameterSet & config )

Constructor.

Definition at line 25 of file BaseTnPEfficiencyTask.cc.

References LogTrace.

     : m_nEvents(0),
       m_muToken(consumes<reco::MuonCollection>(config.getUntrackedParameter<edm::InputTag>("inputTagMuons"))),
       m_borderCut(config.getUntrackedParameter<double>("borderCut")),
       m_dxCut(config.getUntrackedParameter<double>("dx_cut")),
       m_detailedAnalysis(config.getUntrackedParameter<bool>("detailedAnalysis")),
       m_primaryVerticesToken(
           consumes<std::vector<reco::Vertex>>(config.getUntrackedParameter<edm::InputTag>("inputTagPrimaryVertices"))),
       m_triggerResultsToken(
           consumes<edm::TriggerResults>(config.getUntrackedParameter<edm::InputTag>("trigResultsTag"))),
       m_triggerEventToken(consumes<trigger::TriggerEvent>(config.getUntrackedParameter<edm::InputTag>("trigEventTag"))),
       m_trigName(config.getUntrackedParameter<std::string>("trigName")),
       m_probeSelector(config.getUntrackedParameter<std::string>("probeCut")),
       m_dxyCut(config.getUntrackedParameter<double>("probeDxyCut")),
       m_dzCut(config.getUntrackedParameter<double>("probeDzCut")),
       m_tagSelector(config.getUntrackedParameter<std::string>("tagCut")),
       m_lowPairMassCut(config.getUntrackedParameter<double>("lowPairMassCut")),
       m_highPairMassCut(config.getUntrackedParameter<double>("highPairMassCut")) {
   LogTrace("DQMOffline|MuonDPG|BaseTnPEfficiencyTask") << "[BaseTnPEfficiencyTask]: Constructor" << std::endl;
 }

◆ ~BaseTnPEfficiencyTask()

BaseTnPEfficiencyTask::~BaseTnPEfficiencyTask ( )

override

Destructor.

Definition at line 46 of file BaseTnPEfficiencyTask.cc.

References LogTrace, and m_nEvents.

                                               {
   LogTrace("DQMOffline|MuonDPG|BaseTnPEfficiencyTask")
       << "[BaseTnPEfficiencyTask]: analyzed " << m_nEvents << " events" << std::endl;
 }

Member Function Documentation

◆ analyze()

void BaseTnPEfficiencyTask::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	context
	)

overrideprotectedvirtual

Analyze.

Reimplemented from DQMEDAnalyzer.

Reimplemented in RPCTnPEfficiencyTask, DTTnPEfficiencyTask, and CSCTnPEfficiencyTask.

Definition at line 89 of file BaseTnPEfficiencyTask.cc.

References funct::abs(), reco::LeafCandidate::charge(), hasTrigger(), edm::HandleBase::isValid(), m_dxyCut, m_dzCut, m_highPairMassCut, m_histos, m_muToken, m_nEvents, m_primaryVerticesToken, m_probeIndices, m_probeSelector, m_tagIndices, m_tagSelector, m_triggerEventToken, m_triggerResultsToken, m_trigIndices, HLT_2022v12_cff::muon, PDWG_BPHSkim_cff::muons, reco::LeafCandidate::polarP4(), egammaIdentification::pt_max, makeGlobalPositionRcd_cfg::tag, PDWG_DiPhoton_SD_cff::triggerEvent, triggerResults, and bphysicsOniaDQM_cfi::vertex.

Referenced by CSCTnPEfficiencyTask::analyze(), DTTnPEfficiencyTask::analyze(), and RPCTnPEfficiencyTask::analyze().

                                                                                        {
   ++m_nEvents;
 
   edm::Handle<reco::MuonCollection> muons;
   if (!event.getByToken(m_muToken, muons))
     return;
 
   edm::Handle<std::vector<reco::Vertex>> vtxs;
   if (!event.getByToken(m_primaryVerticesToken, vtxs))
     return;
   const reco::Vertex& vertex = vtxs->at(0);
 
   edm::Handle<edm::TriggerResults> triggerResults;
   if (!event.getByToken(m_triggerResultsToken, triggerResults))
     return;
 
   edm::Handle<trigger::TriggerEvent> triggerEvent;
   if (!event.getByToken(m_triggerEventToken, triggerEvent))
     return;
 
   //common tnp variables
   std::vector<unsigned> preSel_tag_indices;
   std::vector<unsigned> tag_indices;
   std::vector<unsigned> preSel_probe_indices;
   std::vector<unsigned> probe_indices;
 
   if (muons.isValid() && vtxs.isValid()) {
     //Is there a better way to initialize two different type variables?
     for (auto [muon, muonColl_index] = std::tuple{std::vector<reco::Muon>::const_iterator{(*muons).begin()}, 0};
          muon != (*muons).end();
          ++muon, ++muonColl_index) {
       bool trigMatch = false;
 
       //Getting trigger infos for tag selection
       if (triggerResults.isValid() && triggerEvent.isValid()) {
         const trigger::TriggerObjectCollection trigObjColl = triggerEvent->getObjects();
         trigMatch = hasTrigger(m_trigIndices, trigObjColl, triggerEvent, *muon);
       }
 
       //Probe selection
       if (m_probeSelector(*muon) && (std::abs(muon->muonBestTrack()->dxy(vertex.position())) < m_dxyCut) &&
           (std::abs(muon->muonBestTrack()->dz(vertex.position())) < m_dzCut)) {
         preSel_probe_indices.push_back(muonColl_index);
       }
       //Tag selection
       if (m_tagSelector(*muon) && trigMatch) {
         preSel_tag_indices.push_back(muonColl_index);
       }
 
     }  //loop over muons
   }
 
   //Probe selection
   for (const auto i_tag : preSel_tag_indices) {
     reco::Muon tag = (*muons).at(i_tag);
     float pt_max = 0.;
     int max_pt_idx;
     bool pair_found = false;
 
     for (const auto i_probe : preSel_probe_indices) {
       //Prevent tag and probe to be the same object
       if (i_probe == i_tag)
         continue;
 
       reco::Muon preSel_probe = (*muons).at(i_probe);
 
       int pair_charge_product = tag.charge() * preSel_probe.charge();
 
       //check if tag+probe pair is compatible with Z decay
       if (pair_charge_product > 0)
         continue;
 
       float pair_mass = (tag.polarP4() + preSel_probe.polarP4()).M();
       m_histos.find("pairMass")->second->Fill(pair_mass);
 
       if (pair_mass < m_lowPairMassCut || pair_mass > m_highPairMassCut)
         continue;
 
       float pair_pt = (tag.polarP4() + preSel_probe.polarP4()).Pt();
       if (pair_pt > pt_max) {
         pair_found = true;
         pt_max = pair_pt;
         max_pt_idx = i_probe;
       }
     }
     if (pair_found) {
       probe_indices.push_back(max_pt_idx);
       tag_indices.push_back(i_tag);
     }
   }
 
   m_probeIndices.push_back(probe_indices);
   m_tagIndices.push_back(tag_indices);
 }

◆ bookHistograms()

void BaseTnPEfficiencyTask::bookHistograms	(	DQMStore::IBooker &	iBooker,
		edm::Run const &	run,
		edm::EventSetup const &	context
	)

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Reimplemented in RPCTnPEfficiencyTask, CSCTnPEfficiencyTask, and DTTnPEfficiencyTask.

Definition at line 68 of file BaseTnPEfficiencyTask.cc.

References ALCARECODTCalibSynchDQM_cff::baseDir, dqm::implementation::IBooker::book1D(), LogTrace, m_detailedAnalysis, m_histos, Pi, dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, and topFolder().

Referenced by CSCTnPEfficiencyTask::bookHistograms(), DTTnPEfficiencyTask::bookHistograms(), and RPCTnPEfficiencyTask::bookHistograms().

                                                                          {
   LogTrace("DQMOffline|MuonDPG|BaseTnPEfficiencyTask") << "[BaseTnPEfficiencyTask]: bookHistograms" << std::endl;
 
   if (m_detailedAnalysis) {
     std::string baseDir = topFolder() + "/detailed/";
     iBooker.setCurrentFolder(baseDir);
 
     LogTrace("DQMOffline|MuonDPG|BaseTnPEfficiencyTask")
         << "[BaseTnPEfficiencyTask]: booking histos in " << baseDir << std::endl;
 
     m_histos["probePt"] = iBooker.book1D("probePt", "probePt;probe p_{T} [GeV];Events", 125, 0., 250.);
     m_histos["probeEta"] = iBooker.book1D("probeEta", "probeEta;probe #eta;Events", 24, -2.4, 2.4);
     m_histos["probePhi"] = iBooker.book1D("probePhi", "probePhi;probe #phi; Events", 36, -TMath::Pi(), TMath::Pi());
     m_histos["probeNumberOfMatchedStations"] = iBooker.book1D(
         "probeNumberOfMatchedStations", "probeNumberOfMatchedStations;Number of matched stations;Events", 5, 0., 5);
     m_histos["pairMass"] = iBooker.book1D("pairMass", "pairMass", 25, 70., 120.);
   }
 }

◆ dqmBeginRun()

void BaseTnPEfficiencyTask::dqmBeginRun	(	const edm::Run &	run,
		const edm::EventSetup &	context
	)

overrideprotectedvirtual

BeginRun.

Reimplemented from DQMEDAnalyzer.

Definition at line 51 of file BaseTnPEfficiencyTask.cc.

References visDQMUpload::context, HLTConfigProvider::init(), m_hltConfig, m_trigIndices, m_trigName, hltMonBTagIPClient_cfi::pathName, writedatasetfile::run, HLTConfigProvider::size(), simpleEdmComparison::tName, and HLTConfigProvider::triggerName().

                                                                                        {
   bool changed = true;
   m_hltConfig.init(run, context, "HLT", changed);
 
   bool enableWildCard = true;
 
   TString tName = TString(m_trigName);
   TRegexp tNamePattern = TRegexp(tName, enableWildCard);
 
   for (unsigned iPath = 0; iPath < m_hltConfig.size(); ++iPath) {
     TString pathName = TString(m_hltConfig.triggerName(iPath));
     if (pathName.Contains(tNamePattern)) {
       m_trigIndices.push_back(static_cast<int>(iPath));
     }
   }
 }

◆ hasTrigger()

bool BaseTnPEfficiencyTask::hasTrigger	(	std::vector< int > &	trigIndices,
		const trigger::TriggerObjectCollection &	trigObjs,
		edm::Handle< trigger::TriggerEvent > &	trigEvent,
		const reco::Muon &	muon
	)

protected

Definition at line 184 of file BaseTnPEfficiencyTask.cc.

References HLTMuonOfflineAnalyzer_cfi::deltaR2, spr::find(), relativeConstraints::keys, m_hltConfig, HLTConfigProvider::moduleLabels(), B2GTnPMonitor_cfi::trigEvent, and TriggerAnalyzer::trigObjs.

Referenced by analyze().

                                                              {
   float dR2match = 999.;
   for (int trigIdx : trigIndices) {
     const std::vector<std::string> trigModuleLabels = m_hltConfig.moduleLabels(trigIdx);
 
     const unsigned trigModuleIndex =
         std::find(trigModuleLabels.begin(), trigModuleLabels.end(), "hltBoolEnd") - trigModuleLabels.begin() - 1;
     const unsigned hltFilterIndex = trigEvent->filterIndex(edm::InputTag(trigModuleLabels[trigModuleIndex], "", "HLT"));
     if (hltFilterIndex < trigEvent->sizeFilters()) {
       const trigger::Keys keys = trigEvent->filterKeys(hltFilterIndex);
       const trigger::Vids vids = trigEvent->filterIds(hltFilterIndex);
       const unsigned nTriggers = vids.size();
 
       for (unsigned iTrig = 0; iTrig < nTriggers; ++iTrig) {
         trigger::TriggerObject trigObj = trigObjs[keys[iTrig]];
         float dR2 = deltaR2(muon, trigObj);
         if (dR2 < dR2match)
           dR2match = dR2;
       }
     }
   }
 
   return dR2match < 0.01;
 }

◆ topFolder()

virtual std::string BaseTnPEfficiencyTask::topFolder ( ) const

protectedpure virtual

Return the top folder.

Implemented in RPCTnPEfficiencyTask, CSCTnPEfficiencyTask, and DTTnPEfficiencyTask.

Referenced by bookHistograms().

Member Data Documentation

◆ m_borderCut

const double BaseTnPEfficiencyTask::m_borderCut

protected

Definition at line 68 of file BaseTnPEfficiencyTask.h.

Referenced by CSCTnPEfficiencyTask::analyze(), DTTnPEfficiencyTask::analyze(), and RPCTnPEfficiencyTask::analyze().

◆ m_detailedAnalysis

const bool BaseTnPEfficiencyTask::m_detailedAnalysis

protected

Definition at line 70 of file BaseTnPEfficiencyTask.h.

Referenced by CSCTnPEfficiencyTask::analyze(), DTTnPEfficiencyTask::analyze(), RPCTnPEfficiencyTask::analyze(), and bookHistograms().

◆ m_dxCut

const double BaseTnPEfficiencyTask::m_dxCut

protected

Definition at line 69 of file BaseTnPEfficiencyTask.h.

Referenced by CSCTnPEfficiencyTask::analyze(), DTTnPEfficiencyTask::analyze(), and RPCTnPEfficiencyTask::analyze().

◆ m_dxyCut

const double BaseTnPEfficiencyTask::m_dxyCut

private

Definition at line 82 of file BaseTnPEfficiencyTask.h.

Referenced by analyze().

◆ m_dzCut

const double BaseTnPEfficiencyTask::m_dzCut

private

Definition at line 83 of file BaseTnPEfficiencyTask.h.

Referenced by analyze().

◆ m_highPairMassCut

const double BaseTnPEfficiencyTask::m_highPairMassCut

private

Definition at line 92 of file BaseTnPEfficiencyTask.h.

Referenced by analyze().

◆ m_histos

std::map<std::string, MonitorElement*> BaseTnPEfficiencyTask::m_histos

protected

Definition at line 62 of file BaseTnPEfficiencyTask.h.

Referenced by CSCTnPEfficiencyTask::analyze(), DTTnPEfficiencyTask::analyze(), RPCTnPEfficiencyTask::analyze(), analyze(), RPCTnPEfficiencyTask::bookEndcapHistos(), CSCTnPEfficiencyTask::bookHistograms(), DTTnPEfficiencyTask::bookHistograms(), RPCTnPEfficiencyTask::bookHistograms(), bookHistograms(), DTTnPEfficiencyTask::bookWheelHistos(), and RPCTnPEfficiencyTask::bookWheelHistos().

◆ m_hltConfig

HLTConfigProvider BaseTnPEfficiencyTask::m_hltConfig

private

Definition at line 78 of file BaseTnPEfficiencyTask.h.

Referenced by dqmBeginRun(), and hasTrigger().

◆ m_lowPairMassCut

const double BaseTnPEfficiencyTask::m_lowPairMassCut

private

Definition at line 91 of file BaseTnPEfficiencyTask.h.

◆ m_muToken

const edm::EDGetTokenT<reco::MuonCollection> BaseTnPEfficiencyTask::m_muToken

protected

Definition at line 66 of file BaseTnPEfficiencyTask.h.

Referenced by CSCTnPEfficiencyTask::analyze(), DTTnPEfficiencyTask::analyze(), RPCTnPEfficiencyTask::analyze(), and analyze().

◆ m_nEvents

int BaseTnPEfficiencyTask::m_nEvents

protected

Definition at line 64 of file BaseTnPEfficiencyTask.h.

Referenced by analyze(), ~BaseTnPEfficiencyTask(), CSCTnPEfficiencyTask::~CSCTnPEfficiencyTask(), DTTnPEfficiencyTask::~DTTnPEfficiencyTask(), and RPCTnPEfficiencyTask::~RPCTnPEfficiencyTask().

◆ m_primaryVerticesToken

const edm::EDGetTokenT<std::vector<reco::Vertex> > BaseTnPEfficiencyTask::m_primaryVerticesToken

private

Definition at line 73 of file BaseTnPEfficiencyTask.h.

Referenced by analyze().

◆ m_probeIndices

std::vector<std::vector<unsigned> > BaseTnPEfficiencyTask::m_probeIndices

protected

Definition at line 59 of file BaseTnPEfficiencyTask.h.

Referenced by CSCTnPEfficiencyTask::analyze(), DTTnPEfficiencyTask::analyze(), RPCTnPEfficiencyTask::analyze(), and analyze().

◆ m_probeSelector

const StringCutObjectSelector<reco::Candidate, true> BaseTnPEfficiencyTask::m_probeSelector

private

Definition at line 81 of file BaseTnPEfficiencyTask.h.

Referenced by analyze().

◆ m_tagIndices

std::vector<std::vector<unsigned> > BaseTnPEfficiencyTask::m_tagIndices

protected

Definition at line 60 of file BaseTnPEfficiencyTask.h.

Referenced by analyze().

◆ m_tagSelector

const StringCutObjectSelector<reco::Muon, true> BaseTnPEfficiencyTask::m_tagSelector

private

Definition at line 86 of file BaseTnPEfficiencyTask.h.

Referenced by analyze().

◆ m_triggerEventToken

const edm::EDGetTokenT<trigger::TriggerEvent> BaseTnPEfficiencyTask::m_triggerEventToken

private

Definition at line 75 of file BaseTnPEfficiencyTask.h.

Referenced by analyze().

◆ m_triggerResultsToken

const edm::EDGetTokenT<edm::TriggerResults> BaseTnPEfficiencyTask::m_triggerResultsToken

private

Definition at line 74 of file BaseTnPEfficiencyTask.h.

Referenced by analyze().

◆ m_trigIndices

std::vector<int> BaseTnPEfficiencyTask::m_trigIndices

private

Definition at line 89 of file BaseTnPEfficiencyTask.h.

Referenced by analyze(), and dqmBeginRun().

◆ m_trigName

const std::string BaseTnPEfficiencyTask::m_trigName

private

Definition at line 77 of file BaseTnPEfficiencyTask.h.

Referenced by dqmBeginRun().

Public Member Functions

Protected Member Functions

Protected Attributes

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ BaseTnPEfficiencyTask()

◆ ~BaseTnPEfficiencyTask()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ dqmBeginRun()

◆ hasTrigger()

◆ topFolder()

Member Data Documentation

◆ m_borderCut

◆ m_detailedAnalysis

◆ m_dxCut

◆ m_dxyCut

◆ m_dzCut

◆ m_highPairMassCut

◆ m_histos

◆ m_hltConfig

◆ m_lowPairMassCut

◆ m_muToken

◆ m_nEvents

◆ m_primaryVerticesToken

◆ m_probeIndices

◆ m_probeSelector

◆ m_tagIndices

◆ m_tagSelector

◆ m_triggerEventToken

◆ m_triggerResultsToken

◆ m_trigIndices

◆ m_trigName