#include <CSCWireDigiValidation.h>

Inheritance diagram for CSCWireDigiValidation:

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

void	bookHistograms (DQMStore::IBooker &)

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

void	plotResolution (const PSimHit &hit, const CSCWireDigi &digi, const CSCLayer *layer, int chamberType)

	~CSCWireDigiValidation () 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
edm::InputTag	inputTag_

MonitorElement *	theNDigisPerEventPlot

MonitorElement *	theNDigisPerLayerPlots [10]

MonitorElement *	theResolutionPlots [10]

MonitorElement *	theTimeBinPlots [10]

edm::EDGetTokenT < CSCWireDigiCollection >	wires_Token_

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

Constructor & Destructor Documentation

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

Definition at line 9 of file CSCWireDigiValidation.cc.

References edm::ParameterSet::getParameterSet(), inputTag_, TrackValidation_cff::pset, and wires_Token_.

     : CSCBaseValidation(ps), theTimeBinPlots(), theNDigisPerLayerPlots() {
   const auto &pset = ps.getParameterSet("cscWireDigi");
   inputTag_ = pset.getParameter<edm::InputTag>("inputTag");
   wires_Token_ = iC.consumes<CSCWireDigiCollection>(inputTag_);
 }

CSCWireDigiValidation::~CSCWireDigiValidation ( )

override

Definition at line 16 of file CSCWireDigiValidation.cc.

16 {}

Member Function Documentation

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

overridevirtual

Implements CSCBaseValidation.

Definition at line 43 of file CSCWireDigiValidation.cc.

References CSCLayer::chamber(), CSCChamberSpecs::chamberType(), CSCBaseValidation::doSim_, edm::InputTag::encode(), dqm::impl::MonitorElement::Fill(), CSCBaseValidation::findLayer(), edm::Event::getByToken(), PSimHitMap::hits(), inputTag_, edm::HandleBase::isValid(), dqmiolumiharvest::j, phase1PixelTopology::layer, plotResolution(), trackerHits::simHits, CSCChamber::specs(), theNDigisPerEventPlot, theNDigisPerLayerPlots, CSCBaseValidation::theSimHitMap, theTimeBinPlots, DigiDM_cff::wires, and wires_Token_.

                                                                             {
   edm::Handle<CSCWireDigiCollection> wires;
 
   e.getByToken(wires_Token_, wires);
 
   if (!wires.isValid()) {
     edm::LogError("CSCWireDigiValidation") << "Cannot get wires by label " << inputTag_.encode();
   }
 
   unsigned nDigisPerEvent = 0;
 
   for (auto j = wires->begin(); j != wires->end(); j++) {
     auto beginDigi = (*j).second.first;
     auto endDigi = (*j).second.second;
     int detId = (*j).first.rawId();
 
     const CSCLayer *layer = findLayer(detId);
     int chamberType = layer->chamber()->specs()->chamberType();
     int nDigis = endDigi - beginDigi;
     nDigisPerEvent += nDigis;
     theNDigisPerLayerPlots[chamberType - 1]->Fill(nDigis);
 
     for (auto digiItr = beginDigi; digiItr != endDigi; ++digiItr) {
       theTimeBinPlots[chamberType - 1]->Fill(digiItr->getTimeBin());
     }
 
     if (doSim_) {
       const edm::PSimHitContainer &simHits = theSimHitMap->hits(detId);
       if (nDigis == 1 && simHits.size() == 1) {
         plotResolution(simHits[0], *beginDigi, layer, chamberType);
       }
     }
   }
 
   theNDigisPerEventPlot->Fill(nDigisPerEvent);
 }

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

Definition at line 18 of file CSCWireDigiValidation.cc.

References dqm::implementation::IBooker::book1D(), CSCDetId::chamberName(), mps_fire::i, dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, theNDigisPerEventPlot, theNDigisPerLayerPlots, theResolutionPlots, and theTimeBinPlots.

                                                                    {
   iBooker.setCurrentFolder("MuonCSCDigisV/CSCDigiTask/Wire/Occupancy/");
   theNDigisPerEventPlot =
       iBooker.book1D("CSCWireDigisPerEvent", "CSC Wire Digis per event;CSC Wire Digis per event;Entries", 100, 0, 100);
   for (int i = 1; i <= 10; ++i) {
     const std::string t1("CSCWireDigiTime_" + CSCDetId::chamberName(i));
     const std::string t2("CSCWireDigisPerLayer_" + CSCDetId::chamberName(i));
     theTimeBinPlots[i - 1] =
         iBooker.book1D(t1, "Wire Time Bin " + CSCDetId::chamberName(i) + ";Wire Time Bin; Entries", 16, 0, 16);
     theNDigisPerLayerPlots[i - 1] = iBooker.book1D(
         t2, "Number of Wire Digis " + CSCDetId::chamberName(i) + ";Number of Wire Digis; Entries", 100, 0, 20);
   }
 
   iBooker.setCurrentFolder("MuonCSCDigisV/CSCDigiTask/Wire/Resolution/");
   for (int i = 1; i <= 10; ++i) {
     const std::string t3("CSCWireDigiResolution_" + CSCDetId::chamberName(i));
     theResolutionPlots[i - 1] = iBooker.book1D(
         t3,
         "Wire Y Position Resolution " + CSCDetId::chamberName(i) + ";Wire Y Position Resolution [cm]; Entries",
         100,
         -10,
         10);
   }
 }

void CSCWireDigiValidation::plotResolution	(	const PSimHit &	hit,
		const CSCWireDigi &	digi,
		const CSCLayer *	layer,
		int	chamberType
	)

Definition at line 80 of file CSCWireDigiValidation.cc.

References dqm::impl::MonitorElement::Fill(), CSCLayer::geometry(), CSCWireDigi::getWireGroup(), PSimHit::localPosition(), theResolutionPlots, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and CSCLayerGeometry::yOfWireGroup().

Referenced by analyze().

                                                             {
   double hitX = hit.localPosition().x();
   double hitY = hit.localPosition().y();
   double digiY = layer->geometry()->yOfWireGroup(digi.getWireGroup(), hitX);
   theResolutionPlots[chamberType - 1]->Fill(digiY - hitY);
 }