JFIT::JRegressor< JPoint4E, JGandalf > Struct Reference

Regressor function object for JPoint4E fit using JGandalf minimiser. More...

#include <JShowerBrightPointRegressor.hh>

Inheritance diagram for JFIT::JRegressor< JPoint4E, JGandalf >:

Public Types
typedef JTOOLS::JSplineFunction1S_t	JFunction1D_t
typedef JTOOLS::JMAPLIST< JTOOLS::JPolint2FunctionalMap, JTOOLS::JPolint1FunctionalGridMap >::maplist	JPDFMapList_t
typedef JPHYSICS::JPDFTable< JFunction1D_t, JPDFMapList_t >	JPDF_t
typedef JGandalf< JPoint4E >	minimiser_type
typedef JRegressor< JPoint4E, JGandalf >	regressor_type
typedef minimiser_type::result_type	result_type
typedef JFIT_LOCAL::JTypedef_t< JModel_t >::parameter_type	parameter_type
	Data type of fit parameter. More...

Public Member Functions
	JRegressor (const std::string &fileDescriptor, const double TTS, const int numberOfPoints=25, const double epsilon=1.0e-10)
	Parameterized constructor. More...
template<class JHit_t >
result_type	operator() (const JPoint4E &vx, const JHit_t &hit) const
	Fit function. More...
JPDF_t::result_type	getH0 (const double R_Hz, const double t1) const
	Get background hypothesis value for time differentiated PDF. More...
JPDF_t::result_type	getH1 (const double D, const double ct, const double t) const
	Get signal hypothesis value per 1 GeV for bright point emission PDF. More...
double	getRmax () const
	Get maximal road width of PDF. More...
result_type	operator() (const JPoint4E &value, T __begin, T __end)
	Global fit. More...
template<class JFunction_t , class T , class ... Args>
result_type	operator() (const JFunction_t &fit, T __begin, T __end, Args ...args)
	Multi-dimensional fit of multiple data sets. More...

Public Attributes
JPDF_t	pdf [NUMBER_OF_PDFS]
	PDF. More...
std::vector< parameter_type >	parameters
	fit parameters More...
int	numberOfIterations
	number of iterations More...
double	lambda
	control parameter More...
JModel_t	value
	value More...
JModel_t	error
	error More...
JMATH::JMatrixNS	V
	Hesse matrix. More...

Static Public Attributes
static JTimeRange	T_ns
	Time window with respect to Cherenkov hypothesis [ns]. More...
static double	Vmax_npe = std::numeric_limits<double>::max()
	Maximal integral of PDF [npe]. More...
static const int	NUMBER_OF_PDFS = 2
static const JPDFType_t	pdf_t [NUMBER_OF_PDFS]
	PDF types. More...
static int	MAXIMUM_ITERATIONS = 1000
	maximal number of iterations More...
static double	EPSILON = 1.0e-3
	maximal distance to minimum More...
static bool	EPSILON_ABSOLUTE = false
	set epsilon to absolute difference instead of relative More...
static double	LAMBDA_MIN = 0.01
	minimal value control parameter More...
static double	LAMBDA_MAX = 100.0
	maximal value control parameter More...
static double	LAMBDA_UP = 10.0
	multiplication factor control parameter More...
static double	LAMBDA_DOWN = 10.0
	multiplication factor control parameter More...
static double	PIVOT = std::numeric_limits<double>::epsilon()
	minimal value diagonal element of Hesse matrix More...
static int	debug = 0
	debug level (default is off). More...

Private Member Functions
void	reset ()
	Reset current parameters. More...
template<class JFunction_t , class T , class ... Args>
void	update (const JFunction_t &fit, T __begin, T __end, Args ...args)
	Recursive method to update current parameters. More...
template<class JFunction_t >
void	update (const JFunction_t &fit)
	Termination method to update current parameters. More...

Static Private Member Functions
static double	get (const JModel_t &model, double JModel_t::*parameter)
	Read/write access to parameter value by data member. More...
static double &	get (JModel_t &model, double JModel_t::*parameter)
	Read/write access to parameter value by data member. More...
static double	get (const JModel_t &model, const size_t index)
	Read/write access to parameter value by index. More...
static double &	get (JModel_t &model, const size_t index)
	Read/write access to parameter value by index. More...
static double	get (const JModel_t &model, const int index)
	Read/write access to parameter value by index. More...
static double &	get (JModel_t &model, const int index)
	Read/write access to parameter value by index. More...

Private Attributes
std::vector< double >	h
JMATH::JVectorND	x
struct {
result_type   result
}	current
struct {
JModel_t   value
result_type   result
}	previous

Detailed Description

Regressor function object for JPoint4E fit using JGandalf minimiser.

Definition at line 44 of file JShowerBrightPointRegressor.hh.

Member Typedef Documentation

JFunction1D_t

typedef JTOOLS::JSplineFunction1S_t JFIT::JRegressor< JPoint4E, JGandalf >::JFunction1D_t

Definition at line 49 of file JShowerBrightPointRegressor.hh.

JPDFMapList_t

typedef JTOOLS::JMAPLIST<JTOOLS::JPolint2FunctionalMap, JTOOLS::JPolint1FunctionalGridMap>::maplist JFIT::JRegressor< JPoint4E, JGandalf >::JPDFMapList_t

Definition at line 51 of file JShowerBrightPointRegressor.hh.

JPDF_t

typedef JPHYSICS::JPDFTable<JFunction1D_t, JPDFMapList_t> JFIT::JRegressor< JPoint4E, JGandalf >::JPDF_t

Definition at line 52 of file JShowerBrightPointRegressor.hh.

minimiser_type

typedef JGandalf <JPoint4E > JFIT::JAbstractRegressor< JPoint4E , JGandalf >::minimiser_type

inherited

Definition at line 80 of file JRegressor.hh.

regressor_type

typedef JRegressor<JPoint4E , JGandalf > JFIT::JAbstractRegressor< JPoint4E , JGandalf >::regressor_type

inherited

Definition at line 81 of file JRegressor.hh.

result_type

typedef minimiser_type::result_type JFIT::JAbstractRegressor< JPoint4E , JGandalf >::result_type

inherited

Definition at line 82 of file JRegressor.hh.

parameter_type

template<class JModel_t >

typedef JFIT_LOCAL::JTypedef_t<JModel_t>::parameter_type JFIT::JGandalf< JModel_t >::parameter_type

inherited

Data type of fit parameter.

Definition at line 95 of file JGandalf.hh.

Constructor & Destructor Documentation

JRegressor()

JFIT::JRegressor< JPoint4E, JGandalf >::JRegressor ( const std::string &  fileDescriptor, const double  TTS, const int  numberOfPoints = 25, const double  epsilon = 1.0e-10  )

inline

Parameterized constructor.

The PDF file descriptor should contain the wild card character JPHYSICS::WILDCARD which will be replaced by the corresponding PDF types.

Parameters

fileDescriptor	PDF file descriptor
TTS	TTS [ns]
numberOfPoints	number of points for Gauss-Hermite integration of TTS
epsilon	precision for Gauss-Hermite integration of TTS

Definition at line 66 of file JShowerBrightPointRegressor.hh.

    {
      using namespace std;
      using namespace JPP;
      
      const JPDF_t::JSupervisor supervisor(new JPDF_t::JDefaultResult(JMATH::zero));
      
      for (int i = 0; i != NUMBER_OF_PDFS; ++i) {
        
        try {
 
          const string file_name = getFilename(fileDescriptor, pdf_t[i]);
          
          NOTICE("loading PDF from file " << file_name << "... " << flush);
 
          pdf[i].load(file_name.c_str());
 
          NOTICE("OK" << endl);
 
          pdf[i].setExceptionHandler(supervisor);
        }
        catch(const JException& error) {
          FATAL(error.what() << endl);
        }
      }
 
      // Add PDFs
      for (int i = 1; i < NUMBER_OF_PDFS; i += 2) {
  
        pdf[ i ].add(pdf[i-1]); 
 
        JPDF_t buffer;
 
        pdf[i-1].swap(buffer);
        
        if        (TTS > 0.0) {
          pdf[i].blur(TTS, numberOfPoints, epsilon);
        } else if (TTS < 0.0) {
          ERROR("Illegal value of TTS [ns]: " << TTS << endl);
        }
      }    
    }

Member Function Documentation

operator()() [1/3]

template<class JHit_t >

result_type JFIT::JRegressor< JPoint4E, JGandalf >::operator() ( const JPoint4E &  vx, const JHit_t &  hit  ) const

inline

Fit function.

This method is used to determine the chi2 and gradient of given hit with respect a bright point emitting isotropically

JHit_t refers to a data structure which should have the following member methods:

• double getX(); // [m]
• double getY(); // [m]
• double getZ(); // [m]
• double getDX(); // [u]
• double getDY(); // [u]
• double getDZ(); // [u]
• double getT(); // [ns]

Parameters

vx	shower vertex
hit	hit

Returns

chi2 and gradient

Definition at line 130 of file JShowerBrightPointRegressor.hh.

    {      
      using namespace JPP;
 
      JPosition3D  D(hit.getPosition());
      JDirection3D U(hit.getDirection());
 
      D.sub(vx.getPosition());
      double length = D.getLength();
      double ct = U.getDot(D) / length;
 
      if (ct > +1.0) { ct = +1.0; }
      if (ct < -1.0) { ct = -1.0; }
 
      const double t  = vx.getT() + (length * getIndexOfRefraction() * getInverseSpeedOfLight());
 
      const double dt = T_ns.constrain(hit.getT() - t);
 
      JPDF_t::result_type H0 = getH0(hit.getR(), dt);                    // getH0 = Get background hypothesis value 
      JPDF_t::result_type H1 = getH1(length, ct, dt);  // getH1 = Get signal hypothesis value / 1 GeV
 
      if (get_value(H1) >= Vmax_npe) {
        H1 *= Vmax_npe / get_value(H1);
      }
 
      double H1_value = get_value(H1);
      double v_H1 = H1.v; //Integral from tmin to t of just H1
      double V_H1 = H1.V; //Integral from tmin to tmax of just H1
      H1 *= vx.getE();
      
      JPDF_t::result_type HT = H1+H0; //now H1 is signal + background
      double HT_value = get_value(HT);
      result_type result;
      result.chi2 = HT.getChi2() - H0.getChi2();  // Likelihood ratio
      
      double exp_V_HT = exp(-HT.V); //V is the integral from tmin to tmax of EH1+H0
 
      double energy_gradient = -1 / HT_value; //dPdE
      energy_gradient *= (H1_value -  HT_value * v_H1) * (1-exp_V_HT)  -  HT_value * exp_V_HT * V_H1; //Numerator
      energy_gradient /= (1-exp_V_HT); // Denominator
 
      /*
       * Here it is evaluated: d(chi2)/d(ct) * d(ct)/d(x0,y0,z0,t0)  + d(chi2)/dE
       */
      result.gradient = JPoint4E(JPoint4D(JVector3D(-getIndexOfRefraction() * D.getX() / length,    // d(ct)/d(x0) 
                                                    -getIndexOfRefraction() * D.getY() / length,    // d(ct)/d(y0) 
                                                    -getIndexOfRefraction() * D.getZ() / length),   // d(ct)/d(z0) 
                                          getSpeedOfLight()),                                               // d(ct)/d(t0)
                                 energy_gradient);                                                              // d(chi2)/d(E)
      
      static_cast<JPoint4D&>(result.gradient).mul(getInverseSpeedOfLight() * (HT.getDerivativeOfChi2() -
                                                                              H0.getDerivativeOfChi2()));                 // x d(chi2)/d(ct1)
 
      return result;
 
    }

getH0()

JPDF_t::result_type JFIT::JRegressor< JPoint4E, JGandalf >::getH0 ( const double  R_Hz, const double  t1  ) const

inline

Get background hypothesis value for time differentiated PDF.

Parameters

R_Hz	rate [Hz]
t1	time [ns]

Returns

hypothesis value

Definition at line 194 of file JShowerBrightPointRegressor.hh.

    {
      using namespace JPP;
 
      return JPDF_t::result_type(R_Hz * 1e-9, t1, T_ns);
    }

getH1()

JPDF_t::result_type JFIT::JRegressor< JPoint4E, JGandalf >::getH1 ( const double  D, const double  ct, const double  t  ) const

inline

Get signal hypothesis value per 1 GeV for bright point emission PDF.

Parameters

D	hit distance from shower vertex [m]
ct	cosine of the HIT angle
t	arrival time of the light

Returns

hypothesis value / GeV

Definition at line 210 of file JShowerBrightPointRegressor.hh.

    {
      using namespace JPP;
 
      JPDF_t::result_type h1 = JMATH::zero;
        
      for (int i = 0; i != NUMBER_OF_PDFS; ++i) {
 
        if (!pdf[i].empty() && D <= pdf[i].getXmax()) {
 
          try {
 
            JPDF_t::result_type y1 = pdf[i](std::max(D, pdf[i].getXmin()), ct, t);
            
            // safety measures
            
            if (y1.f <= 0.0) {
              y1.f  = 0.0;
              y1.fp = 0.0;
            }
            
            if (y1.v <= 0.0) {
              y1.v  = 0.0;
            }
            
            h1 += y1;
          
          }
          catch(JLANG::JException& error) {
            ERROR(error << std::endl);
          }
        }
      }
 
      return h1;
    }

getRmax()

double JFIT::JRegressor< JPoint4E, JGandalf >::getRmax ( ) const

inline

Get maximal road width of PDF.

Returns

road width [m]

Definition at line 254 of file JShowerBrightPointRegressor.hh.

    {
      using namespace JPP;
 
      double xmax = 0.0;
 
      for (int i = 0; i != NUMBER_OF_PDFS; ++i) {
 
        if (!pdf[i].empty() && pdf[i].getXmax() > xmax) {
          xmax = pdf[i].getXmax();
        }
        
      }
 
      return xmax;
    }

operator()() [2/3]

result_type JFIT::JAbstractRegressor< JPoint4E , JGandalf >::operator() ( const JPoint4E &  value, T  __begin, T  __end  )

inlineinherited

Global fit.

Parameters

value	start value
__begin	begin of data set
__end	end of data set

Returns

chi2

Definition at line 94 of file JRegressor.hh.

    {
      static_cast<minimiser_type&>(*this).value = value;
 
      return static_cast<minimiser_type&>(*this)(static_cast<regressor_type&>(*this), __begin, __end);
    }

operator()() [3/3]

template<class JModel_t >

template<class JFunction_t , class T , class ... Args>

result_type JFIT::JGandalf< JModel_t >::operator() ( const JFunction_t &  fit, T  __begin, T  __end, Args ...  args  )

inlineinherited

Multi-dimensional fit of multiple data sets.

The fit function should return the chi2 as well as the partial derivatives for the current value of the model and a given data point.

Parameters

fit	fit function
__begin	begin of data
__end	end of data
args	optional data

Returns

chi2 and gradient

Definition at line 160 of file JGandalf.hh.

    {
      using namespace std;
      using namespace JPP;
 
      // note that all model values should be assigned to the start value of the model before use
      // because the actual list of model parameters can vary from fit to fit
      // (e.g. if model consists of a container).
 
      const size_t N = parameters.size();
 
      V.resize(N);
      h.resize(N);
      x.resize(N);
 
      previous.result.chi2 = numeric_limits<double>::max();
 
      current.result.chi2     = numeric_limits<double>::max();
      current.result.gradient = value;
      current.result.gradient = zero;
 
      error     = value;
      error     = zero;
 
      lambda    = LAMBDA_MIN;
 
      for (numberOfIterations = 0; numberOfIterations != MAXIMUM_ITERATIONS; ++numberOfIterations) {
 
        DEBUG("step:     " << numberOfIterations << endl);
 
        reset();
 
        update(fit, __begin, __end, args...);
 
        DEBUG("lambda:   " << FIXED(12,5) << lambda              << endl);
        DEBUG("chi2:     " << FIXED(12,5) << current.result.chi2 << endl);
 
        if (current.result.chi2 < previous.result.chi2) {
 
          if (numberOfIterations != 0) {
 
            const double tolerance = EPSILON * (EPSILON_ABSOLUTE ? 1.0 : fabs(previous.result.chi2)); 
 
            if (fabs(previous.result.chi2 - current.result.chi2) <= tolerance) {
 
              // normal end
 
              const result_type result = current.result;
 
              reset();
 
              update(fit, __begin, __end, args...);
 
              try {
                V.invert();
              }
              catch (const exception& error) {}
 
              for (size_t i = 0; i != N; ++i) {
                get(error, parameters[i]) = sqrt(V(i,i));
              }
 
              return result;
            }
 
            if (lambda > LAMBDA_MIN) {
              lambda /= LAMBDA_DOWN;
            }
          }
 
          // store current values
 
          previous.value  = value;
          previous.result = current.result;
 
        } else {
 
          value   = previous.value;                                 // restore value
 
          lambda *= LAMBDA_UP;
          
          if (lambda > LAMBDA_MAX) {
            break;
          }
 
          reset();
 
          update(fit, __begin, __end, args...);
        }
 
        DEBUG("Hesse matrix:" << endl << V << endl);
 
        // force definite positiveness
 
        for (size_t i = 0; i != N; ++i) {
 
          if (V(i,i) < PIVOT) {
            V(i,i) = PIVOT;
          }
 
          h[i] = 1.0 / sqrt(V(i,i));
        }
 
        // normalisation
        
        for (size_t row = 0; row != N; ++row) {
          for (size_t col = 0; col != row; ++col) {
            V(row,col) *= h[row] * h[col];
            V(col,row)  = V(row,col);
          }
        }
        
        for (size_t i = 0; i != N; ++i) {
          V(i,i) = 1.0 + lambda;
        }
 
        // solve A x = b
 
        for (size_t col = 0; col != N; ++col) {
          x[col] = h[col] * get(current.result.gradient, parameters[col]);
        }
 
        try {
          V.solve(x);
        }
        catch (const exception& error) {
 
          ERROR("JGandalf: " << error.what() << endl << V << endl);
 
          break;
        }
 
        // update value
 
        for (size_t row = 0; row != N; ++row) {
 
          DEBUG("u[" << noshowpos << setw(3) << row << "] = " << showpos << FIXED(15,5) << get(value, parameters[row]));
 
          get(value, parameters[row]) -= h[row] * x[row];
 
          DEBUG(" -> " << FIXED(15,5) << get(value, parameters[row]) << noshowpos << endl);
        }
 
        model(value);
      }
 
      // abnormal end
 
      const result_type result = previous.result;
 
      value   = previous.value;                                     // restore value
 
      reset();
 
      update(fit, __begin, __end, args...);
 
      try {
        V.invert();
      }
      catch (const exception& error) {}
 
      for (size_t i = 0; i != N; ++i) {
        get(error, parameters[i]) = sqrt(V(i,i));
      }
 
      return result;
    }

reset()

template<class JModel_t >

void JFIT::JGandalf< JModel_t >::reset ( )

inlineprivateinherited

Reset current parameters.

Definition at line 349 of file JGandalf.hh.

    {
      using namespace JPP;
 
      current.result.chi2     = 0.0;
      current.result.gradient = zero;
 
      V.reset();          
    }

update() [1/2]

template<class JModel_t >

template<class JFunction_t , class T , class ... Args>

void JFIT::JGandalf< JModel_t >::update ( const JFunction_t &  fit, T  __begin, T  __end, Args ...  args  )

inlineprivateinherited

Recursive method to update current parameters.

Parameters

fit	fit function
__begin	begin of data
__end	end of data
args	optional data

Definition at line 369 of file JGandalf.hh.

    {
      for (T i = __begin; i != __end; ++i) {
 
        const result_type& result = fit(value, *i);
 
        current.result.chi2     += result.chi2;
        current.result.gradient += result.gradient;
 
        for (size_t row = 0; row != parameters.size(); ++row) {
          for (size_t col = row; col != parameters.size(); ++col) {
            V(row,col) += get(result.gradient, parameters[row]) * get(result.gradient, parameters[col]);
          }
        }
      }
 
      update(fit, args...);
    }

update() [2/2]

template<class JModel_t >

template<class JFunction_t >

void JFIT::JGandalf< JModel_t >::update ( const JFunction_t &  fit )

inlineprivateinherited

Termination method to update current parameters.

Parameters

fit	fit function

Definition at line 395 of file JGandalf.hh.

    {
      for (size_t row = 0; row != parameters.size(); ++row) {
        for (size_t col = 0; col != row; ++col) {
          V(row,col) = V(col,row);
        }
      }
    }

get() [1/6]

template<class JModel_t >

static double JFIT::JGandalf< JModel_t >::get ( const JModel_t &  model, double JModel_t::*  parameter  )

inlinestaticprivateinherited

Read/write access to parameter value by data member.

Parameters

model	model
parameter	parameter

Returns

value

Definition at line 412 of file JGandalf.hh.

    {
      return model.*parameter;
    }

get() [2/6]

template<class JModel_t >

static double& JFIT::JGandalf< JModel_t >::get ( JModel_t &  model, double JModel_t::*  parameter  )

inlinestaticprivateinherited

Read/write access to parameter value by data member.

Parameters

model	model
parameter	parameter

Returns

value

Definition at line 425 of file JGandalf.hh.

    {
      return model.*parameter;
    }

get() [3/6]

template<class JModel_t >

static double JFIT::JGandalf< JModel_t >::get ( const JModel_t &  model, const size_t  index  )

inlinestaticprivateinherited

Read/write access to parameter value by index.

Parameters

model	model
index	index

Returns

value

Definition at line 438 of file JGandalf.hh.

    {
      return model[index];
    }

get() [4/6]

template<class JModel_t >

static double& JFIT::JGandalf< JModel_t >::get ( JModel_t &  model, const size_t  index  )

inlinestaticprivateinherited

Read/write access to parameter value by index.

Parameters

model	model
index	index

Returns

value

Definition at line 451 of file JGandalf.hh.

    {
      return model[index];
    }

get() [5/6]

template<class JModel_t >

static double JFIT::JGandalf< JModel_t >::get ( const JModel_t &  model, const int  index  )

inlinestaticprivateinherited

Read/write access to parameter value by index.

Parameters

model	model
index	index

Returns

value

Definition at line 464 of file JGandalf.hh.

    {
      return model[index];
    }

get() [6/6]

template<class JModel_t >

static double& JFIT::JGandalf< JModel_t >::get ( JModel_t &  model, const int  index  )

inlinestaticprivateinherited

Read/write access to parameter value by index.

Parameters

model	model
index	index

Returns

value

Definition at line 477 of file JGandalf.hh.

    {
      return model[index];
    }

Member Data Documentation

T_ns

JTimeRange JFIT::JRegressor< JPoint4E, JGandalf >::T_ns

static

Time window with respect to Cherenkov hypothesis [ns].

Default values.

Definition at line 271 of file JShowerBrightPointRegressor.hh.

Vmax_npe

double JFIT::JRegressor< JPoint4E, JGandalf >::Vmax_npe = std::numeric_limits<double>::max()

static

Maximal integral of PDF [npe].

Definition at line 272 of file JShowerBrightPointRegressor.hh.

NUMBER_OF_PDFS

const int JFIT::JRegressor< JPoint4E, JGandalf >::NUMBER_OF_PDFS = 2

static

Definition at line 274 of file JShowerBrightPointRegressor.hh.

pdf_t

const JPDFType_t JFIT::JRegressor< JPoint4E, JGandalf >::pdf_t

static

Initial value:

= { DIRECT_LIGHT_FROM_BRIGHT_POINT, 
                                                               SCATTERED_LIGHT_FROM_BRIGHT_POINT }

PDF types.

Definition at line 276 of file JShowerBrightPointRegressor.hh.

pdf

JPDF_t JFIT::JRegressor< JPoint4E, JGandalf >::pdf[NUMBER_OF_PDFS]

PDF.

Definition at line 278 of file JShowerBrightPointRegressor.hh.

MAXIMUM_ITERATIONS

template<class JModel_t >

int JFIT::JGandalf< JModel_t >::MAXIMUM_ITERATIONS = 1000

staticinherited

maximal number of iterations

maximal number of iterations.

Definition at line 329 of file JGandalf.hh.

EPSILON

template<class JModel_t >

double JFIT::JGandalf< JModel_t >::EPSILON = 1.0e-3

staticinherited

maximal distance to minimum

maximal distance to minimum.

Definition at line 330 of file JGandalf.hh.

EPSILON_ABSOLUTE

template<class JModel_t >

bool JFIT::JGandalf< JModel_t >::EPSILON_ABSOLUTE = false

staticinherited

set epsilon to absolute difference instead of relative

set epsilon to absolute difference instead of relative.

Definition at line 331 of file JGandalf.hh.

LAMBDA_MIN

template<class JModel_t >

double JFIT::JGandalf< JModel_t >::LAMBDA_MIN = 0.01

staticinherited

minimal value control parameter

Definition at line 332 of file JGandalf.hh.

LAMBDA_MAX

template<class JModel_t >

double JFIT::JGandalf< JModel_t >::LAMBDA_MAX = 100.0

staticinherited

maximal value control parameter

Definition at line 333 of file JGandalf.hh.

LAMBDA_UP

template<class JModel_t >

double JFIT::JGandalf< JModel_t >::LAMBDA_UP = 10.0

staticinherited

multiplication factor control parameter

Definition at line 334 of file JGandalf.hh.

LAMBDA_DOWN

template<class JModel_t >

double JFIT::JGandalf< JModel_t >::LAMBDA_DOWN = 10.0

staticinherited

multiplication factor control parameter

Definition at line 335 of file JGandalf.hh.

PIVOT

template<class JModel_t >

double JFIT::JGandalf< JModel_t >::PIVOT = std::numeric_limits<double>::epsilon()

staticinherited

minimal value diagonal element of Hesse matrix

minimal value diagonal element of matrix

Definition at line 336 of file JGandalf.hh.

parameters

template<class JModel_t >

std::vector<parameter_type> JFIT::JGandalf< JModel_t >::parameters

inherited

fit parameters

Definition at line 338 of file JGandalf.hh.

numberOfIterations

template<class JModel_t >

int JFIT::JGandalf< JModel_t >::numberOfIterations

inherited

number of iterations

Definition at line 339 of file JGandalf.hh.

lambda

template<class JModel_t >

double JFIT::JGandalf< JModel_t >::lambda

inherited

control parameter

Definition at line 340 of file JGandalf.hh.

value

template<class JModel_t >

JModel_t JFIT::JGandalf< JModel_t >::value

inherited

value

Definition at line 341 of file JGandalf.hh.

error

template<class JModel_t >

JModel_t JFIT::JGandalf< JModel_t >::error

inherited

error

Definition at line 342 of file JGandalf.hh.

V

template<class JModel_t >

JMATH::JMatrixNS JFIT::JGandalf< JModel_t >::V

inherited

Hesse matrix.

Definition at line 343 of file JGandalf.hh.

h

template<class JModel_t >

std::vector<double> JFIT::JGandalf< JModel_t >::h

privateinherited

Definition at line 482 of file JGandalf.hh.

x

template<class JModel_t >

JMATH::JVectorND JFIT::JGandalf< JModel_t >::x

privateinherited

Definition at line 483 of file JGandalf.hh.

result

template<class JModel_t >

result_type JFIT::JGandalf< JModel_t >::result

inherited

Definition at line 486 of file JGandalf.hh.

struct { ... } JFIT::JGandalf< JModel_t >::current

struct { ... } JFIT::JGandalf< JModel_t >::previous

debug

template<class T >

int JEEP::JMessage< T >::debug = 0

staticinherited

debug level (default is off).

Definition at line 45 of file JMessage.hh.

---

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

• JShowerBrightPointRegressor.hh