Performs a linear regression of data. More...

#include "SimpleFits.h"

Inheritance diagram for lar::util::LinearFit< T >:

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
Data_t	Intercept () const
	Returns the intercept of the fit. More...

Data_t	Slope () const
	Returns the slope of the fit. More...

Data_t	InterceptError () const
	Returns the error on intercept of the fit. More...

Data_t	SlopeError () const
	Returns the error in slope of the fit. More...

Data_t	InterceptSlopeCovariance () const
	Returns the covariance between intercept and slope of the fit. More...

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	Y () const
	Aliases. More...

Data_t	XY () 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, 1U >

Detailed Description

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

Performs a linear regression of data.

Template Parameters

T	type of the quantities
W	type of the weight (as T by default)

The linear regression connects measurements $( y_{i} \pm \sigma_{y,i} )$ with a parameter $( x_{i} )$ not affected by uncertainty. The returned parameters describe a straight line $ y = a x + b $ obtained by minimization of $\chi^{2} = \sum_{i} \frac{ \left(y_{i} - a x_{i} - b \right)^{2} }{ \sigma^{2}_{y,i} }$

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

This simple linear 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 833 of file SimpleFits.h.

Member Typedef Documentation

template<typename T>

using lar::util::LinearFit< T >::Base_t = details::SimplePolyFitterBase<T, 1U>

private

Definition at line 834 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::LinearFit< T >::FitMatrix_t = std::array<Data_t, sqr(NParams)>

Definition at line 851 of file SimpleFits.h.

template<typename T>

using lar::util::LinearFit< T >::FitParameters_t = std::array<Data_t, NParams>

type of set of fit parameters

Definition at line 850 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); }

x

Float_t x

Definition: compare.C:6

y

Float_t y

Definition: compare.C:6

lar::util::details::SimplePolyFitterDataBase::stats

Collector_t stats

statistics collected from fit data input

Definition: SimpleFits.h:446

lar::util::details::FitDataCollector::add

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); }

value

double value

Definition: spectrum.C:18

lar::util::details::SimplePolyFitterDataBase::stats

Collector_t stats

statistics collected from fit data input

Definition: SimpleFits.h:446

lar::util::details::FitDataCollector::add

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); }

value

double value

Definition: spectrum.C:18

lar::util::details::SimplePolyFitterDataBase::stats

Collector_t stats

statistics collected from fit data input

Definition: SimpleFits.h:446

lar::util::details::FitDataCollector::add

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.

         {
           return stats.add_with_uncertainty(
             begin_value, end_value, begin_uncertainty, value_extractor, uncertainty_extractor);
         }

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.

         {
           return stats.add_with_uncertainty(begin, end);
         }

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.

         {
           return stats.add_with_uncertainty(cont);
         }

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.

         {
           stats.add_without_uncertainty(begin, end);
         }

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.

         {
           stats.add_without_uncertainty(begin, end, extractor);
         }

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.

         {
           stats.add_without_uncertainty(cont, extractor);
         }

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.

         {
           stats.add_without_uncertainty(cont);
         }

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(); }

lar::util::details::FitDataCollector::AverageUncertainty

Data_t AverageUncertainty() const

Returns an average of the uncertainties.

Definition: SimpleFits.h:280

lar::util::details::SimplePolyFitterDataBase::stats

Collector_t stats

statistics collected from fit data input

Definition: SimpleFits.h:446

template<typename T >

auto lar::util::LinearFit< 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, 1U >.

Definition at line 1690 of file SimpleFits.h.

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

 {
   FitParameters_t fit_params = this->FitParameters();
   const Data_t b = fit_params[0];
   const Data_t a = fit_params[1];
   return Y2() + sqr(a) * X2() + sqr(b) * I() + Data_t(2) * (a * b * X2() - a * XY() - b * Y());
 } // LinearFit<T>::ChiSquare()

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(); }

lar::util::details::FitDataCollector::clear

void clear()

Clears all the statistics.

Definition: SimpleFits.h:1469

lar::util::details::SimplePolyFitterDataBase::stats

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.

         {
           return MatrixOps::Determinant(mat);
         }

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

x	the 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

x	the 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

params	the fitted values of the parameters
Xmat	the matrix of the x^n/s^2 sums
Smat	the covariance matrix
det	the 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

params	the fitted values of the parameters
paramerrors	the uncertainty on the fitted parameters
Xmat	the matrix of the x^n/s^2 sums
Smat	the covariance matrix
det	the 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

params	the fitted values of the parameters
Xmat	the matrix of the x^n/s^2 sums
Smat	the covariance matrix
det	the 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

params	the fitted values of the parameters
paramerrors	the uncertainty on the fitted parameters
Xmat	the matrix of the x^n/s^2 sums
Smat	the covariance matrix
det	the 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

params	the fitted values of the parameters
paramerrors	the 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

n	degree of the parameter; must be no larger than Degree

Returns: the parameter of the fit, in y/x^n units

Exceptions

std::range_error if 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_error if 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

n	degree 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_error if 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_error if 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.

         {
           return ExtractParameterErrors(Smat);
         }

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_error if there is no unique solution

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

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

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::LinearFit< T >::I ( ) const

inlineprotected

Aliases.

Definition at line 902 of file SimpleFits.h.

902 { return Base_t::stats.template XN<0>(); }

lar::util::details::SimplePolyFitterDataBase::stats

Collector_t stats

statistics collected from fit data input

Definition: SimpleFits.h:446

template<typename T>

Data_t lar::util::LinearFit< T >::Intercept ( ) const

inline

Returns the intercept of the fit.

Returns: the intercept of the fit, in y units

Exceptions

std::range_error if there is no unique solution

Definition at line 860 of file SimpleFits.h.

860 { return this->FitParameter(0); }

lar::util::details::SimplePolyFitterBase< T, 1U >::FitParameter

virtual Data_t FitParameter(unsigned int n) const override

Returns the parameter n of the fit result.

template<typename T>

Data_t lar::util::LinearFit< T >::InterceptError ( ) const

inline

Returns the error on intercept of the fit.

Returns: the error on intercept of the fit, in y units

Exceptions

std::range_error if there is no unique solution

Definition at line 874 of file SimpleFits.h.

874 { return this->FitParameterError(0); }

lar::util::details::SimplePolyFitterBase< T, 1U >::FitParameterError

virtual Data_t FitParameterError(unsigned int n) const override

Returns the error on parameter n of the fit result.

template<typename T>

Data_t lar::util::LinearFit< T >::InterceptSlopeCovariance ( ) const

inline

Returns the covariance between intercept and slope of the fit.

Returns: the covariance between intercept and slope of the fit, in y^2 units

Exceptions

std::range_error if there is no unique solution

Definition at line 888 of file SimpleFits.h.

       {
         return this->FitParameterCovariance()[0 * NParams + 1];
       }

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.

         {
           return MatrixOps::InvertSymmetricMatrix(mat, det);
         }

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.

         {
           return MatrixOps::InvertSymmetricMatrix(mat);
         }

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:

insufficient data has been add()ed (no more than the fit Degree)
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.

         {
           return MatrixOps::MatrixVectorProduct(mat, vec);
         }

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(); }

lar::util::details::FitDataCollector::N

int N() const

Returns the number of entries added.

Definition: SimpleFits.h:269

lar::util::details::SimplePolyFitterDataBase::stats

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.

725 { return Base_t::N() - NParams; }

lar::util::details::SimplePolyFitterBase< T, 1U >::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); }

x

Float_t x

Definition: compare.C:6

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

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.

         {
           stats.Print(out);
         }

template<typename T>

Data_t lar::util::LinearFit< T >::Slope ( ) const

inline

Returns the slope of the fit.

Returns: the slope of the fit, in y/x units

Exceptions

std::range_error if there is no unique solution

Definition at line 867 of file SimpleFits.h.

867 { return this->FitParameter(1); }

lar::util::details::SimplePolyFitterBase< T, 1U >::FitParameter

virtual Data_t FitParameter(unsigned int n) const override

Returns the parameter n of the fit result.

template<typename T>

Data_t lar::util::LinearFit< T >::SlopeError ( ) const

inline

Returns the error in slope of the fit.

Returns: the error on slope of the fit, in y/x units

Exceptions

std::range_error if there is no unique solution

Definition at line 881 of file SimpleFits.h.

881 { return this->FitParameterError(1); }

lar::util::details::SimplePolyFitterBase< T, 1U >::FitParameterError

virtual Data_t FitParameterError(unsigned int n) const override

Returns the error on parameter n of the fit result.

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().

         {
           return Collector_t::sqr(v);
         }

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::LinearFit< T >::X ( ) const

inlineprotected

Aliases.

Definition at line 903 of file SimpleFits.h.

903 { return Base_t::stats.template XN<1>(); }

lar::util::details::SimplePolyFitterDataBase::stats

Collector_t stats

statistics collected from fit data input

Definition: SimpleFits.h:446

template<typename T>

Data_t lar::util::LinearFit< T >::X2 ( ) const

inlineprotected

Aliases.

Definition at line 904 of file SimpleFits.h.

904 { return Base_t::stats.template XN<2>(); }

lar::util::details::SimplePolyFitterDataBase::stats

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); }

lar::util::details::FitDataCollector::XN

Data_t XN(unsigned int n) const

Returns the weighted sum of x^n.

Definition: SimpleFits.h:290

n

Char_t n[5]

Definition: comparison_ascii.C:55

lar::util::details::SimplePolyFitterDataBase::stats

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); }

lar::util::details::FitDataCollector::XNY

Data_t XNY(unsigned int n) const

Returns the weighted sum of x^n y.

Definition: SimpleFits.h:293

n

Char_t n[5]

Definition: comparison_ascii.C:55

lar::util::details::SimplePolyFitterDataBase::stats

Collector_t stats

statistics collected from fit data input

Definition: SimpleFits.h:446

template<typename T>

Data_t lar::util::LinearFit< T >::XY ( ) const

inlineprotected

Aliases.

Definition at line 906 of file SimpleFits.h.

906 { return Base_t::stats.template XNY<1>(); }

lar::util::details::SimplePolyFitterDataBase::stats

Collector_t stats

statistics collected from fit data input

Definition: SimpleFits.h:446

template<typename T>

Data_t lar::util::LinearFit< T >::Y ( ) const

inlineprotected

Aliases.

Definition at line 905 of file SimpleFits.h.

905 { return Base_t::stats.template XNY<0>(); }

lar::util::details::SimplePolyFitterDataBase::stats

Collector_t stats

statistics collected from fit data input

Definition: SimpleFits.h:446

template<typename T>

Data_t lar::util::LinearFit< T >::Y2 ( ) const

inlineprotected

Aliases.

Definition at line 907 of file SimpleFits.h.

907 { return Base_t::stats.Y2(); }

lar::util::details::SimplePolyFitterDataBase::stats

Collector_t stats

statistics collected from fit data input

Definition: SimpleFits.h:446

lar::util::details::FitDataCollector::Y2

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:

lardata/v09_16_03/source/lardata/Utilities/SimpleFits.h

Public Types

Public Member Functions

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Protected Attributes

Private Types

Detailed Description

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

Member Typedef Documentation

Member Function Documentation

Member Data Documentation

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