Jpp  17.1.1
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< 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< JFirst_t, JSecond_t > JGEOMETRY3D::JVector3D JMATH::JMath< JFirst_t, JSecond_t >

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...
 
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...
 
JFirst_t & mul (const JSecond_t &object)
 Multiply with object. More...
 
JFirst_t & mul (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...
 
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...
 
JVector3DgetCross (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

Default constructor.

Definition at line 96 of file JLine1ZEstimator.hh.

96  :
97  JLine1Z()
98  {}
JLine1Z()
Default constructor.
Definition: JLine1Z.hh:39
template<class T >
JFIT::JEstimator< JLine1Z >::JEstimator ( T  __begin,
T  __end 
)
inline

Constructor.

Parameters
__beginbegin of data
__endend of data

Definition at line 108 of file JLine1ZEstimator.hh.

108  :
109  JLine1Z()
110  {
111  (*this)(__begin, __end);
112  }
JLine1Z()
Default constructor.
Definition: JLine1Z.hh:39

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

Definition at line 128 of file JLine1ZEstimator.hh.

129  {
130  using namespace std;
131  using namespace JPP;
132 
133  const int N = distance(__begin, __end);
134 
135  if (N >= NUMBER_OF_PARAMETERS) {
136 
137  __x = 0.0;
138  __y = 0.0;
139  __z = 0.0;
140  __t = 0.0;
141 
142  double t0 = 0.0;
143 
144  for (T i = __begin; i != __end; ++i) {
145  __x += i->getX();
146  __y += i->getY();
147  __z += i->getZ();
148  t0 += i->getT();
149  }
150 
151  const double W = 1.0/N;
152 
153  __x *= W;
154  __y *= W;
155  __z *= W;
156  t0 *= W;
157 
158  V.reset();
159 
160  t0 *= getSpeedOfLight();
161 
162  double y0 = 0.0;
163  double y1 = 0.0;
164  double y2 = 0.0;
165 
166  T j = __begin;
167 
168  double xi = j->getX() - getX();
169  double yi = j->getY() - getY();
170  double ti = (j->getT() * getSpeedOfLight() - t0 - j->getZ() + getZ()) / getKappaC();
171 
172  for (bool done = false; !done; ) {
173 
174  if ((done = (++j == __end))) {
175  j = __begin;
176  }
177 
178  double xj = j->getX() - getX();
179  double yj = j->getY() - getY();
180  double tj = (j->getT() * getSpeedOfLight() - t0 - j->getZ() + getZ()) / getKappaC();
181 
182  double dx = xj - xi;
183  double dy = yj - yi;
184  double dt = ti - tj; // opposite sign!
185 
186  const double y = ((xj + xi) * dx +
187  (yj + yi) * dy +
188  (tj + ti) * dt);
189 
190  dx *= 2;
191  dy *= 2;
192  dt *= 2;
193 
194  V.a00 += dx * dx;
195  V.a01 += dx * dy;
196  V.a02 += dx * dt;
197  V.a11 += dy * dy;
198  V.a12 += dy * dt;
199  V.a22 += dt * dt;
200 
201  y0 += dx * y;
202  y1 += dy * y;
203  y2 += dt * y;
204 
205  xi = xj;
206  yi = yj;
207  ti = tj;
208  }
209 
210  t0 *= getInverseSpeedOfLight();
211 
212  V.a10 = V.a01;
213  V.a20 = V.a02;
214  V.a21 = V.a12;
215 
216  svd.decompose(V);
217 
218  if (fabs(svd.S.a11) < MINIMAL_SVD_WEIGHT * fabs(svd.S.a00)) {
219  THROW(JValueOutOfRange, "JEstimator<JLine1Z>::JEstimator(): singular value " << svd.S.a11 << ' ' << svd.S.a00);
220  }
221 
223 
224  __x += V.a00 * y0 + V.a01 * y1 + V.a02 * y2;
225  __y += V.a10 * y0 + V.a11 * y1 + V.a12 * y2;
226  __t = V.a20 * y0 + V.a21 * y1 + V.a22 * y2;
227 
229  __t += t0;
230 
231  } else {
232  throw JValueOutOfRange("JEstimator<JLine1Z>::JEstimator(): Not enough data points.");
233  }
234 
235  return *this;
236  }
double getKappaC()
Get average R-dependence of arrival time of Cherenkov light (a.k.a.
std::vector< T >::difference_type distance(typename std::vector< T >::const_iterator first, typename PhysicsEvent::const_iterator< T > second)
Specialisation of STL distance.
static const int NUMBER_OF_PARAMETERS
number of parameters of fit
#define THROW(JException_t, A)
Marco for throwing exception with std::ostream compatible message.
Definition: JException.hh:696
then JShowerPostfit f $INPUT_FILE o $OUTPUT_FILE N
JMatrix3D & reset()
Set matrix to the null matrix.
do set_variable OUTPUT_DIRECTORY $WORKDIR T
void decompose(const JMatrix3D &A)
Decompose given matrix.
Definition: JSVD3D.hh:52
JMatrix3D S
Definition: JSVD3D.hh:223
double getY() const
Get y position.
Definition: JVector3D.hh:104
JMATH::JMatrix3S V
co-variance matrix of fit parameters
const double getSpeedOfLight()
Get speed of light.
const JMatrix3D & invert(const double precision=1.0e-12) const
Get inverted matrix.
Definition: JSVD3D.hh:192
const double getInverseSpeedOfLight()
Get inverse speed of light.
double getX() const
Get x position.
Definition: JVector3D.hh:94
int j
Definition: JPolint.hh:703
static double MINIMAL_SVD_WEIGHT
minimal SVD weight.
double getZ() const
Get z position.
Definition: JVector3D.hh:115
esac done
Definition: JAddHDE.sh:21
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
__beginbegin of data
__endend of data
Aco-variance matrix of hits

Definition at line 256 of file JLine1ZEstimator.hh.

257  {
258  using namespace std;
259  using namespace JPP;
260 
261  const int N = distance(__begin, __end);
262 
263  if (N != (int) A.size()) {
264  THROW(JValueOutOfRange, "JEstimator<JLine1Z>::update(): Wrong number of points " << N << " != " << A.size());
265  }
266 
267  if (N >= NUMBER_OF_PARAMETERS) {
268 
269  double x1 = 0.0;
270  double y1 = 0.0;
271  double t1 = 0.0;
272 
273  V.reset();
274 
275  T i = __begin;
276 
277  for (size_t row = 0; row != A.size(); ++row, ++i) {
278 
279  const double dx = i->getX() - getX();
280  const double dy = i->getY() - getY();
281 
282  const double rt = sqrt(dx*dx + dy*dy);
283 
284  double xr = getKappaC() * getInverseSpeedOfLight();
285  double yr = getKappaC() * getInverseSpeedOfLight();
286  double tr = 1.0;
287 
288  if (rt != 0.0) {
289  xr *= -dx / rt;
290  yr *= -dy / rt;
291  }
292 
293  T j = __begin;
294 
295  for (size_t col = 0; col != A.size(); ++col, ++j) {
296 
297  const double dx = j->getX() - getX();
298  const double dy = j->getY() - getY();
299  const double dz = j->getZ() - getZ();
300 
301  const double rt = sqrt(dx*dx + dy*dy);
302 
303  double xc = getKappaC() * getInverseSpeedOfLight();
304  double yc = getKappaC() * getInverseSpeedOfLight();
305  double tc = 1.0;
306 
307  if (rt != 0.0) {
308  xc *= -dx / rt;
309  yc *= -dy / rt;
310  }
311 
312  const double ts = j->getT() - (dz + rt * getKappaC()) * getInverseSpeedOfLight();
313 
314  const double vs = A(row,col);
315 
316  x1 += xr * vs * ts;
317  y1 += yr * vs * ts;
318  t1 += tr * vs * ts;
319 
320  V.a00 += xr * vs * xc;
321  V.a01 += xr * vs * yc;
322  V.a02 += xr * vs * tc;
323  V.a11 += yr * vs * yc;
324  V.a12 += yr * vs * tc;
325  V.a22 += tr * vs * tc;
326  }
327  }
328 
329  V.a10 = V.a01;
330  V.a20 = V.a02;
331  V.a21 = V.a12;
332 
333  svd.decompose(V);
334 
335  if (fabs(svd.S.a11) < MINIMAL_SVD_WEIGHT * fabs(svd.S.a00)) {
336  THROW(JValueOutOfRange, "JEstimator<JLine1Z>::update(): singular value " << svd.S.a11 << ' ' << svd.S.a00);
337  }
338 
340 
341  __x += V.a00 * x1 + V.a01 * y1 + V.a02 * t1;
342  __y += V.a10 * x1 + V.a11 * y1 + V.a12 * t1;
343  __t = V.a20 * x1 + V.a21 * y1 + V.a22 * t1;
344 
345  } else {
346  THROW(JValueOutOfRange, "JEstimator<JLine1Z>::update(): Not enough data points " << N);
347  }
348  }
double getKappaC()
Get average R-dependence of arrival time of Cherenkov light (a.k.a.
std::vector< T >::difference_type distance(typename std::vector< T >::const_iterator first, typename PhysicsEvent::const_iterator< T > second)
Specialisation of STL distance.
static const int NUMBER_OF_PARAMETERS
number of parameters of fit
#define THROW(JException_t, A)
Marco for throwing exception with std::ostream compatible message.
Definition: JException.hh:696
then JShowerPostfit f $INPUT_FILE o $OUTPUT_FILE N
JMatrix3D & reset()
Set matrix to the null matrix.
do set_variable OUTPUT_DIRECTORY $WORKDIR T
void decompose(const JMatrix3D &A)
Decompose given matrix.
Definition: JSVD3D.hh:52
JMatrix3D S
Definition: JSVD3D.hh:223
double getY() const
Get y position.
Definition: JVector3D.hh:104
JMATH::JMatrix3S V
co-variance matrix of fit parameters
size_t size() const
Get dimension of matrix.
Definition: JMatrixND.hh:181
const JMatrix3D & invert(const double precision=1.0e-12) const
Get inverted matrix.
Definition: JSVD3D.hh:192
const double getInverseSpeedOfLight()
Get inverse speed of light.
double getX() const
Get x position.
Definition: JVector3D.hh:94
int j
Definition: JPolint.hh:703
source $JPP_DIR setenv csh $JPP_DIR &dev null eval JShellParser o a A
static double MINIMAL_SVD_WEIGHT
minimal SVD weight.
double getZ() const
Get z position.
Definition: JVector3D.hh:115
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  }
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  }
void move(const double step, const double velocity)
Move vertex along this line with given velocity.
Definition: JLine1Z.hh:62
double getZ() const
Get z position.
Definition: JVector3D.hh:115
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  }
double getY() const
Get y position.
Definition: JVector3D.hh:104
double getX() const
Get x position.
Definition: JVector3D.hh:94
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  }
double getDistanceSquared(const JVector3D &pos) const
Get distance squared.
Definition: JLine1Z.hh:87
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 JPP;
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  }
double getT() const
Get time.
double getKappaC()
Get average R-dependence of arrival time of Cherenkov light (a.k.a.
double getY() const
Get y position.
Definition: JVector3D.hh:104
then usage $script[distance] fi case set_variable R
Definition: JDrawLED.sh:43
const double getInverseSpeedOfLight()
Get inverse speed of light.
double getX() const
Get x position.
Definition: JVector3D.hh:94
double getZ() const
Get z position.
Definition: JVector3D.hh:115
double JGEOMETRY3D::JTime::getT ( ) const
inlineinherited

Get time.

Returns
time

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

123  {
124  return __t;
125  }
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 JPP;
137 
138  return pos.getZ() - this->getDistance(pos) / getTanThetaC();
139  }
double getDistance(const JVector3D &pos) const
Get distance.
Definition: JLine1Z.hh:102
double getTanThetaC()
Get average tangent of Cherenkov angle of water corresponding to group velocity.
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  }
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 JPP;
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  }
double getCosThetaC()
Get average cosine of Cherenkov angle of water corresponding to group velocity.
double getY() const
Get y position.
Definition: JVector3D.hh:104
then usage $script[distance] fi case set_variable R
Definition: JDrawLED.sh:43
double getX() const
Get x position.
Definition: JVector3D.hh:94
double getSinThetaC()
Get average sine of Cherenkov angle of water corresponding to group velocity.
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  }
double getDot(const JAxis3D &axis) const
Get cosine angle of impact of Cherenkov light on PMT.
Definition: JLine1Z.hh:172
JVersor3D getDirection(const JVector3D &pos) const
Get photon direction of Cherenkov light on PMT.
Definition: JLine1Z.hh:148
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  }
double getY() const
Get y position.
Definition: JVector3D.hh:104
double getX() const
Get x position.
Definition: JVector3D.hh:94
double getZ() const
Get z position.
Definition: JVector3D.hh:115
double getDX() const
Get x direction.
Definition: JAngle3D.hh:108
double getDZ() const
Get z direction.
Definition: JAngle3D.hh:130
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  }
double getDY() const
Get y direction.
Definition: JVersor3D.hh:106
double getDX() const
Get x direction.
Definition: JVersor3D.hh:95
double getY() const
Get y position.
Definition: JVector3D.hh:104
double getX() const
Get x position.
Definition: JVector3D.hh:94
double getZ() const
Get z position.
Definition: JVector3D.hh:115
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  }
double getDY() const
Get y direction.
Definition: JVersor3Z.hh:158
double getY() const
Get y position.
Definition: JVector3D.hh:104
double getDX() const
Get x direction.
Definition: JVersor3Z.hh:147
double getX() const
Get x position.
Definition: JVector3D.hh:94
double getDZ() const
Get z direction.
Definition: JVersor3Z.hh:169
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  }
double getY() const
Get y position.
Definition: JVector3D.hh:104
double getX() const
Get x position.
Definition: JVector3D.hh:94
double getZ() const
Get z position.
Definition: JVector3D.hh:115
static parameter_type JFIT::JLine1Z::pX ( )
inlinestaticinherited

Definition at line 180 of file JLine1Z.hh.

180 { return &JLine1Z::__x; }
static parameter_type JFIT::JLine1Z::pY ( )
inlinestaticinherited

Definition at line 181 of file JLine1Z.hh.

181 { return &JLine1Z::__y; }
static parameter_type JFIT::JLine1Z::pT ( )
inlinestaticinherited

Definition at line 182 of file JLine1Z.hh.

182 { return &JLine1Z::__t; }
JVertex3D& JGEOMETRY3D::JVertex3D::negate ( )
inlineinherited

Prefix unary minus.

Returns
line

Definition at line 73 of file JVertex3D.hh.

74  {
77 
78  return *this;
79  }
JTime & negate()
Prefix unary minus.
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  }
JTime & add(const JTime &value)
Addition operator.
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  }
double getY() const
Get y position.
Definition: JVector3D.hh:104
double getX() const
Get x position.
Definition: JVector3D.hh:94
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  }
double getT() const
Get time.
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  }
JTime & sub(const JTime &value)
Subtraction operator.
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  }
double getY() const
Get y position.
Definition: JVector3D.hh:104
double getX() const
Get x position.
Definition: JVector3D.hh:94
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  }
double getT() const
Get time.
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  }
JVector3D & mul(const double factor)
Scale vector.
Definition: JVector3D.hh:174
JTime & mul(const double value)
Multiplication operator.
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 357 of file JMath.hh.

358  {
359  return static_cast<JFirst_t&>(*this) = JCalculator<JFirst_t>::calculator.mul(static_cast<const JFirst_t&>(*this), object);
360  }
Auxiliary class for arithmetic operations on objects.
Definition: JCalculator.hh:18
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 357 of file JMath.hh.

358  {
359  return static_cast<JFirst_t&>(*this) = JCalculator<JFirst_t>::calculator.mul(static_cast<const JFirst_t&>(*this), object);
360  }
Auxiliary class for arithmetic operations on objects.
Definition: JCalculator.hh:18
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  }
JTime & div(const double value)
Division operator.
JVector3D & div(const double factor)
Scale vector.
Definition: JVector3D.hh:190
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  }
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  }
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  }
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  }
Data structure for angles in three dimensions.
Definition: JAngle3D.hh:33
double getY() const
Get y position.
Definition: JVector3D.hh:104
double getX() const
Get x position.
Definition: JVector3D.hh:94
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  }
double getY() const
Get y position.
Definition: JVector3D.hh:104
double getX() const
Get x position.
Definition: JVector3D.hh:94
Data structure for normalised vector in three dimensions.
Definition: JVersor3D.hh:26
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  }
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  }
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  }
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  }
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  }
void rotate(double &__x, double &__y, double &__z) const
Rotate.
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  }
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  }
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  }
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  }
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  }
void rotate_back(double &__x, double &__y, double &__z) const
Rotate back.
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  }
JVector3D & sub(const JVector3D &vector)
Subtract vector.
Definition: JVector3D.hh:158
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  }
do set_variable OUTPUT_DIRECTORY $WORKDIR T
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  }
JPosition3D & rotate_back(const JRotation3D &R)
Rotate back.
Definition: JPosition3D.hh:200
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  }
double getY() const
Get y position.
Definition: JVector3D.hh:104
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  }
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  }
double getY() const
Get y position.
Definition: JVector3D.hh:104
double getX() const
Get x position.
Definition: JVector3D.hh:94
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  }
double getY() const
Get y position.
Definition: JVector3D.hh:104
double getX() const
Get x position.
Definition: JVector3D.hh:94
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  }
double getLengthSquared() const
Get length squared.
Definition: JVector3D.hh:235
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  }
double getY() const
Get y position.
Definition: JVector3D.hh:104
double getX() const
Get x position.
Definition: JVector3D.hh:94
double getZ() const
Get z position.
Definition: JVector3D.hh:115

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.

co-variance matrix of fit parameters

Definition at line 352 of file JLine1ZEstimator.hh.

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: