|
|
|
#include <SegmentsTrackAssociator.h>
Public Member Functions | |
| MuonTransientTrackingRecHit::MuonRecHitContainer | associate (const edm::Event &, const edm::EventSetup &, const reco::Track &) |
| Get the analysis. | |
| SegmentsTrackAssociator (const edm::ParameterSet &) | |
| Constructor. | |
| virtual | ~SegmentsTrackAssociator () |
| Destructor. | |
Private Attributes | |
| std::string | metname |
| int | numRecHit |
| int | numRecHitCSC |
| int | numRecHitDT |
| edm::InputTag | theCSCSegmentLabel |
| edm::InputTag | theDTSegmentLabel |
| edm::InputTag | theSegmentContainerName |
tool which take as input a muon track and return a vector with the segments used to fit it
Definition at line 37 of file SegmentsTrackAssociator.h.
| SegmentsTrackAssociator::SegmentsTrackAssociator | ( | const edm::ParameterSet & | iConfig | ) |
Constructor.
Definition at line 44 of file SegmentsTrackAssociator.cc.
References edm::ParameterSet::getUntrackedParameter(), and metname.
{
theDTSegmentLabel = iConfig.getUntrackedParameter<InputTag>("segmentsDt");
theCSCSegmentLabel = iConfig.getUntrackedParameter<InputTag>("segmentsCSC");
theSegmentContainerName = iConfig.getUntrackedParameter<InputTag>("SelectedSegments");
numRecHit=0;
numRecHitDT=0;
numRecHitCSC=0;
metname = "SegmentsTrackAssociator";
}
| SegmentsTrackAssociator::~SegmentsTrackAssociator | ( | ) | [inline, virtual] |
| MuonTransientTrackingRecHit::MuonRecHitContainer SegmentsTrackAssociator::associate | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup, | ||
| const reco::Track & | Track | ||
| ) |
Get the analysis.
Definition at line 60 of file SegmentsTrackAssociator.cc.
References edm::check(), filterCSVwithJSON::copy, CSC(), cscSegments_cfi::cscSegments, DetId::det(), GeomDetEnumerators::DT, Reference_intrackfit_cff::endcap, edm::EventSetup::get(), edm::Event::getByLabel(), LogTrace, metname, DetId::Muon, n, funct::pow(), DetId::rawId(), reco::Track::recHitsBegin(), reco::Track::recHitsEnd(), relativeConstraints::ring, MuonTransientTrackingRecHit::specificBuild(), mathSSE::sqrt(), relativeConstraints::station, DetId::subdetId(), PV3DBase< T, PVType, FrameType >::x(), x, PV3DBase< T, PVType, FrameType >::y(), detailsBasic3DVector::y, PV3DBase< T, PVType, FrameType >::z(), and z.
Referenced by SegmentTrackAnalyzer::analyze().
{
// The segment collections
Handle<DTRecSegment4DCollection> dtSegments;
iEvent.getByLabel(theDTSegmentLabel, dtSegments);
Handle<CSCSegmentCollection> cscSegments;
iEvent.getByLabel(theCSCSegmentLabel, cscSegments);
ESHandle<GlobalTrackingGeometry> theTrackingGeometry;
iSetup.get<GlobalTrackingGeometryRecord>().get(theTrackingGeometry);
DTRecSegment4DCollection::const_iterator segment;
CSCSegmentCollection::const_iterator segment2;
MuonTransientTrackingRecHit::MuonRecHitContainer selectedSegments;
MuonTransientTrackingRecHit::MuonRecHitContainer selectedDtSegments;
MuonTransientTrackingRecHit::MuonRecHitContainer selectedCscSegments;
//loop over recHit
for(trackingRecHit_iterator recHit = Track.recHitsBegin(); recHit != Track.recHitsEnd(); ++recHit){
if(!(*recHit)->isValid()) continue;
numRecHit++;
//get the detector Id
DetId idRivHit = (*recHit)->geographicalId();
// DT recHits
if (idRivHit.det() == DetId::Muon && idRivHit.subdetId() == MuonSubdetId::DT ) {
// get the RecHit Local Position
LocalPoint posTrackRecHit = (*recHit)->localPosition();
DTRecSegment4DCollection::range range;
numRecHitDT++;
// get the chamber Id
DTChamberId chamberId(idRivHit.rawId());
// get the segments of the chamber
range = dtSegments->get(chamberId);
// loop over segments
for (segment = range.first; segment!=range.second; segment++){
DetId id = segment->geographicalId();
const GeomDet* det = theTrackingGeometry->idToDet(id);
vector<const TrackingRecHit*> segments2D = (&(*segment))->recHits();
// container for 4D segment recHits
vector <const TrackingRecHit*> dtRecHits;
for(vector<const TrackingRecHit*>::const_iterator segm2D = segments2D.begin();
segm2D != segments2D.end();
segm2D++) {
vector <const TrackingRecHit*> rHits1D = (*segm2D)->recHits();
for (int hit=0; hit<int(rHits1D.size()); hit++){
dtRecHits.push_back(rHits1D[hit]);
}
}
// loop over the recHit checking if there's the recHit of the track
for (unsigned int hit = 0; hit < dtRecHits.size(); hit++) {
DetId idRivHitSeg = (*dtRecHits[hit]).geographicalId();
LocalPoint posDTSegment= segment->localPosition();
LocalPoint posSegDTRecHit = (*dtRecHits[hit]).localPosition();
double rDT=sqrt(pow((posSegDTRecHit.x()-posTrackRecHit.x()),2) +pow((posSegDTRecHit.y()-posTrackRecHit.y()),2) + pow((posSegDTRecHit.z()-posTrackRecHit.z()),2));
if (idRivHit == idRivHitSeg && rDT<0.0001){
if (selectedDtSegments.empty()){
selectedDtSegments.push_back(MuonTransientTrackingRecHit::specificBuild(det,&*segment));
LogTrace(metname) <<"First selected segment (from DT). Position : "<<posDTSegment<<" Chamber : "<<segment->chamberId();
}
else{
int check=0;
for(int segm=0; segm < int(selectedDtSegments.size()); ++segm) {
double dist=sqrt(pow((((*(selectedDtSegments[segm])).localPosition()).x()-posDTSegment.x()),2) +pow((((*(selectedDtSegments[segm])).localPosition()).y()-posDTSegment.y()),2) + pow((((*(selectedDtSegments[segm])).localPosition()).z()-posDTSegment.z()),2));
if(dist>30) check++;
}
if(check==int(selectedDtSegments.size())){
selectedDtSegments.push_back(MuonTransientTrackingRecHit::specificBuild(det,&*segment));
LogTrace(metname) <<"New DT selected segment. Position : "<<posDTSegment<<" Chamber : "<<segment->chamberId();
}
}
} // check to tag the segment as "selected"
} // loop over segment recHits
} // loop over DT segments
}
// CSC recHits
if (idRivHit.det() == DetId::Muon && idRivHit.subdetId() == MuonSubdetId::CSC ) {
// get the RecHit Local Position
LocalPoint posTrackRecHit = (*recHit)->localPosition();
CSCSegmentCollection::range range;
numRecHitCSC++;
// get the chamber Id
CSCDetId tempchamberId(idRivHit.rawId());
int ring = tempchamberId.ring();
int station = tempchamberId.station();
int endcap = tempchamberId.endcap();
int chamber = tempchamberId.chamber();
CSCDetId chamberId(endcap, station, ring, chamber, 0);
// get the segments of the chamber
range = cscSegments->get(chamberId);
// loop over segments
for(segment2 = range.first; segment2!=range.second; segment2++){
DetId id2 = segment2->geographicalId();
const GeomDet* det2 = theTrackingGeometry->idToDet(id2);
// container for CSC segment recHits
vector<const TrackingRecHit*> cscRecHits = (&(*segment2))->recHits();
// loop over the recHit checking if there's the recHit of the track
for (unsigned int hit = 0; hit < cscRecHits.size(); hit++) {
DetId idRivHitSeg = (*cscRecHits[hit]).geographicalId();
LocalPoint posSegCSCRecHit = (*cscRecHits[hit]).localPosition();
LocalPoint posCSCSegment= segment2->localPosition();
double rCSC=sqrt(pow((posSegCSCRecHit.x()-posTrackRecHit.x()),2) +pow((posSegCSCRecHit.y()-posTrackRecHit.y()),2) + pow((posSegCSCRecHit.z()-posTrackRecHit.z()),2));
if (idRivHit==idRivHitSeg && rCSC < 0.0001){
if (selectedCscSegments.empty()){
selectedCscSegments.push_back(MuonTransientTrackingRecHit::specificBuild(det2,&*segment2));
LogTrace(metname) <<"First selected segment (from CSC). Position: "<<posCSCSegment;
}
else{
int check=0;
for(int n=0; n< int(selectedCscSegments.size()); n++){
double dist = sqrt(pow(((*(selectedCscSegments[n])).localPosition().x()-posCSCSegment.x()),2) +pow(((*(selectedCscSegments[n])).localPosition().y()-posCSCSegment.y()),2) + pow(((*(selectedCscSegments[n])).localPosition().z()-posCSCSegment.z()),2));
if(dist>30) check++;
}
if(check==int(selectedCscSegments.size())){
selectedCscSegments.push_back(MuonTransientTrackingRecHit::specificBuild(det2,&*segment2));
LogTrace(metname) <<"New CSC segment selected. Position : "<<posCSCSegment;
}
}
} // check to tag the segment as "selected"
} // loop over segment recHits
} // loop over DT segments
}
} // loop over track recHits
LogTrace(metname) <<"N recHit:"<<numRecHit;
numRecHit=0;
LogTrace(metname) <<"N recHit DT:"<<numRecHitDT;
numRecHitDT=0;
LogTrace(metname) <<"N recHit CSC:"<<numRecHitCSC;
numRecHitCSC=0;
copy(selectedDtSegments.begin(), selectedDtSegments.end(), back_inserter(selectedSegments));
LogTrace(metname) <<"N selected Dt segments:"<<selectedDtSegments.size();
copy(selectedCscSegments.begin(), selectedCscSegments.end(), back_inserter(selectedSegments));
LogTrace(metname) <<"N selected Csc segments:"<<selectedCscSegments.size();
LogTrace(metname) <<"N selected segments:"<<selectedSegments.size();
return selectedSegments;
}
std::string SegmentsTrackAssociator::metname [private] |
Definition at line 64 of file SegmentsTrackAssociator.h.
int SegmentsTrackAssociator::numRecHit [private] |
Definition at line 55 of file SegmentsTrackAssociator.h.
int SegmentsTrackAssociator::numRecHitCSC [private] |
Definition at line 57 of file SegmentsTrackAssociator.h.
int SegmentsTrackAssociator::numRecHitDT [private] |
Definition at line 56 of file SegmentsTrackAssociator.h.
Definition at line 61 of file SegmentsTrackAssociator.h.
Definition at line 60 of file SegmentsTrackAssociator.h.
Definition at line 62 of file SegmentsTrackAssociator.h.
1.7.3