Jpp
Public Types | Public Member Functions | Static Public Member Functions | Public Attributes | Static Public Attributes | Protected Attributes | List of all members
JFIT::JEstimator< JLine1Z > Class Template Reference

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 >:
JFIT::JLine1Z JGEOMETRY3D::JVertex3D JGEOMETRY3D::JPosition3D JGEOMETRY3D::JTime JMATH::JMath< JVertex3D > JGEOMETRY3D::JVector3D JMATH::JMath< JVector3D >

Public Types

typedef double JLine1Z::* parameter_type
 

Public Member Functions

template<class T >
 JEstimator (T __begin, T __end)
 Fit constructor. More...
 
template<class T >
void update (T __begin, T __end, const JMatrixNZ &A)
 Update track parameters using updated co-variance matrix (e.g. 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 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...
 
double getDot (const JVector3D &vector) const
 Get dot product. More...
 
JVertex3Dnegate ()
 Prefix unary minus. More...
 
JVertex3Dadd (const JVertex3D &value)
 Addition operator. More...
 
JVector3Dadd (const JVector3D &vector)
 Add vector. More...
 
JTimeadd (const JTime &value)
 Addition operator. More...
 
JVertex3Dsub (const JVertex3D &value)
 Subtraction operator. More...
 
JVector3Dsub (const JVector3D &vector)
 Subtract vector. More...
 
JTimesub (const JTime &value)
 Subtraction operator. More...
 
JVertex3Dmul (const double value)
 Multiplication operator. More...
 
JVector3Dmul (const JSecond_t &object)
 Multiply with object. More...
 
JVertex3Dmul (const JSecond_t &object)
 Multiply with object. More...
 
JVertex3Ddiv (const double value)
 Division operator. More...
 
const JPosition3DgetPosition () const
 Get position. More...
 
JPosition3DgetPosition ()
 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...
 
JPosition3Drotate (const JRotation3D &R)
 Rotate. More...
 
JPosition3Drotate (const JRotation3X &R)
 Rotate around X-axis. More...
 
JPosition3Drotate (const JRotation3Y &R)
 Rotate around Y-axis. More...
 
JPosition3Drotate (const JRotation3Z &R)
 Rotate around Z-axis. More...
 
JPosition3Drotate (const JQuaternion3D &Q)
 Rotate. More...
 
JPosition3Drotate_back (const JRotation3D &R)
 Rotate back. More...
 
JPosition3Drotate_back (const JRotation3X &R)
 Rotate back around X-axis. More...
 
JPosition3Drotate_back (const JRotation3Y &R)
 Rotate back around Y-axis. More...
 
JPosition3Drotate_back (const JRotation3Z &R)
 Rotate back around Z-axis. More...
 
JPosition3Drotate_back (const JQuaternion3D &Q)
 Rotate back. More...
 
void transform (const JRotation3D &R, const JVector3D &pos)
 Transform position. More...
 
JVector3Dtransform (const JMatrix3D &T)
 Transform. More...
 
JVector3Dtransform (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...
 
JVector3Dcross (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...
 

Static Public Attributes

static const int NUMBER_OF_PARAMETERS = 3
 number of parameters of fit 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 88 of file JLine1ZEstimator.hh.

Member Typedef Documentation

◆ parameter_type

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

Definition at line 178 of file JLine1Z.hh.

Constructor & Destructor Documentation

◆ JEstimator()

template<class T >
JFIT::JEstimator< JLine1Z >::JEstimator ( __begin,
__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 105 of file JLine1ZEstimator.hh.

105  :
106  JLine1Z()
107  {
108  using namespace std;
109  using namespace JTOOLS;
110 
111  const int N = distance(__begin, __end);
112 
113  if (N >= NUMBER_OF_PARAMETERS) {
114 
115  double t0 = 0.0;
116 
117  for (T i = __begin; i != __end; ++i) {
118  __x += i->getX();
119  __y += i->getY();
120  __z += i->getZ();
121  t0 += i->getT();
122  }
123 
124  div(N);
125  t0 /= N;
126 
127  V.reset();
128 
129  t0 *= getSpeedOfLight();
130 
131  double y0 = 0.0;
132  double y1 = 0.0;
133  double y2 = 0.0;
134 
135  T j = __begin;
136 
137  double xi = j->getX() - getX();
138  double yi = j->getY() - getY();
139  double ti = (j->getT() * getSpeedOfLight() - t0 - j->getZ() + getZ()) / getKappaC();
140 
141  for (bool done = false; !done; ) {
142 
143  if ((done = (++j == __end))) {
144  j = __begin;
145  }
146 
147  double xj = j->getX() - getX();
148  double yj = j->getY() - getY();
149  double tj = (j->getT() * getSpeedOfLight() - t0 - j->getZ() + getZ()) / getKappaC();
150 
151  double dx = xj - xi;
152  double dy = yj - yi;
153  double dt = ti - tj; // opposite sign!
154 
155  const double y = ((xj + xi) * dx +
156  (yj + yi) * dy +
157  (tj + ti) * dt);
158 
159  dx *= 2;
160  dy *= 2;
161  dt *= 2;
162 
163  V.a00 += dx * dx;
164  V.a01 += dx * dy;
165  V.a02 += dx * dt;
166  V.a11 += dy * dy;
167  V.a12 += dy * dt;
168  V.a22 += dt * dt;
169 
170  y0 += dx * y;
171  y1 += dy * y;
172  y2 += dt * y;
173 
174  xi = xj;
175  yi = yj;
176  ti = tj;
177  }
178 
179  t0 *= getInverseSpeedOfLight();
180 
181  V.a10 = V.a01;
182  V.a20 = V.a02;
183  V.a21 = V.a12;
184 
185  V.invert();
186 
187  __x += V.a00 * y0 + V.a01 * y1 + V.a02 * y2;
188  __y += V.a10 * y0 + V.a11 * y1 + V.a12 * y2;
189  __t = V.a20 * y0 + V.a21 * y1 + V.a22 * y2;
190 
192  __t += t0;
193 
194  } else {
195  throw JValueOutOfRange("JEstimator<JLine1Z>::JEstimator(): Not enough data points.");
196  }
197  }

Member Function Documentation

◆ update()

template<class T >
void JFIT::JEstimator< JLine1Z >::update ( __begin,
__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
__beginbegin of data
__endend of data
Aco-variance matrix of hits

Definition at line 217 of file JLine1ZEstimator.hh.

218  {
219  using namespace std;
220  using namespace JTOOLS;
221 
222  const int N = distance(__begin, __end);
223 
224  if (N >= NUMBER_OF_PARAMETERS) {
225 
226  double x1 = 0.0;
227  double y1 = 0.0;
228  double t1 = 0.0;
229 
230  V.reset();
231 
232  T i = __begin;
233 
234  for (JMatrixNZ::const_row_type row = A.begin(); row != A.end(); ++row, ++i) {
235 
236  const double dx = i->getX() - getX();
237  const double dy = i->getY() - getY();
238 
239  const double rt = sqrt(dx*dx + dy*dy);
240 
241  double xr = getKappaC() * getInverseSpeedOfLight();
242  double yr = getKappaC() * getInverseSpeedOfLight();
243  double tr = 1.0;
244 
245  if (rt != 0.0) {
246  xr *= -dx / rt;
247  yr *= -dy / rt;
248  }
249 
250  T j = __begin;
251 
252  for (JMatrixNZ::const_col_type col = row->begin(); col != row->end(); ++col, ++j) {
253 
254  const double dx = j->getX() - getX();
255  const double dy = j->getY() - getY();
256  const double dz = j->getZ() - getZ();
257 
258  const double rt = sqrt(dx*dx + dy*dy);
259 
260  double xc = getKappaC() * getInverseSpeedOfLight();
261  double yc = getKappaC() * getInverseSpeedOfLight();
262  double tc = 1.0;
263 
264  if (rt != 0.0) {
265  xc *= -dx / rt;
266  yc *= -dy / rt;
267  }
268 
269  const double ts = j->getT() - (dz + rt * getKappaC()) * getInverseSpeedOfLight();
270 
271  x1 += xr * (*col) * ts;
272  y1 += yr * (*col) * ts;
273  t1 += tr * (*col) * ts;
274 
275  V.a00 += xr * (*col) * xc;
276  V.a01 += xr * (*col) * yc;
277  V.a02 += xr * (*col) * tc;
278  V.a11 += yr * (*col) * yc;
279  V.a12 += yr * (*col) * tc;
280  V.a22 += tr * (*col) * tc;
281  }
282  }
283 
284  V.a10 = V.a01;
285  V.a20 = V.a02;
286  V.a21 = V.a12;
287 
288  V.invert();
289 
290  __x += V.a00 * x1 + V.a01 * y1 + V.a02 * t1;
291  __y += V.a10 * x1 + V.a11 * y1 + V.a12 * t1;
292  __t = V.a20 * x1 + V.a21 * y1 + V.a22 * t1;
293 
294  } else {
295  throw JValueOutOfRange("JEstimator<JLine1Z>::update(): Not enough data points.");
296  }
297  }

◆ move()

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

Move vertex along this line with given velocity.

Parameters
stepstep
velocityvelocity

Definition at line 62 of file JLine1Z.hh.

63  {
64  __z += step;
65  __t += step / velocity;
66  }

◆ setZ()

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

Set z-position of vertex.

Parameters
zz-position
velocityvelocity

Definition at line 75 of file JLine1Z.hh.

76  {
77  move(z - getZ(), velocity);
78  }

◆ getDistanceSquared()

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

Get distance squared.

Parameters
posposition
Returns
square of distance

Definition at line 87 of file JLine1Z.hh.

88  {
89  const double dx = pos.getX() - this->getX();
90  const double dy = pos.getY() - this->getY();
91 
92  return dx*dx + dy*dy;
93  }

◆ getDistance()

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

Get distance.

Parameters
posposition
Returns
distance

Definition at line 102 of file JLine1Z.hh.

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

◆ getT() [1/3]

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

Get arrival time of Cherenkov light at given position.

Parameters
posposition [m]
Returns
time [ns]

Definition at line 114 of file JLine1Z.hh.

115  {
116  using namespace JTOOLS;
117 
118  const double dx = pos.getX() - this->getX();
119  const double dy = pos.getY() - this->getY();
120  const double dz = pos.getZ() - this->getZ();
121 
122  const double R = sqrt(dx*dx + dy*dy);
123 
124  return this->getT() + (dz + R * getKappaC()) * getInverseSpeedOfLight();
125  }

◆ getT() [2/3]

double JGEOMETRY3D::JVertex3D::getT
inlineinherited

Get arrival time of Cherenkov light at given position.

Parameters
posposition [m]
Returns
time [ns]

Definition at line 144 of file JVertex3D.hh.

145  {
146  using namespace JTOOLS;
147 
148  return this->getT() + this->getDistance(pos) * getInverseSpeedOfLight() * getIndexOfRefraction();
149  }

◆ getT() [3/3]

double JGEOMETRY3D::JTime::getT ( ) const
inlineinherited

Get time.

Returns
time

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

122  {
123  return __t;
124  }

◆ getZ() [1/2]

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

Get point of emission of Cherenkov light along muon path.

Parameters
posposition
Returns
position along muon path

Definition at line 134 of file JLine1Z.hh.

135  {
136  using namespace JTOOLS;
137 
138  return pos.getZ() - this->getDistance(pos) / getTanThetaC();
139  }

◆ getZ() [2/2]

double JGEOMETRY3D::JVector3D::getZ ( ) const
inlineinherited

Get z position.

Returns
z position

Definition at line 114 of file JVector3D.hh.

115  {
116  return __z;
117  }

◆ getDirection()

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

Get photon direction of Cherenkov light on PMT.

Parameters
posPMT position
Returns
direction

Definition at line 148 of file JLine1Z.hh.

149  {
150  using namespace JTOOLS;
151 
152  double dx = pos.getX() - this->getX();
153  double dy = pos.getY() - this->getY();
154 
155  const double R = sqrt(dx*dx + dy*dy);
156 
157  dx *= getSinThetaC() / R;
158  dy *= getSinThetaC() / R;
159 
160  const double dz = getCosThetaC();
161 
162  return JVersor3D(dx,dy,dz);
163  }

◆ getDot() [1/6]

double JFIT::JLine1Z::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 172 of file JLine1Z.hh.

173  {
174  return getDirection(axis.getPosition()).getDot(axis.getDirection());
175  }

◆ getDot() [2/6]

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

Get dot product.

Parameters
angleangle
Returns
dot product

Definition at line 377 of file JPosition3D.hh.

378  {
379  return
380  getX() * angle.getDX() +
381  getY() * angle.getDY() +
382  getZ() * angle.getDZ();
383  }

◆ getDot() [3/6]

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

Get dot product.

Parameters
dirdirection
Returns
dot product

Definition at line 392 of file JPosition3D.hh.

393  {
394  return
395  getX() * dir.getDX() +
396  getY() * dir.getDY() +
397  getZ() * dir.getDZ();
398  }

◆ getDot() [4/6]

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

Get dot product.

Parameters
dirdirection
Returns
dot product

Definition at line 407 of file JPosition3D.hh.

408  {
409  return
410  getX() * dir.getDX() +
411  getY() * dir.getDY() +
412  getZ() * dir.getDZ();
413  }

◆ getDot() [5/6]

double JGEOMETRY3D::JVector3D::getDot
inlineinherited

Get dot product.

Parameters
vectorvector
Returns
dot product

Definition at line 280 of file JVector3D.hh.

281  {
282  return
283  getX() * vector.getX() +
284  getY() * vector.getY() +
285  getZ() * vector.getZ();
286  }

◆ getDot() [6/6]

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

Get dot product.

Parameters
vectorvector
Returns
dot product

Definition at line 280 of file JVector3D.hh.

281  {
282  return
283  getX() * vector.getX() +
284  getY() * vector.getY() +
285  getZ() * vector.getZ();
286  }

◆ pX()

static parameter_type JFIT::JLine1Z::pX ( )
inlinestaticinherited

Definition at line 180 of file JLine1Z.hh.

180 { return &JLine1Z::__x; }

◆ pY()

static parameter_type JFIT::JLine1Z::pY ( )
inlinestaticinherited

Definition at line 181 of file JLine1Z.hh.

181 { return &JLine1Z::__y; }

◆ pT()

static parameter_type JFIT::JLine1Z::pT ( )
inlinestaticinherited

Definition at line 182 of file JLine1Z.hh.

182 { return &JLine1Z::__t; }

◆ negate()

JVertex3D& JGEOMETRY3D::JVertex3D::negate ( )
inlineinherited

Prefix unary minus.

Returns
line

Definition at line 70 of file JVertex3D.hh.

71  {
74 
75  return *this;
76  }

◆ add() [1/3]

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

Addition operator.

Parameters
valueline
Returns
line

Definition at line 84 of file JVertex3D.hh.

85  {
86  JPosition3D::add(value);
87  JTime ::add(value);
88 
89  return *this;
90  }

◆ add() [2/3]

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

Add vector.

Parameters
vectorvector
Returns
this vector

Definition at line 141 of file JVector3D.hh.

142  {
143  __x += vector.getX();
144  __y += vector.getY();
145  __z += vector.getZ();
146 
147  return *this;
148  }

◆ add() [3/3]

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

Addition operator.

Parameters
valuetime
Returns
time

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

67  {
68  __t += value.getT();
69 
70  return *this;
71  }

◆ sub() [1/3]

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

Subtraction operator.

Parameters
valueline
Returns
line

Definition at line 99 of file JVertex3D.hh.

100  {
101  JPosition3D::sub(value);
102  JTime ::sub(value);
103 
104  return *this;
105  }

◆ sub() [2/3]

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

Subtract vector.

Parameters
vectorvector
Returns
this vector

Definition at line 157 of file JVector3D.hh.

158  {
159  __x -= vector.getX();
160  __y -= vector.getY();
161  __z -= vector.getZ();
162 
163  return *this;
164  }

◆ sub() [3/3]

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

Subtraction operator.

Parameters
valuetime
Returns
time

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

81  {
82  __t -= value.getT();
83 
84  return *this;
85  }

◆ mul() [1/3]

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

Multiplication operator.

Parameters
valuemultiplication factor
Returns
line

Definition at line 114 of file JVertex3D.hh.

115  {
116  JPosition3D::mul(value);
117  JTime ::mul(value);
118 
119  return *this;
120  }

◆ mul() [2/3]

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

Multiply with object.

Parameters
objectobject
Returns
result object

Definition at line 273 of file JMath.hh.

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

◆ mul() [3/3]

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

Multiply with object.

Parameters
objectobject
Returns
result object

Definition at line 273 of file JMath.hh.

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

◆ div()

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

Division operator.

Parameters
valuemultiplication factor
Returns
line

Definition at line 129 of file JVertex3D.hh.

130  {
131  JPosition3D::div(value);
132  JTime ::div(value);
133 
134  return *this;
135  }

◆ getPosition() [1/2]

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

Get position.

Returns
position

Definition at line 129 of file JPosition3D.hh.

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

◆ getPosition() [2/2]

JPosition3D& JGEOMETRY3D::JPosition3D::getPosition ( )
inlineinherited

Get position.

Returns
position

Definition at line 140 of file JPosition3D.hh.

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

◆ setPosition()

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

Set position.

Parameters
posposition

Definition at line 151 of file JPosition3D.hh.

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

◆ operator JAngle3D()

JGEOMETRY3D::JPosition3D::operator JAngle3D ( ) const
inlineinherited

Type conversion operator.

Returns
angle

Definition at line 162 of file JPosition3D.hh.

163  {
164  return JAngle3D(getX(), getY(), getZ());
165  }

◆ operator JVersor3D()

JGEOMETRY3D::JPosition3D::operator JVersor3D ( ) const
inlineinherited

Type conversion operator.

Returns
direction

Definition at line 173 of file JPosition3D.hh.

174  {
175  return JVersor3D(getX(), getY(), getZ());
176  }

◆ rotate() [1/5]

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

Rotate.

Parameters
Rrotation matrix
Returns
this position

Definition at line 185 of file JPosition3D.hh.

186  {
187  R.rotate(__x, __y, __z);
188 
189  return *this;
190  }

◆ rotate() [2/5]

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

Rotate around X-axis.

Parameters
Rrotation matrix
Returns
this position

Definition at line 213 of file JPosition3D.hh.

214  {
215  R.rotate(__y, __z);
216 
217  return *this;
218  }

◆ rotate() [3/5]

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

Rotate around Y-axis.

Parameters
Rrotation matrix
Returns
this position

Definition at line 241 of file JPosition3D.hh.

242  {
243  R.rotate(__x, __z);
244 
245  return *this;
246  }

◆ rotate() [4/5]

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

Rotate around Z-axis.

Parameters
Rrotation matrix
Returns
this position

Definition at line 269 of file JPosition3D.hh.

270  {
271  R.rotate(__x, __y);
272 
273  return *this;
274  }

◆ rotate() [5/5]

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

Rotate.

Parameters
Qquaternion
Returns
this position

Definition at line 297 of file JPosition3D.hh.

298  {
299  Q.rotate(__x, __y, __z);
300 
301  return *this;
302  }

◆ rotate_back() [1/5]

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

Rotate back.

Parameters
Rrotation matrix
Returns
this position

Definition at line 199 of file JPosition3D.hh.

200  {
201  R.rotate_back(__x, __y, __z);
202 
203  return *this;
204  }

◆ rotate_back() [2/5]

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

Rotate back around X-axis.

Parameters
Rrotation matrix
Returns
this position

Definition at line 227 of file JPosition3D.hh.

228  {
229  R.rotate_back(__y, __z);
230 
231  return *this;
232  }

◆ rotate_back() [3/5]

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

Rotate back around Y-axis.

Parameters
Rrotation matrix
Returns
this position

Definition at line 255 of file JPosition3D.hh.

256  {
257  R.rotate_back(__x, __z);
258 
259  return *this;
260  }

◆ rotate_back() [4/5]

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

Rotate back around Z-axis.

Parameters
Rrotation matrix
Returns
this position

Definition at line 283 of file JPosition3D.hh.

284  {
285  R.rotate_back(__x, __y);
286 
287  return *this;
288  }

◆ rotate_back() [5/5]

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

Rotate back.

Parameters
Qquaternion
Returns
this position

Definition at line 311 of file JPosition3D.hh.

312  {
313  Q.rotate_back(__x, __y, __z);
314 
315  return *this;
316  }

◆ transform() [1/3]

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 330 of file JPosition3D.hh.

332  {
333  // rotate geometry to system with particle direction along z-axis
334 
335  rotate(R);
336 
337  // offset with respect to origin
338 
339  sub(pos);
340 
341  // rotate geometry to x-z plane
342 
343  __x = sqrt(__x*__x + __y*__y);
344  __y = 0.0;
345  }

◆ transform() [2/3]

JVector3D& JGEOMETRY3D::JVector3D::transform
inlineinherited

Transform.

Parameters
Tmatrix
Returns
this vector

Definition at line 205 of file JVector3D.hh.

206  {
207  T.transform(__x, __y, __z);
208 
209  return *this;
210  }

◆ transform() [3/3]

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

Transform.

Parameters
Tmatrix
Returns
this vector

Definition at line 205 of file JVector3D.hh.

206  {
207  T.transform(__x, __y, __z);
208 
209  return *this;
210  }

◆ transform_back()

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 358 of file JPosition3D.hh.

360  {
361  // offset with respect to origin
362 
363  add(pos);
364 
365  // rotate back geometry to system with particle direction along z-axis
366 
367  rotate_back(R);
368  }

◆ operator JVector2D()

JGEOMETRY3D::JVector3D::operator JVector2D ( ) const
inlineinherited

Type conversion operator.

Returns
JVector2D

Definition at line 82 of file JVector3D.hh.

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

◆ getX()

double JGEOMETRY3D::JVector3D::getX ( ) const
inlineinherited

Get x position.

Returns
x position

Definition at line 93 of file JVector3D.hh.

94  {
95  return __x;
96  }

◆ getY()

double JGEOMETRY3D::JVector3D::getY ( ) const
inlineinherited

Get y position.

Returns
y position

Definition at line 103 of file JVector3D.hh.

104  {
105  return __y;
106  }

◆ equals()

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

Check equality.

Parameters
vectorvector
Returns
true if vectors are equal; else false

Definition at line 219 of file JVector3D.hh.

221  {
222  return (fabs(getX() - vector.getX()) <= precision &&
223  fabs(getY() - vector.getY()) <= precision &&
224  fabs(getZ() - vector.getZ()) <= precision);
225  }

◆ getLengthSquared()

double JGEOMETRY3D::JVector3D::getLengthSquared ( ) const
inlineinherited

Get length squared.

Returns
square of length

Definition at line 233 of file JVector3D.hh.

234  {
235  return getX()*getX() + getY()*getY() + getZ()*getZ();
236  }

◆ getLength()

double JGEOMETRY3D::JVector3D::getLength ( ) const
inlineinherited

Get length.

Returns
length

Definition at line 244 of file JVector3D.hh.

245  {
246  return sqrt(getLengthSquared());
247  }

◆ cross()

JVector3D& JGEOMETRY3D::JVector3D::cross ( 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 297 of file JVector3D.hh.

299  {
300  __x = first .getY() * second.getZ() - second.getY() * first .getZ();
301  __y = second.getX() * first .getZ() - first .getX() * second.getZ();
302  __z = first .getX() * second.getY() - second.getX() * first .getY();
303 
304  return *this;
305  }

Member Data Documentation

◆ NUMBER_OF_PARAMETERS

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

number of parameters of fit

Definition at line 300 of file JLine1ZEstimator.hh.

◆ V

co-variance matrix of fit parameters

Definition at line 301 of file JLine1ZEstimator.hh.

◆ __x

double JGEOMETRY3D::JVector3D::__x
protectedinherited

Definition at line 308 of file JVector3D.hh.

◆ __y

double JGEOMETRY3D::JVector3D::__y
protectedinherited

Definition at line 309 of file JVector3D.hh.

◆ __z

double JGEOMETRY3D::JVector3D::__z
protectedinherited

Definition at line 310 of file JVector3D.hh.

◆ __t

double JGEOMETRY3D::JTime::__t
protectedinherited

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


The documentation for this class was generated from the following file:
JGEOMETRY3D::JVersor3D::getDZ
double getDZ() const
Get z direction.
Definition: JVersor3D.hh:113
JMATH::JMatrix3D::a01
double a01
Definition: JMath/JMatrix3D.hh:385
JFIT::JLine1Z::getDirection
JVersor3D getDirection(const JVector3D &pos) const
Get photon direction of Cherenkov light on PMT.
Definition: JLine1Z.hh:148
JTOOLS::getSinThetaC
double getSinThetaC()
Get average sine of Cherenkov angle of water.
Definition: JConstants.hh:155
JGEOMETRY2D::JRotation2D::rotate_back
void rotate_back(double &__x, double &__y) const
Rotate back.
Definition: JRotation2D.hh:112
JMATH::JMatrix3S::invert
void invert()
Invert matrix.
Definition: JMatrix3S.hh:70
JMATH::JMatrix3D::a21
double a21
Definition: JMath/JMatrix3D.hh:387
JGEOMETRY3D::JRotation3D::rotate
void rotate(double &__x, double &__y, double &__z) const
Rotate.
Definition: JRotation3D.hh:336
JMATH::JMatrix3D::a20
double a20
Definition: JMath/JMatrix3D.hh:387
JTOOLS::getSpeedOfLight
const double getSpeedOfLight()
Number of bytes in a gigabyte.
Definition: JConstants.hh:89
JGEOMETRY3D::JVector3D::__y
double __y
Definition: JVector3D.hh:309
JGEOMETRY3D::JRotation3D::rotate_back
void rotate_back(double &__x, double &__y, double &__z) const
Rotate back.
Definition: JRotation3D.hh:355
JGEOMETRY3D::JVector3D::getLengthSquared
double getLengthSquared() const
Get length squared.
Definition: JVector3D.hh:233
JGEOMETRY3D::JQuaternion3D::rotate_back
void rotate_back(double &__x, double &__y, double &__z) const
Rotate back.
Definition: JQuaternion3D.hh:668
JGEOMETRY3D::JVertex3D::div
JVertex3D & div(const double value)
Division operator.
Definition: JVertex3D.hh:129
JGEOMETRY3D::JVector3D::getZ
double getZ() const
Get z position.
Definition: JVector3D.hh:114
JFIT::JEstimator< JLine1Z >::V
JMATH::JMatrix3S V
co-variance matrix of fit parameters
Definition: JLine1ZEstimator.hh:301
JGEOMETRY3D::JTime::getT
double getT() const
Get time.
Definition: JGeometry3D/JTime.hh:121
JTOOLS::getCosThetaC
double getCosThetaC()
Get average cosine of Cherenkov angle of water.
Definition: JConstants.hh:144
JTOOLS::j
int j
Definition: JPolint.hh:634
JGEOMETRY3D::JVersor3D
Data structure for normalised vector in three dimensions.
Definition: JVersor3D.hh:23
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
JGEOMETRY3D::JVector3D::mul
JVector3D & mul(const double factor)
Scale vector.
Definition: JVector3D.hh:173
JFIT::JLine1Z::getDistanceSquared
double getDistanceSquared(const JVector3D &pos) const
Get distance squared.
Definition: JLine1Z.hh:87
JGEOMETRY3D::JPosition3D::rotate
JPosition3D & rotate(const JRotation3D &R)
Rotate.
Definition: JPosition3D.hh:185
JGEOMETRY3D::JVector3D::negate
JVector3D & negate()
Negate vector.
Definition: JVector3D.hh:125
JGEOMETRY2D::JRotation2D::rotate
void rotate(double &__x, double &__y) const
Rotate.
Definition: JRotation2D.hh:96
JTOOLS::getInverseSpeedOfLight
const double getInverseSpeedOfLight()
Get inverse speed of light.
Definition: JConstants.hh:100
JGEOMETRY3D::JTime::__t
double __t
Definition: JGeometry3D/JTime.hh:188
JMATH::JMatrix3D::a00
double a00
Definition: JMath/JMatrix3D.hh:385
JTOOLS::getIndexOfRefraction
double getIndexOfRefraction()
Get average index of refraction of water.
Definition: JConstants.hh:111
JMATH::JMatrix3D::reset
JMatrix3D & reset()
Set matrix to the null matrix.
Definition: JMath/JMatrix3D.hh:130
JGEOMETRY3D::JPosition3D::rotate_back
JPosition3D & rotate_back(const JRotation3D &R)
Rotate back.
Definition: JPosition3D.hh:199
JGEOMETRY3D::JAngle3D::getDX
double getDX() const
Get x direction.
Definition: JAngle3D.hh:106
JGEOMETRY3D::JVersor3D::getDX
double getDX() const
Get x direction.
Definition: JVersor3D.hh:91
JFIT::JEstimator< JLine1Z >::NUMBER_OF_PARAMETERS
static const int NUMBER_OF_PARAMETERS
number of parameters of fit
Definition: JLine1ZEstimator.hh:300
JMATH::JMatrix3D::a12
double a12
Definition: JMath/JMatrix3D.hh:386
JGEOMETRY3D::JVersor3D::getDY
double getDY() const
Get y direction.
Definition: JVersor3D.hh:102
JMATH::JMatrixND_t::const_row_type
matrix_type::const_iterator const_row_type
Definition: JMatrixND.hh:46
JTOOLS::getTanThetaC
double getTanThetaC()
Get average tangent of Cherenkov angle of water.
Definition: JConstants.hh:133
JGEOMETRY3D::JTime::add
JTime & add(const JTime &value)
Addition operator.
Definition: JGeometry3D/JTime.hh:66
JGEOMETRY3D::JAngle3D
Data structure for angles in three dimensions.
Definition: JAngle3D.hh:31
JGEOMETRY3D::JAngle3D::getDZ
double getDZ() const
Get z direction.
Definition: JAngle3D.hh:128
JFIT::JLine1Z::getDot
double getDot(const JAxis3D &axis) const
Get cosine angle of impact of Cherenkov light on PMT.
Definition: JLine1Z.hh:172
JFIT::JLine1Z::JLine1Z
JLine1Z()
Default constructor.
Definition: JLine1Z.hh:39
JMATH::JMatrixND_t::const_col_type
std::vector< double >::const_iterator const_col_type
Definition: JMatrixND.hh:47
JTOOLS::getKappaC
double getKappaC()
Get average kappa of Cherenkov light in water.
Definition: JConstants.hh:166
std
Definition: jaanetDictionary.h:36
JMATH::JMatrix3D::a22
double a22
Definition: JMath/JMatrix3D.hh:387
JGEOMETRY3D::JVector3D::getY
double getY() const
Get y position.
Definition: JVector3D.hh:103
JGEOMETRY3D::JAngle3D::getDY
double getDY() const
Get y direction.
Definition: JAngle3D.hh:117
JGEOMETRY3D::JVersor3Z::getDY
double getDY() const
Get y direction.
Definition: JVersor3Z.hh:156
JFIT::JLine1Z::move
void move(const double step, const double velocity)
Move vertex along this line with given velocity.
Definition: JLine1Z.hh:62
JGEOMETRY3D::JVector3D::__z
double __z
Definition: JVector3D.hh:310
JMATH::JMatrix3D::a11
double a11
Definition: JMath/JMatrix3D.hh:386
JGEOMETRY3D::JTime::negate
JTime & negate()
Prefix unary minus.
Definition: JGeometry3D/JTime.hh:52
JMATH::JCalculator
Auxiliary class for arithmetic operations on objects.
Definition: JCalculator.hh:18
JMATH::JMatrix3D::a02
double a02
Definition: JMath/JMatrix3D.hh:385
JGEOMETRY3D::JVector3D::__x
double __x
Definition: JVector3D.hh:308
JGEOMETRY3D::JVector3D::add
JVector3D & add(const JVector3D &vector)
Add vector.
Definition: JVector3D.hh:141
JGEOMETRY3D::JVector3D::getX
double getX() const
Get x position.
Definition: JVector3D.hh:93
JMATH::JMatrix3D::a10
double a10
Definition: JMath/JMatrix3D.hh:386
JGEOMETRY3D::JVector3D::div
JVector3D & div(const double factor)
Scale vector.
Definition: JVector3D.hh:189
JTOOLS
Auxiliary classes and methods for multi-dimensional interpolations and histograms.
Definition: JAbstractCollection.hh:9
JGEOMETRY3D::JVersor3Z::getDZ
double getDZ() const
Get z direction.
Definition: JVersor3Z.hh:167
JGEOMETRY3D::JVector3D::sub
JVector3D & sub(const JVector3D &vector)
Subtract vector.
Definition: JVector3D.hh:157
JFIT::JLine1Z::getDistance
double getDistance(const JVector3D &pos) const
Get distance.
Definition: JLine1Z.hh:102
JGEOMETRY3D::JTime::div
JTime & div(const double value)
Division operator.
Definition: JGeometry3D/JTime.hh:108
JGEOMETRY3D::JTime::mul
JTime & mul(const double value)
Multiplication operator.
Definition: JGeometry3D/JTime.hh:94
JGEOMETRY3D::JQuaternion3D::rotate
void rotate(double &__x, double &__y, double &__z) const
Rotate.
Definition: JQuaternion3D.hh:645
JGEOMETRY3D::JVersor3Z::getDX
double getDX() const
Get x direction.
Definition: JVersor3Z.hh:145
JGEOMETRY3D::JTime::sub
JTime & sub(const JTime &value)
Subtraction operator.
Definition: JGeometry3D/JTime.hh:80