LArSoft  v09_90_00
Liquid Argon Software toolkit - https://larsoft.org/
lar::util::QuadraticFit< T > Class Template Referenceabstract

Performs a second-degree fit of data. More...

#include "SimpleFits.h"

Inheritance diagram for lar::util::QuadraticFit< T >:
lar::util::details::SimplePolyFitterBase< T, 2U > lar::util::details::SimpleFitterInterface< T, D+1 > lar::util::details::SimplePolyFitterDataBase< T, D >

Public Types

using FitParameters_t = std::array< Data_t, NParams >
 type of set of fit parameters More...
 
using FitMatrix_t = std::array< Data_t, sqr(NParams)>
 
using Data_t = typename Base_t::Data_t
 type of the data More...
 
using MatrixOps = FastMatrixOperations< Data_t, NParams >
 
using Measurement_t = typename Collector_t::Measurement_t
 type of measurement without uncertainty More...
 
using MeasurementAndUncertainty_t = typename Collector_t::MeasurementAndUncertainty_t
 type of measurement with uncertainty More...
 

Public Member Functions

virtual Data_t ChiSquare () const override
 Returns the $ \chi^{2} $ of the fit. More...
 
virtual bool isValid () const override
 Returns if the fit has valid results. More...
 
bool FillResults (FitParameters_t &params, FitMatrix_t &Xmat, Data_t &det, FitMatrix_t &Smat) const override
 Fills the specified parameters. More...
 
bool FillResults (FitParameters_t &params, FitParameters_t &paramerrors, FitMatrix_t &Xmat, Data_t &det, FitMatrix_t &Smat) const override
 Fills the specified parameters. More...
 
bool FillResults (FitParameters_t &params, FitParameters_t &paramerrors) const override
 Fills the specified parameters. More...
 
virtual bool FillResults (FitParameters_t &params, FitMatrix_t &Xmat, Data_t &det, FitMatrix_t &Smat) const =0
 Fills the specified parameters. More...
 
virtual bool FillResults (FitParameters_t &params, FitParameters_t &paramerrors, FitMatrix_t &Xmat, Data_t &det, FitMatrix_t &Smat) const =0
 Fills the specified parameters. More...
 
virtual Data_t Evaluate (Data_t x) const override
 Evaluates the fitted function at the specified point. More...
 
virtual Data_t Evaluate (Data_t x) const =0
 Evaluates the fitted function at the specified point. More...
 
Data_t operator() (Data_t x) const
 Evaluates the fitted function; alias of Evaluate() More...
 
void clear ()
 Clears all the statistics. More...
 
template<typename Stream >
void PrintStats (Stream &out) const
 Prints the collected statistics into a stream. More...
 
Fitting
virtual FitParameters_t FitParameters () const override
 Computes and returns all the parameters of the fit result. More...
 
virtual FitParameters_t FitParameterErrors () const override
 Computes and returns all the parameter errors of the fit result. More...
 
virtual FitMatrix_t FitParameterCovariance () const override
 Computes and returns all the covariance matrix of the fit result. More...
 
virtual Data_t FitParameter (unsigned int n) const override
 Returns the parameter n of the fit result. More...
 
virtual Data_t FitParameterError (unsigned int n) const override
 Returns the error on parameter n of the fit result. More...
 
virtual int NDF () const override
 Returns the degrees of freedom in the determination of the fit. More...
 
Add elements
See also
FitDataCollector
bool add (Data_t x, Data_t y, Data_t sy=Data_t(1.0))
 
bool add (Measurement_t value, Data_t sy=Data_t(1.0))
 
bool add (MeasurementAndUncertainty_t value)
 
template<typename Iter >
void add_without_uncertainty (Iter begin, Iter end)
 
template<typename Iter , typename Pred >
void add_without_uncertainty (Iter begin, Iter end, Pred extractor)
 
template<typename Cont , typename Pred >
void add_without_uncertainty (Cont cont, Pred extractor)
 
template<typename Cont >
void add_without_uncertainty (Cont cont)
 
template<typename VIter , typename UIter , typename VPred , typename UPred = identity>
unsigned int add_with_uncertainty (VIter begin_value, VIter end_value, UIter begin_uncertainty, VPred value_extractor, UPred uncertainty_extractor=UPred())
 
template<typename Iter >
unsigned int add_with_uncertainty (Iter begin, Iter end)
 
template<typename Cont >
unsigned int add_with_uncertainty (Cont cont)
 
Statistic retrieval
See also
FitDataCollector
int N () const
 
Data_t AverageUncertainty () const
 

Static Public Member Functions

static FitParameters_t ExtractParameterErrors (FitMatrix_t const &Smat)
 Extracts parameter errors from diagonal of the covarriance matrix. More...
 
static constexpr Data_t sqr (Data_t v)
 Returns the square of the specified data value. More...
 
template<typename V >
static constexpr V sqr (V const &v)
 Returns the square of the specified value. More...
 
static constexpr Data_t cube (Data_t v)
 Returns the cube of the specified data value. More...
 

Static Public Attributes

static constexpr unsigned int Degree
 Degree of the fit. More...
 
static constexpr unsigned int NParams
 Number of parameters in the fit. More...
 

Protected Member Functions

virtual FitParameters_t FitParameterErrors (FitMatrix_t const &Smat) const
 Computes and returns all the parameter errors of the fit result. More...
 
virtual FitMatrix_t MakeMatrixX () const
 Fills and returns the matrix of x^n sum coefficients ( { x^(i+j) } ) More...
 
virtual FitParameters_t MakeMatrixY () const
 Fills and returns the matrix (vector) of x^n y sum coefficients. More...
 
virtual Data_t Determinant (FitMatrix_t const &mat) const
 Computes the determinant of a matrix. More...
 
virtual FitMatrix_t InvertMatrix (FitMatrix_t const &mat, Data_t det) const
 Computes the inverse of a matrix (using provided determinant) More...
 
virtual FitMatrix_t InvertMatrix (FitMatrix_t const &mat) const
 Computes the inverse of a matrix. More...
 
virtual FitParameters_t MatrixProduct (FitMatrix_t const &mat, FitParameters_t const &vec) const
 Computes the product of a FitMatrix_t and a FitParameters_t. More...
 
Data_t XN (unsigned int n) const
 Returns the weighted sum of x^n. More...
 
Data_t XNY (unsigned int n) const
 Returns the weighted sum of x^n y. More...
 
Data_t I () const
 Aliases. More...
 
Data_t X () const
 Aliases. More...
 
Data_t X2 () const
 Aliases. More...
 
Data_t X3 () const
 Aliases. More...
 
Data_t X4 () const
 Aliases. More...
 
Data_t Y () const
 Aliases. More...
 
Data_t XY () const
 Aliases. More...
 
Data_t X2Y () const
 Aliases. More...
 
Data_t Y2 () const
 Aliases. More...
 
virtual FitParameters_t FitParameters (FitMatrix_t const &Xmat) const
 Returns the fitted parameters using the provided information. More...
 
virtual FitParameters_t FitParameters (FitMatrix_t const &Smat, Data_t) const
 Returns the fitted parameters using the provided information. More...
 
virtual Data_t Param (unsigned int n, FitMatrix_t const &Xmat) const
 Computes a single fit parameter using the given information. More...
 
virtual Data_t Param (unsigned int n, FitMatrix_t const &Xmat, Data_t detXmat) const
 Computes a single fit parameter using the given information. More...
 

Protected Attributes

Collector_t stats
 statistics collected from fit data input More...
 

Private Types

using Base_t = details::SimplePolyFitterBase< T, 2U >
 

Detailed Description

template<typename T>
class lar::util::QuadraticFit< T >

Performs a second-degree fit of data.


Template Parameters
Ttype of the quantities
Wtype of the weight (as T by default)

The quadratic fit connects measurements $ ( y_{i} \pm \sigma_{y,i} ) $ with a parameter $ ( x_{i} ) $ not affected by uncertainty. The returned parameters describe a quadratic curve $ f(x) = a_{0} + a_{1} x + a_{2} x^{2} $ obtained by minimization of $ \chi^{2} = \sum_{i} \frac{ \left(y_{i} - f(x_{i}) \right)^{2} }{ \sigma^{2}_{y,i} }$

This saves having to link to ROOT for the simplest cases.

This simple quadratic fitter does not store any result: each time a result is requested, it is computed anew. In particular that is true also for ChiSquare(), that requires the full parameters set and therefore reruns the full fit (FitParameters()) and for the covariance matrix of the parameters.

Definition at line 933 of file SimpleFits.h.

Member Typedef Documentation

template<typename T >
using lar::util::QuadraticFit< T >::Base_t = details::SimplePolyFitterBase<T, 2U>
private

Definition at line 934 of file SimpleFits.h.

template<typename T >
using lar::util::details::SimplePolyFitterBase< T, D >::Data_t = typename Base_t::Data_t

type of the data

Definition at line 649 of file SimpleFits.h.

template<typename T >
using lar::util::QuadraticFit< T >::FitMatrix_t = std::array<Data_t, sqr(NParams)>

Definition at line 951 of file SimpleFits.h.

template<typename T >
using lar::util::QuadraticFit< T >::FitParameters_t = std::array<Data_t, NParams>

type of set of fit parameters

Definition at line 950 of file SimpleFits.h.

using lar::util::details::SimpleFitterInterface< T, N >::MatrixOps = FastMatrixOperations<Data_t, NParams>
inherited

Definition at line 465 of file SimpleFits.h.

template<typename T , unsigned int D>
using lar::util::details::SimplePolyFitterDataBase< T, D >::Measurement_t = typename Collector_t::Measurement_t
inherited

type of measurement without uncertainty

Definition at line 352 of file SimpleFits.h.

template<typename T , unsigned int D>
using lar::util::details::SimplePolyFitterDataBase< T, D >::MeasurementAndUncertainty_t = typename Collector_t::MeasurementAndUncertainty_t
inherited

type of measurement with uncertainty

Definition at line 355 of file SimpleFits.h.

Member Function Documentation

template<typename T , unsigned int D>
bool lar::util::details::SimplePolyFitterDataBase< T, D >::add ( Data_t  x,
Data_t  y,
Data_t  sy = Data_t(1.0) 
)
inlineinherited

Definition at line 363 of file SimpleFits.h.

Referenced by cluster::ClusterParamsAlg::IntegrateFitCharge().

363 { return stats.add(x, y, sy); }
Float_t x
Definition: compare.C:6
Float_t y
Definition: compare.C:6
Collector_t stats
statistics collected from fit data input
Definition: SimpleFits.h:446
bool add(Data_t x, Data_t y, Data_t sy=Data_t(1.0))
Adds one entry with specified x, y and uncertainty.
Definition: SimpleFits.h:1412
template<typename T , unsigned int D>
bool lar::util::details::SimplePolyFitterDataBase< T, D >::add ( Measurement_t  value,
Data_t  sy = Data_t(1.0) 
)
inlineinherited

Definition at line 365 of file SimpleFits.h.

365 { return stats.add(value, sy); }
double value
Definition: spectrum.C:18
Collector_t stats
statistics collected from fit data input
Definition: SimpleFits.h:446
bool add(Data_t x, Data_t y, Data_t sy=Data_t(1.0))
Adds one entry with specified x, y and uncertainty.
Definition: SimpleFits.h:1412
template<typename T , unsigned int D>
bool lar::util::details::SimplePolyFitterDataBase< T, D >::add ( MeasurementAndUncertainty_t  value)
inlineinherited

Definition at line 367 of file SimpleFits.h.

367 { return stats.add(value); }
double value
Definition: spectrum.C:18
Collector_t stats
statistics collected from fit data input
Definition: SimpleFits.h:446
bool add(Data_t x, Data_t y, Data_t sy=Data_t(1.0))
Adds one entry with specified x, y and uncertainty.
Definition: SimpleFits.h:1412
template<typename T , unsigned int D>
template<typename VIter , typename UIter , typename VPred , typename UPred = identity>
unsigned int lar::util::details::SimplePolyFitterDataBase< T, D >::add_with_uncertainty ( VIter  begin_value,
VIter  end_value,
UIter  begin_uncertainty,
VPred  value_extractor,
UPred  uncertainty_extractor = UPred() 
)
inlineinherited

Definition at line 394 of file SimpleFits.h.

399  {
401  begin_value, end_value, begin_uncertainty, value_extractor, uncertainty_extractor);
402  }
unsigned int add_with_uncertainty(VIter begin_value, VIter end_value, UIter begin_uncertainty, VPred value_extractor, UPred uncertainty_extractor=UPred())
Adds measurements with uncertainties from a sequence.
Definition: SimpleFits.h:1442
Collector_t stats
statistics collected from fit data input
Definition: SimpleFits.h:446
template<typename T , unsigned int D>
template<typename Iter >
unsigned int lar::util::details::SimplePolyFitterDataBase< T, D >::add_with_uncertainty ( Iter  begin,
Iter  end 
)
inlineinherited

Definition at line 405 of file SimpleFits.h.

406  {
408  }
unsigned int add_with_uncertainty(VIter begin_value, VIter end_value, UIter begin_uncertainty, VPred value_extractor, UPred uncertainty_extractor=UPred())
Adds measurements with uncertainties from a sequence.
Definition: SimpleFits.h:1442
decltype(auto) constexpr end(T &&obj)
ADL-aware version of std::end.
Definition: StdUtils.h:77
decltype(auto) constexpr begin(T &&obj)
ADL-aware version of std::begin.
Definition: StdUtils.h:69
Collector_t stats
statistics collected from fit data input
Definition: SimpleFits.h:446
template<typename T , unsigned int D>
template<typename Cont >
unsigned int lar::util::details::SimplePolyFitterDataBase< T, D >::add_with_uncertainty ( Cont  cont)
inlineinherited

Definition at line 411 of file SimpleFits.h.

412  {
413  return stats.add_with_uncertainty(cont);
414  }
unsigned int add_with_uncertainty(VIter begin_value, VIter end_value, UIter begin_uncertainty, VPred value_extractor, UPred uncertainty_extractor=UPred())
Adds measurements with uncertainties from a sequence.
Definition: SimpleFits.h:1442
Collector_t stats
statistics collected from fit data input
Definition: SimpleFits.h:446
template<typename T , unsigned int D>
template<typename Iter >
void lar::util::details::SimplePolyFitterDataBase< T, D >::add_without_uncertainty ( Iter  begin,
Iter  end 
)
inlineinherited

Definition at line 370 of file SimpleFits.h.

371  {
373  }
void add_without_uncertainty(Iter begin, Iter end)
Adds measurements from a sequence, with no uncertainty.
Definition: SimpleFits.h:136
decltype(auto) constexpr end(T &&obj)
ADL-aware version of std::end.
Definition: StdUtils.h:77
decltype(auto) constexpr begin(T &&obj)
ADL-aware version of std::begin.
Definition: StdUtils.h:69
Collector_t stats
statistics collected from fit data input
Definition: SimpleFits.h:446
template<typename T , unsigned int D>
template<typename Iter , typename Pred >
void lar::util::details::SimplePolyFitterDataBase< T, D >::add_without_uncertainty ( Iter  begin,
Iter  end,
Pred  extractor 
)
inlineinherited

Definition at line 376 of file SimpleFits.h.

377  {
378  stats.add_without_uncertainty(begin, end, extractor);
379  }
void add_without_uncertainty(Iter begin, Iter end)
Adds measurements from a sequence, with no uncertainty.
Definition: SimpleFits.h:136
decltype(auto) constexpr end(T &&obj)
ADL-aware version of std::end.
Definition: StdUtils.h:77
decltype(auto) constexpr begin(T &&obj)
ADL-aware version of std::begin.
Definition: StdUtils.h:69
Collector_t stats
statistics collected from fit data input
Definition: SimpleFits.h:446
template<typename T , unsigned int D>
template<typename Cont , typename Pred >
void lar::util::details::SimplePolyFitterDataBase< T, D >::add_without_uncertainty ( Cont  cont,
Pred  extractor 
)
inlineinherited

Definition at line 382 of file SimpleFits.h.

383  {
384  stats.add_without_uncertainty(cont, extractor);
385  }
void add_without_uncertainty(Iter begin, Iter end)
Adds measurements from a sequence, with no uncertainty.
Definition: SimpleFits.h:136
Collector_t stats
statistics collected from fit data input
Definition: SimpleFits.h:446
template<typename T , unsigned int D>
template<typename Cont >
void lar::util::details::SimplePolyFitterDataBase< T, D >::add_without_uncertainty ( Cont  cont)
inlineinherited

Definition at line 388 of file SimpleFits.h.

389  {
391  }
void add_without_uncertainty(Iter begin, Iter end)
Adds measurements from a sequence, with no uncertainty.
Definition: SimpleFits.h:136
Collector_t stats
statistics collected from fit data input
Definition: SimpleFits.h:446
template<typename T , unsigned int D>
Data_t lar::util::details::SimplePolyFitterDataBase< T, D >::AverageUncertainty ( ) const
inlineinherited

Definition at line 427 of file SimpleFits.h.

427 { return stats.AverageUncertainty(); }
Data_t AverageUncertainty() const
Returns an average of the uncertainties.
Definition: SimpleFits.h:280
Collector_t stats
statistics collected from fit data input
Definition: SimpleFits.h:446
template<typename T >
auto lar::util::QuadraticFit< T >::ChiSquare ( ) const
overridevirtual

Returns the $ \chi^{2} $ of the fit.

Returns
the $ \chi^{2} $ of the fit (not divided by NDF())

Reimplemented from lar::util::details::SimplePolyFitterBase< T, 2U >.

Definition at line 1703 of file SimpleFits.h.

References lar::util::details::FitDataCollector< T, D >::sqr(), X, X2, Y, and lar::util::details::FitDataCollector< T, D >::Y2().

1704 {
1705  FitParameters_t a = this->FitParameters();
1706  return Y2() - Data_t(2) * (a[0] * Y() + a[1] * XY() + a[2] * X2Y()) + sqr(a[0]) * I() +
1707  Data_t(2) * a[0] * (a[1] * X() + a[2] * X2()) + sqr(a[1]) * X2() +
1708  Data_t(2) * a[1] * (a[2] * X3()) + sqr(a[2]) * X4();
1709 } // QuadraticFit<T>::ChiSquare()
Data_t X4() const
Aliases.
Definition: SimpleFits.h:968
Data_t Y2() const
Aliases.
Definition: SimpleFits.h:972
typename Base_t::Data_t Data_t
type of the data
Definition: SimpleFits.h:649
std::array< Data_t, NParams > FitParameters_t
type of set of fit parameters
Definition: SimpleFits.h:950
Data_t Y() const
Aliases.
Definition: SimpleFits.h:969
virtual FitParameters_t FitParameters() const override
Computes and returns all the parameters of the fit result.
static constexpr V sqr(V const &v)
Returns the square of the specified value.
Definition: SimpleFits.h:440
Data_t X2Y() const
Aliases.
Definition: SimpleFits.h:971
Data_t I() const
Aliases.
Definition: SimpleFits.h:964
Data_t XY() const
Aliases.
Definition: SimpleFits.h:970
Data_t X() const
Aliases.
Definition: SimpleFits.h:965
Data_t X2() const
Aliases.
Definition: SimpleFits.h:966
Data_t X3() const
Aliases.
Definition: SimpleFits.h:967
template<typename T , unsigned int D>
void lar::util::details::SimplePolyFitterDataBase< T, D >::clear ( )
inlineinherited

Clears all the statistics.

Definition at line 419 of file SimpleFits.h.

419 { stats.clear(); }
void clear()
Clears all the statistics.
Definition: SimpleFits.h:1469
Collector_t stats
statistics collected from fit data input
Definition: SimpleFits.h:446
static constexpr Data_t lar::util::details::SimpleFitterInterface< T, N >::cube ( Data_t  v)
inlinestaticinherited

Returns the cube of the specified data value.

Definition at line 603 of file SimpleFits.h.

603 { return v * v * v; }
virtual Data_t lar::util::details::SimpleFitterInterface< T, N >::Determinant ( FitMatrix_t const &  mat) const
inlineprotectedvirtualinherited

Computes the determinant of a matrix.

Definition at line 607 of file SimpleFits.h.

608  {
609  return MatrixOps::Determinant(mat);
610  }
Float_t mat
Definition: plot.C:38
virtual Data_t lar::util::details::SimpleFitterInterface< T, N >::Evaluate ( Data_t  x) const
pure virtualinherited

Evaluates the fitted function at the specified point.

Parameters
xthe point where to evaluate the fit function
Returns
the value of the fit function

No check is performed whether the fit is valid.

virtual Data_t lar::util::details::SimplePolyFitterBase< T, D >::Evaluate ( Data_t  x) const
overridevirtualinherited

Evaluates the fitted function at the specified point.

Parameters
xthe point where to evaluate the fit function
Returns
the value of the fit function

No check is performed whether the fit is valid.

static FitParameters_t lar::util::details::SimplePolyFitterBase< T, D >::ExtractParameterErrors ( FitMatrix_t const &  Smat)
staticinherited

Extracts parameter errors from diagonal of the covarriance matrix.

virtual bool lar::util::details::SimpleFitterInterface< T, N >::FillResults ( FitParameters_t params,
FitMatrix_t Xmat,
Data_t det,
FitMatrix_t Smat 
) const
pure virtualinherited

Fills the specified parameters.

Parameters
paramsthe fitted values of the parameters
Xmatthe matrix of the x^n/s^2 sums
Smatthe covariance matrix
detthe determinant of Xmat
Returns
true if the fit is valid (i.e. if a unique solution exists)
virtual bool lar::util::details::SimpleFitterInterface< T, N >::FillResults ( FitParameters_t params,
FitParameters_t paramerrors,
FitMatrix_t Xmat,
Data_t det,
FitMatrix_t Smat 
) const
pure virtualinherited

Fills the specified parameters.

Parameters
paramsthe fitted values of the parameters
paramerrorsthe uncertainty on the fitted parameters
Xmatthe matrix of the x^n/s^2 sums
Smatthe covariance matrix
detthe determinant of Xmat
Returns
true if the fit is valid (i.e. if a unique solution exists)
bool lar::util::details::SimplePolyFitterBase< T, D >::FillResults ( FitParameters_t params,
FitMatrix_t Xmat,
Data_t det,
FitMatrix_t Smat 
) const
overrideinherited

Fills the specified parameters.

Parameters
paramsthe fitted values of the parameters
Xmatthe matrix of the x^n/s^2 sums
Smatthe covariance matrix
detthe determinant of Xmat
Returns
true if the fit is valid (i.e. if a unique solution exists)
bool lar::util::details::SimplePolyFitterBase< T, D >::FillResults ( FitParameters_t params,
FitParameters_t paramerrors,
FitMatrix_t Xmat,
Data_t det,
FitMatrix_t Smat 
) const
overrideinherited

Fills the specified parameters.

Parameters
paramsthe fitted values of the parameters
paramerrorsthe uncertainty on the fitted parameters
Xmatthe matrix of the x^n/s^2 sums
Smatthe covariance matrix
detthe determinant of Xmat
Returns
true if the fit is valid (i.e. if a unique solution exists)
bool lar::util::details::SimplePolyFitterBase< T, D >::FillResults ( FitParameters_t params,
FitParameters_t paramerrors 
) const
overridevirtualinherited

Fills the specified parameters.

Parameters
paramsthe fitted values of the parameters
paramerrorsthe uncertainty on the fitted parameters
Returns
true if the fit is valid (i.e. if a unique solution exists)

Implements lar::util::details::SimpleFitterInterface< T, D+1 >.

virtual Data_t lar::util::details::SimplePolyFitterBase< T, D >::FitParameter ( unsigned int  n) const
overridevirtualinherited

Returns the parameter n of the fit result.

Parameters
ndegree of the parameter; must be no larger than Degree
Returns
the parameter of the fit, in y/x^n units
Exceptions
std::range_errorif there is no unique solution

Reimplemented from lar::util::details::SimpleFitterInterface< T, D+1 >.

virtual FitMatrix_t lar::util::details::SimplePolyFitterBase< T, D >::FitParameterCovariance ( ) const
overridevirtualinherited

Computes and returns all the covariance matrix of the fit result.

Returns
the the covariance matrix of the fit
Exceptions
std::range_errorif there is no unique solution

The matrix is symmetric, and stored in a linear way (say, first row, then second row).

Implements lar::util::details::SimpleFitterInterface< T, D+1 >.

virtual Data_t lar::util::details::SimplePolyFitterBase< T, D >::FitParameterError ( unsigned int  n) const
overridevirtualinherited

Returns the error on parameter n of the fit result.

Parameters
ndegree of the parameter; must be no larger than Degree
Returns
the error on the parameter of the fit, in y/x^n units
Exceptions
std::range_errorif there is no unique solution

Reimplemented from lar::util::details::SimpleFitterInterface< T, D+1 >.

virtual FitParameters_t lar::util::details::SimplePolyFitterBase< T, D >::FitParameterErrors ( ) const
overridevirtualinherited

Computes and returns all the parameter errors of the fit result.

Returns
the full set of parameter errors of the fit
Exceptions
std::range_errorif there is no unique solution

Implements lar::util::details::SimpleFitterInterface< T, D+1 >.

virtual FitParameters_t lar::util::details::SimplePolyFitterBase< T, D >::FitParameterErrors ( FitMatrix_t const &  Smat) const
inlineprotectedvirtualinherited

Computes and returns all the parameter errors of the fit result.

Definition at line 785 of file SimpleFits.h.

References n.

786  {
787  return ExtractParameterErrors(Smat);
788  }
static FitParameters_t ExtractParameterErrors(FitMatrix_t const &Smat)
Extracts parameter errors from diagonal of the covarriance matrix.
virtual FitParameters_t lar::util::details::SimplePolyFitterBase< T, D >::FitParameters ( ) const
overridevirtualinherited

Computes and returns all the parameters of the fit result.

Returns
the full set of parameters of the fit
Exceptions
std::range_errorif there is no unique solution

Implements lar::util::details::SimpleFitterInterface< T, D+1 >.

virtual FitParameters_t lar::util::details::SimplePolyFitterBase< T, D >::FitParameters ( FitMatrix_t const &  Xmat) const
protectedvirtualinherited

Returns the fitted parameters using the provided information.

virtual FitParameters_t lar::util::details::SimplePolyFitterBase< T, D >::FitParameters ( FitMatrix_t const &  Smat,
Data_t   
) const
protectedvirtualinherited

Returns the fitted parameters using the provided information.

template<typename T >
Data_t lar::util::QuadraticFit< T >::I ( ) const
inlineprotected

Aliases.

Definition at line 964 of file SimpleFits.h.

964 { return Base_t::stats.template XN<0>(); }
Collector_t stats
statistics collected from fit data input
Definition: SimpleFits.h:446
virtual FitMatrix_t lar::util::details::SimpleFitterInterface< T, N >::InvertMatrix ( FitMatrix_t const &  mat,
Data_t  det 
) const
inlineprotectedvirtualinherited

Computes the inverse of a matrix (using provided determinant)

Definition at line 613 of file SimpleFits.h.

614  {
615  return MatrixOps::InvertSymmetricMatrix(mat, det);
616  }
Float_t mat
Definition: plot.C:38
virtual FitMatrix_t lar::util::details::SimpleFitterInterface< T, N >::InvertMatrix ( FitMatrix_t const &  mat) const
inlineprotectedvirtualinherited

Computes the inverse of a matrix.

Definition at line 619 of file SimpleFits.h.

620  {
621  return MatrixOps::InvertSymmetricMatrix(mat);
622  }
Float_t mat
Definition: plot.C:38
virtual bool lar::util::details::SimplePolyFitterBase< T, D >::isValid ( ) const
overridevirtualinherited

Returns if the fit has valid results.

Returns
if the fit has valid results

The fit has no valid results if:

  1. insufficient data has been add()ed (no more than the fit Degree)
  2. if input points are vertically aligned

Note that checking point 2 is expensive in terms of time.

Implements lar::util::details::SimpleFitterInterface< T, D+1 >.

virtual FitMatrix_t lar::util::details::SimplePolyFitterBase< T, D >::MakeMatrixX ( ) const
protectedvirtualinherited

Fills and returns the matrix of x^n sum coefficients ( { x^(i+j) } )

virtual FitParameters_t lar::util::details::SimplePolyFitterBase< T, D >::MakeMatrixY ( ) const
protectedvirtualinherited

Fills and returns the matrix (vector) of x^n y sum coefficients.

virtual FitParameters_t lar::util::details::SimpleFitterInterface< T, N >::MatrixProduct ( FitMatrix_t const &  mat,
FitParameters_t const &  vec 
) const
inlineprotectedvirtualinherited

Computes the product of a FitMatrix_t and a FitParameters_t.

Definition at line 625 of file SimpleFits.h.

627  {
628  return MatrixOps::MatrixVectorProduct(mat, vec);
629  }
Float_t mat
Definition: plot.C:38
template<typename T , unsigned int D>
int lar::util::details::SimplePolyFitterDataBase< T, D >::N ( ) const
inlineinherited

Definition at line 425 of file SimpleFits.h.

425 { return stats.N(); }
int N() const
Returns the number of entries added.
Definition: SimpleFits.h:269
Collector_t stats
statistics collected from fit data input
Definition: SimpleFits.h:446
virtual int lar::util::details::SimplePolyFitterBase< T, D >::NDF ( ) const
inlineoverridevirtualinherited

Returns the degrees of freedom in the determination of the fit.

Returns
the degrees of freedom in the determination of the fit

The return value may be 0 or negative if insufficient points have been added.

Implements lar::util::details::SimpleFitterInterface< T, D+1 >.

Definition at line 725 of file SimpleFits.h.

References lar::util::details::FitDataCollector< T, D >::NParams.

725 { return Base_t::N() - NParams; }
static constexpr unsigned int NParams
Number of parameters in the fit.
Definition: SimpleFits.h:647
Data_t lar::util::details::SimpleFitterInterface< T, N >::operator() ( Data_t  x) const
inlineinherited

Evaluates the fitted function; alias of Evaluate()

Definition at line 597 of file SimpleFits.h.

597 { return Evaluate(x); }
Float_t x
Definition: compare.C:6
virtual Data_t Evaluate(Data_t x) const =0
Evaluates the fitted function at the specified point.
virtual Data_t lar::util::details::SimplePolyFitterBase< T, D >::Param ( unsigned int  n,
FitMatrix_t const &  Xmat 
) const
protectedvirtualinherited

Computes a single fit parameter using the given information.

virtual Data_t lar::util::details::SimplePolyFitterBase< T, D >::Param ( unsigned int  n,
FitMatrix_t const &  Xmat,
Data_t  detXmat 
) const
protectedvirtualinherited

Computes a single fit parameter using the given information.

template<typename T , unsigned int D>
template<typename Stream >
void lar::util::details::SimplePolyFitterDataBase< T, D >::PrintStats ( Stream &  out) const
inlineinherited

Prints the collected statistics into a stream.

Definition at line 433 of file SimpleFits.h.

434  {
435  stats.Print(out);
436  }
void Print(Stream &out) const
Prints the statistics into a stream.
Definition: SimpleFits.h:1480
Collector_t stats
statistics collected from fit data input
Definition: SimpleFits.h:446
template<typename T , unsigned int D>
template<typename V >
static constexpr V lar::util::details::SimplePolyFitterDataBase< T, D >::sqr ( V const &  v)
inlinestaticinherited

Returns the square of the specified value.

Definition at line 440 of file SimpleFits.h.

References sqr().

441  {
442  return Collector_t::sqr(v);
443  }
static constexpr V sqr(V const &v)
Returns the square of the specified value.
Definition: SimpleFits.h:284
static constexpr Data_t lar::util::details::SimpleFitterInterface< T, N >::sqr ( Data_t  v)
inlinestaticinherited

Returns the square of the specified data value.

Definition at line 600 of file SimpleFits.h.

600 { return v * v; }
template<typename T >
Data_t lar::util::QuadraticFit< T >::X ( ) const
inlineprotected

Aliases.

Definition at line 965 of file SimpleFits.h.

965 { return Base_t::stats.template XN<1>(); }
Collector_t stats
statistics collected from fit data input
Definition: SimpleFits.h:446
template<typename T >
Data_t lar::util::QuadraticFit< T >::X2 ( ) const
inlineprotected

Aliases.

Definition at line 966 of file SimpleFits.h.

966 { return Base_t::stats.template XN<2>(); }
Collector_t stats
statistics collected from fit data input
Definition: SimpleFits.h:446
template<typename T >
Data_t lar::util::QuadraticFit< T >::X2Y ( ) const
inlineprotected

Aliases.

Definition at line 971 of file SimpleFits.h.

971 { return Base_t::stats.template XNY<2>(); }
Collector_t stats
statistics collected from fit data input
Definition: SimpleFits.h:446
template<typename T >
Data_t lar::util::QuadraticFit< T >::X3 ( ) const
inlineprotected

Aliases.

Definition at line 967 of file SimpleFits.h.

967 { return Base_t::stats.template XN<3>(); }
Collector_t stats
statistics collected from fit data input
Definition: SimpleFits.h:446
template<typename T >
Data_t lar::util::QuadraticFit< T >::X4 ( ) const
inlineprotected

Aliases.

Definition at line 968 of file SimpleFits.h.

968 { return Base_t::stats.template XN<4>(); }
Collector_t stats
statistics collected from fit data input
Definition: SimpleFits.h:446
template<typename T , unsigned int D>
Data_t lar::util::details::SimplePolyFitterDataBase< T, D >::XN ( unsigned int  n) const
inlineprotectedinherited

Returns the weighted sum of x^n.

Definition at line 449 of file SimpleFits.h.

References lar::util::details::FitDataCollector< T, D >::XN().

449 { return stats.XN(n); }
Data_t XN(unsigned int n) const
Returns the weighted sum of x^n.
Definition: SimpleFits.h:290
Char_t n[5]
Collector_t stats
statistics collected from fit data input
Definition: SimpleFits.h:446
template<typename T , unsigned int D>
Data_t lar::util::details::SimplePolyFitterDataBase< T, D >::XNY ( unsigned int  n) const
inlineprotectedinherited

Returns the weighted sum of x^n y.

Definition at line 452 of file SimpleFits.h.

References lar::util::details::FitDataCollector< T, D >::XNY().

452 { return stats.XNY(n); }
Data_t XNY(unsigned int n) const
Returns the weighted sum of x^n y.
Definition: SimpleFits.h:293
Char_t n[5]
Collector_t stats
statistics collected from fit data input
Definition: SimpleFits.h:446
template<typename T >
Data_t lar::util::QuadraticFit< T >::XY ( ) const
inlineprotected

Aliases.

Definition at line 970 of file SimpleFits.h.

970 { return Base_t::stats.template XNY<1>(); }
Collector_t stats
statistics collected from fit data input
Definition: SimpleFits.h:446
template<typename T >
Data_t lar::util::QuadraticFit< T >::Y ( ) const
inlineprotected

Aliases.

Definition at line 969 of file SimpleFits.h.

969 { return Base_t::stats.template XNY<0>(); }
Collector_t stats
statistics collected from fit data input
Definition: SimpleFits.h:446
template<typename T >
Data_t lar::util::QuadraticFit< T >::Y2 ( ) const
inlineprotected

Aliases.

Definition at line 972 of file SimpleFits.h.

972 { return Base_t::stats.Y2(); }
Collector_t stats
statistics collected from fit data input
Definition: SimpleFits.h:446
Data_t Y2() const
Returns the weighted sum of y^2.
Definition: SimpleFits.h:310

Member Data Documentation

constexpr unsigned int lar::util::details::SimplePolyFitterBase< T, D >::Degree
staticinherited

Degree of the fit.

Definition at line 644 of file SimpleFits.h.

constexpr unsigned int lar::util::details::SimplePolyFitterBase< T, D >::NParams
staticinherited

Number of parameters in the fit.

Definition at line 647 of file SimpleFits.h.

template<typename T , unsigned int D>
Collector_t lar::util::details::SimplePolyFitterDataBase< T, D >::stats
protectedinherited

statistics collected from fit data input

Definition at line 446 of file SimpleFits.h.


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