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

#include <Track.h>

Inheritance diagram for genfit::Track:

Collaboration diagram for genfit::Track:

Public Member Functions
	Track ()

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

	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 &)

Track &	operator= (Track)

void	swap (Track &other)

virtual	~Track ()

virtual void	Clear (Option_t *="")

TrackPoint *	getPoint (int id) const

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

unsigned int	getNumPoints () const

TrackPoint *	getPointWithMeasurement (int id) const

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

unsigned int	getNumPointsWithMeasurement () const

TrackPoint *	getPointWithMeasurementAndFitterInfo (int id, const AbsTrackRep *rep) const

const MeasuredStateOnPlane &	getFittedState (int id=0, const AbsTrackRep *rep=NULL, bool biased=true) const
	Shortcut to get FittedStates.

AbsTrackRep *	getTrackRep (int id) const

unsigned int	getNumReps () const

int	getIdForRep (const AbsTrackRep *rep) const
	This is used when streaming TrackPoints.

AbsTrackRep *	getCardinalRep () const
	Get cardinal track representation.

unsigned int	getCardinalRepId () const

int	getMcTrackId () const
	Get the MCT track id, for MC simulations - default value = -1.

bool	hasFitStatus (const AbsTrackRep *rep=NULL) const
	Check if track has a FitStatus for given AbsTrackRep. Per default, check for cardinal rep.

FitStatus *	getFitStatus (const AbsTrackRep *rep=NULL) const
	Get FitStatus for a AbsTrackRep. Per default, return FitStatus for cardinalRep.

bool	hasKalmanFitStatus (const AbsTrackRep *rep=NULL) const
	Check if track has a KalmanFitStatus for given AbsTrackRep. Per default, check for cardinal rep.

KalmanFitStatus *	getKalmanFitStatus (const AbsTrackRep *rep=NULL) const
	If FitStatus is a KalmanFitStatus, return it. Otherwise return NULL.

void	setFitStatus (FitStatus fitStatus, const AbsTrackRep rep)

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.

void	insertPoint (TrackPoint *point, int id=-1)
	Insert TrackPoint BEFORE TrackPoint with position id, if id >= 0.

void	insertPoints (std::vector< genfit::TrackPoint * > points, int id=-1)
	Insert TrackPoints BEFORE TrackPoint with position id, if id >= 0.

void	deletePoint (int id)

void	insertMeasurement (AbsMeasurement *measurement, int id=-1)
	Creates a new TrackPoint containing the measurement, and adds it to the track.

void	mergeTrack (const Track *other, int id=-1)
	Merge two tracks.

void	addTrackRep (AbsTrackRep *trackRep)

void	deleteTrackRep (int id)
	Delete a AbsTrackRep and all corresponding AbsFitterInfo objects in every TrackPoint.

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.

bool	sort ()
	Sort TrackPoint and according to their sorting parameters.

bool	udpateSeed (int id=0, AbsTrackRep *rep=NULL, bool biased=true)

void	reverseTrackPoints ()
	Flip the ordering of the TrackPoints.

void	reverseMomSeed ()
	Flip direction of momentum seed.

void	switchPDGSigns (AbsTrackRep *rep=NULL)
	Switch the pdg signs of specified rep (of all reps if rep == NULL).

void	reverseTrack ()
	Make track ready to be fitted in reverse direction.

void	deleteForwardInfo (int startId=0, int endId=-1, const AbsTrackRep *rep=NULL)

void	deleteBackwardInfo (int startId=0, int endId=-1, const AbsTrackRep *rep=NULL)

void	deleteReferenceInfo (int startId=0, int endId=-1, const AbsTrackRep *rep=NULL)

void	deleteMeasurementInfo (int startId=0, int endId=-1, const AbsTrackRep *rep=NULL)

void	deleteFitterInfo (int startId=0, int endId=-1, const AbsTrackRep *rep=NULL)

double	getTrackLen (AbsTrackRep *rep=NULL, int startId=0, int endId=-1) const
	get TrackLength between to trackPoints (if NULL, for cardinal rep)

double	getTOF (AbsTrackRep *rep=NULL, int startId=0, int endId=-1) const
	get time of flight in ns between to trackPoints (if NULL, for cardinal rep)

TrackCand *	constructTrackCand () const
	Construct a new TrackCand containing the hit IDs of the measurements.

void	fixWeights (AbsTrackRep *rep=NULL, int startId=0, int endId=-1)

void	prune (const Option_t *="CFLWRMIU")
	Delete unneeded information from the Track.

void	Print (const Option_t *="") const

bool	checkConsistency () const

Private Member Functions
void	trackHasChanged ()

void	fillPointsWithMeasurement ()

Private Attributes
std::vector< AbsTrackRep * >	trackReps_

unsigned int	cardinalRep_

std::vector< TrackPoint * >	trackPoints_

std::vector< TrackPoint * >	trackPointsWithMeasurement_

std::map< const AbsTrackRep , FitStatus >	fitStatuses_
	helper

int	mcTrackId_

TVectorD	stateSeed_

TMatrixDSym	covSeed_

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 71 of file Track.h.

Constructor & Destructor Documentation

◆ Track() [1/6]

genfit::Track::Track ( )

Definition at line 42 of file Track.cc.

             :
  cardinalRep_(0), fitStatuses_(), stateSeed_(6), covSeed_(6)
{
  ;
}

◆ Track() [2/6]

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

@ 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).

Definition at line 49 of file Track.cc.

                                                                                                                  :
  cardinalRep_(0), fitStatuses_(), stateSeed_(6), covSeed_(6)
{
 
  if (rep != NULL)
    addTrackRep(rep);
 
  // create the measurements using the factory.
  std::vector <genfit::AbsMeasurement*> factoryHits = factory.createMany(trackCand);
 
  // create TrackPoints
  PlanarMeasurement* lastPlanarMeas(NULL);
  PlanarMeasurement* planarMeas(NULL);
  for (unsigned int i=0; i<factoryHits.size(); ++i){
    planarMeas = dynamic_cast<PlanarMeasurement*>(factoryHits[i]);
 
    if (lastPlanarMeas != NULL && planarMeas != NULL &&
        lastPlanarMeas->getDetId() == planarMeas->getDetId() &&
        planarMeas->getPlaneId() != -1 &&   // -1 is default plane id
        lastPlanarMeas->getPlaneId() == planarMeas->getPlaneId() ) {
      trackPoints_.back()->addRawMeasurement(factoryHits[i]);
    }
    else
      insertPoint(new genfit::TrackPoint(factoryHits[i], this));
 
    lastPlanarMeas = dynamic_cast<PlanarMeasurement*>(factoryHits[i]);
  }
 
  // fill seed state
  stateSeed_ = trackCand.getStateSeed();
 
  // initial guess for cov
  covSeed_ = trackCand.getCovSeed();
 
  // fill cache
  fillPointsWithMeasurement();
 
  // self test
  assert(checkConsistency());
}

◆ Track() [3/6]

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

Definition at line 91 of file Track.cc.

                                                             :
  cardinalRep_(0), fitStatuses_(), stateSeed_(stateSeed),
  covSeed_(TMatrixDSym::kUnit, TMatrixDSym(6))
{
  addTrackRep(trackRep);
}

◆ Track() [4/6]

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

Definition at line 99 of file Track.cc.

                                                                                    :
  cardinalRep_(0), fitStatuses_(), stateSeed_(6),
  covSeed_(TMatrixDSym::kUnit, TMatrixDSym(6))
{
  addTrackRep(trackRep);
  setStateSeed(posSeed, momSeed);
}

◆ Track() [5/6]

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

Definition at line 108 of file Track.cc.

                                                                                         :
  cardinalRep_(0), fitStatuses_(), stateSeed_(stateSeed), covSeed_(covSeed)
{
  addTrackRep(trackRep);
}

◆ Track() [6/6]

genfit::Track::Track ( const Track & rhs )

Definition at line 115 of file Track.cc.

                             :
  TObject(rhs),
  cardinalRep_(rhs.cardinalRep_), stateSeed_(rhs.stateSeed_), covSeed_(rhs.covSeed_)
{
  assert(rhs.checkConsistency());
 
  std::map<const AbsTrackRep*, AbsTrackRep*> oldRepNewRep;
 
  for (std::vector<AbsTrackRep*>::const_iterator it=rhs.trackReps_.begin(); it!=rhs.trackReps_.end(); ++it) {
    AbsTrackRep* newRep = (*it)->clone();
    addTrackRep(newRep);
    oldRepNewRep[(*it)] = newRep;
  }
 
  trackPoints_.reserve(rhs.trackPoints_.size());
  for (std::vector<TrackPoint*>::const_iterator it=rhs.trackPoints_.begin(); it!=rhs.trackPoints_.end(); ++it) {
    trackPoints_.push_back(new TrackPoint(**it, oldRepNewRep));
    trackPoints_.back()->setTrack(this);
  }
 
  for (std::map< const AbsTrackRep*, FitStatus* >::const_iterator it=rhs.fitStatuses_.begin(); it!=rhs.fitStatuses_.end(); ++it) {
    setFitStatus(it->second->clone(), oldRepNewRep[it->first]);
  }
 
  fillPointsWithMeasurement();
 
  // self test
  assert(checkConsistency());
}

◆ ~Track()

genfit::Track::~Track ( )

virtual

Definition at line 171 of file Track.cc.

              {
 // causes problem for python, is it needed ? TR 2014  this->Clear();
 this->Clear();
}

Member Function Documentation

◆ addTrackRep()

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

Definition at line 506 of file Track.cc.

                                             {
  trackReps_.push_back(trackRep);
  fitStatuses_[trackRep] = new FitStatus();
}

◆ checkConsistency()

bool genfit::Track::checkConsistency ( ) const

Definition at line 1192 of file Track.cc.

                                   {
 
  bool retVal(true);
 
  std::map<const AbsTrackRep*, const KalmanFitterInfo*> prevFis;
 
  // check if seed is 6D
  if (stateSeed_.GetNrows() != 6) {
    std::cerr << "Track::checkConsistency(): stateSeed_ dimension != 6" << std::endl;
    retVal = false;
  }
 
  if (covSeed_.GetNrows() != 6) {
    std::cerr << "Track::checkConsistency(): covSeed_ dimension != 6" << std::endl;
    retVal = false;
  }
 
  if (covSeed_.Max() == 0.) {
    std::cerr << "Track::checkConsistency(): Warning: covSeed_ is zero" << std::endl;
    //retVal = false;
  }
 
  // check if cardinalRep_ is in range of trackReps_
  if (trackReps_.size() && cardinalRep_ >= trackReps_.size()) {
    std::cerr << "Track::checkConsistency(): cardinalRep id " << cardinalRep_ << " out of bounds" << std::endl;
    retVal = false;
  }
 
  for (std::vector<AbsTrackRep*>::const_iterator rep = trackReps_.begin(); rep != trackReps_.end(); ++rep) {
    // check for NULL
    if ((*rep) == NULL) {
      std::cerr << "Track::checkConsistency(): TrackRep is NULL" << std::endl;
      retVal = false;
    }
 
    // check for valid pdg code
    TParticlePDG* particle = TDatabasePDG::Instance()->GetParticle((*rep)->getPDG());
    if (particle == NULL) {
      std::cerr << "Track::checkConsistency(): TrackRep pdg ID " << (*rep)->getPDG() << " is not valid" << std::endl;
      retVal = false;
    }
 
  }
 
  // check TrackPoints
  for (std::vector<TrackPoint*>::const_iterator tp = trackPoints_.begin(); tp != trackPoints_.end(); ++tp) {
    // check for NULL
    if ((*tp) == NULL) {
      std::cerr << "Track::checkConsistency(): TrackPoint is NULL" << std::endl;
      retVal = false;
    }
    // check if trackPoint points back to this track
    if ((*tp)->getTrack() != this) {
      std::cerr << "Track::checkConsistency(): TrackPoint does not point back to this track" << std::endl;
      retVal = false;
    }
 
    // check rawMeasurements
    const std::vector<AbsMeasurement*>& rawMeasurements = (*tp)->getRawMeasurements();
    for (std::vector<AbsMeasurement*>::const_iterator m = rawMeasurements.begin(); m != rawMeasurements.end(); ++m) {
      // check for NULL
      if ((*m) == NULL) {
        std::cerr << "Track::checkConsistency(): Measurement is NULL" << std::endl;
        retVal = false;
      }
      // check if measurement points back to TrackPoint
      if ((*m)->getTrackPoint() != *tp) {
        std::cerr << "Track::checkConsistency(): Measurement does not point back to correct TrackPoint" << std::endl;
        retVal = false;
      }
    }
 
    // check fitterInfos
    std::vector<AbsFitterInfo*> fitterInfos = (*tp)->getFitterInfos();
    for (std::vector<AbsFitterInfo*>::const_iterator fi = fitterInfos.begin(); fi != fitterInfos.end(); ++fi) {
      // check for NULL
      if ((*fi) == NULL) {
        std::cerr << "Track::checkConsistency(): FitterInfo is NULL. TrackPoint: " << *tp << std::endl;
        retVal = false;
      }
 
      if (!( (*fi)->checkConsistency() ) ) {
        std::cerr << "Track::checkConsistency(): FitterInfo not consistent. TrackPoint: " << *tp << std::endl;
        retVal = false;
      }
 
      // check if fitterInfos point to valid TrackReps in trackReps_
      int mycount (0);
      for (std::vector<AbsTrackRep*>::const_iterator rep = trackReps_.begin(); rep != trackReps_.end(); ++rep) {
        if ( (*rep) == (*fi)->getRep() ) {
          ++mycount;
        }
      }
      if (mycount ==  0) {
        std::cerr << "Track::checkConsistency(): fitterInfo points to TrackRep which is not in Track" << std::endl;
        retVal = false;
      }
 
      if (dynamic_cast<KalmanFitterInfo*>(*fi) != NULL) {
        if (prevFis[(*fi)->getRep()] != NULL &&
            static_cast<KalmanFitterInfo*>(*fi)->hasReferenceState() &&
            prevFis[(*fi)->getRep()]->hasReferenceState() ) {
          double len = static_cast<KalmanFitterInfo*>(*fi)->getReferenceState()->getForwardSegmentLength();
          double prevLen = prevFis[(*fi)->getRep()]->getReferenceState()->getBackwardSegmentLength();
          if (fabs(prevLen + len) > 1E-10 ) {
            std::cerr << "Track::checkConsistency(): segment lengths of reference states for rep " << (*fi)->getRep() << " (id " << getIdForRep((*fi)->getRep()) << ") at TrackPoint " << (*tp) << " don't match" << std::endl;
            std::cerr << prevLen << " + " << len << " = " << prevLen + len << std::endl;
            std::cerr << "TrackPoint " << *tp << ", FitterInfo " << *fi << ", rep " << getIdForRep((*fi)->getRep()) << std::endl;
            retVal = false;
          }
        }
 
        prevFis[(*fi)->getRep()] = static_cast<KalmanFitterInfo*>(*fi);
      }
      else
        prevFis[(*fi)->getRep()] = NULL;
 
    } // end loop over FitterInfos
 
  } // end loop over TrackPoints
 
 
  // check trackPointsWithMeasurement_
  std::vector<TrackPoint*> trackPointsWithMeasurement;
 
  for (std::vector<TrackPoint*>::const_iterator it = trackPoints_.begin(); it != trackPoints_.end(); ++it) {
    if ((*it)->hasRawMeasurements()) {
      trackPointsWithMeasurement.push_back(*it);
    }
  }
 
  if (trackPointsWithMeasurement.size() != trackPointsWithMeasurement_.size()) {
    std::cerr << "Track::checkConsistency(): trackPointsWithMeasurement_ has incorrect size" << std::endl;
    retVal = false;
  }
 
  for (unsigned int i = 0; i < trackPointsWithMeasurement.size(); ++i) {
    if (trackPointsWithMeasurement[i] != trackPointsWithMeasurement_[i]) {
      std::cerr << "Track::checkConsistency(): trackPointsWithMeasurement_ is not correct" << std::endl;
      std::cerr << "has         id " << i << ", adress " << trackPointsWithMeasurement_[i] << std::endl;
      std::cerr << "should have id " << i << ", adress " << trackPointsWithMeasurement[i] << std::endl;
      retVal = false;
    }
  }
 
  return retVal;
}

◆ Clear()

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

virtual

Definition at line 176 of file Track.cc.

{
  // This function is needed for TClonesArray embedding.
  // FIXME: smarter containers or pointers needed ...
  for (size_t i = 0; i < trackPoints_.size(); ++i)
    delete trackPoints_[i];
 
  trackPoints_.clear();
  trackPointsWithMeasurement_.clear();
 
  for (std::map< const AbsTrackRep*, FitStatus* >::iterator it = fitStatuses_.begin(); it!= fitStatuses_.end(); ++it)
    delete it->second;
  fitStatuses_.clear();
 
  for (size_t i = 0; i < trackReps_.size(); ++i)
    delete trackReps_[i];
  trackReps_.clear();
 
  cardinalRep_ = 0;
 
  stateSeed_ *= 0;
  covSeed_ *= 0;
}

◆ constructTrackCand()

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

inline

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.

Definition at line 260 of file Track.h.

260{ return 0; /* implementation missing FIXME RDM */ }

◆ deleteBackwardInfo()

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

Definition at line 710 of file Track.cc.

                                                                             {
 
  #ifdef DEBUG
  std::cout << "Track::deleteBackwardInfo from position " << startId  << " to " << endId << "\n";
  #endif
 
  trackHasChanged();
 
  if (startId < 0)
    startId += trackPoints_.size();
  if (endId < 0)
    endId += trackPoints_.size();
  endId += 1;
 
  assert (endId >= startId);
 
 
  for (std::vector<TrackPoint*>::const_iterator pointIt = trackPoints_.begin() + startId; pointIt != trackPoints_.begin() + endId; ++pointIt) {
    if (rep != NULL) {
      if ((*pointIt)->hasFitterInfo(rep))
        (*pointIt)->getFitterInfo(rep)->deleteBackwardInfo();
    }
    else {
      const std::vector<AbsFitterInfo*> fitterInfos = (*pointIt)->getFitterInfos();
      for (std::vector<AbsFitterInfo*>::const_iterator fitterInfoIt = fitterInfos.begin(); fitterInfoIt != fitterInfos.end(); ++fitterInfoIt) {
        (*fitterInfoIt)->deleteBackwardInfo();
      }
    }
  }
}

◆ deleteFitterInfo()

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

Definition at line 801 of file Track.cc.

                                                                           {
 
  #ifdef DEBUG
  std::cout << "Track::deleteFitterInfo from position " << startId  << " to " << endId << "\n";
  #endif
 
  trackHasChanged();
 
  if (startId < 0)
    startId += trackPoints_.size();
  if (endId < 0)
    endId += trackPoints_.size();
  endId += 1;
 
  assert (endId >= startId);
 
  for (std::vector<TrackPoint*>::const_iterator pointIt = trackPoints_.begin() + startId; pointIt != trackPoints_.begin() + endId; ++pointIt) {
    if (rep != NULL) {
      if ((*pointIt)->hasFitterInfo(rep))
        (*pointIt)->deleteFitterInfo(rep);
    }
    else {
      for (std::vector<AbsTrackRep*>::const_iterator repIt = trackReps_.begin(); repIt != trackReps_.end(); ++repIt) {
        if ((*pointIt)->hasFitterInfo(*repIt))
          (*pointIt)->deleteFitterInfo(*repIt);
      }
    }
  }
}

◆ deleteForwardInfo()

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

Definition at line 681 of file Track.cc.

                                                                            {
  #ifdef DEBUG
  std::cout << "Track::deleteForwardInfo from position " << startId  << " to " << endId << "\n";
  #endif
 
  trackHasChanged();
 
  if (startId < 0)
    startId += trackPoints_.size();
  if (endId < 0)
    endId += trackPoints_.size();
  endId += 1;
 
  assert (endId >= startId);
 
  for (std::vector<TrackPoint*>::const_iterator pointIt = trackPoints_.begin() + startId; pointIt != trackPoints_.begin() + endId; ++pointIt) {
    if (rep != NULL) {
      if ((*pointIt)->hasFitterInfo(rep))
        (*pointIt)->getFitterInfo(rep)->deleteForwardInfo();
    }
    else {
      const std::vector<AbsFitterInfo*> fitterInfos = (*pointIt)->getFitterInfos();
      for (std::vector<AbsFitterInfo*>::const_iterator fitterInfoIt = fitterInfos.begin(); fitterInfoIt != fitterInfos.end(); ++fitterInfoIt) {
        (*fitterInfoIt)->deleteForwardInfo();
      }
    }
  }
}

◆ deleteMeasurementInfo()

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

Definition at line 771 of file Track.cc.

                                                                                {
 
  #ifdef DEBUG
  std::cout << "Track::deleteMeasurementInfo from position " << startId  << " to " << endId << "\n";
  #endif
 
  trackHasChanged();
 
  if (startId < 0)
    startId += trackPoints_.size();
  if (endId < 0)
    endId += trackPoints_.size();
  endId += 1;
 
  assert (endId >= startId);
 
  for (std::vector<TrackPoint*>::const_iterator pointIt = trackPoints_.begin() + startId; pointIt != trackPoints_.begin() + endId; ++pointIt) {
    if (rep != NULL) {
      if ((*pointIt)->hasFitterInfo(rep))
        (*pointIt)->getFitterInfo(rep)->deleteMeasurementInfo();
    }
    else {
      std::vector<AbsFitterInfo*> fitterInfos = (*pointIt)->getFitterInfos();
      for (std::vector<AbsFitterInfo*>::const_iterator fitterInfoIt = fitterInfos.begin(); fitterInfoIt != fitterInfos.end(); ++fitterInfoIt) {
        (*fitterInfoIt)->deleteMeasurementInfo();
      }
    }
  }
}

◆ deletePoint()

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

Definition at line 417 of file Track.cc.

                              {
 
  #ifdef DEBUG
  std::cout << "Track::deletePoint at position " << id  << "\n";
  #endif
 
  trackHasChanged();
 
  if (id < 0)
    id += trackPoints_.size();
  assert(id>0);
 
 
  // delete forwardInfo after point (backwardInfo before point) if deleted point has a measurement
  if (trackPoints_[id]->hasRawMeasurements()) {
    deleteForwardInfo(id, -1);
    deleteBackwardInfo(0, id-1);
  }
 
  // delete reference info of neighbouring points
  deleteReferenceInfo(std::max(0, id-1), std::min((int)trackPoints_.size()-1, id+1));
 
  // delete point
  std::vector<TrackPoint*>::iterator it = std::find(trackPointsWithMeasurement_.begin(), trackPointsWithMeasurement_.end(), trackPoints_[id]);
  if (it != trackPointsWithMeasurement_.end())
    trackPointsWithMeasurement_.erase(it);
 
  delete trackPoints_[id];
  trackPoints_.erase(trackPoints_.begin()+id);
 
  fillPointsWithMeasurement();
 
}

◆ deleteReferenceInfo()

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

Definition at line 741 of file Track.cc.

                                                                              {
 
  #ifdef DEBUG
  std::cout << "Track::deleteReferenceInfo from position " << startId  << " to " << endId << "\n";
  #endif
 
  trackHasChanged();
 
  if (startId < 0)
    startId += trackPoints_.size();
  if (endId < 0)
    endId += trackPoints_.size();
  endId += 1;
 
  assert (endId >= startId);
 
  for (std::vector<TrackPoint*>::const_iterator pointIt = trackPoints_.begin() + startId; pointIt != trackPoints_.begin() + endId; ++pointIt) {
    if (rep != NULL) {
      if ((*pointIt)->hasFitterInfo(rep))
        (*pointIt)->getFitterInfo(rep)->deleteReferenceInfo();
    }
    else {
      std::vector<AbsFitterInfo*> fitterInfos = (*pointIt)->getFitterInfos();
      for (std::vector<AbsFitterInfo*>::const_iterator fitterInfoIt = fitterInfos.begin(); fitterInfoIt != fitterInfos.end(); ++fitterInfoIt) {
        (*fitterInfoIt)->deleteReferenceInfo();
      }
    }
  }
}

◆ deleteTrackRep()

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

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

Definition at line 512 of file Track.cc.

                                 {
  if (id < 0)
    id += trackReps_.size();
 
  AbsTrackRep* rep = trackReps_.at(id);
 
  // update cardinalRep_
  if (int(cardinalRep_) == id)
    cardinalRep_ = 0; // reset
  else if (int(cardinalRep_) > id)
    --cardinalRep_; // make cardinalRep_ point to the same TrackRep before and after deletion
 
  // delete FitterInfos related to the deleted TrackRep
  for (std::vector<TrackPoint*>::const_iterator pointIt = trackPoints_.begin(); pointIt != trackPoints_.end(); ++pointIt) {
    (*pointIt)->deleteFitterInfo(rep);
  }
 
  // delete fitStatus
  delete fitStatuses_.at(rep);
  fitStatuses_.erase(rep);
 
  // delete rep
  delete rep;
  trackReps_.erase(trackReps_.begin()+id);
}

◆ 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.

Definition at line 553 of file Track.cc.

                                 {
 
  // Todo: test
 
  if (trackReps_.size() <= 1)
    return;
 
  double minChi2(9.E99);
  const AbsTrackRep* bestRep(NULL);
 
  for (std::map< const AbsTrackRep*, FitStatus* >::const_iterator it = fitStatuses_.begin(); it != fitStatuses_.end(); ++it) {
    if (it->second->isFitConverged()) {
      if (it->second->getChi2() < minChi2) {
        minChi2 = it->second->getChi2();
        bestRep = it->first;
      }
    }
  }
 
  if (bestRep != NULL) {
    setCardinalRep(getIdForRep(bestRep));
  }
}

◆ fillPointsWithMeasurement()

void genfit::Track::fillPointsWithMeasurement ( )

private

Definition at line 1356 of file Track.cc.

                                      {
  trackPointsWithMeasurement_.clear();
  trackPointsWithMeasurement_.reserve(trackPoints_.size());
 
  for (std::vector<TrackPoint*>::const_iterator it = trackPoints_.begin(); it != trackPoints_.end(); ++it) {
    if ((*it)->hasRawMeasurements()) {
      trackPointsWithMeasurement_.push_back(*it);
    }
  }
}

◆ fixWeights()

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

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

Definition at line 919 of file Track.cc.

                                                               {
 
  if (startId < 0)
    startId += trackPoints_.size();
  if (endId < 0)
    endId += trackPoints_.size();
 
  assert(startId >= 0);
  assert(startId <= endId);
  assert(endId <= (int)trackPoints_.size());
 
  std::vector< genfit::AbsFitterInfo* > fis;
 
  for (std::vector<TrackPoint*>::iterator tp = trackPoints_.begin() + startId; tp != trackPoints_.begin() + endId; ++tp) {
    fis.clear();
    if (rep == NULL) {
      fis = (*tp)->getFitterInfos();
    }
    else if ((*tp)->hasFitterInfo(rep)) {
      fis.push_back((*tp)->getFitterInfo(rep));
    }
 
    for (std::vector< genfit::AbsFitterInfo* >::iterator fi = fis.begin(); fi != fis.end(); ++fi) {
      KalmanFitterInfo* kfi = dynamic_cast<KalmanFitterInfo*>(*fi);
      if (kfi == NULL)
        continue;
 
      kfi->fixWeights();
    }
  }
}

◆ 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 140 of file Track.h.

140{return trackReps_.at(cardinalRep_);}

◆ getCardinalRepId()

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

inline

Definition at line 141 of file Track.h.

141{return cardinalRep_;}

◆ getCovSeed()

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

inline

Definition at line 162 of file Track.h.

162{return covSeed_;}

◆ getFitStatus()

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

inline

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

Definition at line 149 of file Track.h.

149{if (rep == NULL) rep = getCardinalRep(); return fitStatuses_.at(rep);}

genfit::Track::getCardinalRep

AbsTrackRep * getCardinalRep() const

Get cardinal track representation.

Definition Track.h:140

◆ getFittedState()

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

Shortcut to get FittedStates.

Uses getPointWithMeasurementAndFitterInfo(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 AbsMeasurement and AbsFitterInfo objects is returned. If no AbsTrackRep is specified, the AbsFitterInfo of the cardinal rep will be used.

Definition at line 231 of file Track.cc.

                                                                                                   {
  if (rep == NULL)
    rep = getCardinalRep();
 
  TrackPoint* point = getPointWithMeasurementAndFitterInfo(id, rep);
  if (point == NULL) {
    Exception exc("Track::getFittedState ==> no trackPoint with fitterInfo for rep",__LINE__,__FILE__);
    exc.setFatal();
    throw exc;
  }
  return point->getFitterInfo(rep)->getFittedState(biased);
}

◆ getIdForRep()

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

This is used when streaming TrackPoints.

Definition at line 245 of file Track.cc.

{
  for (size_t i = 0; i < trackReps_.size(); ++i)
    if (trackReps_[i] == rep)
      return i;
 
  // std::cout << "track debug " << trackReps_.size() <<" "<< trackReps_[0] << " " << rep <<"\n";
 
  //assert(0 == 1);  // Cannot happen.
  return 0;  //RDM
}

◆ getKalmanFitStatus()

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

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

Definition at line 280 of file Track.cc.

                                                                       {
  return dynamic_cast<KalmanFitStatus*>(getFitStatus(rep));
}

◆ getMcTrackId()

int genfit::Track::getMcTrackId ( ) const

inline

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

Definition at line 144 of file Track.h.

144{return mcTrackId_;}

genfit::Track::mcTrackId_

int mcTrackId_

Definition Track.h:300

◆ getNumPoints()

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

inline

Definition at line 108 of file Track.h.

108{return trackPoints_.size();}

◆ getNumPointsWithMeasurement()

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

inline

Definition at line 112 of file Track.h.

112{return trackPointsWithMeasurement_.size();}

◆ getNumReps()

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

inline

Definition at line 128 of file Track.h.

128{return trackReps_.size();}

◆ getPoint()

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

Definition at line 201 of file Track.cc.

                                        {
  if (id < 0)
    id += trackPoints_.size();
 
  return trackPoints_.at(id);
}

◆ getPoints()

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

inline

Definition at line 107 of file Track.h.

107{return trackPoints_;}

◆ getPointsWithMeasurement()

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

inline

Definition at line 111 of file Track.h.

111{return trackPointsWithMeasurement_;}

◆ getPointWithMeasurement()

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

Definition at line 209 of file Track.cc.

                                                       {
  if (id < 0)
    id += trackPointsWithMeasurement_.size();
 
  return trackPointsWithMeasurement_.at(id);
}

◆ getPointWithMeasurementAndFitterInfo()

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

Definition at line 217 of file Track.cc.

                                                                                            {
  int i(0);
  for (std::vector<TrackPoint*>::const_iterator it = trackPointsWithMeasurement_.begin(); it != trackPointsWithMeasurement_.end(); ++it) {
    if ((*it)->hasFitterInfo(rep)) {
      if (id == i)
        return (*it);
      ++i;
    }
  }
 
  return NULL;
}

◆ getStateSeed()

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

inline

Definition at line 158 of file Track.h.

158{return stateSeed_;}

◆ getTOF()

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

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

Definition at line 875 of file Track.cc.

                                                                   {
 
  if (startId < 0)
    startId += trackPoints_.size();
  if (endId < 0)
    endId += trackPoints_.size();
 
  bool backwards(false);
  if (startId > endId) {
    double temp = startId;
    startId = endId;
    endId = temp;
    backwards = true;
  }
 
  endId += 1;
 
  if (rep == NULL)
    rep = getCardinalRep();
 
  double tof(0);
  StateOnPlane state;
 
  for (std::vector<TrackPoint*>::const_iterator pointIt = trackPoints_.begin() + startId; pointIt != trackPoints_.begin() + endId; ++pointIt) {
    if (! (*pointIt)->hasFitterInfo(rep)) {
      Exception e("Track::getTOF: trackPoint has no fitterInfo", __LINE__,__FILE__);
      throw e;
    }
 
    if (pointIt != trackPoints_.begin() + startId) {
      rep->extrapolateToPlane(state, (*pointIt)->getFitterInfo(rep)->getPlane());
      tof += rep->getTOF();
    }
 
    state = (*pointIt)->getFitterInfo(rep)->getFittedState();
  }
 
  if (backwards)
    tof *= -1.;
 
  return tof;
}

◆ getTrackLen()

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

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

Definition at line 832 of file Track.cc.

                                                                        {
 
  if (startId < 0)
    startId += trackPoints_.size();
  if (endId < 0)
    endId += trackPoints_.size();
 
  bool backwards(false);
  if (startId > endId) {
    double temp = startId;
    startId = endId;
    endId = temp;
    backwards = true;
  }
 
  endId += 1;
 
  if (rep == NULL)
    rep = getCardinalRep();
 
  double trackLen(0);
  StateOnPlane state;
 
  for (std::vector<TrackPoint*>::const_iterator pointIt = trackPoints_.begin() + startId; pointIt != trackPoints_.begin() + endId; ++pointIt) {
    if (! (*pointIt)->hasFitterInfo(rep)) {
      Exception e("Track::getTracklength: trackPoint has no fitterInfo", __LINE__,__FILE__);
      throw e;
    }
 
    if (pointIt != trackPoints_.begin() + startId) {
      trackLen += rep->extrapolateToPlane(state, (*pointIt)->getFitterInfo(rep)->getPlane());
    }
 
    state = (*pointIt)->getFitterInfo(rep)->getFittedState();
  }
 
  if (backwards)
    trackLen *= -1.;
 
  return trackLen;
}

◆ getTrackRep()

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

inline

Definition at line 127 of file Track.h.

127{return trackReps_.at(id);}

◆ hasFitStatus()

bool genfit::Track::hasFitStatus ( const AbsTrackRep * rep = NULL ) const

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

Definition at line 258 of file Track.cc.

                                                     {
  if (rep == NULL)
    rep = getCardinalRep();
 
  if (fitStatuses_.find(rep) == fitStatuses_.end())
    return false;
 
  return (fitStatuses_.at(rep) != NULL);
}

◆ hasKalmanFitStatus()

bool genfit::Track::hasKalmanFitStatus ( const AbsTrackRep * rep = NULL ) const

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

Definition at line 269 of file Track.cc.

                                                           {
  if (rep == NULL)
    rep = getCardinalRep();
 
  if (fitStatuses_.find(rep) == fitStatuses_.end())
    return false;
 
  return (dynamic_cast<KalmanFitStatus*>(fitStatuses_.at(rep)) != NULL);
}

◆ insertMeasurement()

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

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

Definition at line 452 of file Track.cc.

                                                                 {
  insertPoint(new TrackPoint(measurement, this), id);
}

◆ 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.

Definition at line 307 of file Track.cc.

                                                 {
 
  point->setTrack(this);
 
  #ifdef DEBUG
  std::cout << "Track::insertPoint at position " << id  << "\n";
  #endif
  assert(point!=NULL);
  trackHasChanged();
 
  point->setTrack(this);
 
  if (trackPoints_.size() == 0) {
    trackPoints_.push_back(point);
 
    if (point->hasRawMeasurements())
      trackPointsWithMeasurement_.push_back(point);
 
    return;
  }
 
  if (id == -1 || id == (int)trackPoints_.size()) {
    trackPoints_.push_back(point);
 
    if (point->hasRawMeasurements())
      trackPointsWithMeasurement_.push_back(point);
 
    deleteReferenceInfo(std::max(0, (int)trackPoints_.size()-2), (int)trackPoints_.size()-1);
 
    // delete fitter infos if inserted point has a measurement
    if (point->hasRawMeasurements()) {
      deleteForwardInfo(-1, -1);
      deleteBackwardInfo(0, -2);
    }
 
    return;
  }
 
  // [-size, size-1] is the allowed range
  assert(id < (ssize_t)trackPoints_.size() || -id-1 <= (ssize_t)trackPoints_.size());
 
  if (id < 0)
    id += trackPoints_.size() + 1;
 
  // insert
  trackPoints_.insert(trackPoints_.begin() + id, point);  // insert inserts BEFORE
 
  // delete fitter infos if inserted point has a measurement
  if (point->hasRawMeasurements()) {
    deleteForwardInfo(id, -1);
    deleteBackwardInfo(0, id);
  }
 
  // delete reference info of neighbouring points
  deleteReferenceInfo(std::max(0, id-1), std::min((int)trackPoints_.size()-1, id+1));
 
  fillPointsWithMeasurement();
}

◆ 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.

Definition at line 367 of file Track.cc.

                                                                    {
 
  int nBefore = getNumPoints();
  int n = points.size();
 
  if (n == 0)
    return;
  if (n == 1) {
    insertPoint(points[0], id);
    return;
  }
 
  for (std::vector<genfit::TrackPoint*>::iterator p = points.begin(); p != points.end(); ++p)
    (*p)->setTrack(this);
 
  if (id == -1 || id == (int)trackPoints_.size()) {
    trackPoints_.insert(trackPoints_.end(), points.begin(), points.end());
 
    deleteReferenceInfo(std::max(0, nBefore-1), nBefore);
 
    deleteForwardInfo(nBefore, -1);
    deleteBackwardInfo(0, std::max(0, nBefore-1));
 
    fillPointsWithMeasurement();
 
    return;
  }
 
 
  assert(id < (ssize_t)trackPoints_.size() || -id-1 <= (ssize_t)trackPoints_.size());
 
  if (id < 0)
    id += trackPoints_.size() + 1;
 
 
  // insert
  trackPoints_.insert(trackPoints_.begin() + id, points.begin(), points.end());  // insert inserts BEFORE
 
  // delete fitter infos if inserted point has a measurement
  deleteForwardInfo(id, -1);
  deleteBackwardInfo(0, id+n);
 
  // delete reference info of neighbouring points
  deleteReferenceInfo(std::max(0, id-1), std::min((int)trackPoints_.size()-1, id));
  deleteReferenceInfo(std::max(0, id+n-1), std::min((int)trackPoints_.size()-1, id+n));
 
  fillPointsWithMeasurement();
}

◆ 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.

Definition at line 457 of file Track.cc.

                                                 {
 
  #ifdef DEBUG
  std::cout << "Track::mergeTrack\n";
  #endif
 
  if (other->getNumPoints() == 0)
    return;
 
  std::map<const AbsTrackRep*, AbsTrackRep*> otherRepThisRep;
  std::vector<const AbsTrackRep*> otherRepsToRemove;
 
  for (std::vector<AbsTrackRep*>::const_iterator otherRep=other->trackReps_.begin(); otherRep!=other->trackReps_.end(); ++otherRep) {
    bool found(false);
    for (std::vector<AbsTrackRep*>::const_iterator thisRep=trackReps_.begin(); thisRep!=trackReps_.end(); ++thisRep) {
      if ((*thisRep)->isSame(*otherRep)) {
        otherRepThisRep[*otherRep] = *thisRep;
        #ifdef DEBUG
        std::cout << " map other rep " << *otherRep << " to " << (*thisRep) << "\n";
        #endif
        if (found) {
          Exception exc("Track::mergeTrack ==> more than one matching rep.",__LINE__,__FILE__);
          exc.setFatal();
          throw exc;
        }
        found = true;
        break;
      }
    }
    if (!found) {
      otherRepsToRemove.push_back(*otherRep);
      #ifdef DEBUG
      std::cout << " remove other rep " << *otherRep << "\n";
      #endif
    }
  }
 
 
  std::vector<TrackPoint*> points;
  points.reserve(other->getNumPoints());
 
  for (std::vector<TrackPoint*>::const_iterator otherTp=other->trackPoints_.begin(); otherTp!=other->trackPoints_.end(); ++otherTp) {
    points.push_back(new TrackPoint(**otherTp, otherRepThisRep, &otherRepsToRemove));
  }
 
  insertPoints(points, id);
}

◆ operator=()

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

Definition at line 145 of file Track.cc.

                                   {
  swap(other);
 
  for (std::vector<TrackPoint*>::const_iterator it=trackPoints_.begin(); it!=trackPoints_.end(); ++it) {
    trackPoints_.back()->setTrack(this);
  }
 
  fillPointsWithMeasurement();
 
  // self test
  assert(checkConsistency());
 
  return *this;
}

◆ Print()

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

Definition at line 1026 of file Track.cc.

                                              {
  TString opt = option;
  opt.ToUpper();
  if (opt.Contains("C")) { // compact
 
    std::cout << "\n    ";
    for (unsigned int i=0; i<trackPoints_.size(); ++i) {
 
      int color = 32*(size_t)(trackPoints_[i]) % 15;
      switch (color) {
        case 0:
          std::cout<<"\033[1;30m";
          break;
        case 1:
          std::cout<<"\033[0;34m";
          break;
        case 2:
          std::cout<<"\033[1;34m";
          break;
        case 3:
          std::cout<<"\033[0;32m";
          break;
        case 4:
          std::cout<<"\033[1;32m";
          break;
        case 5:
          std::cout<<"\033[0;36m";
          break;
        case 6:
          std::cout<<"\033[1;36m";
          break;
        case 7:
          std::cout<<"\033[0;31m";
          break;
        case 8:
          std::cout<<"\033[1;31m";
          break;
        case 9:
          std::cout<<"\033[0;35m";
          break;
        case 10:
          std::cout<<"\033[1;35m";
          break;
        case 11:
          std::cout<<"\033[0;33m";
          break;
        case 12:
          std::cout<<"\033[1;33m";
          break;
        case 13:
          std::cout<<"\033[0;37m";
          break;
        default:
          ;
      }
      std::cout << trackPoints_[i] << "\033[00m  ";
    }
    std::cout << "\n";
 
    std::cout << "   ";
    for (unsigned int i=0; i<trackPoints_.size(); ++i) {
      printf("% -9.3g  ", trackPoints_[i]->getSortingParameter());
    }
 
    for (std::vector<AbsTrackRep*>::const_iterator rep = trackReps_.begin(); rep != trackReps_.end(); ++rep) {
      std::cout << "\n" << getIdForRep(*rep) << "   ";
      for (unsigned int i=0; i<trackPoints_.size(); ++i) {
        if (! trackPoints_[i]->hasFitterInfo(*rep)) {
          std::cout << "           ";
          continue;
        }
        AbsFitterInfo* fi = trackPoints_[i]->getFitterInfo(*rep);
        if (fi->hasMeasurements())
          std::cout << "M";
        else
          std::cout << " ";
 
        if (fi->hasReferenceState())
          std::cout << "R";
        else
          std::cout << " ";
 
        std::cout << "         ";
      }
      std::cout << "\n";
 
      std::cout << " -> ";
      for (unsigned int i=0; i<trackPoints_.size(); ++i) {
        if (! trackPoints_[i]->hasFitterInfo(*rep)) {
          std::cout << "           ";
          continue;
        }
        AbsFitterInfo* fi = trackPoints_[i]->getFitterInfo(*rep);
        if (fi->hasForwardPrediction())
          std::cout << "P";
        else
          std::cout << " ";
 
        if (fi->hasForwardUpdate())
          std::cout << "U";
        else
          std::cout << " ";
 
        std::cout << "         ";
      }
      std::cout << "\n";
 
      std::cout << " <- ";
      for (unsigned int i=0; i<trackPoints_.size(); ++i) {
        if (! trackPoints_[i]->hasFitterInfo(*rep)) {
          std::cout << "           ";
          continue;
        }
        AbsFitterInfo* fi = trackPoints_[i]->getFitterInfo(*rep);
        if (fi->hasBackwardPrediction())
          std::cout << "P";
        else
          std::cout << " ";
 
        if (fi->hasBackwardUpdate())
          std::cout << "U";
        else
          std::cout << " ";
 
        std::cout << "         ";
      }
 
      std::cout << "\n";
 
    } //end loop over reps
 
    std::cout << "\n";
    return;
  }
 
 
 
  std::cout << "=======================================================================================\n";
  std::cout << "genfit::Track, containing " << trackPoints_.size() << " TrackPoints and " << trackReps_.size() << " TrackReps.\n";
  std::cout << " Seed state: "; stateSeed_.Print();
 
  for (unsigned int i=0; i<trackReps_.size(); ++i) {
    std::cout << " TrackRep Nr. " << i;
    if (i == cardinalRep_)
      std::cout << " (This is the cardinal rep)";
    std::cout << "\n";
    trackReps_[i]->Print();
  }
 
  std::cout << "---------------------------------------------------------------------------------------\n";
 
  for (unsigned int i=0; i<trackPoints_.size(); ++i) {
    std::cout << "TrackPoint Nr. " << i << "\n";
    trackPoints_[i]->Print();
    std::cout << "..........................................................................\n";
  }
 
  for (std::map< const AbsTrackRep*, FitStatus* >::const_iterator it=fitStatuses_.begin(); it!=fitStatuses_.end(); ++it) {
    it->second->Print();
  }
 
  std::cout << "=======================================================================================\n";
 
}

◆ prune()

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

Delete unneeded information from the Track.

Possible options: C: prune all reps except cardinalRep F: prune all points except first point L: prune all points except last point FL: prune all points except first and last point 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

Definition at line 952 of file Track.cc.

                                        {
 
  TString opt = option;
  opt.ToUpper();
 
  for (std::map< const AbsTrackRep*, FitStatus* >::const_iterator it=fitStatuses_.begin(); it!=fitStatuses_.end(); ++it) {
    it->second->setIsTrackPruned();
  }
 
  // prune trackPoints
  if (opt.Contains("F") || opt.Contains("L")) {
    TrackPoint* firstPoint = getPointWithMeasurement(0);
    TrackPoint* lastPoint = getPointWithMeasurement(-1);
    for (unsigned int i = 0; i<trackPoints_.size(); ++i) {
      if (trackPoints_[i] == firstPoint && opt.Contains("F"))
        continue;
 
      if (trackPoints_[i] == lastPoint && opt.Contains("L"))
        continue;
 
      delete trackPoints_[i];
      trackPoints_.erase(trackPoints_.begin()+i);
      --i;
    }
  }
 
  // prune TrackReps
  if (opt.Contains("C")) {
    for (unsigned int i = 0; i < trackReps_.size(); ++i) {
      if (i != cardinalRep_) {
        deleteTrackRep(i);
        --i;
      }
    }
  }
 
 
  // from remaining trackPoints: prune measurementsOnPlane, unneeded fitterInfoStuff
  for (unsigned int i = 0; i<trackPoints_.size(); ++i) {
    if (opt.Contains("W"))
      trackPoints_[i]->deleteRawMeasurements();
 
    std::vector< genfit::AbsFitterInfo* > fis =  trackPoints_[i]->getFitterInfos();
    for (unsigned int j = 0; j<fis.size(); ++j) {
 
      if (i == 0 && opt.Contains("I") && opt.Contains("F") && opt.Contains("L"))
        fis[j]->deleteForwardInfo();
      else if (i == trackPoints_.size()-1 && opt.Contains("I") && opt.Contains("F") && opt.Contains("L"))
        fis[j]->deleteBackwardInfo();
      else if (opt.Contains("I") && opt.Contains("F"))
        fis[j]->deleteForwardInfo();
      else if (opt.Contains("I") && opt.Contains("L"))
        fis[j]->deleteBackwardInfo();
 
      if (opt.Contains("U") && dynamic_cast<KalmanFitterInfo*>(fis[j]) != NULL) {
        static_cast<KalmanFitterInfo*>(fis[j])->deletePredictions();
      }
 
      if (opt.Contains("R"))
        fis[j]->deleteReferenceInfo();
      if (opt.Contains("M"))
        fis[j]->deleteMeasurementInfo();
    }
  }
 
  fillPointsWithMeasurement();
 
  #ifdef DEBUG
  std::cout << "pruned Track: "; Print();
  #endif
 
}

◆ reverseMomSeed()

void genfit::Track::reverseMomSeed ( )

inline

Flip direction of momentum seed.

Definition at line 225 of file Track.h.

                        {
    stateSeed_(3) *= -1; stateSeed_(4) *= -1; stateSeed_(5) *= -1;
  }

◆ 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.

Definition at line 673 of file Track.cc.

                         {
  udpateSeed(-1); // set fitted state of last hit as new seed
  reverseMomSeed(); // flip momentum direction
  switchPDGSigns();
  reverseTrackPoints(); // also deletes all fitterInfos
}

◆ reverseTrackPoints()

void genfit::Track::reverseTrackPoints ( )

Flip the ordering of the TrackPoints.

Definition at line 649 of file Track.cc.

                               {
 
  std::reverse(trackPoints_.begin(),trackPoints_.end());
 
  deleteForwardInfo(0, -1);
  deleteBackwardInfo(0, -1);
  deleteReferenceInfo(0, -1);
 
  fillPointsWithMeasurement();
}

◆ setCardinalRep()

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

Definition at line 539 of file Track.cc.

                                 {
 
  if (id < 0)
    id += trackReps_.size();
 
  if (id >= 0 && (unsigned int)id < trackReps_.size())
    cardinalRep_ = id;
  else {
    cardinalRep_ = 0;
    std::cerr << "Track::setCardinalRep: Attempted to set cardinalRep_ to a value out of bounds. Resetting  cardinalRep_ to 0." << std::endl;
  }
}

◆ setCovSeed()

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

inline

Definition at line 163 of file Track.h.

163{covSeed_.ResizeTo(c); covSeed_ = c;}

hnl.c

c

Definition hnl.py:100

◆ setFitStatus()

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

Definition at line 285 of file Track.cc.

                                                                     {
  if (fitStatuses_.find(rep) != fitStatuses_.end())
    delete fitStatuses_.at(rep);
 
  fitStatuses_[rep] = fitStatus;
}

◆ setMcTrackId()

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

inline

Set the MCT track id, for MC simulations.

Definition at line 166 of file Track.h.

166{mcTrackId_ = i;}

◆ setStateSeed() [1/2]

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

Definition at line 293 of file Track.cc.

                                                                 {
  stateSeed_.ResizeTo(6);
 
  stateSeed_(0) = pos.X();
  stateSeed_(1) = pos.Y();
  stateSeed_(2) = pos.Z();
 
  stateSeed_(3) = mom.X();
  stateSeed_(4) = mom.Y();
  stateSeed_(5) = mom.Z();
}

◆ setStateSeed() [2/2]

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

inline

Definition at line 159 of file Track.h.

159{stateSeed_.ResizeTo(s); stateSeed_ = s;}

study_GammaConv.s

s

Definition study_GammaConv.py:22

◆ sort()

bool genfit::Track::sort ( )

Sort TrackPoint and according to their sorting parameters.

Returns if the order of the TrackPoint has actually changed.

Definition at line 578 of file Track.cc.

                 {
  #ifdef DEBUG
  std::cout << "Track::sort \n";
  #endif
 
  int nPoints(trackPoints_.size());
  // original order
  std::vector<TrackPoint*> pointsBefore(trackPoints_);
 
  // sort
  std::stable_sort(trackPoints_.begin(), trackPoints_.end(), TrackPointComparator());
 
  // see where order changed
  int equalUntil(-1), equalFrom(nPoints);
  for (int i = 0; i<nPoints; ++i) {
    if (pointsBefore[i] == trackPoints_[i])
      equalUntil = i;
    else
      break;
  }
 
  if (equalUntil == nPoints-1)
    return false; // sorting did not change anything
 
 
  trackHasChanged();
 
  for (int i = nPoints-1; i>equalUntil; --i) {
    if (pointsBefore[i] == trackPoints_[i])
      equalFrom = i;
    else
      break;
  }
 
  #ifdef DEBUG
  std::cout << "Track::sort. Equal up to (including) hit " << equalUntil << " and from (including) hit " << equalFrom << " \n";
  #endif
 
  deleteForwardInfo(equalUntil+1, -1);
  deleteBackwardInfo(0, equalFrom-1);
 
  deleteReferenceInfo(std::max(0, equalUntil+1), std::min((int)trackPoints_.size()-1, equalFrom-1));
 
  fillPointsWithMeasurement();
 
  return true;
}

◆ swap()

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

Definition at line 160 of file Track.cc.

                             {
  std::swap(this->trackReps_, other.trackReps_);
  std::swap(this->cardinalRep_, other.cardinalRep_);
  std::swap(this->trackPoints_, other.trackPoints_);
  std::swap(this->trackPointsWithMeasurement_, other.trackPointsWithMeasurement_);
  std::swap(this->fitStatuses_, other.fitStatuses_);
  std::swap(this->stateSeed_, other.stateSeed_);
  std::swap(this->covSeed_, other.covSeed_);
 
}

◆ switchPDGSigns()

void genfit::Track::switchPDGSigns ( AbsTrackRep * rep = NULL )

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

Definition at line 661 of file Track.cc.

                                           {
  if (rep != NULL) {
    rep->switchPDGSign();
    return;
  }
 
  for (unsigned int i = 0; i<trackReps_.size(); ++i) {
    trackReps_[i]->switchPDGSign();
  }
}

◆ trackHasChanged()

void genfit::Track::trackHasChanged ( )

private

Definition at line 1341 of file Track.cc.

                            {
 
  #ifdef DEBUG
  std::cout << "Track::trackHasChanged \n";
  #endif
 
  if (fitStatuses_.empty())
    return;
 
  for (std::map< const AbsTrackRep*, FitStatus* >::const_iterator it=fitStatuses_.begin(); it!=fitStatuses_.end(); ++it) {
    it->second->setHasTrackChanged();
  }
}

◆ udpateSeed()

bool genfit::Track::udpateSeed	(	int	id = `0`,
		AbsTrackRep *	rep = `NULL`,
		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.

Definition at line 627 of file Track.cc.

                                                            {
  try {
    const MeasuredStateOnPlane& fittedState = getFittedState(id, rep, biased);
    setStateSeed(fittedState.get6DState());
    setCovSeed(fittedState.get6DCov());
 
    double fittedCharge = fittedState.getCharge();
 
    for (unsigned int i = 0; i<trackReps_.size(); ++i) {
      if (trackReps_[i]->getPDGCharge() * fittedCharge < 0) {
        trackReps_[i]->switchPDGSign();
      }
    }
  }
  catch (Exception& e) {
    // in this case the original track seed will be used
    return false;
  }
  return true;
}

Member Data Documentation

◆ cardinalRep_

unsigned int genfit::Track::cardinalRep_

private

Definition at line 293 of file Track.h.

◆ covSeed_

TMatrixDSym genfit::Track::covSeed_

private

Definition at line 302 of file Track.h.

◆ fitStatuses_

std::map< const AbsTrackRep*, FitStatus* > genfit::Track::fitStatuses_

private

helper

Definition at line 298 of file Track.h.

◆ mcTrackId_

int genfit::Track::mcTrackId_

private

if MC simulation, store the mc track id here

Definition at line 300 of file Track.h.

◆ stateSeed_

TVectorD genfit::Track::stateSeed_

private

Definition at line 301 of file Track.h.

◆ trackPoints_

std::vector<TrackPoint*> genfit::Track::trackPoints_

private

Definition at line 295 of file Track.h.

◆ trackPointsWithMeasurement_

std::vector<TrackPoint*> genfit::Track::trackPointsWithMeasurement_

private

Definition at line 296 of file Track.h.

◆ trackReps_

std::vector<AbsTrackRep*> genfit::Track::trackReps_

private

Definition at line 292 of file Track.h.

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

genfit/core/include/Track.h
genfit/core/src/Track.cc

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ Track() [1/6]

◆ Track() [2/6]

◆ Track() [3/6]

◆ Track() [4/6]

◆ Track() [5/6]

◆ Track() [6/6]

◆ ~Track()

Member Function Documentation

◆ addTrackRep()

◆ checkConsistency()

◆ Clear()

◆ constructTrackCand()

◆ deleteBackwardInfo()

◆ deleteFitterInfo()

◆ deleteForwardInfo()

◆ deleteMeasurementInfo()

◆ deletePoint()

◆ deleteReferenceInfo()

◆ deleteTrackRep()

◆ determineCardinalRep()

◆ fillPointsWithMeasurement()

◆ fixWeights()

◆ getCardinalRep()

◆ getCardinalRepId()

◆ getCovSeed()

◆ getFitStatus()

◆ getFittedState()

◆ getIdForRep()

◆ getKalmanFitStatus()

◆ getMcTrackId()

◆ getNumPoints()

◆ getNumPointsWithMeasurement()

◆ getNumReps()

◆ getPoint()

◆ getPoints()

◆ getPointsWithMeasurement()

◆ getPointWithMeasurement()

◆ getPointWithMeasurementAndFitterInfo()

◆ getStateSeed()

◆ getTOF()

◆ getTrackLen()

◆ getTrackRep()

◆ hasFitStatus()

◆ hasKalmanFitStatus()

◆ insertMeasurement()

◆ insertPoint()

◆ insertPoints()

◆ mergeTrack()

◆ operator=()

◆ Print()

◆ prune()

◆ reverseMomSeed()

◆ reverseTrack()

◆ reverseTrackPoints()

◆ setCardinalRep()

◆ setCovSeed()

◆ setFitStatus()

◆ setMcTrackId()

◆ setStateSeed() [1/2]

◆ setStateSeed() [2/2]

◆ sort()

◆ swap()

◆ switchPDGSigns()

◆ trackHasChanged()

◆ udpateSeed()

Member Data Documentation

◆ cardinalRep_

◆ covSeed_

◆ fitStatuses_

◆ mcTrackId_

◆ stateSeed_

◆ trackPoints_

◆ trackPointsWithMeasurement_

◆ trackReps_