PandaRoot
genfit::Track Class Reference

Collection of TrackPoint objects, AbsTrackRep objects and FitStatus objects. More...

#include <Track.h>

Inheritance diagram for genfit::Track:

Public Member Functions

 Track ()
 
 Track (const TrackCand &trackCand, const MeasurementFactory< genfit::AbsMeasurement > &factory, AbsTrackRep *rep=nullptr)
 
 Track (AbsTrackRep *trackRep, const TVectorD &stateSeed)
 
 Track (AbsTrackRep *trackRep, const TVector3 &posSeed, const TVector3 &momSeed)
 
 Track (AbsTrackRep *trackRep, const TVectorD &stateSeed, const TMatrixDSym &covSeed)
 
 Track (const Track &)
 
Trackoperator= (Track)
 
void swap (Track &other)
 
virtual ~Track ()
 
virtual void Clear (Option_t *="")
 
TrackPointgetPoint (int id) const
 
const std::vector< genfit::TrackPoint * > & getPoints () const
 
unsigned int getNumPoints () const
 
TrackPointgetPointWithMeasurement (int id) const
 
const std::vector< genfit::TrackPoint * > & getPointsWithMeasurement () const
 
unsigned int getNumPointsWithMeasurement () const
 
TrackPointgetPointWithMeasurementAndFitterInfo (int id, const AbsTrackRep *rep) const
 
TrackPointgetPointWithFitterInfo (int id, const AbsTrackRep *rep=nullptr) const
 
const MeasuredStateOnPlanegetFittedState (int id=0, const AbsTrackRep *rep=nullptr, bool biased=true) const
 Shortcut to get FittedStates. More...
 
AbsTrackRepgetTrackRep (int id) const
 
unsigned int getNumReps () const
 
int getIdForRep (const AbsTrackRep *rep) const
 This is used when streaming TrackPoints. More...
 
AbsTrackRepgetCardinalRep () const
 Get cardinal track representation. More...
 
unsigned int getCardinalRepId () const
 
int getMcTrackId () const
 Get the MCT track id, for MC simulations - default value = -1. More...
 
bool hasFitStatus (const AbsTrackRep *rep=nullptr) const
 Check if track has a FitStatus for given AbsTrackRep. Per default, check for cardinal rep. More...
 
FitStatusgetFitStatus (const AbsTrackRep *rep=nullptr) const
 Get FitStatus for a AbsTrackRep. Per default, return FitStatus for cardinalRep. More...
 
bool hasKalmanFitStatus (const AbsTrackRep *rep=nullptr) const
 Check if track has a KalmanFitStatus for given AbsTrackRep. Per default, check for cardinal rep. More...
 
KalmanFitStatusgetKalmanFitStatus (const AbsTrackRep *rep=nullptr) const
 If FitStatus is a KalmanFitStatus, return it. Otherwise return nullptr. More...
 
void setFitStatus (FitStatus *fitStatus, const AbsTrackRep *rep)
 
double getTimeSeed () const
 
void setTimeSeed (double time)
 
const TVectorD & getStateSeed () const
 
void setStateSeed (const TVectorD &s)
 
void setStateSeed (const TVector3 &pos, const TVector3 &mom)
 
const TMatrixDSym & getCovSeed () const
 
void setCovSeed (const TMatrixDSym &c)
 
void setMcTrackId (int i)
 Set the MCT track id, for MC simulations. More...
 
void insertPoint (TrackPoint *point, int id=-1)
 Insert TrackPoint BEFORE TrackPoint with position id, if id >= 0. More...
 
void insertPoints (std::vector< genfit::TrackPoint *> points, int id=-1)
 Insert TrackPoints BEFORE TrackPoint with position id, if id >= 0. More...
 
void deletePoint (int id)
 
void insertMeasurement (AbsMeasurement *measurement, int id=-1)
 Creates a new TrackPoint containing the measurement, and adds it to the track. More...
 
void mergeTrack (const Track *other, int id=-1)
 Merge two tracks. More...
 
void addTrackRep (AbsTrackRep *trackRep)
 
void deleteTrackRep (int id)
 Delete a AbsTrackRep and all corresponding AbsFitterInfo objects in every TrackPoint. More...
 
void setCardinalRep (int id)
 
void determineCardinalRep ()
 See with which AbsTrackRep the track was fitted best (converged fit w/ smallest chi2) and set the cardinal rep accordingly. More...
 
bool sort ()
 Sort TrackPoint and according to their sorting parameters. More...
 
bool udpateSeed (int id=0, AbsTrackRep *rep=nullptr, bool biased=true)
 
void reverseTrackPoints ()
 Flip the ordering of the TrackPoints. More...
 
void reverseMomSeed ()
 Flip direction of momentum seed. More...
 
void switchPDGSigns (AbsTrackRep *rep=nullptr)
 Switch the pdg signs of specified rep (of all reps if rep == nullptr). More...
 
void reverseTrack ()
 Make track ready to be fitted in reverse direction. More...
 
void deleteForwardInfo (int startId=0, int endId=-1, const AbsTrackRep *rep=nullptr)
 
void deleteBackwardInfo (int startId=0, int endId=-1, const AbsTrackRep *rep=nullptr)
 
void deleteReferenceInfo (int startId=0, int endId=-1, const AbsTrackRep *rep=nullptr)
 
void deleteMeasurementInfo (int startId=0, int endId=-1, const AbsTrackRep *rep=nullptr)
 
void deleteFitterInfo (int startId=0, int endId=-1, const AbsTrackRep *rep=nullptr)
 
double getTrackLen (AbsTrackRep *rep=nullptr, int startId=0, int endId=-1) const
 get TrackLength between to trackPoints (if nullptr, for cardinal rep) More...
 
double getTOF (AbsTrackRep *rep=nullptr, int startId=0, int endId=-1) const
 get time of flight in ns between to trackPoints (if nullptr, for cardinal rep) More...
 
TrackCandconstructTrackCand () const
 Construct a new TrackCand containing the hit IDs of the measurements. More...
 
void fixWeights (AbsTrackRep *rep=nullptr, int startId=0, int endId=-1)
 
void prune (const Option_t *="CFLWRMIU")
 Delete unneeded information from the Track. More...
 
void Print (const Option_t *="") const
 
bool checkConsistency () const
 

Detailed Description

Collection of TrackPoint objects, AbsTrackRep objects and FitStatus objects.

Holds a number of AbsTrackRep objects, which correspond to the different particle hypotheses or track models which should be fitted. A 6D seed #stateSeed_ (x,y,z,p_x,p_y,p_z) and 6x6 #covSeed_ should be provided as start values for fitting. When fitting the Track with a AbsFitter, a FitStatus object will be created, containing information about the fit. The fitted states will be stored in AbsFitterInfo objects in every TrackPoints.

The fit will be performed for every AbsTrackRep, so after the fit there will be one AbsFitterInfo for each AbsTrackRep in every TrackPoint, as well as one FitStatus for every AbsTrackRep.

Definition at line 67 of file Track.h.

Constructor & Destructor Documentation

◆ Track() [1/6]

genfit::Track::Track ( )

◆ Track() [2/6]

genfit::Track::Track ( const TrackCand trackCand,
const MeasurementFactory< genfit::AbsMeasurement > &  factory,
AbsTrackRep rep = nullptr 
)

@ brief Construct Track from TrackCand, using a MeasurementFactory

The MeasurementFactory will be used to create AbsMeasuremen objects. TrackPoints will be created. If two or more consecutive PlanarMeasurement objects with the same detector- and planeId are created by the factory, they will be put into the same TrackPoint.

Optionally, a AbsTrackRep can be provided.

The stateSeed_ of the Track will be filled with the seed of the TrackCand. A guess for covSeed_ will be made using the largest entry of the cov of the first measurement and the number of measurements (For the covSeed_, it is just important that it will be big enough not to bias the fit too much, but not too big in order to avoid numerical problems).

◆ Track() [3/6]

genfit::Track::Track ( AbsTrackRep trackRep,
const TVectorD &  stateSeed 
)

◆ Track() [4/6]

genfit::Track::Track ( AbsTrackRep trackRep,
const TVector3 &  posSeed,
const TVector3 &  momSeed 
)

◆ Track() [5/6]

genfit::Track::Track ( AbsTrackRep trackRep,
const TVectorD &  stateSeed,
const TMatrixDSym &  covSeed 
)

◆ Track() [6/6]

genfit::Track::Track ( const Track )

◆ ~Track()

virtual genfit::Track::~Track ( )
virtual

Member Function Documentation

◆ addTrackRep()

void genfit::Track::addTrackRep ( AbsTrackRep trackRep)

◆ checkConsistency()

bool genfit::Track::checkConsistency ( ) const

◆ Clear()

virtual void genfit::Track::Clear ( Option_t *  = "")
virtual

◆ constructTrackCand()

TrackCand* genfit::Track::constructTrackCand ( ) const

Construct a new TrackCand containing the hit IDs of the measurements.

The idea is hat you can get a TrackCand for storing the hit IDs after a track has been fitted. His could have been reordered, added or removed, so that the original TackCand no longer represents the Track correctly. You might want to call determineCardinalRep() and/or udpateSeed() before.

◆ deleteBackwardInfo()

void genfit::Track::deleteBackwardInfo ( int  startId = 0,
int  endId = -1,
const AbsTrackRep rep = nullptr 
)

◆ deleteFitterInfo()

void genfit::Track::deleteFitterInfo ( int  startId = 0,
int  endId = -1,
const AbsTrackRep rep = nullptr 
)

◆ deleteForwardInfo()

void genfit::Track::deleteForwardInfo ( int  startId = 0,
int  endId = -1,
const AbsTrackRep rep = nullptr 
)

◆ deleteMeasurementInfo()

void genfit::Track::deleteMeasurementInfo ( int  startId = 0,
int  endId = -1,
const AbsTrackRep rep = nullptr 
)

◆ deletePoint()

void genfit::Track::deletePoint ( int  id)

◆ deleteReferenceInfo()

void genfit::Track::deleteReferenceInfo ( int  startId = 0,
int  endId = -1,
const AbsTrackRep rep = nullptr 
)

◆ deleteTrackRep()

void genfit::Track::deleteTrackRep ( int  id)

Delete a AbsTrackRep and all corresponding AbsFitterInfo objects in every TrackPoint.

◆ determineCardinalRep()

void genfit::Track::determineCardinalRep ( )

See with which AbsTrackRep the track was fitted best (converged fit w/ smallest chi2) and set the cardinal rep accordingly.

◆ fixWeights()

void genfit::Track::fixWeights ( AbsTrackRep rep = nullptr,
int  startId = 0,
int  endId = -1 
)

Helper function: For all KalmanFitterInfos belonging to rep (if nullptr, for all reps), call the fixWeights() function, so that e.g. the DAF will not alter weights anymore.

◆ getCardinalRep()

AbsTrackRep* genfit::Track::getCardinalRep ( ) const
inline

Get cardinal track representation.

The user has to choose which AbsTrackRep should be considered the best one after the fit. E.g. the track representation giving the smallest chi2 could be chosen. By default the first in the list is returned.

See also
determineCardinalRep()

Definition at line 136 of file Track.h.

136 { return trackReps_.at(cardinalRep_); }

◆ getCardinalRepId()

unsigned int genfit::Track::getCardinalRepId ( ) const
inline

Definition at line 137 of file Track.h.

137 { return cardinalRep_; }

◆ getCovSeed()

const TMatrixDSym& genfit::Track::getCovSeed ( ) const
inline

Definition at line 170 of file Track.h.

170 { return covSeed_; }

◆ getFitStatus()

FitStatus* genfit::Track::getFitStatus ( const AbsTrackRep rep = nullptr) const
inline

Get FitStatus for a AbsTrackRep. Per default, return FitStatus for cardinalRep.

Definition at line 145 of file Track.h.

146  {
147  if (rep == nullptr)
148  rep = getCardinalRep();
149  return fitStatuses_.at(rep);
150  }
AbsTrackRep * getCardinalRep() const
Get cardinal track representation.
Definition: Track.h:136

◆ getFittedState()

const MeasuredStateOnPlane& genfit::Track::getFittedState ( int  id = 0,
const AbsTrackRep rep = nullptr,
bool  biased = true 
) const

Shortcut to get FittedStates.

Uses getPointWithFitterInfo(id, rep). Gets the fitted state at trackpoint id for the track representation rep. Per default, the fitted state of the fitterInfo of the first TrackPoint with one or more AbsFitterInfo objects is returned. If no AbsTrackRep is specified, the AbsFitterInfo of the cardinal rep will be used.

◆ getIdForRep()

int genfit::Track::getIdForRep ( const AbsTrackRep rep) const

This is used when streaming TrackPoints.

◆ getKalmanFitStatus()

KalmanFitStatus* genfit::Track::getKalmanFitStatus ( const AbsTrackRep rep = nullptr) const

If FitStatus is a KalmanFitStatus, return it. Otherwise return nullptr.

◆ getMcTrackId()

int genfit::Track::getMcTrackId ( ) const
inline

Get the MCT track id, for MC simulations - default value = -1.

Definition at line 140 of file Track.h.

140 { return mcTrackId_; }

◆ getNumPoints()

unsigned int genfit::Track::getNumPoints ( ) const
inline

Definition at line 103 of file Track.h.

103 { return trackPoints_.size(); }

◆ getNumPointsWithMeasurement()

unsigned int genfit::Track::getNumPointsWithMeasurement ( ) const
inline

Definition at line 107 of file Track.h.

107 { return trackPointsWithMeasurement_.size(); }

◆ getNumReps()

unsigned int genfit::Track::getNumReps ( ) const
inline

Definition at line 124 of file Track.h.

124 { return trackReps_.size(); }

◆ getPoint()

TrackPoint* genfit::Track::getPoint ( int  id) const

◆ getPoints()

const std::vector<genfit::TrackPoint *>& genfit::Track::getPoints ( ) const
inline

Definition at line 102 of file Track.h.

102 { return trackPoints_; }

◆ getPointsWithMeasurement()

const std::vector<genfit::TrackPoint *>& genfit::Track::getPointsWithMeasurement ( ) const
inline

Definition at line 106 of file Track.h.

106 { return trackPointsWithMeasurement_; }

◆ getPointWithFitterInfo()

TrackPoint* genfit::Track::getPointWithFitterInfo ( int  id,
const AbsTrackRep rep = nullptr 
) const

◆ getPointWithMeasurement()

TrackPoint* genfit::Track::getPointWithMeasurement ( int  id) const

◆ getPointWithMeasurementAndFitterInfo()

TrackPoint* genfit::Track::getPointWithMeasurementAndFitterInfo ( int  id,
const AbsTrackRep rep 
) const

◆ getStateSeed()

const TVectorD& genfit::Track::getStateSeed ( ) const
inline

Definition at line 162 of file Track.h.

162 { return stateSeed_; }

◆ getTimeSeed()

double genfit::Track::getTimeSeed ( ) const
inline

Definition at line 159 of file Track.h.

159 { return timeSeed_; }

◆ getTOF()

double genfit::Track::getTOF ( AbsTrackRep rep = nullptr,
int  startId = 0,
int  endId = -1 
) const

get time of flight in ns between to trackPoints (if nullptr, for cardinal rep)

◆ getTrackLen()

double genfit::Track::getTrackLen ( AbsTrackRep rep = nullptr,
int  startId = 0,
int  endId = -1 
) const

get TrackLength between to trackPoints (if nullptr, for cardinal rep)

◆ getTrackRep()

AbsTrackRep* genfit::Track::getTrackRep ( int  id) const
inline

Definition at line 123 of file Track.h.

123 { return trackReps_.at(id); }

◆ hasFitStatus()

bool genfit::Track::hasFitStatus ( const AbsTrackRep rep = nullptr) const

Check if track has a FitStatus for given AbsTrackRep. Per default, check for cardinal rep.

◆ hasKalmanFitStatus()

bool genfit::Track::hasKalmanFitStatus ( const AbsTrackRep rep = nullptr) const

Check if track has a KalmanFitStatus for given AbsTrackRep. Per default, check for cardinal rep.

◆ insertMeasurement()

void genfit::Track::insertMeasurement ( AbsMeasurement measurement,
int  id = -1 
)

Creates a new TrackPoint containing the measurement, and adds it to the track.

◆ insertPoint()

void genfit::Track::insertPoint ( TrackPoint point,
int  id = -1 
)

Insert TrackPoint BEFORE TrackPoint with position id, if id >= 0.

Id -1 means after last TrackPoint. Id -2 means before last TrackPoint. ... Also deletes backwardInfos before new point and forwardInfos after new point. Also sets Track backpointer of point accordingly.

◆ insertPoints()

void genfit::Track::insertPoints ( std::vector< genfit::TrackPoint *>  points,
int  id = -1 
)

Insert TrackPoints BEFORE TrackPoint with position id, if id >= 0.

Id -1 means after last TrackPoint. Id -2 means before last TrackPoint. ... Also deletes backwardInfos before and for new points and forwardInfos after and for new points. Also sets Track backpointers of points accordingly.

◆ mergeTrack()

void genfit::Track::mergeTrack ( const Track other,
int  id = -1 
)

Merge two tracks.

The TrackPoint objects of other will be cloned and inserted after id (per default, they will be appended at the end). The other Track will not be altered, the TrackPoint objects will be (deep) copied. Only copies the TrackPoint objects, NOT the AbsTrackRep, FitStatus, seed state and other objects of the other track.

◆ operator=()

Track& genfit::Track::operator= ( Track  )

◆ Print()

void genfit::Track::Print ( const Option_t *  = "") const

◆ prune()

void genfit::Track::prune ( const Option_t *  = "CFLWRMIU")

Delete unneeded information from the Track.

Possible options: (see also PruneFlags defined in FitStatus.h) C: prune all reps except cardinalRep F: prune all points except first point (also prune referenceInfo from fitterInfos) L: prune all points except last point (also prune referenceInfo from fitterInfos) FL: prune all points except first and last point (also prune referenceInfo from fitterInfos) W: prune rawMeasurements from TrackPoints R: prune referenceInfo from fitterInfos M: prune measurementInfo from fitterInfos I: if F, L, or FL is set, prune forward (backward) info of first (last) point U: if fitterInfo is a KalmanFitterInfo, prune predictions and keep updates

◆ reverseMomSeed()

void genfit::Track::reverseMomSeed ( )
inline

Flip direction of momentum seed.

Definition at line 237 of file Track.h.

238  {
239  stateSeed_(3) *= -1;
240  stateSeed_(4) *= -1;
241  stateSeed_(5) *= -1;
242  }

◆ reverseTrack()

void genfit::Track::reverseTrack ( )

Make track ready to be fitted in reverse direction.

Flip the order of TrackPoints and the momentum direction of the seed state. If possible, take the smoothed state of the last hit as new seed state. Flip charge of the TrackReps.

◆ reverseTrackPoints()

void genfit::Track::reverseTrackPoints ( )

Flip the ordering of the TrackPoints.

◆ setCardinalRep()

void genfit::Track::setCardinalRep ( int  id)

◆ setCovSeed()

void genfit::Track::setCovSeed ( const TMatrixDSym &  c)
inline

Definition at line 171 of file Track.h.

172  {
173  covSeed_.ResizeTo(c);
174  covSeed_ = c;
175  }

◆ setFitStatus()

void genfit::Track::setFitStatus ( FitStatus fitStatus,
const AbsTrackRep rep 
)

◆ setMcTrackId()

void genfit::Track::setMcTrackId ( int  i)
inline

Set the MCT track id, for MC simulations.

Definition at line 178 of file Track.h.

References i.

178 { mcTrackId_ = i; }
unsigned int i
Definition: P4_F32vec4.h:21

◆ setStateSeed() [1/2]

void genfit::Track::setStateSeed ( const TVectorD &  s)
inline

Definition at line 163 of file Track.h.

164  {
165  stateSeed_.ResizeTo(s);
166  stateSeed_ = s;
167  }

◆ setStateSeed() [2/2]

void genfit::Track::setStateSeed ( const TVector3 &  pos,
const TVector3 &  mom 
)

◆ setTimeSeed()

void genfit::Track::setTimeSeed ( double  time)
inline

Definition at line 160 of file Track.h.

160 { timeSeed_ = time; }

◆ sort()

bool genfit::Track::sort ( )

Sort TrackPoint and according to their sorting parameters.

Returns if the order of the TrackPoint has actually changed.

◆ swap()

void genfit::Track::swap ( Track other)

◆ switchPDGSigns()

void genfit::Track::switchPDGSigns ( AbsTrackRep rep = nullptr)

Switch the pdg signs of specified rep (of all reps if rep == nullptr).

◆ udpateSeed()

bool genfit::Track::udpateSeed ( int  id = 0,
AbsTrackRep rep = nullptr,
bool  biased = true 
)

Try to set the fitted state as seed. Return if it was successfull. Adapt the sign of all TrackReps' pdg to the actual fitted charge.


The documentation for this class was generated from the following file: