QMathFunction.cc

Go to the documentation of this file.

#include "QMathFunction.hh"
#include "QMinimizer.hh"
#include  "QVector.hh"
#include "TF1.h"
 
Q_BEGIN_NAMESPACE
 
gsl_function QMathFunction::GetGSLFunction()
{
    gsl_function F;
    F.function =  (double (*)(double, void*))fRoutine;
    F.params = (void*) fParameters;
    return F;
}
 
QMathFunction::QMathFunction(QRoutinePointer funcpointer, unsigned int NParameters, double LowBound, double UpBound)
{
    fRoutine = funcpointer;
    fNParameters = NParameters;
    fParameters = new double[NParameters];
    fParNames.resize(fNParameters);
    for(size_t  p = 0; p < fNParameters; p++) {
        fParameters[p] = 0.0;
        fParNames[p] = "";
    }
    fLowBound = LowBound;
    fUpBound = UpBound;
    fVerbose = false;
}
 
QMathFunction::QMathFunction(const QMathFunction& func)
{
    fNParameters = 0;
    fParameters = NULL;
    *this = func;
}
 
QMathFunction::~QMathFunction()
{
    if(fParameters) delete [] fParameters;
}
 
const QMathFunction& QMathFunction::operator=(const QMathFunction& func)
{
    fRoutine = func.fRoutine;
    if(fNParameters!=func.fNParameters) {
        if(fParameters) delete [] fParameters;
        fNParameters = func.fNParameters;
        fParameters = new double[fNParameters];
        fParNames.resize(fNParameters);
    }
    for(unsigned int p = 0; p < fNParameters; p++) {
        fParameters[p] = func.fParameters[p];
        fParNames[p] = func.fParNames[p];
    }
    fLowBound = func.fLowBound;
    fUpBound = func.fUpBound;
    fVerbose = false;
    return *this;
 
}
 
void QMathFunction::Print()
{
    for(size_t  p = 0; p < fNParameters; p++) {
        if(std::string(fParNames[p]) != "") 
            printf("%s ",fParNames[p]);
        else
            printf("%u ",(unsigned int)p);
        printf("%f\n", fParameters[p]);
    }
 
}
 
void QMathFunction::SetParameter(unsigned int p, double val)
{
    if(p < fNParameters) {
        fParameters[p] = val;
    }
}
 
void QMathFunction::SetParameter(unsigned int p, const char* name, double val)
{
    SetParameter(p,val);
    SetParName(p,name);
}
 
 
QError QMathFunction::GetMinimum(double& min, double& minpos, double epsabs,double epsrel)
{
    QError err;
    QMinimizer* minimizer = new QMinimizer();
    if(fVerbose) minimizer->SetVerbose();
    minpos = minimizer->GetMinimum(*this, fLowBound, fUpBound,epsabs,epsrel);
    if(minpos == Q_DOUBLE_DEFAULT) {
        min = Q_DOUBLE_DEFAULT;
        err = QERR_UNKNOWN_ERR;
    } else {
        min = (*fRoutine)(minpos,fParameters);
        err = QERR_SUCCESS;
    }
    delete minimizer;
    return err;
}
 
QError QMathFunction::GetMaximum(double& max, double& maxpos, double epsabs,double epsrel)
{
    QError err;
    QMinimizer* minimizer = new QMinimizer();
    if(fVerbose) minimizer->SetVerbose();
    maxpos = minimizer->GetMaximum(*this, fLowBound, fUpBound,epsabs,epsrel);
    if(maxpos == Q_DOUBLE_DEFAULT) {
        max = Q_DOUBLE_DEFAULT;
        err = QERR_UNKNOWN_ERR;
    } else {
        max = (*fRoutine)(maxpos,fParameters);
        err = QERR_SUCCESS;
    }
    delete minimizer;
    return err;
}
 
void QMathFunction::SetTF1(TF1& tf1)
{
    for(size_t p = 0; p < fNParameters; p++) {
        tf1.SetParameter(p,fParameters[p]);
        tf1.SetParName(p,fParNames[p]);
    }
}
 
Q_END_NAMESPACE