Trajectory Class Reference

#include <Trajectory.h>

Public Types
using	ConstRecHitContainer = TrackingRecHit::ConstRecHitContainer
typedef std::vector< TrajectoryMeasurement >	DataContainer
using	RecHitContainer = ConstRecHitContainer

Public Member Functions
int	cccBadHits () const
float	cccThreshold () const
float	chiSquared () const
TrajectoryMeasurement const &	closestMeasurement (GlobalPoint) const
DataContainer const &	data () const
	obsolete name, use measurements() instead. More...
PropagationDirection const &	direction () const
float	dPhiCacheForLoopersReconstruction () const
bool	empty () const
	True if trajectory has no measurements. More...
TrajectoryMeasurement const &	firstMeasurement () const
int	foundHits () const
int	foundPixelHits () const
TrajectoryStateOnSurface	geometricalInnermostState () const
void	incrementLoops ()
void	invalidate ()
	Method to invalidate a trajectory. Useful during ambiguity resolution. More...
bool	isLooper () const
bool	isValid () const
const DetLayer *	lastLayer () const
	Redundant method, returns the layer of lastMeasurement() . More...
TrajectoryMeasurement const &	lastMeasurement () const
int	lostHits () const
DataContainer &	measurements ()
DataContainer const &	measurements () const
int	ndof (bool bon=true) const
signed char	nLoops () const
int	numberOfCCCBadHits (float ccc_threshold)
Trajectory &	operator= (Trajectory &&rh)
Trajectory &	operator= (Trajectory const &rh)=default
void	pop ()
void	push (const TrajectoryMeasurement &tm)
void	push (const TrajectoryMeasurement &tm, double chi2Increment)
void	push (TrajectoryMeasurement &&tm)
void	push (TrajectoryMeasurement &&tm, double chi2Increment)
ConstRecHitContainer	recHits () const
void	reserve (unsigned int n)
void	reverse ()
TrajectorySeed const &	seed () const
	Access to the seed used to reconstruct the Trajectory. More...
unsigned int	seedNHits () const
edm::RefToBase< TrajectorySeed >	seedRef (void) const
void	setDPhiCacheForLoopersReconstruction (float dphi)
void	setNLoops (signed char value)
void	setSeedRef (const edm::RefToBase< TrajectorySeed > &seedRef)
void	setSharedSeed (const std::shared_ptr< const TrajectorySeed > &seed)
void	setStopReason (StopReason s)
const std::shared_ptr< const TrajectorySeed > &	sharedSeed () const
StopReason	stopReason () const
int	trailingFoundHits () const
	Trajectory ()
	Trajectory (const std::shared_ptr< const TrajectorySeed > &seed, PropagationDirection dir)
	Trajectory (const TrajectorySeed &seed)
	Trajectory (const TrajectorySeed &seed, PropagationDirection dir)
	Trajectory (PropagationDirection dir)
	Trajectory (Trajectory &&rh)
	Trajectory (Trajectory const &rh)=default
void	validRecHits (ConstRecHitContainer &cont) const

Static Public Member Functions
static bool	inactive ()
static bool	isBad (const TrackingRecHit &hit)
static bool	lost (const TrackingRecHit &hit)
static bool	pixel (const TrackingRecHit &hit)

Private Member Functions
bool	badForCCC (const TrajectoryMeasurement &tm)
void	check () const
void	pushAux (double chi2Increment)
void	updateBadForCCC (float ccc_threshold)

Private Attributes
edm::RefToBase< TrajectorySeed >	seedRef_
StopReason	stopReason_ = StopReason::UNINITIALIZED
float	theCCCThreshold_ = std::numeric_limits<float>::max()
float	theChiSquared = 0
float	theChiSquaredBad = 0
DataContainer	theData
PropagationDirection	theDirection = anyDirection
bool	theDirectionValidity = false
float	theDPhiCache = 0
signed char	theNLoops = 0
signed short	theNumberOfCCCBadHits_ = 0
signed short	theNumberOfFoundHits = 0
signed short	theNumberOfFoundPixelHits = 0
signed short	theNumberOfLostHits = 0
signed short	theNumberOfTrailingFoundHits = 0
std::shared_ptr< const TrajectorySeed >	theSeed
bool	theValid = false

Detailed Description

A class for detailed particle trajectory representation. It is used during trajectory building to "grow" a trajectory. The trajectory is represented as an ordered sequence of TrajectoryMeasurement objects with a stack-like interface. The measurements are added to the Trajectory in the order of increasing precision: each new TrajectoryMeasurement is assumed to improve the precision of the last one, normally by adding a constraint from a new RecHit. However the Trajectory class does not have the means to verify that measurements are added in the correct order, and thus cannot guarantee the order, which is the responsibility of the TrajectoryBuilder. The Trajectory provides some security by allowing to add or remove measurements only on one of it's ends, with push(TM) and pop() methods. The last measurement in a Trajectory can thus be either the innermost (closest to the interaction point) or the outermost, depending on the way the Trajectory was built. The direction of building is represented as a PropagationDirection, which has two possible values: alongMomentum (outwards) and oppositeToMomentum (inwards), and is accessed with the direction() method.

Definition at line 38 of file Trajectory.h.

Member Typedef Documentation

ConstRecHitContainer

using Trajectory::ConstRecHitContainer = TrackingRecHit::ConstRecHitContainer

Definition at line 41 of file Trajectory.h.

DataContainer

typedef std::vector<TrajectoryMeasurement> Trajectory::DataContainer

Definition at line 40 of file Trajectory.h.

RecHitContainer

using Trajectory::RecHitContainer = ConstRecHitContainer

Definition at line 42 of file Trajectory.h.

Constructor & Destructor Documentation

Trajectory() [1/7]

Trajectory::Trajectory ( )

inline

Default constructor of an empty trajectory with undefined seed and undefined direction. This constructor is necessary in order to transiently copy vector<Trajectory> in the edm::Event

Definition at line 49 of file Trajectory.h.

{}

Trajectory() [2/7]

Trajectory::Trajectory ( const TrajectorySeed &  seed )

inlineexplicit

Constructor of an empty trajectory with undefined direction. The direction will be defined at the moment of the push of a second measurement, from the relative radii of the first and second measurements.

Definition at line 57 of file Trajectory.h.

: theSeed(new TrajectorySeed(seed)), theValid(true) {}

Trajectory() [3/7]

Trajectory::Trajectory ( const TrajectorySeed &  seed, PropagationDirection  dir  )

inline

Constructor of an empty trajectory with defined direction. No check is made in the push method that measurements are added in the correct direction.

Definition at line 63 of file Trajectory.h.

      : theSeed(new TrajectorySeed(seed)), theDirection(dir), theDirectionValidity(true), theValid(true) {}

Trajectory() [4/7]

Trajectory::Trajectory ( const std::shared_ptr< const TrajectorySeed > &  seed, PropagationDirection  dir  )

inline

Constructor of an empty trajectory with defined direction. No check is made in the push method that measurements are added in the correct direction.

Definition at line 70 of file Trajectory.h.

      : theSeed(seed), theDirection(dir), theDirectionValidity(true), theValid(true) {}

Trajectory() [5/7]

Trajectory::Trajectory ( PropagationDirection  dir )

inlineexplicit

Constructor of an empty trajectory with defined direction. No check is made in the push method that measurements are added in the correct direction.

Definition at line 77 of file Trajectory.h.

: theDirection(dir), theDirectionValidity(true), theValid(true) {}

Trajectory() [6/7]

Trajectory::Trajectory ( Trajectory const &  rh )

default

Trajectory() [7/7]

Trajectory::Trajectory ( Trajectory &&  rh )

inline

Definition at line 82 of file Trajectory.h.

      : theSeed(std::move(rh.theSeed)),
        seedRef_(std::move(rh.seedRef_)),
        theData(std::move(rh.theData)),
        theChiSquared(rh.theChiSquared),
        theChiSquaredBad(rh.theChiSquaredBad),
        theNumberOfFoundHits(rh.theNumberOfFoundHits),
        theNumberOfFoundPixelHits(rh.theNumberOfFoundPixelHits),
        theNumberOfLostHits(rh.theNumberOfLostHits),
        theNumberOfTrailingFoundHits(rh.theNumberOfTrailingFoundHits),
        theNumberOfCCCBadHits_(rh.theNumberOfCCCBadHits_),
        theDirection(rh.theDirection),
        theDirectionValidity(rh.theDirectionValidity),
        theValid(rh.theValid),
        theDPhiCache(rh.theDPhiCache),
        theCCCThreshold_(rh.theCCCThreshold_),
        theNLoops(rh.theNLoops),
        stopReason_(rh.stopReason_) {}

Member Function Documentation

badForCCC()

bool Trajectory::badForCCC ( const TrajectoryMeasurement &  tm )

private

Definition at line 182 of file Trajectory.cc.

                                                          {
  if (!trackerHitRTTI::isFromDet(*tm.recHit()))
    return false;
  auto const* thit = static_cast<const BaseTrackerRecHit*>(tm.recHit()->hit());
  if (!thit)
    return false;
  if (thit->isPixel() || thit->isPhase2())
    return false;
  if (!tm.updatedState().isValid())
    return false;
  return siStripClusterTools::chargePerCM(thit->rawId(),
                                          thit->firstClusterRef().stripCluster(),
                                          tm.updatedState().localParameters()) < theCCCThreshold_;
}

References siStripClusterTools::chargePerCM(), trackerHitRTTI::isFromDet(), TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localParameters(), TrajectoryMeasurement::recHit(), and TrajectoryMeasurement::updatedState().

cccBadHits()

int Trajectory::cccBadHits ( ) const

inline

Number of hits that are not compatible with the CCC used during patter recognition. Used mainly as a criteria for abandoning a trajectory candidate during trajectory building.

Definition at line 227 of file Trajectory.h.

{ return theNumberOfCCCBadHits_; }

References theNumberOfCCCBadHits_.

cccThreshold()

float Trajectory::cccThreshold ( ) const

inline

Definition at line 322 of file Trajectory.h.

{ return theCCCThreshold_; }

References theCCCThreshold_.

check()

void Trajectory::check ( void  ) const

private

Definition at line 140 of file Trajectory.cc.

                             {
  if (theData.empty())
    throw cms::Exception("TrackingTools/PatternTools",
                         "Trajectory::check() - information requested from empty Trajectory");
}

Referenced by firstMeasurement(), lastLayer(), and lastMeasurement().

chiSquared()

float Trajectory::chiSquared ( ) const

inline

• Trajectories with at least 1 valid hit: value of the raw Chi2 of the trajectory, not normalised to the N.D.F. (evaluated using only valid hits)
• Trajectories with only invalid hits: raw Chi2 (not norm.) of invalid hits w.r.t. the "default" trajectory (traj. containing only the seed information)

Definition at line 241 of file Trajectory.h.

{ return (theNumberOfFoundHits ? theChiSquared : theChiSquaredBad); }

References theChiSquared, theChiSquaredBad, and theNumberOfFoundHits.

Referenced by CosmicTrajectoryBuilder::AddHit(), CRackTrajectoryBuilder::AddHit(), SiStripGainFromData::algoAnalyze(), analyzer::SiPixelLorentzAngle::analyze(), OverlapValidation::analyze(), DAFTrackProducerAlgorithm::buildTrack(), TrackExtenderWithMTDT< TrackCollection >::buildTrack(), MuonTrackLoader::buildTrackAtPCA(), GlobalTrackerMuonAlignment::debugTrajectory(), AlignmentMonitorTracksFromTrajectories::event(), LowPtGsfElectronSeedProducer::lightGsfTracking(), GoodSeedProducer::produce(), and CosmicMuonTrajectoryBuilder::trajectories().

closestMeasurement()

TrajectoryMeasurement const & Trajectory::closestMeasurement

(

GlobalPoint

point

)

const

Definition at line 256 of file Trajectory.cc.

                                                                                   {
  check();
  vector<TrajectoryMeasurement>::const_iterator iter =
      std::min_element(measurements().begin(), measurements().end(), LessMag(point));
 
  return (*iter);
}

References begin, RPCNoise_example::check, end, and point.

Referenced by ConversionSeedFilter::getTSOS().

data()

DataContainer const& Trajectory::data ( ) const

inline

obsolete name, use measurements() instead.

Definition at line 182 of file Trajectory.h.

{ return measurements(); }

References measurements().

direction()

PropagationDirection const & Trajectory::direction

(

)

const

Direction of "growing" of the trajectory. Possible values are alongMomentum (outwards) and oppositeToMomentum (inwards).

Definition at line 133 of file Trajectory.cc.

                                                        {
  if (theDirectionValidity)
    return theDirection;
  else
    throw cms::Exception("TrackingTools/PatternTools", "Trajectory::direction() requested but not set");
}

References Exception.

Referenced by PFTrackTransformer::addPoints(), PFTrackTransformer::addPointsAndBrems(), CkfDebugger::analyseRecHitNotFound(), FWTrackProxyBuilderFullFramework::buildTrack(), DAFTrackProducerAlgorithm::buildTrack(), MuonTrackLoader::buildTrackExtra(), TrackExtenderWithMTDT< TrackCollection >::buildTrackExtra(), GlobalTrackerMuonAlignment::debugTrajectory(), KFSplittingFitter::fitOne(), NuclearTrackCorrector::getNewTrackExtra(), TransientInitialStateEstimator::innerState(), CRackTrajectoryBuilder::innerState(), Traj2TrackHits::operator()(), TrajectoryFactoryBase::orderedTrajectoryMeasurements(), GsfTrackProducerBase::putInEvt(), TrackProducerWithSCAssociation::putInEvt(), GroupedDAFHitCollector::recHits(), CosmicMuonTrajectoryBuilder::reverseTrajectoryPropagationDirection(), NuclearInteractionFinder::run(), DAFTrackProducerAlgorithm::runWithCandidate(), TrajAnnealing::TrajAnnealing(), CosmicMuonTrajectoryBuilder::trajectories(), GsfTrajectorySmoother::trajectory(), KFTrajectorySmoother::trajectory(), and trajectoryToResiduals().

dPhiCacheForLoopersReconstruction()

float Trajectory::dPhiCacheForLoopersReconstruction ( ) const

inline

accessor to the delta phi angle betweem the directions of the two measurements on the last two layers crossed by the trajectory

Definition at line 320 of file Trajectory.h.

{ return theDPhiCache; }

References theDPhiCache.

empty()

bool Trajectory::empty ( void  ) const

inline

True if trajectory has no measurements.

Definition at line 233 of file Trajectory.h.

{ return theData.empty(); }

References theData.

Referenced by BaseCkfTrajectoryBuilder::addToResult(), CosmicMuonTrajectoryBuilder::build(), CosmicMuonTrajectoryBuilder::buildSecondHalf(), GsfTrajectoryFitter::fitOne(), KFSplittingFitter::fitOne(), KFTrajectoryFitter::fitOne(), CosmicMuonTrajectoryBuilder::reverseTrajectoryPropagationDirection(), NuclearInteractionFinder::run(), CosmicMuonTrajectoryBuilder::selfDuplicate(), StandAloneMuonTrajectoryBuilder::trajectories(), CosmicMuonTrajectoryBuilder::trajectories(), GsfTrajectorySmoother::trajectory(), and KFTrajectorySmoother::trajectory().

firstMeasurement()

TrajectoryMeasurement const& Trajectory::firstMeasurement ( ) const

inline

Access to the first measurement. It is the least precise one in a trajectory before smoothing. It is precise in a smoothed trajectory. It's the innermost measurement if direction() == alongMomentum, the outermost one if direction() == oppositeToMomentum.

Definition at line 166 of file Trajectory.h.

                                                        {
    check();
    if (theData.front().recHitR().hit() != nullptr)
      return theData.front();
    else if (theData.size() > 2)
      return *(theData.begin() + 1);
    else
      throw cms::Exception("TrajectoryMeasurement::firstMeasurement - Too few measurements in trajectory");
  }

References check(), Exception, and theData.

Referenced by CosmicTrajectoryBuilder::AddHit(), CkfDebugger::analyseCompatibleMeasurements(), CosmicMuonTrajectoryBuilder::build(), CosmicMuonTrajectoryBuilder::buildSecondHalf(), DAFTrackProducerAlgorithm::buildTrack(), MuonTrackLoader::buildTrackExtra(), TrackExtenderWithMTDT< TrackCollection >::buildTrackExtra(), GlobalTrackerMuonAlignment::debugTrajectory(), CosmicMuonTrajectoryBuilder::estimateDirection(), GsfTrajectoryFitter::fitOne(), KFSplittingFitter::fitOne(), KFTrajectoryFitter::fitOne(), NuclearTrackCorrector::getNewTrackExtra(), TransientInitialStateEstimator::innerState(), LowPtGsfElectronSeedProducer::lightGsfTracking(), TrackerToMuonPropagator::produce(), GoodSeedProducer::produce(), GsfTrackProducerBase::putInEvt(), TrackProducerWithSCAssociation::putInEvt(), CosmicMuonTrajectoryBuilder::reverseTrajectory(), CRackTrajectoryBuilder::run(), and CosmicMuonTrajectoryBuilder::trajectories().

foundHits()

int Trajectory::foundHits ( ) const

inline

Number of valid RecHits used to determine the trajectory. Can be less than the number of measurements in data() since detector layers crossed without using RecHits from them are also stored as measurements.

Definition at line 206 of file Trajectory.h.

{ return theNumberOfFoundHits; }

References theNumberOfFoundHits.

Referenced by CosmicTrajectoryBuilder::AddHit(), BaseCkfTrajectoryBuilder::addToResult(), SiStripGainFromData::algoAnalyze(), OverlapValidation::analyze(), DAFTrackProducerAlgorithm::countingGoodHits(), GlobalTrackerMuonAlignment::debugTrajectory(), GsfTrajectoryFitter::fitOne(), CRackTrajectoryBuilder::innerState(), cms::CkfTrackCandidateMakerBase::produceBase(), ThresholdPtTrajectoryFilter::qualityFilter(), MinPtTrajectoryFilter::qualityFilter(), CosmicTrajectoryBuilder::qualityFilter(), CRackTrajectoryBuilder::qualityFilter(), NuclearInteractionFinder::run(), CRackTrajectoryBuilder::run(), DAFTrackProducerAlgorithm::runWithCandidate(), ThresholdPtTrajectoryFilter::toBeContinued(), MinPtTrajectoryFilter::toBeContinued(), and StandAloneMuonTrajectoryBuilder::trajectories().

foundPixelHits()

int Trajectory::foundPixelHits ( ) const

inline

Number of valid pixel RecHits used to determine the trajectory.

Definition at line 210 of file Trajectory.h.

{ return theNumberOfFoundPixelHits; }

References theNumberOfFoundPixelHits.

geometricalInnermostState()

TrajectoryStateOnSurface Trajectory::geometricalInnermostState

(

)

const

Definition at line 217 of file Trajectory.cc.

                                                                     {
  check();
 
  //if trajectory is in one half, return the end closer to origin point
  if (firstMeasurement().updatedState().globalMomentum().perp() > 1.0 &&
      (firstMeasurement().updatedState().globalPosition().basicVector().dot(
           firstMeasurement().updatedState().globalMomentum().basicVector()) *
           lastMeasurement().updatedState().globalPosition().basicVector().dot(
               lastMeasurement().updatedState().globalMomentum().basicVector()) >
       0)) {
    return (firstMeasurement().updatedState().globalPosition().mag() <
            lastMeasurement().updatedState().globalPosition().mag())
               ? firstMeasurement().updatedState()
               : lastMeasurement().updatedState();
  }
 
  //more complicated in case of traversing and low-pt trajectories with loops
  return closestMeasurement(GlobalPoint(0.0, 0.0, 0.0)).updatedState();
}

References RPCNoise_example::check, dot(), mag(), perp(), and TrackInfoProducer_cfi::updatedState.

Referenced by MuonTrackLoader::buildTrackAtPCA().

inactive()

static bool Trajectory::inactive ( )

inlinestatic

Definition of inactive Det from the Trajectory point of view.

Definition at line 267 of file Trajectory.h.

    {
    return false;
  }  //FIXME

incrementLoops()

void Trajectory::incrementLoops ( )

inline

Definition at line 332 of file Trajectory.h.

{ theNLoops++; }

References theNLoops.

invalidate()

void Trajectory::invalidate ( )

inline

Method to invalidate a trajectory. Useful during ambiguity resolution.

Definition at line 260 of file Trajectory.h.

{ theValid = false; }

References theValid.

Referenced by TrajectoryCleanerBySharedHits::clean().

isBad()

bool Trajectory::isBad ( const TrackingRecHit &  hit )

static

Returns true if the hit type is TrackingRecHit::bad Used in stand-alone trajectory construction

Definition at line 163 of file Trajectory.cc.

                                                {
  if (hit.isValid())
    return false;
  else {
    if (hit.geographicalId().rawId() == 0) {
      return false;
    } else {
      return hit.getType() == TrackingRecHit::bad;
    }
  }
}

References TrackingRecHit::bad.

isLooper()

bool Trajectory::isLooper ( ) const

inline

Definition at line 328 of file Trajectory.h.

{ return (theNLoops > 0); }

References theNLoops.

isValid()

bool Trajectory::isValid ( void  ) const

inline

Returns true if the Trajectory is valid. Trajectories are invalidated e.g. during ambiguity resolution.

Definition at line 257 of file Trajectory.h.

{ return theValid; }

References theValid.

Referenced by ntupleDataFormat._Object::_checkIsValid(), MkFitOutputConverter::backwardFit(), DAFTrackProducerAlgorithm::buildTrack(), DAFTrackProducerAlgorithm::calculateNdof(), GlobalTrackerMuonAlignment::debugTrajectory(), DAFTrackProducerAlgorithm::fit(), TransientInitialStateEstimator::innerState(), LowPtGsfElectronSeedProducer::lightGsfTracking(), GoodSeedProducer::produce(), cms::CkfTrackCandidateMakerBase::produceBase(), ChargeSignificanceTrajectoryFilter::qualityFilter(), core.AutoHandle.AutoHandle::ReallyLoad(), NuclearInteractionFinder::run(), DAFTrackProducerAlgorithm::runWithCandidate(), and CosmicMuonTrajectoryBuilder::trajectories().

lastLayer()

const DetLayer* Trajectory::lastLayer ( ) const

inline

Redundant method, returns the layer of lastMeasurement() .

Definition at line 288 of file Trajectory.h.

                                    {
    check();
    if (theData.back().recHit()->hit() != nullptr)
      return theData.back().layer();
    else if (theData.size() > 2)
      return (theData.end() - 2)->layer();
    else
      throw cms::Exception("TrajectoryMeasurement::lastMeasurement - Too few measurements in trajectory");
  }

References check(), Exception, and theData.

Referenced by CkfDebugger::analyseRecHitNotFound().

lastMeasurement()

TrajectoryMeasurement const& Trajectory::lastMeasurement ( ) const

inline

Access to the last measurement. It's the most precise one in a trajectory before smoothing. It's the outermost measurement if direction() == alongMomentum, the innermost one if direction() == oppositeToMomentum.

Definition at line 150 of file Trajectory.h.

                                                       {
    check();
    if (theData.back().recHitR().hit() != nullptr)
      return theData.back();
    else if (theData.size() > 2)
      return *(theData.end() - 2);
    else
      throw cms::Exception("TrajectoryMeasurement::lastMeasurement - Too few measurements in trajectory");
  }

References check(), Exception, and theData.

Referenced by CosmicTrajectoryBuilder::AddHit(), CRackTrajectoryBuilder::AddHit(), BaseCkfTrajectoryBuilder::addToResult(), CkfDebugger::analyseCompatibleMeasurements(), CkfDebugger::analyseRecHitNotFound(), MkFitOutputConverter::backwardFit(), CosmicMuonTrajectoryBuilder::build(), CosmicMuonTrajectoryBuilder::buildSecondHalf(), DAFTrackProducerAlgorithm::buildTrack(), MuonTrackLoader::buildTrackExtra(), TrackExtenderWithMTDT< TrackCollection >::buildTrackExtra(), GlobalTrackerMuonAlignment::debugTrajectory(), CosmicMuonTrajectoryBuilder::estimateDirection(), CosmicMuonTrajectoryBuilder::flipTrajectory(), NuclearTrackCorrector::getNewTrackExtra(), TransientInitialStateEstimator::innerState(), LowPtGsfElectronSeedProducer::lightGsfTracking(), CkfDebugger::nextCorrectHits(), GoodSeedProducer::produce(), GsfTrackProducerBase::putInEvt(), TrackProducerWithSCAssociation::putInEvt(), ThresholdPtTrajectoryFilter::qualityFilter(), MinPtTrajectoryFilter::qualityFilter(), CosmicMuonTrajectoryBuilder::reverseTrajectory(), StandAloneMuonRefitter::singleRefit(), CRackTrajectoryBuilder::SortHits(), ThresholdPtTrajectoryFilter::toBeContinued(), MinPtTrajectoryFilter::toBeContinued(), StandAloneMuonTrajectoryBuilder::trajectories(), and CosmicMuonTrajectoryBuilder::trajectories().

lost()

bool Trajectory::lost ( const TrackingRecHit &  hit )

static

Definition of what it means for a hit to be "lost". This definition is also used by the TrajectoryBuilder.

Definition at line 146 of file Trajectory.cc.

                                               {
  if (hit.isValid())
    return false;
  else {
    //     // A DetLayer is always inactive in this logic.
    //     // The DetLayer is the Det of an invalid RecHit only if no DetUnit
    //     // is compatible with the predicted state, so we don't really expect
    //     // a hit in this case.
 
    if (hit.geographicalId().rawId() == 0) {
      return false;
    } else {
      return hit.getType() == TrackingRecHit::missing;
    }
  }
}

References TrackingRecHit::missing.

Referenced by MaxConsecLostHitsTrajectoryFilter::TBC().

lostHits()

int Trajectory::lostHits ( ) const

inline

Number of detector layers crossed without valid RecHits. Used mainly as a criteria for abandoning a trajectory candidate during trajectory building.

Definition at line 217 of file Trajectory.h.

{ return theNumberOfLostHits; }

References theNumberOfLostHits.

Referenced by BaseCkfTrajectoryBuilder::addToResult(), and OverlapValidation::analyze().

measurements() [1/2]

DataContainer& Trajectory::measurements ( )

inline

Definition at line 179 of file Trajectory.h.

{ return theData; }

References theData.

measurements() [2/2]

DataContainer const& Trajectory::measurements ( ) const

inline

Return all measurements in a container.

Definition at line 178 of file Trajectory.h.

{ return theData; }

References theData.

Referenced by PFTrackTransformer::addPoints(), PFTrackTransformer::addPointsAndBrems(), BaseCkfTrajectoryBuilder::addToResult(), DeDxDiscriminatorLearner::algoAnalyze(), SiStripGainFromData::algoAnalyze(), CkfDebugger::analyseCompatibleMeasurements(), analyzer::SiPixelLorentzAngle::analyze(), OverlapValidation::analyze(), FWTrackProxyBuilderFullFramework::buildTrack(), TrackExtenderWithMTDT< TrackCollection >::buildTrack(), DAFTrackProducerAlgorithm::calculateNdof(), DAFTrackProducerAlgorithm::checkHits(), MuonTrackResidualAnalyzer::computeResolution(), DAFTrackProducerAlgorithm::countingGoodHits(), StandAloneMuonFilter::createDefaultTrajectory(), data(), AlignmentMonitorTemplate::event(), AlignmentMonitorMuonResiduals::event(), AlignmentMonitorTracksFromTrajectories::event(), TrackerValidationVariables::fillHitQuantities(), ApeEstimator::fillTrackVariables(), MuonKinkFinder::fillTrkKink(), DynamicTruncation::filter(), SiStripFineDelayTLA::findtrackangle(), CosmicMuonSmoother::fit(), TransientInitialStateEstimator::innerState(), CRackTrajectoryBuilder::innerState(), MuonResidualsFromTrack::MuonResidualsFromTrack(), Traj2TrackHits::operator()(), TrajectoryFactoryBase::orderedTrajectoryMeasurements(), ShallowTrackClustersProducer::produce(), TkAlCaOverlapTagger::produce(), HLTTrackClusterRemoverNew::produce(), ShallowGainCalibration::produce(), reco::modules::CosmicTrackSplitter::produce(), reco::modules::TrackerTrackHitFilter::produceFromTrajectory(), GsfTrackProducerBase::putInEvt(), StripSubClusterShapeTrajectoryFilter::qualityFilter(), GroupedDAFHitCollector::recHits(), SimpleDAFHitCollector::recHits(), CosmicMuonTrajectoryBuilder::reverseTrajectory(), CosmicMuonTrajectoryBuilder::reverseTrajectoryPropagationDirection(), MuonMillepedeAlgorithm::run(), MuonDTLocalMillepedeAlgorithm::run(), HIPAlignmentAlgorithm::run(), CSCOverlapsAlignmentAlgorithm::run(), NuclearInteractionFinder::run(), GlobalMuonRefitter::selectMuonHits(), ClusterShapeTrajectoryFilter::toBeContinued(), TrajAnnealing::TrajAnnealing(), CosmicMuonTrajectoryBuilder::trajectories(), GsfTrajectorySmoother::trajectory(), KFTrajectorySmoother::trajectory(), trajectoryToResiduals(), and DAFTrackProducerAlgorithm::updateHits().

ndof()

int Trajectory::ndof

(

bool

bon = true

)

const

Number of dof of the trajectory. The method accepts a bool in order to properly take into account the presence of magnetic field in the dof computation. By default the MF is considered ON.

Definition at line 97 of file Trajectory.cc.

                                   {
  Trajectory::RecHitContainer&& transRecHits = recHits();
 
  int dof = 0;
  int dofBad = 0;
 
  for (auto& rechit : transRecHits) {
    if ((rechit)->isValid())
      dof += (rechit)->dimension();
    else if (isBad(*rechit) && (rechit)->geographicalId().det() == DetId::Muon)
      dofBad += (rechit)->dimension();
  }
 
  // If dof!=0 (there is at least 1 valid hit),
  //    return ndof=ndof(fit)
  // If dof=0 (all rec hits are invalid, only for STA trajectories),
  //    return ndof=ndof(invalid hits)
  if (dof) {
    int constr = bon ? 5 : 4;
    return std::max(dof - constr, 0);
  } else {
    // A STA can have < 5 (invalid) hits
    // if this is the case ==> ndof = 1
    // (to avoid divisions by 0)
    int constr = bon ? 5 : 4;
    return std::max(dofBad - constr, 1);
  }
}

References pat::helper::ParametrizationHelper::dimension(), SiStripPI::max, DetId::Muon, and FastTrackerRecHitMaskProducer_cfi::recHits.

Referenced by analyzer::SiPixelLorentzAngle::analyze(), OverlapValidation::analyze(), TrackExtenderWithMTDT< TrackCollection >::buildTrack(), and MuonTrackLoader::buildTrackAtPCA().

nLoops()

signed char Trajectory::nLoops ( ) const

inline

Definition at line 329 of file Trajectory.h.

{ return theNLoops; }

References theNLoops.

Referenced by TransientInitialStateEstimator::innerState().

numberOfCCCBadHits()

int Trajectory::numberOfCCCBadHits

(

float

ccc_threshold

)

Definition at line 212 of file Trajectory.cc.

                                                      {
  updateBadForCCC(ccc_threshold);
  return theNumberOfCCCBadHits_;
}

Referenced by TempTrajectory::toTrajectory().

operator=() [1/2]

Trajectory& Trajectory::operator= ( Trajectory &&  rh )

inline

Definition at line 101 of file Trajectory.h.

                                         {
    using std::swap;
    swap(theData, rh.theData);
    theChiSquared = rh.theChiSquared;
    theChiSquaredBad = rh.theChiSquaredBad;
    theValid = rh.theValid;
    theDPhiCache = rh.theDPhiCache;
    theCCCThreshold_ = rh.theCCCThreshold_;
    theNLoops = rh.theNLoops;
    theNumberOfFoundHits = rh.theNumberOfFoundHits;
    theNumberOfFoundPixelHits = rh.theNumberOfFoundPixelHits;
    theNumberOfLostHits = rh.theNumberOfLostHits;
    theNumberOfTrailingFoundHits = rh.theNumberOfTrailingFoundHits;
    theNumberOfCCCBadHits_ = rh.theNumberOfCCCBadHits_;
    theDirection = rh.theDirection;
    theDirectionValidity = rh.theDirectionValidity;
    stopReason_ = rh.stopReason_;
    swap(theSeed, rh.theSeed);
    swap(seedRef_, rh.seedRef_);
 
    return *this;
  }

References seedRef_, stopReason_, edm::swap(), std::swap(), theCCCThreshold_, theChiSquared, theChiSquaredBad, theData, theDirection, theDirectionValidity, theDPhiCache, theNLoops, theNumberOfCCCBadHits_, theNumberOfFoundHits, theNumberOfFoundPixelHits, theNumberOfLostHits, theNumberOfTrailingFoundHits, theSeed, and theValid.

operator=() [2/2]

Trajectory& Trajectory::operator= ( Trajectory const &  rh )

default

pixel()

bool Trajectory::pixel ( const TrackingRecHit &  hit )

static

Returns true if the hit type is TrackingRecHit::bad Used in trajectory filtering

Definition at line 175 of file Trajectory.cc.

                                                {
  if (!trackerHitRTTI::isFromDetOrFast(hit))
    return false;
  auto const* thit = static_cast<const BaseTrackerRecHit*>(hit.hit());
  return thit->isPixel();
}

References trackerHitRTTI::isFromDetOrFast().

Referenced by SeedExtensionTrajectoryFilter::looseTBC(), TempTrajectory::pop(), and TempTrajectory::pushAux().

pop()

void Trajectory::pop

(

)

Remove the last measurement from the trajectory.

Definition at line 30 of file Trajectory.cc.

                     {
  if (!empty()) {
    if (theData.back().recHit()->isValid()) {
      theNumberOfFoundHits--;
      theChiSquared -= theData.back().estimate();
      if (badForCCC(theData.back()))
        theNumberOfCCCBadHits_--;
      if (pixel(*(theData.back().recHit())))
        theNumberOfFoundPixelHits--;
    } else if (lost(*(theData.back().recHit()))) {
      theNumberOfLostHits--;
    } else if (isBad(*(theData.back().recHit())) && theData.back().recHit()->geographicalId().det() == DetId::Muon) {
      theChiSquaredBad -= theData.back().estimate();
    }
 
    theData.pop_back();
    theNumberOfTrailingFoundHits = countTrailingValidHits(theData);
  }
}

References std::empty(), DetId::Muon, and muonClassificationByHits_cfi::pixel.

Referenced by BaseCkfTrajectoryBuilder::addToResult(), cms::CkfTrackCandidateMakerBase::produceBase(), CosmicMuonTrajectoryBuilder::reverseTrajectoryPropagationDirection(), and CosmicMuonTrajectoryBuilder::trajectories().

push() [1/4]

void Trajectory::push

(

const TrajectoryMeasurement &

tm

)

Add a new measurement to a Trajectory. The Chi2 of the trajectory is incremented by the value of tm.estimate() .

Definition at line 50 of file Trajectory.cc.

{ push(tm, tm.estimate()); }

References TrajectoryMeasurement::estimate().

Referenced by CosmicTrajectoryBuilder::AddHit(), CRackTrajectoryBuilder::AddHit(), StandAloneMuonFilter::createDefaultTrajectory(), CosmicMuonSmoother::fit(), GsfTrajectoryFitter::fitOne(), NuclearTrackCorrector::newTrajNeeded(), MuonMillepedeTrackRefitter::produce(), TrackerToMuonPropagator::produce(), CSCOverlapsTrackPreparation::produce(), cms::CkfTrackCandidateMakerBase::produceBase(), TempTrajectory::toTrajectory(), TrackTransformer::transform(), and MuonTrajectoryUpdator::update().

push() [2/4]

void Trajectory::push	(	const TrajectoryMeasurement &	tm,
		double	chi2Increment
	)

same as the one-argument push, but the trajectory Chi2 is incremented by chi2Increment. Useful e.g. in trajectory smoothing.

Definition at line 54 of file Trajectory.cc.

                                                                           {
  theData.push_back(tm);
  pushAux(chi2Increment);
}

push() [3/4]

void Trajectory::push

(

TrajectoryMeasurement &&

tm

)

Definition at line 52 of file Trajectory.cc.

{ push(tm, tm.estimate()); }

push() [4/4]

void Trajectory::push	(	TrajectoryMeasurement &&	tm,
		double	chi2Increment
	)

Definition at line 59 of file Trajectory.cc.

                                                                      {
  theData.push_back(tm);
  pushAux(chi2Increment);
}

pushAux()

void Trajectory::pushAux ( double  chi2Increment )

private

Definition at line 64 of file Trajectory.cc.

                                             {
  const TrajectoryMeasurement& tm = theData.back();
  if (tm.recHit()->isValid()) {
    theChiSquared += chi2Increment;
    theNumberOfFoundHits++;
    theNumberOfTrailingFoundHits++;
    if (badForCCC(tm))
      theNumberOfCCCBadHits_++;
    if (pixel(*(tm.recHit())))
      theNumberOfFoundPixelHits++;
  }
  // else if (lost( tm.recHit()) && !inactive(tm.recHit().det())) theNumberOfLostHits++;
  else if (lost(*(tm.recHit()))) {
    theNumberOfLostHits++;
    theNumberOfTrailingFoundHits = 0;
  }
 
  else if (isBad(*(tm.recHit())) && tm.recHit()->geographicalId().det() == DetId::Muon) {
    theChiSquaredBad += chi2Increment;
  }
 
  // in case of a Trajectory constructed without direction,
  // determine direction from the radii of the first two measurements
 
  if (!theDirectionValidity && theData.size() >= 2) {
    if (theData[0].updatedState().globalPosition().perp2() < theData.back().updatedState().globalPosition().perp2())
      theDirection = alongMomentum;
    else
      theDirection = oppositeToMomentum;
    theDirectionValidity = true;
  }
}

References alongMomentum, DetId::Muon, oppositeToMomentum, muonClassificationByHits_cfi::pixel, TrajectoryMeasurement::recHit(), and TrackInfoProducer_cfi::updatedState.

recHits()

ConstRecHitContainer Trajectory::recHits ( void  ) const

inline

Return all RecHits in a container.

Definition at line 186 of file Trajectory.h.

                                       {
    ConstRecHitContainer hits;
    hits.reserve(theData.size());
    for (Trajectory::DataContainer::const_iterator itm = theData.begin(); itm != theData.end(); itm++) {
      hits.push_back((*itm).recHit());
    }
    return hits;
  }

References hfClusterShapes_cfi::hits, and theData.

Referenced by SiStripGainFromData::algoAnalyze(), MuonTrackLoader::buildTrackExtra(), TrackExtenderWithMTDT< TrackCollection >::buildTrackExtra(), CkfDebugger::correctTrajectory(), CosmicMuonTrajectoryBuilder::estimateDirection(), GsfTrajectoryFitter::fitOne(), KFSplittingFitter::fitOne(), KFTrajectoryFitter::fitOne(), CosmicMuonTrajectoryBuilder::flipTrajectory(), CosmicMuonTrajectoryBuilder::getDirectionByTime(), CosmicTrajectoryBuilder::qualityFilter(), CRackTrajectoryBuilder::qualityFilter(), StandAloneMuonRefitter::refit(), CRackTrajectoryBuilder::run(), CosmicMuonTrajectoryBuilder::selfDuplicate(), and StandAloneMuonRefitter::singleRefit().

reserve()

void Trajectory::reserve ( unsigned int  n )

inline

Reserves space in the vector to avoid lots of allocations when push_back-ing measurements

Definition at line 126 of file Trajectory.h.

{ theData.reserve(n); }

References dqmiodumpmetadata::n, and theData.

Referenced by StandAloneMuonFilter::createDefaultTrajectory(), cms::CkfTrackCandidateMakerBase::produceBase(), and TempTrajectory::toTrajectory().

reverse()

void Trajectory::reverse

(

)

Reverse the propagation direction and the order of the trajectory measurements. It doesn't reverse the forward and backward predicted states within each trajectory measurement

Definition at line 264 of file Trajectory.cc.

                         {
  // reverse the direction (without changing it if it's not along or opposite)
  if (theDirection == alongMomentum)
    theDirection = oppositeToMomentum;
  else if (theDirection == oppositeToMomentum)
    theDirection = alongMomentum;
  // reverse the order of the hits
  std::reverse(theData.begin(), theData.end());
}

References alongMomentum, oppositeToMomentum, and groupFilesInBlocks::reverse.

seed()

TrajectorySeed const& Trajectory::seed ( ) const

inline

Access to the seed used to reconstruct the Trajectory.

Definition at line 263 of file Trajectory.h.

{ return *theSeed; }

References theSeed.

Referenced by CosmicTrajectoryBuilder::AddHit(), CosmicMuonTrajectoryBuilder::estimateDirection(), DAFTrackProducerAlgorithm::fit(), GsfTrajectoryFitter::fitOne(), KFSplittingFitter::fitOne(), KFTrajectoryFitter::fitOne(), CosmicMuonTrajectoryBuilder::flipTrajectory(), NuclearTrackCorrector::getTrackFromTrajectory(), CosmicMuonTrajectoryBuilder::reverseTrajectory(), CosmicMuonTrajectoryBuilder::reverseTrajectoryPropagationDirection(), CRackTrajectoryBuilder::run(), DAFTrackProducerAlgorithm::runWithCandidate(), seedNHits(), setSharedSeed(), StandAloneMuonRefitter::singleRefit(), GsfTrajectorySmoother::trajectory(), and KFTrajectorySmoother::trajectory().

seedNHits()

unsigned int Trajectory::seedNHits ( ) const

inline

Definition at line 230 of file Trajectory.h.

{ return seed().nHits(); }

References TrajectorySeed::nHits(), and seed().

seedRef()

edm::RefToBase<TrajectorySeed> Trajectory::seedRef ( void  ) const

inline

return the Reference to the trajectory seed in the original seeds collection. If the collection has been dropped from the Event, the reference may be invalid. Its validity should be tested, before the reference is actually used.

Definition at line 303 of file Trajectory.h.

{ return seedRef_; }

References seedRef_.

Referenced by MuonTrackLoader::buildTrackExtra(), TrackExtenderWithMTDT< TrackCollection >::buildTrackExtra(), GsfTrackProducerBase::putInEvt(), TrackProducerWithSCAssociation::putInEvt(), and setSeedRef().

setDPhiCacheForLoopersReconstruction()

void Trajectory::setDPhiCacheForLoopersReconstruction ( float  dphi )

inline

method to set the delta phi angle betweem the directions of the two measurements on the last two layers crossed by the trajectory

Definition at line 326 of file Trajectory.h.

{ theDPhiCache = dphi; }

References theDPhiCache.

setNLoops()

void Trajectory::setNLoops ( signed char  value )

inline

Definition at line 331 of file Trajectory.h.

{ theNLoops = value; }

References theNLoops, and relativeConstraints::value.

Referenced by cms::CkfTrackCandidateMakerBase::produceBase(), and TempTrajectory::toTrajectory().

setSeedRef()

void Trajectory::setSeedRef ( const edm::RefToBase< TrajectorySeed > &  seedRef )

inline

Definition at line 305 of file Trajectory.h.

{ seedRef_ = seedRef; }

References seedRef(), and seedRef_.

Referenced by cms::CkfTrackCandidateMakerBase::produceBase().

setSharedSeed()

void Trajectory::setSharedSeed ( const std::shared_ptr< const TrajectorySeed > &  seed )

inline

Definition at line 316 of file Trajectory.h.

{ theSeed = seed; }

References seed(), and theSeed.

Referenced by BaseCkfTrajectoryBuilder::addToResult().

setStopReason()

void Trajectory::setStopReason ( StopReason  s )

inline

Definition at line 334 of file Trajectory.h.

{ stopReason_ = s; }

References alignCSCRings::s, and stopReason_.

Referenced by cms::CkfTrackCandidateMakerBase::produceBase(), MinPtTrajectoryFilter::toBeContinued(), and TempTrajectory::toTrajectory().

sharedSeed()

const std::shared_ptr<const TrajectorySeed>& Trajectory::sharedSeed ( ) const

inline

Definition at line 315 of file Trajectory.h.

{ return theSeed; }

References theSeed.

stopReason()

StopReason Trajectory::stopReason ( ) const

inline

Definition at line 335 of file Trajectory.h.

{ return stopReason_; }

References stopReason_.

trailingFoundHits()

int Trajectory::trailingFoundHits ( ) const

inline

Number of valid RecHits at the end of the trajectory after last lost hit.

Definition at line 221 of file Trajectory.h.

{ return theNumberOfTrailingFoundHits; }

References theNumberOfTrailingFoundHits.

updateBadForCCC()

void Trajectory::updateBadForCCC ( float  ccc_threshold )

private

Definition at line 197 of file Trajectory.cc.

                                                    {
  // If the supplied threshold is the same as the currently cached
  // one, then return the current number of bad hits for CCC,
  // otherwise do a new full rescan.
  if (ccc_threshold == theCCCThreshold_)
    return;
 
  theCCCThreshold_ = ccc_threshold;
  theNumberOfCCCBadHits_ = 0;
  for (auto const& h : theData) {
    if (badForCCC(h))
      theNumberOfCCCBadHits_++;
  }
}

validRecHits()

void Trajectory::validRecHits

(

ConstRecHitContainer &

cont

)

const

Just valid hits..

Definition at line 126 of file Trajectory.cc.

                                                              {
  hits.reserve(foundHits());
  for (auto const& tm : theData)
    if (tm.recHit()->isValid())
      hits.push_back(tm.recHit());
}

References hfClusterShapes_cfi::hits.

Referenced by NuclearTrackCorrector::getTrackFromTrajectory().

Member Data Documentation

seedRef_

edm::RefToBase<TrajectorySeed> Trajectory::seedRef_

private

Definition at line 345 of file Trajectory.h.

Referenced by operator=(), seedRef(), and setSeedRef().

stopReason_

StopReason Trajectory::stopReason_ = StopReason::UNINITIALIZED

private

Definition at line 364 of file Trajectory.h.

Referenced by operator=(), setStopReason(), and stopReason().

theCCCThreshold_

float Trajectory::theCCCThreshold_ = std::numeric_limits<float>::max()

private

Definition at line 362 of file Trajectory.h.

Referenced by cccThreshold(), and operator=().

theChiSquared

float Trajectory::theChiSquared = 0

private

Definition at line 348 of file Trajectory.h.

Referenced by chiSquared(), and operator=().

theChiSquaredBad

float Trajectory::theChiSquaredBad = 0

private

Definition at line 349 of file Trajectory.h.

Referenced by chiSquared(), and operator=().

theData

DataContainer Trajectory::theData

private

Definition at line 347 of file Trajectory.h.

Referenced by empty(), firstMeasurement(), lastLayer(), lastMeasurement(), measurements(), operator=(), recHits(), and reserve().

theDirection

PropagationDirection Trajectory::theDirection = anyDirection

private

Definition at line 357 of file Trajectory.h.

Referenced by operator=().

theDirectionValidity

bool Trajectory::theDirectionValidity = false

private

Definition at line 358 of file Trajectory.h.

Referenced by operator=().

theDPhiCache

float Trajectory::theDPhiCache = 0

private

Definition at line 361 of file Trajectory.h.

Referenced by dPhiCacheForLoopersReconstruction(), operator=(), and setDPhiCacheForLoopersReconstruction().

theNLoops

signed char Trajectory::theNLoops = 0

private

Definition at line 363 of file Trajectory.h.

Referenced by incrementLoops(), isLooper(), nLoops(), operator=(), and setNLoops().

theNumberOfCCCBadHits_

signed short Trajectory::theNumberOfCCCBadHits_ = 0

private

Definition at line 355 of file Trajectory.h.

Referenced by cccBadHits(), and operator=().

theNumberOfFoundHits

signed short Trajectory::theNumberOfFoundHits = 0

private

Definition at line 351 of file Trajectory.h.

Referenced by chiSquared(), foundHits(), and operator=().

theNumberOfFoundPixelHits

signed short Trajectory::theNumberOfFoundPixelHits = 0

private

Definition at line 352 of file Trajectory.h.

Referenced by foundPixelHits(), and operator=().

theNumberOfLostHits

signed short Trajectory::theNumberOfLostHits = 0

private

Definition at line 353 of file Trajectory.h.

Referenced by lostHits(), and operator=().

theNumberOfTrailingFoundHits

signed short Trajectory::theNumberOfTrailingFoundHits = 0

private

Definition at line 354 of file Trajectory.h.

Referenced by operator=(), and trailingFoundHits().

theSeed

std::shared_ptr<const TrajectorySeed> Trajectory::theSeed

private

Definition at line 344 of file Trajectory.h.

Referenced by operator=(), seed(), setSharedSeed(), and sharedSeed().

theValid

bool Trajectory::theValid = false

private

Definition at line 359 of file Trajectory.h.

Referenced by invalidate(), isValid(), and operator=().