CSCSegmentValidation Class Reference

#include <CSCSegmentValidation.h>

Inheritance diagram for CSCSegmentValidation:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
	CSCSegmentValidation (DQMStore *dbe, const edm::InputTag &inputTag, edm::ConsumesCollector &&iC)
virtual	~CSCSegmentValidation ()
Public Member Functions inherited from CSCBaseValidation
	CSCBaseValidation (DQMStore *dbe, 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
DQMStore *	dbe_
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 31 of file CSCSegmentValidation.h.

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

private

Definition at line 38 of file CSCSegmentValidation.h.

Constructor & Destructor Documentation

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

Definition at line 7 of file CSCSegmentValidation.cc.

References DQMStore::book1D(), CSCBaseValidation::dbe_, i, segments_Token_, DQMStore::setCurrentFolder(), thedXdZPullPlots, thedXdZResolutionPlots, thedYdZPullPlots, thedYdZResolutionPlots, theRdPhiPullPlots, theRdPhiResolutionPlots, theThetaPullPlots, and theThetaResolutionPlots.

: CSCBaseValidation(dbe, inputTag),
  theLayerHitsPerChamber(),
  theChamberSegmentMap(),
  theShowerThreshold(10),
  theNPerEventPlot( dbe_->book1D("CSCSegmentsPerEvent", "Number of CSC segments per event", 100, 0, 50) ),
  theNRecHitsPlot( dbe_->book1D("CSCRecHitsPerSegment", "Number of CSC rec hits per segment" , 8, 0, 7) ),
  theNPerChamberTypePlot( dbe_->book1D("CSCSegmentsPerChamberType", "Number of CSC segments per chamber type", 11, 0, 10) ),
  theTypePlot4HitsNoShower( dbe_->book1D("CSCSegments4HitsNoShower", "", 100, 0, 10) ),
  theTypePlot4HitsNoShowerSeg( dbe_->book1D("CSCSegments4HitsNoShowerSeg", "", 100, 0, 10) ),
  theTypePlot4HitsShower( dbe_->book1D("CSCSegments4HitsShower", "", 100, 0, 10) ),
  theTypePlot4HitsShowerSeg( dbe_->book1D("CSCSegments4HitsShowerSeg", "", 100, 0, 10) ),
  theTypePlot5HitsNoShower( dbe_->book1D("CSCSegments5HitsNoShower", "", 100, 0, 10) ),
  theTypePlot5HitsNoShowerSeg( dbe_->book1D("CSCSegments5HitsNoShowerSeg", "", 100, 0, 10) ),
  theTypePlot5HitsShower( dbe_->book1D("CSCSegments5HitsShower", "", 100, 0, 10) ),
  theTypePlot5HitsShowerSeg( dbe_->book1D("CSCSegments5HitsShowerSeg", "", 100, 0, 10) ),
  theTypePlot6HitsNoShower( dbe_->book1D("CSCSegments6HitsNoShower", "", 100, 0, 10) ),
  theTypePlot6HitsNoShowerSeg( dbe_->book1D("CSCSegments6HitsNoShowerSeg", "", 100, 0, 10) ),
  theTypePlot6HitsShower( dbe_->book1D("CSCSegments6HitsShower", "", 100, 0, 10) ),
  theTypePlot6HitsShowerSeg( dbe_->book1D("CSCSegments6HitsShowerSeg", "", 100, 0, 10) )
{
   segments_Token_ = iC.consumes<CSCSegmentCollection>(inputTag);

   dbe_->setCurrentFolder("CSCRecHitsV/CSCRecHitTask");
   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] = dbe_->book1D(title1, title1, 100, -0.4, 0.4);
    theRdPhiPullPlots[i] = dbe_->book1D(title2, title2, 100, -5, 5);
    theThetaResolutionPlots[i] = dbe_->book1D(title3, title3, 100, -1, 1);
    theThetaPullPlots[i] = dbe_->book1D(title4, title4, 100, -5, 5);
    thedXdZResolutionPlots[i] = dbe_->book1D(title5, title5, 100, -1, 1);
    thedXdZPullPlots[i] = dbe_->book1D(title6, title6, 100, -5, 5);
    thedYdZResolutionPlots[i] = dbe_->book1D(title7, title7, 100, -1, 1);
    thedYdZPullPlots[i] = dbe_->book1D(title8, title8, 100, -5, 5);

  }
}

virtual CSCSegmentValidation::~CSCSegmentValidation ( )

inlinevirtual

Definition at line 17 of file CSCSegmentValidation.h.

{}

Member Function Documentation

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

virtual

Implements CSCBaseValidation.

Definition at line 57 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 != 0) 
    {
      const CSCLayer * layer = findLayer(hit->detUnitId());
      plotResolution(*hit, *segmentItr, layer, chamberType);
    }  
  }

  theNPerEventPlot->Fill(nPerEvent);

  fillLayerHitsPerChamber();
  fillEfficiencyPlots();
}

void CSCSegmentValidation::fillEfficiencyPlots ( )

private

Definition at line 93 of file CSCSegmentValidation.cc.

References MonitorElement::Fill(), hasSegment(), python.multivaluedict::sort(), theLayerHitsPerChamber, theShowerThreshold, theTypePlot4HitsNoShower, theTypePlot4HitsNoShowerSeg, theTypePlot4HitsShower, theTypePlot4HitsShowerSeg, theTypePlot5HitsNoShower, theTypePlot5HitsNoShowerSeg, theTypePlot5HitsShower, theTypePlot5HitsShowerSeg, theTypePlot6HitsNoShower, theTypePlot6HitsNoShowerSeg, theTypePlot6HitsShower, theTypePlot6HitsShowerSeg, 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 201 of file CSCSegmentValidation.cc.

References CSCDetId::chamberId(), PSimHitMap::detsWithHits(), PSimHitMap::hits(), i, 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 156 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 229 of file CSCSegmentValidation.cc.

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

Referenced by analyze().

{
  const PSimHit * result = 0;
  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 169 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 162 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());
}

Member Data Documentation

edm::EDGetTokenT<CSCSegmentCollection> CSCSegmentValidation::segments_Token_

private

Definition at line 27 of file CSCSegmentValidation.h.

Referenced by analyze(), and CSCSegmentValidation().

ChamberSegmentMap CSCSegmentValidation::theChamberSegmentMap

private

Definition at line 39 of file CSCSegmentValidation.h.

Referenced by analyze(), and hasSegment().

MonitorElement* CSCSegmentValidation::thedXdZPullPlots[10]

private

Definition at line 51 of file CSCSegmentValidation.h.

Referenced by CSCSegmentValidation(), and plotResolution().

MonitorElement* CSCSegmentValidation::thedXdZResolutionPlots[10]

private

Definition at line 50 of file CSCSegmentValidation.h.

Referenced by CSCSegmentValidation(), and plotResolution().

MonitorElement* CSCSegmentValidation::thedYdZPullPlots[10]

private

Definition at line 53 of file CSCSegmentValidation.h.

Referenced by CSCSegmentValidation(), and plotResolution().

MonitorElement* CSCSegmentValidation::thedYdZResolutionPlots[10]

private

Definition at line 52 of file CSCSegmentValidation.h.

Referenced by CSCSegmentValidation(), and plotResolution().

ChamberHitMap CSCSegmentValidation::theLayerHitsPerChamber

private

Definition at line 32 of file CSCSegmentValidation.h.

Referenced by fillEfficiencyPlots(), and fillLayerHitsPerChamber().

MonitorElement* CSCSegmentValidation::theNPerChamberTypePlot

private

Definition at line 45 of file CSCSegmentValidation.h.

Referenced by analyze().

MonitorElement* CSCSegmentValidation::theNPerEventPlot

private

Definition at line 43 of file CSCSegmentValidation.h.

Referenced by analyze().

MonitorElement* CSCSegmentValidation::theNRecHitsPlot

private

Definition at line 44 of file CSCSegmentValidation.h.

Referenced by analyze().

MonitorElement* CSCSegmentValidation::theRdPhiPullPlots[10]

private

Definition at line 47 of file CSCSegmentValidation.h.

Referenced by CSCSegmentValidation(), and plotResolution().

MonitorElement* CSCSegmentValidation::theRdPhiResolutionPlots[10]

private

Definition at line 46 of file CSCSegmentValidation.h.

Referenced by CSCSegmentValidation(), and plotResolution().

int CSCSegmentValidation::theShowerThreshold

private

Definition at line 41 of file CSCSegmentValidation.h.

Referenced by fillEfficiencyPlots().

MonitorElement* CSCSegmentValidation::theThetaPullPlots[10]

private

Definition at line 49 of file CSCSegmentValidation.h.

Referenced by CSCSegmentValidation().

MonitorElement* CSCSegmentValidation::theThetaResolutionPlots[10]

private

Definition at line 48 of file CSCSegmentValidation.h.

Referenced by CSCSegmentValidation(), and plotResolution().

MonitorElement* CSCSegmentValidation::theTypePlot4HitsNoShower

private

Definition at line 57 of file CSCSegmentValidation.h.

Referenced by fillEfficiencyPlots().

MonitorElement* CSCSegmentValidation::theTypePlot4HitsNoShowerSeg

private

Definition at line 58 of file CSCSegmentValidation.h.

Referenced by fillEfficiencyPlots().

MonitorElement* CSCSegmentValidation::theTypePlot4HitsShower

private

Definition at line 59 of file CSCSegmentValidation.h.

Referenced by fillEfficiencyPlots().

MonitorElement* CSCSegmentValidation::theTypePlot4HitsShowerSeg

private

Definition at line 60 of file CSCSegmentValidation.h.

Referenced by fillEfficiencyPlots().

MonitorElement* CSCSegmentValidation::theTypePlot5HitsNoShower

private

Definition at line 61 of file CSCSegmentValidation.h.

Referenced by fillEfficiencyPlots().

MonitorElement* CSCSegmentValidation::theTypePlot5HitsNoShowerSeg

private

Definition at line 62 of file CSCSegmentValidation.h.

Referenced by fillEfficiencyPlots().

MonitorElement* CSCSegmentValidation::theTypePlot5HitsShower

private

Definition at line 63 of file CSCSegmentValidation.h.

Referenced by fillEfficiencyPlots().

MonitorElement* CSCSegmentValidation::theTypePlot5HitsShowerSeg

private

Definition at line 64 of file CSCSegmentValidation.h.

Referenced by fillEfficiencyPlots().

MonitorElement* CSCSegmentValidation::theTypePlot6HitsNoShower

private

Definition at line 65 of file CSCSegmentValidation.h.

Referenced by fillEfficiencyPlots().

MonitorElement* CSCSegmentValidation::theTypePlot6HitsNoShowerSeg

private

Definition at line 66 of file CSCSegmentValidation.h.

Referenced by fillEfficiencyPlots().

MonitorElement* CSCSegmentValidation::theTypePlot6HitsShower

private

Definition at line 67 of file CSCSegmentValidation.h.

Referenced by fillEfficiencyPlots().

MonitorElement* CSCSegmentValidation::theTypePlot6HitsShowerSeg

private

Definition at line 68 of file CSCSegmentValidation.h.

Referenced by fillEfficiencyPlots().