|
|
|
#include <SegmentToTrackAssociator.h>
Public Types | |
| typedef std::vector < std::vector< int > > | intDVector |
Public Member Functions | |
| MuonTransientTrackingRecHit::MuonRecHitContainer | associate (const edm::Event &, const edm::EventSetup &, const reco::Track &, std::string) |
| void | clear () |
| SegmentToTrackAssociator (const edm::ParameterSet &) | |
| virtual | ~SegmentToTrackAssociator () |
Private Attributes | |
| intDVector | indexCollectionCSC |
| intDVector | indexCollectionDT |
| edm::InputTag | theCSCSegmentLabel |
| edm::InputTag | theDTSegmentLabel |
Definition at line 25 of file SegmentToTrackAssociator.h.
| typedef std::vector< std::vector<int> > SegmentToTrackAssociator::intDVector |
Definition at line 28 of file SegmentToTrackAssociator.h.
| SegmentToTrackAssociator::SegmentToTrackAssociator | ( | const edm::ParameterSet & | iConfig | ) |
Definition at line 39 of file SegmentToTrackAssociator.cc.
References edm::ParameterSet::getParameter(), theCSCSegmentLabel, and theDTSegmentLabel.
{
theDTSegmentLabel = iConfig.getParameter<edm::InputTag>( "segmentsDT" );
theCSCSegmentLabel = iConfig.getParameter<edm::InputTag>( "segmentsCSC" );
}
| SegmentToTrackAssociator::~SegmentToTrackAssociator | ( | ) | [inline, virtual] |
Definition at line 34 of file SegmentToTrackAssociator.h.
{}
| MuonTransientTrackingRecHit::MuonRecHitContainer SegmentToTrackAssociator::associate | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup, | ||
| const reco::Track & | track, | ||
| std::string | TrackRefitterType | ||
| ) |
Definition at line 58 of file SegmentToTrackAssociator.cc.
References CSC(), GeomDetEnumerators::DT, TrackingRecHit::geographicalId(), edm::EventSetup::get(), edm::Ref< C, T, F >::get(), edm::Event::getByLabel(), indexCollectionCSC, indexCollectionDT, DetId::rawId(), reco::Track::recHit(), reco::Track::recHitsSize(), MuonTransientTrackingRecHit::specificBuild(), GeomDet::subDetector(), theCSCSegmentLabel, and theDTSegmentLabel.
{
// The segment collections
edm::Handle<DTRecSegment4DCollection> DTSegments;
iEvent.getByLabel(theDTSegmentLabel, DTSegments);
edm::Handle<CSCSegmentCollection> CSCSegments;
iEvent.getByLabel(theCSCSegmentLabel, CSCSegments);
edm::ESHandle<GlobalTrackingGeometry> theTrackingGeometry;
iSetup.get<GlobalTrackingGeometryRecord>().get( theTrackingGeometry );
MuonTransientTrackingRecHit::MuonRecHitContainer SelectedSegments;
DTRecSegment4DCollection::const_iterator segmentDT;
CSCSegmentCollection::const_iterator segmentCSC;
std::vector<int> positionDT;
std::vector<int> positionCSC;
std::vector<TrackingRecHit *> my4DTrack;
//Loop over the hits of the track
for( unsigned int counter = 0; counter != track.recHitsSize()-1; counter++ )
{
TrackingRecHitRef myRef = track.recHit( counter );
const TrackingRecHit *rechit = myRef.get();
const GeomDet* geomDet = theTrackingGeometry->idToDet( rechit->geographicalId() );
//It's a DT Hit
if( geomDet->subDetector() == GeomDetEnumerators::DT )
{
//Take the layer associated to this hit
DTLayerId myLayer( rechit->geographicalId().rawId() );
int NumberOfDTSegment = 0;
//Loop over segments
for( segmentDT = DTSegments->begin(); segmentDT != DTSegments->end(); ++segmentDT ) {
//By default the chamber associated to this Segment is new
bool isNewChamber = true;
//Loop over segments already included in the vector of segments in the actual track
for( std::vector<int>::iterator positionIt = positionDT.begin();
positionIt != positionDT.end(); positionIt++ )
{
//If this segment has been used before isNewChamber = false
if(NumberOfDTSegment == *positionIt) isNewChamber = false;
}
//Loop over vectors of segments associated to previous tracks
for( std::vector<std::vector<int> >::iterator collect = indexCollectionDT.begin();
collect != indexCollectionDT.end(); ++collect)
{
//Loop over segments associated to a track
for( std::vector<int>::iterator positionIt = (*collect).begin();
positionIt != (*collect).end(); positionIt++ )
{
//If this segment was used in a previos track then isNewChamber = false
if( NumberOfDTSegment == *positionIt ) isNewChamber = false;
}
}
//If the chamber is new
if( isNewChamber )
{
DTChamberId myChamber( (*segmentDT).geographicalId().rawId() );
//If the layer of the hit belongs to the chamber of the 4D Segment
if( myLayer.wheel() == myChamber.wheel() &&
myLayer.station() == myChamber.station() &&
myLayer.sector() == myChamber.sector() )
{
//push position of the segment and tracking rechit
positionDT.push_back( NumberOfDTSegment );
const GeomDet* DTgeomDet = theTrackingGeometry->idToDet( myChamber );
SelectedSegments.push_back( MuonTransientTrackingRecHit::specificBuild( DTgeomDet, (TrackingRecHit *) &*segmentDT ) );
//edm::LogWarning("Alignment") << "TagSeg: " << "NumberOfDTSegment " << NumberOfDTSegment << " Wheel " << myChamber.wheel() << " Sector " << myChamber.sector() << " Chamber " << myChamber.station() << std::endl;
}
}
NumberOfDTSegment++;
}
//In case is a CSC
}
else if ( geomDet->subDetector() == GeomDetEnumerators::CSC )
{
//Take the layer associated to this hit
CSCDetId myLayer( rechit->geographicalId().rawId() );
int NumberOfCSCSegment = 0;
//Loop over 4Dsegments
for( segmentCSC = CSCSegments->begin(); segmentCSC != CSCSegments->end(); segmentCSC++ )
{
//By default the chamber associated to the segment is new
bool isNewChamber = true;
//Loop over segments in the current track
for( std::vector<int>::iterator positionIt = positionCSC.begin();
positionIt != positionCSC.end(); positionIt++ )
{
//If this segment has been used then newchamber = false
if( NumberOfCSCSegment == *positionIt ) isNewChamber = false;
}
//Loop over vectors of segments in previous tracks
for( std::vector<std::vector<int> >::iterator collect = indexCollectionCSC.begin();
collect != indexCollectionCSC.end(); ++collect )
{
//Loop over segments in a track
for( std::vector<int>::iterator positionIt = (*collect).begin();
positionIt != (*collect).end(); positionIt++ )
{
//If the segment was used in a previous track isNewChamber = false
if( NumberOfCSCSegment == *positionIt ) isNewChamber = false;
}
}
//If the chamber is new
if( isNewChamber )
{
CSCDetId myChamber( (*segmentCSC).geographicalId().rawId() );
//If the chambers are the same
if( myLayer.chamberId() == myChamber.chamberId() )
{
//push
positionCSC.push_back( NumberOfCSCSegment );
const GeomDet* CSCgeomDet = theTrackingGeometry->idToDet( myChamber );
SelectedSegments.push_back(MuonTransientTrackingRecHit::specificBuild( CSCgeomDet, (TrackingRecHit *) &*segmentCSC ) );
}
}
NumberOfCSCSegment++;
}
}
}
indexCollectionDT.push_back(positionDT);
indexCollectionCSC.push_back(positionCSC);
if ( TrackRefitterType == "CosmicLike" )
{
std::reverse(SelectedSegments.begin(),SelectedSegments.end());
}
return SelectedSegments;
}
| void SegmentToTrackAssociator::clear | ( | void | ) |
Definition at line 51 of file SegmentToTrackAssociator.cc.
References indexCollectionCSC, and indexCollectionDT.
{
indexCollectionDT.clear();
indexCollectionCSC.clear();
}
Definition at line 46 of file SegmentToTrackAssociator.h.
Referenced by associate(), and clear().
Definition at line 45 of file SegmentToTrackAssociator.h.
Referenced by associate(), and clear().
Definition at line 49 of file SegmentToTrackAssociator.h.
Referenced by associate(), and SegmentToTrackAssociator().
Definition at line 48 of file SegmentToTrackAssociator.h.
Referenced by associate(), and SegmentToTrackAssociator().
1.7.3