#include <Landau.h>
Inheritance diagram for Landau:
This class represents a Landau distribution, as approximated by the Moyal formula
See J.E. Moyal, Theory of ionization fluctuations, Phil. Mag. 46 (1955) 263.
Note that this analytical approximation is too low in the tail.
In order to allow for a fit, we define
with x the dataset variable. and are fit parameters PAR(1) and PAR(2) respectively, and correspond to the mode and the width of the distribution. PAR(3) is a normalization factor, so that the concrete implemented formula is: .
Definition at line 47 of file Landau.h.
Public Member Functions | |
virtual void | addToComposite (FunctionBase *) |
Does nothing. | |
virtual FunctionBase * | clone () const |
Creates a new function object by copying an existing one. | |
virtual int | count () |
Returns 0. | |
virtual unsigned int | dimensions () const |
Returns the number of dimensions of the data coordinate. | |
virtual const std::vector< double > & | getParameters () const |
Returns the values of the parameters as a vector. | |
virtual bool | hasDerivatives () const |
Returns true if function can calculate its partial derivatives. | |
virtual void | initialParameters (const FunctionHelper *helper) |
Sets the FunctionHelper so that the function can calculate a reasonable set of initial parameter values. | |
virtual double | integrate (double lower_limit, double upper_limit) const |
virtual bool | isComposite () const |
Returns false . | |
Landau (double peak, double counts, double sigma) | |
This constructor takes the values of the parameters as an argument. | |
Landau () | |
The default constructor. | |
const std::string & | name () const |
Returns the name of the function. | |
virtual double | operator() (const std::vector< double > &v) const |
The function call operator for multi-dimension coordinate variable. | |
virtual double | operator() (double x) const |
The original Fortran. | |
virtual const std::vector< std::string > & | parmNames () const |
Returns a reference to a vector of parameter names. | |
virtual void | removeFromComposite (FunctionBase *) |
Does nothing. | |
virtual std::vector< double >::const_iterator | setParameters (std::vector< double >::const_iterator it) |
Sets the parameter values to the value pointed to by the iterator. | |
virtual void | setParameters (const std::vector< double > &incr) |
Sets the parameter values. | |
void | setParmNames (const std::vector< std::string > &names) |
Sets the names of the parameters. | |
virtual int | size () const |
Returns the number of parameters. | |
Protected Member Functions | |
virtual double | derivByParm (int i, double x) const |
Returns the function's derivative at the coordinate value x with respect to the i-th parameter. | |
virtual void | initialize () |
Initializes the function and parameter names. | |
virtual void | resize () |
Re-sizes the appropriate vectors maintained in this base class. | |
void | setName (const char *) |
Sets the name of the function. | |
Protected Attributes | |
std::string | m_name |
The name of the function. | |
std::vector< std::string > | m_parm_names |
The names of the function parameters. | |
std::vector< double > | m_parms |
The parameter values. | |
Private Types | |
enum | { peak, norm, sigma } |
Indices to the parameter vectors. More... | |
Private Member Functions | |
double | calcY (double x) const |
compute the reduced variable: | |
double | calcZ (double x) const |
Compute the helper variable for derivation: . | |
double | derivByNorm (double x) const |
Calculates the derivative with respect to the counts parameter. | |
double | derivByPeak (double x) const |
Calculates the derivative with respect to the peak parameter. | |
double | derivBySigma (double x) const |
Calculates the derivative with respect to the sigma parameter. |
anonymous enum [private] |
Landau | ( | ) |
The default constructor.
Definition at line 34 of file Landau.cxx.
References Landau::initialize().
Referenced by Landau::clone().
Landau | ( | double | peak, | |
double | counts, | |||
double | sigma | |||
) |
This constructor takes the values of the parameters as an argument.
Definition at line 39 of file Landau.cxx.
References Landau::initialize(), FunctionBase::m_parms, Landau::norm, Landau::peak, and Landau::sigma.
void addToComposite | ( | FunctionBase * | ) | [virtual, inherited] |
Does nothing.
Derived classes that are a composite of functions should override this member function and add functions to the composite.
Reimplemented in LinearSumFunction.
Definition at line 166 of file FunctionBase.cxx.
double calcY | ( | double | x | ) | const [inline, private] |
compute the reduced variable:
Definition at line 84 of file Landau.h.
References sigma.
Referenced by Landau::operator()().
double calcZ | ( | double | x | ) | const [inline, private] |
Compute the helper variable for derivation: .
Definition at line 92 of file Landau.h.
Referenced by Landau::derivByPeak(), and Landau::derivBySigma().
FunctionBase * clone | ( | ) | const [virtual] |
Creates a new function object by copying an existing one.
The implementation raises an assertion as this member function should be reimplemented in s derived class.
Reimplemented from FunctionBase.
Definition at line 59 of file Landau.cxx.
References Landau::Landau().
int count | ( | ) | [virtual, inherited] |
Returns 0.
Derived classes that are a composite of functions should override this member function to return the number of functions in the composite.
Reimplemented in LinearSumFunction.
Definition at line 174 of file FunctionBase.cxx.
double derivByNorm | ( | double | x | ) | const [private] |
Calculates the derivative with respect to the counts parameter.
Definition at line 122 of file Landau.cxx.
References FunctionBase::m_parms, Landau::norm, and Landau::operator()().
Referenced by Landau::derivByParm().
double derivByParm | ( | int | i, | |
double | x | |||
) | const [protected, virtual] |
Returns the function's derivative at the coordinate value x
with respect to the i-th
parameter.
Derived classes must implement this function.
Reimplemented from FunctionBase.
Definition at line 100 of file Landau.cxx.
References Landau::derivByNorm(), Landau::derivByPeak(), Landau::derivBySigma(), Landau::norm, Landau::peak, and Landau::sigma.
double derivByPeak | ( | double | x | ) | const [private] |
Calculates the derivative with respect to the peak parameter.
Definition at line 129 of file Landau.cxx.
References Landau::calcZ(), FunctionBase::m_parms, Landau::operator()(), and Landau::sigma.
Referenced by Landau::derivByParm().
double derivBySigma | ( | double | x | ) | const [private] |
Calculates the derivative with respect to the sigma parameter.
Definition at line 134 of file Landau.cxx.
References Landau::calcZ(), FunctionBase::m_parms, Landau::operator()(), Landau::peak, and Landau::sigma.
Referenced by Landau::derivByParm().
unsigned int dimensions | ( | ) | const [virtual, inherited] |
Returns the number of dimensions of the data coordinate.
This implementation in the base class returns 1
. Derived class should override this function if the data coordinate is multi-dimensional.
Reimplemented in FunctionWrap.
Definition at line 213 of file FunctionBase.cxx.
Referenced by FunctionWrap::dimensions().
const vector< double > & getParameters | ( | ) | const [virtual, inherited] |
Returns the values of the parameters as a vector.
The function is non-const because a derived class may need to create the vector.
Reimplemented in LinearSumFunction.
Definition at line 99 of file FunctionBase.cxx.
References FunctionBase::m_parms.
Referenced by FunctionProjectorXML::createElement(), FuncParmRep::drawProjectedValues(), hippodraw::Python::export_FunctionBase(), FunctionRep::parameters(), FunctionProjector::saveParameters(), and FunctionProjector::setParameters().
bool hasDerivatives | ( | ) | const [virtual, inherited] |
Returns true
if function can calculate its partial derivatives.
This base class implementation returns true
. Derived class that can not calculate its partial derivatives should return false
. This allows an application to not attempt to use a minimizer that requires derivatives to use this function.
Reimplemented in LogNormal, LogParabola, Novosibirsk, and FunctionWrap.
Definition at line 126 of file FunctionBase.cxx.
void initialize | ( | ) | [protected, virtual] |
Initializes the function and parameter names.
Reimplemented from FunctionBase.
Definition at line 48 of file Landau.cxx.
References FunctionBase::m_name, FunctionBase::m_parm_names, and FunctionBase::resize().
Referenced by Landau::Landau().
void initialParameters | ( | const FunctionHelper * | helper | ) | [virtual] |
Sets the FunctionHelper so that the function can calculate a reasonable set of initial parameter values.
Reimplemented from FunctionBase.
Definition at line 93 of file Landau.cxx.
References FunctionBase::m_parms, FunctionHelper::maxValue(), FunctionHelper::meanCoord(), Landau::norm, Landau::peak, Landau::sigma, and FunctionHelper::stdCoord().
double integrate | ( | double | lower_limit, | |
double | upper_limit | |||
) | const [virtual, inherited] |
As of now all the integration is done here numerically using the trapezoidal rule. Derived class may implement an analytical expression in place of this numerical method.
Definition at line 138 of file FunctionBase.cxx.
References FunctionBase::operator()().
bool isComposite | ( | ) | const [virtual, inherited] |
Returns false
.
Derived classes that are a composite of functions should override this member function and return true
.
Reimplemented in LinearSumFunction.
Definition at line 161 of file FunctionBase.cxx.
Referenced by FunctionRep::isComposite(), FunctionController::restoreParameters(), and FunctionController::saveParameters().
const string & name | ( | ) | const [inherited] |
Returns the name of the function.
Definition at line 81 of file FunctionBase.cxx.
References FunctionBase::m_name.
Referenced by FunctionProjectorXML::createElement(), FunctionController::createNTuple(), hippodraw::Python::export_FunctionBase(), FunctionController::functionNames(), FunctionProjector::getTitle(), FunctionProjector::getYLabel(), and LinearSumFunction::parmNames().
double operator() | ( | const std::vector< double > & | v | ) | const [virtual, inherited] |
The function call operator for multi-dimension coordinate variable.
Returns the value of the function at coordinate represented by v. The implementation in this base class uses the first element of the vector and calls operator () ( double ).
Reimplemented in FunctionWrap.
Definition at line 181 of file FunctionBase.cxx.
References FunctionBase::operator()().
double operator() | ( | double | x | ) | const [virtual] |
The original Fortran.
Reimplemented from FunctionBase.
Definition at line 81 of file Landau.cxx.
References Landau::calcY(), FunctionBase::m_parms, and Landau::norm.
Referenced by Landau::derivByNorm(), Landau::derivByPeak(), and Landau::derivBySigma().
const vector< string > & parmNames | ( | ) | const [virtual, inherited] |
Returns a reference to a vector of parameter names.
Reimplemented in LinearSumFunction.
Definition at line 86 of file FunctionBase.cxx.
References FunctionBase::m_parm_names.
Referenced by FuncParmRep::drawProjectedValues(), hippodraw::Python::export_FunctionBase(), and FunctionRep::parmNames().
void removeFromComposite | ( | FunctionBase * | ) | [virtual, inherited] |
Does nothing.
Derived classes that are a composite of functions should override this member function and remove functions from the composite.
Reimplemented in LinearSumFunction.
Definition at line 170 of file FunctionBase.cxx.
void resize | ( | ) | [protected, virtual, inherited] |
Re-sizes the appropriate vectors maintained in this base class.
Derived classes should call this member function after initializing the m_parm_names data member, as the new size is taken from that data member's size.
Definition at line 71 of file FunctionBase.cxx.
References FunctionBase::m_parm_names, FunctionBase::m_parms, and FunctionBase::size().
Referenced by Quadratic2::initialize(), Quadratic::initialize(), PowerLaw::initialize(), Novosibirsk::initialize(), LogParabola::initialize(), LogNormal::initialize(), Linear::initialize(), Landau::initialize(), Gaussian::initialize(), FunctionWrap::initialize(), Exponential::initialize(), Erfc::initialize(), ConstantF::initialize(), Chi2Dist::initialize(), BrokenPowerLaw::initialize(), FunctionWrap::setParmNames(), and FunctionBase::setParmNames().
void setName | ( | const char * | ) | [protected, inherited] |
Sets the name of the function.
Definition at line 66 of file FunctionBase.cxx.
References FunctionBase::m_name.
Referenced by FunctionWrap::setName().
vector< double >::const_iterator setParameters | ( | std::vector< double >::const_iterator | it | ) | [virtual, inherited] |
Sets the parameter values to the value pointed to by the iterator.
Returns an iterator to a data element one past the last used data element.
Reimplemented in LinearSumFunction.
Definition at line 114 of file FunctionBase.cxx.
References FunctionBase::m_parms, and FunctionBase::size().
void setParameters | ( | const std::vector< double > & | incr | ) | [virtual, inherited] |
Sets the parameter values.
Derived classes should implement this function by invoking incrementValues ( vector< double >::const_iterator ) so that they correctly work with the LinearSumFunction class.
Definition at line 106 of file FunctionBase.cxx.
Referenced by hippodraw::Python::export_FunctionBase(), FunctionProjector::restoreParameters(), FunctionProjectorXML::setAttributes(), and FunctionProjector::setParameters().
void setParmNames | ( | const std::vector< std::string > & | names | ) | [inherited] |
Sets the names of the parameters.
Reimplemented in FunctionWrap.
Definition at line 93 of file FunctionBase.cxx.
References FunctionBase::m_parm_names, and FunctionBase::resize().
int size | ( | ) | const [virtual, inherited] |
Returns the number of parameters.
Reimplemented in LinearSumFunction.
Definition at line 156 of file FunctionBase.cxx.
References FunctionBase::m_parm_names.
Referenced by Novosibirsk::initialParameters(), LogParabola::initialParameters(), Gaussian::initialParameters(), FunctionWrap::operator()(), FunctionBase::resize(), FunctionProjectorXML::setAttributes(), and FunctionBase::setParameters().
std::string m_name [protected, inherited] |
The name of the function.
Definition at line 103 of file FunctionBase.h.
Referenced by FunctionBase::FunctionBase(), Quadratic2::initialize(), Quadratic::initialize(), PowerLaw::initialize(), Novosibirsk::initialize(), LogParabola::initialize(), LogNormal::initialize(), LinearSumFunction::initialize(), Linear::initialize(), Landau::initialize(), Gaussian::initialize(), Exponential::initialize(), Erfc::initialize(), ConstantF::initialize(), Chi2Dist::initialize(), BrokenPowerLaw::initialize(), FunctionBase::name(), and FunctionBase::setName().
std::vector<std::string> m_parm_names [protected, inherited] |
The names of the function parameters.
Definition at line 106 of file FunctionBase.h.
Referenced by Quadratic2::initialize(), Quadratic::initialize(), PowerLaw::initialize(), Novosibirsk::initialize(), LogParabola::initialize(), LogNormal::initialize(), Linear::initialize(), Landau::initialize(), Gaussian::initialize(), Exponential::initialize(), Erfc::initialize(), ConstantF::initialize(), Chi2Dist::initialize(), BrokenPowerLaw::initialize(), LinearSumFunction::parmNames(), FunctionBase::parmNames(), FunctionBase::resize(), FunctionWrap::setParmNames(), FunctionBase::setParmNames(), and FunctionBase::size().
std::vector< double > m_parms [protected, inherited] |
The parameter values.
Definition at line 109 of file FunctionBase.h.
Referenced by BrokenPowerLaw::BrokenPowerLaw(), Erfc::calcRed(), Chi2Dist::Chi2Dist(), ConstantF::ConstantF(), Gaussian::derivByMean(), Landau::derivByNorm(), Gaussian::derivByNorm(), Quadratic2::derivByParm(), PowerLaw::derivByParm(), LogParabola::derivByParm(), Exponential::derivByParm(), Erfc::derivByParm(), Chi2Dist::derivByParm(), BrokenPowerLaw::derivByParm(), Landau::derivByPeak(), Landau::derivBySigma(), Gaussian::derivBySigma(), Erfc::Erfc(), Exponential::Exponential(), Gaussian::Gaussian(), LinearSumFunction::getParameters(), FunctionBase::getParameters(), Quadratic2::initialParameters(), Quadratic::initialParameters(), PowerLaw::initialParameters(), Novosibirsk::initialParameters(), LogParabola::initialParameters(), LogNormal::initialParameters(), Linear::initialParameters(), Landau::initialParameters(), Gaussian::initialParameters(), Exponential::initialParameters(), Erfc::initialParameters(), ConstantF::initialParameters(), Chi2Dist::initialParameters(), BrokenPowerLaw::initialParameters(), Landau::Landau(), Linear::Linear(), LogNormal::LogNormal(), LogParabola::LogParabola(), Novosibirsk::Novosibirsk(), Quadratic2::operator()(), Quadratic::operator()(), PowerLaw::operator()(), Novosibirsk::operator()(), LogParabola::operator()(), LogNormal::operator()(), Linear::operator()(), Landau::operator()(), Gaussian::operator()(), Exponential::operator()(), ConstantF::operator()(), Chi2Dist::operator()(), BrokenPowerLaw::operator()(), PowerLaw::PowerLaw(), Quadratic::Quadratic(), Quadratic2::Quadratic2(), FunctionBase::resize(), and FunctionBase::setParameters().