QMathFunction.hh

Go to the documentation of this file.

#ifndef _Q_MATH_FUNCTION_
#define _Q_MATH_FUNCTION_
#include <gsl/gsl_math.h>
#include "QDiana.hh"
#include "QError.hh"
#include <vector>
 
class TF1;
 
Q_BEGIN_NAMESPACE
 
class QMathFunction 
{
    public:
        #ifndef __CINT__
        typedef double(*QRoutinePointer)(double,double* params);
        
        QMathFunction(QRoutinePointer funcpointer, unsigned int NParameters, double LowBound, double UpBound);
        #endif
 
        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; }         
        #ifndef __CINT__
        QRoutinePointer fRoutine;
        #endif
 
        size_t fNParameters;
        
        double* fParameters;
 
        std::vector<const char*> fParNames;
 
        double fLowBound;
        double fUpBound;
 
        bool fVerbose;
        
        friend class QMinimizer;
};
 
Q_END_NAMESPACE
 
#endif