Inheritance diagram for CSCTnPEfficiencyTask:

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

	~CSCTnPEfficiencyTask () override
	Destructor. More...

Public Member Functions inherited from BaseTnPEfficiencyTask
	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

std::string	topFolder () const override
	Return the top folder. More...

Protected Member Functions inherited from BaseTnPEfficiencyTask
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)

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

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

Protected Attributes inherited from BaseTnPEfficiencyTask
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_

Detailed Description

Definition at line 14 of file CSCTnPEfficiencyTask.cc.

Constructor & Destructor Documentation

◆ CSCTnPEfficiencyTask()

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

Constructor.

Definition at line 31 of file CSCTnPEfficiencyTask.cc.

References LogTrace.

                                                                         : BaseTnPEfficiencyTask(config) {
   LogTrace("DQMOffline|MuonDPG|CSCTnPEfficiencyTask") << "[CSCTnPEfficiencyTask]: Constructor" << std::endl;
 }

◆ ~CSCTnPEfficiencyTask()

CSCTnPEfficiencyTask::~CSCTnPEfficiencyTask ( )

override

Destructor.

Definition at line 35 of file CSCTnPEfficiencyTask.cc.

References LogTrace, and BaseTnPEfficiencyTask::m_nEvents.

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

Member Function Documentation

◆ analyze()

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

overrideprotectedvirtual

Analyze.

Reimplemented from BaseTnPEfficiencyTask.

Definition at line 115 of file CSCTnPEfficiencyTask.cc.

References funct::abs(), BaseTnPEfficiencyTask::analyze(), visDQMUpload::context, MuonSubdetId::CSC, PVValHelper::dx, mps_fire::i, dqmiolumiharvest::j, BaseTnPEfficiencyTask::m_borderCut, BaseTnPEfficiencyTask::m_detailedAnalysis, BaseTnPEfficiencyTask::m_dxCut, BaseTnPEfficiencyTask::m_histos, BaseTnPEfficiencyTask::m_muToken, BaseTnPEfficiencyTask::m_probeIndices, DiMuonV_cfg::muons, relativeConstraints::ring, relativeConstraints::station, and CSCDetId::zendcap().

                                                                                       {
   BaseTnPEfficiencyTask::analyze(event, context);
 
   edm::Handle<reco::MuonCollection> muons;
   event.getByToken(m_muToken, muons);
 
   //CSC variables
   std::vector<std::vector<int>> probe_coll_CSC_zend;
   std::vector<std::vector<int>> probe_coll_CSC_ring;
   std::vector<std::vector<int>> probe_coll_CSC_sta;
   std::vector<std::vector<float>> probe_coll_CSC_dx;
   std::vector<uint8_t> probe_coll_CSC_staMatch;
 
   std::vector<unsigned> probe_indices;
   if (!m_probeIndices.empty())
     probe_indices = m_probeIndices.back();
 
   //Fill probe dx + subdetector coordinates
   for (const auto i : probe_indices) {
     //CSC variables
     std::vector<int> probe_CSC_zend;
     std::vector<int> probe_CSC_ring;
     std::vector<int> probe_CSC_sta;
     std::vector<float> probe_CSC_dx;
     uint8_t CSC_stationMatching = 0;
 
     float csc_matched = false;  // fill detailed plots only for probes matching CSC
 
     for (const auto& chambMatch : (*muons).at(i).matches()) {
       // look in CSCs
       if (chambMatch.detector() == MuonSubdetId::CSC) {
         if (chambMatch.edgeX < m_borderCut && chambMatch.edgeY < m_borderCut) {
           csc_matched = true;  //fill detailed plots if at least one CSC match
 
           CSCDetId chId(chambMatch.id.rawId());
 
           int zendcap = chId.zendcap();
           int ring = chId.ring();
           int station = chId.station();
 
           reco::MuonSegmentMatch closest_matchedSegment;
           double smallestDx = 99999.;
           for (auto& seg : chambMatch.segmentMatches) {
             float dx = std::abs(chambMatch.x - seg.x);
             if (dx < smallestDx) {
               smallestDx = dx;
               closest_matchedSegment = seg;
             }
           }
 
           CSC_stationMatching = CSC_stationMatching | (1 << (station - 1));
 
           if (station == 1 && ring == 4 && chambMatch.y < -31.5) {
             probe_CSC_zend.push_back(zendcap);
             probe_CSC_ring.push_back(ring);
             probe_CSC_sta.push_back(station);
             probe_CSC_dx.push_back(smallestDx);
           } else if (station == 1 && ring == 1 && chambMatch.y > -31.5) {
             probe_CSC_zend.push_back(zendcap);
             probe_CSC_ring.push_back(ring);
             probe_CSC_sta.push_back(station);
             probe_CSC_dx.push_back(smallestDx);
           } else if (station > 1 || ring == 2 || ring == 3) {
             probe_CSC_zend.push_back(zendcap);
             probe_CSC_ring.push_back(ring);
             probe_CSC_sta.push_back(station);
             probe_CSC_dx.push_back(smallestDx);
           }
         }
       } else
         continue;
     }  //loop over chamber matches
 
     //Fill detailed plots
     if (m_detailedAnalysis && csc_matched) {
       m_histos.find("probeEta")->second->Fill((*muons).at(i).eta());
       m_histos.find("probePhi")->second->Fill((*muons).at(i).phi());
       m_histos.find("probeNumberOfMatchedStations")->second->Fill((*muons).at(i).numberOfMatchedStations());
       m_histos.find("probePt")->second->Fill((*muons).at(i).pt());
     }
 
     //Fill CSC variables
     probe_coll_CSC_zend.push_back(probe_CSC_zend);
     probe_coll_CSC_ring.push_back(probe_CSC_ring);
     probe_coll_CSC_sta.push_back(probe_CSC_sta);
     probe_coll_CSC_dx.push_back(probe_CSC_dx);
     probe_coll_CSC_staMatch.push_back(CSC_stationMatching);
   }  //loop over probe collection
 
   //Loop over probes
   for (unsigned i = 0; i < probe_indices.size(); ++i) {
     uint8_t CSC_matchPatt = probe_coll_CSC_staMatch.at(i);
 
     //Loop over CSC matches
     unsigned nCSC_matches = probe_coll_CSC_zend.at(i).size();
     for (unsigned j = 0; j < nCSC_matches; ++j) {
       int CSC_zendcap = probe_coll_CSC_zend.at(i).at(j);
       int CSC_sta = probe_coll_CSC_sta.at(i).at(j);
       int CSC_ring = probe_coll_CSC_ring.at(i).at(j);
       float CSC_dx = probe_coll_CSC_dx.at(i).at(j);
 
       //Fill CSC plots
       if ((CSC_matchPatt & (1 << (CSC_sta - 1))) != 0 &&  //avoids 0 station matching
           (CSC_matchPatt & (1 << (CSC_sta - 1))) !=
               CSC_matchPatt)  //avoids matching with the station under consideration only
       {
         if (CSC_dx < m_dxCut) {
           m_histos.find("CSC_nPassingProbe_allCh")->second->Fill(CSC_zendcap * CSC_sta, CSC_ring);
           m_histos.find("CSC_nPassingProbe_allCh_1D")->second->Fill(CSC_zendcap * CSC_sta);
         } else {
           m_histos.find("CSC_nFailingProbe_allCh")->second->Fill(CSC_zendcap * CSC_sta, CSC_ring);
           m_histos.find("CSC_nFailingProbe_allCh_1D")->second->Fill(CSC_zendcap * CSC_sta);
         }
       }
     }
   }
 }

◆ bookHistograms()

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

overrideprotectedvirtual

Reimplemented from BaseTnPEfficiencyTask.

Definition at line 40 of file CSCTnPEfficiencyTask.cc.

References ALCARECODTCalibSynchDQM_cff::baseDir, dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), BaseTnPEfficiencyTask::bookHistograms(), visDQMUpload::context, mps_fire::i, LogTrace, BaseTnPEfficiencyTask::m_histos, writedatasetfile::run, dqm::impl::MonitorElement::setAxisTitle(), dqm::impl::MonitorElement::setBinLabel(), dqm::implementation::NavigatorBase::setCurrentFolder(), to_string(), and topFolder().

                                                                         {
   BaseTnPEfficiencyTask::bookHistograms(iBooker, run, context);
 
   LogTrace("DQMOffline|MuonDPG|CSCTnPEfficiencyTask") << "[CSCTnPEfficiencyTask]: bookHistograms" << std::endl;
 
   auto baseDir = topFolder() + "Task/";
   iBooker.setCurrentFolder(baseDir);
 
   MonitorElement* me_CSC_pass_allCh =
       iBooker.book2D("CSC_nPassingProbe_allCh", "CSC_nPassingProbe_allCh", 9, -4., 5., 4, 0., 4.5);
   MonitorElement* me_CSC_fail_allCh =
       iBooker.book2D("CSC_nFailingProbe_allCh", "CSC_nFailingProbe_allCh", 9, -4., 5., 4, 0., 4.5);
 
   MonitorElement* me_CSC_pass_allCh_1D =
       iBooker.book1D("CSC_nPassingProbe_allCh_1D", "CSC_nPassingProbe_allCh_1D", 9, -4., 5.);
   MonitorElement* me_CSC_fail_allCh_1D =
       iBooker.book1D("CSC_nFailingProbe_allCh_1D", "CSC_nFailingProbe_allCh_1D", 9, -4., 5.);
 
   me_CSC_pass_allCh->setBinLabel(1, "ME-4", 1);
   me_CSC_pass_allCh->setBinLabel(2, "ME-3", 1);
   me_CSC_pass_allCh->setBinLabel(3, "ME-2", 1);
   me_CSC_pass_allCh->setBinLabel(4, "ME-1", 1);
   me_CSC_pass_allCh->setBinLabel(6, "ME1", 1);
   me_CSC_pass_allCh->setBinLabel(7, "ME2", 1);
   me_CSC_pass_allCh->setBinLabel(8, "ME3", 1);
   me_CSC_pass_allCh->setBinLabel(9, "ME4", 1);
   for (int i = 1; i < 5; ++i) {
     me_CSC_pass_allCh->setBinLabel(i, std::to_string(i), 2);
   }
   me_CSC_pass_allCh->setAxisTitle("Ring", 2);
   me_CSC_pass_allCh->setAxisTitle("Number of passing probes", 3);
 
   me_CSC_fail_allCh->setBinLabel(1, "ME-4", 1);
   me_CSC_fail_allCh->setBinLabel(2, "ME-3", 1);
   me_CSC_fail_allCh->setBinLabel(3, "ME-2", 1);
   me_CSC_fail_allCh->setBinLabel(4, "ME-1", 1);
   me_CSC_fail_allCh->setBinLabel(6, "ME1", 1);
   me_CSC_fail_allCh->setBinLabel(7, "ME2", 1);
   me_CSC_fail_allCh->setBinLabel(8, "ME3", 1);
   me_CSC_fail_allCh->setBinLabel(9, "ME4", 1);
   for (int i = 1; i < 5; ++i) {
     me_CSC_fail_allCh->setBinLabel(i, std::to_string(i), 2);
   }
   me_CSC_fail_allCh->setAxisTitle("Ring", 2);
   me_CSC_fail_allCh->setAxisTitle("Number of failing probes", 3);
 
   me_CSC_pass_allCh_1D->setBinLabel(1, "ME-4", 1);
   me_CSC_pass_allCh_1D->setBinLabel(2, "ME-3", 1);
   me_CSC_pass_allCh_1D->setBinLabel(3, "ME-2", 1);
   me_CSC_pass_allCh_1D->setBinLabel(4, "ME-1", 1);
   me_CSC_pass_allCh_1D->setBinLabel(6, "ME1", 1);
   me_CSC_pass_allCh_1D->setBinLabel(7, "ME2", 1);
   me_CSC_pass_allCh_1D->setBinLabel(8, "ME3", 1);
   me_CSC_pass_allCh_1D->setBinLabel(9, "ME4", 1);
   me_CSC_pass_allCh_1D->setAxisTitle("Number of passing probes", 2);
 
   me_CSC_fail_allCh_1D->setBinLabel(1, "ME-4", 1);
   me_CSC_fail_allCh_1D->setBinLabel(2, "ME-3", 1);
   me_CSC_fail_allCh_1D->setBinLabel(3, "ME-2", 1);
   me_CSC_fail_allCh_1D->setBinLabel(4, "ME-1", 1);
   me_CSC_fail_allCh_1D->setBinLabel(6, "ME1", 1);
   me_CSC_fail_allCh_1D->setBinLabel(7, "ME2", 1);
   me_CSC_fail_allCh_1D->setBinLabel(8, "ME3", 1);
   me_CSC_fail_allCh_1D->setBinLabel(9, "ME4", 1);
   me_CSC_fail_allCh_1D->setAxisTitle("Number of failing probes", 2);
 
   m_histos["CSC_nPassingProbe_allCh"] = me_CSC_pass_allCh;
   m_histos["CSC_nFailingProbe_allCh"] = me_CSC_fail_allCh;
 
   m_histos["CSC_nPassingProbe_allCh_1D"] = me_CSC_pass_allCh_1D;
   m_histos["CSC_nFailingProbe_allCh_1D"] = me_CSC_fail_allCh_1D;
 }