#include "SurfPlane.h"

Inheritance diagram for trkf::SurfPlane:

Public Types
enum	TrackDirection { FORWARD, BACKWARD, UNKNOWN }
	Track direction enum. More...

Public Member Functions
	SurfPlane ()
	Default constructor. More...

virtual	~SurfPlane ()
	Destructor. More...

double	PointingError (const TrackVector &vec, const TrackError &err) const
	Get pointing error of track. More...

void	getStartingError (TrackError &err) const
	Get starting error matrix for Kalman filter. More...

virtual Surface *	clone () const =0
	Clone method. More...

virtual bool	isTrackValid (const TrackVector &vec) const =0
	Surface-specific tests of validity of track parameters. More...

virtual void	toLocal (const double xyz[3], double uvw[3]) const =0
	Transform global to local coordinates. More...

virtual void	toGlobal (const double uvw[3], double xyz[3]) const =0
	Transform local to global coordinates. More...

virtual TrackVector	getDiff (const TrackVector &vec1, const TrackVector &vec2) const
	Calculate difference of two track parameter vectors. More...

virtual void	getPosition (const TrackVector &vec, double xyz[3]) const =0
	Get position of track. More...

virtual TrackDirection	getDirection (const TrackVector &, TrackDirection dir=UNKNOWN) const
	Get direction of track (default UNKNOWN). More...

virtual void	getMomentum (const TrackVector &vec, double mom[3], TrackDirection dir=UNKNOWN) const =0
	Get momentum vector of track. More...

virtual bool	isParallel (const Surface &surf) const =0
	Test whether two surfaces are parallel, within tolerance. More...

virtual double	distanceTo (const Surface &surf) const =0
	Find perpendicular forward distance to a parallel surface. More...

virtual bool	isEqual (const Surface &surf) const =0
	Test two surfaces for equality, within tolerance. More...

virtual std::ostream &	Print (std::ostream &out) const =0
	Printout. More...

Detailed Description

Definition at line 28 of file SurfPlane.h.

Member Enumeration Documentation

enum trkf::Surface::TrackDirection

inherited

Track direction enum.

Enumerator
FORWARD
BACKWARD
UNKNOWN

Definition at line 54 of file Surface.h.

54 { FORWARD, BACKWARD, UNKNOWN };

trkf::Surface::BACKWARD

Definition: Surface.h:54

trkf::Surface::UNKNOWN

Definition: Surface.h:54

trkf::Surface::FORWARD

Definition: Surface.h:54

Constructor & Destructor Documentation

trkf::SurfPlane::SurfPlane ( )

Default constructor.

Definition at line 17 of file SurfPlane.cxx.

17 {}

trkf::SurfPlane::~SurfPlane ( )

virtual

Destructor.

Definition at line 20 of file SurfPlane.cxx.

20 {}

Member Function Documentation

virtual Surface* trkf::Surface::clone ( ) const

pure virtualinherited

Clone method.

Implemented in trkf::SurfXYZPlane, trkf::SurfYZLine, and trkf::SurfYZPlane.

virtual double trkf::Surface::distanceTo ( const Surface & surf ) const

pure virtualinherited

Find perpendicular forward distance to a parallel surface.

Implemented in trkf::SurfXYZPlane, trkf::SurfYZLine, and trkf::SurfYZPlane.

Referenced by trkf::Surface::getDirection().

TrackVector trkf::Surface::getDiff	(	const TrackVector &	vec1,
		const TrackVector &	vec2
	)		const

virtualinherited

Calculate difference of two track parameter vectors.

Calculate difference of two track parameter vectors. This method has a default implementation which is just the numeric difference. Surfaces that require a more sophisticated difference (e.g. phi-wrap difference) should override this method.

Reimplemented in trkf::SurfYZLine.

Definition at line 31 of file Surface.cxx.

   {
     return vec1 - vec2;
   }

virtual TrackDirection trkf::Surface::getDirection	(	const TrackVector &	,
		TrackDirection	dir = `UNKNOWN`
	)		const

inlinevirtualinherited

Get direction of track (default UNKNOWN).

Definition at line 83 of file Surface.h.

References dir, trkf::Surface::distanceTo(), trkf::Surface::getMomentum(), trkf::Surface::getStartingError(), trkf::Surface::isEqual(), trkf::Surface::isParallel(), trkf::operator<<(), trkf::Surface::PointingError(), and trkf::Surface::Print().

Referenced by trkf::SurfYZPlane::getMomentum(), and trkf::SurfXYZPlane::getMomentum().

     {
       return dir;
     }

virtual void trkf::Surface::getMomentum	(	const TrackVector &	vec,
		double	mom[3],
		TrackDirection	dir = `UNKNOWN`
	)		const

pure virtualinherited

Get momentum vector of track.

Implemented in trkf::SurfXYZPlane, trkf::SurfYZLine, and trkf::SurfYZPlane.

Referenced by trkf::Surface::getDirection().

virtual void trkf::Surface::getPosition	(	const TrackVector &	vec,
		double	xyz[3]
	)		const

pure virtualinherited

Get position of track.

Implemented in trkf::SurfXYZPlane, trkf::SurfYZLine, and trkf::SurfYZPlane.

void trkf::SurfPlane::getStartingError ( TrackError & err ) const

virtual

Get starting error matrix for Kalman filter.

Arguments:

err - Error matrix.

Implements trkf::Surface.

Definition at line 83 of file SurfPlane.cxx.

   {
     err.resize(5, false);
     err.clear();
     err(0, 0) = 1000.;
     err(1, 1) = 1000.;
     err(2, 2) = 0.25;
     err(3, 3) = 0.25;
     err(4, 4) = 10.;
   }

virtual bool trkf::Surface::isEqual ( const Surface & surf ) const

pure virtualinherited

Test two surfaces for equality, within tolerance.

Implemented in trkf::SurfXYZPlane, trkf::SurfYZLine, and trkf::SurfYZPlane.

Referenced by trkf::Surface::getDirection().

virtual bool trkf::Surface::isParallel ( const Surface & surf ) const

pure virtualinherited

Test whether two surfaces are parallel, within tolerance.

Implemented in trkf::SurfXYZPlane, trkf::SurfYZLine, and trkf::SurfYZPlane.

Referenced by trkf::Surface::getDirection().

virtual bool trkf::Surface::isTrackValid ( const TrackVector & vec ) const

pure virtualinherited

Surface-specific tests of validity of track parameters.

Implemented in trkf::SurfXYZPlane, trkf::SurfYZLine, and trkf::SurfYZPlane.

double trkf::SurfPlane::PointingError	(	const TrackVector &	vec,
		const TrackError &	err
	)		const

virtual

Get pointing error of track.

Arguments:

vec - Track parameters. err - Track error matrix.

Returns: Pointing error.

This method calculates the track pointing error based on the slope track paramers and errors (parameters 2 and 3).

Implements trkf::Surface.

Definition at line 34 of file SurfPlane.cxx.

References den.

   {
     // Get slope parameters and error matrix.
 
     double xp = vec(2);
     double yp = vec(3);
     double exx = err(2, 2);
     double eyy = err(3, 3);
     double exy = err(3, 2);
 
     // Calculate error matrix of pointing unit vector in some coordinate system.
 
     double den = 1. + xp * xp + yp * yp;
     double den3 = den * den * den;
 
     double vxx = ((1. + yp * yp) * (1. + yp * yp) * exx + xp * xp * yp * yp * eyy -
                   2. * xp * yp * (1. + yp * yp) * exy) /
                  den3;
     double vyy = (xp * xp * yp * yp * exx + (1. + xp * xp) * (1. + xp * xp) * eyy -
                   2. * xp * yp * (1. + xp * xp) * exy) /
                  den3;
     double vzz = (xp * xp * exx + yp * yp * eyy + 2. * xp * yp * exy) / den3;
 
     double vxy = (-xp * yp * (1. + yp * yp) * exx - xp * yp * (1. + xp * xp) * eyy +
                   (1. + xp * xp + yp * yp + 2. * xp * xp * yp * yp) * exy) /
                  den3;
     double vyz =
       (xp * xp * yp * exx - yp * (1. + xp * xp) * eyy - xp * (1. + xp * xp - yp * yp) * exy) / den3;
     double vxz =
       (-xp * (1. + yp * yp) * exx + xp * yp * yp * eyy - yp * (1. - xp * xp + yp * yp) * exy) /
       den3;
 
     // Calculate square root of the largest eigenvalue of error matrix.
 
     double ddd2 = vxx * vxx + vyy * vyy + vzz * vzz - 2. * vxx * vyy - 2. * vxx * vzz -
                   2. * vyy * vzz + 4. * vxy * vxy + 4. * vyz * vyz + 4. * vxz * vxz;
     double ddd = sqrt(ddd2 > 0. ? ddd2 : 0.);
     double lambda2 = 0.5 * (vxx + vyy + vzz + ddd);
     double lambda = sqrt(lambda2 > 0. ? lambda2 : 0.);
 
     return lambda;
   }