#include <CSCComparatorDigiValidation.h>

Inheritance diagram for CSCComparatorDigiValidation:

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

	CSCComparatorDigiValidation (DQMStore *dbe, const edm::InputTag &inputTag, const edm::InputTag &stripDigiInputTag, edm::ConsumesCollector &&iC)

	~CSCComparatorDigiValidation ()

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
edm::EDGetTokenT < CSCComparatorDigiCollection >	comparators_Token_

MonitorElement *	the3StripPlots [10]

MonitorElement *	theNDigisPerEventPlot

MonitorElement *	theNDigisPerLayerPlots [10]

edm::EDGetTokenT < CSCStripDigiCollection >	theStripDigi_Token_

MonitorElement *	theStripDigiPlots [10]

MonitorElement *	theTimeBinPlots [10]

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 13 of file CSCComparatorDigiValidation.h.

Constructor & Destructor Documentation

CSCComparatorDigiValidation::CSCComparatorDigiValidation	(	DQMStore *	dbe,
		const edm::InputTag &	inputTag,
		const edm::InputTag &	stripDigiInputTag,
		edm::ConsumesCollector &&	iC
	)

Definition at line 7 of file CSCComparatorDigiValidation.cc.

References DQMStore::book1D(), comparators_Token_, CSCBaseValidation::dbe_, i, the3StripPlots, theNDigisPerLayerPlots, theStripDigiPlots, and theTimeBinPlots.

 : CSCBaseValidation(dbe, inputTag),
   theStripDigi_Token_(iC.consumes<CSCStripDigiCollection>(stripDigiInputTag)),
   theTimeBinPlots(),
   theNDigisPerLayerPlots(),
   theStripDigiPlots(),
   the3StripPlots(),
   theNDigisPerEventPlot( dbe_->book1D("CSCComparatorDigisPerEvent", "CSC Comparator Digis per event", 100, 0, 100) )
 {
   comparators_Token_ = iC.consumes<CSCComparatorDigiCollection>(inputTag);
 
   for(int i = 0; i < 10; ++i)
   {
     char title1[200], title2[200], title3[200], title4[200];
     sprintf(title1, "CSCComparatorDigiTimeType%d", i+1);
     sprintf(title2, "CSCComparatorDigisPerLayerType%d", i+1);
     sprintf(title3, "CSCComparatorStripAmplitudeType%d", i+1);
     sprintf(title4, "CSCComparator3StripAmplitudeType%d", i+1);
     theTimeBinPlots[i] = dbe_->book1D(title1, title1, 9, 0, 8);
     theNDigisPerLayerPlots[i] = dbe_->book1D(title2, title2, 100, 0, 20);
     theStripDigiPlots[i] = dbe_->book1D(title3, title3, 100, 0, 1000);
     the3StripPlots[i] = dbe_->book1D(title4, title4, 100, 0, 1000);
   }
 }

CSCComparatorDigiValidation::~CSCComparatorDigiValidation ( )

Definition at line 37 of file CSCComparatorDigiValidation.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") << "Mean charge of " << the3StripPlots[i]->getName()
 //       << " is " << the3StripPlots[i]->getMean();
 //   }
 }

Member Function Documentation

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

virtual

Implements CSCBaseValidation.

Definition at line 50 of file CSCComparatorDigiValidation.cc.

References ecalMGPA::adc(), CSCLayer::chamber(), CSCChamberSpecs::chamberType(), comparators_Token_, edm::InputTag::encode(), MonitorElement::Fill(), CSCBaseValidation::findLayer(), edm::Event::getByToken(), edm::HandleBase::isValid(), j, prof2calltree::last, DetId::rawId(), CSCChamber::specs(), the3StripPlots, CSCBaseValidation::theInputTag, theNDigisPerEventPlot, theNDigisPerLayerPlots, theStripDigi_Token_, theStripDigiPlots, and theTimeBinPlots.

Referenced by CSCDigiValidation::analyze().

 {
   edm::Handle<CSCComparatorDigiCollection> comparators;
   edm::Handle<CSCStripDigiCollection> stripDigis;
 
   e.getByToken(comparators_Token_, comparators);
   if (!comparators.isValid()) {
     edm::LogError("CSCDigiDump") << "Cannot get comparators by label " << theInputTag.encode();
   }
   e.getByToken(theStripDigi_Token_, stripDigis);
   if (!stripDigis.isValid()) {
     edm::LogError("CSCDigiDump") << "Cannot get comparators by label " << theInputTag.encode();
   }
 
   unsigned nDigisPerEvent = 0;
 
   for (CSCComparatorDigiCollection::DigiRangeIterator j=comparators->begin(); j!=comparators->end(); j++) {
     std::vector<CSCComparatorDigi>::const_iterator digiItr = (*j).second.first;
     std::vector<CSCComparatorDigi>::const_iterator last = (*j).second.second;
 
     CSCDetId detId((*j).first);
     const CSCLayer * layer = findLayer(detId.rawId());
     int chamberType = layer->chamber()->specs()->chamberType();
 
     CSCStripDigiCollection::Range stripRange = stripDigis->get(detId);
 
     theNDigisPerLayerPlots[chamberType-1]->Fill(last-digiItr);
 
     for( ; digiItr != last; ++digiItr) {
       ++nDigisPerEvent;
       theTimeBinPlots[chamberType-1]->Fill(digiItr->getTimeBin());
 
       int strip = digiItr->getStrip();
       for(std::vector<CSCStripDigi>::const_iterator stripItr = stripRange.first;
           stripItr != stripRange.second; ++stripItr)
       {
         if(stripItr->getStrip() == strip)
         {
           std::vector<int> adc = stripItr->getADCCounts();
           float pedc = 0.5*(adc[0]+adc[1]);
           float amp = adc[4] - pedc;
           theStripDigiPlots[chamberType-1]->Fill(amp);
           // check neighbors
           if(stripItr != stripRange.first && stripItr != stripRange.second-1)
           {
             std::vector<int> adcl = (stripItr-1)->getADCCounts();
             std::vector<int> adcr = (stripItr+1)->getADCCounts();
             float pedl = 0.5*(adcl[0]+adcl[1]);
             float pedr = 0.5*(adcr[0]+adcr[1]);
             float three = adcl[4]-pedl
                         + adcr[4]-pedr
                         + amp;
             the3StripPlots[chamberType-1]->Fill(three);
           }
         }
       }
     }
 
   }
 
   theNDigisPerEventPlot->Fill(nDigisPerEvent);
 }