#include <TrkHelixRep.h>

Inheritance diagram for TrkHelixRep:

Public Member Functions
	TrkHelixRep (const TrkExchangePar &inPar, TrkRecoTrk myTrack, PdtPid::PidType hypo, const TrkHotList =0)

	TrkHelixRep (const TrkExchangePar &inPar, TrkRecoTrk *trk, PdtPid::PidType hypo, int nact, int nsv, int ndc, double chi2, double startFoundRange, double endFoundRange)

virtual	~TrkHelixRep ()

virtual TrkHelixRep *	clone (TrkRecoTrk *newTrack) const

virtual TrkHelixRep *	cloneNewHypo (PdtPid::PidType hypo)

virtual TrkDifTraj &	traj ()

virtual const TrkDifTraj &	traj () const

const IfdKey &	myKey () const

virtual TrkExchangePar	helix (double fltLen) const

virtual bool	resid (const TrkHitOnTrk *h, double &residual, double &residErr, bool exclude) const

virtual void	printAll (std::ostream &ostr) const

virtual void	print (std::ostream &ostr) const

virtual TrkSimpTraj &	simpTraj ()

virtual const TrkSimpTraj &	simpTraj () const

Public Member Functions inherited from TrkSimpleRep
	TrkSimpleRep (const TrkHotList , TrkRecoTrk , PdtPid::PidType)

	TrkSimpleRep (TrkRecoTrk *, PdtPid::PidType, int nActive, int nSvt, int nDch, double chi2, double startFoundRange, double endFoundRange)

	TrkSimpleRep (const TrkSimpleRep &right, TrkRecoTrk *trk, PdtPid::PidType)

virtual	~TrkSimpleRep ()

int	nDof () const

double	chisq () const

void	setChisq (double c)

virtual int	charge () const

virtual Hep3Vector	momentum (double fltL=0.) const

virtual double	pt (double fltL=0.) const

virtual BesVectorErr	momentumErr (double fltL) const

virtual HepMatrix	posmomCov (double fltL) const

virtual void	getAllCovs (double fltL, HepSymMatrix &xxCov, HepSymMatrix &ppCov, HepMatrix &xpCov) const

virtual void	getAllWeights (double fltL, HepVector &pos, HepVector &mom, HepSymMatrix &xxWeight, HepSymMatrix &ppWeight, HepMatrix &xpWeight) const

virtual void	getAllWeights (const HepPoint3D &pt, HepVector &pos, HepVector &mom, HepSymMatrix &xxWeight, HepSymMatrix &ppWeight, HepMatrix &xpWeight) const

virtual TrkErrCode	fit ()

TrkHelixFitter &	fitter ()

Public Member Functions inherited from TrkRep
	TrkRep (const TrkHotList &inHots, TrkRecoTrk *trk, PdtPid::PidType hypo)

	TrkRep (TrkHotList &inHots, TrkRecoTrk *trk, PdtPid::PidType hypo, bool stealHots=false)

	TrkRep (const TrkHotList inHots, TrkRecoTrk trk, PdtPid::PidType hypo)

	TrkRep (TrkHotList inHots, TrkRecoTrk trk, PdtPid::PidType hypo, bool takeownership=false)

	TrkRep (TrkRecoTrk *trk, PdtPid::PidType hypo, int nAct, int nSvt, int nMdc, double startFoundRange, double endFoundRange)

	TrkRep (TrkRecoTrk *trk, PdtPid::PidType hypo, bool createHotList=false)

	TrkRep (const TrkRep &oldRep, TrkRecoTrk *trk, PdtPid::PidType hypo)

virtual	~TrkRep ()

virtual TrkRep *	clone (TrkRecoTrk *newTrack) const =0

virtual TrkRep *	cloneNewHypo (PdtPid::PidType hypo)=0

bool	operator== (const TrkRep &)

virtual ChisqConsistency	chisqConsistency () const

virtual int	nActive () const

virtual int	nSvt () const

virtual int	nMdc () const

virtual PdtPid::PidType	particleType () const

TrkRecoTrk *	parentTrack ()

const TrkRecoTrk *	parentTrack () const

double	startValidRange () const

double	endValidRange () const

virtual double	startFoundRange () const

virtual double	endFoundRange () const

virtual const IfdKey &	myKey () const

double	trackT0 () const

virtual HepPoint3D	position (double fltL) const

virtual Hep3Vector	direction (double fltL) const

virtual double	arrivalTime (double fltL) const

virtual BesPointErr	positionErr (double fltL) const

virtual BesVectorErr	directionErr (double fltL) const

virtual void	addHot (TrkHitOnTrk *theHot)

virtual void	removeHot (TrkHitOnTrk *theHot)

virtual void	activateHot (TrkHitOnTrk *theHot)

virtual void	deactivateHot (TrkHitOnTrk *theHot)

virtual TrkHotList *	hotList ()

virtual const TrkHotList *	hotList () const

virtual void	updateHots ()

virtual bool	resid (const TrkHitOnTrk *theHot, double &residual, double &residErr, bool exclude=false) const

bool	hitCapable () const

virtual TrkErrCode	fit ()=0

Public Member Functions inherited from TrkFitStatus
bool	is2d () const

bool	fitCurrent () const

bool	fitValid () const

bool	multScat () const

std::ostream &	printStatus (std::ostream &os=std::cout) const

void	setValid (bool v)

void	setCurrent (bool c)

void	set2d (bool d)

void	setMultScat (bool m)

const TrkErrCode &	fitStatus () const

history_iterator	beginHistory () const

history_iterator	endHistory () const

history_riterator	reverseBeginHistory () const

history_riterator	reverseEndHistory () const

std::pair< history_iterator, history_iterator >	history () const

const std::vector< TrkHistory > &	historyVector () const

virtual void	addHistory (const TrkErrCode &status, const char *modulename)

template<class T >
void	addHistory (T begin, T end)

template<class T >
void	addHistory (std::pair< T, T > p)

std::ostream &	printHistory (std::ostream &os=std::cout) const

Public Member Functions inherited from TrkFit
virtual ChisqConsistency	chisqConsistency () const =0

virtual bool	validFlightLength (double fltL, double tolerance=0.0) const

void	printType (std::ostream &ostr) const

virtual int	nActive () const =0

virtual int	nSvt () const =0

virtual int	nMdc () const =0

virtual PdtPid::PidType	particleType () const =0

virtual TrkExchangePar	helix (double fltL) const =0

virtual double	arrivalTime (double fltL) const =0

virtual double	startFoundRange () const =0

virtual double	endFoundRange () const =0

virtual int	charge () const =0

virtual double	chisq () const =0

virtual int	nDof () const =0

virtual const TrkDifTraj &	traj () const =0

virtual HepPoint3D	position (double fltL) const =0

virtual Hep3Vector	direction (double fltL) const =0

virtual Hep3Vector	momentum (double fltL=0.) const =0

virtual double	pt (double fltL=0.) const =0

virtual BesPointErr	positionErr (double fltL) const =0

virtual BesVectorErr	directionErr (double fltL) const =0

virtual BesVectorErr	momentumErr (double fltL) const =0

virtual HepMatrix	posmomCov (double fltL) const =0

virtual void	getAllCovs (double fltL, HepSymMatrix &xxCov, HepSymMatrix &ppCov, HepMatrix &xpCov) const =0

virtual void	getAllWeights (double fltL, HepVector &pos, HepVector &mom, HepSymMatrix &xxWeight, HepSymMatrix &ppWeight, HepMatrix &xpWeight) const =0

virtual double	startValidRange () const =0

virtual double	endValidRange () const =0

virtual void	printAll (std::ostream &ostr) const =0

virtual void	print (std::ostream &ostr) const =0

Public Member Functions inherited from TrkHitOnTrkUpdater
virtual	~TrkHitOnTrkUpdater ()=0

Protected Member Functions
	TrkHelixRep (const TrkHelixRep &right, TrkRecoTrk *trk, PdtPid::PidType)

Protected Member Functions inherited from TrkSimpleRep
TrkSimpleRep &	operator= (const TrkSimpleRep &right)

virtual TrkSimpTraj &	simpTraj ()=0

virtual const TrkSimpTraj &	simpTraj () const =0

Protected Member Functions inherited from TrkRep
TrkRep &	operator= (const TrkRep &)

const TrkContext &	trkContext () const

Protected Member Functions inherited from TrkFitStatus
virtual	~TrkFitStatus ()

	TrkFitStatus ()

TrkFitStatus &	operator= (const TrkFitStatus &)

	TrkFitStatus (const TrkFitStatus &)

Protected Member Functions inherited from TrkFit
	TrkFit ()

virtual	~TrkFit ()

Protected Member Functions inherited from TrkAbsFit
	TrkAbsFit ()

virtual	~TrkAbsFit ()

Protected Member Functions inherited from TrkHitOnTrkUpdater
TrkErrCode	updateMeasurement (TrkHitOnTrk &hot, const TrkDifTraj *traj=0, bool maintainAmbiguity=false) const

void	setActivity (TrkHitOnTrk &hot, bool active) const

void	setParent (TrkHitOnTrk &hot, TrkRep *parent) const

TrkBase::Functors::updateMeasurement	updateMeasurement (const TrkDifTraj *traj=0, bool maintainAmbiguity=false) const

TrkBase::Functors::setParent	setParent (TrkRep *parent) const

TrkBase::Functors::setActive	setActive (bool active) const

Additional Inherited Members
Public Types inherited from TrkFitStatus
typedef std::unary_function< TrkFitStatus, bool >	predicate_type

Protected Attributes inherited from TrkRep
std::auto_ptr< TrkHotList >	_hotList

Detailed Description

Definition at line 22 of file TrkHelixRep.h.

Constructor & Destructor Documentation

◆ TrkHelixRep() [1/3]

TrkHelixRep::TrkHelixRep	(	const TrkExchangePar &	inPar,
		TrkRecoTrk *	myTrack,
		PdtPid::PidType	hypo,
		const TrkHotList *	hots = `0`
	)

Definition at line 21 of file TrkHelixRep.cxx.

                            :
  TrkSimpleRep(hots, trk, hypo),
  _traj(inPar)
{}

◆ TrkHelixRep() [2/3]

TrkHelixRep::TrkHelixRep	(	const TrkExchangePar &	inPar,
		TrkRecoTrk *	trk,
		PdtPid::PidType	hypo,
		int	nact,
		int	nsv,
		int	ndc,
		double	chi2,
		double	startFoundRange,
		double	endFoundRange
	)

Definition at line 31 of file TrkHelixRep.cxx.

                                       :
  TrkSimpleRep(trk, hypo, nact, nsv, ndc, chi2, stFndRng, endFndRng),
  _traj(inPar.params(), inPar.covariance())
{}

◆ ~TrkHelixRep()

TrkHelixRep::~TrkHelixRep ( )

virtual

Definition at line 43 of file TrkHelixRep.cxx.

44{}

◆ TrkHelixRep() [3/3]

TrkHelixRep::TrkHelixRep	(	const TrkHelixRep &	right,
		TrkRecoTrk *	trk,
		PdtPid::PidType	hypo
	)

protected

Definition at line 39 of file TrkHelixRep.cxx.

                                                                                    :
  TrkSimpleRep(old, trk, hypo), _traj(old._traj)
{}

Member Function Documentation

◆ clone()

TrkHelixRep * TrkHelixRep::clone ( TrkRecoTrk * newTrack ) const

virtual

Implements TrkRep.

Definition at line 47 of file TrkHelixRep.cxx.

{
  TrkHelixRep* newRep = new TrkHelixRep(*this, theTrack, this->particleType());
  newRep->setCurrent(fitCurrent());
  return newRep;
}

◆ cloneNewHypo()

TrkHelixRep * TrkHelixRep::cloneNewHypo ( PdtPid::PidType hypo )

virtual

Implements TrkRep.

Definition at line 55 of file TrkHelixRep.cxx.

{
  TrkHelixRep* newRep = new TrkHelixRep(*this, parentTrack(), hypo);
  newRep->setValid(fitValid());
  newRep->setCurrent(fitCurrent());   // true for helices, anyway
  return newRep;
}

◆ helix()

TrkExchangePar TrkHelixRep::helix ( double fltLen ) const

virtual

Implements TrkFit.

Definition at line 64 of file TrkHelixRep.cxx.

{
  HepVector outPar = _traj.parameters()->parameter();
  outPar[HelixTraj::phi0Index] = BesAngle(outPar[HelixTraj::phi0Index]).rad();
  return TrkExchangePar(outPar, _traj.parameters()->covariance());
}

◆ myKey()

const IfdKey & TrkHelixRep::myKey ( ) const

virtual

Reimplemented from TrkRep.

Definition at line 118 of file TrkHelixRep.cxx.

{
  static IfdIntKey _theKey(3133);
  return _theKey;
}

◆ print()

void TrkHelixRep::print ( std::ostream & ostr ) const

virtual

Implements TrkAbsFit.

Definition at line 72 of file TrkHelixRep.cxx.

{
  ostr << "TrkHelixRep "
    << " q: "<<charge()
    << " phi0: " << BesAngle(_traj.phi0()).rad()
    << " om: "<<_traj.omega()
    << " pt: "<<pt()
    << " p: "<<momentum()
    << " d0: " << _traj.d0() 
    << " z0: " << _traj.z0()
    << " ct: " << _traj.tanDip() 
    << "parent track: "<<parentTrack()->id() << endl;
}

Referenced by printAll().

◆ printAll()

void TrkHelixRep::printAll ( std::ostream & ostr ) const

virtual

Implements TrkAbsFit.

Definition at line 87 of file TrkHelixRep.cxx.

{
  print(ostr);
}

◆ resid()

bool TrkHelixRep::resid	(	const TrkHitOnTrk *	h,
		double &	residual,
		double &	residErr,
		bool	exclude
	)		const

virtual

Reimplemented from TrkRep.

Definition at line 126 of file TrkHelixRep.cxx.

{
  if (this != h->getParentRep() ) return false;
  static HepVector derivs; double deltaChi;
  TrkErrCode errCode = h->getFitStuff( derivs, deltaChi );
  // TODO: add caching of derivatives here...
  if (errCode.failure()) return false;
  if (!exclude || !h->isActive()) {
    HepSymMatrix c = _traj.parameters()->covariance().similarityT(derivs);
    double tErr2 = c(1,1);
    double rms = h->hitRms();
    residual=h->residual();
    double e2=rms*rms+(h->isActive()?-1:1)*tErr2;
    if (e2<0) return false;
    residErr=sqrt(e2);
  } else {
    assert(h->isActive());
    // FIXME: never really checked this is really OK...
    // FIXME: should compare to HelixRefiner...
    HepSymMatrix w=_traj.parameters()->covariance();
    int ier;
    w.invert(ier);
    if (ier!=0) return false;
    HepVector p0 = _traj.parameters()->parameter();
    HepVector p1 = w*p0;
    p1 -= deltaChi*derivs;
    for (unsigned i=0; i<p0.num_row(); ++i) for (unsigned j=0; j<i+1 ;++j) {
      w.fast(i+1,j+1) -= derivs[i]*derivs[j];
    }
    w.invert(ier);
    if (ier!=0) return false;
    residual=h->residual()+h->hitRms()*dot(derivs,w*p1-p0);
    residErr=h->hitRms()*sqrt(w(1,1)+1);
  }
  return true;
}

◆ simpTraj() [1/2]

TrkSimpTraj & TrkHelixRep::simpTraj ( )

virtual

Implements TrkSimpleRep.

Definition at line 93 of file TrkHelixRep.cxx.

{
  return _traj;
}

◆ simpTraj() [2/2]

const TrkSimpTraj & TrkHelixRep::simpTraj ( ) const

virtual

Implements TrkSimpleRep.

Definition at line 99 of file TrkHelixRep.cxx.

{
  return _traj;
}

◆ traj() [1/2]

TrkDifTraj & TrkHelixRep::traj ( )

virtual

Definition at line 105 of file TrkHelixRep.cxx.

{
  return _traj;
}

◆ traj() [2/2]

const TrkDifTraj & TrkHelixRep::traj ( ) const

virtual

Implements TrkAbsFit.

Definition at line 111 of file TrkHelixRep.cxx.

{
  return _traj;
}

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

Public Member Functions

Protected Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ TrkHelixRep() [1/3]

◆ TrkHelixRep() [2/3]

◆ ~TrkHelixRep()

◆ TrkHelixRep() [3/3]

Member Function Documentation

◆ clone()

◆ cloneNewHypo()

◆ helix()

◆ myKey()

◆ print()

◆ printAll()

◆ resid()

◆ simpTraj() [1/2]

◆ simpTraj() [2/2]

◆ traj() [1/2]

◆ traj() [2/2]