Base Class for genfit track representations. Defines interface for track parameterizations. More...

#include <GFAbsTrackRep.h>

Inheritance diagram for genf::GFAbsTrackRep:

Public Member Functions
virtual GFAbsTrackRep *	clone () const =0

virtual GFAbsTrackRep *	prototype () const =0

virtual double	extrapolate (const GFDetPlane &plane, TMatrixT< Double_t > &statePred)
	returns the tracklength spanned in this extrapolation More...

	GFAbsTrackRep ()

	GFAbsTrackRep (int)

virtual	~GFAbsTrackRep ()

virtual void	extrapolateToPoint (const TVector3 &point, TVector3 &poca, TVector3 &normVec)
	This method is to extrapolate the track to point of closest approach to a point in space. More...

virtual void	extrapolateToLine (const TVector3 &point1, const TVector3 &point2, TVector3 &poca, TVector3 &normVec, TVector3 &poca_onwire)
	This method extrapolates to the point of closest approach to a line. More...

virtual void	stepalong (double h)
	make step of h cm along the track More...

virtual double	extrapolate (const GFDetPlane &plane, TMatrixT< Double_t > &statePred, TMatrixT< Double_t > &covPred)=0
	Extrapolates the track to the given detectorplane. More...

double	extrapolate (const GFDetPlane &plane)
	This changes the state and cov and plane of the rep. More...

unsigned int	getDim () const
	returns dimension of state vector More...

virtual void	Print (std::ostream &out=std::cout) const

const TMatrixT< Double_t > &	getState () const

const TMatrixT< Double_t > &	getCov () const

double	getStateElem (int i) const

double	getCovElem (int i, int j) const

virtual TVector3	getPos (const GFDetPlane &pl)=0

virtual TVector3	getMom (const GFDetPlane &pl)=0

virtual void	getPosMom (const GFDetPlane &pl, TVector3 &pos, TVector3 &mom)=0

virtual void	getPosMomCov (const GFDetPlane &pl, TVector3 &pos, TVector3 &mom, TMatrixT< Double_t > &cov)
	method which gets position, momentum and 6x6 covariance matrix More...

virtual double	getCharge () const =0

TVector3	getPos ()

TVector3	getMom ()

void	getPosMomCov (TVector3 &pos, TVector3 &mom, TMatrixT< Double_t > &c)

TMatrixT< Double_t >	getFirstState () const

TMatrixT< Double_t >	getFirstCov () const

GFDetPlane	getFirstPlane () const

TMatrixT< Double_t >	getLastState () const

TMatrixT< Double_t >	getLastCov () const

GFDetPlane	getLastPlane () const

double	getChiSqu () const

double	getRedChiSqu () const
	returns chi2/ndf More...

unsigned int	getNDF () const

virtual void	setData (const TMatrixT< Double_t > &st, const GFDetPlane &pl, const TMatrixT< Double_t > *cov=NULL)

void	setCov (const TMatrixT< Double_t > &aCov)

void	setFirstState (const TMatrixT< Double_t > &aState)

void	setFirstCov (const TMatrixT< Double_t > &aCov)

void	setFirstPlane (const GFDetPlane &aPlane)

void	setLastState (const TMatrixT< Double_t > &aState)

void	setLastCov (const TMatrixT< Double_t > &aCov)

void	setLastPlane (const GFDetPlane &aPlane)

const GFDetPlane &	getReferencePlane () const

void	setChiSqu (double aChiSqu)

void	setNDF (unsigned int n)

void	addChiSqu (double aChiSqu)

void	addNDF (unsigned int n)

void	setStatusFlag (int _val)

virtual void	switchDirection ()=0

bool	setInverted (bool f=true)
	Deprecated. Should be removed soon. More...

bool	getStatusFlag ()

virtual void	reset ()

Protected Attributes
unsigned int	fDimension
	Dimensionality of track representation. More...

TMatrixT< Double_t >	fState
	The vector of track parameters. More...

TMatrixT< Double_t >	fCov
	The covariance matrix. More...

double	fChiSqu
	chiSqu of the track fit More...

unsigned int	fNdf

int	fStatusFlag
	status of track representation: 0 means everything's OK More...

bool	fInverted
	specifies the direction of flight of the particle More...

TMatrixT< Double_t >	fFirstState
	state, cov and plane for first and last point in fit More...

TMatrixT< Double_t >	fFirstCov

TMatrixT< Double_t >	fLastState

TMatrixT< Double_t >	fLastCov

GFDetPlane	fFirstPlane

GFDetPlane	fLastPlane

GFDetPlane	fRefPlane

Private Member Functions
void	Abort (std::string method)

virtual void	Print (Option_t *) const

Detailed Description

Base Class for genfit track representations. Defines interface for track parameterizations.

Author: Christian Höppner (Technische Universität München, original author); Sebastian Neubert (Technische Universität München, original author)

It is important to understand the difference between a track and a track representation in genfit:

A track representation is a specific parameterization of a trajectory. It contains the parameters that describe the track at some point and code for the extrapolation of the track parameters through space. The actual extrapolation code is not part of genfit but has to be supplied in some additional package (e.g. GEANE). LSLTrackRep is a very basic example of a track representation.
A Track is a collection of RecoHits (see GFAbsRecoHit) plus a collection of track representation objects. The hits can be from different detectors. There can be several representations of the same track. This makes it possible to perform several fits in parallel, for example to compare different parameterizations or to fit different particle hypotheses.

All track tepresentations must inherit GFAbsTrackRep to be available in genfit. Algorithms in genfit use this class as interface to access track parameters

Provides:

Matrix objects to store track parameters
... and covariances
interface to track extrapolation code

The track extrapolation engine can be exchanged in genfit. Or one can even use more than one engine in parallel! In order to use a track extrapolation engine (like e.g. GEANE) with genfit one has to write a TrackRep class that inherits from GFAbsTrackRep. This makes it possible to uses different track extrapolation codes within a unified framework without major changes in the detector code.

There is only one thing one has to do to use a specific track representation together with the hits from a detector: add the respective code in the GFAbsRecoHit::getHMatrix method implementation of the RecoHit in question.

Definition at line 83 of file GFAbsTrackRep.h.

Constructor & Destructor Documentation

genf::GFAbsTrackRep::GFAbsTrackRep ( )

Definition at line 23 of file GFAbsTrackRep.cxx.

                                  : fDimension(5),fState(5,1), fCov(5,5), fChiSqu(0), fNdf(0), fStatusFlag(0), fInverted(false), fFirstState(5,1), fFirstCov(5,5), fLastState(5,1), fLastCov(5,5)
 {
 }

genf::GFAbsTrackRep::GFAbsTrackRep ( int dim )

Definition at line 27 of file GFAbsTrackRep.cxx.

                                         : fDimension(dim), fState(dim,1), fCov(dim,dim), fChiSqu(0), fNdf(0), fStatusFlag(0), fInverted(false), fFirstState(dim,1), fFirstCov(dim,dim), fLastState(dim,1), fLastCov(dim,dim)
 {
 }

genf::GFAbsTrackRep::~GFAbsTrackRep ( )

virtual

Definition at line 31 of file GFAbsTrackRep.cxx.

31 {}

Member Function Documentation

void genf::GFAbsTrackRep::Abort ( std::string method )

private

Definition at line 47 of file GFAbsTrackRep.cxx.

                                              {
   std::cerr << method <<  " as implemented in " << __FILE__
             << " was called. This means that this feature was used "
             << "in a track rep which didnt overwrite this method. "
             << std::endl << "C++ throw;" << std::endl;
   //system call abort
   throw GFException("genf::GFAbsTrackRep: " + method + "() not implemented", __LINE__, __FILE__).setFatal();
 }

void genf::GFAbsTrackRep::addChiSqu ( double aChiSqu )

inline

Definition at line 291 of file GFAbsTrackRep.h.

                                         {
     fChiSqu += aChiSqu;
   }

void genf::GFAbsTrackRep::addNDF ( unsigned int n )

inline

Definition at line 294 of file GFAbsTrackRep.h.

                                      {
     fNdf += n;
   }

virtual GFAbsTrackRep* genf::GFAbsTrackRep::clone ( ) const

pure virtual

Implemented in genf::RKTrackRep, and genf::SlTrackRep.

double genf::GFAbsTrackRep::extrapolate	(	const GFDetPlane &	plane,
		TMatrixT< Double_t > &	statePred
	)

virtual

returns the tracklength spanned in this extrapolation

Reimplemented in genf::RKTrackRep.

Definition at line 42 of file GFAbsTrackRep.cxx.

                                                                                            {
   TMatrixT<Double_t> cov(fDimension,fDimension);
   return extrapolate(plane,statePred,cov);
 }

virtual double genf::GFAbsTrackRep::extrapolate	(	const GFDetPlane &	plane,
		TMatrixT< Double_t > &	statePred,
		TMatrixT< Double_t > &	covPred
	)

pure virtual

Extrapolates the track to the given detectorplane.

Results are put into statePred and covPred This method does NOT alter the state of the object!

Implemented in genf::RKTrackRep.

double genf::GFAbsTrackRep::extrapolate ( const GFDetPlane & plane )

This changes the state and cov and plane of the rep.

This method extrapolates to to the plane and sets the results of state, cov and also plane in itself.

Definition at line 33 of file GFAbsTrackRep.cxx.

                                                             {
   TMatrixT<Double_t> statePred(fDimension,1);
   TMatrixT<Double_t> covPred(fDimension,fDimension);
   double retVal = extrapolate(plane,statePred,covPred);
   setData(statePred,plane,&covPred);
   return retVal;
 }

void genf::GFAbsTrackRep::extrapolateToLine	(	const TVector3 &	point1,
		const TVector3 &	point2,
		TVector3 &	poca,
		TVector3 &	normVec,
		TVector3 &	poca_onwire
	)

virtual

This method extrapolates to the point of closest approach to a line.

This method extrapolates to the POCA to a line, i.e. a wire. There is a default implementation just like for the extrapolateToPoca for trackReps which do not need this feature, which will abort the execution if it is ever called.

Reimplemented in genf::RKTrackRep, and genf::SlTrackRep.

Definition at line 62 of file GFAbsTrackRep.cxx.

                                                                                     {
   Abort("extrapolateToLine()");
 }

void genf::GFAbsTrackRep::extrapolateToPoint	(	const TVector3 &	point,
		TVector3 &	poca,
		TVector3 &	normVec
	)

virtual

This method is to extrapolate the track to point of closest approach to a point in space.

Reimplemented in genf::RKTrackRep, and genf::SlTrackRep.

Definition at line 56 of file GFAbsTrackRep.cxx.

                                                 {
   Abort("extrapolateToPoca()");
 }

virtual double genf::GFAbsTrackRep::getCharge ( ) const

pure virtual

Implemented in genf::RKTrackRep, and genf::SlTrackRep.

double genf::GFAbsTrackRep::getChiSqu ( ) const

inline

Definition at line 243 of file GFAbsTrackRep.h.

                                   {
     return fChiSqu;
   }

const TMatrixT<Double_t>& genf::GFAbsTrackRep::getCov ( ) const

inline

Definition at line 197 of file GFAbsTrackRep.h.

197 { return fCov; }

genf::GFAbsTrackRep::fCov

TMatrixT< Double_t > fCov

The covariance matrix.

Definition: GFAbsTrackRep.h:94

double genf::GFAbsTrackRep::getCovElem	(	int	i,
		int	j
	)		const

inline

Definition at line 200 of file GFAbsTrackRep.h.

200 {return fCov(i,j);}

genf::GFAbsTrackRep::fCov

TMatrixT< Double_t > fCov

The covariance matrix.

Definition: GFAbsTrackRep.h:94

unsigned int genf::GFAbsTrackRep::getDim ( ) const

inline

returns dimension of state vector

Definition at line 192 of file GFAbsTrackRep.h.

192 {return fDimension;}

genf::GFAbsTrackRep::fDimension

unsigned int fDimension

Dimensionality of track representation.

Definition: GFAbsTrackRep.h:88

TMatrixT<Double_t> genf::GFAbsTrackRep::getFirstCov ( ) const

inline

Definition at line 228 of file GFAbsTrackRep.h.

                                                 {
     return fFirstCov;
   }

GFDetPlane genf::GFAbsTrackRep::getFirstPlane ( ) const

inline

Definition at line 231 of file GFAbsTrackRep.h.

                                           {
     return fFirstPlane;
   }

TMatrixT<Double_t> genf::GFAbsTrackRep::getFirstState ( ) const

inline

Definition at line 225 of file GFAbsTrackRep.h.

                                                   {
     return fFirstState;
   }

TMatrixT<Double_t> genf::GFAbsTrackRep::getLastCov ( ) const

inline

Definition at line 237 of file GFAbsTrackRep.h.

                                                {
     return fLastCov;
   }

GFDetPlane genf::GFAbsTrackRep::getLastPlane ( ) const

inline

Definition at line 240 of file GFAbsTrackRep.h.

                                          {
     return fLastPlane;
   }

TMatrixT<Double_t> genf::GFAbsTrackRep::getLastState ( ) const

inline

Definition at line 234 of file GFAbsTrackRep.h.

                                                  {
     return fLastState;
   }

virtual TVector3 genf::GFAbsTrackRep::getMom ( const GFDetPlane & pl )

pure virtual

Implemented in genf::RKTrackRep, and genf::SlTrackRep.

TVector3 genf::GFAbsTrackRep::getMom ( )

inline

Definition at line 219 of file GFAbsTrackRep.h.

219 {return getMom(fRefPlane);}

genf::GFAbsTrackRep::getMom

TVector3 getMom()

Definition: GFAbsTrackRep.h:219

genf::GFAbsTrackRep::fRefPlane

GFDetPlane fRefPlane

Definition: GFAbsTrackRep.h:115

unsigned int genf::GFAbsTrackRep::getNDF ( ) const

inline

Definition at line 251 of file GFAbsTrackRep.h.

                                      {
     if(fNdf>getDim())  return fNdf-getDim();
     return 0;
   }

virtual TVector3 genf::GFAbsTrackRep::getPos ( const GFDetPlane & pl )

pure virtual

Implemented in genf::RKTrackRep, and genf::SlTrackRep.

TVector3 genf::GFAbsTrackRep::getPos ( )

inline

Definition at line 218 of file GFAbsTrackRep.h.

218 {return getPos(fRefPlane);}

genf::GFAbsTrackRep::getPos

TVector3 getPos()

Definition: GFAbsTrackRep.h:218

genf::GFAbsTrackRep::fRefPlane

GFDetPlane fRefPlane

Definition: GFAbsTrackRep.h:115

virtual void genf::GFAbsTrackRep::getPosMom	(	const GFDetPlane &	pl,
		TVector3 &	pos,
		TVector3 &	mom
	)

pure virtual

Implemented in genf::RKTrackRep, and genf::SlTrackRep.

void genf::GFAbsTrackRep::getPosMomCov	(	const GFDetPlane &	pl,
		TVector3 &	pos,
		TVector3 &	mom,
		TMatrixT< Double_t > &	cov
	)

virtual

method which gets position, momentum and 6x6 covariance matrix

default implementation in cxx file, if a ConcreteTrackRep can not implement this functionality

Definition at line 75 of file GFAbsTrackRep.cxx.

                                                                                                      {
   Abort("getPosMomCov()");
 }

void genf::GFAbsTrackRep::getPosMomCov	(	TVector3 &	pos,
		TVector3 &	mom,
		TMatrixT< Double_t > &	c
	)

inline

Definition at line 221 of file GFAbsTrackRep.h.

                                                                       {
     getPosMomCov(fRefPlane,pos,mom,c);
   }

double genf::GFAbsTrackRep::getRedChiSqu ( ) const

inline

returns chi2/ndf

Definition at line 247 of file GFAbsTrackRep.h.

                                      {
     if(getNDF()>0) return getChiSqu()/getNDF();
     return 0;
   }

const GFDetPlane& genf::GFAbsTrackRep::getReferencePlane ( ) const

inline

Definition at line 283 of file GFAbsTrackRep.h.

283 {return fRefPlane;}

genf::GFAbsTrackRep::fRefPlane

GFDetPlane fRefPlane

Definition: GFAbsTrackRep.h:115

const TMatrixT<Double_t>& genf::GFAbsTrackRep::getState ( ) const

inline

Definition at line 196 of file GFAbsTrackRep.h.

196 { return fState; }

genf::GFAbsTrackRep::fState

TMatrixT< Double_t > fState

The vector of track parameters.

Definition: GFAbsTrackRep.h:91

double genf::GFAbsTrackRep::getStateElem ( int i ) const

inline

Definition at line 199 of file GFAbsTrackRep.h.

199 {return fState(i,0);}

genf::GFAbsTrackRep::fState

TMatrixT< Double_t > fState

The vector of track parameters.

Definition: GFAbsTrackRep.h:91

bool genf::GFAbsTrackRep::getStatusFlag ( )

inline

Definition at line 306 of file GFAbsTrackRep.h.

                               {
     return fStatusFlag;
   }

void genf::GFAbsTrackRep::Print ( std::ostream & out = std::cout ) const

virtual

Definition at line 93 of file GFAbsTrackRep.cxx.

                                                 {
   out << "++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++" << std::endl;
   out <<"GFAbsTrackRep::Parameters at reference plane ";
   fRefPlane.Print(out);
   out <<"GFAbsTrackRep::State"<<std::endl;
   PrintROOTmatrix(out, fState);
   out <<"GFAbsTrackRep::Covariances"<<std::endl;
   PrintROOTmatrix(out, fCov);
   out <<"GFAbsTrackRep::chi^2"<<std::endl;
   out <<fChiSqu<<std::endl;
   out << "++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++" << std::endl;
 }

virtual void genf::GFAbsTrackRep::Print ( Option_t * ) const

inlineprivatevirtual

Definition at line 315 of file GFAbsTrackRep.h.

316 { throw std::logic_error(std::string(__func__) + "::Print(Option_t*) not available"); }

virtual GFAbsTrackRep* genf::GFAbsTrackRep::prototype ( ) const

pure virtual

Implemented in genf::RKTrackRep, and genf::SlTrackRep.

void genf::GFAbsTrackRep::reset ( )

virtual

Definition at line 80 of file GFAbsTrackRep.cxx.

                         {
   std::cout<<"GFAbsTrackRep::reset"<<std::endl;
   TVector3 nullVec(0.,0.,0.);
   fRefPlane.set(nullVec,nullVec,nullVec);
   fState.Zero();
   fCov.Zero();
   fFirstState.Zero();
   fFirstCov.Zero();
   fLastState.Zero();
   fLastCov.Zero();
 }

void genf::GFAbsTrackRep::setChiSqu ( double aChiSqu )

inline

Definition at line 285 of file GFAbsTrackRep.h.

                                         {
     fChiSqu = aChiSqu;
   }

void genf::GFAbsTrackRep::setCov ( const TMatrixT< Double_t > & aCov )

inline

Definition at line 261 of file GFAbsTrackRep.h.

                                                      {
     fCov=aCov;
   }

virtual void genf::GFAbsTrackRep::setData	(	const TMatrixT< Double_t > &	st,
		const GFDetPlane &	pl,
		const TMatrixT< Double_t > *	cov = `NULL`
	)

inlinevirtual

Definition at line 256 of file GFAbsTrackRep.h.

                                                                                                               {
     fState=st;
     fRefPlane=pl;
     if(cov!=NULL) fCov=*cov;
   }

void genf::GFAbsTrackRep::setFirstCov ( const TMatrixT< Double_t > & aCov )

inline

Definition at line 267 of file GFAbsTrackRep.h.

                                                           {
     fFirstCov = aCov;
   }

void genf::GFAbsTrackRep::setFirstPlane ( const GFDetPlane & aPlane )

inline

Definition at line 270 of file GFAbsTrackRep.h.

                                                       {
     fFirstPlane = aPlane;;
   }

void genf::GFAbsTrackRep::setFirstState ( const TMatrixT< Double_t > & aState )

inline

Definition at line 264 of file GFAbsTrackRep.h.

                                                               {
     fFirstState = aState;
   }

bool genf::GFAbsTrackRep::setInverted ( bool f = true )

inline

Deprecated. Should be removed soon.

Definition at line 304 of file GFAbsTrackRep.h.

304 {fInverted=f; return true;}

genf::GFAbsTrackRep::fInverted

bool fInverted

specifies the direction of flight of the particle

Definition: GFAbsTrackRep.h:103

void genf::GFAbsTrackRep::setLastCov ( const TMatrixT< Double_t > & aCov )

inline

Definition at line 276 of file GFAbsTrackRep.h.

                                                          {
     fLastCov = aCov;
   }

void genf::GFAbsTrackRep::setLastPlane ( const GFDetPlane & aPlane )

inline

Definition at line 279 of file GFAbsTrackRep.h.

                                                      {
     fLastPlane = aPlane;;
   }

void genf::GFAbsTrackRep::setLastState ( const TMatrixT< Double_t > & aState )

inline

Definition at line 273 of file GFAbsTrackRep.h.

                                                              {
     fLastState = aState;
   }

void genf::GFAbsTrackRep::setNDF ( unsigned int n )

inline

Definition at line 288 of file GFAbsTrackRep.h.

                                      {
     fNdf = n;
   }

void genf::GFAbsTrackRep::setStatusFlag ( int _val )

inline

Definition at line 297 of file GFAbsTrackRep.h.

                                       {
     fStatusFlag = _val;
   }

void genf::GFAbsTrackRep::stepalong ( double h )

virtual

make step of h cm along the track

There is an emply implementation in GFAbsTrackRep.cxx which will abort (see one of the extrapolate methods above). This can be overwritten, if this feature is needed.

Definition at line 71 of file GFAbsTrackRep.cxx.

                                          {
   Abort("stepalong()");
 }

virtual void genf::GFAbsTrackRep::switchDirection ( )

pure virtual

Implemented in genf::RKTrackRep, and genf::SlTrackRep.

Member Data Documentation

double genf::GFAbsTrackRep::fChiSqu

protected

chiSqu of the track fit

Definition at line 97 of file GFAbsTrackRep.h.

TMatrixT<Double_t> genf::GFAbsTrackRep::fCov

protected

The covariance matrix.

Definition at line 94 of file GFAbsTrackRep.h.

unsigned int genf::GFAbsTrackRep::fDimension

protected

Dimensionality of track representation.

Definition at line 88 of file GFAbsTrackRep.h.

TMatrixT<Double_t> genf::GFAbsTrackRep::fFirstCov

protected

Definition at line 107 of file GFAbsTrackRep.h.

GFDetPlane genf::GFAbsTrackRep::fFirstPlane

protected

Definition at line 111 of file GFAbsTrackRep.h.

TMatrixT<Double_t> genf::GFAbsTrackRep::fFirstState

protected

state, cov and plane for first and last point in fit

Definition at line 106 of file GFAbsTrackRep.h.

bool genf::GFAbsTrackRep::fInverted

protected

specifies the direction of flight of the particle

Definition at line 103 of file GFAbsTrackRep.h.

TMatrixT<Double_t> genf::GFAbsTrackRep::fLastCov

protected

Definition at line 110 of file GFAbsTrackRep.h.

GFDetPlane genf::GFAbsTrackRep::fLastPlane

protected

Definition at line 112 of file GFAbsTrackRep.h.

TMatrixT<Double_t> genf::GFAbsTrackRep::fLastState

protected

Definition at line 109 of file GFAbsTrackRep.h.

unsigned int genf::GFAbsTrackRep::fNdf

protected

Definition at line 98 of file GFAbsTrackRep.h.

GFDetPlane genf::GFAbsTrackRep::fRefPlane

protected

Definition at line 115 of file GFAbsTrackRep.h.

TMatrixT<Double_t> genf::GFAbsTrackRep::fState

protected

The vector of track parameters.

Definition at line 91 of file GFAbsTrackRep.h.

int genf::GFAbsTrackRep::fStatusFlag

protected

status of track representation: 0 means everything's OK

Definition at line 101 of file GFAbsTrackRep.h.

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

Public Member Functions

Protected Attributes

Private Member Functions

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation