#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::ParameterSet &, edm::ConsumesCollector &&iC)

	~CSCSegmentValidation () 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 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)

Private Attributes
edm::InputTag	inputTag_

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 *	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
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 8 of file CSCSegmentValidation.h.

Member Typedef Documentation

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

private

Definition at line 26 of file CSCSegmentValidation.h.

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

private

Definition at line 33 of file CSCSegmentValidation.h.

Constructor & Destructor Documentation

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

Definition at line 5 of file CSCSegmentValidation.cc.

References edm::ParameterSet::getParameterSet(), inputTag_, TrackValidation_cff::pset, and segments_Token_.

     : CSCBaseValidation(ps), theLayerHitsPerChamber(), theChamberSegmentMap(), theShowerThreshold(10) {
   const auto &pset = ps.getParameterSet("cscSegment");
   inputTag_ = pset.getParameter<edm::InputTag>("inputTag");
   segments_Token_ = iC.consumes<CSCSegmentCollection>(inputTag_);
 }

CSCSegmentValidation::~CSCSegmentValidation ( )

override

Definition at line 12 of file CSCSegmentValidation.cc.

12 {}

Member Function Documentation

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

overridevirtual

Implements CSCBaseValidation.

Definition at line 53 of file CSCSegmentValidation.cc.

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

                                                                                      {
   // 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 (auto segmentItr = cscRecHits->begin(); segmentItr != cscRecHits->end(); segmentItr++) {
     ++nPerEvent;
     int detId = segmentItr->geographicalId().rawId();
     int chamberType = segmentItr->cscDetId().iChamberType();
 
     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 14 of file CSCSegmentValidation.cc.

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

                                                                   {
   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 = 1; i <= 10; ++i) {
     const std::string cn(CSCDetId::chamberName(i));
 
     const std::string t1("CSCSegmentRdPhiResolution_" + cn);
     const std::string t2("CSCSegmentRdPhiPull_" + cn);
     const std::string t3("CSCSegmentThetaResolution_" + cn);
     const std::string t5("CSCSegmentdXdZResolution_" + cn);
     const std::string t6("CSCSegmentdXdZPull_" + cn);
     const std::string t7("CSCSegmentdYdZResolution_" + cn);
     const std::string t8("CSCSegmentdYdZPull_" + cn);
 
     theRdPhiResolutionPlots[i - 1] = iBooker.book1D(t1, t1, 100, -0.4, 0.4);
     theRdPhiPullPlots[i - 1] = iBooker.book1D(t2, t2, 100, -5, 5);
     theThetaResolutionPlots[i - 1] = iBooker.book1D(t3, t3, 100, -1, 1);
     thedXdZResolutionPlots[i - 1] = iBooker.book1D(t5, t5, 100, -1, 1);
     thedXdZPullPlots[i - 1] = iBooker.book1D(t6, t6, 100, -5, 5);
     thedYdZResolutionPlots[i - 1] = iBooker.book1D(t7, t7, 100, -1, 1);
     thedYdZPullPlots[i - 1] = iBooker.book1D(t8, t8, 100, -5, 5);
   }
 }

void CSCSegmentValidation::fillEfficiencyPlots ( )

private

Definition at line 84 of file CSCSegmentValidation.cc.

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

Referenced by analyze().

                                                {
   // now plot efficiency by looping over all chambers with hits
   for (auto 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 = CSCDetId(chamberId).iChamberType();
     // 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 178 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 (auto 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 192 of file CSCSegmentValidation.cc.

References PSimHitMap::hits(), fireworks::p1, fireworks::p2, PSimHit::pabs(), mps_fire::result, and CSCBaseValidation::theSimHitMap.

Referenced by analyze().

                                                                {
   auto SimHitPabsLessThan = [](const PSimHit &p1, const PSimHit &p2) -> bool { return p1.pabs() < p2.pabs(); };
 
   const PSimHit *result = nullptr;
   int layerId = chamberId + 3;
   const auto &layerHits = theSimHitMap->hits(layerId);
 
   if (!layerHits.empty()) {
     // pick the hit with maximum energy
     auto hitItr = std::max_element(layerHits.begin(), layerHits.end(), SimHitPabsLessThan);
     result = &(*hitItr);
   }
   return result;
 }

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

private

Definition at line 148 of file CSCSegmentValidation.cc.

References dqm::impl::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 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);
 
   thedXdZResolutionPlots[chamberType - 1]->Fill(ddxdz);
   thedXdZPullPlots[chamberType - 1]->Fill(ddxdz / sigmadxdz);
   thedYdZResolutionPlots[chamberType - 1]->Fill(ddydz);
   thedYdZPullPlots[chamberType - 1]->Fill(ddydz / sigmadydz);
 }