|
|
|
#include <CSCSegmentValidation.h>
Definition at line 10 of file CSCSegmentValidation.h.
typedef std::map<int, std::vector<int> > CSCSegmentValidation::ChamberHitMap [private] |
Definition at line 27 of file CSCSegmentValidation.h.
typedef std::map<int, std::vector<CSCSegment> > CSCSegmentValidation::ChamberSegmentMap [private] |
Definition at line 34 of file CSCSegmentValidation.h.
| CSCSegmentValidation::CSCSegmentValidation | ( | DQMStore * | dbe, |
| const edm::InputTag & | inputTag | ||
| ) |
Definition at line 8 of file CSCSegmentValidation.cc.
References DQMStore::book1D(), CSCBaseValidation::dbe_, i, 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) ) { 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 | ( | ) | [inline, virtual] |
Definition at line 15 of file CSCSegmentValidation.h.
{}
| void CSCSegmentValidation::analyze | ( | const edm::Event & | e, |
| const edm::EventSetup & | eventSetup | ||
| ) | [virtual] |
Implements CSCBaseValidation.
Definition at line 56 of file CSCSegmentValidation.cc.
References PSimHit::detUnitId(), MonitorElement::Fill(), fillEfficiencyPlots(), fillLayerHitsPerChamber(), CSCBaseValidation::findLayer(), edm::Event::getByLabel(), keyHit(), plotResolution(), edm::Handle< T >::product(), theChamberSegmentMap, CSCBaseValidation::theInputTag, theNPerChamberTypePlot, theNPerEventPlot, theNRecHitsPlot, and whatChamberType().
Referenced by CSCRecHitValidation::analyze().
{
// get the collection of CSCRecHsegmentItrD
edm::Handle<CSCSegmentCollection> hRecHits;
e.getByLabel(theInputTag, 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 92 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, 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 200 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 155 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 228 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 168 of file CSCSegmentValidation.cc.
References MonitorElement::Fill(), PSimHit::localDirection(), CSCSegment::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 | ) | [static, private] |
Definition at line 161 of file CSCSegmentValidation.cc.
References CSCDetId::ring(), and CSCDetId::station().
Referenced by analyze(), and fillEfficiencyPlots().
{
CSCDetId cscDetId(detId);
return CSCChamberSpecs::whatChamberType(cscDetId.station(), cscDetId.ring());
}
Definition at line 35 of file CSCSegmentValidation.h.
Referenced by analyze(), and hasSegment().
MonitorElement* CSCSegmentValidation::thedXdZPullPlots[10] [private] |
Definition at line 47 of file CSCSegmentValidation.h.
Referenced by CSCSegmentValidation(), and plotResolution().
MonitorElement* CSCSegmentValidation::thedXdZResolutionPlots[10] [private] |
Definition at line 46 of file CSCSegmentValidation.h.
Referenced by CSCSegmentValidation(), and plotResolution().
MonitorElement* CSCSegmentValidation::thedYdZPullPlots[10] [private] |
Definition at line 49 of file CSCSegmentValidation.h.
Referenced by CSCSegmentValidation(), and plotResolution().
MonitorElement* CSCSegmentValidation::thedYdZResolutionPlots[10] [private] |
Definition at line 48 of file CSCSegmentValidation.h.
Referenced by CSCSegmentValidation(), and plotResolution().
Definition at line 28 of file CSCSegmentValidation.h.
Referenced by fillEfficiencyPlots(), and fillLayerHitsPerChamber().
Definition at line 41 of file CSCSegmentValidation.h.
Referenced by analyze().
Definition at line 39 of file CSCSegmentValidation.h.
Referenced by analyze().
Definition at line 40 of file CSCSegmentValidation.h.
Referenced by analyze().
MonitorElement* CSCSegmentValidation::theRdPhiPullPlots[10] [private] |
Definition at line 43 of file CSCSegmentValidation.h.
Referenced by CSCSegmentValidation(), and plotResolution().
MonitorElement* CSCSegmentValidation::theRdPhiResolutionPlots[10] [private] |
Definition at line 42 of file CSCSegmentValidation.h.
Referenced by CSCSegmentValidation(), and plotResolution().
int CSCSegmentValidation::theShowerThreshold [private] |
Definition at line 37 of file CSCSegmentValidation.h.
Referenced by fillEfficiencyPlots().
MonitorElement* CSCSegmentValidation::theThetaPullPlots[10] [private] |
Definition at line 45 of file CSCSegmentValidation.h.
Referenced by CSCSegmentValidation().
MonitorElement* CSCSegmentValidation::theThetaResolutionPlots[10] [private] |
Definition at line 44 of file CSCSegmentValidation.h.
Referenced by CSCSegmentValidation(), and plotResolution().
Definition at line 53 of file CSCSegmentValidation.h.
Referenced by fillEfficiencyPlots().
Definition at line 54 of file CSCSegmentValidation.h.
Referenced by fillEfficiencyPlots().
Definition at line 55 of file CSCSegmentValidation.h.
Referenced by fillEfficiencyPlots().
Definition at line 56 of file CSCSegmentValidation.h.
Referenced by fillEfficiencyPlots().
Definition at line 57 of file CSCSegmentValidation.h.
Referenced by fillEfficiencyPlots().
Definition at line 58 of file CSCSegmentValidation.h.
Referenced by fillEfficiencyPlots().
Definition at line 59 of file CSCSegmentValidation.h.
Referenced by fillEfficiencyPlots().
Definition at line 60 of file CSCSegmentValidation.h.
Referenced by fillEfficiencyPlots().
Definition at line 61 of file CSCSegmentValidation.h.
Referenced by fillEfficiencyPlots().
Definition at line 62 of file CSCSegmentValidation.h.
Referenced by fillEfficiencyPlots().
Definition at line 63 of file CSCSegmentValidation.h.
Referenced by fillEfficiencyPlots().
Definition at line 64 of file CSCSegmentValidation.h.
Referenced by fillEfficiencyPlots().
1.7.3