TrackTransformer Class Reference

#include <TrackTransformer.h>

Inheritance diagram for TrackTransformer:

Public Member Functions
TransientTrackingRecHit::ConstRecHitContainer	getTransientRecHits (const reco::TransientTrack &track) const
const MagneticField *	magneticField () const
	the magnetic field More...
std::unique_ptr < TrajectoryFitter > const &	refitter () const
	the refitter used to refit the reco::Track More...
virtual void	setServices (const edm::EventSetup &)
	set the services needed by the TrackTransformer More...
std::unique_ptr < TrajectorySmoother > const &	smoother () const
	the smoother used to smooth the trajectory which came from the refitting step More...
edm::ESHandle < GlobalTrackingGeometry >	trackingGeometry () const
	the tracking geometry More...
	TrackTransformer (const edm::ParameterSet &)
	Constructor. More...
virtual std::vector< Trajectory >	transform (const reco::Track &) const
	Convert a reco::Track into Trajectory. More...
std::vector< Trajectory >	transform (const reco::TrackRef &) const
	Convert a reco::TrackRef into Trajectory. More...
std::vector< Trajectory >	transform (const reco::TransientTrack &, const TransientTrackingRecHit::ConstRecHitContainer &) const
	Convert a reco::TrackRef into Trajectory, refit with a new set of hits. More...
virtual	~TrackTransformer ()
	Destructor. More...
Public Member Functions inherited from TrackTransformerBase
	TrackTransformerBase ()
	Constructor. More...
virtual	~TrackTransformerBase ()
	Destructor. More...

Private Member Functions
RefitDirection::GeometricalDirection	checkRecHitsOrdering (TransientTrackingRecHit::ConstRecHitContainer &) const
edm::ESHandle< Propagator >	propagator () const

Private Attributes
TkClonerImpl	hitCloner
unsigned long long	theCacheId_GTG
unsigned long long	theCacheId_MG
unsigned long long	theCacheId_TC
unsigned long long	theCacheId_TRH
bool	theDoPredictionsOnly
std::unique_ptr< TrajectoryFitter >	theFitter
std::string	theFitterName
edm::ESHandle< MagneticField >	theMGField
edm::ESHandle < TransientTrackingRecHitBuilder >	theMuonRecHitBuilder
std::string	theMuonRecHitBuilderName
edm::ESHandle< Propagator >	thePropagator
std::string	thePropagatorName
RefitDirection	theRefitDirection
bool	theRPCInTheFit
std::unique_ptr < TrajectorySmoother >	theSmoother
std::string	theSmootherName
edm::ESHandle < TransientTrackingRecHitBuilder >	theTrackerRecHitBuilder
std::string	theTrackerRecHitBuilderName
edm::ESHandle < GlobalTrackingGeometry >	theTrackingGeometry

Detailed Description

This class takes a reco::Track and refits the rechits inside it. The final result is a Trajectory refitted and smoothed. To make the refitting (and the smoothing) the usual KF tools are used.

CAVEAT: till now (it will be changed in the near future) the class stores the pointers to the services, therefore EACH event the setServices(const edm::EventSetup&) method MUST be called in the code in which the TrackTransformer is used.

Author

R. Bellan - INFN Torino ricca.nosp@m.rdo..nosp@m.bella.nosp@m.n@ce.nosp@m.rn.ch

This class takes a reco::Track and refits the rechits inside it. The final result is a Trajectory refitted and smoothed. To make the refitting (and the smoothing) the usual KF tools are used.

CAVEAT: till now (it will be changed in the near future) the class stores the pointers to the services, therefore EACH event the setServices(const edm::EventSetup&) method MUST be called in the code in which the TrackTransformer is used.

Rec hits are ordered in rank of increasing |z| assuing muons are coming from the top of the detector.

$Date: 2009/01/15 author R. Bellan - CERN ricca.nosp@m.rdo..nosp@m.bella.nosp@m.n@ce.nosp@m.rn.ch modified for zed ordering by N. Kypreos - UF nicho.nosp@m.las..nosp@m.theod.nosp@m.ore..nosp@m.kypre.nosp@m.os@c.nosp@m.ern.c.nosp@m.h

This class takes a reco::Track and refits the rechits inside it. The final result is a Trajectory refitted and smoothed. To make the refitting (and the smoothing) the usual KF tools are used.

CAVEAT: till now (it will be changed in the near future) the class stores the pointers to the services, therefore EACH event the setServices(const edm::EventSetup&) method MUST be called in the code in which the TrackTransformer is used.

Author

R. Bellan - CERN ricca.nosp@m.rdo..nosp@m.bella.nosp@m.n@ce.nosp@m.rn.ch

Definition at line 38 of file TrackTransformer.h.

Constructor & Destructor Documentation

TrackTransformer::TrackTransformer

(

const edm::ParameterSet &

parameterSet

)

Constructor.

Definition at line 29 of file TrackTransformer.cc.

References edm::ParameterSet::getParameter().

                                                                  {
  
  // Refit direction
  string refitDirectionName = parameterSet.getParameter<string>("RefitDirection");
  theRefitDirection = RefitDirection(refitDirectionName);
  
  theFitterName = parameterSet.getParameter<string>("Fitter");  
  theSmootherName = parameterSet.getParameter<string>("Smoother");  
  thePropagatorName = parameterSet.getParameter<string>("Propagator");

  theTrackerRecHitBuilderName = parameterSet.getParameter<string>("TrackerRecHitBuilder");
  theMuonRecHitBuilderName = parameterSet.getParameter<string>("MuonRecHitBuilder");

  theRPCInTheFit = parameterSet.getParameter<bool>("RefitRPCHits");
  theDoPredictionsOnly = parameterSet.getParameter<bool>("DoPredictionsOnly");

  theCacheId_TC = theCacheId_GTG = theCacheId_MG = theCacheId_TRH = 0;
}

TrackTransformer::~TrackTransformer ( )

virtual

Destructor.

Definition at line 49 of file TrackTransformer.cc.

{}

Member Function Documentation

RefitDirection::GeometricalDirection TrackTransformer::checkRecHitsOrdering ( TransientTrackingRecHit::ConstRecHitContainer &  recHits ) const

private

Definition at line 133 of file TrackTransformer.cc.

References first, RefitDirection::insideOut, prof2calltree::last, LogDebug, PV3DBase< T, PVType, FrameType >::mag(), RefitDirection::outsideIn, position, and RefitDirection::undetermined.

                                                                                                 {
  
  if (!recHits.empty()){
    GlobalPoint first = trackingGeometry()->idToDet(recHits.front()->geographicalId())->position();
    GlobalPoint last = trackingGeometry()->idToDet(recHits.back()->geographicalId())->position();
    
    double rFirst = first.mag();
    double rLast  = last.mag();
    if(rFirst < rLast) return RefitDirection::insideOut;
    else if(rFirst > rLast) return RefitDirection::outsideIn;
    else{
      LogDebug("Reco|TrackingTools|TrackTransformer") << "Impossible to determine the rechits order" <<endl;
      return RefitDirection::undetermined;
    }
  }
  else{
    LogDebug("Reco|TrackingTools|TrackTransformer") << "Impossible to determine the rechits order" <<endl;
    return RefitDirection::undetermined;
  }
}

TransientTrackingRecHit::ConstRecHitContainer TrackTransformer::getTransientRecHits

(

const reco::TransientTrack &

track

)

const

Definition at line 108 of file TrackTransformer.cc.

References LogTrace, DetId::Muon, reco::TransientTrack::recHitsBegin(), reco::TransientTrack::recHitsEnd(), query::result, and DetId::Tracker.

                                                                           {

  TransientTrackingRecHit::ConstRecHitContainer result;
  auto tkbuilder = static_cast<TkTransientTrackingRecHitBuilder const *>(theTrackerRecHitBuilder.product());

  
  for (trackingRecHit_iterator hit = track.recHitsBegin(); hit != track.recHitsEnd(); ++hit) {
    if((*hit)->isValid()) {
      if ( (*hit)->geographicalId().det() == DetId::Tracker ) {
    result.emplace_back((**hit).cloneForFit(*tkbuilder->geometry()->idToDet( (**hit).geographicalId() ) ) );
      } else if ( (*hit)->geographicalId().det() == DetId::Muon ){
    if( (*hit)->geographicalId().subdetId() == 3 && !theRPCInTheFit){
      LogTrace("Reco|TrackingTools|TrackTransformer") << "RPC Rec Hit discarged"; 
      continue;
    }
    result.push_back(theMuonRecHitBuilder->build(&**hit));
      }
    }
  }
  
  return result;
}

const MagneticField* TrackTransformer::magneticField ( ) const

inline

the magnetic field

Definition at line 61 of file TrackTransformer.h.

References theMGField.

{return &*theMGField;}

edm::ESHandle<Propagator> TrackTransformer::propagator ( void  ) const

inlineprivate

Definition at line 83 of file TrackTransformer.h.

References thePropagator.

{return thePropagator;}

std::unique_ptr<TrajectoryFitter> const& TrackTransformer::refitter ( ) const

inline

the refitter used to refit the reco::Track

Definition at line 70 of file TrackTransformer.h.

References theFitter.

{return theFitter;}

void TrackTransformer::setServices ( const edm::EventSetup &  setup )

virtual

set the services needed by the TrackTransformer

Implements TrackTransformerBase.

Definition at line 52 of file TrackTransformer.cc.

References edm::EventSetup::get(), LogTrace, metname, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by MuonKinkFinder::init(), MuonReSeeder::produce(), TrackerToMuonPropagator::produce(), CSCOverlapsAlignmentAlgorithm::run(), and GlobalTrajectoryBuilderBase::setEvent().

                                                         {
  
  const std::string metname = "Reco|TrackingTools|TrackTransformer";

  edm::ESHandle<TrajectoryFitter> aFitter;
  edm::ESHandle<TrajectorySmoother> aSmoother;
  setup.get<TrajectoryFitter::Record>().get(theFitterName,aFitter);
  setup.get<TrajectoryFitter::Record>().get(theSmootherName,aSmoother);
  theFitter = aFitter->clone();
  theSmoother.reset(aSmoother->clone());
  
  unsigned long long newCacheId_TC = setup.get<TrackingComponentsRecord>().cacheIdentifier();

  if ( newCacheId_TC != theCacheId_TC ){
    LogTrace(metname) << "Tracking Component changed!";
    theCacheId_TC = newCacheId_TC;
    setup.get<TrackingComponentsRecord>().get(thePropagatorName,thePropagator);
  }

  // Global Tracking Geometry
  unsigned long long newCacheId_GTG = setup.get<GlobalTrackingGeometryRecord>().cacheIdentifier();
  if ( newCacheId_GTG != theCacheId_GTG ) {
    LogTrace(metname) << "GlobalTrackingGeometry changed!";
    theCacheId_GTG = newCacheId_GTG;
    setup.get<GlobalTrackingGeometryRecord>().get(theTrackingGeometry); 
  }
  
  // Magfield Field
  unsigned long long newCacheId_MG = setup.get<IdealMagneticFieldRecord>().cacheIdentifier();
  if ( newCacheId_MG != theCacheId_MG ) {
    LogTrace(metname) << "Magnetic Field changed!";
    theCacheId_MG = newCacheId_MG;
    setup.get<IdealMagneticFieldRecord>().get(theMGField);
  }
  
  // Transient Rechit Builders
  unsigned long long newCacheId_TRH = setup.get<TransientRecHitRecord>().cacheIdentifier();
  if ( newCacheId_TRH != theCacheId_TRH ) {
    theCacheId_TRH = newCacheId_TRH;
    LogTrace(metname) << "TransientRecHitRecord changed!";
    setup.get<TransientRecHitRecord>().get(theTrackerRecHitBuilderName,theTrackerRecHitBuilder);
    setup.get<TransientRecHitRecord>().get(theMuonRecHitBuilderName,theMuonRecHitBuilder);
    hitCloner = static_cast<TkTransientTrackingRecHitBuilder const *>(theTrackerRecHitBuilder.product())->cloner();
  }
  theFitter->setHitCloner(&hitCloner);
  theSmoother->setHitCloner(&hitCloner);

}

std::unique_ptr<TrajectorySmoother> const& TrackTransformer::smoother ( ) const

inline

the smoother used to smooth the trajectory which came from the refitting step

Definition at line 73 of file TrackTransformer.h.

References theSmoother.

{return theSmoother;}

edm::ESHandle<GlobalTrackingGeometry> TrackTransformer::trackingGeometry ( ) const

inline

the tracking geometry

Definition at line 64 of file TrackTransformer.h.

References theTrackingGeometry.

{return theTrackingGeometry;}

vector< Trajectory > TrackTransformer::transform ( const reco::Track &  newTrack ) const

virtual

Convert a reco::Track into Trajectory.

Convert Tracks into Trajectories.

Implements TrackTransformerBase.

Definition at line 165 of file TrackTransformer.cc.

References metname, AlCaHLTBitMon_QueryRunRegistry::string, and create_public_lumi_plots::transform.

Referenced by GlobalTrajectoryBuilderBase::build(), CSCPairResidualsConstraint::dphidzFromTrack(), MuonKinkFinder::fillTrkKink(), MuonReSeeder::produce(), and TrackerToMuonPropagator::produce().

                                                                              {

  const std::string metname = "Reco|TrackingTools|TrackTransformer";
  
  reco::TransientTrack track(newTrack,magneticField(),trackingGeometry());   

  // Build the transient Rechits
  TransientTrackingRecHit::ConstRecHitContainer recHitsForReFit = getTransientRecHits(track);

  return transform(track, recHitsForReFit);
}

vector< Trajectory > TrackTransformer::transform

(

const reco::TrackRef &

track

)

const

Convert a reco::TrackRef into Trajectory.

Definition at line 102 of file TrackTransformer.cc.

References create_public_lumi_plots::transform.

                                                                              {
  return transform(*track);
}

vector< Trajectory > TrackTransformer::transform	(	const reco::TransientTrack &	track,
		const TransientTrackingRecHit::ConstRecHitContainer &	_recHitsForReFit
	)		const

Convert a reco::TrackRef into Trajectory, refit with a new set of hits.

Convert Tracks into Trajectories with a given set of hits.

Definition at line 179 of file TrackTransformer.cc.

References alongMomentum, anyDirection, TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), reco::Track::innerDetId(), reco::TransientTrack::innermostMeasurementState(), RefitDirection::insideOut, TrajectoryStateOnSurface::isValid(), LogTrace, metname, oppositeToMomentum, reco::Track::outerDetId(), reco::TransientTrack::outermostMeasurementState(), RefitDirection::outsideIn, AlCaHLTBitMon_ParallelJobs::p, position, Trajectory::push(), fileCollector::seed, AlCaHLTBitMon_QueryRunRegistry::string, TrajectoryStateOnSurface::surface(), reco::TransientTrack::track(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

                                                                                                                          {
  
  TransientTrackingRecHit::ConstRecHitContainer recHitsForReFit =  _recHitsForReFit;
  const std::string metname = "Reco|TrackingTools|TrackTransformer";

  if(recHitsForReFit.size() < 2) return vector<Trajectory>();
 
  // 8 cases are foreseen: 
  // [RH = rec hit order, P = momentum dir, FD = fit direction. IO/OI = inside-out/outside-in, AM/OM = along momentum/opposite to momentum]
  // (1) RH IO | P IO | FD AM  ---> Start from IN
  // (2) RH IO | P IO | FD OM  ---> Reverse RH and start from OUT
  // (3) RH IO | P OI | FD AM  ---> Reverse RH and start from IN
  // (4) RH IO | P OI | FD OM  ---> Start from OUT
  // (5) RH OI | P IO | FD AM  ---> Reverse RH and start from IN
  // (6) RH OI | P IO | FD OM  ---> Start from OUT
  // (7) RH OI | P OI | FD AM  ---> Start from IN
  // (8) RH OI | P OI | FD OM  ---> Reverse RH and start from OUT
  //
  // *** Rules: ***
  // -A- If RH-FD agree (IO-AM,OI-OM) do not reverse the RH
  // -B- If FD along momentum start from innermost state, otherwise use outermost
  
  // Other special cases can be handled:
  // (1 bis) RH IO | P IO | GFD IO => FD AM  ---> Start from IN
  // (2 bis) RH IO | P IO | GFD OI => FD OM  ---> Reverse RH and start from OUT
  // (3 bis) RH IO | P OI | GFD OI => FD AM  ---> Reverse RH and start from OUT
  // (4 bis) RH IO | P OI | GFD IO => FD OM  ---> Start from IN
  // (5 bis) RH OI | P IO | GFD IO => FD AM  ---> Reverse RH and start from IN
  // (6 bis) RH OI | P IO | GFD OI => FD OM  ---> Start from OUT
  // (7 bis) RH OI | P OI | GFD OI => FD AM  ---> Start from OUT
  // (8 bis) RH OI | P OI | GFD IO => FD OM  ---> Reverse RH and start from IN
  // 
  // *** Additional rule: ***
  // -A0- If P and GFD agree, then FD is AM otherwise is OM
  // -A00- rechit must be ordered as GFD in order to handle the case of cosmics
  // -B0- The starting state is decided by GFD

  // Determine the RH order
  RefitDirection::GeometricalDirection recHitsOrder = checkRecHitsOrdering(recHitsForReFit); // FIXME change nome of the *type*  --> RecHit order!
  LogTrace(metname) << "RH order (0-insideOut, 1-outsideIn): " << recHitsOrder;

  PropagationDirection propagationDirection = theRefitDirection.propagationDirection();

  // Apply rule -A0-
  if(propagationDirection == anyDirection){
    GlobalVector momentum = track.innermostMeasurementState().globalMomentum();
    GlobalVector position = track.innermostMeasurementState().globalPosition() - GlobalPoint(0,0,0);
    RefitDirection::GeometricalDirection p = (momentum.x()*position.x() > 0 || momentum.y()*position.y() > 0) ? RefitDirection::insideOut : RefitDirection::outsideIn;

    propagationDirection = p == theRefitDirection.geometricalDirection() ? alongMomentum : oppositeToMomentum;
    LogTrace(metname) << "P  (0-insideOut, 1-outsideIn): " << p;
    LogTrace(metname) << "FD (0-OM, 1-AM, 2-ANY): " << propagationDirection;
  }
  // -A0-

  // Apply rule -A-
  if(theRefitDirection.propagationDirection() != anyDirection){
    if((recHitsOrder == RefitDirection::insideOut && propagationDirection == oppositeToMomentum) ||
       (recHitsOrder == RefitDirection::outsideIn && propagationDirection == alongMomentum) ) 
      reverse(recHitsForReFit.begin(),recHitsForReFit.end());}
  // -A-
  // Apply rule -A00-
  else{
    // reorder the rechit as defined in theRefitDirection.geometricalDirection(); 
    if(theRefitDirection.geometricalDirection() != recHitsOrder) reverse(recHitsForReFit.begin(),recHitsForReFit.end()); 
  }
  // -A00-
    
  // Apply rule -B-
  TrajectoryStateOnSurface firstTSOS = track.innermostMeasurementState();
  unsigned int innerId = track.track().innerDetId();
  if(theRefitDirection.propagationDirection() != anyDirection){
    if(propagationDirection == oppositeToMomentum){
      innerId   = track.track().outerDetId();
      firstTSOS = track.outermostMeasurementState();
    }
  }
  else { // if(theRefitDirection.propagationDirection() == anyDirection)
    // Apply rule -B0-
    if(theRefitDirection.geometricalDirection() == RefitDirection::outsideIn){
      innerId   = track.track().outerDetId();
      firstTSOS = track.outermostMeasurementState();
    }
    // -B0-
  }
  // -B-

  if(!firstTSOS.isValid()){
    LogTrace(metname)<<"Error wrong initial state!"<<endl;
    return vector<Trajectory>();
  }

  TrajectorySeed seed(PTrajectoryStateOnDet(),TrajectorySeed::recHitContainer(),propagationDirection);

  if(recHitsForReFit.front()->geographicalId() != DetId(innerId)){
    LogTrace(metname)<<"Propagation occured"<<endl;
    firstTSOS = propagator()->propagate(firstTSOS, recHitsForReFit.front()->det()->surface());
    if(!firstTSOS.isValid()){
      LogTrace(metname)<<"Propagation error!"<<endl;
      return vector<Trajectory>();
    }
  }

  if(theDoPredictionsOnly){
    Trajectory aTraj(seed,propagationDirection);
    TrajectoryStateOnSurface predTSOS = firstTSOS;
    for(TransientTrackingRecHit::ConstRecHitContainer::const_iterator ihit = recHitsForReFit.begin(); 
    ihit != recHitsForReFit.end(); ++ihit ) {
      predTSOS = propagator()->propagate(predTSOS, (*ihit)->det()->surface());
      if (predTSOS.isValid()) aTraj.push(TrajectoryMeasurement(predTSOS, *ihit));
    }
    return vector<Trajectory>(1, aTraj);
  }


  vector<Trajectory> trajectories = theFitter->fit(seed,recHitsForReFit,firstTSOS);
  
  if(trajectories.empty()){
    LogTrace(metname)<<"No Track refitted!"<<endl;
    return vector<Trajectory>();
  }
  
  Trajectory trajectoryBW = trajectories.front();
    
  vector<Trajectory> trajectoriesSM = theSmoother->trajectories(trajectoryBW);

  if(trajectoriesSM.empty()){
    LogTrace(metname)<<"No Track smoothed!"<<endl;
    return vector<Trajectory>();
  }
  
  return trajectoriesSM;

}

Member Data Documentation

TkClonerImpl TrackTransformer::hitCloner

private

Definition at line 112 of file TrackTransformer.h.

unsigned long long TrackTransformer::theCacheId_GTG

private

Definition at line 87 of file TrackTransformer.h.

unsigned long long TrackTransformer::theCacheId_MG

private

Definition at line 88 of file TrackTransformer.h.

unsigned long long TrackTransformer::theCacheId_TC

private

Definition at line 86 of file TrackTransformer.h.

unsigned long long TrackTransformer::theCacheId_TRH

private

Definition at line 89 of file TrackTransformer.h.

bool TrackTransformer::theDoPredictionsOnly

private

Definition at line 93 of file TrackTransformer.h.

std::unique_ptr<TrajectoryFitter> TrackTransformer::theFitter

private

Definition at line 100 of file TrackTransformer.h.

Referenced by refitter().

std::string TrackTransformer::theFitterName

private

Definition at line 99 of file TrackTransformer.h.

edm::ESHandle<MagneticField> TrackTransformer::theMGField

private

Definition at line 97 of file TrackTransformer.h.

Referenced by magneticField().

edm::ESHandle<TransientTrackingRecHitBuilder> TrackTransformer::theMuonRecHitBuilder

private

Definition at line 115 of file TrackTransformer.h.

std::string TrackTransformer::theMuonRecHitBuilderName

private

Definition at line 114 of file TrackTransformer.h.

edm::ESHandle<Propagator> TrackTransformer::thePropagator

private

Definition at line 84 of file TrackTransformer.h.

Referenced by propagator().

std::string TrackTransformer::thePropagatorName

private

Definition at line 82 of file TrackTransformer.h.

RefitDirection TrackTransformer::theRefitDirection

private

Definition at line 94 of file TrackTransformer.h.

bool TrackTransformer::theRPCInTheFit

private

Definition at line 91 of file TrackTransformer.h.

std::unique_ptr<TrajectorySmoother> TrackTransformer::theSmoother

private

Definition at line 103 of file TrackTransformer.h.

Referenced by smoother().

std::string TrackTransformer::theSmootherName

private

Definition at line 102 of file TrackTransformer.h.

edm::ESHandle<TransientTrackingRecHitBuilder> TrackTransformer::theTrackerRecHitBuilder

private

Definition at line 111 of file TrackTransformer.h.

std::string TrackTransformer::theTrackerRecHitBuilderName

private

Definition at line 110 of file TrackTransformer.h.

edm::ESHandle<GlobalTrackingGeometry> TrackTransformer::theTrackingGeometry

private

Definition at line 96 of file TrackTransformer.h.

Referenced by trackingGeometry().