66 GFunc(
const TF1*
function ,
double y ):fFunction(
function), fY0(y) {}
67 double operator()(
double x)
const {
68 return fFunction->
Eval(x) - fY0;
76 GInverseFunc(
const TF1*
function):fFunction(
function) {}
78 double operator()(
double x)
const {
79 return - fFunction->
Eval(x);
83 class GInverseFuncNdim {
86 GInverseFuncNdim(
TF1*
function):fFunction(
function) {}
88 double operator()(
const double*
x)
const {
105 fFunc->InitArgs(fX, fPar);
106 if (par) fFunc->SetParameters(par);
111 TF1_EvalWrapper * f =
new TF1_EvalWrapper( *
this);
112 f->fFunc->InitArgs(f->fX, f->fPar);
119 Double_t fval = fFunc->EvalPar( fX, 0);
120 if (fAbsVal && fval < 0)
return -fval;
126 return fX[0] *
TMath::Abs( fFunc->EvalPar( fX, 0) );
398 fXmin(0), fXmax(0), fNpar(0), fNdim(0),
400 fNpfits(0), fNDF(0), fChisquare(0),
401 fMinimum(-1111), fMaximum(-1111),
402 fParent(0), fHistogram(0),
403 fMethodCall(0), fNormalized(false), fNormIntegral(0),
404 fFormula(0), fParams(0)
449 if (
fNdim > 1 && xmin < xmax) {
450 Error(
"TF1",
"function: %s/%s has dimension %d instead of 1",name,formula,
fNdim);
487 Info(
"TF1",
"TF1 has name * - it is not well defined");
491 Error(
"TF1",
"requires a proper function name!");
499 Error(
"TF1",
"No function found with the signature %s(Double_t*,Double_t*)",name);
619 if (doAdd &&
gROOT) {
624 gROOT->GetListOfFunctions()->Remove(f1old);
625 gROOT->GetListOfFunctions()->Add(
this);
657 if (!
gROOT)
return false;
663 assert (
gROOT->GetListOfFunctions()->FindObject(
this) !=
nullptr);
670 gROOT->GetListOfFunctions()->Add(
this);
673 else if (prevStatus) {
679 Warning(
"AddToGlobalList",
"Function is supposed to be in the global list but it is not present");
682 gROOT->GetListOfFunctions()->Remove(
this);
711 if (
gROOT)
gROOT->GetListOfFunctions()->Remove(
this);
734 ((
TF1&)f1).Copy(*
this);
767 delete ((
TF1&)obj).fHistogram;
768 delete ((
TF1&)obj).fMethodCall;
792 ((
TF1&)obj).fHistogram = 0;
793 ((
TF1&)obj).fMethodCall = 0;
796 ((
TF1&)obj).fFormula = 0;
805 ((
TF1&)obj).fMethodCall =
m;
810 if (formulaToCopy)
delete formulaToCopy;
813 ((
TF1&)obj).fFormula = formulaToCopy;
817 if (paramsToCopy) *paramsToCopy = *
fParams;
861 Warning(
"Derivative",
"Function dimension is larger than one");
868 double h = eps* std::abs(xmax-xmin);
869 if ( h <= 0 ) h = 0.001;
927 Warning(
"Derivative2",
"Function dimension is larger than one");
934 double h = eps* std::abs(xmax-xmin);
935 if ( h <= 0 ) h = 0.001;
993 Warning(
"Derivative3",
"Function dimension is larger than one");
1000 double h = eps* std::abs(xmax-xmin);
1001 if ( h <= 0 ) h = 0.001;
1044 Int_t distance = 9999;
1047 if (distance <= 1)
return distance;
1054 xx[0] =
gPad->PadtoX(x);
1055 if (xx[0] <
fXmin || xx[0] >
fXmax)
return distance;
1111 TF1 *newf1 = (
TF1*)this->IsA()->New();
1139 if (padsav) padsav->
cd();
1164 if (padsav) padsav->
cd();
1197 if (
fType==2) ((
TF1*)
this)->InitArgs(xx,pp);
1198 return ((
TF1*)
this)->EvalPar(xx,pp);
1271 if (!
gPad->GetView()) {
1283 if (ipar < 0 || ipar >
GetNpar()-1)
return;
1295 ::Warning(
"TF1::GetCurrent",
"This function is obsolete and is working only for the current painted functions");
1309 const_cast<TF1*
>(
this)->
fHistogram = const_cast<TF1*>(
this)->CreateHistogram();
1335 if (xmin >= xmax) {xmin =
fXmin; xmax =
fXmax;}
1337 if (!logx &&
gPad != 0) logx =
gPad->GetLogx();
1340 GInverseFunc g(
this);
1345 bm.
Minimize(maxiter, epsilon, epsilon );
1373 if (xmin >= xmax) {xmin =
fXmin; xmax =
fXmax;}
1375 if (!logx &&
gPad != 0) logx =
gPad->GetLogx();
1378 GInverseFunc g(
this);
1383 bm.
Minimize(maxiter, epsilon, epsilon );
1411 if (xmin >= xmax) {xmin =
fXmin; xmax =
fXmax;}
1413 if (!logx &&
gPad != 0) logx =
gPad->GetLogx();
1420 bm.
Minimize(maxiter, epsilon, epsilon );
1439 Error(
"GetMinimumNDim",
"Function of dimension 0 - return Eval(x)");
1440 return (const_cast<TF1&>(*
this))(
x);
1449 Error(
"GetMinimumNDim",
"Error creating minimizer %s",minimName);
1460 GInverseFuncNdim invFunc(const_cast<TF1*>(
this));
1467 std::vector<double> rmin(ndim);
1468 std::vector<double> rmax(ndim);
1470 for (
int i = 0; i < ndim; ++i) {
1471 const char * xname = 0;
1472 double stepSize = 0.1;
1474 if (rmax[i] > rmin[i])
1475 stepSize = (rmax[i] - rmin[i])/100;
1476 else if (std::abs(x[i]) > 1.)
1477 stepSize = 0.1*x[i];
1496 if (rmin[i] < rmax[i] ) {
1507 Error(
"GetMinimumNDim",
"Error minimizing function %s",
GetName() );
1509 if (min->
X() ) std::copy (min->
X(), min->
X()+ndim,
x);
1513 return (findmax) ? -fmin : fmin;
1539 if (xmin >= xmax) {xmin =
fXmin; xmax =
fXmax;}
1546 bm.
Minimize(maxiter, epsilon, epsilon );
1576 if (xmin >= xmax) {xmin =
fXmin; xmax =
fXmax;}
1578 if (!logx &&
gPad != 0) logx =
gPad->GetLogx();
1586 brf.
Solve(maxiter, epsilon, epsilon);
1612 for (
Int_t i=0;i<nfree;i++) {
1613 ((
TF1*)
this)->GetParLimits(i,al,bl);
1614 if (al*bl != 0 && al >= bl) nfree--;
1627 static char info[64];
1639 if (ipar < 0 || ipar >
GetNpar()-1)
return 0;
1653 if (ipar < 0 || ipar > n-1)
return;
1664 if (
fNDF <= 0)
return 0;
1708 const Double_t dx = (xMax-xMin)/npx;
1716 Int_t intNegative = 0;
1718 for (i = 0; i < npx; i++) {
1720 if (integ < 0) {intNegative++; integ = -integ;}
1721 integral[i+1] = integral[i] + integ;
1724 if (intNegative > 0)
1725 Warning(
"GetQuantiles",
"function:%s has %d negative values: abs assumed",
1727 if (integral[npx] == 0) {
1728 Error(
"GetQuantiles",
"Integral of function is zero");
1733 for (i = 1; i <= npx; i++) integral[i] /= total;
1737 for (i = 0; i < npx; i++) {
1739 const Double_t r2 = integral[i+1]-integral[i];
1741 gamma[i] = (2*r2-4*
r1)/(dx*dx);
1742 beta[i] = r2/dx-gamma[i]*dx;
1749 for (i = 0; i < nprobSum; i++) {
1758 const Double_t rr = r-integral[bin];
1761 const Double_t fac = -2.*gamma[bin]*rr/beta[bin]/beta[bin];
1762 if (fac != 0 && fac <= 1)
1763 xx = (-beta[bin]+
TMath::Sqrt(beta[bin]*beta[bin]+2*gamma[bin]*rr))/gamma[bin];
1764 else if (beta[bin] != 0.)
1766 q[i] = alpha[bin]+xx;
1769 if (integral[bin+1] == r) q[i] += dx;
1808 Int_t intNegative = 0;
1814 if (xmin > 0 && xmax/xmin>
fNpx) {
1823 for (i=0;i<
fNpx;i++) {
1827 for (i=0;i<
fNpx;i++) {
1833 if (integ < 0) {intNegative++; integ = -integ;}
1836 if (intNegative > 0) {
1837 Warning(
"GetRandom",
"function:%s has %d negative values: abs assumed",
GetName(),intNegative);
1841 Error(
"GetRandom",
"Integral of function is zero");
1845 for (i=1;i<=
fNpx;i++) {
1852 for (i=0;i<
fNpx;i++) {
1921 Int_t intNegative = 0;
1923 for (i=0;i<
fNpx;i++) {
1925 if (integ < 0) {intNegative++; integ = -integ;}
1928 if (intNegative > 0) {
1929 Warning(
"GetRandom",
"function:%s has %d negative values: abs assumed",
GetName(),intNegative);
1932 Error(
"GetRandom",
"Integral of function is zero");
1936 for (i=1;i<=
fNpx;i++) {
1943 for (i=0;i<
fNpx;i++) {
1977 }
while(x<xmin || x>xmax);
1990 for (
int i = 0; i < ndim; ++i) {
1998 rmin[2] = zmin; rmax[2] = zmax;
2049 if (
fSave.size() == 0)
return 0;
2051 int fNsave =
fSave.size();
2057 xmin =
fSave[fNsave-3];
2058 xmax =
fSave[fNsave-2];
2059 if (
fSave[fNsave-1] == xmax) {
2068 ylow =
fSave[bin-bin1];
2069 yup =
fSave[bin-bin1+1];
2073 ylow =
fSave[bin-bin1-1];
2074 yup =
fSave[bin-bin1];
2077 y = ((xup*ylow-xlow*yup) + x*(yup-ylow))/dx;
2081 Int_t np = fNsave - 3;
2084 dx = (xmax-
xmin)/np;
2085 if (x < xmin || x > xmax)
return 0;
2088 if (dx <= 0)
return 0;
2091 xlow = xmin + bin*dx;
2095 y = ((xup*ylow-xlow*yup) + x*(yup-ylow))/dx;
2152 if(eps< 1e-10 || eps > 1) {
2153 Warning(
"Derivative",
"parameter esp=%g out of allowed range[1e-10,1], reset to 0.01",eps);
2168 ((
TF1*)
this)->GetParLimits(ipar,al,bl);
2169 if (al*bl != 0 && al >= bl) {
2182 parameters[ipar] = par0 +
h; f1 = func->
EvalPar(x,parameters);
2183 parameters[ipar] = par0 -
h; f2 = func->
EvalPar(x,parameters);
2184 parameters[ipar] = par0 + h/2; g1 = func->
EvalPar(x,parameters);
2185 parameters[ipar] = par0 - h/2; g2 = func->
EvalPar(x,parameters);
2195 parameters[ipar] =
par0;
2216 if(eps< 1e-10 || eps > 1) {
2217 Warning(
"Derivative",
"parameter esp=%g out of allowed range[1e-10,1], reset to 0.01",eps);
2234 if (params) args[1] = (
Long_t)params;
2248 if (!
gROOT->GetListOfFunctions()->FindObject(
"gaus")) {
2254 for (
Int_t i=0;i<10;i++) {
2255 f1 =
new TF1(
Form(
"pol%d",i),
Form(
"pol%d",i),-1,1);
2285 Warning(
"analytical integral not available for %s - with number %d compute numerical integral",
GetName(),
GetNumber());
2380 error = iod.
Error();
2393 error = iod.
Error();
2398 Warning(
"IntegralOneDim",
"Error found in integrating function %s in [%f,%f] using %s. Result = %f +/- %f - status = %d",
GetName(),a,b,igName.c_str(),
result,error,status);
2399 TString msg(
"\t\tFunction Parameters = {");
2400 for (
int ipar = 0; ipar <
GetNpar(); ++ipar) {
2405 Info(
"IntegralOneDim",
"%s",msg.Data());
2463 x1[0] =
a, x2[0] =
b;
2544 Warning(
"IntegralMultiple",
"failed code=%d, ",ifail);
2601 nfnevl = aimd.
NEval();
2608 relerr = (result != 0) ? imd.
Error()/ std::abs(result) : imd.
Error();
2638 printf(
"Formula based function: %s \n",
GetName());
2642 else if (
fType > 0) {
2644 printf(
"Interpreted based function: %s(double *x, double *p). Ndim = %d, Npar = %d \n",
GetName(),
GetNpar(),
GetNdim());
2647 printf(
"Compiled based function: %s based on a functor object. Ndim = %d, Npar = %d\n",
GetName(),
GetNpar(),
GetNdim());
2649 printf(
"Function based on a list of points from a compiled based function: %s. Ndim = %d, Npar = %d, Npx = %d\n",
GetName(),
GetNpar(),
GetNdim(),
int(
fSave.size()));
2651 Warning(
"Print",
"Function %s is based on a list of points but list is empty",
GetName());
2659 printf(
"List of Parameters: \n");
2660 for (
int i = 0; i <
fNpar; ++i)
2663 if (!
fSave.empty() ) {
2665 printf(
"List of Saved points (N=%d): \n",
int(
fSave.size()));
2667 printf(
"( %10f ) ",
x);
2673 printf(
"Contained histogram\n");
2695 pmin =
gPad->PadtoX(
gPad->GetUxmin());
2696 pmax =
gPad->PadtoX(
gPad->GetUxmax());
2699 if (xmax < pmin)
return;
2700 if (xmin > pmax)
return;
2701 if (xmin < pmin) xmin = pmin;
2702 if (xmax > pmax) xmax = pmax;
2711 if (minimum <= 0 && gPad && gPad->GetLogy()) minimum = -1111;
2714 if (minimum == -1111) {
2721 if (minimum == -1111) {
2723 if (optSAME &&
gPad) hmin =
gPad->GetUymin();
2728 if (optSAME &&
gPad) hmax =
gPad->GetUymax();
2730 hmin -= 0.05*(hmax-hmin);
2731 if (hmin < 0) hmin = 0;
2732 if (hmin <= 0 && gPad && gPad->GetLogy()) hmin = hminpos;
2739 if (maximum == -1111) {
2767 TH1 * histogram = 0;
2774 char *semicol = (
char*)strstr(
GetTitle(),
";");
2776 Int_t nxt = strlen(semicol);
2777 char *ctemp =
new char[nxt];
2778 strlcpy(ctemp,semicol+1,nxt);
2779 semicol = (
char*)strstr(ctemp,
";");
2801 if (xmin > 0 &&
gPad &&
gPad->GetLogx()) {
2806 for (i=0;i<=
fNpx;i++) {
2807 xbins[i] =
gPad->PadtoX(xlogmin+ i*dlogx);
2827 for (i=1;i<=
fNpx;i++) {
2857 if (ipar < 0 || ipar >
GetNpar()-1)
return;
2875 int fNsave = bin2-bin1+4;
2877 fSave.resize(fNsave);
2881 for (
Int_t i=bin1;i<=bin2;i++) {
2891 int fNsave =
fNpx+3;
2892 if (fNsave <= 3) {
return;}
2894 fSave.resize(fNsave);
2899 xmin =
fXmin +0.5*dx;
2900 xmax =
fXmax -0.5*dx;
2905 xv[0] = xmin + dx*i;
2922 if (strstr(option,
"cc")) {
2923 out <<
"double " <<
GetName() <<
"(double xv) {"<<std::endl;
2925 out <<
" double x[" <<
fNpx <<
"] = {" << std::endl;
2928 for (i=0; i<
fNpx; i++) {
2929 out <<
fXmin + dx*i ;
2930 if (i<fNpx-1) out <<
", ";
2931 if (n++ == 10) { out << std::endl; out <<
" "; n = 0;}
2934 out <<
" };" << std::endl;
2935 out <<
" double y[" << fNpx <<
"] = {" << std::endl;
2938 for (i=0; i<
fNpx; i++) {
2940 if (i<fNpx-1) out <<
", ";
2941 if (n++ == 10) { out << std::endl; out <<
" "; n = 0;}
2944 out <<
" };" << std::endl;
2945 out <<
" if (xv<x[0]) return y[0];" << std::endl;
2946 out <<
" if (xv>x[" << fNpx-1 <<
"]) return y[" << fNpx-1 <<
"];" << std::endl;
2947 out <<
" int i, j=0;" << std::endl;
2948 out <<
" for (i=1; i<" << fNpx <<
"; i++) { if (xv < x[i]) break; j++; }" << std::endl;
2949 out <<
" return y[j] + (y[j + 1] - y[j]) / (x[j + 1] - x[j]) * (xv - x[j]);" << std::endl;
2950 out <<
"}"<<std::endl;
2954 out<<
" "<<std::endl;
2957 static Int_t f1Number = 0;
2959 const char *
l = strstr(option,
"#");
2961 sscanf(&l[1],
"%d",&f1Number);
2968 out<<
" TF1 *"<<f1Name.
Data()<<
" = new TF1("<<quote<<
GetName()<<quote<<
","<<quote<<
GetTitle()<<quote<<
","<<
fXmin<<
","<<
fXmax<<
");"<<std::endl;
2970 out<<
" "<<f1Name.
Data()<<
"->SetNpx("<<
fNpx<<
");"<<std::endl;
2974 out<<
" //The original function : "<<
GetTitle()<<
" had originally been created by:" <<std::endl;
2976 out<<
" "<<f1Name.
Data()<<
"->SetRange("<<
fXmin<<
","<<
fXmax<<
");"<<std::endl;
2977 out<<
" "<<f1Name.
Data()<<
"->SetName("<<quote<<
GetName()<<quote<<
");"<<std::endl;
2978 out<<
" "<<f1Name.
Data()<<
"->SetTitle("<<quote<<
GetTitle()<<quote<<
");"<<std::endl;
2980 out<<
" "<<f1Name.
Data()<<
"->SetNpx("<<
fNpx<<
");"<<std::endl;
2986 for (i=0;i<=
fNpx;i++) {
2987 xv[0] =
fXmin + dx*i;
2989 out<<
" "<<f1Name.
Data()<<
"->SetSavedPoint("<<i<<
","<<save<<
");"<<std::endl;
2991 out<<
" "<<f1Name.
Data()<<
"->SetSavedPoint("<<
fNpx+1<<
","<<
fXmin<<
");"<<std::endl;
2992 out<<
" "<<f1Name.
Data()<<
"->SetSavedPoint("<<
fNpx+2<<
","<<
fXmax<<
");"<<std::endl;
2996 out<<
" "<<f1Name.
Data()<<
"->SetBit(TF1::kNotDraw);"<<std::endl;
3001 out<<
" "<<f1Name.
Data()<<
"->SetFillColor(ci);" << std::endl;
3003 out<<
" "<<f1Name.
Data()<<
"->SetFillColor("<<
GetFillColor()<<
");"<<std::endl;
3006 out<<
" "<<f1Name.
Data()<<
"->SetFillStyle("<<
GetFillStyle()<<
");"<<std::endl;
3011 out<<
" "<<f1Name.
Data()<<
"->SetMarkerColor(ci);" << std::endl;
3024 out<<
" "<<f1Name.
Data()<<
"->SetLineColor(ci);" << std::endl;
3026 out<<
" "<<f1Name.
Data()<<
"->SetLineColor("<<
GetLineColor()<<
");"<<std::endl;
3029 out<<
" "<<f1Name.
Data()<<
"->SetLineWidth("<<
GetLineWidth()<<
");"<<std::endl;
3032 out<<
" "<<f1Name.
Data()<<
"->SetLineStyle("<<
GetLineStyle()<<
");"<<std::endl;
3035 out<<
" "<<f1Name.
Data()<<
"->SetChisquare("<<
GetChisquare()<<
");"<<std::endl;
3036 out<<
" "<<f1Name.
Data()<<
"->SetNDF("<<
GetNDF()<<
");"<<std::endl;
3044 out<<
" "<<f1Name.
Data()<<
"->SetParameter("<<i<<
","<<
GetParameter(i)<<
");"<<std::endl;
3045 out<<
" "<<f1Name.
Data()<<
"->SetParError("<<i<<
","<<
GetParError(i)<<
");"<<std::endl;
3047 out<<
" "<<f1Name.
Data()<<
"->SetParLimits("<<i<<
","<<parmin<<
","<<parmax<<
");"<<std::endl;
3049 if (!strstr(option,
"nodraw")) {
3050 out<<
" "<<f1Name.
Data()<<
"->Draw(" 3051 <<quote<<option<<quote<<
");"<<std::endl;
3077 Warning(
"SetFitResult",
"Empty Fit result - nothing is set in TF1");
3080 if (indpar == 0 && npar != (
int) result.
NPar() ) {
3081 Error(
"SetFitResult",
"Invalid Fit result passed - number of parameter is %d , different than TF1::GetNpar() = %d",npar,result.
NPar());
3084 if (result.
Chi2() > 0)
3093 for (
Int_t i = 0; i < npar; ++i) {
3094 Int_t ipar = (indpar != 0) ? indpar[i] : i;
3095 if (ipar < 0)
continue;
3098 if (ipar < (
int) result.
Errors().size() )
3153 const Int_t minPx = 4;
3154 Int_t maxPx = 10000000;
3155 if (
GetNdim() > 1) maxPx = 10000;
3156 if (npx >= minPx && npx <= maxPx) {
3160 if(npx < minPx)
fNpx = minPx;
3161 if(npx > maxPx)
fNpx = maxPx;
3162 Warning(
"SetNpx",
"Number of points must be >=%d && <= %d, fNpx set to %d",minPx,maxPx,
fNpx);
3172 if (ipar <0 || ipar >=
GetNpar())
return;
3182 void TF1::SetParNames(
const char*name0,
const char*name1,
const char*name2,
const char*name3,
const char*name4,
3183 const char*name5,
const char*name6,
const char*name7,
const char*name8,
const char*name9,
const char*name10)
3186 fFormula->
SetParNames(name0,name1,name2,name3,name4,name5,name6,name7,name8,name9,name10);
3188 fParams->
SetParNames(name0,name1,name2,name3,name4,name5,name6,name7,name8,name9,name10);
3195 if (ipar < 0 || ipar >
GetNpar()-1)
return;
3206 if (!errors)
return;
3221 if (ipar < 0 || ipar > npar-1)
return;
3248 if (
fSave.size() == 0) {
3251 if (point < 0 || point >=
int(
fSave.size()))
return;
3252 fSave[point] = value;
3275 void TF1::Streamer(
TBuffer &b)
3289 gROOT->GetListOfFunctions()->Add(
this);
3300 if (fold.
fType == 0) {
3340 for (
int ibit = 0; ibit < 24; ++ibit)
3345 fOldLine.
Copy(*
this);
3347 fOldFill.
Copy(*
this);
3349 fOldMarker.
Copy(*
this);
3363 if (saved) {
fSave.clear(); }
3426 TF1_EvalWrapper
func(
this, params,
kTRUE, n);
3435 Error(
"Moment",
"Integral zero over range");
3460 TF1_EvalWrapper
func(
this, params,
kTRUE, n);
3469 Error(
"Moment",
"Integral zero over range");
3567 for (
unsigned int i = 0; i < fParNames.size(); ++i) {
3568 if (fParNames[i] == std::string(name) )
return i;
3580 unsigned int npar = fParameters.size();
3581 if (npar > 0) fParameters[0] = p0;
3582 if (npar > 1) fParameters[1] =
p1;
3583 if (npar > 2) fParameters[2] =
p2;
3584 if (npar > 3) fParameters[3] =
p3;
3585 if (npar > 4) fParameters[4] = p4;
3586 if (npar > 5) fParameters[5] = p5;
3587 if (npar > 6) fParameters[6] = p6;
3588 if (npar > 7) fParameters[7] = p7;
3589 if (npar > 8) fParameters[8] = p8;
3590 if (npar > 9) fParameters[9] = p9;
3591 if (npar >10) fParameters[10]= p10;
3598 const char *name4,
const char *name5,
const char *name6,
const char *name7,
3599 const char *name8,
const char *name9,
const char *name10)
3601 unsigned int npar = fParNames.size();
3602 if (npar > 0) fParNames[0] = name0;
3603 if (npar > 1) fParNames[1] = name1;
3604 if (npar > 2) fParNames[2] = name2;
3605 if (npar > 3) fParNames[3] = name3;
3606 if (npar > 4) fParNames[4] = name4;
3607 if (npar > 5) fParNames[5] = name5;
3608 if (npar > 6) fParNames[6] = name6;
3609 if (npar > 7) fParNames[7] = name7;
3610 if (npar > 8) fParNames[8] = name8;
3611 if (npar > 9) fParNames[9] = name9;
3612 if (npar >10) fParNames[10]= name10;
virtual const double * X() const =0
return pointer to X values at the minimum
ROOT::Math::ParamFunctor fFunctor
virtual const char * GetTitle() const
Returns title of object.
int NEval() const
return number of function evaluations in calculating the integral
double Error() const
Return the estimate of the absolute Error of the last Integral calculation.
virtual Int_t GetQuantiles(Int_t nprobSum, Double_t *q, const Double_t *probSum)
Compute Quantiles for density distribution of this function.
void Copy(TAttMarker &attmarker) const
Copy this marker attributes to a new TAttMarker.
virtual Style_t GetLineStyle() const
Return the line style.
virtual Style_t GetFillStyle() const
Return the fill area style.
virtual void SetLineWidth(Width_t lwidth)
Set the line width.
double Derivative3(double x)
Returns the third derivative of the function at point x, computed by Richardson's extrapolation metho...
virtual void SaveAttributes(std::ostream &out, const char *name, const char *subname)
Save axis attributes as C++ statement(s) on output stream out.
virtual void Paint(Option_t *option="")
Control routine to paint any kind of histograms.
virtual void SetParameters(const Double_t *params)
virtual Double_t GetRandom()
Return a random number following this function shape.
Bool_t fNormalized
Pointer to MethodCall in case of interpreted function.
double IntegralUp(double a)
Returns Integral of function on an upper semi-infinite interval.
virtual Int_t WriteClassBuffer(const TClass *cl, void *pointer)=0
Interface (abstract class) for generic functions objects of one-dimension Provides a method to evalua...
virtual void SetNpx(Int_t npx=100)
Set the number of points used to draw the function.
virtual void ReleaseParameter(Int_t ipar)
Release parameter number ipar If used in a fit, the parameter can vary freely.
virtual void SetMaximum(Double_t maximum=-1111)
bool Solve(int maxIter=100, double absTol=1E-8, double relTol=1E-10)
Returns the X value corresponding to the function value fy for (xmin<x<xmax).
static double p3(double t, double a, double b, double c, double d)
virtual void SetParName(Int_t ipar, const char *name)
Set name of parameter number ipar.
virtual void SetLimits(Double_t xmin, Double_t xmax)
This namespace contains pre-defined functions to be used in conjuction with TExecutor::Map and TExecu...
int Status() const
return status of integration
Option_t * GetDrawOption() const
Get option used by the graphics system to draw this object.
virtual Bool_t InheritsFrom(const char *classname) const
Returns kTRUE if object inherits from class "classname".
virtual void SetParNames(const char *name0="p0", const char *name1="p1", const char *name2="p2", const char *name3="p3", const char *name4="p4", const char *name5="p5", const char *name6="p6", const char *name7="p7", const char *name8="p8", const char *name9="p9", const char *name10="p10")
Set up to 10 parameter names.
const std::vector< double > & Errors() const
parameter errors (return st::vector)
virtual void SetDirectory(TDirectory *dir)
By default when an histogram is created, it is added to the list of histogram objects in the current ...
double RelError() const
return relative error
virtual Int_t GetNumberFreeParameters() const
Return the number of free parameters.
static std::atomic< Bool_t > fgAbsValue
void SetLogScan(bool on)
Set a log grid scan (default is equidistant bins) will work only if xlow > 0.
virtual const char * GetParName(Int_t ipar) const
R__EXTERN TStyle * gStyle
unsigned int Ndf() const
Number of degree of freedom.
void DoInitialize(EAddToList addToGlobList)
Common initialization of the TF1.
virtual Double_t IntegralFast(Int_t num, Double_t *x, Double_t *w, Double_t a, Double_t b, Double_t *params=0, Double_t epsilon=1e-12)
Gauss-Legendre integral, see CalcGaussLegendreSamplingPoints.
static TF1 * GetCurrent()
Static function returning the current function being processed.
static void SaveColor(std::ostream &out, Int_t ci)
Save a color with index > 228 as a C++ statement(s) on output stream out.
virtual void Info(const char *method, const char *msgfmt,...) const
Issue info message.
TObject * fParent
Array gamma.
virtual void SetRange(Double_t xmin, Double_t xmax)
Initialize the upper and lower bounds to draw the function.
void Copy(TAttLine &attline) const
Copy this line attributes to a new TAttLine.
virtual void ExecuteEvent(Int_t event, Int_t px, Int_t py)
Execute action corresponding to one event.
void ToUpper()
Change string to upper case.
Buffer base class used for serializing objects.
virtual void Save(Double_t xmin, Double_t xmax, Double_t ymin, Double_t ymax, Double_t zmin, Double_t zmax)
Save values of function in array fSave.
virtual Double_t IntegralMultiple(Int_t n, const Double_t *a, const Double_t *b, Int_t maxpts, Double_t epsrel, Double_t epsabs, Double_t &relerr, Int_t &nfnevl, Int_t &ifail)
This function computes, to an attempted specified accuracy, the value of the integral.
virtual void SetMinimum(Double_t minimum=-1111)
double Integral(const double *xmin, const double *xmax)
evaluate the integral with the previously given function between xmin[] and xmax[] ...
static ROOT::Math::Minimizer * CreateMinimizer(const std::string &minimizerType="", const std::string &algoType="")
static method to create the corrisponding Minimizer given the string Supported Minimizers types are: ...
Class to Wrap a ROOT Function class (like TF1) in a IParamFunction interface of one dimensions to be ...
virtual void SetNumberFitPoints(Int_t npfits)
static Bool_t DefaultAddToGlobalList(Bool_t on=kTRUE)
Static method to add/avoid to add automatically functions to the global list (gROOT->GetListOfFunctio...
static Double_t DerivativeError()
Static function returning the error of the last call to the of Derivative's functions.
void Copy(TAttFill &attfill) const
Copy this fill attributes to a new TAttFill.
Short_t Min(Short_t a, Short_t b)
void ToLower()
Change string to lower-case.
virtual Double_t GetParError(Int_t ipar) const
Return value of parameter number ipar.
R__EXTERN TVirtualMutex * gROOTMutex
virtual void Draw(Option_t *option="")
Draw this function with its current attributes.
virtual void SetFillStyle(Style_t fstyle)
Set the fill area style.
virtual Double_t Derivative2(Double_t x, Double_t *params=0, Double_t epsilon=0.001) const
Returns the second derivative of the function at point x, computed by Richardson's extrapolation meth...
double Integral(double a, double b)
Returns Integral of function between a and b.
Width_t GetFuncWidth() const
virtual void SetMinimum(Double_t minimum=-1111)
Set the minimum value along Y for this function In case the function is already drawn, set also the minimum in the helper histogram.
double Derivative2(double x)
Returns the second derivative of the function at point x, computed by Richardson's extrapolation meth...
const Double_t * GetParameters() const
virtual Double_t Integral(Double_t a, Double_t b, Double_t epsrel=1.e-12)
IntegralOneDim or analytical integral.
User class for performing function integration.
virtual TObject * DrawIntegral(Option_t *option="al")
Draw integral of this function.
virtual Int_t GetNDF() const
Return the number of degrees of freedom in the fit the fNDF parameter has been previously computed du...
virtual Double_t GetXmin() const
virtual void Draw(Option_t *chopt="")
Draw this graph with its current attributes.
virtual TVirtualPad * cd(Int_t subpadnumber=0)=0
bit set when zooming on Y axis
Double_t Prob(Double_t chi2, Int_t ndf)
Computation of the probability for a certain Chi-squared (chi2) and number of degrees of freedom (ndf...
virtual TH1 * GetHistogram() const
Return a pointer to the histogram used to visualise the function.
LongDouble_t Power(LongDouble_t x, LongDouble_t y)
virtual bool SetLimitedVariable(unsigned int ivar, const std::string &name, double val, double step, double lower, double upper)
set a new upper/lower limited variable (override if minimizer supports them ) otherwise as default se...
void SetBit(UInt_t f, Bool_t set)
Set or unset the user status bits as specified in f.
double IntegralLow(double b)
Returns Integral of function on a lower semi-infinite interval.
double beta(double x, double y)
Calculates the beta function.
static void SetCurrent(TF1 *f1)
Static function setting the current function.
virtual Double_t GetMinMaxNDim(Double_t *x, Bool_t findmax, Double_t epsilon=0, Int_t maxiter=0) const
Find the minimum of a function of whatever dimension.
if object in a list can be deleted
virtual void AppendPad(Option_t *option="")
Append graphics object to current pad.
virtual void RecursiveRemove(TObject *obj)
Recursively remove this object from a list.
double Error() const
return the estimate of the absolute Error of the last Integral calculation
Template class to wrap any C++ callable object which takes one argument i.e.
void SetParamPtrs(void *paramArr, Int_t nparam=-1)
ParamArr is an array containing the function argument values.
const char * Data() const
virtual Int_t DistancetoPrimitive(Int_t px, Int_t py)
Compute distance from point px,py to a function.
virtual Double_t Moment(Double_t n, Double_t a, Double_t b, const Double_t *params=0, Double_t epsilon=0.000001)
Return nth moment of function between a and b.
virtual void SetRelTolerance(double eps)
Set the desired relative Error.
void SetParameters(const double *p)
set parameter values
TFormula * fFormula
Functor object to wrap any C++ callable object.
static const double x2[5]
Fill Area Attributes class.
virtual Bool_t IsValid() const
Return kTRUE if the function is valid.
static TString Format(const char *fmt,...)
Static method which formats a string using a printf style format descriptor and return a TString...
Double_t GetChisquare() const
int Status() const
return the Error Status of the last Integral calculation
The TNamed class is the base class for all named ROOT classes.
virtual Double_t GetMaximumStored() const
unsigned int r3[N_CITIES]
virtual void DrawF1(Double_t xmin, Double_t xmax, Option_t *option="")
Draw function between xmin and xmax.
void GetWeightVectors(double *x, double *w) const
Returns the arrays x and w containing the abscissa and weight of the Gauss-Legendre n-point quadratur...
virtual void Copy(TObject &named) const
Copy this to obj.
Color_t GetFuncColor() const
Abstract Minimizer class, defining the interface for the various minimizer (like Minuit2, Minuit, GSL, etc..) Plug-in's exist in ROOT to be able to instantiate the derived classes like ROOT::Math::GSLMinimizer or ROOT::Math::Minuit2Minimizer via the plug-in manager.
Double_t Log10(Double_t x)
static double p2(double t, double a, double b, double c)
static double DefaultRelTolerance()
virtual void Copy(TObject &f1) const
Copy this F1 to a new F1.
TF1()
TF1 default constructor.
virtual bool Minimize()=0
method to perform the minimization
void SetParameters(const Double_t *params)
virtual void SetMarkerColor(Color_t mcolor=1)
Set the marker color.
double Derivative1(double x)
Returns the first derivative of the function at point x, computed by Richardson's extrapolation metho...
double Error() const
return integration error
static IntegrationOneDim::Type DefaultIntegratorType()
you should not use this method at all Int_t Int_t Double_t bm
virtual void GetRange(Double_t *xmin, Double_t *xmax) const
Return range of a generic N-D function.
Method or function calling interface.
virtual void Error(const char *method, const char *msgfmt,...) const
Issue error message.
User class for performing function minimization.
TVirtualPad is an abstract base class for the Pad and Canvas classes.
int Status() const
return the Error Status of the last Integral calculation
virtual Double_t GetXmax() const
virtual double MinValue() const =0
return minimum function value
std::vector< Double_t > fIntegral
double Integral(const double *xmin, const double *xmax)
evaluate the integral with the previously given function between xmin[] and xmax[] ...
virtual void GetParLimits(Int_t ipar, Double_t &parmin, Double_t &parmax) const
Return limits for parameter ipar.
virtual Double_t Rndm()
Machine independent random number generator.
virtual Double_t GradientPar(Int_t ipar, const Double_t *x, Double_t eps=0.01)
Compute the gradient (derivative) wrt a parameter ipar.
double IntegralLow(const IGenFunction &f, double b)
evaluate the Integral of a function f over the over the semi-infinite interval (-inf,b)
double Root() const
Returns root value.
Bool_t AreEqualRel(Double_t af, Double_t bf, Double_t relPrec)
virtual Int_t DistancetoPrimitive(Int_t px, Int_t py)
Compute distance from point px,py to a line.
unsigned int NFreeParameters() const
get total number of free parameters
virtual void SetLineColor(Color_t lcolor)
Set the line color.
double Integral(Function &f, double a, double b)
evaluate the Integral of a function f over the defined interval (a,b)
Using a TBrowser one can browse all ROOT objects.
virtual void ExecuteEvent(Int_t event, Int_t px, Int_t py)
Execute action corresponding to one event.
virtual void SetChisquare(Double_t chi2)
virtual Double_t GetBinCenter(Int_t bin) const
Return bin center for 1D histogram.
virtual Size_t GetMarkerSize() const
Return the marker size.
virtual void SetParLimits(Int_t ipar, Double_t parmin, Double_t parmax)
Set limits for parameter ipar.
double IntegralUp(const IGenFunction &f, double a)
evaluate the Integral of a function f over the semi-infinite interval (a,+inf)
static Bool_t fgRejectPoint
virtual void SetFunction(const ROOT::Math::IMultiGenFunction &func)=0
set the function to minimize
unsigned int NPar() const
total number of parameters (abbreviation)
virtual void Print(Option_t *option="") const
Print some global quantities for this histogram.
TH1 * fHistogram
Parent object hooking this function (if one)
Class to manage histogram axis.
static const std::string & DefaultMinimizerType()
virtual Double_t GetMinimumStored() const
virtual void SetFillColor(Color_t fcolor)
Set the fill area color.
void SetLogScan(bool on)
Set a log grid scan (default is equidistant bins) will work only if xlow > 0.
double IntegralError(TF1 *func, Int_t ndim, const double *a, const double *b, const double *params, const double *covmat, double epsilon)
virtual void SetParError(Int_t ipar, Double_t error)
Set error for parameter number ipar.
unsigned int r1[N_CITIES]
void SetFunction(const ROOT::Math::IGenFunction &f, double xlow, double xup)
Sets function to be minimized.
virtual void SetBinContent(Int_t bin, Double_t content)
Set bin content see convention for numbering bins in TH1::GetBin In case the bin number is greater th...
virtual Color_t GetFillColor() const
Return the fill area color.
virtual double XMinimum() const
Return current estimate of the position of the minimum.
Bool_t TestBit(UInt_t f) const
virtual Int_t GetNdim() const
char * Form(const char *fmt,...)
User class for performing function integration.
virtual void FixParameter(Int_t ipar, Double_t value)
Fix the value of a parameter The specified value will be used in a fit operation. ...
User Class for performing numerical integration of a function in one dimension.
static void RejectPoint(Bool_t reject=kTRUE)
Static function to set the global flag to reject points the fgRejectPoint global flag is tested by al...
virtual const char * GetName() const
Returns name of object.
static void InitStandardFunctions()
Create the basic function objects.
virtual void SetMarkerStyle(Style_t mstyle=1)
Set the marker style.
const char * GetTitle() const
Returns title of object.
static double p1(double t, double a, double b)
virtual void Update()
Called by functions such as SetRange, SetNpx, SetParameters to force the deletion of the associated h...
virtual ~TF1()
TF1 default destructor.
virtual Double_t GetProb() const
Return the fit probability.
R__EXTERN TRandom * gRandom
virtual Color_t GetLineColor() const
Return the line color.
virtual TH1 * DoCreateHistogram(Double_t xmin, Double_t xmax, Bool_t recreate=kFALSE)
Create histogram with bin content equal to function value computed at the bin center This histogram w...
virtual void SetMarkerSize(Size_t msize=1)
Set the marker size.
static IntegrationMultiDim::Type DefaultIntegratorType()
TAxis * GetYaxis() const
Get y axis of the function.
const Double_t * GetArray() const
void InitWithPrototype(TClass *cl, const char *method, const char *proto, Bool_t objectIsConst=kFALSE, ROOT::EFunctionMatchMode mode=ROOT::kConversionMatch)
Initialize the method invocation environment.
#define R__LOCKGUARD2(mutex)
static void CalcGaussLegendreSamplingPoints(Int_t num, Double_t *x, Double_t *w, Double_t eps=3.0e-11)
Type safe interface (static method) The number of sampling points are taken from the TGraph...
Class for finding the root of a one dimensional function using the Brent algorithm.
virtual Double_t GetMinimum(Double_t xmin=0, Double_t xmax=0, Double_t epsilon=1.E-10, Int_t maxiter=100, Bool_t logx=false) const
Returns the minimum value of the function on the (xmin, xmax) interval.
virtual void SetMaximum(Double_t maximum=-1111)
Set the maximum value along Y for this function In case the function is already drawn, set also the maximum in the helper histogram.
virtual Int_t FindBin(Double_t x)
Find bin number corresponding to abscissa x.
virtual void SetTitle(const char *title="")
Set function title if title has the form "fffffff;xxxx;yyyy", it is assumed that the function title i...
static unsigned int total
virtual Int_t ReadClassBuffer(const TClass *cl, void *pointer, const TClass *onfile_class=0)=0
virtual Double_t Derivative(Double_t x, Double_t *params=0, Double_t epsilon=0.001) const
Returns the first derivative of the function at point x, computed by Richardson's extrapolation metho...
static const std::string & DefaultMinimizerAlgo()
void SetMaxFunctionCalls(unsigned int maxfcn)
set maximum of function calls
class containg the result of the fit and all the related information (fitted parameter values...
static const double x1[5]
void SetFunction(const IGenFunction &)
Set integration function (flag control if function must be copied inside).
std::vector< Double_t > fSave
virtual Double_t GetX(Double_t y, Double_t xmin=0, Double_t xmax=0, Double_t epsilon=1.E-10, Int_t maxiter=100, Bool_t logx=false) const
Returns the X value corresponding to the function value fy for (xmin<x<xmax).
void SetTolerance(double tol)
set the tolerance
std::vector< Double_t > fParErrors
virtual Double_t GetMaximum(Double_t xmin=0, Double_t xmax=0, Double_t epsilon=1.E-10, Int_t maxiter=100, Bool_t logx=false) const
Returns the maximum value of the function.
virtual Double_t Derivative3(Double_t x, Double_t *params=0, Double_t epsilon=0.001) const
Returns the third derivative of the function at point x, computed by Richardson's extrapolation metho...
virtual void Print(Option_t *option="") const
Print TNamed name and title.
void SetParName(Int_t iparam, const char *name)
double Error(unsigned int i) const
parameter error by index
Double_t AnalyticalIntegral(TF1 *f, Double_t a, Double_t b)
double func(double *x, double *p)
std::vector< Double_t > fGamma
Array beta. is approximated by x = alpha +beta*r *gamma*r**2.
void Streamer(TBuffer &b, Int_t version, UInt_t start, UInt_t count, const TClass *onfile_class=0)
Stream a class object.
you should not use this method at all Int_t Int_t Double_t Double_t Double_t e
static std::string GetName(IntegrationOneDim::Type)
static function to get a string from the enumeration
virtual void SavePrimitive(std::ostream &out, Option_t *option="")
Save primitive as a C++ statement(s) on output stream out.
virtual Double_t Uniform(Double_t x1=1)
Returns a uniform deviate on the interval (0, x1).
std::vector< Double_t > fAlpha
Integral of function binned on fNpx bins.
Style_t GetFuncStyle() const
virtual Int_t GetNumber() const
virtual void SetLineStyle(Style_t lstyle)
Set the line style.
Template class to wrap any C++ callable object implementing operator() (const double * x) in a multi-...
Int_t GetParNumber(const char *name) const
Returns the parameter number given a name not very efficient but list of parameters is typically smal...
Array of doubles (64 bits per element).
Namespace for new Math classes and functions.
virtual void InitArgs(const Double_t *x, const Double_t *params)
Initialize parameters addresses.
std::vector< Double_t > fBeta
Array alpha. for each bin in x the deconvolution r of fIntegral.
TAxis * GetZaxis() const
Get z axis of the function. (In case this object is a TF2 or TF3)
void SetParNames(const char *name0="p0", const char *name1="p1", const char *name2="p2", const char *name3="p3", const char *name4="p4", const char *name5="p5", const char *name6="p6", const char *name7="p7", const char *name8="p8", const char *name9="p9", const char *name10="p10")
Set parameter names.
virtual Double_t GetSave(const Double_t *x)
Get value corresponding to X in array of fSave values.
virtual Double_t GetParameter(Int_t ipar) const
TAxis * GetXaxis() const
Get x axis of the function.
virtual Double_t GetBinCenter(Int_t bin) const
Return center of bin.
Mother of all ROOT objects.
virtual Double_t * GetParameters() const
you should not use this method at all Int_t Int_t z
double MinFcnValue() const
Return value of the objective function (chi2 or likelihood) used in the fit.
double Error() const
Returns the estimate of the absolute Error of the last derivative calculation.
TMethodCall * fMethodCall
Pointer to histogram used for visualisation.
virtual Double_t IntegralOneDim(Double_t a, Double_t b, Double_t epsrel, Double_t epsabs, Double_t &err)
Return Integral of function between a and b using the given parameter values and relative and absolut...
virtual Double_t GetMaximumX(Double_t xmin=0, Double_t xmax=0, Double_t epsilon=1.E-10, Int_t maxiter=100, Bool_t logx=false) const
Returns the X value corresponding to the maximum value of the function.
virtual Color_t GetMarkerColor() const
Return the marker color.
virtual bool SetVariable(unsigned int ivar, const std::string &name, double val, double step)=0
set a new free variable
virtual double FValMinimum() const
Return function value at current estimate of the minimum.
User class for performing multidimensional integration.
static Bool_t RejectedPoint()
See TF1::RejectPoint above.
double f2(const double *x)
virtual void SetParErrors(const Double_t *errors)
Set errors for all active parameters when calling this function, the array errors must have at least ...
void Execute(const char *, const char *, int *=0)
Execute method on this object with the given parameter string, e.g.
void SetNpx(int npx)
Set the number of point used to bracket root using a grid.
virtual Bool_t AddToGlobalList(Bool_t on=kTRUE)
Add to global list of functions (gROOT->GetListOfFunctions() ) return previous status (true if the fu...
Short_t Max(Short_t a, Short_t b)
Bool_t Contains(const char *pat, ECaseCompare cmp=kExact) const
virtual void SetSavedPoint(Int_t point, Double_t value)
Restore value of function saved at point.
A Graph is a graphics object made of two arrays X and Y with npoints each.
virtual void SetFitResult(const ROOT::Fit::FitResult &result, const Int_t *indpar=0)
Set the result from the fit parameter values, errors, chi2, etc...
virtual Double_t Eval(Double_t x, Double_t y=0, Double_t z=0, Double_t t=0) const
Evaluate this function.
Param Functor class for Multidimensional functions.
you should not use this method at all Int_t Int_t Double_t Double_t Double_t Int_t Double_t Double_t Double_t Double_t b
static std::atomic< Bool_t > fgAddToGlobList
THist< 1, double, THistStatContent, THistStatUncertainty > TH1D
static void AbsValue(Bool_t reject=kTRUE)
Static function: set the fgAbsValue flag.
virtual void SetNDF(Int_t ndf)
Set the number of degrees of freedom ndf should be the number of points used in a fit - the number of...
virtual Double_t GetMinimumX(Double_t xmin=0, Double_t xmax=0, Double_t epsilon=1.E-10, Int_t maxiter=100, Bool_t logx=false) const
Returns the X value corresponding to the minimum value of the function on the (xmin, xmax) interval.
double Chi2() const
Chi2 fit value in case of likelihood must be computed ?
virtual Double_t IntegralError(Double_t a, Double_t b, const Double_t *params=0, const Double_t *covmat=0, Double_t epsilon=1.E-2)
Return Error on Integral of a parametric function between a and b due to the parameter uncertainties...
int Status() const
return the status of the last integration - 0 in case of success
TF1 & operator=(const TF1 &rhs)
Operator =.
virtual void SetParameter(Int_t param, Double_t value)
virtual void SetTitle(const char *title)
Change (i.e.
virtual bool Minimize(int maxIter, double absTol=1.E-8, double relTol=1.E-10)
Find minimum position iterating until convergence specified by the absolute and relative tolerance or...
bool SetFunction(const ROOT::Math::IGenFunction &f, double xlow, double xup)
Sets the function for the rest of the algorithms.
void SetNpx(int npx)
Set the number of point used to bracket root using a grid.
Functor1D class for one-dimensional functions.
Double_t Sqrt(Double_t x)
std::vector< Double_t > fParMax
double Parameter(unsigned int i) const
parameter value by index
virtual Style_t GetMarkerStyle() const
Return the marker style.
virtual char * GetObjectInfo(Int_t px, Int_t py) const
Redefines TObject::GetObjectInfo.
virtual Double_t CentralMoment(Double_t n, Double_t a, Double_t b, const Double_t *params=0, Double_t epsilon=0.000001)
Return nth central moment of function between a and b (i.e the n-th moment around the mean value) ...
virtual Double_t GetMaximum(Double_t maxval=FLT_MAX) const
Return maximum value smaller than maxval of bins in the range, unless the value has been overridden b...
bool IsEmpty() const
True if a fit result does not exist (even invalid) with parameter values.
virtual Width_t GetLineWidth() const
Return the line width.
virtual Double_t EvalPar(const Double_t *x, const Double_t *params=0)
Evaluate function with given coordinates and parameters.
virtual void SetTitle(const char *title="")
Set the title of the TNamed.
virtual void Browse(TBrowser *b)
Browse.
double norm(double *x, double *p)
std::vector< Double_t > fParMin
virtual Int_t GetNpar() const
Long64_t BinarySearch(Long64_t n, const T *array, T value)
Bool_t IsValid() const
Return true if the method call has been properly initialized and is usable.
static Double_t gErrorTF1
virtual Double_t GetMinimum(Double_t minval=-FLT_MAX) const
Return minimum value larger than minval of bins in the range, unless the value has been overridden by...
User class for calculating the derivatives of a function.
unsigned int r2[N_CITIES]
class for adaptive quadrature integration in multi-dimensions using rectangular regions.
virtual Version_t ReadVersion(UInt_t *start=0, UInt_t *bcnt=0, const TClass *cl=0)=0
virtual void Paint(Option_t *option="")
Paint this function with its current attributes.
virtual TObject * DrawDerivative(Option_t *option="al")
Draw derivative of this function.
virtual TF1 * DrawCopy(Option_t *option="") const
Draw a copy of this function with its current attributes.
virtual void Warning(const char *method, const char *msgfmt,...) const
Issue warning message.