Jpp  18.1.0
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Public Types | Public Member Functions | Static Public Member Functions | Public Attributes | Static Public Attributes | Protected Attributes | List of all members
JFIT::JEstimator< JPoint4D > Class Template Reference

Linear fit of bright point (position and time) between hits (objects with position and time). More...

#include <JPoint4DEstimator.hh>

Inheritance diagram for JFIT::JEstimator< JPoint4D >:
JFIT::JPoint4D JGEOMETRY3D::JVertex3D JGEOMETRY3D::JPosition3D JGEOMETRY3D::JTime JMATH::JMath< JFirst_t, JSecond_t > JGEOMETRY3D::JVector3D JMATH::JMath< JFirst_t, JSecond_t >

Public Types

typedef double JPoint4D::* parameter_type
 

Public Member Functions

 JEstimator ()
 Fit constructor. More...
 
template<class T >
 JEstimator (T __begin, T __end)
 Fit constructor. More...
 
template<class T >
void operator() (T __begin, T __end)
 Fit function. 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...
 
double getT (const JVector3D &pos) const
 Get arrival time of Cherenkov light at given position. More...
 
double getT () const
 Get time. 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...
 
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...
 
double getZ () const
 Get z 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...
 
double getDistanceSquared (const JVector3D &pos) const
 Get squared of distance to point. More...
 
double getDistance (const JVector3D &pos) const
 Get distance to point. 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 pZ ()
 
static parameter_type pT ()
 

Public Attributes

JMATH::JMatrix4S V
 co-variance matrix of fit parameters More...
 

Static Public Attributes

static const int NUMBER_OF_PARAMETERS = 4
 number of parameters of fit More...
 

Protected Attributes

double __x
 
double __y
 
double __z
 
double __t
 

Detailed Description

template<>
class JFIT::JEstimator< JPoint4D >

Linear fit of bright point (position and time) between hits (objects with position and time).

\begin{center}\setlength{\unitlength}{0.6cm}\begin{picture}(12,7) \put( 6.0, 1.0){\circle*{0.3}} \put( 6.0, 0.0){\makebox(0,0)[b]{$(x_{0},y_{0},z_{0})$}} \multiput(6.0, 1.0)(-0.5, 0.5){10}{\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( 1.0, 6.0){\circle*{0.2}} \put( 1.0, 6.5){\makebox(0,0)[b]{$(x_i,y_i,z_i,t_i)$}} \multiput(6.0, 1.0)( 0.5, 0.5){10}{\qbezier(0.0,0.0)(0.35,-0.1)( 0.25,0.25)\qbezier( 0.25,0.25)(0.15, 0.6)( 0.5,0.5)} \put(11.0, 6.0){\circle*{0.2}} \put(11.0, 6.5){\makebox(0,0)[b]{$(x_j,y_j,z_j,t_j)$}} \end{picture} \end{center}

\[ t_j = t_0 + \frac{c}{n} \times \sqrt((x_j - x_0)^2 + (y_j - y_0)^2 + (z_j - z_0)^2) \]

where:

\begin{eqnarray*} x_0 & = & \textrm{x position of vertex (fit parameter)} \\ y_0 & = & \textrm{y position of vertex (fit parameter)} \\ z_0 & = & \textrm{z position of vertex (fit parameter)} \\ t_0 & = & \textrm{time at vertex (fit parameter)} \\ \\ c & = & \textrm{speed of light (in vacuum)} \\ n & = & \textrm{index of refraction corresponding to the group velocity of light} \\ \end{eqnarray*}

Defining:

\begin{eqnarray*} t_j' & \equiv & nct_j \\ t_0' & \equiv & nct_0 \\ \end{eqnarray*}

\[ \Rightarrow (t_j' - t_0')^2 = (x_j - x_0)^2 + (y_j - y_0)^2 + (z_j - z_0)^2 \]

The parameters $ \{x_0, y_0, z_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 $, $ z_0 $ and $ t_0 $ are eliminated.

Definition at line 69 of file JPoint4DEstimator.hh.

Member Typedef Documentation

typedef double JPoint4D::* JFIT::JPoint4D::parameter_type
inherited

Definition at line 56 of file JPoint4D.hh.

Constructor & Destructor Documentation

JFIT::JEstimator< JPoint4D >::JEstimator ( )
inline

Fit constructor.

Definition at line 77 of file JPoint4DEstimator.hh.

77  :
78  JPoint4D()
79  {}
JFIT::JPoint4D::JPoint4D
JPoint4D()
Default constructor.
Definition: JPoint4D.hh:29
template<class T >
JFIT::JEstimator< JPoint4D >::JEstimator ( T  __begin,
T  __end 
)
inline

Fit constructor.

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
__beginbegin of data
__endend of data

Definition at line 95 of file JPoint4DEstimator.hh.

95  :
96  JPoint4D()
97  {
98  (*this)(__begin, __end);
99  }
JFIT::JPoint4D::JPoint4D
JPoint4D()
Default constructor.
Definition: JPoint4D.hh:29

Member Function Documentation

template<class T >
void JFIT::JEstimator< JPoint4D >::operator() ( T  __begin,
T  __end 
)
inline

Fit function.

This method is used to find the vertex of a given set of hits

Parameters
__beginbegin of data
__endend of data

Definition at line 110 of file JPoint4DEstimator.hh.

111  {
112  using namespace std;
113  using namespace JPP;
114 
115  const int N = distance(__begin, __end);
116 
117  if (N >= NUMBER_OF_PARAMETERS) {
118 
119  double t0 = 0.0;
120 
121  __x = 0.0;
122  __y = 0.0;
123  __z = 0.0;
124 
125  for (T i = __begin; i != __end; ++i) {
126  __x += i->getX();
127  __y += i->getY();
128  __z += i->getZ();
129  t0 += i->getT();
130  }
131 
132  div(N);
133  t0 /= N;
134 
135  V.reset();
136 
137  t0 *= getSpeedOfLight();
138 
139  double y0 = 0.0;
140  double y1 = 0.0;
141  double y2 = 0.0;
142  double y3 = 0.0;
143 
144  T j = __begin;
145 
146  double xi = j->getX() - getX();
147  double yi = j->getY() - getY();
148  double zi = j->getZ() - getZ();
149  double ti = (j->getT() * getSpeedOfLight() - t0) / getIndexOfRefraction();
150 
151  for (bool done = false; !done; ) {
152 
153  if ((done = (++j == __end))) {
154  j = __begin;
155  }
156 
157  double xj = j->getX() - getX();
158  double yj = j->getY() - getY();
159  double zj = j->getZ() - getZ();
160  double tj = (j->getT() * getSpeedOfLight() - t0) / getIndexOfRefraction();
161 
162  double dx = xj - xi;
163  double dy = yj - yi;
164  double dz = zj - zi;
165  double dt = ti - tj; // opposite sign!
166 
167  const double y = ((xj + xi) * dx +
168  (yj + yi) * dy +
169  (zj + zi) * dz +
170  (tj + ti) * dt);
171 
172  dx *= 2;
173  dy *= 2;
174  dz *= 2;
175  dt *= 2;
176 
177  V.a00 += dx * dx;
178  V.a01 += dx * dy;
179  V.a02 += dx * dz;
180  V.a03 += dx * dt;
181  V.a11 += dy * dy;
182  V.a12 += dy * dz;
183  V.a13 += dy * dt;
184  V.a22 += dz * dz;
185  V.a23 += dz * dt;
186  V.a33 += dt * dt;
187 
188  y0 += dx * y;
189  y1 += dy * y;
190  y2 += dz * y;
191  y3 += dt * y;
192 
193  xi = xj;
194  yi = yj;
195  zi = zj;
196  ti = tj;
197  }
198 
199  t0 *= getInverseSpeedOfLight();
200 
201  V.a10 = V.a01;
202  V.a20 = V.a02;
203  V.a30 = V.a03;
204  V.a21 = V.a12;
205  V.a31 = V.a13;
206  V.a32 = V.a23;
207 
208  V.invert();
209 
210  __x += V.a00 * y0 + V.a01 * y1 + V.a02 * y2 + V.a03 * y3;
211  __y += V.a10 * y0 + V.a11 * y1 + V.a12 * y2 + V.a13 * y3;
212  __z += V.a20 * y0 + V.a21 * y1 + V.a22 * y2 + V.a23 * y3;
213  __t = V.a30 * y0 + V.a31 * y1 + V.a32 * y2 + V.a33 * y3;
214 
215  __t *= getInverseSpeedOfLight() * getIndexOfRefraction();
216  __t += t0;
217 
218  } else {
219  throw JValueOutOfRange("JEstimator<JPoint4D>::operator(): Not enough data points.");
220  }
221  }
JFIT::JEstimator< JPoint4D >::NUMBER_OF_PARAMETERS
static const int NUMBER_OF_PARAMETERS
number of parameters of fit
Definition: JPoint4DEstimator.hh:223
JFIT::JEstimator< JPoint4D >::V
JMATH::JMatrix4S V
co-variance matrix of fit parameters
Definition: JPoint4DEstimator.hh:224
JMATH::JMatrix4S::invert
void invert()
Invert matrix.
Definition: JMatrix4S.hh:77
distance
std::vector< T >::difference_type distance(typename std::vector< T >::const_iterator first, typename PhysicsEvent::const_iterator< T > second)
Specialisation of STL distance.
Definition: PhysicsEvent.hh:434
JPHYSICS::getIndexOfRefraction
double getIndexOfRefraction()
Get average index of refraction of water corresponding to group velocity.
Definition: JPhysics/JConstants.hh:134
JMATH::JMatrix4D::reset
JMatrix4D & reset()
Set matrix to the null matrix.
Definition: JMatrix4D.hh:142
JGEOMETRY3D::JVertex3D::div
JVertex3D & div(const double value)
Division operator.
Definition: JVertex3D.hh:132
T
do set_variable OUTPUT_DIRECTORY $WORKDIR T
Definition: JCalibrateHeight.sh:61
JGEOMETRY3D::JVector3D::getY
double getY() const
Get y position.
Definition: JVector3D.hh:104
N
then usage $script< input file >[option[primary[working directory]]] nWhere option can be N
Definition: JMuonPostfit.sh:40
JPHYSICS::getSpeedOfLight
const double getSpeedOfLight()
Get speed of light.
Definition: JPhysics/JConstants.hh:112
JPHYSICS::getInverseSpeedOfLight
const double getInverseSpeedOfLight()
Get inverse speed of light.
Definition: JPhysics/JConstants.hh:123
JGEOMETRY3D::JVector3D::getX
double getX() const
Get x position.
Definition: JVector3D.hh:94
JTOOLS::j
int j
Definition: JPolint.hh:703
JGEOMETRY3D::JVector3D::getZ
double getZ() const
Get z position.
Definition: JVector3D.hh:115
done
esac done
Definition: JAddHDE.sh:21
static parameter_type JFIT::JPoint4D::pX ( )
inlinestaticinherited

Definition at line 58 of file JPoint4D.hh.

58 { return &JPoint4D::__x; }
static parameter_type JFIT::JPoint4D::pY ( )
inlinestaticinherited

Definition at line 59 of file JPoint4D.hh.

59 { return &JPoint4D::__y; }
static parameter_type JFIT::JPoint4D::pZ ( )
inlinestaticinherited

Definition at line 60 of file JPoint4D.hh.

60 { return &JPoint4D::__z; }
static parameter_type JFIT::JPoint4D::pT ( )
inlinestaticinherited

Definition at line 61 of file JPoint4D.hh.

61 { return &JPoint4D::__t; }
JVertex3D& JGEOMETRY3D::JVertex3D::negate ( )
inlineinherited

Prefix unary minus.

Returns
line

Definition at line 73 of file JVertex3D.hh.

74  {
75  JPosition3D::negate();
76  JTime ::negate();
77 
78  return *this;
79  }
JGEOMETRY3D::JTime::negate
JTime & negate()
Prefix unary minus.
Definition: JGeometry3D/JTime.hh:53
JGEOMETRY3D::JVector3D::negate
JVector3D & negate()
Negate vector.
Definition: JVector3D.hh:126
JVertex3D& JGEOMETRY3D::JVertex3D::add ( const JVertex3D &  value)
inlineinherited

Addition operator.

Parameters
valueline
Returns
line

Definition at line 87 of file JVertex3D.hh.

88  {
89  JPosition3D::add(value);
90  JTime ::add(value);
91 
92  return *this;
93  }
JGEOMETRY3D::JTime::add
JTime & add(const JTime &value)
Addition operator.
Definition: JGeometry3D/JTime.hh:67
JGEOMETRY3D::JVector3D::add
JVector3D & add(const JVector3D &vector)
Add vector.
Definition: JVector3D.hh:142
JVector3D& JGEOMETRY3D::JVector3D::add ( const JVector3D &  vector)
inlineinherited

Add vector.

Parameters
vectorvector
Returns
this vector

Definition at line 142 of file JVector3D.hh.

143  {
144  __x += vector.getX();
145  __y += vector.getY();
146  __z += vector.getZ();
147 
148  return *this;
149  }
JGEOMETRY3D::JVector3D::getY
double getY() const
Get y position.
Definition: JVector3D.hh:104
JGEOMETRY3D::JVector3D::getX
double getX() const
Get x position.
Definition: JVector3D.hh:94
JGEOMETRY3D::JVector3D::getZ
double getZ() const
Get z position.
Definition: JVector3D.hh:115
JTime& JGEOMETRY3D::JTime::add ( const JTime &  value)
inlineinherited

Addition operator.

Parameters
valuetime
Returns
time

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

68  {
69  __t += value.getT();
70 
71  return *this;
72  }
JGEOMETRY3D::JTime::getT
double getT() const
Get time.
Definition: JGeometry3D/JTime.hh:122
JVertex3D& JGEOMETRY3D::JVertex3D::sub ( const JVertex3D &  value)
inlineinherited

Subtraction operator.

Parameters
valueline
Returns
line

Definition at line 102 of file JVertex3D.hh.

103  {
104  JPosition3D::sub(value);
105  JTime ::sub(value);
106 
107  return *this;
108  }
JGEOMETRY3D::JTime::sub
JTime & sub(const JTime &value)
Subtraction operator.
Definition: JGeometry3D/JTime.hh:81
JGEOMETRY3D::JVector3D::sub
JVector3D & sub(const JVector3D &vector)
Subtract vector.
Definition: JVector3D.hh:158
JVector3D& JGEOMETRY3D::JVector3D::sub ( const JVector3D &  vector)
inlineinherited

Subtract vector.

Parameters
vectorvector
Returns
this vector

Definition at line 158 of file JVector3D.hh.

159  {
160  __x -= vector.getX();
161  __y -= vector.getY();
162  __z -= vector.getZ();
163 
164  return *this;
165  }
JGEOMETRY3D::JVector3D::getY
double getY() const
Get y position.
Definition: JVector3D.hh:104
JGEOMETRY3D::JVector3D::getX
double getX() const
Get x position.
Definition: JVector3D.hh:94
JGEOMETRY3D::JVector3D::getZ
double getZ() const
Get z position.
Definition: JVector3D.hh:115
JTime& JGEOMETRY3D::JTime::sub ( const JTime &  value)
inlineinherited

Subtraction operator.

Parameters
valuetime
Returns
time

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

82  {
83  __t -= value.getT();
84 
85  return *this;
86  }
JGEOMETRY3D::JTime::getT
double getT() const
Get time.
Definition: JGeometry3D/JTime.hh:122
JVertex3D& JGEOMETRY3D::JVertex3D::mul ( const double  value)
inlineinherited

Multiplication operator.

Parameters
valuemultiplication factor
Returns
line

Definition at line 117 of file JVertex3D.hh.

118  {
119  JPosition3D::mul(value);
120  JTime ::mul(value);
121 
122  return *this;
123  }
JGEOMETRY3D::JVector3D::mul
JVector3D & mul(const double factor)
Scale vector.
Definition: JVector3D.hh:174
JGEOMETRY3D::JTime::mul
JTime & mul(const double value)
Multiplication operator.
Definition: JGeometry3D/JTime.hh:95
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
objectobject
Returns
result object

Definition at line 354 of file JMath.hh.

355  {
356  return static_cast<JFirst_t&>(*this) = JFirst_t().mul(static_cast<const JFirst_t&>(*this), object);
357  }
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
objectobject
Returns
result object

Definition at line 354 of file JMath.hh.

355  {
356  return static_cast<JFirst_t&>(*this) = JFirst_t().mul(static_cast<const JFirst_t&>(*this), object);
357  }
JVertex3D& JGEOMETRY3D::JVertex3D::div ( const double  value)
inlineinherited

Division operator.

Parameters
valuemultiplication factor
Returns
line

Definition at line 132 of file JVertex3D.hh.

133  {
134  JPosition3D::div(value);
135  JTime ::div(value);
136 
137  return *this;
138  }
JGEOMETRY3D::JTime::div
JTime & div(const double value)
Division operator.
Definition: JGeometry3D/JTime.hh:109
JGEOMETRY3D::JVector3D::div
JVector3D & div(const double factor)
Scale vector.
Definition: JVector3D.hh:190
double JGEOMETRY3D::JVertex3D::getT ( const JVector3D &  pos) const
inlineinherited

Get arrival time of Cherenkov light at given position.

Parameters
posposition [m]
Returns
time [ns]

Definition at line 147 of file JVertex3D.hh.

148  {
149  using namespace JPHYSICS;
150 
151  return this->getT() + this->getDistance(pos) * getInverseSpeedOfLight() * getIndexOfRefraction();
152  }
JGEOMETRY3D::JTime::getT
double getT() const
Get time.
Definition: JGeometry3D/JTime.hh:122
JPHYSICS::getIndexOfRefraction
double getIndexOfRefraction()
Get average index of refraction of water corresponding to group velocity.
Definition: JPhysics/JConstants.hh:134
JGEOMETRY3D::JVector3D::getDistance
double getDistance(const JVector3D &pos) const
Get distance to point.
Definition: JVector3D.hh:270
JPHYSICS::getInverseSpeedOfLight
const double getInverseSpeedOfLight()
Get inverse speed of light.
Definition: JPhysics/JConstants.hh:123
double JGEOMETRY3D::JTime::getT ( ) const
inlineinherited

Get time.

Returns
time

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

123  {
124  return __t;
125  }
JVersor3D JGEOMETRY3D::JVertex3D::getDirection ( const JVector3D &  pos) const
inlineinherited

Get photon direction of Cherenkov light on PMT.

Parameters
posPMT position
Returns
direction

Definition at line 161 of file JVertex3D.hh.

162  {
163  JPosition3D D(pos);
164 
165  D.sub(this->getPosition());
166 
167  return JVersor3D(D);
168  }
JGEOMETRY3D::JPosition3D::getPosition
const JPosition3D & getPosition() const
Get position.
Definition: JPosition3D.hh:130
JGEOMETRY3D::JPosition3D
Data structure for position in three dimensions.
Definition: JPosition3D.hh:36
JGEOMETRY3D::JVersor3D
Data structure for normalised vector in three dimensions.
Definition: JVersor3D.hh:26
D
do echo Generating $dir eval D
Definition: JDrawLED.sh:53
double JGEOMETRY3D::JVertex3D::getDot ( const JAxis3D &  axis) const
inlineinherited

Get cosine angle of impact of Cherenkov light on PMT.

Parameters
axisPMT axis
Returns
cosine angle of impact

Definition at line 177 of file JVertex3D.hh.

178  {
179  return getDirection(axis.getPosition()).getDot(axis.getDirection());
180  }
JGEOMETRY3D::JDirection3D::getDirection
const JDirection3D & getDirection() const
Get direction.
Definition: JDirection3D.hh:107
JGEOMETRY3D::JVertex3D::getDot
double getDot(const JAxis3D &axis) const
Get cosine angle of impact of Cherenkov light on PMT.
Definition: JVertex3D.hh:177
JGEOMETRY3D::JPosition3D::getPosition
const JPosition3D & getPosition() const
Get position.
Definition: JPosition3D.hh:130
JGEOMETRY3D::JVertex3D::getDirection
JVersor3D getDirection(const JVector3D &pos) const
Get photon direction of Cherenkov light on PMT.
Definition: JVertex3D.hh:161
double JGEOMETRY3D::JPosition3D::getDot ( const JAngle3D &  angle) const
inlineinherited

Get dot product.

Parameters
angleangle
Returns
dot product

Definition at line 378 of file JPosition3D.hh.

379  {
380  return
381  getX() * angle.getDX() +
382  getY() * angle.getDY() +
383  getZ() * angle.getDZ();
384  }
JGEOMETRY3D::JVector3D::getY
double getY() const
Get y position.
Definition: JVector3D.hh:104
JGEOMETRY3D::JVector3D::getX
double getX() const
Get x position.
Definition: JVector3D.hh:94
JGEOMETRY3D::JVector3D::getZ
double getZ() const
Get z position.
Definition: JVector3D.hh:115
JGEOMETRY3D::JAngle3D::getDX
double getDX() const
Get x direction.
Definition: JAngle3D.hh:108
JGEOMETRY3D::JAngle3D::getDZ
double getDZ() const
Get z direction.
Definition: JAngle3D.hh:130
JGEOMETRY3D::JAngle3D::getDY
double getDY() const
Get y direction.
Definition: JAngle3D.hh:119
double JGEOMETRY3D::JPosition3D::getDot ( const JVersor3D &  dir) const
inlineinherited

Get dot product.

Parameters
dirdirection
Returns
dot product

Definition at line 393 of file JPosition3D.hh.

394  {
395  return
396  getX() * dir.getDX() +
397  getY() * dir.getDY() +
398  getZ() * dir.getDZ();
399  }
JGEOMETRY3D::JVersor3D::getDY
double getDY() const
Get y direction.
Definition: JVersor3D.hh:106
JGEOMETRY3D::JVersor3D::getDX
double getDX() const
Get x direction.
Definition: JVersor3D.hh:95
JGEOMETRY3D::JVector3D::getY
double getY() const
Get y position.
Definition: JVector3D.hh:104
JGEOMETRY3D::JVector3D::getX
double getX() const
Get x position.
Definition: JVector3D.hh:94
JGEOMETRY3D::JVector3D::getZ
double getZ() const
Get z position.
Definition: JVector3D.hh:115
JGEOMETRY3D::JVersor3D::getDZ
double getDZ() const
Get z direction.
Definition: JVersor3D.hh:117
double JGEOMETRY3D::JPosition3D::getDot ( const JVersor3Z &  dir) const
inlineinherited

Get dot product.

Parameters
dirdirection
Returns
dot product

Definition at line 408 of file JPosition3D.hh.

409  {
410  return
411  getX() * dir.getDX() +
412  getY() * dir.getDY() +
413  getZ() * dir.getDZ();
414  }
JGEOMETRY3D::JVersor3Z::getDY
double getDY() const
Get y direction.
Definition: JVersor3Z.hh:158
JGEOMETRY3D::JVector3D::getY
double getY() const
Get y position.
Definition: JVector3D.hh:104
JGEOMETRY3D::JVersor3Z::getDX
double getDX() const
Get x direction.
Definition: JVersor3Z.hh:147
JGEOMETRY3D::JVector3D::getX
double getX() const
Get x position.
Definition: JVector3D.hh:94
JGEOMETRY3D::JVersor3Z::getDZ
double getDZ() const
Get z direction.
Definition: JVersor3Z.hh:169
JGEOMETRY3D::JVector3D::getZ
double getZ() const
Get z position.
Definition: JVector3D.hh:115
double JGEOMETRY3D::JVector3D::getDot ( const JVector3D &  vector) const
inlineinherited

Get dot product.

Parameters
vectorvector
Returns
dot product

Definition at line 282 of file JVector3D.hh.

283  {
284  return
285  getX() * vector.getX() +
286  getY() * vector.getY() +
287  getZ() * vector.getZ();
288  }
JGEOMETRY3D::JVector3D::getY
double getY() const
Get y position.
Definition: JVector3D.hh:104
JGEOMETRY3D::JVector3D::getX
double getX() const
Get x position.
Definition: JVector3D.hh:94
JGEOMETRY3D::JVector3D::getZ
double getZ() const
Get z position.
Definition: JVector3D.hh:115
const JPosition3D& JGEOMETRY3D::JPosition3D::getPosition ( ) const
inlineinherited

Get position.

Returns
position

Definition at line 130 of file JPosition3D.hh.

131  {
132  return static_cast<const JPosition3D&>(*this);
133  }
JGEOMETRY3D::JPosition3D
Data structure for position in three dimensions.
Definition: JPosition3D.hh:36
JPosition3D& JGEOMETRY3D::JPosition3D::getPosition ( )
inlineinherited

Get position.

Returns
position

Definition at line 141 of file JPosition3D.hh.

142  {
143  return static_cast<JPosition3D&>(*this);
144  }
JGEOMETRY3D::JPosition3D
Data structure for position in three dimensions.
Definition: JPosition3D.hh:36
void JGEOMETRY3D::JPosition3D::setPosition ( const JVector3D &  pos)
inlineinherited

Set position.

Parameters
posposition

Definition at line 152 of file JPosition3D.hh.

153  {
154  static_cast<JVector3D&>(*this) = pos;
155  }
JGEOMETRY3D::JVector3D
Data structure for vector in three dimensions.
Definition: JVector3D.hh:34
JGEOMETRY3D::JPosition3D::operator JAngle3D ( ) const
inlineinherited

Type conversion operator.

Returns
angle

Definition at line 163 of file JPosition3D.hh.

164  {
165  return JAngle3D(getX(), getY(), getZ());
166  }
JGEOMETRY3D::JAngle3D
Data structure for angles in three dimensions.
Definition: JAngle3D.hh:33
JGEOMETRY3D::JVector3D::getY
double getY() const
Get y position.
Definition: JVector3D.hh:104
JGEOMETRY3D::JVector3D::getX
double getX() const
Get x position.
Definition: JVector3D.hh:94
JGEOMETRY3D::JVector3D::getZ
double getZ() const
Get z position.
Definition: JVector3D.hh:115
JGEOMETRY3D::JPosition3D::operator JVersor3D ( ) const
inlineinherited

Type conversion operator.

Returns
direction

Definition at line 174 of file JPosition3D.hh.

175  {
176  return JVersor3D(getX(), getY(), getZ());
177  }
JGEOMETRY3D::JVector3D::getY
double getY() const
Get y position.
Definition: JVector3D.hh:104
JGEOMETRY3D::JVector3D::getX
double getX() const
Get x position.
Definition: JVector3D.hh:94
JGEOMETRY3D::JVersor3D
Data structure for normalised vector in three dimensions.
Definition: JVersor3D.hh:26
JGEOMETRY3D::JVector3D::getZ
double getZ() const
Get z position.
Definition: JVector3D.hh:115
JPosition3D& JGEOMETRY3D::JPosition3D::rotate ( const JRotation3D &  R)
inlineinherited

Rotate.

Parameters
Rrotation matrix
Returns
this position

Definition at line 186 of file JPosition3D.hh.

187  {
188  R.rotate(__x, __y, __z);
189 
190  return *this;
191  }
JGEOMETRY3D::JRotation3D::rotate
void rotate(double &__x, double &__y, double &__z) const
Rotate.
Definition: JRotation3D.hh:336
JPosition3D& JGEOMETRY3D::JPosition3D::rotate ( const JRotation3X &  R)
inlineinherited

Rotate around X-axis.

Parameters
Rrotation matrix
Returns
this position

Definition at line 214 of file JPosition3D.hh.

215  {
216  R.rotate(__y, __z);
217 
218  return *this;
219  }
JGEOMETRY2D::JRotation2D::rotate
void rotate(double &__x, double &__y) const
Rotate.
Definition: JRotation2D.hh:96
JPosition3D& JGEOMETRY3D::JPosition3D::rotate ( const JRotation3Y &  R)
inlineinherited

Rotate around Y-axis.

Parameters
Rrotation matrix
Returns
this position

Definition at line 242 of file JPosition3D.hh.

243  {
244  R.rotate(__x, __z);
245 
246  return *this;
247  }
JGEOMETRY2D::JRotation2D::rotate
void rotate(double &__x, double &__y) const
Rotate.
Definition: JRotation2D.hh:96
JPosition3D& JGEOMETRY3D::JPosition3D::rotate ( const JRotation3Z &  R)
inlineinherited

Rotate around Z-axis.

Parameters
Rrotation matrix
Returns
this position

Definition at line 270 of file JPosition3D.hh.

271  {
272  R.rotate(__x, __y);
273 
274  return *this;
275  }
JGEOMETRY2D::JRotation2D::rotate
void rotate(double &__x, double &__y) const
Rotate.
Definition: JRotation2D.hh:96
JPosition3D& JGEOMETRY3D::JPosition3D::rotate ( const JQuaternion3D &  Q)
inlineinherited

Rotate.

Parameters
Qquaternion
Returns
this position

Definition at line 298 of file JPosition3D.hh.

299  {
300  Q.rotate(__x, __y, __z);
301 
302  return *this;
303  }
JGEOMETRY3D::JQuaternion3D::rotate
void rotate(double &__x, double &__y, double &__z) const
Rotate.
Definition: JQuaternion3D.hh:805
JPosition3D& JGEOMETRY3D::JPosition3D::rotate_back ( const JRotation3D &  R)
inlineinherited

Rotate back.

Parameters
Rrotation matrix
Returns
this position

Definition at line 200 of file JPosition3D.hh.

201  {
202  R.rotate_back(__x, __y, __z);
203 
204  return *this;
205  }
JGEOMETRY3D::JRotation3D::rotate_back
void rotate_back(double &__x, double &__y, double &__z) const
Rotate back.
Definition: JRotation3D.hh:355
JPosition3D& JGEOMETRY3D::JPosition3D::rotate_back ( const JRotation3X &  R)
inlineinherited

Rotate back around X-axis.

Parameters
Rrotation matrix
Returns
this position

Definition at line 228 of file JPosition3D.hh.

229  {
230  R.rotate_back(__y, __z);
231 
232  return *this;
233  }
JGEOMETRY2D::JRotation2D::rotate_back
void rotate_back(double &__x, double &__y) const
Rotate back.
Definition: JRotation2D.hh:112
JPosition3D& JGEOMETRY3D::JPosition3D::rotate_back ( const JRotation3Y &  R)
inlineinherited

Rotate back around Y-axis.

Parameters
Rrotation matrix
Returns
this position

Definition at line 256 of file JPosition3D.hh.

257  {
258  R.rotate_back(__x, __z);
259 
260  return *this;
261  }
JGEOMETRY2D::JRotation2D::rotate_back
void rotate_back(double &__x, double &__y) const
Rotate back.
Definition: JRotation2D.hh:112
JPosition3D& JGEOMETRY3D::JPosition3D::rotate_back ( const JRotation3Z &  R)
inlineinherited

Rotate back around Z-axis.

Parameters
Rrotation matrix
Returns
this position

Definition at line 284 of file JPosition3D.hh.

285  {
286  R.rotate_back(__x, __y);
287 
288  return *this;
289  }
JGEOMETRY2D::JRotation2D::rotate_back
void rotate_back(double &__x, double &__y) const
Rotate back.
Definition: JRotation2D.hh:112
JPosition3D& JGEOMETRY3D::JPosition3D::rotate_back ( const JQuaternion3D &  Q)
inlineinherited

Rotate back.

Parameters
Qquaternion
Returns
this position

Definition at line 312 of file JPosition3D.hh.

313  {
314  Q.rotate_back(__x, __y, __z);
315 
316  return *this;
317  }
JGEOMETRY3D::JQuaternion3D::rotate_back
void rotate_back(double &__x, double &__y, double &__z) const
Rotate back.
Definition: JQuaternion3D.hh:824
void JGEOMETRY3D::JPosition3D::transform ( const JRotation3D &  R,
const JVector3D &  pos 
)
inlineinherited

Transform position.

The final position is obtained as follows:

  1. rotation of the position according matrix R;
  2. offset position with pos;
  3. rotation of position around z-axis, such that final position lies in x-z plane;
Parameters
Rrotation matrix
posposition of origin (after rotation)

Definition at line 331 of file JPosition3D.hh.

333  {
334  // rotate geometry to system with particle direction along z-axis
335 
336  rotate(R);
337 
338  // offset with respect to origin
339 
340  sub(pos);
341 
342  // rotate geometry to x-z plane
343 
344  __x = sqrt(__x*__x + __y*__y);
345  __y = 0.0;
346  }
JGEOMETRY3D::JVector3D::sub
JVector3D & sub(const JVector3D &vector)
Subtract vector.
Definition: JVector3D.hh:158
JGEOMETRY3D::JPosition3D::rotate
JPosition3D & rotate(const JRotation3D &R)
Rotate.
Definition: JPosition3D.hh:186
JVector3D& JGEOMETRY3D::JVector3D::transform ( const JMatrix3D &  T)
inlineinherited

Transform.

Parameters
Tmatrix
Returns
this vector

Definition at line 206 of file JVector3D.hh.

207  {
208  T.transform(__x, __y, __z);
209 
210  return *this;
211  }
T
do set_variable OUTPUT_DIRECTORY $WORKDIR T
Definition: JCalibrateHeight.sh:61
void JGEOMETRY3D::JPosition3D::transform_back ( const JRotation3D &  R,
const JVector3D &  pos 
)
inlineinherited

Transform back position.

The final position is obtained as follows:

  1. offset position with position pos;
  2. rotation of postion according matrix R;
Parameters
Rrotation matrix
posposition of origin (before rotation)

Definition at line 359 of file JPosition3D.hh.

361  {
362  // offset with respect to origin
363 
364  add(pos);
365 
366  // rotate back geometry to system with particle direction along z-axis
367 
368  rotate_back(R);
369  }
JGEOMETRY3D::JPosition3D::rotate_back
JPosition3D & rotate_back(const JRotation3D &R)
Rotate back.
Definition: JPosition3D.hh:200
JGEOMETRY3D::JVector3D::add
JVector3D & add(const JVector3D &vector)
Add vector.
Definition: JVector3D.hh:142
JGEOMETRY3D::JVector3D::operator JVector2D ( ) const
inlineinherited

Type conversion operator.

Returns
JVector2D

Definition at line 83 of file JVector3D.hh.

84  {
85  return JVector2D(this->getX(), this->getY());
86  }
JGEOMETRY3D::JVector3D::getY
double getY() const
Get y position.
Definition: JVector3D.hh:104
JGEOMETRY3D::JVector3D::getX
double getX() const
Get x position.
Definition: JVector3D.hh:94
double JGEOMETRY3D::JVector3D::getX ( ) const
inlineinherited

Get x position.

Returns
x position

Definition at line 94 of file JVector3D.hh.

95  {
96  return __x;
97  }
double JGEOMETRY3D::JVector3D::getY ( ) const
inlineinherited

Get y position.

Returns
y position

Definition at line 104 of file JVector3D.hh.

105  {
106  return __y;
107  }
double JGEOMETRY3D::JVector3D::getZ ( ) const
inlineinherited

Get z position.

Returns
z position

Definition at line 115 of file JVector3D.hh.

116  {
117  return __z;
118  }
bool JGEOMETRY3D::JVector3D::equals ( const JVector3D &  vector,
const double  precision = std::numeric_limits<double>::min() 
) const
inlineinherited

Check equality.

Parameters
vectorvector
precisionprecision
Returns
true if vectors are equal; else false

Definition at line 221 of file JVector3D.hh.

223  {
224  return (fabs(getX() - vector.getX()) <= precision &&
225  fabs(getY() - vector.getY()) <= precision &&
226  fabs(getZ() - vector.getZ()) <= precision);
227  }
JGEOMETRY3D::JVector3D::getY
double getY() const
Get y position.
Definition: JVector3D.hh:104
JGEOMETRY3D::JVector3D::getX
double getX() const
Get x position.
Definition: JVector3D.hh:94
JGEOMETRY3D::JVector3D::getZ
double getZ() const
Get z position.
Definition: JVector3D.hh:115
double JGEOMETRY3D::JVector3D::getLengthSquared ( ) const
inlineinherited

Get length squared.

Returns
square of length

Definition at line 235 of file JVector3D.hh.

236  {
237  return getX()*getX() + getY()*getY() + getZ()*getZ();
238  }
JGEOMETRY3D::JVector3D::getY
double getY() const
Get y position.
Definition: JVector3D.hh:104
JGEOMETRY3D::JVector3D::getX
double getX() const
Get x position.
Definition: JVector3D.hh:94
JGEOMETRY3D::JVector3D::getZ
double getZ() const
Get z position.
Definition: JVector3D.hh:115
double JGEOMETRY3D::JVector3D::getLength ( ) const
inlineinherited

Get length.

Returns
length

Definition at line 246 of file JVector3D.hh.

247  {
248  return sqrt(getLengthSquared());
249  }
JGEOMETRY3D::JVector3D::getLengthSquared
double getLengthSquared() const
Get length squared.
Definition: JVector3D.hh:235
double JGEOMETRY3D::JVector3D::getDistanceSquared ( const JVector3D &  pos) const
inlineinherited

Get squared of distance to point.

Parameters
posposition
Returns
square of distance

Definition at line 258 of file JVector3D.hh.

259  {
260  return JVector3D(pos).sub(*this).getLengthSquared();
261  }
JGEOMETRY3D::JVector3D::JVector3D
JVector3D()
Default constructor.
Definition: JVector3D.hh:41
double JGEOMETRY3D::JVector3D::getDistance ( const JVector3D &  pos) const
inlineinherited

Get distance to point.

Parameters
posposition
Returns
distance

Definition at line 270 of file JVector3D.hh.

271  {
272  return sqrt(getDistanceSquared(pos));
273  }
JGEOMETRY3D::JVector3D::getDistanceSquared
double getDistanceSquared(const JVector3D &pos) const
Get squared of distance to point.
Definition: JVector3D.hh:258
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
firstfirst vector
secondsecond vector
Returns
this vector

Definition at line 299 of file JVector3D.hh.

301  {
302  __x = first .getY() * second.getZ() - second.getY() * first .getZ();
303  __y = second.getX() * first .getZ() - first .getX() * second.getZ();
304  __z = first .getX() * second.getY() - second.getX() * first .getY();
305 
306  return *this;
307  }
JGEOMETRY3D::JVector3D::getY
double getY() const
Get y position.
Definition: JVector3D.hh:104
JGEOMETRY3D::JVector3D::getX
double getX() const
Get x position.
Definition: JVector3D.hh:94
JGEOMETRY3D::JVector3D::getZ
double getZ() const
Get z position.
Definition: JVector3D.hh:115

Member Data Documentation

const int JFIT::JEstimator< JPoint4D >::NUMBER_OF_PARAMETERS = 4
static

number of parameters of fit

Definition at line 223 of file JPoint4DEstimator.hh.

JMATH::JMatrix4S JFIT::JEstimator< JPoint4D >::V

co-variance matrix of fit parameters

Definition at line 224 of file JPoint4DEstimator.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:
Generated by   doxygen 1.8.5