TrkDifLineTraj Class Reference

#include <TrkDifLineTraj.h>

Inheritance diagram for TrkDifLineTraj:

Public Types
enum	{ d0Ind =0 , phi0Ind , z0Ind , tanDipInd }
enum	{ NLINPRM = 4 }

Public Member Functions
	TrkDifLineTraj (const HepVector &, const HepSymMatrix &, double lowlim=-99999., double hilim=99999., const HepPoint3D &refpoint=_theOrigin)
	TrkDifLineTraj (const TrkExchangePar &, double lowlim=-99999., double hilim=99999., const HepPoint3D &refpoint=_theOrigin)
	TrkDifLineTraj (const TrkDifLineTraj &)
TrkDifLineTraj *	clone () const
	~TrkDifLineTraj ()
TrkDifLineTraj &	operator= (const TrkDifLineTraj &)
virtual HepPoint3D	position (double fltLen) const
virtual Hep3Vector	direction (double fltLen) const
virtual Hep3Vector	delDirect (double) const
virtual void	getInfo (double fltLen, HepPoint3D &pos, Hep3Vector &dir) const
virtual void	getInfo (double fltLen, HepPoint3D &, Hep3Vector &dir, Hep3Vector &delDir) const
virtual void	getDFInfo (double fltLen, DifPoint &, DifVector &dir, DifVector &delDir) const
virtual double	curvature (double fltLen) const
double	cosDip () const
double	d0 () const
double	z0 () const
double	tanDip () const
double	phi0 () const
int	nPar () const
virtual double	distTo1stError (double flt, double tol, int pathDir) const
virtual double	distTo2ndError (double flt, double tol, int pathDir) const
HepMatrix	derivDeflect (double fltlen, deflectDirection) const
HepMatrix	derivDisplace (double fltlen, deflectDirection) const
HepMatrix	derivPFract (double fltlen) const
TranslateParams	paramFunction () const
void	invertParams (TrkParams *params, std::vector< bool > &flags) const
virtual void	visitAccept (TrkVisitor *vis) const
Public Member Functions inherited from TrkSimpTraj
	TrkSimpTraj (const HepVector &params, const HepSymMatrix &cov, const double startRange=-99999., const double endRange=99999., const HepPoint3D &refpoint=_theOrigin)
	TrkSimpTraj (const TrkParams &params, const double startRange=-99999., const double endRange=99999., const HepPoint3D &refpoint=_theOrigin)
virtual	~TrkSimpTraj ()
virtual TrkSimpTraj *	clone () const =0
TrkParams *	parameters ()
const TrkParams *	parameters () const
virtual const TrkSimpTraj *	localTrajectory (double fltLen, double &localFlt) const
const HepPoint3D &	referencePoint () const
virtual void	print (std::ostream &os) const
virtual void	printAll (std::ostream &os) const
virtual int	nPar () const
void	changePoint (const HepPoint3D &newpoint, double &fltlen)
void	setPoint (const HepPoint3D &newpoint)
TrkSimpTraj &	invert ()
virtual void	invertParams (TrkParams *newparams, std::vector< bool > &flags) const =0
virtual TranslateParams	paramFunction () const =0
virtual void	visitAccept (TrkVisitor *vis) const =0
bool	operator== (const TrkSimpTraj &) const
Public Member Functions inherited from TrkDifTraj
	TrkDifTraj (const double lowlim=-99999., const double hilim=99999.)
virtual	~TrkDifTraj ()
virtual void	getDFInfo (double fltLen, DifPoint &pos, DifVector &direction, DifVector &delDirect) const =0
virtual void	getDFInfo2 (double fltLen, DifPoint &pos, DifVector &direction) const
virtual const TrkSimpTraj *	localTrajectory (double fltLen, double &localFlt) const =0
Public Member Functions inherited from Trajectory
	Trajectory (double lowlim, double hilim)
virtual	~Trajectory ()
virtual HepPoint3D	position (double) const =0
virtual Hep3Vector	direction (double) const =0
virtual Hep3Vector	delDirect (double) const =0
virtual double	curvature (double) const =0
virtual void	getInfo (double fltLen, HepPoint3D &pos, Hep3Vector &direction) const =0
virtual void	getInfo (double fltLen, HepPoint3D &pos, Hep3Vector &direction, Hep3Vector &delDirect) const =0
double	distTo0thError (double s, double tol, int pathDir) const
virtual double	distTo1stError (double s, double tol, int pathDir) const =0
virtual double	distTo2ndError (double s, double tol, int pathDir) const =0
virtual Trajectory *	clone () const =0
bool	validFlightDistance (double f, double tolerance=0.0) const
virtual void	setFlightRange (double newrange[2])
double	lowRange () const
double	hiRange () const
double	range () const
virtual void	print (std::ostream &os) const
virtual void	printAll (std::ostream &os) const
virtual HepMatrix	derivDeflect (double fltlen, deflectDirection idir) const =0
virtual HepMatrix	derivDisplace (double fltlen, deflectDirection idir) const =0
virtual HepMatrix	derivPFract (double fltlen) const =0

Additional Inherited Members
Static Public Attributes inherited from TrkSimpTraj
static HepPoint3D	_theOrigin
Protected Member Functions inherited from Trajectory
Trajectory &	operator= (const Trajectory &)
Protected Attributes inherited from TrkSimpTraj
TrkParams	_dtparams
HepPoint3D	_refpoint
Protected Attributes inherited from Trajectory
double	flightrange [2]

Detailed Description

Definition at line 25 of file TrkDifLineTraj.h.

Member Enumeration Documentation

anonymous enum

anonymous enum

Enumerator
d0Ind
phi0Ind
z0Ind
tanDipInd

Definition at line 29 of file TrkDifLineTraj.h.

{d0Ind=0, phi0Ind, z0Ind, tanDipInd};

anonymous enum

anonymous enum

Enumerator
NLINPRM

Definition at line 30 of file TrkDifLineTraj.h.

{NLINPRM = 4};

Constructor & Destructor Documentation

TrkDifLineTraj() [1/3]

TrkDifLineTraj::TrkDifLineTraj	(	const HepVector &	pvec,
		const HepSymMatrix &	pcov,
		double	lowlim = -99999.,
		double	hilim = 99999.,
		const HepPoint3D &	refpoint = _theOrigin
	)

Definition at line 25 of file TrkDifLineTraj.cxx.

                                                                                        :
  TrkSimpTraj(pvec, pcov, lowlim,hilim,refpoint)
{
}

TrkDifLineTraj() [2/3]

TrkDifLineTraj::TrkDifLineTraj	(	const TrkExchangePar &	inpar,
		double	lowlim = -99999.,
		double	hilim = 99999.,
		const HepPoint3D &	refpoint = _theOrigin
	)

Definition at line 31 of file TrkDifLineTraj.cxx.

                                                                                        :
  TrkSimpTraj(HepVector(NLINPRM,1),HepSymMatrix(NLINPRM,1), lowlim,hilim,refpoint)
{
// the above constructor actually screws up z0 and tandip.  I can't fix it in
// the initializer, so I'll just overwrite the track here.  Ugly and inefficient
  HepVector subvect(NLINPRM,1);
  HepSymMatrix submat(NLINPRM,1);
  const HepSymMatrix covar = inpar.covariance();
  subvect[d0Ind] = inpar.d0();
  subvect[phi0Ind] = inpar.phi0();
  subvect[z0Ind] = inpar.z0();
  subvect[tanDipInd] = inpar.tanDip();
 
  submat.fast(d0Ind+1,d0Ind+1) = covar.fast(TrkExchangePar::ex_d0+1,TrkExchangePar::ex_d0+1);
  submat.fast(d0Ind+1,phi0Ind+1) = covar.fast(TrkExchangePar::ex_d0+1,TrkExchangePar::ex_phi0+1);
  submat.fast(d0Ind+1,z0Ind+1) = covar.fast(TrkExchangePar::ex_d0+1,TrkExchangePar::ex_z0+1);
  submat.fast(d0Ind+1,tanDipInd+1) = covar.fast(TrkExchangePar::ex_d0+1,TrkExchangePar::ex_tanDip+1);
  submat.fast(phi0Ind+1,phi0Ind+1) = covar.fast(TrkExchangePar::ex_phi0+1,TrkExchangePar::ex_phi0+1);
  submat.fast(phi0Ind+1,z0Ind+1) = covar.fast(TrkExchangePar::ex_phi0+1,TrkExchangePar::ex_z0+1);
  submat.fast(phi0Ind+1,tanDipInd+1) = covar.fast(TrkExchangePar::ex_phi0+1,TrkExchangePar::ex_tanDip+1);
  submat.fast(z0Ind+1,z0Ind+1) = covar.fast(TrkExchangePar::ex_z0+1,TrkExchangePar::ex_z0+1);
  submat.fast(z0Ind+1,tanDipInd+1) = covar.fast(TrkExchangePar::ex_z0+1,TrkExchangePar::ex_tanDip+1);
  submat.fast(tanDipInd+1,tanDipInd+1) = covar.fast(TrkExchangePar::ex_tanDip+1,TrkExchangePar::ex_tanDip+1);
// reset the track
  (*parameters()) = TrkParams(subvect,submat);
}

TrkDifLineTraj() [3/3]

TrkDifLineTraj::TrkDifLineTraj

(

const TrkDifLineTraj &

h

)

Definition at line 59 of file TrkDifLineTraj.cxx.

  : TrkSimpTraj(h.parameters()->parameter(), h.parameters()->covariance(),
                h.lowRange(),h.hiRange(),h.referencePoint())
{
}

~TrkDifLineTraj()

TrkDifLineTraj::~TrkDifLineTraj

(

)

Definition at line 82 of file TrkDifLineTraj.cxx.

{
}

Member Function Documentation

clone()

TrkDifLineTraj * TrkDifLineTraj::clone ( ) const

virtual

Implements TrkSimpTraj.

Definition at line 66 of file TrkDifLineTraj.cxx.

{
  return new TrkDifLineTraj(*this);
}

cosDip()

double TrkDifLineTraj::cosDip ( ) const

inline

Definition at line 54 of file TrkDifLineTraj.h.

{return 1./sqrt(1.+tanDip()*tanDip());}

Referenced by derivDeflect(), derivDisplace(), direction(), and position().

curvature()

double TrkDifLineTraj::curvature ( double  fltLen ) const

virtual

Implements Trajectory.

Definition at line 202 of file TrkDifLineTraj.cxx.

{
  return 0.;
}

d0()

double TrkDifLineTraj::d0 ( ) const

inline

Definition at line 55 of file TrkDifLineTraj.h.

{return parameters()->parameter()[d0Index()];}

Referenced by derivDeflect(), getDFInfo(), TrkLineRep::helix(), position(), and TrkLineRep::print().

delDirect()

Hep3Vector TrkDifLineTraj::delDirect ( double  ) const

virtual

Implements Trajectory.

Definition at line 127 of file TrkDifLineTraj.cxx.

{
  return Hep3Vector(0.0, 0.0, 0.0);
}

Referenced by getInfo().

derivDeflect()

HepMatrix TrkDifLineTraj::derivDeflect ( double  fltlen, deflectDirection  idirect  ) const

virtual

Implements TrkKalDeriv.

Definition at line 214 of file TrkDifLineTraj.cxx.

{
//  This function computes the column matrix of derivatives for the change
//  in parameters for a change in the direction of a track at a point along
//  its flight, holding the momentum and position constant.  The effects for
//  changes in 2 perpendicular directions (theta1 = dip and
//  theta2 = phi*cos(dip)) are uncorrelated.
 
  HepMatrix ddflct(nLinPrm(),1);
//  Compute some common things
  double cosd = cosDip();
//  Go through the parameters
  switch (idirect) {
  case theta1:
    ddflct[tanDipIndex()][0] = 1.0/(cosd*cosd);
    ddflct[d0Index()][0] = 0.0;
    ddflct[phi0Index()][0] =  0.0;
    ddflct[z0Index()][0] = -fltlen/cosd;;
    break;
  case theta2:
    ddflct[tanDipIndex()][0] = 0;
    ddflct[d0Index()][0] = -fltlen;
    ddflct[phi0Index()][0] = 1.0/cosd;
    ddflct[z0Index()][0] =  -(tanDip()/cosd)*d0();
    break;
  }
  return ddflct;
}

derivDisplace()

HepMatrix TrkDifLineTraj::derivDisplace ( double  fltlen, deflectDirection  idirect  ) const

virtual

Implements TrkKalDeriv.

Definition at line 244 of file TrkDifLineTraj.cxx.

{
//  This function computes the column matrix of derivatives for the change
//  in parameters for a change in the position of a track at a point along
//  its flight, holding the momentum and direction constant.  The effects for
//  changes in 2 perpendicular directions (theta1 = dip and
//  theta2 = phi*cos(dip)) are uncorrelated.
 
  HepMatrix ddflct(nLinPrm(),1);
//  Compute some common things
  double cosd = cosDip();
//  Go through the parameters
  switch (idirect) {
  case theta1:
    ddflct[tanDipIndex()][0] = 0.0;
    ddflct[d0Index()][0] = 0.0;
    ddflct[phi0Index()][0] =  0.0;
    ddflct[z0Index()][0] = 1.0/cosd;
    break;
  case theta2:
    ddflct[tanDipIndex()][0] = 0;
    ddflct[d0Index()][0] = 1.0;
    ddflct[phi0Index()][0] = 0.0;
    ddflct[z0Index()][0] =  0.0;
    break;
  }
  return ddflct;
}

derivPFract()

HepMatrix TrkDifLineTraj::derivPFract ( double  fltlen ) const

virtual

Implements TrkKalDeriv.

Definition at line 274 of file TrkDifLineTraj.cxx.

{
//  This function computes the column matrix of derrivatives for the change
//  in parameters from a (fractional) change in the track momentum,
//  holding the direction and position constant.  The momentum change can
//  come from energy loss or bfield inhomogeneities.  Line Trajs have no
//  geometric change from energy loss.
 
  return HepMatrix(nLinPrm(),1, 0);
}

direction()

Hep3Vector TrkDifLineTraj::direction ( double  fltLen ) const

virtual

Implements Trajectory.

Definition at line 118 of file TrkDifLineTraj.cxx.

{
  double cdip = cosDip();
  return Hep3Vector ( cos(phi0()) * cdip,
                      sin(phi0()) * cdip,
                         tanDip() * cdip);
}

Referenced by getInfo().

distTo1stError()

double TrkDifLineTraj::distTo1stError ( double  flt, double  tol, int  pathDir  ) const

virtual

Implements Trajectory.

Definition at line 133 of file TrkDifLineTraj.cxx.

{
  return 999.e4;
}

distTo2ndError()

double TrkDifLineTraj::distTo2ndError ( double  flt, double  tol, int  pathDir  ) const

virtual

Implements Trajectory.

Definition at line 139 of file TrkDifLineTraj.cxx.

{
  return 999.e4;
}

getDFInfo()

void TrkDifLineTraj::getDFInfo ( double  fltLen, DifPoint &  pos, DifVector &  dir, DifVector &  delDir  ) const

virtual

Implements TrkDifTraj.

Definition at line 163 of file TrkDifLineTraj.cxx.

{
  //Provides difNum version of information for calculation of derivatives.
 
  // Create difNumber versions of parameters
  //enum index (phi0Index(), etc) is from TrkLineParams.hh
  DifNumber phi0Df(phi0(), phi0Index()+1, nLinPrm());
  DifNumber d0Df(d0(), d0Index()+1, nLinPrm());
  DifNumber z0Df(z0(), z0Index()+1, nLinPrm());
  DifNumber tanDipDf(tanDip(), tanDipIndex()+1, nLinPrm());
  DifNumber zero(0.0, nLinPrm());
  phi0Df.setIndepPar( parameters() );
  d0Df.setIndepPar( parameters() );
  z0Df.setIndepPar( parameters() );
  tanDipDf.setIndepPar( parameters() );
  zero.setIndepPar( parameters() );
 
  DifNumber sphi0, cphi0;
  phi0Df.cosAndSin(cphi0, sphi0);
 
  DifNumber px(referencePoint().x());
  DifNumber py(referencePoint().y());
  DifNumber pz(referencePoint().z());
 
  DifNumber cdip =  1. / sqrt(1. + tanDipDf*tanDipDf);
  DifNumber xx = -d0Df * sphi0  +  flt * cphi0 * cdip + px;
  DifNumber yy =  d0Df * cphi0  +  flt * sphi0 * cdip + py;
  DifNumber zz =  z0Df  +  flt * tanDipDf * cdip + pz;
 
  pos = DifPoint(xx, yy, zz);
  dir = DifVector( cphi0 * cdip,
                   sphi0 * cdip,
                   tanDipDf * cdip);
 
  delDir = DifVector(zero, zero, zero);
}

getInfo() [1/2]

void TrkDifLineTraj::getInfo ( double  fltLen, HepPoint3D &  pos, Hep3Vector &  dir, Hep3Vector &  delDir  ) const

virtual

Implements Trajectory.

Definition at line 145 of file TrkDifLineTraj.cxx.

{
  // This could be made much more efficient!!!!!!
  pos = position(fltLen);
  dir = direction(fltLen);
  delDir = delDirect(fltLen);
}

getInfo() [2/2]

void TrkDifLineTraj::getInfo ( double  fltLen, HepPoint3D &  pos, Hep3Vector &  dir  ) const

virtual

Implements Trajectory.

Definition at line 155 of file TrkDifLineTraj.cxx.

{
  //  This could be made much more efficient!!!!!
  pos = position(fltLen);
  dir = direction(fltLen);
}

invertParams()

void TrkDifLineTraj::invertParams ( TrkParams *  params, std::vector< bool > &  flags  ) const

virtual

Implements TrkSimpTraj.

Definition at line 296 of file TrkDifLineTraj.cxx.

{
  // Inverts parameters and returns true if the parameter inversion
  // requires a change in sign of elements in the covariance matrix
  for (unsigned iparam = 0; iparam < NLINPRM; iparam++) {
    switch ( iparam ) {
    case d0Ind:  // changes sign
    case tanDipInd:  // changes sign
      params->parameter()[iparam] *= -1.0;
      flags[iparam] = true;
      break;
    case phi0Ind:  // changes by pi, but covariance matrix shouldn't change
      params->parameter()[iparam] =
        BesAngle(params->parameter()[iparam] + Constants::pi);
      flags[iparam] = false;
      break;
    case z0Ind:  // no change
      flags[iparam] = false;
    }
  }
  return;
}

nPar()

int TrkDifLineTraj::nPar ( ) const

inlinevirtual

Reimplemented from TrkSimpTraj.

Definition at line 60 of file TrkDifLineTraj.h.

{return NLINPRM;}

operator=()

TrkDifLineTraj & TrkDifLineTraj::operator=

(

const TrkDifLineTraj &

h

)

Definition at line 72 of file TrkDifLineTraj.cxx.

{
  if( &h != this ){
    Trajectory::operator=(h);
    _dtparams = *h.parameters();
    _refpoint = h._refpoint;
  }
  return *this;
}

paramFunction()

TranslateParams TrkDifLineTraj::paramFunction ( ) const

inlinevirtual

Implements TrkSimpTraj.

Definition at line 72 of file TrkDifLineTraj.h.

{ return TrkDifLineTraj::paramFunc; }

phi0()

double TrkDifLineTraj::phi0

(

)

const

Definition at line 208 of file TrkDifLineTraj.cxx.

{
 return BesAngle(parameters()->parameter()[phi0Index()]).rad();
}

Referenced by direction(), getDFInfo(), TrkLineRep::helix(), position(), and TrkLineRep::print().

position()

HepPoint3D TrkDifLineTraj::position ( double  fltLen ) const

virtual

Implements Trajectory.

Definition at line 107 of file TrkDifLineTraj.cxx.

{
  double cosd = cosDip();
  double cp = cos(phi0());
  double sp = sin(phi0());
  return HepPoint3D(-d0()*sp + f*cp*cosd + referencePoint().x(),
                  d0()*cp  + f*sp*cosd + referencePoint().y(),
                  z0() + f*tanDip()*cosd + referencePoint().z());
}

Referenced by getInfo().

tanDip()

double TrkDifLineTraj::tanDip ( ) const

inline

Definition at line 57 of file TrkDifLineTraj.h.

{return parameters()->parameter()[tanDipIndex()];}

Referenced by cosDip(), derivDeflect(), direction(), getDFInfo(), TrkLineRep::helix(), position(), and TrkLineRep::print().

visitAccept()

void TrkDifLineTraj::visitAccept ( TrkVisitor *  vis ) const

virtual

Implements TrkSimpTraj.

Definition at line 320 of file TrkDifLineTraj.cxx.

{
// Visitor access--just use the TrkVisitor class member function
  vis->trkVisitLineTraj(this);
}

z0()

double TrkDifLineTraj::z0 ( ) const

inline

Definition at line 56 of file TrkDifLineTraj.h.

{return parameters()->parameter()[z0Index()];}

Referenced by getDFInfo(), TrkLineRep::helix(), position(), and TrkLineRep::print().

---

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

• TrkDifLineTraj.h
• TrkDifLineTraj.cxx