HelixTraj Class Reference

#include <HelixTraj.h>

Inheritance diagram for HelixTraj:

Public Types
enum	ParIndex {   d0Index =0 , phi0Index , omegaIndex , z0Index ,   tanDipIndex , d0Index =0 , phi0Index , omegaIndex ,   z0Index , tanDipIndex }
enum	{ NHLXPRM = 5 }
enum	ParIndex {   d0Index =0 , phi0Index , omegaIndex , z0Index ,   tanDipIndex , d0Index =0 , phi0Index , omegaIndex ,   z0Index , tanDipIndex }
enum	{ NHLXPRM = 5 }

Public Member Functions
	HelixTraj (const HepVector &, const HepSymMatrix &, double lowlim=-99999., double hilim=99999., const HepPoint3D &refpoint=_theOrigin)
	HelixTraj (const TrkExchangePar &, double lowlim=-99999., double hilim=99999., const HepPoint3D &refpoint=_theOrigin)
	HelixTraj (const TrkParams &, double lowlim=-99999., double hilim=99999., const HepPoint3D &refpoint=_theOrigin)
	HelixTraj (const HelixTraj &)
HelixTraj *	clone () const
virtual	~HelixTraj ()
HelixTraj &	operator= (const HelixTraj &)
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	distTo1stError (double s, double tol, int pathDir) const
virtual double	distTo2ndError (double s, double tol, int pathDir) const
double	curvature (double fltLen) const
HepMatrix	derivDeflect (double fltlen, deflectDirection) const
HepMatrix	derivDisplace (double fltlen, deflectDirection idir) const
HepMatrix	derivPFract (double fltlen) const
TranslateParams	paramFunction () const
void	invertParams (TrkParams *params, std::vector< bool > &flags) const
int	nPar () const
double	d0 () const
double	phi0 () const
double	omega () const
double	z0 () const
double	tanDip () const
virtual void	visitAccept (TrkVisitor *vis) const
virtual void	print (std::ostream &os) const
virtual void	printAll (std::ostream &os) const
	HelixTraj (const HepVector &, const HepSymMatrix &, double lowlim=-99999., double hilim=99999., const HepPoint3D &refpoint=_theOrigin)
	HelixTraj (const TrkExchangePar &, double lowlim=-99999., double hilim=99999., const HepPoint3D &refpoint=_theOrigin)
	HelixTraj (const TrkParams &, double lowlim=-99999., double hilim=99999., const HepPoint3D &refpoint=_theOrigin)
	HelixTraj (const HelixTraj &)
HelixTraj *	clone () const
virtual	~HelixTraj ()
HelixTraj &	operator= (const HelixTraj &)
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	distTo1stError (double s, double tol, int pathDir) const
virtual double	distTo2ndError (double s, double tol, int pathDir) const
double	curvature (double fltLen) const
HepMatrix	derivDeflect (double fltlen, deflectDirection) const
HepMatrix	derivDisplace (double fltlen, deflectDirection idir) const
HepMatrix	derivPFract (double fltlen) const
TranslateParams	paramFunction () const
void	invertParams (TrkParams *params, std::vector< bool > &flags) const
int	nPar () const
double	d0 () const
double	phi0 () const
double	omega () const
double	z0 () const
double	tanDip () const
virtual void	visitAccept (TrkVisitor *vis) const
virtual void	print (std::ostream &os) const
virtual void	printAll (std::ostream &os) 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
	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
	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
	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
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 &)
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 InstallArea/include/TrkBase/TrkBase/HelixTraj.h.

Member Enumeration Documentation

anonymous enum

anonymous enum

Enumerator
NHLXPRM

Definition at line 31 of file InstallArea/include/TrkBase/TrkBase/HelixTraj.h.

{NHLXPRM = 5};

anonymous enum

anonymous enum

Enumerator
NHLXPRM

Definition at line 31 of file Reconstruction/MdcPatRec/TrkBase/TrkBase-00-01-12/TrkBase/HelixTraj.h.

{NHLXPRM = 5};

ParIndex [1/2]

enum HelixTraj::ParIndex

Enumerator
d0Index
phi0Index
omegaIndex
z0Index
tanDipIndex
d0Index
phi0Index
omegaIndex
z0Index
tanDipIndex

Definition at line 30 of file InstallArea/include/TrkBase/TrkBase/HelixTraj.h.

{d0Index=0, phi0Index, omegaIndex, z0Index, tanDipIndex}; 

ParIndex [2/2]

enum HelixTraj::ParIndex

Enumerator
d0Index
phi0Index
omegaIndex
z0Index
tanDipIndex
d0Index
phi0Index
omegaIndex
z0Index
tanDipIndex

Definition at line 30 of file Reconstruction/MdcPatRec/TrkBase/TrkBase-00-01-12/TrkBase/HelixTraj.h.

{d0Index=0, phi0Index, omegaIndex, z0Index, tanDipIndex}; 

Constructor & Destructor Documentation

HelixTraj() [1/8]

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

Definition at line 28 of file HelixTraj.cxx.

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

HelixTraj() [2/8]

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

Definition at line 44 of file HelixTraj.cxx.

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

HelixTraj() [3/8]

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

Definition at line 51 of file HelixTraj.cxx.

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

HelixTraj() [4/8]

HelixTraj::HelixTraj

(

const HelixTraj &

h

)

Definition at line 59 of file HelixTraj.cxx.

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

~HelixTraj() [1/2]

HelixTraj::~HelixTraj ( )

virtual

Definition at line 82 of file HelixTraj.cxx.

{
}

HelixTraj() [5/8]

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

HelixTraj() [6/8]

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

HelixTraj() [7/8]

HelixTraj::HelixTraj	(	const TrkParams &	,
		double	lowlim = -99999.,
		double	hilim = 99999.,
		const HepPoint3D &	refpoint = _theOrigin
	)

HelixTraj() [8/8]

HelixTraj::HelixTraj

(

const HelixTraj &

)

~HelixTraj() [2/2]

virtual HelixTraj::~HelixTraj ( )

virtual

Member Function Documentation

clone() [1/2]

HelixTraj * HelixTraj::clone ( ) const

virtual

Implements TrkSimpTraj.

Definition at line 66 of file HelixTraj.cxx.

{
  return new HelixTraj(*this);
}

clone() [2/2]

HelixTraj * HelixTraj::clone ( ) const

virtual

Implements TrkSimpTraj.

curvature() [1/2]

double HelixTraj::curvature ( double  fltLen ) const

virtual

Implements Trajectory.

Definition at line 499 of file HelixTraj.cxx.

{
//  Compute the curvature as the magnitude of the 2nd derivative
//  of the position function with respect to the 3-d flight distance
//
  double cosd = cosDip();
//  return sqr(cosd)*fabs(omega());
  return (cosd*cosd)*fabs(omega());
}

curvature() [2/2]

double HelixTraj::curvature ( double  fltLen ) const

virtual

Implements Trajectory.

d0() [1/2]

double HelixTraj::d0 ( ) const

inline

Definition at line 79 of file InstallArea/include/TrkBase/TrkBase/HelixTraj.h.

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

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

d0() [2/2]

double HelixTraj::d0 ( ) const

inline

Definition at line 79 of file Reconstruction/MdcPatRec/TrkBase/TrkBase-00-01-12/TrkBase/HelixTraj.h.

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

delDirect() [1/2]

Hep3Vector HelixTraj::delDirect ( double  fltLen ) const

virtual

Implements Trajectory.

Definition at line 123 of file HelixTraj.cxx.

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

delDirect() [2/2]

virtual Hep3Vector HelixTraj::delDirect ( double  ) const

virtual

Implements Trajectory.

derivDeflect() [1/2]

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

virtual

Implements TrkKalDeriv.

Definition at line 365 of file HelixTraj.cxx.

{
//
//  This function computes the column matrix of derrivatives 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(NHLXPRM,1);
//
//  Compute some common things
//
  double omeg = omega();
  double tand = tanDip();
  double arcl = arc(fltlen);
  double dx = cos(arcl);
  double dy = sin(arcl);
  double cosd = cosDip();
  double darc = omeg*d0();
//
//  Go through the parameters
//
  switch (idirect) {
  case theta1:
    ddflct(omegaIndex+1,1) = omeg*tand;
//    ddflct(tanDipIndex+1,1) = 1.0/sqr(cosd);
    ddflct(tanDipIndex+1,1) = 1.0/(cosd*cosd);
    ddflct(d0Index+1,1) = (1-dx)*tand/omeg;
    ddflct(phi0Index+1,1) =  -dy*tand/(1+darc);
//    ddflct(z0Index+1,1) = - translen(fltlen) - sqr(tand)*dy/(omeg*(1+darc));
    ddflct(z0Index+1,1) = - translen(fltlen) - (tand*tand)*dy/(omeg*(1+darc));
    break;
  case theta2:
    ddflct(omegaIndex+1,1) = 0;
    ddflct(tanDipIndex+1,1) = 0;
    ddflct(d0Index+1,1) = -dy/(cosd*omeg);
    ddflct(phi0Index+1,1) = dx/(cosd*(1+darc));
    ddflct(z0Index+1,1) = -tand*(1- dx/(1+darc))/(cosd*omeg);
    break;
  }
 
  return ddflct;
}

derivDeflect() [2/2]

HepMatrix HelixTraj::derivDeflect ( double  fltlen, deflectDirection    ) const

virtual

Implements TrkKalDeriv.

derivDisplace() [1/2]

HepMatrix HelixTraj::derivDisplace ( double  fltlen, deflectDirection  idir  ) const

virtual

Implements TrkKalDeriv.

Definition at line 413 of file HelixTraj.cxx.

{
//
//  This function computes the column matrix of derrivatives 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' = (-sin(l)cos(p),-sin(l)sin(p),cos(l)) and
//  'theta2' = (-sin(p),cos(p),0).  These are by definition orthogonal and uncorrelated.
//  these displacements are correlated with the angular change above
//
  HepMatrix ddflct(NHLXPRM,1);
//
//  Compute some common things
//
  double omeg = omega();
  double tand = tanDip();
  double arcl = arc(fltlen);
  double dx = cos(arcl);
  double dy = sin(arcl);
  double cosd = cosDip();
  double sind = sinDip();
  double darc_1 = 1.0+omeg*d0();
//
//  Go through the parameters
//
  switch (idirect) {
  case theta1:
    ddflct(omegaIndex+1,1) = 0.0;
    ddflct(tanDipIndex+1,1) = 0.0;
    ddflct(d0Index+1,1) = -sind*dy;
    ddflct(phi0Index+1,1) = sind*dx*omeg/darc_1;
    ddflct(z0Index+1,1) = sind*tand*dx/darc_1 + cosd;
    break;
  case theta2:
    ddflct(omegaIndex+1,1) = 0;
    ddflct(tanDipIndex+1,1) = 0;
    ddflct(d0Index+1,1) = dx;
    ddflct(phi0Index+1,1) = dy*omeg/darc_1;
    ddflct(z0Index+1,1) = tand*dy/darc_1;
    break;
  }
 
  return ddflct;
}

derivDisplace() [2/2]

HepMatrix HelixTraj::derivDisplace ( double  fltlen, deflectDirection  idir  ) const

virtual

Implements TrkKalDeriv.

derivPFract() [1/2]

HepMatrix HelixTraj::derivPFract ( double  fltlen ) const

virtual

Implements TrkKalDeriv.

Definition at line 460 of file HelixTraj.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(NHLXPRM,1);
//
//  Compute some common things
 
  double omeg = omega();
  double tand = tanDip();
  double tranl = translen(fltlen);
  double arcl = tranl*omeg;
  double dx = cos(arcl);
  double dy = sin(arcl);
  double darc = omeg*d0();
 
//  Go through the parameters
// omega
  dmomfrac(omegaIndex+1,1) = -omeg;
// tanDip
  dmomfrac(tanDipIndex+1,1) = 0.0;
// d0
  dmomfrac(d0Index+1,1) = -(1-dx)/omeg;
// phi0
  dmomfrac(phi0Index+1,1) = dy/(1+darc);
// z0
  dmomfrac(z0Index+1,1) = -tand*(tranl-dy/((1+darc)*omeg));
//
  return dmomfrac;
}

derivPFract() [2/2]

HepMatrix HelixTraj::derivPFract ( double  fltlen ) const

virtual

Implements TrkKalDeriv.

direction() [1/2]

Hep3Vector HelixTraj::direction ( double  fltLen ) const

virtual

Implements Trajectory.

Definition at line 110 of file HelixTraj.cxx.

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

direction() [2/2]

virtual Hep3Vector HelixTraj::direction ( double  fltLen ) const

virtual

Implements Trajectory.

distTo1stError() [1/2]

double HelixTraj::distTo1stError ( double  s, double  tol, int  pathDir  ) const

virtual

Implements Trajectory.

Definition at line 133 of file HelixTraj.cxx.

{
//  return sqrt(2.*tol/fabs(omega())*(1.+sqr(tanDip())));
  return sqrt(2.*tol/fabs(omega())*(1.+tanDip()*tanDip()));
}

distTo1stError() [2/2]

virtual double HelixTraj::distTo1stError ( double  s, double  tol, int  pathDir  ) const

virtual

Implements Trajectory.

distTo2ndError() [1/2]

double HelixTraj::distTo2ndError ( double  s, double  tol, int  pathDir  ) const

virtual

Implements Trajectory.

Definition at line 140 of file HelixTraj.cxx.

{
//  return sqrt(1.+sqr(tanDip()))*cbrt(6.*tol/sqr(omega()));
  return sqrt(1.+tanDip()*tanDip())*cbrt(6.*tol/omega()*omega());
}

distTo2ndError() [2/2]

virtual double HelixTraj::distTo2ndError ( double  s, double  tol, int  pathDir  ) const

virtual

Implements Trajectory.

getDFInfo() [1/2]

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

virtual

Implements TrkDifTraj.

Definition at line 280 of file HelixTraj.cxx.

{
  //Provides difNum version of information for calculation of derivatives.
  //  All arithmetic operations have been replaced by +=, etc. versions 
  //  for speed.
 
  // Create difNumber versions of parameters
  DifNumber phi0Df(phi0(), phi0Index+1, NHLXPRM);
  DifNumber d0Df(d0(), d0Index+1, NHLXPRM);
  DifNumber z0Df(z0(), z0Index+1, NHLXPRM);
  DifNumber tanDipDf(tanDip(), tanDipIndex+1, NHLXPRM);
  DifNumber omegaDf(omega(), omegaIndex+1, NHLXPRM);
  phi0Df.setIndepPar( parameters() );
  d0Df.setIndepPar( parameters() );
  z0Df.setIndepPar( parameters() );
  tanDipDf.setIndepPar( parameters() );
  omegaDf.setIndepPar( parameters() );
  DifNumber dipDf = atan(tanDipDf);
 
  static DifNumber cDip;
  dipDf.cosAndSin(cDip, dir.z);
  static DifNumber sinPhi0, cosPhi0;
  phi0Df.cosAndSin(cosPhi0, sinPhi0);
 
  bool lref = (referencePoint().x() != 0. || referencePoint().y() != 0. ||
           referencePoint().z() != 0.);
 
  DifNumber alphaDf = cDip;
  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(dir.x, dir.y);
 
  //  DifNumber x =   (sinAlpha - sinPhi0) / omegaDf - d0Df * sinPhi0 + px;
  //  DifNumber y =  -(cosAlpha - cosPhi0) / omegaDf + d0Df * cosPhi0 + py;
 
  pos.x = dir.y;
  pos.x -= sinPhi0;
  pos.x /= omegaDf;
  DifNumber temp = d0Df;
  temp *= sinPhi0;
  pos.x -= temp;
 
  pos.y = cosPhi0;
  pos.y -= dir.x;
  pos.y /= omegaDf;
  temp = d0Df;
  temp *= cosPhi0;
  pos.y += temp;
 
  pos.z = flt;
  pos.z *= dir.z;
  pos.z += z0Df;
 
  if (lref) {
    DifNumber px(referencePoint().x());
    DifNumber py(referencePoint().y());
    DifNumber pz(referencePoint().z());
    pos.x += px;
    pos.y += py;
    pos.z += pz;
  }
 
  delDir.x = -omegaDf;
  delDir.x *= cDip;
  delDir.x *= cDip;
  delDir.x *= dir.y;
 
  delDir.y =  omegaDf;
  delDir.y *= cDip;
  delDir.y *= cDip;
  delDir.y *= dir.x;
 
  delDir.z = 0.;
 
  dir.x *= cDip;
  dir.y *= cDip;
}

getDFInfo() [2/2]

virtual void HelixTraj::getDFInfo ( double  fltLen, DifPoint &  , DifVector &  dir, DifVector &  delDir  ) const

virtual

Implements TrkDifTraj.

getDFInfo2() [1/2]

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

virtual

Reimplemented from TrkDifTraj.

Definition at line 207 of file HelixTraj.cxx.

{
  //Provides difNum version of information for calculation of derivatives.
  //  All arithmetic operations have been replaced by +=, etc. versions 
  //  for speed.
 
  // Create difNumber versions of parameters
  DifNumber phi0Df(phi0(), phi0Index+1, NHLXPRM);
  phi0Df.setIndepPar( parameters() );
  DifNumber d0Df(d0(), d0Index+1, NHLXPRM);
  d0Df.setIndepPar( parameters() );
  DifNumber z0Df(z0(), z0Index+1, NHLXPRM);
  z0Df.setIndepPar( parameters() );
  DifNumber tanDipDf(tanDip(), tanDipIndex+1, NHLXPRM);
  tanDipDf.setIndepPar( parameters() );
  DifNumber omegaDf(omega(), omegaIndex+1, NHLXPRM);
  omegaDf.setIndepPar( parameters() );
 
  DifNumber dipDf = atan(tanDipDf);
 
  static DifNumber cDip;
  dipDf.cosAndSin(cDip, dir.z);
  static DifNumber sinPhi0, cosPhi0;
  phi0Df.cosAndSin(cosPhi0, sinPhi0);
 
  bool lref = (referencePoint().x() != 0. || referencePoint().y() != 0. ||
               referencePoint().z() != 0.);
 
  DifNumber alphaDf = cDip;
  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(dir.x, dir.y);
 
  //  DifNumber x =   (sinAlpha - sinPhi0) / omegaDf - d0Df * sinPhi0 + px;
  //  DifNumber y =  -(cosAlpha - cosPhi0) / omegaDf + d0Df * cosPhi0 + py;
 
  pos.x = dir.y;
  pos.x -= sinPhi0;
  pos.x /= omegaDf;
  DifNumber temp = d0Df;
  temp *= sinPhi0;
  pos.x -= temp;
 
  pos.y = cosPhi0;
  pos.y -= dir.x;
  pos.y /= omegaDf;
  temp = d0Df;
  temp *= cosPhi0;
  pos.y += temp;
 
  pos.z = flt;
  pos.z *= dir.z;
  pos.z += z0Df;
 
  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 *= cDip;
  dir.y *= cDip;
}

getDFInfo2() [2/2]

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

virtual

Reimplemented from TrkDifTraj.

getInfo() [1/4]

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

virtual

Implements Trajectory.

Definition at line 147 of file HelixTraj.cxx.

{
  //  double ang = angle(fltLen);
  double cDip = cosDip();
  double sDip = tanDip() * cDip;
  double phi00 = parameters()->parameter()[phi0Index];  // Don't normalize
  double ang = phi00 + cDip*fltLen*omega(); 
  double cang = cos(ang);
  double sang = sin(ang);
  double sphi0 = sin(phi00);
  double cphi0 = cos(phi00);
 
  double xt = (sang - sphi0)/omega() - d0()*sphi0 +
    referencePoint().x();
  double yt = -(cang - cphi0)/omega() + d0()*cphi0 +
    referencePoint().y();
  double zt = z0() + fltLen*sDip + referencePoint().z();
  pos.setX(xt);
  pos.setY(yt);
  pos.setZ(zt);
 
  dir.setX(cang * cDip);
  dir.setY(sang * cDip);
  dir.setZ(sDip);
 
  double delX = -omega() * cDip * cDip * sang;
  double delY =  omega() * cDip * cDip * cang;
  delDir.setX(delX);
  delDir.setY(delY);
  delDir.setZ(0.0);
}

getInfo() [2/4]

virtual void HelixTraj::getInfo ( double  fltLen, HepPoint3D &  , Hep3Vector &  dir, Hep3Vector &  delDir  ) const

virtual

Implements Trajectory.

getInfo() [3/4]

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

virtual

Implements Trajectory.

Definition at line 181 of file HelixTraj.cxx.

{
  //  double ang = angle(fltLen);
  double cDip = cosDip();
  double sDip = tanDip() * cDip;
  double phi00 = parameters()->parameter()[phi0Index];  // Don't normalize
  double ang = phi00 + cDip*fltLen*omega(); 
  double cang = cos(ang);
  double sang = sin(ang);
  double sphi0 = sin(phi00);
  double cphi0 = cos(phi00);
 
  double xt = (sang - sphi0)/omega() - d0()*sphi0 +
    referencePoint().x();
  double yt = -(cang - cphi0)/omega() + d0()*cphi0 +
    referencePoint().y();
  double zt = z0() + fltLen*sDip + referencePoint().z();
  pos.setX(xt);
  pos.setY(yt);
  pos.setZ(zt);
 
  dir.setX(cang * cDip);
  dir.setY(sang * cDip);
  dir.setZ(sDip);
}

getInfo() [4/4]

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

virtual

Implements Trajectory.

invertParams() [1/2]

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

virtual

Implements TrkSimpTraj.

Definition at line 598 of file HelixTraj.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 < NHLXPRM; iparam++) {
    switch ( iparam ) {
    case d0Index:  // changes sign
    case omegaIndex:  // changes sign
    case tanDipIndex:  // changes sign
      params->parameter()[iparam] *= -1.0;
      flags[iparam] = true;
      break;
    case phi0Index:  // changes by pi, but covariance matrix shouldn't change
      params->parameter()[iparam] =
        BesAngle(params->parameter()[iparam] + Constants::pi);
      flags[iparam] = false;
      break;
    case z0Index:  // nochange
      flags[iparam] = false;
    }
  }
  return;
}

invertParams() [2/2]

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

virtual

Implements TrkSimpTraj.

nPar() [1/2]

int HelixTraj::nPar ( ) const

inlinevirtual

Reimplemented from TrkSimpTraj.

Definition at line 75 of file InstallArea/include/TrkBase/TrkBase/HelixTraj.h.

{ return NHLXPRM;}

nPar() [2/2]

int HelixTraj::nPar ( ) const

inlinevirtual

Reimplemented from TrkSimpTraj.

Definition at line 75 of file Reconstruction/MdcPatRec/TrkBase/TrkBase-00-01-12/TrkBase/HelixTraj.h.

{ return NHLXPRM;}

omega() [1/2]

double HelixTraj::omega ( ) const

inline

Definition at line 81 of file InstallArea/include/TrkBase/TrkBase/HelixTraj.h.

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

Referenced by curvature(), delDirect(), derivDeflect(), derivDisplace(), derivPFract(), distTo1stError(), distTo2ndError(), getDFInfo(), getDFInfo2(), getInfo(), HelixTraj(), and position().

omega() [2/2]

double HelixTraj::omega ( ) const

inline

Definition at line 81 of file Reconstruction/MdcPatRec/TrkBase/TrkBase-00-01-12/TrkBase/HelixTraj.h.

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

operator=() [1/2]

HelixTraj & HelixTraj::operator=

(

const HelixTraj &

h

)

Definition at line 72 of file HelixTraj.cxx.

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

operator=() [2/2]

HelixTraj & HelixTraj::operator=

(

const HelixTraj &

)

paramFunction() [1/2]

TranslateParams HelixTraj::paramFunction ( ) const

inlinevirtual

Implements TrkSimpTraj.

Definition at line 69 of file InstallArea/include/TrkBase/TrkBase/HelixTraj.h.

{ return HelixTraj::paramFunc; }

paramFunction() [2/2]

TranslateParams HelixTraj::paramFunction ( ) const

inlinevirtual

Implements TrkSimpTraj.

Definition at line 69 of file Reconstruction/MdcPatRec/TrkBase/TrkBase-00-01-12/TrkBase/HelixTraj.h.

{ return HelixTraj::paramFunc; }

phi0() [1/2]

double HelixTraj::phi0

(

)

const

Definition at line 510 of file HelixTraj.cxx.

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

Referenced by getDFInfo(), and getDFInfo2().

phi0() [2/2]

double HelixTraj::phi0

(

)

const

position() [1/2]

HepPoint3D HelixTraj::position ( double  fltLen ) const

virtual

Implements Trajectory.

Definition at line 93 of file HelixTraj.cxx.

{
  double cDip = cosDip();
  double sDip = tanDip() * cDip;
  double phi00 = parameters()->parameter()[phi0Index];  // Don't normalize
  double ang = phi00 + cDip*f*omega();
  double cang = cos(ang);
  double sang = sin(ang);
  double sphi0 = sin(phi00);
  double cphi0 = cos(phi00);
 
  return HepPoint3D((sang - sphi0)/omega() - d0()*sphi0+referencePoint().x(),
                 -(cang - cphi0)/omega() + d0()*cphi0+referencePoint().y(),
                 z0() + f*sDip                       +referencePoint().z());
}

position() [2/2]

virtual HepPoint3D HelixTraj::position ( double  fltLen ) const

virtual

Implements Trajectory.

print() [1/2]

virtual void HelixTraj::print ( std::ostream &  os ) const

virtual

Reimplemented from TrkSimpTraj.

print() [2/2]

virtual void HelixTraj::print ( std::ostream &  os ) const

virtual

Reimplemented from TrkSimpTraj.

printAll() [1/2]

virtual void HelixTraj::printAll ( std::ostream &  os ) const

virtual

Reimplemented from TrkSimpTraj.

printAll() [2/2]

virtual void HelixTraj::printAll ( std::ostream &  os ) const

virtual

Reimplemented from TrkSimpTraj.

tanDip() [1/2]

double HelixTraj::tanDip ( ) const

inline

Definition at line 83 of file InstallArea/include/TrkBase/TrkBase/HelixTraj.h.

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

Referenced by derivDeflect(), derivDisplace(), derivPFract(), direction(), distTo1stError(), distTo2ndError(), getDFInfo(), getDFInfo2(), getInfo(), and position().

tanDip() [2/2]

double HelixTraj::tanDip ( ) const

inline

Definition at line 83 of file Reconstruction/MdcPatRec/TrkBase/TrkBase-00-01-12/TrkBase/HelixTraj.h.

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

visitAccept() [1/2]

void HelixTraj::visitAccept ( TrkVisitor *  vis ) const

virtual

Implements TrkSimpTraj.

Definition at line 591 of file HelixTraj.cxx.

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

visitAccept() [2/2]

virtual void HelixTraj::visitAccept ( TrkVisitor *  vis ) const

virtual

Implements TrkSimpTraj.

z0() [1/2]

double HelixTraj::z0 ( ) const

inline

Definition at line 82 of file InstallArea/include/TrkBase/TrkBase/HelixTraj.h.

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

Referenced by getDFInfo(), getDFInfo2(), getInfo(), and position().

z0() [2/2]

double HelixTraj::z0 ( ) const

inline

Definition at line 82 of file Reconstruction/MdcPatRec/TrkBase/TrkBase-00-01-12/TrkBase/HelixTraj.h.

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

---

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

• InstallArea/include/TrkBase/TrkBase/HelixTraj.h
• Reconstruction/MdcPatRec/TrkBase/TrkBase-00-01-12/TrkBase/HelixTraj.h
• HelixTraj.cxx