CMS 3D CMS Logo

List of all members | Public Member Functions | Private Member Functions | Private Attributes
CSCStripDigiValidation Class Reference

#include <CSCStripDigiValidation.h>

Inheritance diagram for CSCStripDigiValidation:
CSCBaseValidation

Public Member Functions

void analyze (const edm::Event &e, const edm::EventSetup &) override
 
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 () override
 
- 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< CSCStripDigiCollectionstrips_Token_
 
MonitorElementtheAmplitudePlot
 
MonitorElementtheNDigisPerChamberPlot
 
MonitorElementtheNDigisPerEventPlot
 
MonitorElementtheNDigisPerLayerPlot
 
int thePedestalCount
 
float thePedestalCovarianceSum
 
MonitorElementthePedestalNeighborCorrelationPlot
 
MonitorElementthePedestalPlot
 
float thePedestalSum
 
MonitorElementthePedestalTimeCorrelationPlot
 
MonitorElementtheRatio4to5Plot
 
MonitorElementtheRatio6to5Plot
 
MonitorElementtheResolutionPlots [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 CSCLayerfindLayer (int detId) const
 
- Protected Attributes inherited from CSCBaseValidation
const CSCGeometrytheCSCGeometry
 
edm::InputTag theInputTag
 
const PSimHitMaptheSimHitMap
 

Detailed Description

Definition at line 11 of file CSCStripDigiValidation.h.

Constructor & Destructor Documentation

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

Definition at line 6 of file CSCStripDigiValidation.cc.

References SimL1EmulatorRepack_Full_cff::inputTag, and strips_Token_.

7  : CSCBaseValidation(inputTag),
13  theNDigisPerChamberPlot(nullptr) {
15 }
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
CSCBaseValidation(const edm::InputTag &inputTag)
edm::EDGetTokenT< CSCStripDigiCollection > strips_Token_
MonitorElement * theNDigisPerChamberPlot
MonitorElement * thePedestalTimeCorrelationPlot
MonitorElement * thePedestalNeighborCorrelationPlot
CSCStripDigiValidation::~CSCStripDigiValidation ( )
override

Definition at line 17 of file CSCStripDigiValidation.cc.

17 {}

Member Function Documentation

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

Implements CSCBaseValidation.

Definition at line 35 of file CSCStripDigiValidation.cc.

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

Referenced by CSCDigiValidation::analyze().

35  {
37  e.getByToken(strips_Token_, strips);
38  if (!strips.isValid()) {
39  edm::LogError("CSCDigiValidation") << "Cannot get strips by label " << theInputTag.encode();
40  }
41 
42  unsigned nDigisPerEvent = 0;
43 
44  for (CSCStripDigiCollection::DigiRangeIterator j = strips->begin(); j != strips->end(); j++) {
45  std::vector<CSCStripDigi>::const_iterator digiItr = (*j).second.first;
46  std::vector<CSCStripDigi>::const_iterator last = (*j).second.second;
47  int nDigis = last - digiItr;
48  nDigisPerEvent += nDigis;
49  theNDigisPerLayerPlot->Fill(nDigis);
50 
51  double maxAmplitude = 0.;
52  // int maxStrip = 0;
53 
54  for (; digiItr != last; ++digiItr) {
55  // average up the pedestals
56  std::vector<int> adcCounts = digiItr->getADCCounts();
57  thePedestalSum += adcCounts[0];
58  thePedestalSum += adcCounts[1];
59  thePedestalCount += 2;
61  if (adcCounts[4] - pedestal > maxAmplitude) {
62  // maxStrip = digiItr->getStrip();
63  maxAmplitude = adcCounts[4] - pedestal;
64  }
65 
66  // if we have enough pedestal statistics
67  if (thePedestalCount > 100) {
68  fillPedestalPlots(*digiItr);
69 
70  // see if it's big enough to count as "signal"
71  if (adcCounts[5] > (thePedestalSum / thePedestalCount + 100)) {
72  fillSignalPlots(*digiItr);
73  }
74  }
75  }
76  } // loop over digis
77 
78  theNDigisPerEventPlot->Fill(nDigisPerEvent);
79 }
edm::InputTag theInputTag
MonitorElement * theNDigisPerEventPlot
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:525
void fillSignalPlots(const CSCStripDigi &digi)
std::string encode() const
Definition: InputTag.cc:159
void Fill(long long x)
edm::EDGetTokenT< CSCStripDigiCollection > strips_Token_
MonitorElement * theNDigisPerLayerPlot
bool isValid() const
Definition: HandleBase.h:70
void fillPedestalPlots(const CSCStripDigi &digi)
strips
#turn off noise in all subdetectors simHcalUnsuppressedDigis.doNoise = False mix.digitizers.hcal.doNoise = False simEcalUnsuppressedDigis.doNoise = False mix.digitizers.ecal.doNoise = False simEcalUnsuppressedDigis.doESNoise = False simSiPixelDigis.AddNoise = False mix.digitizers.pixel.AddNoise = False simSiStripDigis.Noise = False mix.digitizers.strip.AddNoise = False
Definition: DigiDM_cff.py:32
void CSCStripDigiValidation::bookHistograms ( DQMStore::IBooker iBooker,
bool  doSim 
)

Definition at line 19 of file CSCStripDigiValidation.cc.

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

Referenced by CSCDigiValidation::bookHistograms().

19  {
20  thePedestalPlot = iBooker.book1D("CSCPedestal", "CSC Pedestal ", 400, 550, 650);
21  theAmplitudePlot = iBooker.book1D("CSCStripAmplitude", "CSC Strip Amplitude", 200, 0, 2000);
22  theRatio4to5Plot = iBooker.book1D("CSCStrip4to5", "CSC Strip Ratio tbin 4 to tbin 5", 100, 0, 1);
23  theRatio6to5Plot = iBooker.book1D("CSCStrip6to5", "CSC Strip Ratio tbin 6 to tbin 5", 120, 0, 1.2);
24  theNDigisPerLayerPlot = iBooker.book1D("CSCStripDigisPerLayer", "Number of CSC Strip Digis per layer", 48, 0, 48);
25  theNDigisPerEventPlot = iBooker.book1D("CSCStripDigisPerEvent", "Number of CSC Strip Digis per event", 100, 0, 500);
26  if (doSim) {
27  for (int i = 0; i < 10; ++i) {
28  char title1[200];
29  sprintf(title1, "CSCStripDigiResolution%d", i + 1);
30  theResolutionPlots[i] = iBooker.book1D(title1, title1, 100, -5, 5);
31  }
32  }
33 }
MonitorElement * theNDigisPerEventPlot
MonitorElement * theRatio4to5Plot
MonitorElement * theRatio6to5Plot
MonitorElement * thePedestalPlot
MonitorElement * theNDigisPerLayerPlot
MonitorElement * theAmplitudePlot
MonitorElement * theResolutionPlots[10]
void CSCStripDigiValidation::fillPedestalPlots ( const CSCStripDigi digi)
private

Definition at line 81 of file CSCStripDigiValidation.cc.

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

Referenced by analyze(), and setGeometry().

81  {
82  std::vector<int> adcCounts = digi.getADCCounts();
83  thePedestalPlot->Fill(adcCounts[0]);
84  thePedestalPlot->Fill(adcCounts[1]);
85 }
std::vector< int > const & getADCCounts() const
Get ADC readings.
Definition: CSCStripDigi.h:44
void Fill(long long x)
MonitorElement * thePedestalPlot
void CSCStripDigiValidation::fillSignalPlots ( const CSCStripDigi digi)
private

Definition at line 87 of file CSCStripDigiValidation.cc.

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

Referenced by analyze(), and setGeometry().

87  {
88  std::vector<int> adcCounts = digi.getADCCounts();
90  theAmplitudePlot->Fill(adcCounts[4] - pedestal);
91  theRatio4to5Plot->Fill((adcCounts[3] - pedestal) / (adcCounts[4] - pedestal));
92  theRatio6to5Plot->Fill((adcCounts[5] - pedestal) / (adcCounts[4] - pedestal));
93 }
MonitorElement * theRatio4to5Plot
std::vector< int > const & getADCCounts() const
Get ADC readings.
Definition: CSCStripDigi.h:44
MonitorElement * theRatio6to5Plot
void Fill(long long x)
MonitorElement * theAmplitudePlot
void CSCStripDigiValidation::plotResolution ( const PSimHit hit,
int  strip,
const CSCLayer layer,
int  chamberType 
)

Definition at line 95 of file CSCStripDigiValidation.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 setGeometry().

95  {
96  double hitX = hit.localPosition().x();
97  double hitY = hit.localPosition().y();
98  double digiX = layer->geometry()->xOfStrip(strip, hitY);
99  theResolutionPlots[chamberType - 1]->Fill(digiX - hitX);
100 }
T y() const
Definition: PV3DBase.h:60
float xOfStrip(int strip, float y=0.) const
void Fill(long long x)
Local3DPoint localPosition() const
Definition: PSimHit.h:52
MonitorElement * theResolutionPlots[10]
T x() const
Definition: PV3DBase.h:59
const CSCLayerGeometry * geometry() const
Definition: CSCLayer.h:44
void CSCStripDigiValidation::setGeometry ( const CSCGeometry geom)
inline

Member Data Documentation

edm::EDGetTokenT<CSCStripDigiCollection> CSCStripDigiValidation::strips_Token_
private

Definition at line 24 of file CSCStripDigiValidation.h.

Referenced by analyze(), and CSCStripDigiValidation().

MonitorElement* CSCStripDigiValidation::theAmplitudePlot
private

Definition at line 31 of file CSCStripDigiValidation.h.

Referenced by bookHistograms(), and fillSignalPlots().

MonitorElement* CSCStripDigiValidation::theNDigisPerChamberPlot
private

Definition at line 35 of file CSCStripDigiValidation.h.

MonitorElement* CSCStripDigiValidation::theNDigisPerEventPlot
private

Definition at line 36 of file CSCStripDigiValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CSCStripDigiValidation::theNDigisPerLayerPlot
private

Definition at line 34 of file CSCStripDigiValidation.h.

Referenced by analyze(), and bookHistograms().

int CSCStripDigiValidation::thePedestalCount
private

Definition at line 27 of file CSCStripDigiValidation.h.

Referenced by analyze(), and fillSignalPlots().

float CSCStripDigiValidation::thePedestalCovarianceSum
private

Definition at line 26 of file CSCStripDigiValidation.h.

MonitorElement* CSCStripDigiValidation::thePedestalNeighborCorrelationPlot
private

Definition at line 30 of file CSCStripDigiValidation.h.

MonitorElement* CSCStripDigiValidation::thePedestalPlot
private

Definition at line 28 of file CSCStripDigiValidation.h.

Referenced by bookHistograms(), and fillPedestalPlots().

float CSCStripDigiValidation::thePedestalSum
private

Definition at line 25 of file CSCStripDigiValidation.h.

Referenced by analyze(), and fillSignalPlots().

MonitorElement* CSCStripDigiValidation::thePedestalTimeCorrelationPlot
private

Definition at line 29 of file CSCStripDigiValidation.h.

MonitorElement* CSCStripDigiValidation::theRatio4to5Plot
private

Definition at line 32 of file CSCStripDigiValidation.h.

Referenced by bookHistograms(), and fillSignalPlots().

MonitorElement* CSCStripDigiValidation::theRatio6to5Plot
private

Definition at line 33 of file CSCStripDigiValidation.h.

Referenced by bookHistograms(), and fillSignalPlots().

MonitorElement* CSCStripDigiValidation::theResolutionPlots[10]
private

Definition at line 37 of file CSCStripDigiValidation.h.

Referenced by bookHistograms(), and plotResolution().