CSCStripDigiValidation Class Reference

#include <CSCStripDigiValidation.h>

Inheritance diagram for CSCStripDigiValidation:

Public Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &)
void	bookHistograms (DQMStore::IBooker &, bool doSim)
	CSCStripDigiValidation (const edm::InputTag &inputTag, edm::ConsumesCollector &&iC)
void	plotResolution (const PSimHit &hit, int strip, const CSCLayer *layer, int chamberType)
void	setGeometry (const CSCGeometry *geom)
	~CSCStripDigiValidation ()
Public Member Functions inherited from CSCBaseValidation
	CSCBaseValidation (const edm::InputTag &inputTag)
void	setGeometry (const CSCGeometry *geom)
void	setSimHitMap (const PSimHitMap *simHitMap)
virtual	~CSCBaseValidation ()

Private Member Functions
void	fillPedestalPlots (const CSCStripDigi &digi)
void	fillSignalPlots (const CSCStripDigi &digi)

Private Attributes
edm::EDGetTokenT < CSCStripDigiCollection >	strips_Token_
MonitorElement *	theAmplitudePlot
MonitorElement *	theNDigisPerChamberPlot
MonitorElement *	theNDigisPerEventPlot
MonitorElement *	theNDigisPerLayerPlot
int	thePedestalCount
float	thePedestalCovarianceSum
MonitorElement *	thePedestalNeighborCorrelationPlot
MonitorElement *	thePedestalPlot
float	thePedestalSum
MonitorElement *	thePedestalTimeCorrelationPlot
MonitorElement *	theRatio4to5Plot
MonitorElement *	theRatio6to5Plot
MonitorElement *	theResolutionPlots [10]

Additional Inherited Members
Protected Member Functions inherited from CSCBaseValidation
const CSCLayer *	findLayer (int detId) const
Protected Attributes inherited from CSCBaseValidation
const CSCGeometry *	theCSCGeometry
edm::InputTag	theInputTag
const PSimHitMap *	theSimHitMap

Detailed Description

Definition at line 12 of file CSCStripDigiValidation.h.

Constructor & Destructor Documentation

CSCStripDigiValidation::CSCStripDigiValidation	(	const edm::InputTag &	inputTag,
		edm::ConsumesCollector &&	iC
	)

Definition at line 7 of file CSCStripDigiValidation.cc.

References beam_dqm_sourceclient-live_cfg::inputTag, and strips_Token_.

                                                                          :
  CSCBaseValidation(inputTag),
  thePedestalSum(0),
  thePedestalCovarianceSum(0),
  thePedestalCount(0),
  thePedestalTimeCorrelationPlot(0),
  thePedestalNeighborCorrelationPlot(0),
  theNDigisPerChamberPlot(0)
{
  strips_Token_ = iC.consumes<CSCStripDigiCollection>(inputTag);
}

CSCStripDigiValidation::~CSCStripDigiValidation

(

)

Definition at line 20 of file CSCStripDigiValidation.cc.

{
}

Member Function Documentation

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

virtual

Implements CSCBaseValidation.

Definition at line 43 of file CSCStripDigiValidation.cc.

References edm::InputTag::encode(), MonitorElement::Fill(), fillPedestalPlots(), fillSignalPlots(), edm::Event::getByToken(), edm::HandleBase::isValid(), j, prof2calltree::last, EcalCondDBWriter_cfi::pedestal, RecoTauPiZeroBuilderPlugins_cfi::strips, strips_Token_, CSCBaseValidation::theInputTag, theNDigisPerEventPlot, theNDigisPerLayerPlot, thePedestalCount, and thePedestalSum.

Referenced by CSCDigiValidation::analyze().

{
  edm::Handle<CSCStripDigiCollection> strips;
  e.getByToken(strips_Token_, strips);
  if (!strips.isValid()) {
    edm::LogError("CSCDigiValidation") << "Cannot get strips by label "
                                       << theInputTag.encode();
  }

 unsigned nDigisPerEvent = 0;

 for (CSCStripDigiCollection::DigiRangeIterator j=strips->begin(); j!=strips->end(); j++) {
    std::vector<CSCStripDigi>::const_iterator digiItr = (*j).second.first;
    std::vector<CSCStripDigi>::const_iterator last = (*j).second.second;
    int nDigis = last-digiItr;
    nDigisPerEvent += nDigis;
    theNDigisPerLayerPlot->Fill(nDigis);

    double maxAmplitude = 0.;
    // int maxStrip = 0;

    for( ; digiItr != last; ++digiItr) {
      // average up the pedestals
      std::vector<int> adcCounts = digiItr->getADCCounts();
      thePedestalSum += adcCounts[0];
      thePedestalSum += adcCounts[1];
      thePedestalCount += 2;
      float pedestal = thePedestalSum/thePedestalCount;
      if(adcCounts[4]-pedestal > maxAmplitude)
      {
      //  maxStrip = digiItr->getStrip();
        maxAmplitude = adcCounts[4]-pedestal;
      }

      // if we have enough pedestal statistics
      if(thePedestalCount > 100)
      {
        fillPedestalPlots(*digiItr);

        // see if it's big enough to count as "signal"
        if(adcCounts[5] > (thePedestalSum/thePedestalCount + 100))
        {
          fillSignalPlots(*digiItr);
        }
      }
    }
  } // loop over digis

  theNDigisPerEventPlot->Fill(nDigisPerEvent);
}

void CSCStripDigiValidation::bookHistograms	(	DQMStore::IBooker &	iBooker,
		bool	doSim
	)

Definition at line 24 of file CSCStripDigiValidation.cc.

References DQMStore::IBooker::book1D(), i, theAmplitudePlot, theNDigisPerEventPlot, theNDigisPerLayerPlot, thePedestalPlot, theRatio4to5Plot, theRatio6to5Plot, and theResolutionPlots.

Referenced by CSCDigiValidation::bookHistograms().

{
  thePedestalPlot = iBooker.book1D("CSCPedestal", "CSC Pedestal ", 400, 550, 650);
  theAmplitudePlot = iBooker.book1D("CSCStripAmplitude", "CSC Strip Amplitude", 200, 0, 2000);
  theRatio4to5Plot = iBooker.book1D("CSCStrip4to5", "CSC Strip Ratio tbin 4 to tbin 5", 100, 0, 1);
  theRatio6to5Plot = iBooker.book1D("CSCStrip6to5", "CSC Strip Ratio tbin 6 to tbin 5", 120, 0, 1.2);
  theNDigisPerLayerPlot = iBooker.book1D("CSCStripDigisPerLayer", "Number of CSC Strip Digis per layer", 48, 0, 48);
  theNDigisPerEventPlot = iBooker.book1D("CSCStripDigisPerEvent", "Number of CSC Strip Digis per event", 100, 0, 500); 
  if(doSim) {
    for(int i = 0; i < 10; ++i)
    {
      char title1[200];
      sprintf(title1, "CSCStripDigiResolution%d", i+1);
      theResolutionPlots[i] = iBooker.book1D(title1, title1, 100, -5, 5);
    }
  }
}

void CSCStripDigiValidation::fillPedestalPlots ( const CSCStripDigi &  digi )

private

Definition at line 96 of file CSCStripDigiValidation.cc.

References MonitorElement::Fill(), CSCStripDigi::getADCCounts(), and thePedestalPlot.

Referenced by analyze().

{
  std::vector<int> adcCounts = digi.getADCCounts();
  thePedestalPlot->Fill(adcCounts[0]);
  thePedestalPlot->Fill(adcCounts[1]);
}

void CSCStripDigiValidation::fillSignalPlots ( const CSCStripDigi &  digi )

private

Definition at line 105 of file CSCStripDigiValidation.cc.

References MonitorElement::Fill(), CSCStripDigi::getADCCounts(), EcalCondDBWriter_cfi::pedestal, theAmplitudePlot, thePedestalCount, thePedestalSum, theRatio4to5Plot, and theRatio6to5Plot.

Referenced by analyze().

{
  std::vector<int> adcCounts = digi.getADCCounts();
  float pedestal = thePedestalSum/thePedestalCount;
  theAmplitudePlot->Fill(adcCounts[4] - pedestal);
  theRatio4to5Plot->Fill( (adcCounts[3]-pedestal) / (adcCounts[4]-pedestal) );
  theRatio6to5Plot->Fill( (adcCounts[5]-pedestal) / (adcCounts[4]-pedestal) );
}

void CSCStripDigiValidation::plotResolution	(	const PSimHit &	hit,
		int	strip,
		const CSCLayer *	layer,
		int	chamberType
	)

Definition at line 115 of file CSCStripDigiValidation.cc.

References MonitorElement::Fill(), CSCLayer::geometry(), PSimHit::localPosition(), theResolutionPlots, PV3DBase< T, PVType, FrameType >::x(), CSCLayerGeometry::xOfStrip(), and PV3DBase< T, PVType, FrameType >::y().

{
  double hitX = hit.localPosition().x();
  double hitY = hit.localPosition().y();
  double digiX = layer->geometry()->xOfStrip(strip, hitY);
  theResolutionPlots[chamberType-1]->Fill(digiX - hitX);
}

void CSCStripDigiValidation::setGeometry ( const CSCGeometry *  geom )

inline

Definition at line 20 of file CSCStripDigiValidation.h.

References relativeConstraints::geom, and CSCBaseValidation::theCSCGeometry.

Referenced by CSCDigiValidation::analyze().

{theCSCGeometry = geom;}

Member Data Documentation

edm::EDGetTokenT<CSCStripDigiCollection> CSCStripDigiValidation::strips_Token_

private

Definition at line 28 of file CSCStripDigiValidation.h.

Referenced by analyze(), and CSCStripDigiValidation().

MonitorElement* CSCStripDigiValidation::theAmplitudePlot

private

Definition at line 35 of file CSCStripDigiValidation.h.

Referenced by bookHistograms(), and fillSignalPlots().

MonitorElement* CSCStripDigiValidation::theNDigisPerChamberPlot

private

Definition at line 39 of file CSCStripDigiValidation.h.

MonitorElement* CSCStripDigiValidation::theNDigisPerEventPlot

private

Definition at line 40 of file CSCStripDigiValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CSCStripDigiValidation::theNDigisPerLayerPlot

private

Definition at line 38 of file CSCStripDigiValidation.h.

Referenced by analyze(), and bookHistograms().

int CSCStripDigiValidation::thePedestalCount

private

Definition at line 31 of file CSCStripDigiValidation.h.

Referenced by analyze(), and fillSignalPlots().

float CSCStripDigiValidation::thePedestalCovarianceSum

private

Definition at line 30 of file CSCStripDigiValidation.h.

MonitorElement* CSCStripDigiValidation::thePedestalNeighborCorrelationPlot

private

Definition at line 34 of file CSCStripDigiValidation.h.

MonitorElement* CSCStripDigiValidation::thePedestalPlot

private

Definition at line 32 of file CSCStripDigiValidation.h.

Referenced by bookHistograms(), and fillPedestalPlots().

float CSCStripDigiValidation::thePedestalSum

private

Definition at line 29 of file CSCStripDigiValidation.h.

Referenced by analyze(), and fillSignalPlots().

MonitorElement* CSCStripDigiValidation::thePedestalTimeCorrelationPlot

private

Definition at line 33 of file CSCStripDigiValidation.h.

MonitorElement* CSCStripDigiValidation::theRatio4to5Plot

private

Definition at line 36 of file CSCStripDigiValidation.h.

Referenced by bookHistograms(), and fillSignalPlots().

MonitorElement* CSCStripDigiValidation::theRatio6to5Plot

private

Definition at line 37 of file CSCStripDigiValidation.h.

Referenced by bookHistograms(), and fillSignalPlots().

MonitorElement* CSCStripDigiValidation::theResolutionPlots[10]

private

Definition at line 41 of file CSCStripDigiValidation.h.

Referenced by bookHistograms(), and plotResolution().