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 ()
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
void	changePoint (const HepPoint3D &newpoint, double &fltlen)
void	setPoint (const HepPoint3D &newpoint)
TrkSimpTraj &	invert ()
bool	operator== (const TrkSimpTraj &) const
Public Member Functions inherited from TrkDifTraj
	TrkDifTraj (const double lowlim=-99999., const double hilim=99999.)
virtual	~TrkDifTraj ()
virtual void	getDFInfo2 (double fltLen, DifPoint &pos, DifVector &direction) const
Public Member Functions inherited from Trajectory
	Trajectory (double lowlim, double hilim)
virtual	~Trajectory ()
double	distTo0thError (double s, double tol, int pathDir) const
bool	validFlightDistance (double f, double tolerance=0.0) const
virtual void	setFlightRange (double newrange[2])
double	lowRange () const
double	hiRange () const
double	range () const
Public Member Functions inherited from TrkKalDeriv

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)
{
}

Referenced by clone().

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(), and 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(), and 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