|
|
|
#include <TrackerMuonHitExtractor.h>
Public Member Functions | |
| std::vector< const TrackingRecHit * > | getMuonHits (const reco::Muon &mu) const |
| void | init (const edm::Event &, const edm::EventSetup &) |
| TrackerMuonHitExtractor (const edm::ParameterSet &) | |
| ~TrackerMuonHitExtractor () | |
Private Attributes | |
| edm::Handle< CSCSegmentCollection > | cscSegmentCollectionH_ |
| edm::Handle < DTRecSegment4DCollection > | dtSegmentCollectionH_ |
| edm::InputTag | inputCSCSegmentCollection_ |
| edm::InputTag | inputDTRecSegment4DCollection_ |
Definition at line 16 of file TrackerMuonHitExtractor.h.
| TrackerMuonHitExtractor::TrackerMuonHitExtractor | ( | const edm::ParameterSet & | parset | ) | [explicit] |
Definition at line 15 of file TrackerMuonHitExtractor.cc.
: inputDTRecSegment4DCollection_(parset.getParameter<edm::InputTag>("inputDTRecSegment4DCollection")), inputCSCSegmentCollection_(parset.getParameter<edm::InputTag>("inputCSCSegmentCollection")) { }
| TrackerMuonHitExtractor::~TrackerMuonHitExtractor | ( | ) |
Definition at line 21 of file TrackerMuonHitExtractor.cc.
{
}
| std::vector< const TrackingRecHit * > TrackerMuonHitExtractor::getMuonHits | ( | const reco::Muon & | mu | ) | const |
Definition at line 93 of file TrackerMuonHitExtractor.cc.
References reco::MuonSegmentMatch::BelongsToTrackByDR, reco::MuonSegmentMatch::BestInChamberByDR, CSCDetId::chamber(), MuonSubdetId::CSC, CSCDetId, MuonSubdetId::DT, DTChamberId, Reference_intrackfit_cff::endcap, CSCDetId::endcap(), edm::Ref< C, T, F >::get(), DTRecSegment4D::hasPhi(), DTRecSegment4D::hasZed(), edm::Ref< C, T, F >::key(), reco::Muon::matches(), reco::Muon::numberOfMatches(), DTRecSegment4D::phiSegment(), CSCSegment::recHits(), DTRecSegment2D::recHits(), run_regression::ret, CSCDetId::ring(), relativeConstraints::ring, DTChamberId::sector(), reco::Muon::SegmentAndTrackArbitration, relativeConstraints::station, DTChamberId::station(), CSCDetId::station(), and DTChamberId::wheel().
Referenced by MuonAssociatorByHits::associateMuons().
{
std::vector<const TrackingRecHit *> ret;
int wheel, station, sector;
int endcap, /*station, */ ring, chamber;
edm::LogVerbatim("TrackerMuonHitExtractor") <<"Number of chambers: "<<mu.matches().size()
<<", arbitrated: "<<mu.numberOfMatches(reco::Muon::SegmentAndTrackArbitration);
unsigned int index_chamber = 0;
for(std::vector<reco::MuonChamberMatch>::const_iterator chamberMatch = mu.matches().begin();
chamberMatch != mu.matches().end(); ++chamberMatch, index_chamber++) {
std::stringstream chamberStr;
chamberStr <<"\nchamber index: "<<index_chamber;
int subdet = chamberMatch->detector();
DetId did = chamberMatch->id;
if (subdet == MuonSubdetId::DT) {
DTChamberId dtdetid = DTChamberId(did);
wheel = dtdetid.wheel();
station = dtdetid.station();
sector = dtdetid.sector();
chamberStr << ", DT chamber Wh:"<<wheel<<",St:"<<station<<",Se:"<<sector;
}
else if (subdet == MuonSubdetId::CSC) {
CSCDetId cscdetid = CSCDetId(did);
endcap = cscdetid.endcap();
station = cscdetid.station();
ring = cscdetid.ring();
chamber = cscdetid.chamber();
chamberStr << ", CSC chamber End:"<<endcap<<",St:"<<station<<",Ri:"<<ring<<",Ch:"<<chamber;
}
chamberStr << ", Number of segments: "<<chamberMatch->segmentMatches.size();
edm::LogVerbatim("TrackerMuonHitExtractor") << chamberStr.str();
unsigned int index_segment = 0;
for(std::vector<reco::MuonSegmentMatch>::const_iterator segmentMatch = chamberMatch->segmentMatches.begin();
segmentMatch != chamberMatch->segmentMatches.end(); ++segmentMatch, index_segment++) {
float segmentX = segmentMatch->x;
float segmentY = segmentMatch->y ;
float segmentdXdZ = segmentMatch->dXdZ;
float segmentdYdZ = segmentMatch->dYdZ;
float segmentXerr = segmentMatch->xErr;
float segmentYerr = segmentMatch->yErr;
float segmentdXdZerr = segmentMatch->dXdZErr;
float segmentdYdZerr = segmentMatch->dYdZErr;
CSCSegmentRef segmentCSC = segmentMatch->cscSegmentRef;
DTRecSegment4DRef segmentDT = segmentMatch->dtSegmentRef;
bool segment_arbitrated_Ok = (segmentMatch->isMask(reco::MuonSegmentMatch::BestInChamberByDR) &&
segmentMatch->isMask(reco::MuonSegmentMatch::BelongsToTrackByDR));
std::string ARBITRATED(" ***Arbitrated Off*** ");
if (segment_arbitrated_Ok) ARBITRATED = " ***ARBITRATED OK*** ";
if (subdet == MuonSubdetId::DT) {
edm::LogVerbatim("TrackerMuonHitExtractor")
<<"\n\t segment index: "<<index_segment << ARBITRATED
<<"\n\t Local Position (X,Y)=("<<segmentX<<","<<segmentY<<") +/- ("<<segmentXerr<<","<<segmentYerr<<"), "
<<"\n\t Local Direction (dXdZ,dYdZ)=("<<segmentdXdZ<<","<<segmentdYdZ<<") +/- ("<<segmentdXdZerr<<","<<segmentdYdZerr<<")";
if (!segment_arbitrated_Ok) continue;
if (segmentDT.get() != 0) {
const DTRecSegment4D* segment = segmentDT.get();
edm::LogVerbatim("TrackerMuonHitExtractor")<<"\t ===> MATCHING with DT segment with index = "<<segmentDT.key();
if(segment->hasPhi()) {
const DTChamberRecSegment2D* phiSeg = segment->phiSegment();
std::vector<const TrackingRecHit*> phiHits = phiSeg->recHits();
for(std::vector<const TrackingRecHit*>::const_iterator ihit = phiHits.begin();
ihit != phiHits.end(); ++ihit) {
ret.push_back(*ihit);
}
}
if(segment->hasZed()) {
const DTSLRecSegment2D* zSeg = (*segment).zSegment();
std::vector<const TrackingRecHit*> zedHits = zSeg->recHits();
for(std::vector<const TrackingRecHit*>::const_iterator ihit = zedHits.begin();
ihit != zedHits.end(); ++ihit) {
ret.push_back(*ihit);
}
}
} else edm::LogWarning("TrackerMuonHitExtractor")<<"\n***WARNING: UNMATCHED DT segment ! \n";
} // if (subdet == MuonSubdetId::DT)
else if (subdet == MuonSubdetId::CSC) {
edm::LogVerbatim("TrackerMuonHitExtractor")
<<"\n\t segment index: "<<index_segment << ARBITRATED
<<"\n\t Local Position (X,Y)=("<<segmentX<<","<<segmentY<<") +/- ("<<segmentXerr<<","<<segmentYerr<<"), "
<<"\n\t Local Direction (dXdZ,dYdZ)=("<<segmentdXdZ<<","<<segmentdYdZ<<") +/- ("<<segmentdXdZerr<<","<<segmentdYdZerr<<")";
if (!segment_arbitrated_Ok) continue;
if (segmentCSC.get() != 0) {
const CSCSegment* segment = segmentCSC.get();
edm::LogVerbatim("TrackerMuonHitExtractor")<<"\t ===> MATCHING with CSC segment with index = "<<segmentCSC.key();
std::vector<const TrackingRecHit*> hits = segment->recHits();
for(std::vector<const TrackingRecHit*>::const_iterator ihit = hits.begin();
ihit != hits.end(); ++ihit) {
ret.push_back(*ihit);
}
} else edm::LogWarning("TrackerMuonHitExtractor")<<"\n***WARNING: UNMATCHED CSC segment ! \n";
} // else if (subdet == MuonSubdetId::CSC)
} // loop on vector<MuonSegmentMatch>
} // loop on vector<MuonChamberMatch>
return ret;
}
| void TrackerMuonHitExtractor::init | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) |
Definition at line 24 of file TrackerMuonHitExtractor.cc.
References CSCDetId::chamber(), CSCDetId, cscSegmentCollectionH_, DTChamberId, dtSegmentCollectionH_, Reference_intrackfit_cff::endcap, CSCDetId::endcap(), edm::Event::getByLabel(), inputCSCSegmentCollection_, inputDTRecSegment4DCollection_, CSCDetId::ring(), relativeConstraints::ring, DTChamberId::sector(), mathSSE::sqrt(), relativeConstraints::station, CSCDetId::station(), DTChamberId::station(), DTChamberId::wheel(), PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), PV3DBase< T, PVType, FrameType >::y(), LocalError::yy(), and PV3DBase< T, PVType, FrameType >::z().
Referenced by MuonAssociatorByHits::associateMuons().
{
iEvent.getByLabel(inputDTRecSegment4DCollection_, dtSegmentCollectionH_);
iEvent.getByLabel(inputCSCSegmentCollection_, cscSegmentCollectionH_);
edm::LogVerbatim("TrackerMuonHitExtractor") <<"\nThere are "<< dtSegmentCollectionH_->size()<<" DT segments.";
unsigned int index_dt_segment = 0;
for(DTRecSegment4DCollection::const_iterator segment = dtSegmentCollectionH_->begin();
segment != dtSegmentCollectionH_->end(); ++segment , index_dt_segment++) {
LocalPoint segmentLocalPosition = segment->localPosition();
LocalVector segmentLocalDirection = segment->localDirection();
LocalError segmentLocalPositionError = segment->localPositionError();
LocalError segmentLocalDirectionError = segment->localDirectionError();
DetId geoid = segment->geographicalId();
DTChamberId dtdetid = DTChamberId(geoid);
int wheel = dtdetid.wheel();
int station = dtdetid.station();
int sector = dtdetid.sector();
float segmentX = segmentLocalPosition.x();
float segmentY = segmentLocalPosition.y();
float segmentdXdZ = segmentLocalDirection.x()/segmentLocalDirection.z();
float segmentdYdZ = segmentLocalDirection.y()/segmentLocalDirection.z();
float segmentXerr = sqrt(segmentLocalPositionError.xx());
float segmentYerr = sqrt(segmentLocalPositionError.yy());
float segmentdXdZerr = sqrt(segmentLocalDirectionError.xx());
float segmentdYdZerr = sqrt(segmentLocalDirectionError.yy());
edm::LogVerbatim("TrackerMuonHitExtractor")
<<"\nDT segment index :"<<index_dt_segment
<<"\nchamber Wh:"<<wheel<<",St:"<<station<<",Se:"<<sector
<<"\nLocal Position (X,Y)=("<<segmentX<<","<<segmentY<<") +/- ("<<segmentXerr<<","<<segmentYerr<<"), "
<<"Local Direction (dXdZ,dYdZ)=("<<segmentdXdZ<<","<<segmentdYdZ<<") +/- ("<<segmentdXdZerr<<","<<segmentdYdZerr<<")";
}
edm::LogVerbatim("TrackerMuonHitExtractor") <<"\nThere are "<< cscSegmentCollectionH_->size()<<" CSC segments.";
unsigned int index_csc_segment = 0;
for(CSCSegmentCollection::const_iterator segment = cscSegmentCollectionH_->begin();
segment != cscSegmentCollectionH_->end(); ++segment , index_csc_segment++) {
LocalPoint segmentLocalPosition = segment->localPosition();
LocalVector segmentLocalDirection = segment->localDirection();
LocalError segmentLocalPositionError = segment->localPositionError();
LocalError segmentLocalDirectionError = segment->localDirectionError();
DetId geoid = segment->geographicalId();
CSCDetId cscdetid = CSCDetId(geoid);
int endcap = cscdetid.endcap();
int station = cscdetid.station();
int ring = cscdetid.ring();
int chamber = cscdetid.chamber();
float segmentX = segmentLocalPosition.x();
float segmentY = segmentLocalPosition.y();
float segmentdXdZ = segmentLocalDirection.x()/segmentLocalDirection.z();
float segmentdYdZ = segmentLocalDirection.y()/segmentLocalDirection.z();
float segmentXerr = sqrt(segmentLocalPositionError.xx());
float segmentYerr = sqrt(segmentLocalPositionError.yy());
float segmentdXdZerr = sqrt(segmentLocalDirectionError.xx());
float segmentdYdZerr = sqrt(segmentLocalDirectionError.yy());
edm::LogVerbatim("TrackerMuonHitExtractor")
<<"\nCSC segment index :"<<index_csc_segment
<<"\nchamber Endcap:"<<endcap<<",St:"<<station<<",Ri:"<<ring<<",Ch:"<<chamber
<<"\nLocal Position (X,Y)=("<<segmentX<<","<<segmentY<<") +/- ("<<segmentXerr<<","<<segmentYerr<<"), "
<<"Local Direction (dXdZ,dYdZ)=("<<segmentdXdZ<<","<<segmentdYdZ<<") +/- ("<<segmentdXdZerr<<","<<segmentdYdZerr<<")";
}
}
Definition at line 25 of file TrackerMuonHitExtractor.h.
Referenced by init().
Definition at line 24 of file TrackerMuonHitExtractor.h.
Referenced by init().
Definition at line 28 of file TrackerMuonHitExtractor.h.
Referenced by init().
Definition at line 27 of file TrackerMuonHitExtractor.h.
Referenced by init().
1.7.3