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 63 of file GFTrack.h.

Constructor & Destructor Documentation

genf::GFTrack::GFTrack

(

)

Default constructor – needed for compatibility with ROOT.

Definition at line 30 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 51 of file GFTrack.cxx.

References addTrackRep(), fBookkeeping, fCand, fCardinal_rep, fHits, fNextHitToFit, fRepAtHit, fTrackReps, getHit(), getNumHits(), getNumReps(), and getTrackRep().

                                      : 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 24 of file GFTrack.cxx.

References addTrackRep().

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

genf::GFTrack::~GFTrack ( )

virtual

Definition at line 36 of file GFTrack.cxx.

References fBookkeeping, fHits, fTrackReps, and getNumReps().

                      {
  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 293 of file GFTrack.h.

Referenced by genf::GFKalman::fittingPass(), and genf::GFDaf::processTrack().

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

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

inline

deprecated!

Definition at line 299 of file GFTrack.h.

Referenced by mergeHits(), trkf::Track3DKalman::produce(), and trkf::Track3DKalmanSPS::produce().

                                           { 
    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 305 of file GFTrack.h.

References genf::GFTrackCand::addHit().

                                     {
    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 318 of file GFTrack.h.

References hits().

                                                   {
    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 326 of file GFTrack.h.

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

                                               {
    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 416 of file GFTrack.h.

References getNumReps().

Referenced by genf::GFDaf::processTrack().

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

void genf::GFTrack::clearFailedHits ( )

inline

Definition at line 422 of file GFTrack.h.

References getHitsByPlane(), getNumReps(), and lar::dump::vector().

Referenced by genf::GFDaf::processTrack().

                        {
    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 404 of file GFTrack.h.

References getNumReps(), Print(), and printBookkeeping().

Referenced by genf::GFKalman::processTrack().

                      {
    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 283 of file GFTrack.h.

References fillGeoTrack().

Referenced by fillGeoTrack().

{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 151 of file GFTrack.cxx.

References fCand, fHits, genf::GFTrackCand::getNHits(), genf::GFAbsTrackRep::getPos(), genf::GFAbsTrackRep::getState(), getTrackRep(), n, and genf::PrintROOTobject().

{
  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 334 of file GFTrack.h.

References getNumReps(), setCandidate(), and GFException::setFatal().

Referenced by genf::GFKalman::fittingPass(), and genf::GFDaf::processTrack().

                                    {
    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 157 of file GFTrack.h.

References fCand.

Referenced by mergeHits().

{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 211 of file GFTrack.h.

Referenced by getCharge(), getChiSqu(), getMom(), getNDF(), getPos(), getPosMomCov(), and getRedChiSqu().

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

double genf::GFTrack::getCharge ( ) const

inline

Get charge from fit.

Cardinal representation is used.

Definition at line 279 of file GFTrack.h.

References getCardinalRep(), and genf::GFAbsTrackRep::getCharge().

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

double genf::GFTrack::getChiSqu ( ) const

inline

Get chi2.

Cardinal representation is used.

Definition at line 261 of file GFTrack.h.

References getCardinalRep(), and genf::GFAbsTrackRep::getChiSqu().

{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 148 of file GFTrack.h.

References fCardinal_rep.

Referenced by trkf::Track3DKalman::produce(), and trkf::Track3DKalmanSPS::produce().

                                 {
        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 159 of file GFTrack.h.

Referenced by getResiduals(), GFTrack(), mergeHits(), operator=(), genf::GFKalman::processHit(), and genf::GFDaf::processTrack().

                                     {
    return fHits.at(id);
  }

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

inline

Definition at line 363 of file GFTrack.h.

References fHitChi2.

Referenced by trkf::Track3DKalmanSPS::produce().

{return fHitChi2;}

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

inline

Definition at line 365 of file GFTrack.h.

References fHitCov.

Referenced by trkf::Track3DKalmanSPS::produce().

{return fHitCov;}

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

inline

Definition at line 366 of file GFTrack.h.

References fHitCov.

Referenced by trkf::Track3DKalmanSPS::produce().

{return fHitCov;}

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

inline

Definition at line 361 of file GFTrack.h.

References fHitMeasCov.

Referenced by trkf::Track3DKalmanSPS::produce().

{return fHitMeasCov;}

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

inline

Definition at line 369 of file GFTrack.h.

References fHitPlaneU.

Referenced by trkf::Track3DKalmanSPS::produce().

{return fHitPlaneU;}

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

inline

Definition at line 368 of file GFTrack.h.

References fHitPlaneUxUyUz.

Referenced by trkf::Track3DKalmanSPS::produce().

{return fHitPlaneUxUyUz;}

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

inline

Definition at line 370 of file GFTrack.h.

References fHitPlaneV.

Referenced by trkf::Track3DKalmanSPS::produce().

{return fHitPlaneV;}

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

inline

Definition at line 367 of file GFTrack.h.

References fHitPlaneXYZ.

Referenced by trkf::Track3DKalmanSPS::produce().

{return fHitPlaneXYZ;}

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

inline

Definition at line 155 of file GFTrack.h.

References fHits.

{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 219 of file GFTrack.cxx.

References fCand, fHits, genf::GFTrackCand::getHitWithPlane(), genf::GFTrackCand::getNHits(), and GFException::setFatal().

Referenced by clearFailedHits(), and genf::GFDaf::processTrack().

                                                                  {
  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 364 of file GFTrack.h.

References fHitState.

Referenced by trkf::Track3DKalmanSPS::produce().

{return fHitState;}

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

inline

Definition at line 362 of file GFTrack.h.

References fHitUpdate.

Referenced by trkf::Track3DKalmanSPS::produce().

{return fHitUpdate;}

TVector3 genf::GFTrack::getMom ( ) const

inline

Get momentum at the present position.

Cardinal representation is used.

Definition at line 218 of file GFTrack.h.

References getCardinalRep(), and genf::GFAbsTrackRep::getMom().

{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 225 of file GFTrack.h.

References getCardinalRep(), and genf::GFAbsTrackRep::getMom().

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

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

inline

Get NDF.

Cardinal representation is used.

Definition at line 267 of file GFTrack.h.

References getCardinalRep(), and genf::GFAbsTrackRep::getNDF().

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

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

inline

Accessor for fNextHitToFit.

Definition at line 187 of file GFTrack.h.

References fNextHitToFit.

Referenced by genf::GFKalman::fittingPass().

{return fNextHitToFit;}

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

inline

Definition at line 163 of file GFTrack.h.

References mergeHits().

Referenced by genf::GFKalman::fittingPass(), getResiduals(), GFTrack(), mergeHits(), operator=(), Print(), genf::GFKalman::processHit(), genf::GFDaf::processTrack(), and genf::GFKalman::processTrack().

                                  {
    return fHits.size();
  }

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

inline

Get number of track represenatations.

Definition at line 201 of file GFTrack.h.

Referenced by genf::GFDaf::blowUpCovs(), genf::GFKalman::blowUpCovs(), genf::GFKalman::blowUpCovsDiag(), clearBookkeeping(), clearFailedHits(), clearRepAtHit(), genf::GFKalman::fittingPass(), getBK(), getRepAtHit(), getResiduals(), GFTrack(), operator=(), Print(), printBookkeeping(), genf::GFDaf::processTrack(), genf::GFKalman::processTrack(), reset(), setCandidate(), setRepAtHit(), genf::GFKalman::switchDirection(), and ~GFTrack().

                                  {
    return fTrackReps->GetEntriesFast();
  }

int genf::GFTrack::getPDG ( )

inline

Definition at line 372 of file GFTrack.h.

References fPDG, and getResiduals().

Referenced by genf::GFKalman::processHit().

{return fPDG;}

TVector3 genf::GFTrack::getPos ( ) const

inline

Get present position.

Cardinal representation is used.

Definition at line 231 of file GFTrack.h.

References getCardinalRep(), and genf::GFAbsTrackRep::getPos().

{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 238 of file GFTrack.h.

References getCardinalRep(), and genf::GFAbsTrackRep::getPos().

{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 244 of file GFTrack.h.

References getCardinalRep(), and genf::GFAbsTrackRep::getPosMomCov().

                                                                        {
    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 253 of file GFTrack.h.

References getCardinalRep(), and genf::GFAbsTrackRep::getPosMomCov().

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

double genf::GFTrack::getRedChiSqu ( ) const

inline

Get chi2/NDF.

Cardinal representation is used.

Definition at line 273 of file GFTrack.h.

References getCardinalRep(), and genf::GFAbsTrackRep::getRedChiSqu().

{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 396 of file GFTrack.h.

References getNumReps(), and GFException::setFatal().

Referenced by genf::GFKalman::processHit().

                                    {
    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 168 of file GFTrack.cxx.

References genf::GFAbsTrackRep::extrapolate(), fCand, genf::GFAbsRecoHit::getDetPlane(), genf::GFAbsTrackRep::getDim(), genf::GFTrackCand::getHit(), getHit(), getNumHits(), getNumReps(), getTrackRep(), and genf::GFAbsRecoHit::residualVector().

Referenced by getPDG().

{
  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 195 of file GFTrack.h.

Referenced by genf::GFDaf::blowUpCovs(), genf::GFKalman::blowUpCovs(), genf::GFKalman::blowUpCovsDiag(), fillGeoTrack(), genf::GFKalman::fittingPass(), getResiduals(), GFTrack(), operator=(), Print(), genf::GFKalman::processHit(), genf::GFDaf::processTrack(), genf::GFKalman::processTrack(), and genf::GFKalman::switchDirection().

                                           {
    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 125 of file GFTrack.cxx.

References addHit(), fHits, getCand(), genf::GFTrackCand::getHit(), getHit(), and getNumHits().

Referenced by getNumHits().

                                  {
  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 71 of file GFTrack.cxx.

References addTrackRep(), genf::GFAbsTrackRep::clone(), fBookkeeping, fCand, fCardinal_rep, fHits, fNextHitToFit, fRepAtHit, fTrackReps, getHit(), getNumHits(), getNumReps(), and getTrackRep().

                                                      {
  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 209 of file GFTrack.cxx.

References fBookkeeping, getNumHits(), getNumReps(), getTrackRep(), and genf::GFAbsTrackRep::Print().

Referenced by clearRepAtHit().

                                               {
  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 434 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 200 of file GFTrack.cxx.

References fBookkeeping, and getNumReps().

Referenced by clearRepAtHit().

                                                          {
  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 183 of file GFTrack.h.

{fHits.clear();} 

void genf::GFTrack::reset

(

)

Resets the GFTrack – deletes RecoHits!

Definition at line 110 of file GFTrack.cxx.

References fBookkeeping, fHits, fRepAtHit, fTrackReps, and getNumReps().

Referenced by setCandidate().

                  {
  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 139 of file GFTrack.cxx.

References fCand, fTrackReps, getNumReps(), and reset().

Referenced by getBK().

{
  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 348 of file GFTrack.h.

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

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

inline

Definition at line 352 of file GFTrack.h.

Referenced by genf::GFKalman::processHit().

{fHitChi2.push_back(mat);}

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

inline

Definition at line 354 of file GFTrack.h.

Referenced by genf::GFKalman::processHit().

{fHitCov.push_back(mat);}

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

inline

Definition at line 355 of file GFTrack.h.

Referenced by genf::GFKalman::processHit().

{fHitCov7x7.push_back(mat);}

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

inline

Definition at line 350 of file GFTrack.h.

Referenced by genf::GFKalman::processHit().

{fHitMeasCov.push_back(mat);}

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

inline

Definition at line 358 of file GFTrack.h.

Referenced by genf::GFKalman::processHit().

{ fHitPlaneU.push_back(pl);}

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

inline

Definition at line 357 of file GFTrack.h.

Referenced by genf::GFKalman::processHit().

{ fHitPlaneUxUyUz.push_back(pl);}

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

inline

Definition at line 359 of file GFTrack.h.

Referenced by genf::GFKalman::processHit().

{ fHitPlaneV.push_back(pl);}

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

inline

Definition at line 356 of file GFTrack.h.

Referenced by genf::GFKalman::processHit().

{ fHitPlaneXYZ.push_back(pl);}

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

inline

Definition at line 353 of file GFTrack.h.

Referenced by genf::GFKalman::processHit().

{fHitState.push_back(mat);}

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

inline

Definition at line 351 of file GFTrack.h.

Referenced by genf::GFKalman::processHit().

{fHitUpdate.push_back(mat);}

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

inline

Set next hit to be used in a fit.

Definition at line 191 of file GFTrack.h.

Referenced by genf::GFKalman::fittingPass(), and genf::GFKalman::processTrack().

{fNextHitToFit=i;}

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

inline

Definition at line 360 of file GFTrack.h.

Referenced by trkf::Track3DKalmanSPS::produce().

{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 388 of file GFTrack.h.

References getNumReps(), and GFException::setFatal().

Referenced by genf::GFKalman::processHit().

                                              {
    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 86 of file GFTrack.h.

Referenced by GFTrack(), operator=(), Print(), printBookkeeping(), reset(), and ~GFTrack().

GFTrackCand genf::GFTrack::fCand

private

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

Definition at line 96 of file GFTrack.h.

Referenced by fillGeoTrack(), getCand(), getHitsByPlane(), getResiduals(), GFTrack(), operator=(), and setCandidate().

unsigned int genf::GFTrack::fCardinal_rep

private

Definition at line 113 of file GFTrack.h.

Referenced by getFailedHits(), GFTrack(), and operator=().

const int genf::GFTrack::fDefNumTrackReps = 10

staticprivate

Definition at line 112 of file GFTrack.h.

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

private

Definition at line 102 of file GFTrack.h.

Referenced by getHitChi2().

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

private

Definition at line 104 of file GFTrack.h.

Referenced by getHitCov(), and getHitCov7x7().

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

private

Definition at line 103 of file GFTrack.h.

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

private

Definition at line 99 of file GFTrack.h.

Referenced by getHitMeasuredCov().

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

private

Definition at line 107 of file GFTrack.h.

Referenced by getHitPlaneU().

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

private

Definition at line 106 of file GFTrack.h.

Referenced by getHitPlaneUxUyUz().

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

private

Definition at line 108 of file GFTrack.h.

Referenced by getHitPlaneV().

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

private

Definition at line 105 of file GFTrack.h.

Referenced by getHitPlaneXYZ().

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

private

Collection of RecoHits.

Definition at line 80 of file GFTrack.h.

Referenced by fillGeoTrack(), getHits(), getHitsByPlane(), GFTrack(), mergeHits(), operator=(), reset(), and ~GFTrack().

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

private

Definition at line 101 of file GFTrack.h.

Referenced by getHitState().

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

private

Definition at line 100 of file GFTrack.h.

Referenced by getHitUpdate().

unsigned int genf::GFTrack::fNextHitToFit

private

Definition at line 115 of file GFTrack.h.

Referenced by getNextHitToFit(), GFTrack(), and operator=().

int genf::GFTrack::fPDG

private

Definition at line 110 of file GFTrack.h.

Referenced by getPDG().

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 91 of file GFTrack.h.

Referenced by GFTrack(), operator=(), and reset().

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 76 of file GFTrack.h.

Referenced by GFTrack(), operator=(), reset(), setCandidate(), and ~GFTrack().

---

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

• /cvmfs/larsoft.opensciencegrid.org/products/larreco/v07_10_02/source/larreco/Genfit/GFTrack.h
• /cvmfs/larsoft.opensciencegrid.org/products/larreco/v07_10_02/source/larreco/Genfit/GFTrack.cxx