QMathFunction.hh

Go to the documentation of this file.

#ifndef _Q_MATH_FUNCTION_
#define _Q_MATH_FUNCTION_

#include <gsl/gsl_math.h>
#include "QCuore.hh"
#include "QError.hh"
#include <vector>

class TF1;

Q_BEGIN_NAMESPACE

class QMathFunction 
{
    public:
        typedef double(*QRoutinePointer)(double,double* params);
        
        QMathFunction(QRoutinePointer funcpointer, unsigned int NParameters, double LowBound, double UpBound);

        QMathFunction(const QMathFunction& func);
        
        virtual ~QMathFunction();

        const QMathFunction& operator=(const QMathFunction& func);
        
        void SetParameter(unsigned int p, double val);
        
        void SetParameter(unsigned int p, const char* name, double val);
        
        void SetParameters(const double* param) {memcpy(fParameters,param,sizeof(double)*fNParameters); }

        double GetParameter(unsigned int p) const {return fParameters[p]; }

        const double* GetParameters() const {return fParameters; }
        
        size_t GetNParameters() const {return fNParameters; }
        
        double Eval(double x) const { return (*fRoutine)(x,fParameters);}
        
        double Eval(double x, double* params) const { return (*fRoutine)(x,params);}
        
        virtual void Print();
        
        void SetVerbose(bool verbose = true) { fVerbose = verbose; }

        virtual void SetTF1(TF1& tf1);

        QError GetMinimum(double& min, double &minpos, double epsabs = 0.001, double epsrel = 0.0);

        QError GetMaximum(double& max, double &maxpos, double epsabs = 0.001, double epsrel = 0.0);

    protected:
        gsl_function GetGSLFunction();

        void SetParName(unsigned int  p,const char* name) { if(p < fNParameters) fParNames[p] = name; }         
        QRoutinePointer fRoutine;

        size_t fNParameters;
        
        double* fParameters;

        std::vector<const char*> fParNames;

        double fLowBound;
        double fUpBound;

        bool fVerbose;
        
        friend class QMinimizer;
};

Q_END_NAMESPACE

#endif

---