#include <CSCRecHit2DValidation.h>

Inheritance diagram for CSCRecHit2DValidation:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
	CSCRecHit2DValidation (DQMStore *dbe, const edm::InputTag &inputTag)
virtual	~CSCRecHit2DValidation ()
Private Member Functions
void	plotResolution (const PSimHit &simHit, const CSCRecHit2D &recHit, const CSCLayer *layer, int chamberType)
Private Attributes
MonitorElement *	theNPerEventPlot
MonitorElement *	thePullPlots [10]
MonitorElement *	theRecHitPosInStrip [10]
MonitorElement *	theResolutionPlots [10]
MonitorElement *	theScatterPlots [10]
MonitorElement *	theSimHitPosInStrip [10]
MonitorElement *	theSimHitScatterPlots [10]
MonitorElement *	theTPeaks [10]
MonitorElement *	theYPullPlots [10]
MonitorElement *	theYResolutionPlots [10]

Detailed Description

Definition at line 10 of file CSCRecHit2DValidation.h.

Constructor & Destructor Documentation

CSCRecHit2DValidation::CSCRecHit2DValidation	(	DQMStore *	dbe,
		const edm::InputTag &	inputTag
	)

Definition at line 8 of file CSCRecHit2DValidation.cc.

References DQMStore::book1D(), DQMStore::book2D(), CSCBaseValidation::dbe_, i, DQMStore::setCurrentFolder(), thePullPlots, theRecHitPosInStrip, theResolutionPlots, theScatterPlots, theSimHitPosInStrip, theSimHitScatterPlots, theTPeaks, theYPullPlots, and theYResolutionPlots.

: CSCBaseValidation(dbe, inputTag),
  theNPerEventPlot( dbe_->book1D("CSCRecHitsPerEvent", "Number of CSC Rec Hits per event", 100, 0, 500) )
{
   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], title9[200];
    sprintf(title1, "CSCRecHitResolution%d", i+1);
    sprintf(title2, "CSCRecHitPull%d", i+1);
    sprintf(title3, "CSCRecHitYResolution%d", i+1);
    sprintf(title4, "CSCRecHitYPull%d", i+1);
    sprintf(title5, "CSCRecHitPosInStrip%d", i+1);
    sprintf(title6, "CSCSimHitPosInStrip%d", i+1);
    sprintf(title7, "CSCRecHit%d", i+1);
    sprintf(title8, "CSCSimHit%d", i+1);
    sprintf(title9, "CSCTPeak%d", i+1);


    theResolutionPlots[i] = dbe_->book1D(title1, title1, 100, -0.2, 0.2);
    thePullPlots[i] = dbe_->book1D(title2, title2, 100, -3, 3);
    theYResolutionPlots[i] = dbe_->book1D(title3, title3, 100, -5, 5);
    theYPullPlots[i] = dbe_->book1D(title4, title4, 100, -3, 3);
    theRecHitPosInStrip[i] = dbe_->book1D(title5, title5, 100, -2, 2);
    theSimHitPosInStrip[i] = dbe_->book1D(title6, title6, 100, -2, 2);
    theScatterPlots[i] = dbe->book2D(title7, title7, 200, -20, 20, 200, -250, 250);
    theSimHitScatterPlots[i] = dbe->book2D(title8, title8, 200, -20, 20, 200, -250, 250);
    theTPeaks[i] =  dbe->book1D(title9, title9, 200, 0, 400);
  }

}

CSCRecHit2DValidation::~CSCRecHit2DValidation ( ) [virtual]

Definition at line 41 of file CSCRecHit2DValidation.cc.

References MonitorElement::getMean(), MonitorElement::getName(), MonitorElement::getRMS(), i, theResolutionPlots, and theTPeaks.

{
  for(int i = 0; i < 10; ++i)
  {
    edm::LogInfo("CSCRecHitValidation") << "Resolution of " << theResolutionPlots[i]->getName() << " is " << theResolutionPlots[i]->getRMS();
    edm::LogInfo("CSCRecHitValidation") << "Peak Time is " << theTPeaks[i]->getMean();
  }
}

Member Function Documentation

void CSCRecHit2DValidation::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	eventSetup
	)		`[virtual]`

Implements CSCBaseValidation.

Definition at line 51 of file CSCRecHit2DValidation.cc.

References CSCLayer::chamber(), CSCChamberSpecs::chamberType(), Geom::Phi< T >::degrees(), PSimHitMap::detsWithHits(), MonitorElement::Fill(), CSCBaseValidation::findLayer(), edm::Event::getByLabel(), PSimHitMap::hits(), i, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), plotResolution(), edm::Handle< T >::product(), trackerHits::simHits, CSCChamber::specs(), CSCBaseValidation::theInputTag, theNPerEventPlot, theScatterPlots, CSCBaseValidation::theSimHitMap, theSimHitScatterPlots, theTPeaks, and GeomDet::toGlobal().

Referenced by CSCRecHitValidation::analyze().

{
  // get the collection of CSCRecHrecHitItrD
  edm::Handle<CSCRecHit2DCollection> hRecHits;
  e.getByLabel(theInputTag, hRecHits);
  const CSCRecHit2DCollection * cscRecHits = hRecHits.product();

  unsigned nPerEvent = 0;

  for(CSCRecHit2DCollection::const_iterator recHitItr = cscRecHits->begin(); 
      recHitItr != cscRecHits->end(); recHitItr++) 
  {
    ++nPerEvent;
    int detId = (*recHitItr).cscDetId().rawId();
    edm::PSimHitContainer simHits = theSimHitMap->hits(detId);
    const CSCLayer * layer = findLayer(detId);
    int chamberType = layer->chamber()->specs()->chamberType();
    theTPeaks[chamberType-1]->Fill(recHitItr->tpeak());
    if(simHits.size() == 1)
    {
      plotResolution(simHits[0], *recHitItr, layer, chamberType);
    }
    float localX = recHitItr->localPosition().x();
    float localY = recHitItr->localPosition().y();
    //theYPlots[chamberType-1]->Fill(localY);
    // find a local phi
    float globalR = layer->toGlobal(LocalPoint(0.,0.,0.)).perp();
    GlobalPoint axisThruChamber(globalR+localY, localX, 0.);
    float localPhi = axisThruChamber.phi().degrees();
    //thePhiPlots[chamberType-1]->Fill(axisThruChamber.phi().degrees());
    theScatterPlots[chamberType-1]->Fill( localPhi, localY);
  }    
  theNPerEventPlot->Fill(nPerEvent);
return;
  // fill sim hits
  std::vector<int> layersWithSimHits = theSimHitMap->detsWithHits();
  for(unsigned i = 0; i < layersWithSimHits.size(); ++i)
   {
    edm::PSimHitContainer simHits = theSimHitMap->hits(layersWithSimHits[i]);
    for(edm::PSimHitContainer::const_iterator hitItr = simHits.begin(); hitItr != simHits.end(); ++hitItr)
    {
    const CSCLayer * layer = findLayer(layersWithSimHits[i]);
    int chamberType = layer->chamber()->specs()->chamberType();
      float localX = hitItr->localPosition().x();
      float localY = hitItr->localPosition().y();
      //theYPlots[chamberType-1]->Fill(localY);
      // find a local phi
      float globalR = layer->toGlobal(LocalPoint(0.,0.,0.)).perp();
      GlobalPoint axisThruChamber(globalR+localY, localX, 0.);
      float localPhi = axisThruChamber.phi().degrees();
      //thePhiPlots[chamberType-1]->Fill(axisThruChamber.phi().degrees());
      theSimHitScatterPlots[chamberType-1]->Fill( localPhi, localY);
    }
  }

}

void CSCRecHit2DValidation::plotResolution	(	const PSimHit &	simHit,
		const CSCRecHit2D &	recHit,
		const CSCLayer *	layer,
		int	chamberType
	)		`[private]`

Definition at line 109 of file CSCRecHit2DValidation.cc.

References MonitorElement::Fill(), CSCLayer::geometry(), PSimHit::localPosition(), CSCRecHit2D::localPosition(), CSCRecHit2D::localPositionError(), PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), mathSSE::sqrt(), CSCLayerGeometry::strip(), thePullPlots, theRecHitPosInStrip, theResolutionPlots, theSimHitPosInStrip, theYPullPlots, theYResolutionPlots, GeomDet::toGlobal(), LocalError::xx(), PV3DBase< T, PVType, FrameType >::y(), and LocalError::yy().

Referenced by analyze().

{
  GlobalPoint simHitPos = layer->toGlobal(simHit.localPosition());
  GlobalPoint recHitPos = layer->toGlobal(recHit.localPosition());

  double dphi = recHitPos.phi() - simHitPos.phi();
  double rdphi = recHitPos.perp() * dphi;
  theResolutionPlots[chamberType-1]->Fill( rdphi );
  thePullPlots[chamberType-1]->Fill( rdphi/ sqrt(recHit.localPositionError().xx()) );
  double dy = recHit.localPosition().y() - simHit.localPosition().y();
  theYResolutionPlots[chamberType-1]->Fill( dy );
  theYPullPlots[chamberType-1]->Fill( dy/ sqrt(recHit.localPositionError().yy()) );

  const CSCLayerGeometry * layerGeometry = layer->geometry();
  float recStrip = layerGeometry->strip(recHit.localPosition());
  float simStrip = layerGeometry->strip(simHit.localPosition());
  theRecHitPosInStrip[chamberType-1]->Fill( recStrip - int(recStrip) );
  theSimHitPosInStrip[chamberType-1]->Fill( simStrip - int(simStrip) );
}