// ---------------------------------- F1.h
#ifndef ROOT_TF1
#define ROOT_TF1
//+SEQ,CopyRight,T=NOINCLUDE.
//////////////////////////////////////////////////////////////////////////
// //
// TF1 //
// //
// The Parametric 1-D function //
// //
//////////////////////////////////////////////////////////////////////////
#ifndef ROOT_TFormula
//*KEEP,TFormula.
#include "TFormula.h"
//*KEND.
#endif
#ifndef ROOT_TAttLine
//*KEEP,TAttLine.
#include "TAttLine.h"
//*KEND.
#endif
#ifndef ROOT_TAttFill
//*KEEP,TAttFill.
#include "TAttFill.h"
//*KEND.
#endif
#ifndef ROOT_TAttMarker
//*KEEP,TAttMarker.
#include "TAttMarker.h"
//*KEND.
#endif
#ifndef ROOT_TMethodCall
//*KEEP,TMethodCall,T=C++.
#include "TMethodCall.h"
//*KEND.
#endif
class TH1F;
class TF1;
class TF1 : public TFormula, public TAttLine, public TAttFill, public TAttMarker {
protected:
Float_t fXmin; //Lower bounds for the range
Float_t fXmax; //Upper bounds for the range
Int_t fNpx; //Number of points used for the graphical representation
Int_t fType; //(=0 for standard functions, 1 if pointer to function)
Int_t fNpfits; //Number of points used in the fit
Int_t fNsave; //Number of points used to fill array fSave
Double_t fChisquare; //Function fit chisquare
Double_t *fIntegral; //Integral of function binned on fNpx bins
Double_t *fParErrors; //Array of errors of the fNpar parameters
Double_t *fParMin; //Array of lower limits of the fNpar parameters
Double_t *fParMax; //Array of upper limits of the fNpar parameters
Double_t *fSave; //Array of fNsave function values
TObject *fParent; //Parent object hooking this function (if one)
TH1F *fHistogram; //Pointer to histogram used for visualisation
Float_t fMaximum; //Maximum value for plotting
Float_t fMinimum; //Minimum value for plotting
TMethodCall *fMethodCall; //Pointer to MethodCall in case of interpreted function
Double_t (*fFunction) (Double_t *, Double_t *); //Pointer to function
public:
TF1();
TF1(const char *name, const char *formula, Float_t xmin=0, Float_t xmax=1);
TF1(const char *name, void *fcn, Float_t xmin, Float_t xmax, Int_t npar);
TF1(const char *name, Double_t (*fcn)(Double_t *, Double_t *), Float_t xmin=0, Float_t xmax=1, Int_t npar=0);
TF1(const TF1 &f1);
virtual ~TF1();
virtual void Browse(TBrowser *b);
virtual void Copy(TObject &f1);
virtual Double_t Derivative(Double_t x, Double_t *params=0, Float_t epsilon=0);
virtual Int_t DistancetoPrimitive(Int_t px, Int_t py);
virtual void Draw(Option_t *option="");
virtual TF1 *DrawCopy(Option_t *option="");
virtual void DrawF1(const char *formula, Float_t xmin, Float_t xmax, Option_t *option="");
virtual void DrawPanel(); // *MENU*
virtual Double_t Eval(Double_t x, Double_t y=0, Double_t z=0);
virtual Double_t EvalPar(Double_t *x, Double_t *params=0);
virtual void ExecuteEvent(Int_t event, Int_t px, Int_t py);
Double_t GetChisquare() {return fChisquare;}
TH1F *GetHistogram();
Int_t GetNDF() {return fNpfits-fNpar;}
Int_t GetNpx() {return fNpx;}
TMethodCall *GetMethodCall() {return fMethodCall;}
Int_t GetNumberFitPoints() {return fNpfits;}
virtual Text_t *GetObjectInfo(Int_t px, Int_t py);
TObject *GetParent() {return fParent;}
Double_t GetParError(Int_t ipar) {return fParErrors[ipar];}
virtual void GetParLimits(Int_t ipar, Double_t &parmin, Double_t &parmax);
virtual Double_t GetProb() {return TMath::Prob(fChisquare,fNpfits-fNpar);}
virtual Double_t GetRandom();
virtual void GetRange(Float_t &xmin, Float_t &xmax);
virtual void GetRange(Float_t &xmin, Float_t &ymin, Float_t &xmax, Float_t &ymax);
virtual void GetRange(Float_t &xmin, Float_t &ymin, Float_t &zmin, Float_t &xmax, Float_t &ymax, Float_t &zmax);
virtual Double_t GetSave(Double_t *x);
Float_t GetXmin() {return fXmin;}
Float_t GetXmax() {return fXmax;}
virtual void InitArgs(Double_t *x, Double_t *params);
static void InitStandardFunctions();
virtual Double_t Integral(Float_t a, Float_t b, Double_t *params=0, Float_t epsilon=0.000001);
virtual Double_t Integral(Float_t ax, Float_t bx, Float_t ay, Float_t by, Float_t epsilon=0.000001);
virtual Double_t Integral(Float_t ax, Float_t bx, Float_t ay, Float_t by, Float_t az, Float_t bz, Float_t epsilon=0.000001);
virtual Double_t IntegralMultiple(Int_t n, Float_t *a, Float_t *b, Float_t epsilon, Float_t &relerr);
virtual void Paint(Option_t *option="");
virtual void Print(Option_t *option="");
virtual void Save(Float_t xmin, Float_t xmax);
virtual void SavePrimitive(ofstream &out, Option_t *option);
virtual void SetChisquare(Double_t chi2) {fChisquare = chi2;}
virtual void SetMaximum(Float_t maximum=-1111) {fMaximum=maximum;} // *MENU*
virtual void SetMinimum(Float_t minimum=-1111) {fMinimum=minimum;} // *MENU*
virtual void SetNumberFitPoints(Int_t npfits) {fNpfits = npfits;}
virtual void SetNpx(Int_t npx=100); // *MENU*
virtual void SetParError(Int_t ipar, Double_t error) {fParErrors[ipar] = error;}
virtual void SetParLimits(Int_t ipar, Double_t parmin, Double_t parmax);
virtual void SetParent(TObject *p=0) {fParent = p;}
virtual void SetRange(Float_t xmin, Float_t xmax); // *MENU*
virtual void SetRange(Float_t xmin, Float_t ymin, Float_t xmax, Float_t ymax);
virtual void SetRange(Float_t xmin, Float_t ymin, Float_t zmin, Float_t xmax, Float_t ymax, Float_t zmax);
virtual void Update();
ClassDef(TF1,3) //The Parametric 1-D function
};
inline void TF1::SetRange(Float_t xmin, Float_t, Float_t xmax, Float_t)
{ TF1::SetRange(xmin, xmax); }
inline void TF1::SetRange(Float_t xmin, Float_t, Float_t, Float_t xmax, Float_t, Float_t)
{ TF1::SetRange(xmin, xmax); }
#endif