NeutTraj Class Reference

#include <NeutTraj.h>

Inheritance diagram for NeutTraj:

Public Member Functions
	NeutTraj (const NeutParams &, double lowlim=-99999., double hilim=99999., const HepPoint3D &refpoint=_theOrigin)
	NeutTraj (const NeutTraj &)
NeutTraj *	clone () const
virtual	~NeutTraj ()
NeutTraj &	operator= (const NeutTraj &)
virtual HepPoint3D	position (double fltLen) const
virtual Hep3Vector	direction (double fltLen) const
virtual Hep3Vector	delDirect (double) const
void	getInfo (double fltLen, HepPoint3D &pos, Hep3Vector &dir) const
void	getInfo (double fltLen, HepPoint3D &, Hep3Vector &dir, Hep3Vector &delDir) const
void	getDFInfo (double fltLen, DifPoint &, DifVector &dir, DifVector &delDir) const
virtual double	distTo1stError (double s, double tol, int pathDir) const
virtual double	distTo2ndError (double s, double tol, int pathDir) const
NeutParams &	params ()
const NeutParams &	params () const
int	nPar () 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 *newparams, std::vector< bool > &flags) const
TrkSimpTraj &	invert ()
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 NeutTraj.h.

Constructor & Destructor Documentation

NeutTraj() [1/2]

NeutTraj::NeutTraj	(	const NeutParams &	np,
		double	lowlim = -99999.,
		double	hilim = 99999.,
		const HepPoint3D &	refpoint = _theOrigin
	)

Definition at line 16 of file NeutTraj.cxx.

                                               : 
  TrkSimpTraj(np.parameter(), np.covariance(), lowlim,hilim)
{
}

NeutTraj() [2/2]

NeutTraj::NeutTraj

(

const NeutTraj &

n

)

Definition at line 22 of file NeutTraj.cxx.

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

~NeutTraj()

NeutTraj::~NeutTraj ( )

virtual

Definition at line 47 of file NeutTraj.cxx.

                    {
}

Member Function Documentation

clone()

NeutTraj * NeutTraj::clone ( ) const

virtual

Implements TrkSimpTraj.

Definition at line 30 of file NeutTraj.cxx.

{
  return new NeutTraj(*this);
}

curvature()

double NeutTraj::curvature ( double  fltLen ) const

inlinevirtual

Implements Trajectory.

Definition at line 62 of file NeutTraj.h.

{return 0.;}

delDirect()

Hep3Vector NeutTraj::delDirect ( double  fltLen ) const

virtual

Implements Trajectory.

Definition at line 93 of file NeutTraj.cxx.

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

Referenced by getInfo().

derivDeflect()

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

virtual

Implements TrkKalDeriv.

Definition at line 172 of file NeutTraj.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(NeutParams::_nneutprm,1); 
//
//  Compute some common things
//
  double tand = tanDip();
  double cosd = cosDip();
//
//  Go through the parameters
//
  switch (idirect) {
  case theta1:
    ddflct[NeutParams::_p][0] = 0.;
    ddflct[NeutParams::_tanDip][0] = 1.0/(cosd*cosd);
    ddflct[NeutParams::_d0][0] = 0.;
    ddflct[NeutParams::_phi0][0] =  0.;
    ddflct[NeutParams::_s0][0] =  0.;
    ddflct[NeutParams::_z0][0] = translen(fltlen) * (-1. - (tand*tand));
    break;
  case theta2:
    ddflct[NeutParams::_p][0] = 0;
    ddflct[NeutParams::_tanDip][0] = 0;
    ddflct[NeutParams::_d0][0] = -translen(fltlen)/cosd;
    ddflct[NeutParams::_phi0][0] = 1./cosd;
    ddflct[NeutParams::_z0][0] = 0.;
    ddflct[NeutParams::_s0][0] =  0.;
    break;
  }
  return ddflct;
}

derivDisplace()

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

virtual

Implements TrkKalDeriv.

Definition at line 213 of file NeutTraj.cxx.

{
//
//  This function computes the column matrix of derivatives for the change
//  in parameters for a change in the transvers 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.  The sign convention has been
//  chosen so that correlated scattering and displacement effects may be added
//
  HepMatrix ddflct(NeutParams::_nneutprm,1); 
//
//  Compute some common things
//
  double cosd = cosDip();
//
//  Go through the parameters
//
  switch (idirect) {
  case theta1:
    ddflct[NeutParams::_p][0] = 0.;
    ddflct[NeutParams::_tanDip][0] = 0.0;
    ddflct[NeutParams::_d0][0] = 0.;
    ddflct[NeutParams::_phi0][0] =  0.;
    ddflct[NeutParams::_s0][0] =  0.;
    ddflct[NeutParams::_z0][0] = 1.0/cosd;
    break;
  case theta2:
    ddflct[NeutParams::_p][0] = 0;
    ddflct[NeutParams::_tanDip][0] = 0;
    ddflct[NeutParams::_d0][0] = 1.0;
    ddflct[NeutParams::_phi0][0] = 0.0;
    ddflct[NeutParams::_z0][0] = 0.;
    ddflct[NeutParams::_s0][0] =  0.;
    break;
  }
  return ddflct;
}

derivPFract()

HepMatrix NeutTraj::derivPFract ( double  fltlen ) const

virtual

Implements TrkKalDeriv.

Definition at line 253 of file NeutTraj.cxx.

{
//  This function computes the column matrix of derivatives 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(NeutParams::_nneutprm,1);
 
//  Go through the parameters
  dmomfrac[NeutParams::_p][0] = p(); // momentum
  dmomfrac[NeutParams::_tanDip][0] = 0.0; // tanDip
  dmomfrac[NeutParams::_d0][0] = 0.0; // d0
  dmomfrac[NeutParams::_phi0][0] = 0.0; // phi0
  dmomfrac[NeutParams::_z0][0] = 0.0; // z0
  dmomfrac[NeutParams::_s0][0] = 0.0; // s0
  return dmomfrac;
}

direction()

Hep3Vector NeutTraj::direction ( double  fltLen ) const

virtual

Implements Trajectory.

Definition at line 86 of file NeutTraj.cxx.

                                   {
//-----------------------------------------------------------------
  return Hep3Vector ( cos(phi0())*cosDip(), sin(phi0())*cosDip(),
                       sinDip() );
}

Referenced by getInfo().

distTo1stError()

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

virtual

Implements Trajectory.

Definition at line 100 of file NeutTraj.cxx.

                                                       {
//-----------------------------------------------------------------
  return 9999.;
}

distTo2ndError()

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

virtual

Implements Trajectory.

Definition at line 107 of file NeutTraj.cxx.

                                                       {
//-----------------------------------------------------------------
  return 9999.;
}

getDFInfo()

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

virtual

Implements TrkDifTraj.

Definition at line 132 of file NeutTraj.cxx.

{
  //Provides difNum version of information for calculation of derivatives.
  //  Some of this duplicate code could be reduced if we had member 
  //  function templates.  Maybe.
 
  // Create difNumber versions of parameters
  //enum index (phi0(), etc) is from HelixParams.hh
  DifNumber phi0Df(phi0(), NeutParams::_phi0+1, NeutParams::_nneutprm);
  DifNumber d0Df(d0(), NeutParams::_d0+1, NeutParams::_nneutprm);
  DifNumber z0Df(z0(), NeutParams::_z0+1, NeutParams::_nneutprm);
  DifNumber tanDipDf(tanDip(), NeutParams::_tanDip+1, NeutParams::_nneutprm);
  DifNumber pDf(p(), NeutParams::_p+1, NeutParams::_nneutprm);
  // RF 03/25/99
  // the flight lenght has an error, which is the error along the trajectory 
  // in the creation point
  DifNumber s0Df(flt, NeutParams::_s0+1, NeutParams::_nneutprm);
  phi0Df.setIndepPar( parameters() );
  d0Df.setIndepPar( parameters() );
  z0Df.setIndepPar( parameters() );
  tanDipDf.setIndepPar( parameters() );
  pDf.setIndepPar( parameters() );
  s0Df.setIndepPar( parameters() );
  DifNumber dipDf = atan(tanDipDf);
 
  DifNumber sDip, cDip;
  dipDf.cosAndSin(cDip, sDip);
  DifNumber sinPhi0, cosPhi0;
  phi0Df.cosAndSin(cosPhi0, sinPhi0);
 
  DifNumber x =  cDip*s0Df*cosPhi0 - d0Df*sinPhi0;
  DifNumber y =  cDip*s0Df*sinPhi0 + d0Df*cosPhi0;
  DifNumber z =  z0Df + s0Df*sDip ;
  pos =  DifPoint(x, y, z);
  dir = DifVector( cosPhi0*cDip, sinPhi0*cDip, sDip );
  delDir = DifVector(0., 0., 0.);
}

getInfo() [1/2]

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

virtual

Implements Trajectory.

Definition at line 113 of file NeutTraj.cxx.

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

getInfo() [2/2]

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

virtual

Implements Trajectory.

Definition at line 123 of file NeutTraj.cxx.

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

invert()

TrkSimpTraj & NeutTraj::invert

(

)

invertParams()

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

virtual

Implements TrkSimpTraj.

Definition at line 338 of file NeutTraj.cxx.

{
  assert(1==0);
 
//  // 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 < NeutParams::_nneutprm; iparam++) {
//    switch ( iparam ) {
//    case NeutParams::_d0:  // changes sign
//    case NeutParams::_tanDip:  // changes sign
//      newparams->parameter()[iparam] *= -1.0;
//      flags[iparam] = true;
//      break;
//    case NeutParams::_phi0:  // changes by pi, but covariance
//                             // matrix shouldn't changed
//      newparams->parameter()[iparam] =
//  BesAngle(newparams->parameter()[iparam] + Constants::pi);
//      flags[iparam] = false;
//      break;
//    case NeutParams::_p:  // no change sign
//    case NeutParams::_z0:  // no change
//      flags[iparam] = false;
//    }
//  }
//  return;
}

nPar()

int NeutTraj::nPar ( ) const

inlinevirtual

Reimplemented from TrkSimpTraj.

Definition at line 57 of file NeutTraj.h.

{return params().nPar();}

operator=()

NeutTraj & NeutTraj::operator=

(

const NeutTraj &

n

)

Definition at line 36 of file NeutTraj.cxx.

{
  if( &n != this ){
    for(int iend=0;iend<2;iend++)
      flightrange[iend] = n.flightrange[iend];
    _dtparams = *n._np();
  }
  return *this;
}

paramFunction()

TranslateParams NeutTraj::paramFunction ( ) const

inlinevirtual

Implements TrkSimpTraj.

Definition at line 67 of file NeutTraj.h.

{ return NeutTraj::paramFunc; }

params() [1/2]

NeutParams & NeutTraj::params ( )

inline

Definition at line 54 of file NeutTraj.h.

{return *_np();}

Referenced by nPar(), and TrkCompTrk::print().

params() [2/2]

const NeutParams & NeutTraj::params ( ) const

inline

Definition at line 55 of file NeutTraj.h.

{return *_np();}

position()

HepPoint3D NeutTraj::position ( double  fltLen ) const

virtual

Implements Trajectory.

Definition at line 73 of file NeutTraj.cxx.

                                  {
//-----------------------------------------------------------------
  double cdip = cosDip();
  double sphi0 = sin(phi0());
  double cphi0 = cos(phi0());
 
  return HepPoint3D(cdip * f * cphi0 - d0()*sphi0, 
          cdip * f * sphi0 + d0()*cphi0, 
          z0() + f * tanDip() * cdip);
}

Referenced by getInfo().

visitAccept()

void NeutTraj::visitAccept ( TrkVisitor *  vis ) const

virtual

Implements TrkSimpTraj.

Definition at line 367 of file NeutTraj.cxx.

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

---

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

• NeutTraj.h
• NeutTraj.cxx