#include <CSCStubEfficiencyValidation.h>

Inheritance diagram for CSCStubEfficiencyValidation:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	bookHistograms (DQMStore::IBooker &)

	CSCStubEfficiencyValidation (const edm::ParameterSet &pset, edm::ConsumesCollector &&iC)

	~CSCStubEfficiencyValidation () override

Public Member Functions inherited from CSCBaseValidation
	CSCBaseValidation (const edm::ParameterSet &ps)

void	setGeometry (const CSCGeometry *geom)

void	setSimHitMap (const PSimHitMap *simHitMap)

virtual	~CSCBaseValidation ()

Private Attributes
std::unique_ptr< CSCStubMatcher >	cscStubMatcher_

MonitorElement *	etaALCTDenom [10]

MonitorElement *	etaALCTNum [10]

MonitorElement *	etaCLCTDenom [10]

MonitorElement *	etaCLCTNum [10]

MonitorElement *	etaLCTDenom [10]

MonitorElement *	etaLCTNum [10]

std::vector< double >	etaMaxs_

std::vector< double >	etaMins_

edm::EDGetTokenT< edm::SimTrackContainer >	simTrackInput_

edm::EDGetTokenT< edm::SimVertexContainer >	simVertexInput_

Additional Inherited Members
Public Types inherited from CSCBaseValidation
typedef dqm::legacy::DQMStore	DQMStore

typedef dqm::legacy::MonitorElement	MonitorElement

Protected Member Functions inherited from CSCBaseValidation
const CSCLayer *	findLayer (int detId) const

bool	isSimTrackGood (const SimTrack &t) const

Protected Attributes inherited from CSCBaseValidation
bool	doSim_

double	simTrackMaxEta_

double	simTrackMinEta_

double	simTrackMinPt_

const CSCGeometry *	theCSCGeometry

const PSimHitMap *	theSimHitMap

Detailed Description

Definition at line 10 of file CSCStubEfficiencyValidation.h.

Constructor & Destructor Documentation

◆ CSCStubEfficiencyValidation()

CSCStubEfficiencyValidation::CSCStubEfficiencyValidation	(	const edm::ParameterSet &	pset,
		edm::ConsumesCollector &&	iC
	)

Definition at line 6 of file CSCStubEfficiencyValidation.cc.

References cscStubMatcher_, etaMaxs_, etaMins_, eostools::move(), muonDTDigis_cfi::pset, cscDigiValidation_cfi::simTrack, simTrackInput_, ecalDetailedTimeRecHit_cfi::simVertex, and simVertexInput_.

     : CSCBaseValidation(pset) {
   const auto& simVertex = pset.getParameter<edm::ParameterSet>("simVertex");
   simVertexInput_ = iC.consumes<edm::SimVertexContainer>(simVertex.getParameter<edm::InputTag>("inputTag"));
   const auto& simTrack = pset.getParameter<edm::ParameterSet>("simTrack");
   simTrackInput_ = iC.consumes<edm::SimTrackContainer>(simTrack.getParameter<edm::InputTag>("inputTag"));
 
   // Initialize stub matcher
   cscStubMatcher_ = std::make_unique<CSCStubMatcher>(pset, std::move(iC));
 
   // get the eta ranges
   etaMins_ = pset.getParameter<std::vector<double>>("etaMins");
   etaMaxs_ = pset.getParameter<std::vector<double>>("etaMaxs");
 }

◆ ~CSCStubEfficiencyValidation()

CSCStubEfficiencyValidation::~CSCStubEfficiencyValidation ( )

override

Definition at line 21 of file CSCStubEfficiencyValidation.cc.

21 {}

Member Function Documentation

◆ analyze()

void CSCStubEfficiencyValidation::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	eventSetup
	)

overridevirtual

Implements CSCBaseValidation.

Definition at line 55 of file CSCStubEfficiencyValidation.cc.

References nano_mu_digi_cff::chamberType, cscStubMatcher_, MillePedeFileConverter_cfg::e, etaALCTDenom, etaALCTNum, etaCLCTDenom, etaCLCTNum, etaLCTDenom, etaLCTNum, options_cfi::eventSetup, dqm::impl::MonitorElement::Fill(), mps_fire::i, CSCDetId::iChamberType(), CSCBaseValidation::isSimTrackGood(), edm::Handle< T >::product(), simTrackInput_, simVertexInput_, and submitPVValidationJobs::t.

                                                                                             {
   // Define handles
   edm::Handle<edm::SimTrackContainer> sim_tracks;
   edm::Handle<edm::SimVertexContainer> sim_vertices;
 
   // Use token to retreive event information
   e.getByToken(simTrackInput_, sim_tracks);
   e.getByToken(simVertexInput_, sim_vertices);
 
   // Initialize StubMatcher
   cscStubMatcher_->init(e, eventSetup);
 
   const edm::SimTrackContainer& sim_track = *sim_tracks.product();
   const edm::SimVertexContainer& sim_vert = *sim_vertices.product();
 
   // select simtracks for true muons
   edm::SimTrackContainer sim_track_selected;
   for (const auto& t : sim_track) {
     if (!isSimTrackGood(t))
       continue;
     sim_track_selected.push_back(t);
   }
 
   // Skip events with no selected simtracks
   if (sim_track_selected.empty())
     return;
 
   // Loop through good tracks, use corresponding vetrex to match stubs, then fill hists of chambers where the stub appears.
   for (const auto& t : sim_track_selected) {
     std::vector<bool> hitALCT(10);
     std::vector<bool> hitCLCT(10);
     std::vector<bool> hitLCT(10);
 
     // Match track to stubs with appropriate vertex
     cscStubMatcher_->match(t, sim_vert[t.vertIndex()]);
 
     // Store matched stubs.
     // Key: ChamberID, Value : CSCStubDigiContainer
     const auto& alcts = cscStubMatcher_->alcts();
     const auto& clcts = cscStubMatcher_->clcts();
     const auto& lcts = cscStubMatcher_->lcts();
 
     // denominator histograms
     for (int i = 0; i < 10; ++i) {
       etaALCTDenom[i]->Fill(t.momentum().eta());
       etaCLCTDenom[i]->Fill(t.momentum().eta());
       etaLCTDenom[i]->Fill(t.momentum().eta());
     }
 
     for (auto& [id, container] : alcts) {
       const CSCDetId cscId(id);
       const unsigned chamberType(cscId.iChamberType());
       hitALCT[chamberType - 1] = true;
     }
 
     for (auto& [id, container] : clcts) {
       const CSCDetId cscId(id);
       const unsigned chamberType(cscId.iChamberType());
       hitCLCT[chamberType - 1] = true;
     }
 
     for (auto& [id, container] : lcts) {
       const CSCDetId cscId(id);
       const unsigned chamberType(cscId.iChamberType());
       hitLCT[chamberType - 1] = true;
     }
 
     // numerator histograms
     for (int i = 0; i < 10; ++i) {
       if (hitALCT[i])
         etaALCTNum[i]->Fill(t.momentum().eta());
       if (hitCLCT[i])
         etaCLCTNum[i]->Fill(t.momentum().eta());
       if (hitLCT[i])
         etaLCTNum[i]->Fill(t.momentum().eta());
     }
   }
 }

◆ bookHistograms()

void CSCStubEfficiencyValidation::bookHistograms ( DQMStore::IBooker & iBooker )

Definition at line 23 of file CSCStubEfficiencyValidation.cc.

References dqm::implementation::IBooker::book1D(), CSCDetId::chamberName(), etaALCTDenom, etaALCTNum, etaCLCTDenom, etaCLCTNum, etaLCTDenom, etaLCTNum, etaMaxs_, etaMins_, dqm::legacy::MonitorElement::getTH1(), mps_fire::i, dqmiolumiharvest::j, dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, RandomServiceHelper::t1, RandomServiceHelper::t2, and RandomServiceHelper::t3.

                                                                          {
   for (int i = 1; i <= 10; ++i) {
     int j = i - 1;
     const std::string cn(CSCDetId::chamberName(i));
 
     std::string t1 = "ALCTEtaDenom_" + cn;
     std::string t2 = "CLCTEtaDenom_" + cn;
     std::string t3 = "LCTEtaDenom_" + cn;
 
     iBooker.setCurrentFolder("MuonCSCDigisV/CSCDigiTask/ALCT/Occupancy/");
     etaALCTDenom[j] = iBooker.book1D(t1, t1 + ";True Muon |#eta|; Entries", 50, etaMins_[j], etaMaxs_[j]);
     etaALCTDenom[j]->getTH1()->SetMinimum(0);
     t1 = "ALCTEtaNum_" + cn;
     etaALCTNum[j] = iBooker.book1D(t1, t1 + ";True Muon |#eta|; Entries", 50, etaMins_[j], etaMaxs_[j]);
     etaALCTNum[j]->getTH1()->SetMinimum(0);
 
     iBooker.setCurrentFolder("MuonCSCDigisV/CSCDigiTask/CLCT/Occupancy/");
     etaCLCTDenom[j] = iBooker.book1D(t2, t2 + ";True Muon |#eta|; Entries", 50, etaMins_[j], etaMaxs_[j]);
     etaCLCTDenom[j]->getTH1()->SetMinimum(0);
     t2 = "CLCTEtaNum_" + cn;
     etaCLCTNum[j] = iBooker.book1D(t2, t2 + ";True Muon |#eta|; Entries", 50, etaMins_[j], etaMaxs_[j]);
     etaCLCTNum[j]->getTH1()->SetMinimum(0);
 
     iBooker.setCurrentFolder("MuonCSCDigisV/CSCDigiTask/LCT/Occupancy/");
     etaLCTDenom[j] = iBooker.book1D(t3, t3 + ";True Muon |#eta|; Entries", 50, etaMins_[j], etaMaxs_[j]);
     etaLCTDenom[j]->getTH1()->SetMinimum(0);
     t3 = "LCTEtaNum_" + cn;
     etaLCTNum[j] = iBooker.book1D(t3, t3 + ";True Muon |#eta|; Entries", 50, etaMins_[j], etaMaxs_[j]);
     etaLCTNum[j]->getTH1()->SetMinimum(0);
   }
 }