Linear fit of straight line parallel to z-axis to set of hits (objects with position and time). More...

#include <JLine1ZEstimator.hh>

Inheritance diagram for JFIT::JEstimator< JLine1Z >:

Public Types
typedef double JLine1Z::*	parameter_type

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

template<class T >
	JEstimator (T __begin, T __end)
	Constructor. More...

template<class T >
const JEstimator< JLine1Z > &	operator() (T __begin, T __end)
	Fit. More...

template<class T >
void	update (T __begin, T __end, const JMatrixNZ &A)
	Update track parameters using updated co-variance matrix (e.g. matrix with one hit switched off). More...

void	move (const double step, const double velocity)
	Move vertex along this line with given velocity. More...

void	setZ (const double z, const double velocity)
	Set z-position of vertex. More...

double	getDistanceSquared (const JVector3D &pos) const
	Get distance squared. More...

double	getDistance (const JVector3D &pos) const
	Get distance. More...

double	getT (const JVector3D &pos) const
	Get arrival time of Cherenkov light at given position. More...

double	getT () const
	Get time. More...

double	getZ (const JPosition3D &pos) const
	Get point of emission of Cherenkov light along muon path. More...

double	getZ () const
	Get z position. More...

JVersor3D	getDirection (const JVector3D &pos) const
	Get photon direction of Cherenkov light on PMT. More...

double	getDot (const JAxis3D &axis) const
	Get cosine angle of impact of Cherenkov light on PMT. More...

double	getDot (const JAngle3D &angle) const
	Get dot product. More...

double	getDot (const JVersor3D &dir) const
	Get dot product. More...

double	getDot (const JVersor3Z &dir) const
	Get dot product. More...

double	getDot (const JVector3D &vector) const
	Get dot product. More...

JVertex3D &	negate ()
	Prefix unary minus. More...

JVertex3D &	add (const JVertex3D &value)
	Addition operator. More...

JVector3D &	add (const JVector3D &vector)
	Add vector. More...

JTime &	add (const JTime &value)
	Addition operator. More...

JVertex3D &	sub (const JVertex3D &value)
	Subtraction operator. More...

JVector3D &	sub (const JVector3D &vector)
	Subtract vector. More...

JTime &	sub (const JTime &value)
	Subtraction operator. More...

JVertex3D &	mul (const double value)
	Multiplication operator. More...

JFirst_t &	mul (const JSecond_t &object)
	Multiply with object. More...

JFirst_t &	mul (const JSecond_t &object)
	Multiply with object. More...

JVertex3D &	div (const double value)
	Division operator. More...

const JPosition3D &	getPosition () const
	Get position. More...

JPosition3D &	getPosition ()
	Get position. More...

void	setPosition (const JVector3D &pos)
	Set position. More...

	operator JAngle3D () const
	Type conversion operator. More...

	operator JVersor3D () const
	Type conversion operator. More...

JPosition3D &	rotate (const JRotation3D &R)
	Rotate. More...

JPosition3D &	rotate (const JRotation3X &R)
	Rotate around X-axis. More...

JPosition3D &	rotate (const JRotation3Y &R)
	Rotate around Y-axis. More...

JPosition3D &	rotate (const JRotation3Z &R)
	Rotate around Z-axis. More...

JPosition3D &	rotate (const JQuaternion3D &Q)
	Rotate. More...

JPosition3D &	rotate_back (const JRotation3D &R)
	Rotate back. More...

JPosition3D &	rotate_back (const JRotation3X &R)
	Rotate back around X-axis. More...

JPosition3D &	rotate_back (const JRotation3Y &R)
	Rotate back around Y-axis. More...

JPosition3D &	rotate_back (const JRotation3Z &R)
	Rotate back around Z-axis. More...

JPosition3D &	rotate_back (const JQuaternion3D &Q)
	Rotate back. More...

void	transform (const JRotation3D &R, const JVector3D &pos)
	Transform position. More...

JVector3D &	transform (const JMatrix3D &T)
	Transform. More...

void	transform_back (const JRotation3D &R, const JVector3D &pos)
	Transform back position. More...

	operator JVector2D () const
	Type conversion operator. More...

double	getX () const
	Get x position. More...

double	getY () const
	Get y position. More...

bool	equals (const JVector3D &vector, const double precision=std::numeric_limits< double >::min()) const
	Check equality. More...

double	getLengthSquared () const
	Get length squared. More...

double	getLength () const
	Get length. More...

JVector3D &	getCross (const JVector3D &first, const JVector3D &second)
	Get cross product. More...

Static Public Member Functions
static parameter_type	pX ()

static parameter_type	pY ()

static parameter_type	pT ()

Public Attributes
JMATH::JMatrix3S	V
	co-variance matrix of fit parameters More...

JMATH::JSVD3D	svd

Static Public Attributes
static const int	NUMBER_OF_PARAMETERS = 3
	number of parameters of fit More...

static double	MINIMAL_SVD_WEIGHT = 1.0e-4
	minimal SVD weight. More...

Protected Attributes
double	__x

double	__y

double	__z

double	__t

Detailed Description

template<>
class JFIT::JEstimator< JLine1Z >

Linear fit of straight line parallel to z-axis to set of hits (objects with position and time).

$\begin{center}\setlength{\unitlength}{0.7cm}\begin{picture}(8,12) \put( 1.0, 0.5){\vector(0,1){9}} \put( 1.0,10.0){\makebox(0,0){$z$}} \put( 1.0, 0.0){\makebox(0,0){muon}} \put( 1.0, 8.0){\line(1,0){6}} \put( 4.0, 8.5){\makebox(0,0)[c]{$R$}} \multiput( 1.0, 2.0)(0.5, 0.5){12}{\qbezier(0.0,0.0)(-0.1,0.35)(0.25,0.25)\qbezier(0.25,0.25)(0.6,0.15)(0.5,0.5)} \put( 4.5, 4.5){\makebox(0,0)[l]{photon}} \put( 1.0, 2.0){\circle*{0.2}} \put( 0.5, 2.0){\makebox(0,0)[r]{$(x_{0},y_{0},z_{0},t_{0})$}} \put( 1.0, 8.0){\circle*{0.2}} \put( 0.5, 8.0){\makebox(0,0)[r]{$(x_{0},y_{0},z_{j})$}} \put( 7.0, 8.0){\circle*{0.2}} \put( 7.0, 9.0){\makebox(0,0)[c]{PMT}} \put( 7.5, 8.0){\makebox(0,0)[l]{$(x_{j},y_{j},z_{j},t_{j})$}} \put( 1.1, 2.1){ \put(0.0,1.00){\vector(-1,0){0.1}} \qbezier(0.0,1.0)(0.25,1.0)(0.5,0.75) \put(0.5,0.75){\vector(1,-1){0.1}} \put(0.4,1.5){\makebox(0,0){$\theta_{c}$}} } \end{picture} \end{center}$

$ct_j = ct_0 + (z_j - z_0) + \tan(\theta_{c}) \sqrt((x_j - x_0)^2 + (y_j - y_0)^2)$

where:

$\begin{eqnarray*} x_0 & = & \textrm{x position of track (fit parameter)} \\ y_0 & = & \textrm{y position of track (fit parameter)} \\ z_0 & = & \textrm{z position of track} \\ t_0 & = & \textrm{time of track at } z = z_0 \textrm{ (fit parameter)} \\ \\ c & = & \textrm{speed of light in vacuum} \\ \theta_{c} & = & \textrm{Cherenkov angle} \\ \end{eqnarray*}$

Defining:

$\begin{eqnarray*} t_j' & \equiv & ct_j / \tan(\theta_{c}) - (z_j - z_0)/\tan(\theta_{c}) \\ t_0' & \equiv & ct_0 / \tan(\theta_{c}) \\ \end{eqnarray*}$

$\Rightarrow (t_j' - t_0')^2 = (x_j - x_0)^2 + (y_j - y_0)^2$

The parameters $\{x_0, y_0, t_0\}$ are estimated in the constructor of this class based on consecutive pairs of equations by which the quadratic terms in $ x_0 $ , $ y_0 $ and $ t_0 $ are eliminated.

Definition at line 89 of file JLine1ZEstimator.hh.

Member Typedef Documentation

typedef double JLine1Z::* JFIT::JLine1Z::parameter_type

inherited

Definition at line 178 of file JLine1Z.hh.

Constructor & Destructor Documentation

JFIT::JEstimator< JLine1Z >::JEstimator ( )

inline

Default constructor.

Definition at line 96 of file JLine1ZEstimator.hh.

                  :
       JLine1Z()
     {}

template<class T >

JFIT::JEstimator< JLine1Z >::JEstimator	(	T	__begin,
		T	__end
	)

inline

Constructor.

Parameters

__begin	begin of data
__end	end of data

Definition at line 108 of file JLine1ZEstimator.hh.

                                    :
       JLine1Z()
     {
       (*this)(__begin, __end);
     }

Member Function Documentation

template<class T >

const JEstimator<JLine1Z>& JFIT::JEstimator< JLine1Z >::operator()	(	T	__begin,
		T	__end
	)

inline

Fit.

The template argument T refers to an iterator of a data structure which should have the following member methods:

double getX(); // [m]
double getY(); // [m]
double getZ(); // [m]
double getT(); // [ns]

Parameters

__begin	begin of data
__end	end of data

Definition at line 128 of file JLine1ZEstimator.hh.

     {
       using namespace std;
       using namespace JPP;
 
       const int N = distance(__begin, __end);
 
       if (N >= NUMBER_OF_PARAMETERS) {
 
         __x = 0.0;
         __y = 0.0;
         __z = 0.0;
         __t = 0.0;
 
         double t0 = 0.0;
 
         for (T i = __begin; i != __end; ++i) {
           __x += i->getX();
           __y += i->getY();
           __z += i->getZ();
           t0  += i->getT();
         }
         
         const double W = 1.0/N;
 
         __x *= W;
         __y *= W;
         __z *= W;
         t0  *= W;
 
         V.reset();
 
         t0 *= getSpeedOfLight();
 
         double y0 = 0.0;
         double y1 = 0.0;
         double y2 = 0.0;
 
         T j = __begin;
 
         double xi =  j->getX() - getX();
         double yi =  j->getY() - getY();
         double ti = (j->getT() * getSpeedOfLight()  -  t0  -  j->getZ() + getZ()) / getKappaC();
 
         for (bool done = false; !done; ) {
 
           if ((done = (++j == __end))) {
             j = __begin;
           }
 
           double xj =  j->getX() - getX();
           double yj =  j->getY() - getY();
           double tj = (j->getT() * getSpeedOfLight()  -  t0  -  j->getZ() + getZ()) / getKappaC();
             
           double dx = xj - xi;
           double dy = yj - yi;
           double dt = ti - tj;    // opposite sign!
             
           const double y  = ((xj + xi) * dx +
                              (yj + yi) * dy +
                              (tj + ti) * dt);
 
           dx *= 2;
           dy *= 2;
           dt *= 2;
             
           V.a00 += dx * dx;
           V.a01 += dx * dy;
           V.a02 += dx * dt;
           V.a11 += dy * dy;
           V.a12 += dy * dt;
           V.a22 += dt * dt;
 
           y0 += dx * y;
           y1 += dy * y;
           y2 += dt * y;
 
           xi = xj;
           yi = yj;
           ti = tj;
         }
 
         t0 *= getInverseSpeedOfLight();
 
         V.a10 = V.a01;
         V.a20 = V.a02;
         V.a21 = V.a12;
 
         svd.decompose(V);
 
         if (fabs(svd.S.a11) < MINIMAL_SVD_WEIGHT * fabs(svd.S.a00)) {
           THROW(JValueOutOfRange, "JEstimator<JLine1Z>::JEstimator(): singular value " << svd.S.a11 << ' ' << svd.S.a00);
         }
 
         V = svd.invert(MINIMAL_SVD_WEIGHT);
         
         __x += V.a00 * y0  +  V.a01 * y1  +  V.a02 * y2;
         __y += V.a10 * y0  +  V.a11 * y1  +  V.a12 * y2;
         __t  = V.a20 * y0  +  V.a21 * y1  +  V.a22 * y2;
         
         __t *= getKappaC() * getInverseSpeedOfLight();
         __t += t0;
 
       } else {
         throw JValueOutOfRange("JEstimator<JLine1Z>::JEstimator(): Not enough data points.");
       }
 
       return *this;
     }

template<class T >

void JFIT::JEstimator< JLine1Z >::update	(	T	__begin,
		T	__end,
		const JMatrixNZ &	A
	)

inline

Update track parameters using updated co-variance matrix (e.g. matrix with one hit switched off).

In this, it is assumed that the changes in x and y positions are small compared to the actual distances between the track and the hits.

The template argument T refers to an iterator of a data structure which should have the following member methods:

double getX(); // [m]
double getY(); // [m]
double getZ(); // [m]
double getT(); // [ns]

Parameters

__begin	begin of data
__end	end of data
A	co-variance matrix of hits

Definition at line 256 of file JLine1ZEstimator.hh.

     {
       using namespace std;
       using namespace JPP;
 
       const int N = distance(__begin, __end);
 
       if (N != (int) A.size()) {
         THROW(JValueOutOfRange, "JEstimator<JLine1Z>::update(): Wrong number of points " << N << " != " << A.size());
       }
 
       if (N >= NUMBER_OF_PARAMETERS) {
 
         double x1 = 0.0;
         double y1 = 0.0;
         double t1 = 0.0;
 
         V.reset();
         
         T i = __begin;
         
         for (size_t row = 0; row != A.size(); ++row, ++i) {
           
           const double dx = i->getX() - getX();
           const double dy = i->getY() - getY();
           
           const double rt = sqrt(dx*dx + dy*dy);
           
           double xr = getKappaC() * getInverseSpeedOfLight();
           double yr = getKappaC() * getInverseSpeedOfLight();
           double tr = 1.0;
           
           if (rt != 0.0) {
             xr *= -dx / rt;
             yr *= -dy / rt;
           }
 
           T j = __begin;
 
           for (size_t col = 0; col != A.size(); ++col, ++j) {
             
             const double dx = j->getX() - getX();
             const double dy = j->getY() - getY();
             const double dz = j->getZ() - getZ();
             
             const double rt = sqrt(dx*dx + dy*dy);
             
             double xc = getKappaC() * getInverseSpeedOfLight();
             double yc = getKappaC() * getInverseSpeedOfLight();
             double tc = 1.0;
 
             if (rt != 0.0) {
               xc *= -dx / rt;
               yc *= -dy / rt;
             }
             
             const double ts = j->getT() - (dz + rt * getKappaC()) * getInverseSpeedOfLight();
             
             const double vs = A(row,col);
 
             x1 += xr * vs * ts;
             y1 += yr * vs * ts;
             t1 += tr * vs * ts;
                     
             V.a00 += xr * vs * xc;
             V.a01 += xr * vs * yc;
             V.a02 += xr * vs * tc;
             V.a11 += yr * vs * yc;
             V.a12 += yr * vs * tc;
             V.a22 += tr * vs * tc;
           }
         }
         
         V.a10 = V.a01;
         V.a20 = V.a02;
         V.a21 = V.a12;
         
         svd.decompose(V);
 
         if (fabs(svd.S.a11) < MINIMAL_SVD_WEIGHT * fabs(svd.S.a00)) {
           THROW(JValueOutOfRange, "JEstimator<JLine1Z>::update(): singular value " << svd.S.a11 << ' ' << svd.S.a00);
         }
 
         V = svd.invert(MINIMAL_SVD_WEIGHT);
         
         __x += V.a00 * x1  +  V.a01 * y1  +  V.a02 * t1;
         __y += V.a10 * x1  +  V.a11 * y1  +  V.a12 * t1;
         __t  = V.a20 * x1  +  V.a21 * y1  +  V.a22 * t1;
 
       } else {
         THROW(JValueOutOfRange, "JEstimator<JLine1Z>::update(): Not enough data points " << N);
       }
     }

void JFIT::JLine1Z::move	(	const double	step,
		const double	velocity
	)

inlineinherited

Move vertex along this line with given velocity.

Parameters

step	step
velocity	velocity

Definition at line 62 of file JLine1Z.hh.

     {
       __z += step;
       __t += step / velocity;
     }

void JFIT::JLine1Z::setZ	(	const double	z,
		const double	velocity
	)

inlineinherited

Set z-position of vertex.

Parameters

z	z-position
velocity	velocity

Definition at line 75 of file JLine1Z.hh.

     {
       move(z - getZ(), velocity);
     }

double JFIT::JLine1Z::getDistanceSquared ( const JVector3D & pos ) const

inlineinherited

Get distance squared.

Parameters

pos position

Returns: square of distance

Definition at line 87 of file JLine1Z.hh.

     {
       const double dx = pos.getX() - this->getX();
       const double dy = pos.getY() - this->getY();
 
       return dx*dx + dy*dy;
     }

double JFIT::JLine1Z::getDistance ( const JVector3D & pos ) const

inlineinherited

Get distance.

Parameters

pos position

Returns: distance

Definition at line 102 of file JLine1Z.hh.

     {
       return sqrt(this->getDistanceSquared(pos));
     }

double JFIT::JLine1Z::getT ( const JVector3D & pos ) const

inlineinherited

Get arrival time of Cherenkov light at given position.

Parameters

pos	position [m]

Returns: time [ns]

Definition at line 114 of file JLine1Z.hh.

     {
       using namespace JPP;
 
       const double dx = pos.getX() - this->getX();
       const double dy = pos.getY() - this->getY();
       const double dz = pos.getZ() - this->getZ();
 
       const double R  = sqrt(dx*dx + dy*dy);
       
       return this->getT() + (dz + R * getKappaC()) * getInverseSpeedOfLight();
     }

double JGEOMETRY3D::JTime::getT ( ) const

inlineinherited

Get time.

Returns: time

Definition at line 122 of file JGeometry3D/JTime.hh.

     {
       return __t;
     }

double JFIT::JLine1Z::getZ ( const JPosition3D & pos ) const

inlineinherited

Get point of emission of Cherenkov light along muon path.

Parameters

pos position

Returns: position along muon path

Definition at line 134 of file JLine1Z.hh.

     {
       using namespace JPP;
 
       return pos.getZ()  -  this->getDistance(pos) / getTanThetaC();
     }

double JGEOMETRY3D::JVector3D::getZ ( ) const

inlineinherited

Get z position.

Returns: z position

Definition at line 115 of file JVector3D.hh.

     { 
       return __z;
     }

JVersor3D JFIT::JLine1Z::getDirection ( const JVector3D & pos ) const

inlineinherited

Get photon direction of Cherenkov light on PMT.

Parameters

pos	PMT position

Returns: direction

Definition at line 148 of file JLine1Z.hh.

     {
       using namespace JPP;
 
       double dx = pos.getX() - this->getX();
       double dy = pos.getY() - this->getY();
 
       const double R = sqrt(dx*dx + dy*dy);
 
       dx *= getSinThetaC() / R;
       dy *= getSinThetaC() / R;
 
       const double dz = getCosThetaC();
       
       return JVersor3D(dx,dy,dz);
     }

double JFIT::JLine1Z::getDot ( const JAxis3D & axis ) const

inlineinherited

Get cosine angle of impact of Cherenkov light on PMT.

Parameters

axis PMT axis

Returns: cosine angle of impact

Definition at line 172 of file JLine1Z.hh.

     {
       return getDirection(axis.getPosition()).getDot(axis.getDirection());
     }

double JGEOMETRY3D::JPosition3D::getDot ( const JAngle3D & angle ) const

inlineinherited

Get dot product.

Parameters

angle angle

Returns: dot product

Definition at line 378 of file JPosition3D.hh.

     {
       return
         getX() * angle.getDX() +
         getY() * angle.getDY() +
         getZ() * angle.getDZ();
     }

double JGEOMETRY3D::JPosition3D::getDot ( const JVersor3D & dir ) const

inlineinherited

Get dot product.

Parameters

dir direction

Returns: dot product

Definition at line 393 of file JPosition3D.hh.

     {
       return
         getX() * dir.getDX() +
         getY() * dir.getDY() +
         getZ() * dir.getDZ();
     }

double JGEOMETRY3D::JPosition3D::getDot ( const JVersor3Z & dir ) const

inlineinherited

Get dot product.

Parameters

dir direction

Returns: dot product

Definition at line 408 of file JPosition3D.hh.

     {
       return
         getX() * dir.getDX() +
         getY() * dir.getDY() +
         getZ() * dir.getDZ();
     }

double JGEOMETRY3D::JVector3D::getDot ( const JVector3D & vector ) const

inlineinherited

Get dot product.

Parameters

vector vector

Returns: dot product

Definition at line 282 of file JVector3D.hh.

     {
       return
         getX() * vector.getX() +
         getY() * vector.getY() +
         getZ() * vector.getZ();
     }

static parameter_type JFIT::JLine1Z::pX ( )

inlinestaticinherited

Definition at line 180 of file JLine1Z.hh.

180 { return &JLine1Z::__x; }

JGEOMETRY3D::JVector3D::__x

double __x

Definition: JVector3D.hh:310

static parameter_type JFIT::JLine1Z::pY ( )

inlinestaticinherited

Definition at line 181 of file JLine1Z.hh.

181 { return &JLine1Z::__y; }

JGEOMETRY3D::JVector3D::__y

double __y

Definition: JVector3D.hh:311

static parameter_type JFIT::JLine1Z::pT ( )

inlinestaticinherited

Definition at line 182 of file JLine1Z.hh.

182 { return &JLine1Z::__t; }

JGEOMETRY3D::JTime::__t

double __t

Definition: JGeometry3D/JTime.hh:191

JVertex3D& JGEOMETRY3D::JVertex3D::negate ( )

inlineinherited

Prefix unary minus.

Returns: line

Definition at line 73 of file JVertex3D.hh.

     {
       JPosition3D::negate();
       JTime      ::negate();
 
       return *this;
     }

JVertex3D& JGEOMETRY3D::JVertex3D::add ( const JVertex3D & value )

inlineinherited

Addition operator.

Parameters

value line

Returns: line

Definition at line 87 of file JVertex3D.hh.

     {
       JPosition3D::add(value);
       JTime      ::add(value);
 
       return *this;
     }

JVector3D& JGEOMETRY3D::JVector3D::add ( const JVector3D & vector )

inlineinherited

Add vector.

Parameters

vector vector

Returns: this vector

Definition at line 142 of file JVector3D.hh.

     {
       __x += vector.getX();
       __y += vector.getY();
       __z += vector.getZ();
 
       return *this;
     }

JTime& JGEOMETRY3D::JTime::add ( const JTime & value )

inlineinherited

Addition operator.

Parameters

value time

Returns: time

Definition at line 67 of file JGeometry3D/JTime.hh.

     {
       __t += value.getT();
 
       return *this;
     }

JVertex3D& JGEOMETRY3D::JVertex3D::sub ( const JVertex3D & value )

inlineinherited

Subtraction operator.

Parameters

value line

Returns: line

Definition at line 102 of file JVertex3D.hh.

     {
       JPosition3D::sub(value);
       JTime      ::sub(value);
 
       return *this;
     }

JVector3D& JGEOMETRY3D::JVector3D::sub ( const JVector3D & vector )

inlineinherited

Subtract vector.

Parameters

vector vector

Returns: this vector

Definition at line 158 of file JVector3D.hh.

     {
       __x -= vector.getX();
       __y -= vector.getY();
       __z -= vector.getZ();
 
       return *this;
     }

JTime& JGEOMETRY3D::JTime::sub ( const JTime & value )

inlineinherited

Subtraction operator.

Parameters

value time

Returns: time

Definition at line 81 of file JGeometry3D/JTime.hh.

     {
       __t -= value.getT();
 
       return *this;
     }

JVertex3D& JGEOMETRY3D::JVertex3D::mul ( const double value )

inlineinherited

Multiplication operator.

Parameters

value multiplication factor

Returns: line

Definition at line 117 of file JVertex3D.hh.

     {
       JPosition3D::mul(value);
       JTime      ::mul(value);
 
       return *this;
     }

template<class JFirst_t, class JSecond_t = JNullType>

JFirst_t& JMATH::JMath< JFirst_t, JSecond_t >::mul ( const JSecond_t & object )

inlineinherited

Multiply with object.

Parameters

object object

Returns: result object

Definition at line 357 of file JMath.hh.

     {
       return static_cast<JFirst_t&>(*this) = JCalculator<JFirst_t>::calculator.mul(static_cast<const JFirst_t&>(*this), object);
     }

template<class JFirst_t, class JSecond_t = JNullType>

JFirst_t& JMATH::JMath< JFirst_t, JSecond_t >::mul ( const JSecond_t & object )

inlineinherited

Multiply with object.

Parameters

object object

Returns: result object

Definition at line 357 of file JMath.hh.

     {
       return static_cast<JFirst_t&>(*this) = JCalculator<JFirst_t>::calculator.mul(static_cast<const JFirst_t&>(*this), object);
     }

JVertex3D& JGEOMETRY3D::JVertex3D::div ( const double value )

inlineinherited

Division operator.

Parameters

value multiplication factor

Returns: line

Definition at line 132 of file JVertex3D.hh.

     {
       JPosition3D::div(value);
       JTime      ::div(value);
 
       return *this;
     }

const JPosition3D& JGEOMETRY3D::JPosition3D::getPosition ( ) const

inlineinherited

Get position.

Returns: position

Definition at line 130 of file JPosition3D.hh.

     {
       return static_cast<const JPosition3D&>(*this);
     }

JPosition3D& JGEOMETRY3D::JPosition3D::getPosition ( )

inlineinherited

Get position.

Returns: position

Definition at line 141 of file JPosition3D.hh.

     {
       return static_cast<JPosition3D&>(*this);
     }

void JGEOMETRY3D::JPosition3D::setPosition ( const JVector3D & pos )

inlineinherited

Set position.

Parameters

pos position

Definition at line 152 of file JPosition3D.hh.

     {
       static_cast<JVector3D&>(*this) = pos;
     }

JGEOMETRY3D::JPosition3D::operator JAngle3D ( ) const

inlineinherited

Type conversion operator.

Returns: angle

Definition at line 163 of file JPosition3D.hh.

     {
       return JAngle3D(getX(), getY(), getZ());
     }

JGEOMETRY3D::JPosition3D::operator JVersor3D ( ) const

inlineinherited

Type conversion operator.

Returns: direction

Definition at line 174 of file JPosition3D.hh.

     {
       return JVersor3D(getX(), getY(), getZ());
     }

JPosition3D& JGEOMETRY3D::JPosition3D::rotate ( const JRotation3D & R )

inlineinherited

Rotate.

Parameters

R	rotation matrix

Returns: this position

Definition at line 186 of file JPosition3D.hh.

     {
       R.rotate(__x, __y, __z);
 
       return *this;
     }

JPosition3D& JGEOMETRY3D::JPosition3D::rotate ( const JRotation3X & R )

inlineinherited

Rotate around X-axis.

Parameters

R	rotation matrix

Returns: this position

Definition at line 214 of file JPosition3D.hh.

     {
       R.rotate(__y, __z);
 
       return *this;
     }

JPosition3D& JGEOMETRY3D::JPosition3D::rotate ( const JRotation3Y & R )

inlineinherited

Rotate around Y-axis.

Parameters

R	rotation matrix

Returns: this position

Definition at line 242 of file JPosition3D.hh.

     {
       R.rotate(__x, __z);
 
       return *this;
     }

JPosition3D& JGEOMETRY3D::JPosition3D::rotate ( const JRotation3Z & R )

inlineinherited

Rotate around Z-axis.

Parameters

R	rotation matrix

Returns: this position

Definition at line 270 of file JPosition3D.hh.

     {
       R.rotate(__x, __y);
 
       return *this;
     }

JPosition3D& JGEOMETRY3D::JPosition3D::rotate ( const JQuaternion3D & Q )

inlineinherited

Rotate.

Parameters

Q	quaternion

Returns: this position

Definition at line 298 of file JPosition3D.hh.

     {
       Q.rotate(__x, __y, __z);
       
       return *this;
     }

JPosition3D& JGEOMETRY3D::JPosition3D::rotate_back ( const JRotation3D & R )

inlineinherited

Rotate back.

Parameters

R	rotation matrix

Returns: this position

Definition at line 200 of file JPosition3D.hh.

     {
       R.rotate_back(__x, __y, __z);
 
       return *this;
     }

JPosition3D& JGEOMETRY3D::JPosition3D::rotate_back ( const JRotation3X & R )

inlineinherited

Rotate back around X-axis.

Parameters

R	rotation matrix

Returns: this position

Definition at line 228 of file JPosition3D.hh.

     {
       R.rotate_back(__y, __z);
 
       return *this;
     }

JPosition3D& JGEOMETRY3D::JPosition3D::rotate_back ( const JRotation3Y & R )

inlineinherited

Rotate back around Y-axis.

Parameters

R	rotation matrix

Returns: this position

Definition at line 256 of file JPosition3D.hh.

     {
       R.rotate_back(__x, __z);
 
       return *this;
     }

JPosition3D& JGEOMETRY3D::JPosition3D::rotate_back ( const JRotation3Z & R )

inlineinherited

Rotate back around Z-axis.

Parameters

R	rotation matrix

Returns: this position

Definition at line 284 of file JPosition3D.hh.

     {
       R.rotate_back(__x, __y);
 
       return *this;
     }

JPosition3D& JGEOMETRY3D::JPosition3D::rotate_back ( const JQuaternion3D & Q )

inlineinherited

Rotate back.

Parameters

Q	quaternion

Returns: this position

Definition at line 312 of file JPosition3D.hh.

     {
       Q.rotate_back(__x, __y, __z);
 
       return *this;
     }

void JGEOMETRY3D::JPosition3D::transform	(	const JRotation3D &	R,
		const JVector3D &	pos
	)

inlineinherited

Transform position.

The final position is obtained as follows:

rotation of the position according matrix R;
offset position with pos;
rotation of position around z-axis, such that final position lies in x-z plane;

Parameters

R	rotation matrix
pos	position of origin (after rotation)

Definition at line 331 of file JPosition3D.hh.

     {
       // rotate geometry to system with particle direction along z-axis
 
       rotate(R);
 
       // offset with respect to origin
 
       sub(pos);
 
       // rotate geometry to x-z plane
 
       __x = sqrt(__x*__x + __y*__y);
       __y = 0.0;
     }

JVector3D& JGEOMETRY3D::JVector3D::transform ( const JMatrix3D & T )

inlineinherited

Transform.

Parameters

T matrix

Returns: this vector

Definition at line 206 of file JVector3D.hh.

     {
       T.transform(__x, __y, __z);
 
       return *this;
     }

void JGEOMETRY3D::JPosition3D::transform_back	(	const JRotation3D &	R,
		const JVector3D &	pos
	)

inlineinherited

Transform back position.

The final position is obtained as follows:

offset position with position pos;
rotation of postion according matrix R;

Parameters

R	rotation matrix
pos	position of origin (before rotation)

Definition at line 359 of file JPosition3D.hh.

     {
       // offset with respect to origin
 
       add(pos);
 
       // rotate back geometry to system with particle direction along z-axis
 
       rotate_back(R);
     }

JGEOMETRY3D::JVector3D::operator JVector2D ( ) const

inlineinherited

Type conversion operator.

Returns: JVector2D

Definition at line 83 of file JVector3D.hh.

     {
       return JVector2D(this->getX(), this->getY());
     }

double JGEOMETRY3D::JVector3D::getX ( ) const

inlineinherited

Get x position.

Returns: x position

Definition at line 94 of file JVector3D.hh.

     { 
       return __x;
     }

double JGEOMETRY3D::JVector3D::getY ( ) const

inlineinherited

Get y position.

Returns: y position

Definition at line 104 of file JVector3D.hh.

     { 
       return __y;
     }

bool JGEOMETRY3D::JVector3D::equals	(	const JVector3D &	vector,
		const double	precision = `std::numeric_limits<double>::min()`
	)		const

inlineinherited

Check equality.

Parameters

vector	vector
precision	precision

Returns: true if vectors are equal; else false

Definition at line 221 of file JVector3D.hh.

     {
       return (fabs(getX() - vector.getX()) <= precision &&
               fabs(getY() - vector.getY()) <= precision &&
               fabs(getZ() - vector.getZ()) <= precision);
     }

double JGEOMETRY3D::JVector3D::getLengthSquared ( ) const

inlineinherited

Get length squared.

Returns: square of length

Definition at line 235 of file JVector3D.hh.

     { 
       return getX()*getX() + getY()*getY() + getZ()*getZ();
     }

double JGEOMETRY3D::JVector3D::getLength ( ) const

inlineinherited

Get length.

Returns: length

Definition at line 246 of file JVector3D.hh.

     { 
       return sqrt(getLengthSquared());
     }

JVector3D& JGEOMETRY3D::JVector3D::getCross	(	const JVector3D &	first,
		const JVector3D &	second
	)

inlineinherited

Get cross product.

Note that this vector should not overlap with the first or second vector,

Parameters

first	first vector
second	second vector

Returns: this vector

Definition at line 299 of file JVector3D.hh.

     {
       __x = first .getY() * second.getZ()  -  second.getY() * first .getZ();
       __y = second.getX() * first .getZ()  -  first .getX() * second.getZ();
       __z = first .getX() * second.getY()  -  second.getX() * first .getY();
 
       return *this;
     }

Member Data Documentation

const int JFIT::JEstimator< JLine1Z >::NUMBER_OF_PARAMETERS = 3

static

number of parameters of fit

Definition at line 351 of file JLine1ZEstimator.hh.

JMATH::JMatrix3S JFIT::JEstimator< JLine1Z >::V

co-variance matrix of fit parameters

Definition at line 352 of file JLine1ZEstimator.hh.

JMATH::JSVD3D JFIT::JEstimator< JLine1Z >::svd

Definition at line 353 of file JLine1ZEstimator.hh.

double JFIT::JEstimator< JLine1Z >::MINIMAL_SVD_WEIGHT = 1.0e-4

static

minimal SVD weight.

Set default minimal SVD weight.

Definition at line 354 of file JLine1ZEstimator.hh.

double JGEOMETRY3D::JVector3D::__x

protectedinherited

Definition at line 310 of file JVector3D.hh.

double JGEOMETRY3D::JVector3D::__y

protectedinherited

Definition at line 311 of file JVector3D.hh.

double JGEOMETRY3D::JVector3D::__z

protectedinherited

Definition at line 312 of file JVector3D.hh.

double JGEOMETRY3D::JTime::__t

protectedinherited

Definition at line 191 of file JGeometry3D/JTime.hh.

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

JLine1ZEstimator.hh

Public Types

Public Member Functions

Static Public Member Functions

Public Attributes

Static Public Attributes

Protected Attributes

Detailed Description

template<> class JFIT::JEstimator< JLine1Z >

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation

template<>
class JFIT::JEstimator< JLine1Z >