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.
std::vector< T >::difference_type distance(typename std::vector< T >::const_iterator first, typename PhysicsEvent::const_iterator< T > second)
Specialisation of STL distance.
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.
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.
result_type< T > operator()(T __begin, T __end)
Fit.
JKatoomba< JEstimator > estimator
JKatoomba< JGandalf > gandalf
JGlobalfit(const JGeometry &geometry, const JSoundVelocity &V, const JFitParameters ¶meters)
Constructor.
JKatoomba< JAbstractMinimiser > evaluator
JFitParameters parameters
Template specialisation of fit function of acoustic model based on JAbstractMinimiser minimiser.
Template specialisation of fit function of acoustic model based on linear approximation.
Template specialisation of fit function of acoustic model based on JGandalf minimiser.
Model for fit to acoustics data.
size_t getN() const
Get number of fit parameters.
Implementation for depth dependend velocity of sound.