JACOUSTICS::JKatoomba_t Struct Reference

Worker class for complete fit procedure of acoustic model. More...

#include <JKatoomba_t.hh>

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

Public Member Functions
	JKatoomba_t (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

Private Attributes
JModel	result

Detailed Description

Worker class for complete fit procedure of acoustic model.

Definition at line 29 of file JKatoomba_t.hh.

Constructor & Destructor Documentation

JACOUSTICS::JKatoomba_t::JKatoomba_t ( 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 JKatoomba_t.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::JKatoomba_t::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 JKatoomba_t.hh.

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

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

        for (double chi2 = evaluator(result, __begin, __end), chi2_old = chi2; ; ) {

          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) {

            iter_swap(out, --__end);

            result = estimator(__begin, __end);
            chi2   = evaluator(result, __begin, __end);

            if (chi2_old - chi2  >  parameters.stdev * parameters.stdev) {

              chi2_old = chi2;

              continue;

            } else {

              iter_swap(out, __end++);

              result = estimator(__begin, __end);
              chi2   = chi2_old;
            }
          }

          break;
        }
      }

      gandalf.value = result;                                         // start value

      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::JKatoomba_t::parameters

Definition at line 154 of file JKatoomba_t.hh.

JKatoomba<JEstimator> JACOUSTICS::JKatoomba_t::estimator

Definition at line 155 of file JKatoomba_t.hh.

JKatoomba<JAbstractMinimiser> JACOUSTICS::JKatoomba_t::evaluator

Definition at line 156 of file JKatoomba_t.hh.

JKatoomba<JGandalf> JACOUSTICS::JKatoomba_t::gandalf

Definition at line 157 of file JKatoomba_t.hh.

JModel JACOUSTICS::JKatoomba_t::result

private

Definition at line 160 of file JKatoomba_t.hh.

---

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

• JKatoomba_t.hh