TrkCircleTraj Class Reference

#include <TrkCircleTraj.h>

Inheritance diagram for TrkCircleTraj:

Public Member Functions
	TrkCircleTraj (const HepVector &, const HepSymMatrix &, double lowlim=-99999., double hilim=99999., const HepPoint3D &refpoint=_theOrigin)
	TrkCircleTraj (const TrkExchangePar &, double lowlim=-99999., double hilim=99999., const HepPoint3D &refpoint=_theOrigin)
	TrkCircleTraj (const TrkCircleTraj &)
TrkCircleTraj *	clone () const
	~TrkCircleTraj ()
TrkCircleTraj &	operator= (const TrkCircleTraj &)
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 void	getDFInfo2 (double fltLen, DifPoint &pos, DifVector &dir) const
virtual double	curvature (double fltLen) 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
int	nPar () const
virtual void	visitAccept (TrkVisitor *vis) const
double	d0 () const
double	phi0 () const
double	omega () 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

Static Public Member Functions
static int	d0Index ()
static int	phi0Index ()
static int	omegaIndex ()
static int	nCirPrm ()

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 27 of file TrkCircleTraj.h.

Constructor & Destructor Documentation

TrkCircleTraj() [1/3]

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

Definition at line 32 of file TrkCircleTraj.cxx.

                                                                          :
  TrkSimpTraj(pvec, pcov, lowlim,hilim,refpoint)
//-----------------------------------------------------------------
{ 
  //  Make sure the dimensions of the input matrix and vector are correct
  if( pvec.num_row() != nCirPrm() ||
      pcov.num_row() != nCirPrm()  ){
    std::cout<<"ErrMsg(fatal)" << 
      "CircleTraj: incorrect constructor vector/matrix dimension" << std::endl;
  }
 
  if (omega() == 0.0) parameters()->parameter()[omegaIndex()] = 1.e-9; 
}

TrkCircleTraj() [2/3]

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

Definition at line 49 of file TrkCircleTraj.cxx.

                                                                          :
  TrkSimpTraj(inpar.params(), inpar.covariance(), lowlim,hilim,refpoint) {
//-----------------------------------------------------------------
 
  if (omega() == 0.0) parameters()->parameter()[omegaIndex()] = 1.e-9;   
}

TrkCircleTraj() [3/3]

TrkCircleTraj::TrkCircleTraj

(

const TrkCircleTraj &

h

)

Definition at line 57 of file TrkCircleTraj.cxx.

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

~TrkCircleTraj()

TrkCircleTraj::~TrkCircleTraj

(

)

Definition at line 81 of file TrkCircleTraj.cxx.

{
}

Member Function Documentation

clone()

TrkCircleTraj * TrkCircleTraj::clone ( ) const

virtual

Implements TrkSimpTraj.

Definition at line 65 of file TrkCircleTraj.cxx.

{
  return new TrkCircleTraj(*this);
}

curvature()

double TrkCircleTraj::curvature ( double  fltLen ) const

virtual

Implements Trajectory.

Definition at line 421 of file TrkCircleTraj.cxx.

{
//  Compute the curvature as the magnitude of the 2nd derrivative
//  of the position function with respect to the 3-d flight distance
 
  return fabs(omega());
}

d0()

double TrkCircleTraj::d0 ( ) const

inline

Definition at line 88 of file TrkCircleTraj.h.

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

Referenced by derivDeflect(), derivDisplace(), derivPFract(), getDFInfo(), getDFInfo2(), getInfo(), TrkCircleRep::helix(), position(), and TrkCircleRep::print().

d0Index()

static int TrkCircleTraj::d0Index ( )

inlinestatic

Definition at line 82 of file TrkCircleTraj.h.

{return d0Ind;}

Referenced by d0(), derivDeflect(), derivDisplace(), derivPFract(), getDFInfo(), and getDFInfo2().

delDirect()

Hep3Vector TrkCircleTraj::delDirect ( double  fltLen ) const

virtual

Implements Trajectory.

Definition at line 123 of file TrkCircleTraj.cxx.

{
  double delX = -omega() * sin(angle(fltLen));
  double delY =  omega() * cos(angle(fltLen));
  return Hep3Vector(delX, delY, 0.0);
}

derivDeflect()

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

virtual

Implements TrkKalDeriv.

Definition at line 322 of file TrkCircleTraj.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)) can sometimes be added, as scattering in these
//  are uncorrelated.
 
  HepMatrix ddflct(nCirPrm(),1,0); // initialize with zeros
 
//  Compute some common things
 
  double omeg = omega();
  double arcl = arc(fltlen);
  double dx = cos(arcl);
  double dy = sin(arcl);
  double darc = omeg*d0();
 
//  Go through the parameters
 
  switch (idirect) {
  case theta1:
    break;
  case theta2:
    ddflct[d0Index()][0] = -dy/(omeg);
    ddflct[phi0Index()][0] = dx/(1+darc);
  }
 
  return ddflct;
}

derivDisplace()

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

virtual

Implements TrkKalDeriv.

Definition at line 355 of file TrkCircleTraj.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)) can sometimes be added, as scattering in these
//  are uncorrelated.
 
  HepMatrix ddflct(nCirPrm(),1,0); // initialize with zeros
 
//  Compute some common things
 
  double omeg = omega();
  double arcl = arc(fltlen);
  double dx = cos(arcl);
  double dy = sin(arcl);
  double darc = omeg*d0();
 
//  Go through the parameters
 
  switch (idirect) {
  case theta1:
    break;
  case theta2:
    ddflct[d0Index()][0] = dx;
    ddflct[phi0Index()][0] = dy*omeg/(1+darc);
  }
 
  return ddflct;
}

derivPFract()

HepMatrix TrkCircleTraj::derivPFract ( double  fltlen ) const

virtual

Implements TrkKalDeriv.

Definition at line 388 of file TrkCircleTraj.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.
//
//  For a helix, dp/P = -domega/omega,
//  dParam/d(domega/omega) = omega*dParam/ddomega
 
  HepMatrix dmomfrac(nCirPrm(),1);
 
//  Compute some common things
 
  double omeg = omega();
  double arcl = arc(fltlen);
  double dx = cos(arcl);
  double dy = sin(arcl);
  double darc = omeg*d0();
 
//  Go through the parameters
// omega
  dmomfrac[omegaIndex()][0] = -omeg;
// d0
  dmomfrac[d0Index()][0] = -(1-dx)/omeg;
// phi0
  dmomfrac[phi0Index()][0] = dy/(1+darc);
 
  return dmomfrac;
}

direction()

Hep3Vector TrkCircleTraj::direction ( double  fltLen ) const

virtual

Implements Trajectory.

Definition at line 111 of file TrkCircleTraj.cxx.

{
  // Angle formed by tangent vector after
  // being rotated 'arclength' around orbit.
  double alpha = angle( f );
  // Construct 3-D tangent vector of unit magnitude.
  return Hep3Vector ( cos(alpha),
              sin(alpha),
                     0.0);
}

distTo1stError()

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

virtual

Implements Trajectory.

Definition at line 131 of file TrkCircleTraj.cxx.

{
  double arg = 2. * tol / fabs(omega());
  assert (arg >= 0.);
  return sqrt(arg);
}

distTo2ndError()

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

virtual

Implements Trajectory.

Definition at line 139 of file TrkCircleTraj.cxx.

{
  //return pow(6.*tol / sqr(omega()), 0.33333333);//yzhang changed sqr
  return pow(6.*tol / (omega()*omega()), 0.33333333);
}

getDFInfo()

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

virtual

Implements TrkDifTraj.

Definition at line 253 of file TrkCircleTraj.cxx.

{
  //Provides difNum version of information for calculation of derivatives.
 
  // Create difNumber versions of parameters
  DifNumber phi0Df(phi0(), phi0Index()+1, nCirPrm());
  DifNumber d0Df(d0(), d0Index()+1, nCirPrm());
  DifNumber omegaDf(omega(), omegaIndex()+1, nCirPrm());
  phi0Df.setIndepPar( parameters() );
  d0Df.setIndepPar( parameters() );
  omegaDf.setIndepPar( parameters() );
 
  DifNumber sinPhi0, cosPhi0;
  phi0Df.cosAndSin(cosPhi0, sinPhi0);
 
  DifNumber alphaDf = omegaDf;
  alphaDf *= flt;
  alphaDf += phi0Df;
 
  // This is not the prettiest line imaginable for this operation:
  alphaDf.mod(-Constants::pi, Constants::pi);
  DifNumber sinAlpha, cosAlpha;
  alphaDf.cosAndSin(cosAlpha, sinAlpha);
 
  //  DifNumber x =  (sinAlpha - sinPhi0) / omegaDf - d0Df * sinPhi0 + px;
  //  DifNumber y = -(cosAlpha - cosPhi0) / omegaDf + d0Df * cosPhi0 + py;
 
  DifNumber x = sinAlpha;
  x -= sinPhi0;
  x /= omegaDf;
  x -= (d0Df * sinPhi0);
 
  DifNumber y = -cosAlpha;
  y += cosPhi0;
  y /= omegaDf;
  y += (d0Df * cosPhi0);
 
 
  static DifNumber zNull(0.);
 
  pos.x = x;
  pos.y = y;
  pos.z = zNull;
 
  bool lref = (referencePoint().x() != 0. || referencePoint().y() != 0. ||
           referencePoint().z() != 0.);
  if (lref) {
    DifNumber px(referencePoint().x());
    DifNumber py(referencePoint().y());
    DifNumber pz(referencePoint().z());
    pos.x += px;
    pos.y += py;
    pos.z += pz;
  }
 
  dir.x = cosAlpha;
  dir.y = sinAlpha;
  dir.z = 0.;
 
  delDir.x = -omegaDf;
  delDir.x *= sinAlpha;
 
  delDir.y = omegaDf;
  delDir.y *= cosAlpha;
 
  delDir.z = 0.;
}

getDFInfo2()

void TrkCircleTraj::getDFInfo2 ( double  fltLen, DifPoint &  pos, DifVector &  dir  ) const

virtual

Reimplemented from TrkDifTraj.

Definition at line 191 of file TrkCircleTraj.cxx.

{
  //Provides difNum version of information for calculation of derivatives.
 
  // Create difNumber versions of parameters
  DifNumber phi0Df(phi0(), phi0Index()+1, nCirPrm());
  DifNumber d0Df(d0(), d0Index()+1, nCirPrm());
  DifNumber omegaDf(omega(), omegaIndex()+1, nCirPrm());
  phi0Df.setIndepPar( parameters() );
  d0Df.setIndepPar( parameters() );
  omegaDf.setIndepPar( parameters() );
 
  DifNumber sinPhi0, cosPhi0;
  phi0Df.cosAndSin(cosPhi0, sinPhi0);
 
  DifNumber alphaDf = omegaDf;
  alphaDf *= fltLen;
  alphaDf += phi0Df;
 
  // This is not the prettiest line imaginable for this operation:
  alphaDf.mod(-Constants::pi, Constants::pi);
  DifNumber sinAlpha, cosAlpha;
  alphaDf.cosAndSin(cosAlpha, sinAlpha);
 
  //  DifNumber x =  (sinAlpha - sinPhi0) / omegaDf - d0Df * sinPhi0 + px;
  //  DifNumber y = -(cosAlpha - cosPhi0) / omegaDf + d0Df * cosPhi0 + py;
 
  DifNumber x = sinAlpha;
  x -= sinPhi0;
  x /= omegaDf;
  x -= (d0Df * sinPhi0);
 
  DifNumber y = -cosAlpha;
  y += cosPhi0;
  y /= omegaDf;
  y += (d0Df * cosPhi0);
 
 
  static DifNumber zNull(0.);
 
  pos.x = x;
  pos.y = y;
  pos.z = zNull;
 
  bool lref = (referencePoint().x() != 0. || referencePoint().y() != 0. ||
           referencePoint().z() != 0.);
  if (lref) {
    DifNumber px(referencePoint().x());
    DifNumber py(referencePoint().y());
    DifNumber pz(referencePoint().z());
    pos.x += px;
    pos.y += py;
    pos.z += pz;
  }
 
  dir.x = cosAlpha;
  dir.y = sinAlpha;
  dir.z = 0.;
}

getInfo() [1/2]

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

virtual

Implements Trajectory.

Definition at line 146 of file TrkCircleTraj.cxx.

{
  double sphi0 = sin(phi0());
  double cphi0 = cos(phi0());
  double ang = angle(flt);
  double cang = cos(ang);
  double sang = sin(ang);
  double xt =   (sang-sphi0)/omega() - d0()*sphi0 + referencePoint().x();
  double yt =  -(cang-cphi0)/omega() + d0()*cphi0 + referencePoint().y();
 
  pos.setX(xt);
  pos.setY(yt);
  pos.setZ(referencePoint().z());
 
  dir.setX(cang);
  dir.setY(sang);
  dir.setZ(0.0);
 
  delDir.setX(-omega()*sang);
  delDir.setY( omega()*cang);
  delDir.setZ(0.);
}

getInfo() [2/2]

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

virtual

Implements Trajectory.

Definition at line 171 of file TrkCircleTraj.cxx.

{
  double sphi0 = sin(phi0());
  double cphi0 = cos(phi0());
  double ang = angle(fltLen);
  double cang = cos(ang);
  double sang = sin(ang);
  double xt =   (sang-sphi0)/omega() - d0()*sphi0 + referencePoint().x();
  double yt =  -(cang-cphi0)/omega() + d0()*cphi0 + referencePoint().y();
 
  pos.setX(xt);
  pos.setY(yt);
  pos.setZ(referencePoint().z());
 
  dir.setX(cang);
  dir.setY(sang);
  dir.setZ(0.0);
}

invertParams()

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

virtual

Implements TrkSimpTraj.

Definition at line 495 of file TrkCircleTraj.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 < NCIRPAR; iparam++) {
    switch ( iparam ) {
    case d0Ind:  // changes sign
    case omegaInd:  // 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;
    }
  }
}

nCirPrm()

static int TrkCircleTraj::nCirPrm ( )

inlinestatic

Definition at line 85 of file TrkCircleTraj.h.

{return NCIRPAR;}

Referenced by derivDeflect(), derivDisplace(), derivPFract(), getDFInfo(), getDFInfo2(), and TrkCircleTraj().

nPar()

int TrkCircleTraj::nPar ( ) const

inlinevirtual

Reimplemented from TrkSimpTraj.

Definition at line 69 of file TrkCircleTraj.h.

{return NCIRPAR;}

omega()

double TrkCircleTraj::omega ( ) const

inline

Definition at line 90 of file TrkCircleTraj.h.

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

Referenced by curvature(), delDirect(), derivDeflect(), derivDisplace(), derivPFract(), distTo1stError(), distTo2ndError(), getDFInfo(), getDFInfo2(), getInfo(), TrkCircleRep::helix(), position(), TrkCircleRep::print(), and TrkCircleTraj().

omegaIndex()

static int TrkCircleTraj::omegaIndex ( )

inlinestatic

Definition at line 84 of file TrkCircleTraj.h.

{return omegaInd;}

Referenced by derivPFract(), getDFInfo(), getDFInfo2(), omega(), and TrkCircleTraj().

operator=()

TrkCircleTraj & TrkCircleTraj::operator=

(

const TrkCircleTraj &

h

)

Definition at line 71 of file TrkCircleTraj.cxx.

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

paramFunction()

TranslateParams TrkCircleTraj::paramFunction ( ) const

inlinevirtual

Implements TrkSimpTraj.

Definition at line 64 of file TrkCircleTraj.h.

{ return TrkCircleTraj::paramFunc; }

phi0()

double TrkCircleTraj::phi0

(

)

const

Definition at line 430 of file TrkCircleTraj.cxx.

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

Referenced by getDFInfo(), getDFInfo2(), getInfo(), TrkCircleRep::helix(), position(), and TrkCircleRep::print().

phi0Index()

static int TrkCircleTraj::phi0Index ( )

inlinestatic

Definition at line 83 of file TrkCircleTraj.h.

{return phi0Ind;}

Referenced by derivDeflect(), derivDisplace(), derivPFract(), getDFInfo(), getDFInfo2(), and phi0().

position()

HepPoint3D TrkCircleTraj::position ( double  fltLen ) const

virtual

Implements Trajectory.

Definition at line 100 of file TrkCircleTraj.cxx.

{
  double sphi0 = sin(phi0());
  double cphi0 = cos(phi0());
  double ang = angle(f);
  double xt =  ( sin(ang)-sphi0)/omega() - d0()*sphi0 + referencePoint().x();
  double yt =  -(cos(ang)-cphi0)/omega() + d0()*cphi0 + referencePoint().y();
  return HepPoint3D(xt, yt, referencePoint().z());
}

visitAccept()

void TrkCircleTraj::visitAccept ( TrkVisitor *  vis ) const

virtual

Implements TrkSimpTraj.

Definition at line 522 of file TrkCircleTraj.cxx.

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

---

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

• TrkCircleTraj.h
• TrkCircleTraj.cxx