#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)
	~CSCComparatorDigiValidation ()
Private Attributes
MonitorElement *	the3StripPlots [10]
MonitorElement *	theNDigisPerEventPlot
MonitorElement *	theNDigisPerLayerPlots [10]
edm::InputTag	theStripDigiInputTag
MonitorElement *	theStripDigiPlots [10]
MonitorElement *	theTimeBinPlots [10]

Detailed Description

Definition at line 8 of file CSCComparatorDigiValidation.h.

Constructor & Destructor Documentation

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

Definition at line 9 of file CSCComparatorDigiValidation.cc.

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

: CSCBaseValidation(dbe, inputTag),
  theStripDigiInputTag(stripDigiInputTag),
  theTimeBinPlots(),
  theNDigisPerLayerPlots(),
  theStripDigiPlots(),
  the3StripPlots(),
  theNDigisPerEventPlot( dbe_->book1D("CSCComparatorDigisPerEvent", "CSC Comparator Digis per event", 100, 0, 100) )
{
  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 36 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 49 of file CSCComparatorDigiValidation.cc.

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

Referenced by CSCDigiValidation::analyze().

{
  edm::Handle<CSCComparatorDigiCollection> comparators;
  edm::Handle<CSCStripDigiCollection> stripDigis;

  e.getByLabel(theInputTag, comparators);
  if (!comparators.isValid()) {
    edm::LogError("CSCDigiDump") << "Cannot get comparators by label " << theInputTag.encode();
  }
  e.getByLabel(theStripDigiInputTag, 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);
}