genf::GFTrack Class Reference

#include <GFTrack.h>

Inheritance diagram for genf::GFTrack:

Public Member Functions
	GFTrack ()
	Default constructor – needed for compatibility with ROOT. More...
	GFTrack (const GFTrack &)
	Copy constructor. More...
GFTrack &	operator= (const GFTrack &)
	assignement operator More...
	GFTrack (GFAbsTrackRep *)
	Initializing constructor. More...
virtual	~GFTrack ()
void	reset ()
	Resets the GFTrack – deletes RecoHits! More...
int	getFailedHits (int repId=-1)
	return the number of failed Hits in track fit repId == -1 will use cardinal rep More...
std::vector< GFAbsRecoHit * >	getHits ()
const GFTrackCand &	getCand () const
GFAbsRecoHit *	getHit (int id) const
unsigned int	getNumHits () const
void	mergeHits (GFTrack *trk)
	Merge two GFTracks. Only hits will be merged. More...
void	releaseHits ()
	Clear hit vector. Note that hits will not be deleted! More...
unsigned int	getNextHitToFit () const
	Accessor for fNextHitToFit. More...
void	setNextHitToFit (unsigned int i)
	Set next hit to be used in a fit. More...
GFAbsTrackRep *	getTrackRep (int id) const
	Accessor for track representations. More...
unsigned int	getNumReps () const
	Get number of track represenatations. More...
GFAbsTrackRep *	getCardinalRep () const
	Get cardinal track representation. More...
TVector3	getMom () const
	Get momentum at the present position. More...
TVector3	getMom (const GFDetPlane &pl) const
	Get momentum at GFDetPlane. More...
TVector3	getPos () const
	Get present position. More...
TVector3	getPos (const GFDetPlane &pl) const
	Get position at GFDetPlane. More...
void	getPosMomCov (TVector3 &pos, TVector3 &mom, TMatrixT< Double_t > &cov)
	Get position, momentum, and 6x6 covariance at current position. More...
void	getPosMomCov (const GFDetPlane &pl, TVector3 &pos, TVector3 &mom, TMatrixT< Double_t > &cov)
	Get position, momentum, and 6x6 covariance at GFDetPlane. More...
double	getChiSqu () const
	Get chi2. More...
unsigned int	getNDF () const
	Get NDF. More...
double	getRedChiSqu () const
	Get chi2/NDF. More...
double	getCharge () const
	Get charge from fit. More...
void	fillGeoTrack (TVirtualGeoTrack *tr) const
	Fill TVirtualGeoTrack object Cardinal representation is used. More...
void	fillGeoTrack (TVirtualGeoTrack *tr, unsigned int repid) const
	Fill TVirtualGeoTrack object with data from specified track rep. More...
void	addFailedHit (unsigned int irep, unsigned int id)
void	addHit (GFAbsRecoHit *theHit)
	deprecated! More...
void	addHit (GFAbsRecoHit *theHit, unsigned int detId, unsigned int hitId, double rho=0., unsigned int planeId=0)
	Add single hit. Updates the GFTrackCand. More...
void	addHitVector (std::vector< GFAbsRecoHit * > hits)
	Add collection of hits. More...
void	addTrackRep (GFAbsTrackRep *theTrackRep)
	Add track represenation. More...
GFBookkeeping *	getBK (int index=-1)
	get GFBookKeeping object for particular track rep (default is cardinal rep) More...
void	setCandidate (const GFTrackCand &cand, bool doreset=false)
	set track candidate More...
void	setCardinalRep (unsigned int r)
	Choose cardinal track represenatation. More...
void	setHitMeasuredCov (TMatrixT< Double_t > mat)
void	setHitUpdate (TMatrixT< Double_t > mat)
void	setHitChi2 (Double_t mat)
void	setHitState (TMatrixT< Double_t > mat)
void	setHitCov (TMatrixT< Double_t > mat)
void	setHitCov7x7 (TMatrixT< Double_t > mat)
void	setHitPlaneXYZ (TVector3 pl)
void	setHitPlaneUxUyUz (TVector3 pl)
void	setHitPlaneU (TVector3 pl)
void	setHitPlaneV (TVector3 pl)
void	setPDG (int pdgt)
std::vector< TMatrixT< Double_t > >	getHitMeasuredCov ()
std::vector< TMatrixT< Double_t > >	getHitUpdate ()
std::vector< Double_t >	getHitChi2 ()
std::vector< TMatrixT< Double_t > >	getHitState ()
std::vector< TMatrixT< Double_t > >	getHitCov ()
std::vector< TMatrixT< Double_t > >	getHitCov7x7 ()
std::vector< TVector3 >	getHitPlaneXYZ ()
std::vector< TVector3 >	getHitPlaneUxUyUz ()
std::vector< TVector3 >	getHitPlaneU ()
std::vector< TVector3 >	getHitPlaneV ()
int	getPDG ()
void	getResiduals (unsigned int detId, unsigned int dim, unsigned int rep, std::vector< double > &result)
	Get residuals. More...
void	setRepAtHit (unsigned int irep, int ihit)
	set the hit index at which plane,state&cov of rep irep is defined More...
int	getRepAtHit (unsigned int irep)
	get the hit index at which plane,state&cov of rep irep is defined More...
void	clearRepAtHit ()
	clear the hit indices at which plane,state&cov of reps are defined More...
void	printBookkeeping (std::ostream &out=std::cout) const
	print bookkeeping More...
void	Print (std::ostream &out=std::cout) const
void	clearBookkeeping ()
void	clearFailedHits ()
void	getHitsByPlane (std::vector< std::vector< int > * > &retVal)

Private Member Functions
virtual void	Print (Option_t *) const

Private Attributes
TObjArray *	fTrackReps
	Collection of track representations. More...
std::vector< GFAbsRecoHit * >	fHits
	Collection of RecoHits. More...
std::vector< GFBookkeeping * >	fBookkeeping
	Collection of Bookeeping objects for failed hits in every trackrep. More...
std::vector< int >	fRepAtHit
	repAtHit keeps track of at which hit index which rep is currently defined, to avoid null extrapolations More...
GFTrackCand	fCand
	Helper to store the indices of the hits in the track. See GFTrackCand for details. More...
std::vector< TMatrixT< Double_t > >	fHitMeasCov
std::vector< TMatrixT< Double_t > >	fHitUpdate
std::vector< TMatrixT< Double_t > >	fHitState
std::vector< Double_t >	fHitChi2
std::vector< TMatrixT< Double_t > >	fHitCov7x7
std::vector< TMatrixT< Double_t > >	fHitCov
std::vector< TVector3 >	fHitPlaneXYZ
std::vector< TVector3 >	fHitPlaneUxUyUz
std::vector< TVector3 >	fHitPlaneU
std::vector< TVector3 >	fHitPlaneV
int	fPDG
unsigned int	fCardinal_rep
unsigned int	fNextHitToFit

Static Private Attributes
static const int	fDefNumTrackReps = 10

Detailed Description

Definition at line 64 of file GFTrack.h.

Constructor & Destructor Documentation

genf::GFTrack::GFTrack

(

)

Default constructor – needed for compatibility with ROOT.

Definition at line 31 of file GFTrack.cxx.

  : fTrackReps(NULL), fCardinal_rep(0), fNextHitToFit(0)
{
  //trackReps = new TObjArray(defNumTrackReps);
}

genf::GFTrack::GFTrack

(

const GFTrack &

_tr

)

Copy constructor.

Definition at line 52 of file GFTrack.cxx.

                                      : TObject() {
  fCand=_tr.fCand;
  fCardinal_rep=_tr.fCardinal_rep;
  fNextHitToFit=_tr.fNextHitToFit;
  for(unsigned int i=0;i<_tr.getNumHits();i++) {
    fHits.push_back((_tr.getHit(i))->clone());
  }
  fTrackReps = NULL;
  for(unsigned int i=0; i<_tr.getNumReps();i++) {
    addTrackRep( (_tr.getTrackRep(i))->clone() );
  }
  for(unsigned int i=0; i<fBookkeeping.size(); ++i) delete fBookkeeping[i];
  fBookkeeping.clear();

  for(unsigned int i=0;i<_tr.fBookkeeping.size();++i){
    fBookkeeping.push_back(new GFBookkeeping(*(_tr.fBookkeeping.at(i))));
  }
  fRepAtHit = _tr.fRepAtHit;
}

genf::GFTrack::GFTrack

(

GFAbsTrackRep *

defaultRep

)

Initializing constructor.

A track needs at least one track representation to be useable. The given track representation has to contain starting values for the fit!

Definition at line 25 of file GFTrack.cxx.

  : fTrackReps(NULL),fPDG(2112), fCardinal_rep(0), fNextHitToFit(0)
{
  addTrackRep(defaultRep);
}

genf::GFTrack::~GFTrack ( )

virtual

Definition at line 37 of file GFTrack.cxx.

                      {
  if(fTrackReps!=NULL){
    for(unsigned int i=0;i<getNumReps();i++) {
      delete fTrackReps->At(i);
    }
    delete fTrackReps;
  }
  for(unsigned int i=0;i<fHits.size();i++) {
    delete fHits[i];
  }
  for(unsigned int i=0;i<fBookkeeping.size();++i){
    if(fBookkeeping.at(i)!=NULL) delete fBookkeeping.at(i);
  }
}

Member Function Documentation

void genf::GFTrack::addFailedHit ( unsigned int  irep, unsigned int  id  )

inline

Definition at line 294 of file GFTrack.h.

                                                      {
    fBookkeeping.at(irep)->addFailedHit(id);
  }

void genf::GFTrack::addHit ( GFAbsRecoHit *  theHit )

inline

deprecated!

Definition at line 300 of file GFTrack.h.

                                           {
    fHits.push_back(theHit);
  }

void genf::GFTrack::addHit ( GFAbsRecoHit *  theHit, unsigned int  detId, unsigned int  hitId, double  rho = 0., unsigned int  planeId = 0  )

inline

Add single hit. Updates the GFTrackCand.

Definition at line 306 of file GFTrack.h.

                                     {
    fHits.push_back(theHit);
    fCand.addHit(detId,hitId,rho,planeId);
  }

void genf::GFTrack::addHitVector ( std::vector< GFAbsRecoHit * >  hits )

inline

Add collection of hits.

This is the standard way to fill the track with hit data

Definition at line 319 of file GFTrack.h.

                                                   {
    fHits = hits;
  }

void genf::GFTrack::addTrackRep ( GFAbsTrackRep *  theTrackRep )

inline

Add track represenation.

The given track representation has to contain starting values for fit!

Definition at line 327 of file GFTrack.h.

                                               {
    if(fTrackReps==NULL)fTrackReps=new TObjArray(fDefNumTrackReps);
    fTrackReps->Add(theTrackRep);
    fBookkeeping.push_back( new GFBookkeeping() );
    fRepAtHit.push_back(-1);
  }

void genf::GFTrack::clearBookkeeping ( )

inline

Definition at line 417 of file GFTrack.h.

                         {
    for(unsigned int i=0;i<getNumReps();++i){
      fBookkeeping.at(i)->clearAll();
    }
  }

void genf::GFTrack::clearFailedHits ( )

inline

Definition at line 423 of file GFTrack.h.

                        {
    for(unsigned int i=0;i<getNumReps();++i){
      fBookkeeping.at(i)->clearFailedHits();
    }
  }

void genf::GFTrack::clearRepAtHit ( )

inline

clear the hit indices at which plane,state&cov of reps are defined

Definition at line 405 of file GFTrack.h.

                      {
    for(unsigned int i=0;i<getNumReps();++i){
      fRepAtHit.at(i)=-1;
    }
  }

void genf::GFTrack::fillGeoTrack ( TVirtualGeoTrack *  tr ) const

inline

Fill TVirtualGeoTrack object Cardinal representation is used.

Definition at line 284 of file GFTrack.h.

{fillGeoTrack(tr,fCardinal_rep);}

void genf::GFTrack::fillGeoTrack	(	TVirtualGeoTrack *	tr,
		unsigned int	repid
	)		const

Fill TVirtualGeoTrack object with data from specified track rep.

Definition at line 152 of file GFTrack.cxx.

{
  GFAbsTrackRep* rep=getTrackRep(repid);
  unsigned int n=fCand.getNHits();
  PrintROOTobject(std::cout, rep->getState());
  for(unsigned int i=0; i<n; ++i){// loop over hits
    GFDetPlane pl=fHits[i]->getDetPlane(rep);
    TVector3 pos=rep->getPos(pl);
    #ifdef NDEBUG
      std::cout<<pos.X()<<","<<pos.Y()<<","<<pos.Z()<<std::endl;
    #endif // NDEBUG
    geotrk->AddPoint(pos.X(),pos.Y(),pos.Z(),0);
  }// end loop over hits
}

GFBookkeeping* genf::GFTrack::getBK ( int  index = -1 )

inline

get GFBookKeeping object for particular track rep (default is cardinal rep)

Definition at line 335 of file GFTrack.h.

                                    {
    if(index==-1) return fBookkeeping.at(fCardinal_rep);
    if ((unsigned int)index >= getNumReps())
      throw GFException("genf::GFTrack::getBK(): index out of range", __LINE__, __FILE__).setFatal();
    return fBookkeeping.at(index);
  }

const GFTrackCand& genf::GFTrack::getCand ( ) const

inline

Definition at line 158 of file GFTrack.h.

{return fCand;}

GFAbsTrackRep* genf::GFTrack::getCardinalRep ( ) const

inline

Get cardinal track representation.

The user has to choose which track rep should be considered the best one after the fit. Usually the track representation giving the smallest chi2 is choosen. By default the first in the list is returned.

Definition at line 212 of file GFTrack.h.

{return ((GFAbsTrackRep*)fTrackReps->At(fCardinal_rep));}

double genf::GFTrack::getCharge ( ) const

inline

Get charge from fit.

Cardinal representation is used.

Definition at line 280 of file GFTrack.h.

{return getCardinalRep()->getCharge();}

double genf::GFTrack::getChiSqu ( ) const

inline

Get chi2.

Cardinal representation is used.

Definition at line 262 of file GFTrack.h.

{return getCardinalRep()->getChiSqu();}

int genf::GFTrack::getFailedHits ( int  repId = -1 )

inline

return the number of failed Hits in track fit repId == -1 will use cardinal rep

Definition at line 149 of file GFTrack.h.

                                 {
        int theRep;
        if(repId==-1) theRep=fCardinal_rep;
        else theRep = repId;
        return fBookkeeping.at(theRep)->getNumFailed();
  }

GFAbsRecoHit* genf::GFTrack::getHit ( int  id ) const

inline

Definition at line 160 of file GFTrack.h.

                                     {
    return fHits.at(id);
  }

std::vector< Double_t > genf::GFTrack::getHitChi2 ( )

inline

Definition at line 364 of file GFTrack.h.

{return fHitChi2;}

std::vector< TMatrixT<Double_t> > genf::GFTrack::getHitCov ( )

inline

Definition at line 366 of file GFTrack.h.

{return fHitCov;}

std::vector< TMatrixT<Double_t> > genf::GFTrack::getHitCov7x7 ( )

inline

Definition at line 367 of file GFTrack.h.

{return fHitCov;}

std::vector< TMatrixT<Double_t> > genf::GFTrack::getHitMeasuredCov ( )

inline

Definition at line 362 of file GFTrack.h.

{return fHitMeasCov;}

std::vector< TVector3 > genf::GFTrack::getHitPlaneU ( )

inline

Definition at line 370 of file GFTrack.h.

{return fHitPlaneU;}

std::vector< TVector3 > genf::GFTrack::getHitPlaneUxUyUz ( )

inline

Definition at line 369 of file GFTrack.h.

{return fHitPlaneUxUyUz;}

std::vector< TVector3 > genf::GFTrack::getHitPlaneV ( )

inline

Definition at line 371 of file GFTrack.h.

{return fHitPlaneV;}

std::vector< TVector3 > genf::GFTrack::getHitPlaneXYZ ( )

inline

Definition at line 368 of file GFTrack.h.

{return fHitPlaneXYZ;}

std::vector<GFAbsRecoHit*> genf::GFTrack::getHits ( )

inline

Definition at line 156 of file GFTrack.h.

{return fHits;}

void genf::GFTrack::getHitsByPlane

(

std::vector< std::vector< int > * > &

retVal

)

use planeId information of GFTrackCand and return by ref groups of hit ids which are in the same planes.

Definition at line 220 of file GFTrack.cxx.

                                                                  {
  for(int i=0;retVal.size();++i){
    delete retVal.at(i);
  }
  retVal.clear();
  //this method can only be called when all hits have been loaded
  if (fHits.size() != fCand.getNHits())
    throw GFException("genf::GFTrack::getResiduals(): inconsistent hits", __LINE__, __FILE__).setFatal();
  if (fHits.size() < 2)
    throw GFException("genf::GFTrack::getResiduals(): less than 2 hits", __LINE__, __FILE__).setFatal();
  unsigned int detId,hitId,planeId;
  fCand.getHitWithPlane(0,detId,hitId,planeId);
  //  std::cout << "$$$ " << 0 << " " << detId << " " << hitId << " " << planeId << std::endl;
  unsigned int lastPlane=planeId;
  retVal.push_back(new std::vector<int>);
  retVal.at(0)->push_back(0);
  for(unsigned int i=1;i<fCand.getNHits();++i){
    fCand.getHitWithPlane(i,detId,hitId,planeId);
    //std::cout << "$$$ " << i << " " << detId << " " << hitId << " " << planeId << std::endl;
    if(lastPlane==planeId){
      retVal.at(retVal.size()-1)->push_back(i);
    }
    else{
      lastPlane=planeId;
      retVal.push_back(new std::vector<int>);
      retVal.at(retVal.size()-1)->push_back(i);
    }
  }
}

std::vector< TMatrixT<Double_t> > genf::GFTrack::getHitState ( )

inline

Definition at line 365 of file GFTrack.h.

{return fHitState;}

std::vector< TMatrixT<Double_t> > genf::GFTrack::getHitUpdate ( )

inline

Definition at line 363 of file GFTrack.h.

{return fHitUpdate;}

TVector3 genf::GFTrack::getMom ( ) const

inline

Get momentum at the present position.

Cardinal representation is used.

Definition at line 219 of file GFTrack.h.

{return getCardinalRep()->getMom();}

TVector3 genf::GFTrack::getMom ( const GFDetPlane &  pl ) const

inline

Get momentum at GFDetPlane.

The track will be extrapolated to GFDetPlane to get the momentum there. The track will not be modified. Cardinal representation is used.

Definition at line 226 of file GFTrack.h.

{return getCardinalRep()->getMom(pl);}

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

inline

Get NDF.

Cardinal representation is used.

Definition at line 268 of file GFTrack.h.

{return getCardinalRep()->getNDF();}

unsigned int genf::GFTrack::getNextHitToFit ( ) const

inline

Accessor for fNextHitToFit.

Definition at line 188 of file GFTrack.h.

{return fNextHitToFit;}

unsigned int genf::GFTrack::getNumHits ( ) const

inline

Definition at line 164 of file GFTrack.h.

                                  {
    return fHits.size();
  }

unsigned int genf::GFTrack::getNumReps ( ) const

inline

Get number of track represenatations.

Definition at line 202 of file GFTrack.h.

                                  {
    return fTrackReps->GetEntriesFast();
  }

int genf::GFTrack::getPDG ( )

inline

Definition at line 373 of file GFTrack.h.

{return fPDG;}

TVector3 genf::GFTrack::getPos ( ) const

inline

Get present position.

Cardinal representation is used.

Definition at line 232 of file GFTrack.h.

{return getCardinalRep()->getPos();}

TVector3 genf::GFTrack::getPos ( const GFDetPlane &  pl ) const

inline

Get position at GFDetPlane.

The track will be extrapolated to GFDetPlane to get the position there. The track will not be modified. Cardinal representation is used.

Definition at line 239 of file GFTrack.h.

{return getCardinalRep()->getPos(pl);}

void genf::GFTrack::getPosMomCov ( TVector3 &  pos, TVector3 &  mom, TMatrixT< Double_t > &  cov  )

inline

Get position, momentum, and 6x6 covariance at current position.

Cardinal representation is used.

Definition at line 245 of file GFTrack.h.

                                                                        {
    getCardinalRep()->getPosMomCov(pos,mom,cov);
  }

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

inline

Get position, momentum, and 6x6 covariance at GFDetPlane.

The track will be extrapolated to GFDetPlane to get everything there. The track will not be modified. Cardinal representation is used.

Definition at line 254 of file GFTrack.h.

                                                                                             {
    getCardinalRep()->getPosMomCov(pl,pos,mom,cov);
  }

double genf::GFTrack::getRedChiSqu ( ) const

inline

Get chi2/NDF.

Cardinal representation is used.

Definition at line 274 of file GFTrack.h.

{return getCardinalRep()->getRedChiSqu();}

int genf::GFTrack::getRepAtHit ( unsigned int  irep )

inline

get the hit index at which plane,state&cov of rep irep is defined

Definition at line 397 of file GFTrack.h.

                                    {
    if (irep >= getNumReps())
      throw GFException("genf::GFTrack::getRepAtHit(): index out of range", __LINE__, __FILE__).setFatal();
    return fRepAtHit.at(irep);
  }

void genf::GFTrack::getResiduals	(	unsigned int	detId,
		unsigned int	dim,
		unsigned int	rep,
		std::vector< double > &	result
	)

Get residuals.

Parameters

detId	which detector?
dim	= index of coordinate to choose from resiudal vector
rep	which track representation?
result	results are written to this vector

Definition at line 169 of file GFTrack.cxx.

{
  unsigned int nhits=getNumHits();
  if(repid>=getNumReps())return;
  GFAbsTrackRep* rep=getTrackRep(repid);//->clone();
  for(unsigned int ih=0; ih<nhits; ++ih){// loop over hits
    unsigned int anid;
    unsigned int dummy;
    fCand.getHit(ih,anid,dummy); // check if this is a hit we want to look at
    if(anid==detId){
      GFAbsRecoHit* hit=getHit(ih);
      // extrapolate trackrep there
      int repDim=rep->getDim();
      TMatrixT<Double_t> state(repDim,1);
      GFDetPlane pl=hit->getDetPlane(rep);

      rep->extrapolate(pl,state);
      //rep->setState(state);
      //rep->setReferencePlane(pl);
      double res=hit->residualVector(rep,state,pl)[dim][0];

      //std::cout<<res<<std::endl;

      result.push_back(res);
    }
  }
}

GFAbsTrackRep* genf::GFTrack::getTrackRep ( int  id ) const

inline

Accessor for track representations.

Definition at line 196 of file GFTrack.h.

                                           {
    return reinterpret_cast<GFAbsTrackRep*>(fTrackReps->At(id));
  }

void genf::GFTrack::mergeHits

(

GFTrack *

trk

)

Merge two GFTracks. Only hits will be merged.

All hits from trk will be merged into this GFTrack. trk will be empty afterwards.

Kalman::continueTrack can be used to include the newly added hits in the fit.

Note that the new hits are inserted at the end of the present track!

Definition at line 126 of file GFTrack.cxx.

                                  {
  unsigned int nhits=trk->getNumHits();
  for(unsigned int i=0;i<nhits;++i){
    unsigned int detId;
    unsigned int hitId;
    trk->getCand().getHit(i,detId,hitId);
    GFAbsRecoHit* hit=trk->getHit(i);
    addHit(hit,detId,hitId);
  }
  trk->fHits.clear();
}

genf::GFTrack & genf::GFTrack::operator=

(

const GFTrack &

_tr

)

assignement operator

Definition at line 72 of file GFTrack.cxx.

                                                      {
  if(fTrackReps!=NULL){
    for(unsigned int i=0;i<getNumReps();i++) {
      delete fTrackReps->At(i);
    }
    delete fTrackReps;
    fTrackReps=NULL;
  }
  for(unsigned int i=0;i<fHits.size();i++) {
    delete fHits[i];
  }
  for(unsigned int i=0;i<fBookkeeping.size();++i){
    if(fBookkeeping.at(i)!=NULL) delete fBookkeeping.at(i);
  }

  for(unsigned int i=0;i<_tr.getNumReps();++i){
    addTrackRep(_tr.getTrackRep(i)->clone());
  }
  fCand=_tr.fCand;
  fCardinal_rep=_tr.fCardinal_rep;
  fNextHitToFit=_tr.fNextHitToFit;
  for(unsigned int i=0;i<_tr.getNumHits();i++) {
    fHits.push_back((_tr.getHit(i))->clone());
  }

  //clear the empty bookeeping objs made by addTrackRep and copy the others
  for(unsigned int i=0; i<fBookkeeping.size(); ++i) delete fBookkeeping[i];
  fBookkeeping.clear();
  for(unsigned int i=0;i<_tr.fBookkeeping.size();++i){
    fBookkeeping.push_back(new GFBookkeeping(*(_tr.fBookkeeping.at(i))));
  }
  fRepAtHit = _tr.fRepAtHit;


  return *this;
}

void genf::GFTrack::Print

(

std::ostream &

out = std::cout

)

const

Definition at line 210 of file GFTrack.cxx.

                                               {
  for(unsigned int i=0;i<getNumReps();++i){
    getTrackRep(i)->Print(out);
    fBookkeeping.at(i)->Print(out);
  }
  out << "GFTrack has " << getNumHits() << " detector hits." << std::endl;

}

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

inlineprivatevirtual

Definition at line 435 of file GFTrack.h.

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

void genf::GFTrack::printBookkeeping

(

std::ostream &

out = std::cout

)

const

print bookkeeping

Definition at line 201 of file GFTrack.cxx.

                                                          {
  out << "GFTrack::printBookkeeping()" << std::endl;
  for(unsigned int i=0;i<getNumReps();++i){
    out << "trackRep " << i << ":" << std::endl;
    fBookkeeping.at(i)->Print(out);
  }

}

void genf::GFTrack::releaseHits ( )

inline

Clear hit vector. Note that hits will not be deleted!

Be carefull not to create memory leaks here.

Definition at line 184 of file GFTrack.h.

{fHits.clear();}

void genf::GFTrack::reset

(

)

Resets the GFTrack – deletes RecoHits!

Definition at line 111 of file GFTrack.cxx.

                  {
  if(fTrackReps!=NULL){
    for(unsigned int i=0;i<getNumReps();i++) {
      delete fTrackReps->At(i);
      delete fBookkeeping.at(i);
    }
  }
  for(unsigned int i=0;i<fHits.size();i++) {
    delete fHits[i];
  }
  fHits.clear();
  fRepAtHit.clear();
}

void genf::GFTrack::setCandidate	(	const GFTrackCand &	cand,
		bool	doreset = false
	)

set track candidate

Definition at line 140 of file GFTrack.cxx.

{
  fCand=cand;
  // reset fits
  if(doreset) {
    for(unsigned int i=0;i<getNumReps();i++) {
      ((GFAbsTrackRep*)fTrackReps->At(i))->reset();
    }
  }
}

void genf::GFTrack::setCardinalRep ( unsigned int  r )

inline

Choose cardinal track represenatation.

See also

getCardinalRep

Definition at line 349 of file GFTrack.h.

{if((int)r<fTrackReps->GetEntriesFast())fCardinal_rep=r;}

void genf::GFTrack::setHitChi2 ( Double_t  mat )

inline

Definition at line 353 of file GFTrack.h.

{fHitChi2.push_back(mat);}

void genf::GFTrack::setHitCov ( TMatrixT< Double_t >  mat )

inline

Definition at line 355 of file GFTrack.h.

{fHitCov.push_back(mat);}

void genf::GFTrack::setHitCov7x7 ( TMatrixT< Double_t >  mat )

inline

Definition at line 356 of file GFTrack.h.

{fHitCov7x7.push_back(mat);}

void genf::GFTrack::setHitMeasuredCov ( TMatrixT< Double_t >  mat )

inline

Definition at line 351 of file GFTrack.h.

{fHitMeasCov.push_back(mat);}

void genf::GFTrack::setHitPlaneU ( TVector3  pl )

inline

Definition at line 359 of file GFTrack.h.

{ fHitPlaneU.push_back(pl);}

void genf::GFTrack::setHitPlaneUxUyUz ( TVector3  pl )

inline

Definition at line 358 of file GFTrack.h.

{ fHitPlaneUxUyUz.push_back(pl);}

void genf::GFTrack::setHitPlaneV ( TVector3  pl )

inline

Definition at line 360 of file GFTrack.h.

{ fHitPlaneV.push_back(pl);}

void genf::GFTrack::setHitPlaneXYZ ( TVector3  pl )

inline

Definition at line 357 of file GFTrack.h.

{ fHitPlaneXYZ.push_back(pl);}

void genf::GFTrack::setHitState ( TMatrixT< Double_t >  mat )

inline

Definition at line 354 of file GFTrack.h.

{fHitState.push_back(mat);}

void genf::GFTrack::setHitUpdate ( TMatrixT< Double_t >  mat )

inline

Definition at line 352 of file GFTrack.h.

{fHitUpdate.push_back(mat);}

void genf::GFTrack::setNextHitToFit ( unsigned int  i )

inline

Set next hit to be used in a fit.

Definition at line 192 of file GFTrack.h.

{fNextHitToFit=i;}

void genf::GFTrack::setPDG ( int  pdgt )

inline

Definition at line 361 of file GFTrack.h.

{fPDG = pdgt;}

void genf::GFTrack::setRepAtHit ( unsigned int  irep, int  ihit  )

inline

set the hit index at which plane,state&cov of rep irep is defined

Definition at line 389 of file GFTrack.h.

                                              {
    if (irep >= getNumReps())
      throw GFException("genf::GFTrack::setRepAtHit(): index out of range", __LINE__, __FILE__).setFatal();
    fRepAtHit.at(irep) = ihit;
  }

Member Data Documentation

std::vector< GFBookkeeping* > genf::GFTrack::fBookkeeping

private

Collection of Bookeeping objects for failed hits in every trackrep.

Definition at line 87 of file GFTrack.h.

GFTrackCand genf::GFTrack::fCand

private

Helper to store the indices of the hits in the track. See GFTrackCand for details.

Definition at line 97 of file GFTrack.h.

unsigned int genf::GFTrack::fCardinal_rep

private

Definition at line 114 of file GFTrack.h.

const int genf::GFTrack::fDefNumTrackReps = 10

staticprivate

Definition at line 113 of file GFTrack.h.

std::vector< Double_t > genf::GFTrack::fHitChi2

private

Definition at line 103 of file GFTrack.h.

std::vector< TMatrixT<Double_t> > genf::GFTrack::fHitCov

private

Definition at line 105 of file GFTrack.h.

std::vector< TMatrixT<Double_t> > genf::GFTrack::fHitCov7x7

private

Definition at line 104 of file GFTrack.h.

std::vector< TMatrixT<Double_t> > genf::GFTrack::fHitMeasCov

private

Definition at line 100 of file GFTrack.h.

std::vector< TVector3 > genf::GFTrack::fHitPlaneU

private

Definition at line 108 of file GFTrack.h.

std::vector< TVector3 > genf::GFTrack::fHitPlaneUxUyUz

private

Definition at line 107 of file GFTrack.h.

std::vector< TVector3 > genf::GFTrack::fHitPlaneV

private

Definition at line 109 of file GFTrack.h.

std::vector< TVector3 > genf::GFTrack::fHitPlaneXYZ

private

Definition at line 106 of file GFTrack.h.

std::vector<GFAbsRecoHit*> genf::GFTrack::fHits

private

Collection of RecoHits.

Definition at line 81 of file GFTrack.h.

std::vector< TMatrixT<Double_t> > genf::GFTrack::fHitState

private

Definition at line 102 of file GFTrack.h.

std::vector< TMatrixT<Double_t> > genf::GFTrack::fHitUpdate

private

Definition at line 101 of file GFTrack.h.

unsigned int genf::GFTrack::fNextHitToFit

private

Definition at line 116 of file GFTrack.h.

int genf::GFTrack::fPDG

private

Definition at line 111 of file GFTrack.h.

std::vector<int> genf::GFTrack::fRepAtHit

private

repAtHit keeps track of at which hit index which rep is currently defined, to avoid null extrapolations

Definition at line 92 of file GFTrack.h.

TObjArray* genf::GFTrack::fTrackReps

private

Collection of track representations.

this array is only to be added to in the addTrackRep method because the synchronized construction of bookkeeping objects and repAtHit array is ensured there. NEVER delete elements from this array! If this functionality will be need, it has to be done synchronized with bookkeeping!!

Definition at line 77 of file GFTrack.h.

---

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

• larreco/larreco/Genfit/GFTrack.h
• larreco/larreco/Genfit/GFTrack.cxx