1#ifndef __JACOUSTICS__JGLOBALFIT__
2#define __JACOUSTICS__JGLOBALFIT__
108 double chi2[] = { numeric_limits<double>::max(),
109 numeric_limits<double>::max() };
117 for (T hit = __begin; hit != __end; ++hit) {
119 const double x = fabs(hit->getValue() -
estimator.getToA(
result, *hit)) / hit->getSigma();
129 if (chi2[0] == numeric_limits<double>::max()) {
135 iter_swap(out, --__end);
148 iter_swap(out, __end++);
158 const double chi2 =
gandalf (__begin, __end) /
gandalf.estimator->getRho(1.0);
Model for fit to acoutsics data.
Fit functions of acoustic model.
Auxiliary classes and methods for acoustic position calibration.
double getWeight(T __begin, T __end)
Get total weight of data points.
JMEstimator * getMEstimator(const int type)
Get M-Estimator.
void model(JModel_t &value)
Auxiliary function to constrain model during fit.
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
int mestimator
M-estimator.
double stdev
standard deviation for outlier removal
Auxiliary data structure for return value of fit.
JTimeRange getTimeRange() const
Get time range of data.
int size() const
Get number of hits used in fit (after outlier removal).
Global fit of prameterised detector geometry to acoustics data.
JKatoomba< JEstimator > estimator
JKatoomba< JGandalf > gandalf
JGlobalfit(const JGeometry &geometry, const JSoundVelocity &V, const JFitParameters ¶meters)
Constructor.
JKatoomba< JAbstractMinimiser > evaluator
result_type< T > operator()(T __begin, T __end)
Fit.
JFitParameters parameters
Template definition of fit function of acoustic model.
Model for fit to acoustics data.
size_t getN() const
Get number of fit parameters.
Implementation for depth dependend velocity of sound.