#include <CSCWireDigiValidation.h>

Inheritance diagram for CSCWireDigiValidation:

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

	CSCWireDigiValidation (DQMStore *dbe, const edm::InputTag &inputTag, edm::ConsumesCollector &&iC, bool doSim)

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

	~CSCWireDigiValidation ()

Public Member Functions inherited from CSCBaseValidation
	CSCBaseValidation (DQMStore *dbe, const edm::InputTag &inputTag)

void	setGeometry (const CSCGeometry *geom)

void	setSimHitMap (const PSimHitMap *simHitMap)

virtual	~CSCBaseValidation ()

Private Attributes
bool	theDoSimFlag

MonitorElement *	theNDigisPerEventPlot

MonitorElement *	theNDigisPerLayerPlots [10]

MonitorElement *	theResolutionPlots [10]

MonitorElement *	theTimeBinPlots [10]

edm::EDGetTokenT < CSCWireDigiCollection >	wires_Token_

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

Protected Attributes inherited from CSCBaseValidation
DQMStore *	dbe_

const CSCGeometry *	theCSCGeometry

edm::InputTag	theInputTag

const PSimHitMap *	theSimHitMap

Detailed Description

Definition at line 12 of file CSCWireDigiValidation.h.

Constructor & Destructor Documentation

CSCWireDigiValidation::CSCWireDigiValidation	(	DQMStore *	dbe,
		const edm::InputTag &	inputTag,
		edm::ConsumesCollector &&	iC,
		bool	doSim
	)

Definition at line 11 of file CSCWireDigiValidation.cc.

References DQMStore::book1D(), CSCBaseValidation::dbe_, i, theNDigisPerLayerPlots, theResolutionPlots, theTimeBinPlots, and wires_Token_.

 : CSCBaseValidation(dbe, inputTag),
   theDoSimFlag(doSim),
   theTimeBinPlots(),
   theNDigisPerLayerPlots(),
   theNDigisPerEventPlot( dbe_->book1D("CSCWireDigisPerEvent", "CSC Wire Digis per event", 100, 0, 100) )
 {
   wires_Token_ = iC.consumes<CSCWireDigiCollection>(inputTag);
 
   for(int i = 0; i < 10; ++i)
   {
     char title1[200], title2[200], title3[200];
     sprintf(title1, "CSCWireDigiTimeType%d", i+1);
     sprintf(title2, "CSCWireDigisPerLayerType%d", i+1);
     sprintf(title3, "CSCWireDigiResolution%d", i+1);
     theTimeBinPlots[i] = dbe_->book1D(title1, title1, 9, 0, 8);
     theNDigisPerLayerPlots[i] = dbe_->book1D(title2, title2, 100, 0, 20);
     theResolutionPlots[i] = dbe_->book1D(title3, title3, 100, -10, 10);
   }
 }

CSCWireDigiValidation::~CSCWireDigiValidation ( )

Definition at line 37 of file CSCWireDigiValidation.cc.

 {
 //   for(int i = 0; i < 10; ++i)
 //   {
 //     edm::LogInfo("CSCDigiValidation") << "Mean of " << theTimeBinPlots[i]->getName() 
 //       << " is " << theTimeBinPlots[i]->getMean() 
 //       << " +/- " << theTimeBinPlots[i]->getRMS();
 //     edm::LogInfo("CSCDigiValidation") << "RMS of " << theResolutionPlots[i]->getName() 
 //       << " is " << theResolutionPlots[i]->getRMS();
 //   }
 }

Member Function Documentation

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

virtual

Implements CSCBaseValidation.

Definition at line 50 of file CSCWireDigiValidation.cc.

References CSCLayer::chamber(), CSCChamberSpecs::chamberType(), edm::InputTag::encode(), MonitorElement::Fill(), CSCBaseValidation::findLayer(), edm::Event::getByToken(), PSimHitMap::hits(), edm::HandleBase::isValid(), j, plotResolution(), trackerHits::simHits, CSCChamber::specs(), theDoSimFlag, CSCBaseValidation::theInputTag, theNDigisPerEventPlot, theNDigisPerLayerPlots, CSCBaseValidation::theSimHitMap, theTimeBinPlots, and wires_Token_.

Referenced by CSCDigiValidation::analyze().

 {
   edm::Handle<CSCWireDigiCollection> wires;
 
   e.getByToken(wires_Token_, wires);
 
   if (!wires.isValid()) {
     edm::LogError("CSCDigiDump") << "Cannot get wires by label " << theInputTag.encode();
   }
 
   unsigned nDigisPerEvent = 0;
 
   for (CSCWireDigiCollection::DigiRangeIterator j=wires->begin(); j!=wires->end(); j++) {
     std::vector<CSCWireDigi>::const_iterator beginDigi = (*j).second.first;
     std::vector<CSCWireDigi>::const_iterator 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( std::vector<CSCWireDigi>::const_iterator digiItr = beginDigi;
          digiItr != endDigi; ++digiItr) 
     {
       theTimeBinPlots[chamberType-1]->Fill(digiItr->getTimeBin());
     }
 
     if(theDoSimFlag) 
     {
       const edm::PSimHitContainer simHits = theSimHitMap->hits(detId);
       if(nDigis == 1 && simHits.size() == 1)
       {
         plotResolution(simHits[0], *beginDigi, layer, chamberType);
       }
     }
   }
 
   theNDigisPerEventPlot->Fill(nDigisPerEvent);
 }

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

Definition at line 93 of file CSCWireDigiValidation.cc.

References 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);
 }