SegmentToTrackAssociator Class Reference

#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 &, const GlobalTrackingGeometry *GlobalTrackingGeometry, edm::ConsumesCollector &)
virtual	~SegmentToTrackAssociator ()=default

Private Attributes
const GlobalTrackingGeometry *	globalTrackingGeometry_
intDVector	indexCollectionCSC
intDVector	indexCollectionDT
const edm::InputTag	theCSCSegmentLabel
const edm::InputTag	theDTSegmentLabel
const edm::EDGetTokenT< CSCSegmentCollection >	tokenCSCSegment_
const edm::EDGetTokenT< DTRecSegment4DCollection >	tokenDTSegment_

Detailed Description

Definition at line 32 of file SegmentToTrackAssociator.h.

Member Typedef Documentation

intDVector

typedef std::vector<std::vector<int> > SegmentToTrackAssociator::intDVector

Definition at line 34 of file SegmentToTrackAssociator.h.

Constructor & Destructor Documentation

SegmentToTrackAssociator()

SegmentToTrackAssociator::SegmentToTrackAssociator	(	const edm::ParameterSet &	iConfig,
		const GlobalTrackingGeometry *	GlobalTrackingGeometry,
		edm::ConsumesCollector &	iC
	)

Definition at line 33 of file SegmentToTrackAssociator.cc.

    : globalTrackingGeometry_(GlobalTrackingGeometry),
      theDTSegmentLabel(iConfig.getParameter<edm::InputTag>("segmentsDT")),
      theCSCSegmentLabel(iConfig.getParameter<edm::InputTag>("segmentsCSC")),
      tokenDTSegment_(iC.consumes<DTRecSegment4DCollection>(theDTSegmentLabel)),
      tokenCSCSegment_(iC.consumes<CSCSegmentCollection>(theCSCSegmentLabel)) {}

~SegmentToTrackAssociator()

virtual SegmentToTrackAssociator::~SegmentToTrackAssociator ( )

virtualdefault

Member Function Documentation

associate()

MuonTransientTrackingRecHit::MuonRecHitContainer SegmentToTrackAssociator::associate	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		const reco::Track &	track,
		std::string	TrackRefitterType
	)

Definition at line 47 of file SegmentToTrackAssociator.cc.

References GeomDetEnumerators::CSC, GeomDetEnumerators::DT, TrackingRecHit::geographicalId(), edm::Ref< C, T, F >::get(), globalTrackingGeometry_, GlobalTrackingGeometry::idToDet(), iEvent, indexCollectionCSC, indexCollectionDT, DetId::rawId(), groupFilesInBlocks::reverse, segmentTrackAnalyzer_cfi::SelectedSegments, MuonTransientTrackingRecHit::specificBuild(), GeomDet::subDetector(), tokenCSCSegment_, tokenDTSegment_, HLT_2023v11_cff::track, and trackerHitRTTI::vector.

Referenced by SequenceTypes._ModuleSequenceType::remove(), and SequenceTypes._ModuleSequenceType::replace().

                                                                                                                  {
  // The segment collections
  const edm::Handle<DTRecSegment4DCollection>& DTSegments = iEvent.getHandle(tokenDTSegment_);

  const edm::Handle<CSCSegmentCollection>& CSCSegments = iEvent.getHandle(tokenCSCSegment_);

  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 = globalTrackingGeometry_->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 = globalTrackingGeometry_->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 = globalTrackingGeometry_->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;
}

clear()

void SegmentToTrackAssociator::clear

(

)

Definition at line 42 of file SegmentToTrackAssociator.cc.

References indexCollectionCSC, and indexCollectionDT.

                                     {
  indexCollectionDT.clear();
  indexCollectionCSC.clear();
}

Member Data Documentation

globalTrackingGeometry_

const GlobalTrackingGeometry* SegmentToTrackAssociator::globalTrackingGeometry_

private

Definition at line 57 of file SegmentToTrackAssociator.h.

Referenced by associate().

indexCollectionCSC

intDVector SegmentToTrackAssociator::indexCollectionCSC

private

Definition at line 55 of file SegmentToTrackAssociator.h.

Referenced by associate(), and clear().

indexCollectionDT

intDVector SegmentToTrackAssociator::indexCollectionDT

private

Definition at line 54 of file SegmentToTrackAssociator.h.

Referenced by associate(), and clear().

theCSCSegmentLabel

const edm::InputTag SegmentToTrackAssociator::theCSCSegmentLabel

private

Definition at line 60 of file SegmentToTrackAssociator.h.

theDTSegmentLabel

const edm::InputTag SegmentToTrackAssociator::theDTSegmentLabel

private

Definition at line 59 of file SegmentToTrackAssociator.h.

tokenCSCSegment_

const edm::EDGetTokenT<CSCSegmentCollection> SegmentToTrackAssociator::tokenCSCSegment_

private

Definition at line 62 of file SegmentToTrackAssociator.h.

Referenced by associate().

tokenDTSegment_

const edm::EDGetTokenT<DTRecSegment4DCollection> SegmentToTrackAssociator::tokenDTSegment_

private

Definition at line 61 of file SegmentToTrackAssociator.h.

Referenced by associate().