JACOUSTICS::JGlobalfit Struct Reference

Global fit of prameterised detector geometry to acoustics data. More...

#include <JGlobalfit.hh>

Classes
struct	result_type
	Auxiliary data structure for return value of fit. More...

Public Member Functions
	JGlobalfit (const JGeometry &geometry, const JSoundVelocity &V, const JFitParameters &parameters)
	Constructor. More...
template<class T >
result_type< T >	operator() (T __begin, T __end)
	Fit. More...

Public Attributes
JFitParameters	parameters
JKatoomba< JEstimator >	estimator
JKatoomba< JAbstractMinimiser >	evaluator
JKatoomba< JGandalf >	gandalf

Detailed Description

Global fit of prameterised detector geometry to acoustics data.

Definition at line 29 of file JGlobalfit.hh.

Constructor & Destructor Documentation

JACOUSTICS::JGlobalfit::JGlobalfit ( const JGeometry &  geometry, const JSoundVelocity &  V, const JFitParameters &  parameters  )

inline

Constructor.

Parameters

geometry	detector geometry
V	sound velocity
parameters	parameters

Definition at line 37 of file JGlobalfit.hh.

                                                                                                     :
      parameters(parameters),
      estimator(geometry, V, parameters.option),
      evaluator(geometry, V, parameters.option),
      gandalf  (geometry, V, parameters.option)
    {
      using namespace JPP;

      evaluator.estimator.reset(getMEstimator(parameters.mestimator));
      gandalf  .estimator.reset(getMEstimator(parameters.mestimator));
    }

Member Function Documentation

template<class T >

result_type<T> JACOUSTICS::JGlobalfit::operator() ( T  __begin, T  __end  )

inline

Fit.

Parameters

__begin	begin of data
__end	end of data

Returns

fit result

Definition at line 97 of file JGlobalfit.hh.

    {
      using namespace std;

      JModel  model;
      JModel& result = gandalf.value;                                 // start value

      result = estimator(__begin, __end);                             // pre-fit

      if (parameters.stdev > 0.0) {                                   // remove outliers

        double chi2[] = { numeric_limits<double>::max(),
                          numeric_limits<double>::max() };

        for ( ; ; ) {

          T out = __end;

          double xmax = 0.0;

          for (T hit = __begin; hit != __end; ++hit) {

            const double x = fabs(hit->getValue() - estimator.getToA(result, *hit)) / hit->getSigma();

            if (x > xmax) {
              xmax = x;
              out  = hit;
            }
          }

          if (xmax > parameters.stdev) {

            if (chi2[0] == numeric_limits<double>::max()) {
              chi2[0] = evaluator(result, __begin, __end);
            }

            model = result;

            iter_swap(out, --__end);

            result  = estimator(__begin, __end);
            chi2[1] = evaluator(result, __begin, __end);

            if (chi2[0] - chi2[1]  >  parameters.stdev * parameters.stdev) {

              chi2[0] = chi2[1];

              continue;

            } else {

              iter_swap(out, __end++);

              result = model;
            }
          }

          break;
        }
      }

      const double chi2  =  gandalf  (__begin, __end)  /  gandalf.estimator->getRho(1.0);
      const double ndf   =  getWeight(__begin, __end)  -  gandalf.value.getN();

      return { gandalf.value, chi2, ndf, __begin, __end };            // return values
    }

Member Data Documentation

JFitParameters JACOUSTICS::JGlobalfit::parameters

Definition at line 165 of file JGlobalfit.hh.

JKatoomba<JEstimator> JACOUSTICS::JGlobalfit::estimator

Definition at line 166 of file JGlobalfit.hh.

JKatoomba<JAbstractMinimiser> JACOUSTICS::JGlobalfit::evaluator

Definition at line 167 of file JGlobalfit.hh.

JKatoomba<JGandalf> JACOUSTICS::JGlobalfit::gandalf

Definition at line 168 of file JGlobalfit.hh.

---

The documentation for this struct was generated from the following file:

• JGlobalfit.hh