#include <TrkRep.h>

Inheritance diagram for TrkRep:

Public Member Functions
	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 bool	validFlightLength (double fltL, double tolerance=0.0) const

void	printType (std::ostream &ostr) const

virtual TrkExchangePar	helix (double fltL) const =0

Public Member Functions inherited from TrkAbsFit
virtual int	charge () const =0

virtual double	chisq () const =0

virtual int	nDof () const =0

virtual const TrkDifTraj &	traj () const =0

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

virtual double	pt (double fltL=0.) 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 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
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

Protected Attributes
std::auto_ptr< TrkHotList >	_hotList

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

Detailed Description

Definition at line 43 of file TrkRep.h.

Constructor & Destructor Documentation

◆ TrkRep() [1/7]

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

Definition at line 47 of file TrkRep.cxx.

  :_hotList( hotlist.clone(TrkBase::Functors::cloneHot(this)) )
{
  init(trk, hypo);
}

◆ TrkRep() [2/7]

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

Definition at line 54 of file TrkRep.cxx.

{
  init(trk, hypo);
  if (!stealHots) {
    _hotList.reset( hotlist.clone(TrkBase::Functors::cloneHot(this)) );
  } else {
    _hotList.reset( new TrkHotListFull(hotlist,setParent(this)) );
  }
}

◆ TrkRep() [3/7]

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

Definition at line 65 of file TrkRep.cxx.

{
    //yzhang SegGrouperAx::storePar newTrack come here
    //yzhang SegGrouperSt::storePar addZValue come here too and hotlist!=0 do clone()
  init(trk, hypo);
  _hotList.reset( hotlist!=0?
                  hotlist->clone(TrkBase::Functors::cloneHot(this)):
                  new TrkHotListFull );
}

◆ TrkRep() [4/7]

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

Definition at line 76 of file TrkRep.cxx.

{
  init(trk, hypo);
  if (!takeownership) {
    _hotList.reset( hotlist!=0?
                    hotlist->clone(TrkBase::Functors::cloneHot(this)):
                    new TrkHotListFull );
  } else {
    assert(hotlist!=0);
    _hotList.reset( hotlist->resetParent(setParent(this)) );
  }
}

◆ TrkRep() [5/7]

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

Definition at line 91 of file TrkRep.cxx.

  :_hotList(new TrkHotListEmpty(nact, nsv, ndc, stFndRng, endFndRng))
{
//        cout << " in TrkRep copy ctor 1"   << endl;//yzhang debug
    
  init(trk, hypo);
}

◆ TrkRep() [6/7]

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

Definition at line 41 of file TrkRep.cxx.

  : _hotList( createHotList?new TrkHotListFull:0 )
{
  init(trk, hypo);
}

◆ TrkRep() [7/7]

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

Definition at line 101 of file TrkRep.cxx.

                                                                        :
  TrkFitStatus(oldRep)
{
//  cout << " in TrkRep copy ctor 2"   << endl;//yzhang debug
    
  init(trk, hypo);
  // Hots and hotlist have to be cloned in the derived classes
}

◆ ~TrkRep()

TrkRep::~TrkRep ( )

virtual

Definition at line 129 of file TrkRep.cxx.

130{

131}

Member Function Documentation

◆ activateHot()

void TrkRep::activateHot ( TrkHitOnTrk * theHot )

virtual

Definition at line 154 of file TrkRep.cxx.

{
  if(!hot->isActive()){
// make sure this is my hot we're talking about
    if(this == hot->getParentRep()){
      setCurrent(false);
// actually activate the hot; this is now the rep's job
      hot->setActive(true);
    }
  }
}

Referenced by TrkHitOnTrk::setActivity(), and TrkHitOnTrk::setUsability().

◆ addHot()

void TrkRep::addHot ( TrkHitOnTrk * theHot )

virtual

Definition at line 140 of file TrkRep.cxx.

{
  if (newHot->isActive()) setCurrent(false);
  hotList()->append(newHot);
}

◆ arrivalTime()

double TrkRep::arrivalTime ( double fltL ) const

virtual

Implements TrkFit.

Reimplemented in TrkLineRep.

Definition at line 192 of file TrkRep.cxx.

{
  static double cinv = 1./Constants::c;
  // Initialize cache
  if (_betainv < 0.0) {
    double mass2 = Pdt::lookup(particleType())->mass();
    mass2 = mass2 * mass2;
    double ptot2 = momentum(0.).mag2();
    assert(ptot2 != 0.0);
    _betainv = sqrt( (ptot2 +  mass2)/ ptot2);
  }
  double tof = fltL * _betainv * cinv;
  return trackT0() + tof;
}

Referenced by MdcTrkRecon::fillTrackList(), MdcTrack::storeTrack(), HoughFinder::storeTracks(), and CgemMdcFitAlg::updateTracks().

◆ chisqConsistency()

ChisqConsistency TrkRep::chisqConsistency ( ) const

virtual

Implements TrkFit.

Definition at line 363 of file TrkRep.cxx.

                               {
  if(fitValid())
    return ChisqConsistency(chisq(),nDof());
  else
    return ChisqConsistency();
}

◆ clone()

virtual TrkRep * TrkRep::clone ( TrkRecoTrk * newTrack ) const

pure virtual

Implemented in TrkCircleRep, TrkHelixRep, and TrkLineRep.

◆ cloneNewHypo()

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

pure virtual

Implemented in TrkCircleRep, TrkHelixRep, and TrkLineRep.

◆ deactivateHot()

void TrkRep::deactivateHot ( TrkHitOnTrk * theHot )

virtual

Definition at line 167 of file TrkRep.cxx.

{
  if(hot->isActive()){
// make sure this is my hot we're talking about
    if(this == hot->getParentRep()){
      setCurrent(false);
// actually deactivate the hot; this is now the rep's job
      hot->setActive(false);
    }
  }
}

Referenced by TrkHitOnTrk::setActivity(), and TrkHitOnTrk::setUsability().

◆ direction()

Hep3Vector TrkRep::direction ( double fltL ) const

virtual

Implements TrkAbsFit.

Definition at line 186 of file TrkRep.cxx.

{
  return traj().direction(fltL);
}

◆ directionErr()

BesVectorErr TrkRep::directionErr ( double fltL ) const

virtual

Implements TrkAbsFit.

Definition at line 272 of file TrkRep.cxx.

{
  static DifPoint posD;
  static DifVector dirD;
  traj().getDFInfo2(fltL, posD, dirD);
  BesError symErr(3);
  symErr.assign( dirD.errorMatrix( dirD.x.indepPar()->covariance() ));
  Hep3Vector dir(dirD.x.number(), dirD.y.number(), dirD.z.number());
  return BesVectorErr(dir, symErr);
}

◆ endFoundRange()

double TrkRep::endFoundRange ( ) const

virtual

Implements TrkFit.

Definition at line 302 of file TrkRep.cxx.

{
  return hotList()->endFoundRange();
}

◆ endValidRange()

double TrkRep::endValidRange ( ) const

virtual

Implements TrkAbsFit.

Definition at line 290 of file TrkRep.cxx.

{
  return traj().hiRange();
}

◆ fit()

virtual TrkErrCode TrkRep::fit ( )

pure virtual

Implemented in TrkSimpleRep.

Referenced by TrkHitList::fit().

◆ hitCapable()

bool TrkRep::hitCapable ( ) const

inline

Definition at line 118 of file TrkRep.h.

118{return hotList()->hitCapable();}

TrkHotList::hitCapable

virtual bool hitCapable() const =0

Referenced by TrkSimpleRep::fit().

◆ hotList() [1/2]

virtual TrkHotList * TrkRep::hotList ( )

inlinevirtual

Definition at line 109 of file TrkRep.h.

109{return _hotList.get();}

Referenced by addHot(), TrkHelixMaker::addZValues(), TrkRecoTrk::changeDefault(), endFoundRange(), TrkHitList::fit(), TrkSimpleRep::fit(), hitCapable(), TrkHitList::hotList(), TrkRecoTrk::hots(), TrkRecoTrk::hots(), nActive(), nMdc(), nSvt(), removeHot(), startFoundRange(), TrkSimpleRep::TrkSimpleRep(), and updateHots().

◆ hotList() [2/2]

virtual const TrkHotList * TrkRep::hotList ( ) const

inlinevirtual

Definition at line 110 of file TrkRep.h.

110{return _hotList.get();}

◆ myKey()

const IfdKey & TrkRep::myKey ( ) const

virtual

Reimplemented in TrkHelixRep.

Definition at line 314 of file TrkRep.cxx.

{
  // This provides a default key (used to provide Rep-specific interfaces
  //   to TrkRecoTrk consumers).
  static IfdIntKey _theKey(0);
  return _theKey;
}

Referenced by TrkExtInterface::attach(), and TrkExtInterface::attach().

◆ nActive()

int TrkRep::nActive ( ) const

virtual

Implements TrkFit.

Definition at line 330 of file TrkRep.cxx.

{
  return hotList()->nActive();
}

Referenced by TrkSimpleRep::nDof().

◆ nMdc()

int TrkRep::nMdc ( ) const

virtual

Implements TrkFit.

Definition at line 342 of file TrkRep.cxx.

{
  return hotList()->nMdc();
}

◆ nSvt()

int TrkRep::nSvt ( ) const

virtual

Implements TrkFit.

Definition at line 336 of file TrkRep.cxx.

{
  return hotList()->nSvt();
}

◆ operator=()

TrkRep & TrkRep::operator= ( const TrkRep & right )

protected

Definition at line 111 of file TrkRep.cxx.

{
  if(&right != this){
    init(right._parentTrack, right._partHypo);
    _hotList.reset( right._hotList->clone(this) );
    TrkFitStatus::operator=(right);
  }
  return *this;
}

Referenced by TrkSimpleRep::operator=().

◆ operator==()

bool TrkRep::operator== ( const TrkRep & rhs )

Definition at line 134 of file TrkRep.cxx.

{
  return (&rhs == this);
}

◆ parentTrack() [1/2]

TrkRecoTrk * TrkRep::parentTrack ( )

inline

Definition at line 82 of file TrkRep.h.

82{return _parentTrack;}

Referenced by TrkSimpleRep::charge(), TrkCircleRep::cloneNewHypo(), TrkHelixRep::cloneNewHypo(), TrkLineRep::cloneNewHypo(), TrkSimpleRep::getAllCovs(), TrkSimpleRep::getAllWeights(), TrkHitOnTrk::getParentTrack(), TrkSimpleRep::momentum(), TrkSimpleRep::momentumErr(), TrkHitOnTrk::parentTrack(), TrkSimpleRep::posmomCov(), TrkCircleRep::print(), TrkHelixRep::print(), TrkSimpleRep::pt(), CgemHitOnTrack::timeAbsolute(), MdcHitOnTrack::timeAbsolute(), and trackT0().

◆ parentTrack() [2/2]

const TrkRecoTrk * TrkRep::parentTrack ( ) const

inline

Definition at line 83 of file TrkRep.h.

83{return _parentTrack;}

◆ particleType()

PdtPid::PidType TrkRep::particleType ( ) const

virtual

Implements TrkFit.

Definition at line 308 of file TrkRep.cxx.

{
  return _partHypo;
}

Referenced by arrivalTime(), TrkCircleRep::clone(), TrkHelixRep::clone(), TrkLineRep::clone(), and TrkHitOnTrk::particleType().

◆ position()

HepPoint3D TrkRep::position ( double fltL ) const

virtual

Implements TrkAbsFit.

Definition at line 180 of file TrkRep.cxx.

{
  return traj().position(fltL);
}

Referenced by MdcHitOnTrack::wireAmbig().

◆ positionErr()

BesPointErr TrkRep::positionErr ( double fltL ) const

virtual

Implements TrkAbsFit.

Definition at line 214 of file TrkRep.cxx.

{
  static DifPoint posD;
  static DifVector dirD;
  traj().getDFInfo2(fltL, posD, dirD);
  HepMatrix err = posD.errorMatrix( posD.x.indepPar()->covariance() );
  HepPoint3D point(posD.x.number(), posD.y.number(), posD.z.number());
  BesError symErr(3);
  symErr.assign(err);
 
  if (false) {
#ifdef MDCPATREC_ROUTINE 
    cout<< "ErrMsg(routine) " << "Pos " 
        << err.num_row() << " " << err.num_col() << endl
        << "output:" << endl
      //    << err(1,1) << endl
      //    << err(2,1) << "  " << err(2,2) << endl
      //    << err(3,1) << "  " << err(3,2) << "  " << err(3,3) << endl
                    << "x deriv: " << endl
                    << posD.x.derivatives() << endl
                    << "y deriv: " << endl
                    << posD.y.derivatives() << endl
                    << endl;
#endif
    //  }
 
    Hep3Vector pointDir(point.x(), point.y());
    double dirMag = pointDir.mag();
    double dist = 5.e-3;
    double delx = dist * point.x() / dirMag;
    double dely = dist * point.y() / dirMag;
    int ierr = 0;
    HepMatrix weight = err.inverse(ierr);
    double chisq =     weight(1,1) * delx * delx +
      2 * weight(2,1) * delx * dely +
      weight(2,2) * dely * dely;
#ifdef MDCPATREC_DEBUG
    cout << point << endl;
    cout << symErr << endl;
    cout << "delta: " << delx << "  " << dely << endl;
    cout << "chisq: " << chisq << endl;
#endif
    double phi0 = helix(fltL).phi0();
    delx = dist * cos(phi0);
    dely = dist * sin(phi0);
    chisq =            weight(1,1) * delx * delx +
                   2 * weight(2,1) * delx * dely +
                       weight(2,2) * dely * dely;
#ifdef MDCPATREC_DEBUG
    cout << "delta: " << delx << "  " << dely << endl;
    cout << "chisq: " << chisq << endl;
    cout << endl << endl;
#endif
  }
  return BesPointErr(point, symErr);
}

◆ removeHot()

void TrkRep::removeHot ( TrkHitOnTrk * theHot )

virtual

Definition at line 147 of file TrkRep.cxx.

{
  if(theHot->isActive()) setCurrent(false);     // fit no longer current
  hotList()->remove(theHot);
}

◆ resid()

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

virtual

Reimplemented in TrkHelixRep.

Definition at line 348 of file TrkRep.cxx.

{
  assert (h != 0);
  if (h->parentRep() != this) return false; // HOT must belong to Rep
  if (!h->hasResidual()) return false;      // Residual must be available
  if (exclude) return false;            // FIXME: Can't do unbiased residuals (yet!)
 
  residual=h->residual();
  residErr=h->hitRms();
  return true;
}

Referenced by TrkHitOnTrk::resid(), and TrkHitOnTrk::resid().

◆ startFoundRange()

double TrkRep::startFoundRange ( ) const

virtual

Implements TrkFit.

Definition at line 296 of file TrkRep.cxx.

{
  return hotList()->startFoundRange();
}

◆ startValidRange()

double TrkRep::startValidRange ( ) const

virtual

Implements TrkAbsFit.

Definition at line 284 of file TrkRep.cxx.

{
  return traj().lowRange();
}

◆ trackT0()

double TrkRep::trackT0 ( ) const

Definition at line 208 of file TrkRep.cxx.

{
  return parentTrack()->trackT0();
}

Referenced by TrkLineRep::arrivalTime(), and arrivalTime().

◆ trkContext()

const TrkContext & TrkRep::trkContext ( ) const

protected

◆ updateHots()

void TrkRep::updateHots ( )

virtual

Definition at line 323 of file TrkRep.cxx.

{
  setCurrent(false);
  hotList()->updateHots();
}

Referenced by TrkRecoTrk::updateReps().

Member Data Documentation

◆ _hotList

std::auto_ptr<TrkHotList> TrkRep::_hotList

protected

Definition at line 131 of file TrkRep.h.

Referenced by hotList(), hotList(), operator=(), TrkRep(), TrkRep(), TrkRep(), and TrkSimpleRep::TrkSimpleRep().

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

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ TrkRep() [1/7]

◆ TrkRep() [2/7]

◆ TrkRep() [3/7]

◆ TrkRep() [4/7]

◆ TrkRep() [5/7]

◆ TrkRep() [6/7]

◆ TrkRep() [7/7]

◆ ~TrkRep()

Member Function Documentation

◆ activateHot()

◆ addHot()

◆ arrivalTime()

◆ chisqConsistency()

◆ clone()

◆ cloneNewHypo()

◆ deactivateHot()

◆ direction()

◆ directionErr()

◆ endFoundRange()

◆ endValidRange()

◆ fit()

◆ hitCapable()

◆ hotList() [1/2]

◆ hotList() [2/2]

◆ myKey()

◆ nActive()

◆ nMdc()

◆ nSvt()

◆ operator=()

◆ operator==()

◆ parentTrack() [1/2]

◆ parentTrack() [2/2]

◆ particleType()

◆ position()

◆ positionErr()

◆ removeHot()

◆ resid()

◆ startFoundRange()

◆ startValidRange()

◆ trackT0()

◆ trkContext()

◆ updateHots()

Member Data Documentation

◆ _hotList