TrackMerger Class Reference

#include <TrackMerger.h>

Classes
class	GlobalMomentumSort
class	MomentumSort

Public Member Functions
void	init (const edm::EventSetup &iSetup)
TrackCandidate	merge (const reco::Track &inner, const reco::Track &outer, DuplicateTrackType duplicateType) const
	TrackMerger (const edm::ParameterSet &iConfig, edm::ConsumesCollector)
	~TrackMerger ()

Private Member Functions
void	addSecondTrackHits (std::vector< const TrackingRecHit * > &hits, const reco::Track &outer) const
void	sortByHitPosition (const GlobalVector &v, const std::vector< const TrackingRecHit * > &hits, TrackCandidate::RecHitContainer &ownHits) const

Private Attributes
edm::ESGetToken< TransientTrackingRecHitBuilder, TransientRecHitRecord >	builderToken_
bool	debug_
edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	geometryToken_
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	magFieldToken_
edm::ESHandle< TransientTrackingRecHitBuilder >	theBuilder
std::string	theBuilderName
edm::ESHandle< TrackerGeometry >	theGeometry
edm::ESHandle< MagneticField >	theMagField
edm::ESHandle< TrackerTopology >	theTrkTopo
edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	trackerTopoToken_
bool	useInnermostState_

Detailed Description

Definition at line 19 of file TrackMerger.h.

Constructor & Destructor Documentation

TrackMerger()

TrackMerger::TrackMerger	(	const edm::ParameterSet &	iConfig,
		edm::ConsumesCollector	cc
	)

Definition at line 22 of file TrackMerger.cc.

    : useInnermostState_(iConfig.getParameter<bool>("useInnermostState")),
      theBuilderName(iConfig.getParameter<std::string>("ttrhBuilderName")),
      geometryToken_(cc.esConsumes<TrackerGeometry, TrackerDigiGeometryRecord>()),
      magFieldToken_(cc.esConsumes<MagneticField, IdealMagneticFieldRecord>()),
      builderToken_(
          cc.esConsumes<TransientTrackingRecHitBuilder, TransientRecHitRecord>(edm::ESInputTag("", theBuilderName))),
      trackerTopoToken_(cc.esConsumes<TrackerTopology, TrackerTopologyRcd>()) {}

~TrackMerger()

TrackMerger::~TrackMerger

(

)

Definition at line 31 of file TrackMerger.cc.

{}

Member Function Documentation

addSecondTrackHits()

void TrackMerger::addSecondTrackHits ( std::vector< const TrackingRecHit * > &  hits, const reco::Track &  outer  ) const

private

Definition at line 146 of file TrackMerger.cc.

                                                                                                          {
  size_t nHitsFirstTrack = hits.size();
  for (auto it = outer.recHitsBegin(), ed = outer.recHitsEnd(); it != ed; ++it) {
    const TrackingRecHit *hit = &**it;
    DetId id(hit->geographicalId());
    const auto lay = theTrkTopo->layer(id);
    bool shared = false;
    bool valid = hit->isValid();
    PRINT << "   subdet " << id.subdetId() << "  layer " << theTrkTopo->layer(id) << " valid " << valid
          << "   detid: " << id() << std::endl;
    size_t iHit = 0;
    for (auto hit2 : hits) {
      ++iHit;
      if (iHit > nHitsFirstTrack)
        break;
      DetId id2 = hit2->geographicalId();
      if (id.subdetId() != id2.subdetId())
        continue;
      if (theTrkTopo->layer(id2) != lay)
        continue;
      if (hit->sharesInput(hit2, TrackingRecHit::all)) {
        PRINT << "        discared as duplicate of other hit" << id() << std::endl;
        shared = true;
        break;
      }
      if (hit2->isValid() && !valid) {
        PRINT << "        replacing old invalid hit on detid " << id2() << std::endl;
        hit = hit2;
        shared = true;
        break;
      }
      PRINT << "        discared as additional hit on layer that already contains hit with detid " << id() << std::endl;
      shared = true;
      break;
    }
    if (shared)
      continue;
    hits.push_back(hit);
  }
}

References TrackingRecHit::all, hfClusterShapes_cfi::hits, triggerObjects_cff::id, globals_cff::id2, TrackerTopology::layer(), SurfaceOrientation::outer, PRINT, theTrkTopo, and RunInfoPI::valid.

Referenced by merge().

init()

void TrackMerger::init

(

const edm::EventSetup &

iSetup

)

Definition at line 33 of file TrackMerger.cc.

                                                  {
  theGeometry = iSetup.getHandle(geometryToken_);
  theMagField = iSetup.getHandle(magFieldToken_);
  theBuilder = iSetup.getHandle(builderToken_);
  theTrkTopo = iSetup.getHandle(trackerTopoToken_);
}

References builderToken_, geometryToken_, edm::EventSetup::getHandle(), magFieldToken_, theBuilder, theGeometry, theMagField, theTrkTopo, and trackerTopoToken_.

merge()

TrackCandidate TrackMerger::merge	(	const reco::Track &	inner,
		const reco::Track &	outer,
		DuplicateTrackType	duplicateType
	)		const

Definition at line 40 of file TrackMerger.cc.

                                                                          {
  DPRINT("TrackMerger") << std::abs(inner.eta()) << " merging " << inner.algo() << '/' << outer.algo() << ' '
                        << inner.eta() << '/' << outer.eta() << std::endl;
 
  std::vector<const TrackingRecHit *> hits;
  hits.reserve(inner.recHitsSize() + outer.recHitsSize());
  DPRINT("TrackMerger") << "Inner track hits: " << std::endl;
  for (auto it = inner.recHitsBegin(), ed = inner.recHitsEnd(); it != ed; ++it) {
    hits.push_back(&**it);
    if (debug_) {
      DetId id(hits.back()->geographicalId());
      PRINT << "   subdet " << id.subdetId() << "  layer " << theTrkTopo->layer(id) << " valid "
            << hits.back()->isValid() << "   detid: " << id() << std::endl;
    }
  }
  DPRINT("TrackMerger") << "Outer track hits: " << std::endl;
 
  if (duplicateType == DuplicateTrackType::Disjoint) {
#if TRACK_SORT
    DetId lastId(hits.back()->geographicalId());
    int lastSubdet = lastId.subdetId();
    unsigned int lastLayer = theTrkTopo->layer(lastId);
    for (auto it = outer.recHitsBegin(), ed = outer.recHitsEnd(); it != ed; ++it) {
      const TrackingRecHit *hit = &**it;
      DetId id(hit->geographicalId());
      int thisSubdet = id.subdetId();
      if (thisSubdet > lastSubdet || (thisSubdet == lastSubdet && theTrkTopo->layer(id) > lastLayer)) {
        hits.push_back(hit);
        PRINT << "  adding   subdet " << id.subdetId() << "  layer " << theTrkTopo->layer(id) << " valid "
              << hit->isValid() << "   detid: " << id() << std::endl;
      } else {
        PRINT << "  skipping subdet " << thisSubdet << "  layer " << theTrkTopo->layer(id) << " valid "
              << hit->isValid() << "   detid: " << id() << std::endl;
      }
    }
#else
    addSecondTrackHits(hits, outer);
#endif
  } else if (duplicateType == DuplicateTrackType::Overlapping) {
    addSecondTrackHits(hits, outer);
  }
 
  math::XYZVector p = (inner.innerMomentum() + outer.outerMomentum());
  GlobalVector v(p.x(), p.y(), p.z());
  if (!useInnermostState_)
    v *= -1;
 
  TrackCandidate::RecHitContainer ownHits;
  unsigned int nhits = hits.size();
  ownHits.reserve(nhits);
 
  if (duplicateType == DuplicateTrackType::Disjoint) {
#if TRACK_SORT
    if (!useInnermostState_)
      std::reverse(hits.begin(), hits.end());
    for (auto hit : hits)
      ownHits.push_back(*hit);
#elif DET_SORT
    // OLD METHOD, sometimes fails
    std::sort(hits.begin(), hits.end(), MomentumSort(v, &*theGeometry));
    for (auto hit : hits)
      ownHits.push_back(*hit);
#else
    sortByHitPosition(v, hits, ownHits);
#endif
  } else if (duplicateType == DuplicateTrackType::Overlapping) {
    sortByHitPosition(v, hits, ownHits);
  }
 
  PTrajectoryStateOnDet state;
  PropagationDirection pdir;
  if (useInnermostState_) {
    pdir = alongMomentum;
    if ((inner.outerPosition() - inner.innerPosition()).Dot(inner.momentum()) >= 0) {
      // use inner state
      TrajectoryStateOnSurface originalTsosIn(
          trajectoryStateTransform::innerStateOnSurface(inner, *theGeometry, &*theMagField));
      state = trajectoryStateTransform::persistentState(originalTsosIn, DetId(inner.innerDetId()));
    } else {
      // use outer state
      TrajectoryStateOnSurface originalTsosOut(
          trajectoryStateTransform::outerStateOnSurface(inner, *theGeometry, &*theMagField));
      state = trajectoryStateTransform::persistentState(originalTsosOut, DetId(inner.outerDetId()));
    }
  } else {
    pdir = oppositeToMomentum;
    if ((outer.outerPosition() - inner.innerPosition()).Dot(inner.momentum()) >= 0) {
      // use outer state
      TrajectoryStateOnSurface originalTsosOut(
          trajectoryStateTransform::outerStateOnSurface(outer, *theGeometry, &*theMagField));
      state = trajectoryStateTransform::persistentState(originalTsosOut, DetId(outer.outerDetId()));
    } else {
      // use inner state
      TrajectoryStateOnSurface originalTsosIn(
          trajectoryStateTransform::innerStateOnSurface(outer, *theGeometry, &*theMagField));
      state = trajectoryStateTransform::persistentState(originalTsosIn, DetId(outer.innerDetId()));
    }
  }
  TrajectorySeed seed(state, TrackCandidate::RecHitContainer(), pdir);
  TrackCandidate ret(ownHits, seed, state, (useInnermostState_ ? inner : outer).seedRef());
  ret.setStopReason((uint8_t)(useInnermostState_ ? inner : outer).stopReason());
  return ret;
}

References funct::abs(), addSecondTrackHits(), alongMomentum, debug_, Disjoint, DPRINT, hfClusterShapes_cfi::hits, triggerObjects_cff::id, SurfaceOrientation::inner, trajectoryStateTransform::innerStateOnSurface(), TrackerTopology::layer(), nhits, oppositeToMomentum, SurfaceOrientation::outer, trajectoryStateTransform::outerStateOnSurface(), Overlapping, AlCaHLTBitMon_ParallelJobs::p, trajectoryStateTransform::persistentState(), PRINT, edm::OwnVector< T, P >::push_back(), edm::OwnVector< T, P >::reserve(), runTheMatrix::ret, groupFilesInBlocks::reverse, fileCollector::seed, sortByHitPosition(), DetId::subdetId(), theGeometry, theMagField, theTrkTopo, useInnermostState_, and findQualityFiles::v.

sortByHitPosition()

void TrackMerger::sortByHitPosition ( const GlobalVector &  v, const std::vector< const TrackingRecHit * > &  hits, TrackCandidate::RecHitContainer &  ownHits  ) const

private

Definition at line 187 of file TrackMerger.cc.

                                                                                  {
  // NEW sort, more accurate
  unsigned int nhits = hits.size();
  std::vector<TransientTrackingRecHit::RecHitPointer> ttrh(nhits);
  for (unsigned int i = 0; i < nhits; ++i)
    ttrh[i] = theBuilder->build(hits[i]);
  std::sort(ttrh.begin(), ttrh.end(), GlobalMomentumSort(v));
  for (const auto &hit : ttrh)
    ownHits.push_back(*hit);
}

References TransientTrackingRecHitBuilder::build(), hfClusterShapes_cfi::hits, mps_fire::i, nhits, edm::OwnVector< T, P >::push_back(), theBuilder, and findQualityFiles::v.

Referenced by merge().

Member Data Documentation

builderToken_

edm::ESGetToken<TransientTrackingRecHitBuilder, TransientRecHitRecord> TrackMerger::builderToken_

private

Definition at line 38 of file TrackMerger.h.

Referenced by init().

debug_

bool TrackMerger::debug_

private

Definition at line 32 of file TrackMerger.h.

Referenced by merge().

geometryToken_

edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> TrackMerger::geometryToken_

private

Definition at line 36 of file TrackMerger.h.

Referenced by init().

magFieldToken_

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> TrackMerger::magFieldToken_

private

Definition at line 37 of file TrackMerger.h.

Referenced by init().

theBuilder

edm::ESHandle<TransientTrackingRecHitBuilder> TrackMerger::theBuilder

private

Definition at line 34 of file TrackMerger.h.

Referenced by init(), and sortByHitPosition().

theBuilderName

std::string TrackMerger::theBuilderName

private

Definition at line 33 of file TrackMerger.h.

theGeometry

edm::ESHandle<TrackerGeometry> TrackMerger::theGeometry

private

Definition at line 29 of file TrackMerger.h.

Referenced by init(), and merge().

theMagField

edm::ESHandle<MagneticField> TrackMerger::theMagField

private

Definition at line 30 of file TrackMerger.h.

Referenced by init(), and merge().

theTrkTopo

edm::ESHandle<TrackerTopology> TrackMerger::theTrkTopo

private

Definition at line 35 of file TrackMerger.h.

Referenced by addSecondTrackHits(), init(), and merge().

trackerTopoToken_

edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> TrackMerger::trackerTopoToken_

private

Definition at line 39 of file TrackMerger.h.

Referenced by init().

useInnermostState_

bool TrackMerger::useInnermostState_

private

Definition at line 31 of file TrackMerger.h.

Referenced by merge().