#include <CSCComparatorDigiValidation.h>

Inheritance diagram for CSCComparatorDigiValidation:

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

void	bookHistograms (DQMStore::IBooker &)

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

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

	~CSCComparatorDigiValidation () 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::EDGetTokenT < CSCComparatorDigiCollection >	comparators_Token_

edm::InputTag	inputTagComp_

edm::InputTag	inputTagStrip_

edm::EDGetTokenT < CSCStripDigiCollection >	strips_Token_

MonitorElement *	the3StripPlots [10]

MonitorElement *	theNDigisPerEventPlot

MonitorElement *	theNDigisPerLayerPlots [10]

MonitorElement *	theResolutionPlots [10]

MonitorElement *	theStripDigiPlots [10]

MonitorElement *	theTimeBinPlots [10]

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

Constructor & Destructor Documentation

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

Definition at line 7 of file CSCComparatorDigiValidation.cc.

References comparators_Token_, edm::ParameterSet::getParameterSet(), inputTagComp_, inputTagStrip_, DigiDM_cff::strips, and strips_Token_.

     : CSCBaseValidation(ps), theTimeBinPlots(), theNDigisPerLayerPlots(), theStripDigiPlots(), the3StripPlots() {
   const auto &comps = ps.getParameterSet("cscComparatorDigi");
   inputTagComp_ = comps.getParameter<edm::InputTag>("inputTag");
   comparators_Token_ = iC.consumes<CSCComparatorDigiCollection>(inputTagComp_);
 
   const auto &strips = ps.getParameterSet("cscStripDigi");
   inputTagStrip_ = strips.getParameter<edm::InputTag>("inputTag");
   strips_Token_ = iC.consumes<CSCStripDigiCollection>(inputTagStrip_);
 }

CSCComparatorDigiValidation::~CSCComparatorDigiValidation ( )

override

Definition at line 18 of file CSCComparatorDigiValidation.cc.

18 {}

Member Function Documentation

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

overridevirtual

Implements CSCBaseValidation.

Definition at line 65 of file CSCComparatorDigiValidation.cc.

References gpuClustering::adc, CSCLayer::chamber(), CSCChamberSpecs::chamberType(), comparators_Token_, CSCBaseValidation::doSim_, edm::InputTag::encode(), dqm::impl::MonitorElement::Fill(), CSCBaseValidation::findLayer(), edm::Event::getByToken(), PSimHitMap::hits(), inputTagComp_, edm::HandleBase::isValid(), dqmiolumiharvest::j, dqmdumpme::last, phase1PixelTopology::layer, CSCConstants::NUM_STRIPS_PER_CFEB, plotResolution(), DetId::rawId(), trackerHits::simHits, CSCChamber::specs(), digitizers_cfi::strip, strips_Token_, the3StripPlots, theNDigisPerEventPlot, theNDigisPerLayerPlots, CSCBaseValidation::theSimHitMap, theStripDigiPlots, and theTimeBinPlots.

                                                                                   {
   edm::Handle<CSCComparatorDigiCollection> comparators;
   edm::Handle<CSCStripDigiCollection> stripDigis;
 
   e.getByToken(comparators_Token_, comparators);
   if (!comparators.isValid()) {
     edm::LogError("CSCComparatorDigiValidation") << "Cannot get comparators by label " << inputTagComp_.encode();
   }
   e.getByToken(strips_Token_, stripDigis);
   if (!stripDigis.isValid()) {
     edm::LogError("CSCComparatorDigiValidation") << "Cannot get strips by label " << inputTagComp_.encode();
   }
 
   unsigned nDigisPerEvent = 0;
 
   for (auto j = comparators->begin(); j != comparators->end(); j++) {
     auto digiItr = (*j).second.first;
     auto last = (*j).second.second;
 
     CSCDetId detId((*j).first);
     const CSCLayer *layer = findLayer(detId.rawId());
     int chamberType = layer->chamber()->specs()->chamberType();
     int nDigis = last - digiItr;
 
     theNDigisPerLayerPlots[chamberType - 1]->Fill(last - digiItr);
 
     // require exactly one comparator digi per layer
     if (doSim_) {
       const edm::PSimHitContainer &simHits = theSimHitMap->hits(detId);
       if (nDigis == 1 && simHits.size() == 1) {
         const int fullStrip(digiItr->getStrip() + digiItr->getCFEB() * CSCConstants::NUM_STRIPS_PER_CFEB);
         plotResolution(simHits[0], fullStrip, layer, chamberType);
       }
     }
 
     for (auto stripRange = stripDigis->get(detId); digiItr != last; ++digiItr) {
       ++nDigisPerEvent;
       theTimeBinPlots[chamberType - 1]->Fill(digiItr->getTimeBin());
 
       int strip = digiItr->getStrip();
       for (auto 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);
 }

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

Definition at line 20 of file CSCComparatorDigiValidation.cc.

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

                                                                          {
   iBooker.setCurrentFolder("MuonCSCDigisV/CSCDigiTask/Comparator/Occupancy/");
 
   theNDigisPerEventPlot = iBooker.book1D("CSCComparatorDigisPerEvent",
                                          "CSC Comparator Digis per event;CSC Comparator Digis per event;Entries",
                                          100,
                                          0,
                                          100);
   // 10 chamber types, if you consider ME1/a and ME1/b separate
   for (int i = 1; i <= 10; ++i) {
     const std::string t1("CSCComparatorDigiTime_" + CSCDetId::chamberName(i));
     const std::string t2("CSCComparatorDigisPerLayer_" + CSCDetId::chamberName(i));
     const std::string t3("CSCComparatorStripAmplitude_" + CSCDetId::chamberName(i));
     const std::string t4("CSCComparator3StripAmplitude_" + CSCDetId::chamberName(i));
     theTimeBinPlots[i - 1] = iBooker.book1D(
         t1, "Comparator Time Bin " + CSCDetId::chamberName(i) + " ;Comparator Time Bin; Entries", 16, 0, 16);
     theNDigisPerLayerPlots[i - 1] = iBooker.book1D(
         t2,
         "Number of Comparator Digis " + CSCDetId::chamberName(i) + " ;Number of Comparator Digis; Entries",
         100,
         0,
         20);
     theStripDigiPlots[i - 1] = iBooker.book1D(
         t3, "Comparator Amplitude " + CSCDetId::chamberName(i) + " ;Comparator Amplitude; Entries", 100, 0, 1000);
     the3StripPlots[i - 1] = iBooker.book1D(
         t4,
         "Comparator-triplet Amplitude " + CSCDetId::chamberName(i) + " ;Comparator-triplet Amplitude; Entries",
         100,
         0,
         1000);
   }
   if (doSim_) {
     iBooker.setCurrentFolder("MuonCSCDigisV/CSCDigiTask/Comparator/Resolution/");
     for (int i = 1; i <= 10; ++i) {
       const std::string t1("CSCComparatorPosResolution_" + CSCDetId::chamberName(i));
       theResolutionPlots[i - 1] = iBooker.book1D(t1,
                                                  "Comparator X Position Resolution " + CSCDetId::chamberName(i) +
                                                      ";Comparator X Position Resolution [cm]; Entries",
                                                  100,
                                                  -5,
                                                  5);
     }
   }
 }

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

Definition at line 127 of file CSCComparatorDigiValidation.cc.

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

Referenced by analyze().

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