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.
template<class T >
result_type< T >	operator() (T __begin, T __end)
	Fit.

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

JGlobalfit()

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

operator()()

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

parameters

JFitParameters JACOUSTICS::JGlobalfit::parameters

Definition at line 165 of file JGlobalfit.hh.

estimator

JKatoomba<JEstimator> JACOUSTICS::JGlobalfit::estimator

Definition at line 166 of file JGlobalfit.hh.

evaluator

JKatoomba<JAbstractMinimiser> JACOUSTICS::JGlobalfit::evaluator

Definition at line 167 of file JGlobalfit.hh.

gandalf

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