Base class for result of ROOT Fit. More...

#include <JRootfit.hh>

Inheritance diagram for JROOT::JRootfit_t< JFs_t >:

Public Types
typedef JFIT_LOCAL::JTypedef_t< JModel_t >::parameter_type	parameter_type
	Data type of fit parameter.

Public Member Functions
const JFs_t &	getFunction () const
	Get function.

size_t	getNumberOfParameters () const
	Get number of parameters.

size_t	getNumberOfFreeParameters () const
	Get number of free parameters.

size_t	getN () const
	Get number of data points.

double	getChi2 () const
	Get chi2.

int	getNDF () const
	Get number of degrees of freedom.

double	getValue (size_t i) const
	Get value of parameter at given index.

double	getError (size_t i) const
	Get error of parameter at given index.

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.

Public Attributes
std::vector< parameter_type >	parameters
	fit parameters

int	numberOfIterations
	number of iterations

double	lambda
	control parameter

JModel_t	value
	value

JModel_t	error
	error

JMATH::JMatrixNS	V
	Hesse matrix.

Static Public Attributes
static int	MAXIMUM_ITERATIONS = 1000
	maximal number of iterations

static double	EPSILON = 1.0e-3
	maximal distance to minimum

static bool	EPSILON_ABSOLUTE = false
	set epsilon to absolute difference instead of relative

static double	LAMBDA_MIN = 0.01
	minimal value control parameter

static double	LAMBDA_MAX = 100.0
	maximal value control parameter

static double	LAMBDA_UP = 10.0
	multiplication factor control parameter

static double	LAMBDA_DOWN = 10.0
	multiplication factor control parameter

static double	PIVOT = std::numeric_limits<double>::epsilon()
	minimal value diagonal element of Hesse matrix

static int	debug = 0
	debug level (default is off).

Protected Member Functions
	JRootfit_t ()
	Default constructor.

Protected Attributes
size_t	npx
	number of data points

double	chi2
	chi2

Private Member Functions
void	reset ()
	Reset current parameters.

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.

template<class JFunction_t >
void	update (const JFunction_t &fit)
	Termination method to update current parameters.

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

static double &	get (JModel_t &model, double JModel_t::*parameter)
	Read/write access to parameter value by data member.

static double	get (const JModel_t &model, const size_t index)
	Read/write access to parameter value by index.

static double &	get (JModel_t &model, const size_t index)
	Read/write access to parameter value by index.

static double	get (const JModel_t &model, const int index)
	Read/write access to parameter value by index.

static double &	get (JModel_t &model, const int index)
	Read/write access to parameter value by index.

Private Attributes
std::vector< double >	h

JMATH::JVectorND	x

struct {

result_type result

}	current

struct {

JModel_t value

result_type result

}	previous

Detailed Description

template<class JFs_t>
class JROOT::JRootfit_t< JFs_t >

Base class for result of ROOT Fit.

Definition at line 828 of file JRootfit.hh.

Member Typedef Documentation

◆ parameter_type

template<class JModel_t >

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

inherited

Data type of fit parameter.

Definition at line 96 of file JGandalf.hh.

Constructor & Destructor Documentation

◆ JRootfit_t()

template<class JFs_t >

JROOT::JRootfit_t< JFs_t >::JRootfit_t ( )

inlineprotected

Default constructor.

Definition at line 835 of file JRootfit.hh.

                 :
      npx (0),
      chi2(0.0)
    {}

Member Function Documentation

◆ getFunction()

template<class JFs_t >

const JFs_t & JROOT::JRootfit_t< JFs_t >::getFunction ( ) const

inline

Get function.

Returns: function.

Definition at line 850 of file JRootfit.hh.

    {
      return this->value;
    }

◆ getNumberOfParameters()

template<class JFs_t >

size_t JROOT::JRootfit_t< JFs_t >::getNumberOfParameters ( ) const

inline

Get number of parameters.

Returns: number of parameters

Definition at line 861 of file JRootfit.hh.

    {
      return JFs_t::parameters.size();
    }

◆ getNumberOfFreeParameters()

template<class JFs_t >

size_t JROOT::JRootfit_t< JFs_t >::getNumberOfFreeParameters ( ) const

inline

Get number of free parameters.

Returns: number of free parameters

Definition at line 872 of file JRootfit.hh.

    {
      return this->parameters.size();
    }

◆ getN()

template<class JFs_t >

size_t JROOT::JRootfit_t< JFs_t >::getN ( ) const

inline

Get number of data points.

Returns: number of data points

Definition at line 883 of file JRootfit.hh.

    {
      return npx;
    }

◆ getChi2()

template<class JFs_t >

double JROOT::JRootfit_t< JFs_t >::getChi2 ( ) const

inline

Get chi2.

Returns: chi2

Definition at line 894 of file JRootfit.hh.

    {
      return chi2;
    }

◆ getNDF()

template<class JFs_t >

int JROOT::JRootfit_t< JFs_t >::getNDF ( ) const

inline

Get number of degrees of freedom.

Returns: number of degrees of freedom

Definition at line 905 of file JRootfit.hh.

    {
      return (int) getN() - (int) getNumberOfFreeParameters();
    }

◆ getValue()

template<class JFs_t >

double JROOT::JRootfit_t< JFs_t >::getValue ( size_t i ) const

inline

Get value of parameter at given index.

Parameters

i index

Definition at line 916 of file JRootfit.hh.

    {
      return getParameter(this->value, i);
    }

◆ getError()

template<class JFs_t >

double JROOT::JRootfit_t< JFs_t >::getError ( size_t i ) const

inline

Get error of parameter at given index.

Parameters

i index

Definition at line 927 of file JRootfit.hh.

    {
      return getParameter(this->error, i);
    }

◆ operator()()

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 161 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 350 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 370 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 396 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 413 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 426 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 439 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 452 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 465 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 478 of file JGandalf.hh.

    {
      return model[index];
    }

Member Data Documentation

◆ npx

template<class JFs_t >

size_t JROOT::JRootfit_t< JFs_t >::npx

protected

number of data points

Definition at line 841 of file JRootfit.hh.

◆ chi2

template<class JFs_t >

double JROOT::JRootfit_t< JFs_t >::chi2

protected

chi2

Definition at line 842 of file JRootfit.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 330 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 331 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 332 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 333 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 334 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 335 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 336 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 337 of file JGandalf.hh.

◆ parameters

template<class JModel_t >

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

inherited

fit parameters

Definition at line 339 of file JGandalf.hh.

◆ numberOfIterations

template<class JModel_t >

int JFIT::JGandalf< JModel_t >::numberOfIterations

inherited

number of iterations

Definition at line 340 of file JGandalf.hh.

◆ lambda

template<class JModel_t >

double JFIT::JGandalf< JModel_t >::lambda

inherited

control parameter

Definition at line 341 of file JGandalf.hh.

◆ value

template<class JModel_t >

JModel_t JFIT::JGandalf< JModel_t >::value

inherited

value

Definition at line 342 of file JGandalf.hh.

◆ error

template<class JModel_t >

JModel_t JFIT::JGandalf< JModel_t >::error

inherited

error

Definition at line 343 of file JGandalf.hh.

◆ V

template<class JModel_t >

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

inherited

Hesse matrix.

Definition at line 344 of file JGandalf.hh.

◆ h

template<class JModel_t >

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

privateinherited

Definition at line 483 of file JGandalf.hh.

◆ x

template<class JModel_t >

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

privateinherited

Definition at line 484 of file JGandalf.hh.

◆ result

template<class JModel_t >

result_type JFIT::JGandalf< JModel_t >::result

inherited

Definition at line 487 of file JGandalf.hh.

◆ [struct]

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

◆ [struct]

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 class was generated from the following file:

JRootfit.hh

Public Types

Public Member Functions

Public Attributes

Static Public Attributes

Protected Member Functions

Protected Attributes

Private Member Functions

Static Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ parameter_type

Constructor & Destructor Documentation

◆ JRootfit_t()

Member Function Documentation

◆ getFunction()

◆ getNumberOfParameters()

◆ getNumberOfFreeParameters()

◆ getN()

◆ getChi2()

◆ getNDF()

◆ getValue()

◆ getError()

◆ operator()()

◆ reset()

◆ update() [1/2]

◆ update() [2/2]

◆ get() [1/6]

◆ get() [2/6]

◆ get() [3/6]

◆ get() [4/6]

◆ get() [5/6]

◆ get() [6/6]

Member Data Documentation

◆ npx

◆ chi2

◆ MAXIMUM_ITERATIONS

◆ EPSILON

◆ EPSILON_ABSOLUTE

◆ LAMBDA_MIN

◆ LAMBDA_MAX

◆ LAMBDA_UP

◆ LAMBDA_DOWN

◆ PIVOT

◆ parameters

◆ numberOfIterations

◆ lambda

◆ value

◆ error

◆ V

◆ h

◆ x

◆ result

◆ [struct]

◆ [struct]

◆ debug