#include <CSCALCTDigiValidation.h>

Inheritance diagram for CSCALCTDigiValidation:

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

void	bookHistograms (DQMStore::IBooker &)

	CSCALCTDigiValidation (const edm::ParameterSet &ps, edm::ConsumesCollector &&iC)

	~CSCALCTDigiValidation () 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::vector< double >	alctMaxBin_

std::vector< double >	alctMinBin_

std::vector< unsigned >	alctNBin_

edm::EDGetTokenT< CSCALCTDigiCollection >	alcts_Token_

std::vector< std::string >	alctVars_

std::map< uint32_t, std::map< std::string, MonitorElement * > >	chamberHistos

std::vector< std::string >	chambers_

edm::InputTag	inputTag_

MonitorElement *	theNDigisPerChamberPlots [10]

MonitorElement *	theNDigisPerEventPlot

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 CSCALCTDigiValidation.h.

Constructor & Destructor Documentation

◆ CSCALCTDigiValidation()

CSCALCTDigiValidation::CSCALCTDigiValidation	(	const edm::ParameterSet &	ps,
		edm::ConsumesCollector &&	iC
	)

Definition at line 9 of file CSCALCTDigiValidation.cc.

References alcts_Token_, edm::ParameterSet::getParameterSet(), inputTag_, and muonDTDigis_cfi::pset.

     : CSCBaseValidation(ps),
       theNDigisPerChamberPlots(),
       chambers_(ps.getParameter<std::vector<std::string>>("chambers")),
       // variables
       alctVars_(ps.getParameter<std::vector<std::string>>("alctVars")),
       // binning
       alctNBin_(ps.getParameter<std::vector<unsigned>>("alctNBin")),
       alctMinBin_(ps.getParameter<std::vector<double>>("alctMinBin")),
       alctMaxBin_(ps.getParameter<std::vector<double>>("alctMaxBin")) {
   const auto &pset = ps.getParameterSet("cscALCT");
   inputTag_ = pset.getParameter<edm::InputTag>("inputTag");
   alcts_Token_ = iC.consumes<CSCALCTDigiCollection>(inputTag_);
 }

◆ ~CSCALCTDigiValidation()

CSCALCTDigiValidation::~CSCALCTDigiValidation ( )

override

Definition at line 24 of file CSCALCTDigiValidation.cc.

24 {}

Member Function Documentation

◆ analyze()

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

overridevirtual

Implements CSCBaseValidation.

Definition at line 56 of file CSCALCTDigiValidation.cc.

References alcts_Token_, chamberHistos, chambers_, MillePedeFileConverter_cfg::e, edm::InputTag::encode(), CSCDetId::endcap(), dqm::impl::MonitorElement::Fill(), CSCDetId::iChamberType(), inputTag_, edm::HandleBase::isValid(), dqmiolumiharvest::j, FastTimerService_cff::range, theNDigisPerChamberPlots, and theNDigisPerEventPlot.

                                                                             {
   edm::Handle<CSCALCTDigiCollection> alcts;
   e.getByToken(alcts_Token_, alcts);
   if (!alcts.isValid()) {
     edm::LogError("CSCALCTDigiValidation") << "Cannot get ALCTs by label " << inputTag_.encode();
   }
   unsigned nDigisPerEvent = 0;
 
   for (auto j = alcts->begin(); j != alcts->end(); j++) {
     auto beginDigi = (*j).second.first;
     auto endDigi = (*j).second.second;
     const CSCDetId &detId((*j).first);
     int chamberType = detId.iChamberType();
     int nDigis = endDigi - beginDigi;
     nDigisPerEvent += nDigis;
     theNDigisPerChamberPlots[chamberType - 1]->Fill(nDigis);
 
     auto range = alcts->get((*j).first);
     // 1=forward (+Z); 2=backward (-Z)
     const unsigned typeCorrected(detId.endcap() == 1 ? chamberType - 2 : chamberType - 2 + chambers_.size());
     for (auto alct = range.first; alct != range.second; alct++) {
       if (alct->isValid()) {
         chamberHistos[typeCorrected]["alct_quality"]->Fill(alct->getQuality());
         chamberHistos[typeCorrected]["alct_wiregroup"]->Fill(alct->getKeyWG());
         chamberHistos[typeCorrected]["alct_bx"]->Fill(alct->getBX());
       }
     }
   }
   theNDigisPerEventPlot->Fill(nDigisPerEvent);
 }

◆ bookHistograms()

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

Definition at line 26 of file CSCALCTDigiValidation.cc.

References alctMaxBin_, alctMinBin_, alctNBin_, alctVars_, dqm::implementation::IBooker::book1D(), chamberHistos, CSCDetId::chamberName(), chambers_, B2GTnPMonitor_cfi::histTitle, mps_fire::i, crabWrapper::key, dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, RandomServiceHelper::t2, theNDigisPerChamberPlots, and theNDigisPerEventPlot.

                                                                    {
   iBooker.setCurrentFolder("MuonCSCDigisV/CSCDigiTask/ALCT/Occupancy/");
 
   theNDigisPerEventPlot = iBooker.book1D(
       "CSCALCTDigisPerEvent", "ALCT trigger primitives per event; Number of ALCTs; Entries", 100, 0, 100);
   for (int i = 1; i <= 10; ++i) {
     const std::string t2("CSCALCTDigisPerChamber_" + CSCDetId::chamberName(i));
     theNDigisPerChamberPlots[i - 1] = iBooker.book1D(
         t2, "Number of ALCTs per chamber " + CSCDetId::chamberName(i) + ";Number of ALCTs per chamber;Entries", 4, 0, 4);
   }
 
   // chamber type
   for (unsigned iType = 0; iType < chambers_.size(); iType++) {
     // consider CSC+ and CSC- separately
     for (unsigned iEndcap = 0; iEndcap < 2; iEndcap++) {
       const std::string eSign(iEndcap == 0 ? "+" : "-");
       // alct variable
       for (unsigned iVar = 0; iVar < alctVars_.size(); iVar++) {
         const std::string key("alct_" + alctVars_[iVar]);
         const std::string histName(key + "_" + chambers_[iType] + eSign);
         const std::string histTitle(chambers_[iType] + eSign + " ALCT " + alctVars_[iVar]);
         const unsigned iTypeCorrected(iEndcap == 0 ? iType : iType + chambers_.size());
         chamberHistos[iTypeCorrected][key] =
             iBooker.book1D(histName, histTitle, alctNBin_[iVar], alctMinBin_[iVar], alctMaxBin_[iVar]);
         chamberHistos[iTypeCorrected][key]->getTH1()->SetMinimum(0);
       }
     }
   }
 }