1 #ifndef __JFIT__JPOINT4DESTIMATOR__
2 #define __JFIT__JPOINT4DESTIMATOR__
15 namespace JPP {
using namespace JFIT; }
98 (*this)(__begin, __end);
110 void operator ()(T __begin,T __end)
115 const int N =
distance(__begin, __end);
117 if (N >= NUMBER_OF_PARAMETERS) {
125 for (T i = __begin; i != __end; ++i) {
146 double xi =
j->getX() - getX();
147 double yi =
j->getY() - getY();
148 double zi =
j->getZ() - getZ();
151 for (
bool done =
false; !done; ) {
153 if ((done = (++
j == __end))) {
157 double xj =
j->getX() - getX();
158 double yj =
j->getY() - getY();
159 double zj =
j->getZ() - getZ();
167 const double y = ((xj + xi) * dx +
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;
219 throw JValueOutOfRange(
"JEstimator<JPoint4D>::operator(): Not enough data points.");
223 static const int NUMBER_OF_PARAMETERS = 4;
std::vector< T >::difference_type distance(typename std::vector< T >::const_iterator first, typename PhysicsEvent::const_iterator< T > second)
Specialisation of STL distance.
JEstimator()
Fit constructor.
JMATH::JMatrix4S V
co-variance matrix of fit parameters
JEstimator(T __begin, T __end)
Fit constructor.
Template definition of linear fit.
Data structure for vertex fit.
Exception for accessing a value in a collection that is outside of its range.
Auxiliary classes and methods for linear and iterative data regression.
double getIndexOfRefraction()
Get average index of refraction of water corresponding to group velocity.
const double getInverseSpeedOfLight()
Get inverse speed of light.
const double getSpeedOfLight()
Get speed of light.
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).