#include <CSCSegmentValidation.h>

Inheritance diagram for CSCSegmentValidation:

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

void	bookHistograms (DQMStore::IBooker &)

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

	~CSCSegmentValidation () 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 Types
typedef std::map< int, std::vector< int > >	ChamberHitMap

typedef std::map< int, std::vector< CSCSegment > >	ChamberSegmentMap

Private Member Functions
void	fillEfficiencyPlots ()

void	fillLayerHitsPerChamber ()

bool	hasSegment (int chamberId) const

const PSimHit *	keyHit (int chamberId) const

void	plotResolution (const PSimHit &simHit, const CSCSegment &recHit, const CSCLayer *layer, int chamberType)

Static Private Member Functions
static int	whatChamberType (int detId)

Private Attributes
edm::EDGetTokenT< CSCSegmentCollection >	segments_Token_

ChamberSegmentMap	theChamberSegmentMap

MonitorElement *	thedXdZPullPlots [10]

MonitorElement *	thedXdZResolutionPlots [10]

MonitorElement *	thedYdZPullPlots [10]

MonitorElement *	thedYdZResolutionPlots [10]

ChamberHitMap	theLayerHitsPerChamber

MonitorElement *	theNPerChamberTypePlot

MonitorElement *	theNPerEventPlot

MonitorElement *	theNRecHitsPlot

MonitorElement *	theRdPhiPullPlots [10]

MonitorElement *	theRdPhiResolutionPlots [10]

int	theShowerThreshold

MonitorElement *	theThetaPullPlots [10]

MonitorElement *	theThetaResolutionPlots [10]

MonitorElement *	theTypePlot4HitsNoShower

MonitorElement *	theTypePlot4HitsNoShowerSeg

MonitorElement *	theTypePlot4HitsShower

MonitorElement *	theTypePlot4HitsShowerSeg

MonitorElement *	theTypePlot5HitsNoShower

MonitorElement *	theTypePlot5HitsNoShowerSeg

MonitorElement *	theTypePlot5HitsShower

MonitorElement *	theTypePlot5HitsShowerSeg

MonitorElement *	theTypePlot6HitsNoShower

MonitorElement *	theTypePlot6HitsNoShowerSeg

MonitorElement *	theTypePlot6HitsShower

MonitorElement *	theTypePlot6HitsShowerSeg

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 CSCSegmentValidation.h.

Member Typedef Documentation

typedef std::map<int, std::vector<int> > CSCSegmentValidation::ChamberHitMap

private

Definition at line 29 of file CSCSegmentValidation.h.

typedef std::map<int, std::vector<CSCSegment> > CSCSegmentValidation::ChamberSegmentMap

private

Definition at line 36 of file CSCSegmentValidation.h.

Constructor & Destructor Documentation

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

Definition at line 5 of file CSCSegmentValidation.cc.

References HiClusterCompatibility_cfi::inputTag, and segments_Token_.

     : CSCBaseValidation(inputTag), theLayerHitsPerChamber(), theChamberSegmentMap(), theShowerThreshold(10) {
   segments_Token_ = iC.consumes<CSCSegmentCollection>(inputTag);
 }

CSCSegmentValidation::~CSCSegmentValidation ( )

override

Definition at line 10 of file CSCSegmentValidation.cc.

10 {}

Member Function Documentation

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

overridevirtual

Implements CSCBaseValidation.

Definition at line 51 of file CSCSegmentValidation.cc.

References PSimHit::detUnitId(), MonitorElement::Fill(), fillEfficiencyPlots(), fillLayerHitsPerChamber(), CSCBaseValidation::findLayer(), edm::Event::getByToken(), keyHit(), plotResolution(), edm::Handle< T >::product(), segments_Token_, theChamberSegmentMap, theNPerChamberTypePlot, theNPerEventPlot, theNRecHitsPlot, and whatChamberType().

Referenced by CSCRecHitValidation::analyze().

                                                                                      {
   // get the collection of CSCRecHsegmentItrD
   edm::Handle<CSCSegmentCollection> hRecHits;
   e.getByToken(segments_Token_, hRecHits);
   const CSCSegmentCollection *cscRecHits = hRecHits.product();
 
   theChamberSegmentMap.clear();
   unsigned nPerEvent = 0;
   for (CSCSegmentCollection::const_iterator segmentItr = cscRecHits->begin(); segmentItr != cscRecHits->end();
        segmentItr++) {
     ++nPerEvent;
     int detId = segmentItr->geographicalId().rawId();
     int chamberType = whatChamberType(detId);
 
     theNRecHitsPlot->Fill(segmentItr->nRecHits());
     theNPerChamberTypePlot->Fill(chamberType);
     theChamberSegmentMap[detId].push_back(*segmentItr);
 
     // do the resolution plots
     const PSimHit *hit = keyHit(detId);
     if (hit != nullptr) {
       const CSCLayer *layer = findLayer(hit->detUnitId());
       plotResolution(*hit, *segmentItr, layer, chamberType);
     }
   }
 
   theNPerEventPlot->Fill(nPerEvent);
 
   fillLayerHitsPerChamber();
   fillEfficiencyPlots();
 }

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

Definition at line 12 of file CSCSegmentValidation.cc.

References DQMStore::IBooker::book1D(), mps_fire::i, thedXdZPullPlots, thedXdZResolutionPlots, thedYdZPullPlots, thedYdZResolutionPlots, theNPerChamberTypePlot, theNPerEventPlot, theNRecHitsPlot, theRdPhiPullPlots, theRdPhiResolutionPlots, theThetaPullPlots, theThetaResolutionPlots, theTypePlot4HitsNoShower, theTypePlot4HitsNoShowerSeg, theTypePlot4HitsShower, theTypePlot4HitsShowerSeg, theTypePlot5HitsNoShower, theTypePlot5HitsNoShowerSeg, theTypePlot5HitsShower, theTypePlot5HitsShowerSeg, theTypePlot6HitsNoShower, theTypePlot6HitsNoShowerSeg, theTypePlot6HitsShower, and theTypePlot6HitsShowerSeg.

Referenced by CSCRecHitValidation::bookHistograms().

                                                                   {
   theNPerEventPlot = iBooker.book1D("CSCSegmentsPerEvent", "Number of CSC segments per event", 100, 0, 50);
   theNRecHitsPlot = iBooker.book1D("CSCRecHitsPerSegment", "Number of CSC rec hits per segment", 8, 0, 7);
   theNPerChamberTypePlot =
       iBooker.book1D("CSCSegmentsPerChamberType", "Number of CSC segments per chamber type", 11, 0, 10);
   theTypePlot4HitsNoShower = iBooker.book1D("CSCSegments4HitsNoShower", "", 100, 0, 10);
   theTypePlot4HitsNoShowerSeg = iBooker.book1D("CSCSegments4HitsNoShowerSeg", "", 100, 0, 10);
   theTypePlot4HitsShower = iBooker.book1D("CSCSegments4HitsShower", "", 100, 0, 10);
   theTypePlot4HitsShowerSeg = iBooker.book1D("CSCSegments4HitsShowerSeg", "", 100, 0, 10);
   theTypePlot5HitsNoShower = iBooker.book1D("CSCSegments5HitsNoShower", "", 100, 0, 10);
   theTypePlot5HitsNoShowerSeg = iBooker.book1D("CSCSegments5HitsNoShowerSeg", "", 100, 0, 10);
   theTypePlot5HitsShower = iBooker.book1D("CSCSegments5HitsShower", "", 100, 0, 10);
   theTypePlot5HitsShowerSeg = iBooker.book1D("CSCSegments5HitsShowerSeg", "", 100, 0, 10);
   theTypePlot6HitsNoShower = iBooker.book1D("CSCSegments6HitsNoShower", "", 100, 0, 10);
   theTypePlot6HitsNoShowerSeg = iBooker.book1D("CSCSegments6HitsNoShowerSeg", "", 100, 0, 10);
   theTypePlot6HitsShower = iBooker.book1D("CSCSegments6HitsShower", "", 100, 0, 10);
   theTypePlot6HitsShowerSeg = iBooker.book1D("CSCSegments6HitsShowerSeg", "", 100, 0, 10);
   for (int i = 0; i < 10; ++i) {
     char title1[200], title2[200], title3[200], title4[200], title5[200], title6[200], title7[200], title8[200];
     sprintf(title1, "CSCSegmentRdPhiResolution%d", i + 1);
     sprintf(title2, "CSCSegmentRdPhiPull%d", i + 1);
     sprintf(title3, "CSCSegmentThetaResolution%d", i + 1);
     sprintf(title4, "CSCSegmentThetaPull%d", i + 1);
     sprintf(title5, "CSCSegmentdXdZResolution%d", i + 1);
     sprintf(title6, "CSCSegmentdXdZPull%d", i + 1);
     sprintf(title7, "CSCSegmentdYdZResolution%d", i + 1);
     sprintf(title8, "CSCSegmentdYdZPull%d", i + 1);
 
     theRdPhiResolutionPlots[i] = iBooker.book1D(title1, title1, 100, -0.4, 0.4);
     theRdPhiPullPlots[i] = iBooker.book1D(title2, title2, 100, -5, 5);
     theThetaResolutionPlots[i] = iBooker.book1D(title3, title3, 100, -1, 1);
     theThetaPullPlots[i] = iBooker.book1D(title4, title4, 100, -5, 5);
     thedXdZResolutionPlots[i] = iBooker.book1D(title5, title5, 100, -1, 1);
     thedXdZPullPlots[i] = iBooker.book1D(title6, title6, 100, -5, 5);
     thedYdZResolutionPlots[i] = iBooker.book1D(title7, title7, 100, -1, 1);
     thedYdZPullPlots[i] = iBooker.book1D(title8, title8, 100, -5, 5);
   }
 }

void CSCSegmentValidation::fillEfficiencyPlots ( )

private

Definition at line 83 of file CSCSegmentValidation.cc.

References MonitorElement::Fill(), hasSegment(), theLayerHitsPerChamber, theShowerThreshold, theTypePlot4HitsNoShower, theTypePlot4HitsNoShowerSeg, theTypePlot4HitsShower, theTypePlot4HitsShowerSeg, theTypePlot5HitsNoShower, theTypePlot5HitsNoShowerSeg, theTypePlot5HitsShower, theTypePlot5HitsShowerSeg, theTypePlot6HitsNoShower, theTypePlot6HitsNoShowerSeg, theTypePlot6HitsShower, theTypePlot6HitsShowerSeg, tier0::unique(), findQualityFiles::v, and whatChamberType().

Referenced by analyze().

                                                {
   // now plot efficiency by looping over all chambers with hits
   for (ChamberHitMap::const_iterator mapItr = theLayerHitsPerChamber.begin(), mapEnd = theLayerHitsPerChamber.end();
        mapItr != mapEnd;
        ++mapItr) {
     int chamberId = mapItr->first;
     int nHitsInChamber = mapItr->second.size();
     bool isShower = (nHitsInChamber > theShowerThreshold);
     bool hasSeg = hasSegment(chamberId);
     int chamberType = whatChamberType(chamberId);
     // find how many layers were hit in this chamber
     std::vector<int> v = mapItr->second;
     std::sort(v.begin(), v.end());
     // maybe can just count
     v.erase(std::unique(v.begin(), v.end()), v.end());
     int nLayersHit = v.size();
 
     if (nLayersHit == 4) {
       if (isShower)
         theTypePlot4HitsShower->Fill(chamberType);
       else
         theTypePlot4HitsNoShower->Fill(chamberType);
 
       if (hasSeg) {
         if (isShower)
           theTypePlot4HitsShowerSeg->Fill(chamberType);
         else
           theTypePlot4HitsNoShowerSeg->Fill(chamberType);
       }
     }
 
     if (nLayersHit == 5) {
       if (isShower)
         theTypePlot5HitsShower->Fill(chamberType);
       else
         theTypePlot5HitsNoShower->Fill(chamberType);
 
       if (hasSeg) {
         if (isShower)
           theTypePlot5HitsShowerSeg->Fill(chamberType);
         else
           theTypePlot5HitsNoShowerSeg->Fill(chamberType);
       }
     }
 
     if (nLayersHit == 6) {
       if (isShower)
         theTypePlot6HitsShower->Fill(chamberType);
       else
         theTypePlot6HitsNoShower->Fill(chamberType);
 
       if (hasSeg) {
         if (isShower)
           theTypePlot6HitsShowerSeg->Fill(chamberType);
         else
           theTypePlot6HitsNoShowerSeg->Fill(chamberType);
       }
     }
   }
 }

void CSCSegmentValidation::fillLayerHitsPerChamber ( )

private

Definition at line 185 of file CSCSegmentValidation.cc.

References CSCDetId::chamberId(), PSimHitMap::detsWithHits(), PSimHitMap::hits(), mps_fire::i, nhits, DetId::rawId(), theLayerHitsPerChamber, and CSCBaseValidation::theSimHitMap.

Referenced by analyze().

                                                    {
   theLayerHitsPerChamber.clear();
   std::vector<int> layersHit = theSimHitMap->detsWithHits();
   for (std::vector<int>::const_iterator layerItr = layersHit.begin(), layersHitEnd = layersHit.end();
        layerItr != layersHitEnd;
        ++layerItr) {
     CSCDetId layerId(*layerItr);
     CSCDetId chamberId = layerId.chamberId();
     int nhits = theSimHitMap->hits(*layerItr).size();
     // multiple entries, so we can see showers
     for (int i = 0; i < nhits; ++i) {
       theLayerHitsPerChamber[chamberId.rawId()].push_back(*layerItr);
     }
   }
 }

bool CSCSegmentValidation::hasSegment ( int chamberId ) const

private

Definition at line 144 of file CSCSegmentValidation.cc.

References theChamberSegmentMap.

Referenced by fillEfficiencyPlots().

                                                          {
   return (theChamberSegmentMap.find(chamberId) != theChamberSegmentMap.end());
 }

const PSimHit * CSCSegmentValidation::keyHit ( int chamberId ) const

private

Definition at line 205 of file CSCSegmentValidation.cc.

References PSimHitMap::hits(), mps_fire::result, CSCSegmentValidationUtils::SimHitPabsLessThan(), and CSCBaseValidation::theSimHitMap.

Referenced by analyze().

                                                                {
   const PSimHit *result = nullptr;
   int layerId = chamberId + 3;
   const edm::PSimHitContainer &layerHits = theSimHitMap->hits(layerId);
 
   if (!layerHits.empty()) {
     // pick the hit with maximum energy
     edm::PSimHitContainer::const_iterator hitItr =
         std::max_element(layerHits.begin(), layerHits.end(), CSCSegmentValidationUtils::SimHitPabsLessThan);
     result = &(*hitItr);
   }
   return result;
 }

void CSCSegmentValidation::plotResolution	(	const PSimHit &	simHit,
		const CSCSegment &	recHit,
		const CSCLayer *	layer,
		int	chamberType
	)

private

Definition at line 153 of file CSCSegmentValidation.cc.

References MonitorElement::Fill(), CSCSegment::localDirection(), PSimHit::localDirection(), CSCSegment::localDirectionError(), CSCSegment::localPosition(), PSimHit::localPosition(), CSCSegment::localPositionError(), PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), mathSSE::sqrt(), thedXdZPullPlots, thedXdZResolutionPlots, thedYdZPullPlots, thedYdZResolutionPlots, theRdPhiPullPlots, theRdPhiResolutionPlots, PV3DBase< T, PVType, FrameType >::theta(), theThetaResolutionPlots, GeomDet::toGlobal(), PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), PV3DBase< T, PVType, FrameType >::y(), LocalError::yy(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by analyze().

                                                            {
   GlobalPoint simHitPos = layer->toGlobal(simHit.localPosition());
   GlobalPoint segmentPos = layer->toGlobal(segment.localPosition());
   LocalVector simHitDir = simHit.localDirection();
   LocalVector segmentDir = segment.localDirection();
 
   double dphi = segmentPos.phi() - simHitPos.phi();
   double rdphi = segmentPos.perp() * dphi;
   double dtheta = segmentPos.theta() - simHitPos.theta();
 
   double sigmax = sqrt(segment.localPositionError().xx());
   // double sigmay = sqrt(segment.localPositionError().yy());
 
   double ddxdz = segmentDir.x() / segmentDir.z() - simHitDir.x() / simHitDir.z();
   double ddydz = segmentDir.y() / segmentDir.z() - simHitDir.y() / simHitDir.z();
   double sigmadxdz = sqrt(segment.localDirectionError().xx());
   double sigmadydz = sqrt(segment.localDirectionError().yy());
 
   theRdPhiResolutionPlots[chamberType - 1]->Fill(rdphi);
   theRdPhiPullPlots[chamberType - 1]->Fill(rdphi / sigmax);
   theThetaResolutionPlots[chamberType - 1]->Fill(dtheta);
   // theThetaPullPlots[chamberType-1]->Fill( dy/sigmay );
 
   thedXdZResolutionPlots[chamberType - 1]->Fill(ddxdz);
   thedXdZPullPlots[chamberType - 1]->Fill(ddxdz / sigmadxdz);
   thedYdZResolutionPlots[chamberType - 1]->Fill(ddydz);
   thedYdZPullPlots[chamberType - 1]->Fill(ddydz / sigmadydz);
 }

int CSCSegmentValidation::whatChamberType ( int detId )

staticprivate

Definition at line 148 of file CSCSegmentValidation.cc.

References CSCDetId::ring(), CSCDetId::station(), and CSCChamberSpecs::whatChamberType().

Referenced by analyze(), and fillEfficiencyPlots().

                                                    {
   CSCDetId cscDetId(detId);
   return CSCChamberSpecs::whatChamberType(cscDetId.station(), cscDetId.ring());
 }