#include <CSCCorrelatedLCTDigiValidation.h>

Inheritance diagram for CSCCorrelatedLCTDigiValidation:

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

void	bookHistograms (DQMStore::IBooker &)

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

	~CSCCorrelatedLCTDigiValidation () 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::map< uint32_t, std::map< std::string, MonitorElement * > >	chamberHistos

std::vector< std::string >	chambers_

std::vector< unsigned >	chambersRun3_

edm::InputTag	inputTag_

bool	isRun3_

std::vector< double >	lctMaxBin_

std::vector< double >	lctMinBin_

std::vector< unsigned >	lctNBin_

edm::EDGetTokenT< CSCCorrelatedLCTDigiCollection >	lcts_Token_

std::vector< std::string >	lctVars_

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

Constructor & Destructor Documentation

◆ CSCCorrelatedLCTDigiValidation()

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

Definition at line 9 of file CSCCorrelatedLCTDigiValidation.cc.

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

     : CSCBaseValidation(ps),
       theNDigisPerChamberPlots(),
       chambers_(ps.getParameter<std::vector<std::string>>("chambers")),
       chambersRun3_(ps.getParameter<std::vector<unsigned>>("chambersRun3")),
       // variables
       lctVars_(ps.getParameter<std::vector<std::string>>("lctVars")),
       // binning
       lctNBin_(ps.getParameter<std::vector<unsigned>>("lctNBin")),
       lctMinBin_(ps.getParameter<std::vector<double>>("lctMinBin")),
       lctMaxBin_(ps.getParameter<std::vector<double>>("lctMaxBin")),
       isRun3_(ps.getParameter<bool>("isRun3")) {
   const auto &pset = ps.getParameterSet("cscMPLCT");
   inputTag_ = pset.getParameter<edm::InputTag>("inputTag");
   lcts_Token_ = iC.consumes<CSCCorrelatedLCTDigiCollection>(inputTag_);
 }

◆ ~CSCCorrelatedLCTDigiValidation()

CSCCorrelatedLCTDigiValidation::~CSCCorrelatedLCTDigiValidation ( )

override

Definition at line 26 of file CSCCorrelatedLCTDigiValidation.cc.

26 {}

Member Function Documentation

◆ analyze()

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

overridevirtual

Implements CSCBaseValidation.

Definition at line 88 of file CSCCorrelatedLCTDigiValidation.cc.

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

                                                                                      {
   edm::Handle<CSCCorrelatedLCTDigiCollection> lcts;
   e.getByToken(lcts_Token_, lcts);
   if (!lcts.isValid()) {
     edm::LogError("CSCDigiDump") << "Cannot get CorrelatedLCTs by label " << inputTag_.encode();
   }
 
   unsigned nDigisPerEvent = 0;
 
   for (auto j = lcts->begin(); j != lcts->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 = lcts->get((*j).first);
     // 1=forward (+Z); 2=backward (-Z)
     const unsigned typeCorrected(detId.endcap() == 1 ? chamberType - 2 : chamberType - 2 + chambers_.size());
     for (auto lct = range.first; lct != range.second; lct++) {
       if (lct->isValid()) {
         chamberHistos[typeCorrected]["lct_pattern"]->Fill(lct->getPattern());
         chamberHistos[typeCorrected]["lct_quality"]->Fill(lct->getQuality());
         chamberHistos[typeCorrected]["lct_wiregroup"]->Fill(lct->getKeyWG());
         chamberHistos[typeCorrected]["lct_halfstrip"]->Fill(lct->getStrip());
         chamberHistos[typeCorrected]["lct_bend"]->Fill(lct->getBend());
         chamberHistos[typeCorrected]["lct_bx"]->Fill(lct->getBX());
         chamberHistos[typeCorrected]["lct_type"]->Fill(lct->getType());
         if (isRun3_) {
           // ignore these fields for chambers that do not enable the Run-3 algorithm
           if (std::find(chambersRun3_.begin(), chambersRun3_.end(), chamberType - 2) == chambersRun3_.end())
             continue;
           chamberHistos[typeCorrected]["lct_run3pattern"]->Fill(lct->getRun3Pattern());
           chamberHistos[typeCorrected]["lct_slope"]->Fill(lct->getSlope());
           chamberHistos[typeCorrected]["lct_quartstrip"]->Fill(lct->getStrip(4));
           chamberHistos[typeCorrected]["lct_eighthstrip"]->Fill(lct->getStrip(8));
         }
       }
     }
   }
   theNDigisPerEventPlot->Fill(nDigisPerEvent);
 }

◆ bookHistograms()

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

Definition at line 28 of file CSCCorrelatedLCTDigiValidation.cc.

References dqm::implementation::IBooker::book1D(), chamberHistos, CSCDetId::chamberName(), chambers_, chambersRun3_, spr::find(), B2GTnPMonitor_cfi::histTitle, mps_fire::i, isRun3_, crabWrapper::key, lctMaxBin_, lctMinBin_, lctNBin_, lctVars_, dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, RandomServiceHelper::t2, theNDigisPerChamberPlots, and theNDigisPerEventPlot.

                                                                             {
   iBooker.setCurrentFolder("MuonCSCDigisV/CSCDigiTask/LCT/Occupancy/");
 
   theNDigisPerEventPlot =
       iBooker.book1D("CSCCorrelatedLCTDigisPerEvent",
                      "CorrelatedLCT trigger primitives per event; Number of CorrelatedLCTs; Entries",
                      100,
                      0,
                      100);
   for (int i = 1; i <= 10; ++i) {
     const std::string t2("CSCCorrelatedLCTDigisPerChamber_" + CSCDetId::chamberName(i));
     theNDigisPerChamberPlots[i - 1] =
         iBooker.book1D(t2,
                        "Number of CorrelatedLCTs per chamber " + CSCDetId::chamberName(i) +
                            ";Number of CorrelatedLCTs per chamber;Entries",
                        4,
                        0,
                        4);
   }
 
   // do not analyze Run-3 properties in Run-1 and Run-2 eras
   if (!isRun3_) {
     lctVars_.resize(7);
   }
 
   // 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 ? "+" : "-");
       // lct variable
       for (unsigned iVar = 0; iVar < lctVars_.size(); iVar++) {
         if (std::find(chambersRun3_.begin(), chambersRun3_.end(), iType) == chambersRun3_.end()) {
           if (iVar > 6)
             continue;
         }
         const std::string key("lct_" + lctVars_[iVar]);
         const std::string histName(key + "_" + chambers_[iType] + eSign);
         const std::string histTitle(chambers_[iType] + eSign + " LCT " + lctVars_[iVar]);
         const unsigned iTypeCorrected(iEndcap == 0 ? iType : iType + chambers_.size());
         chamberHistos[iTypeCorrected][key] =
             iBooker.book1D(histName, histTitle, lctNBin_[iVar], lctMinBin_[iVar], lctMaxBin_[iVar]);
         chamberHistos[iTypeCorrected][key]->getTH1()->SetMinimum(0);
         // set bin labels for the "type" plot. very useful in ME1/1 and ME2/1
         // when the GEM-CSC ILTs will be running
         if (lctVars_[iVar] == "type") {
           chamberHistos[iTypeCorrected][key]->setBinLabel(1, "CLCTALCT", 1);
           chamberHistos[iTypeCorrected][key]->setBinLabel(2, "ALCTCLCT", 1);
           chamberHistos[iTypeCorrected][key]->setBinLabel(3, "ALCTCLCTGEM", 1);
           chamberHistos[iTypeCorrected][key]->setBinLabel(4, "ALCTCLCT2GEM", 1);
           chamberHistos[iTypeCorrected][key]->setBinLabel(5, "ALCT2GEM", 1);
           chamberHistos[iTypeCorrected][key]->setBinLabel(6, "CLCT2GEM", 1);
           chamberHistos[iTypeCorrected][key]->setBinLabel(7, "CLCTONLY", 1);
           chamberHistos[iTypeCorrected][key]->setBinLabel(8, "ALCTONLY", 1);
         }
       }
     }
   }
 }