50 :
TF2(name,formula,xmin,xmax,ymax,ymin)
56 if (ndim != 3 && xmin < xmax) {
57 Error(
"TF3",
"function: %s/%s has %d parameters instead of 3",name,formula,ndim);
75 TF3::TF3(
const char *
name,
Double_t (*fcn)(
Double_t *,
Double_t *),
Double_t xmin,
Double_t xmax,
Double_t ymin,
Double_t ymax,
Double_t zmin,
Double_t zmax,
Int_t npar,
Int_t ndim)
76 :
TF2(name,fcn,xmin,xmax,ymin,ymax,npar,ndim)
96 TF3::TF3(
const char *name,
Double_t (*fcn)(
const Double_t *,
const Double_t *),
Double_t xmin,
Double_t xmax,
Double_t ymin,
Double_t ymax,
Double_t zmin,
Double_t zmax,
Int_t npar,
Int_t ndim)
97 :
TF2(name,fcn,xmin,xmax,ymin,ymax,npar,ndim),
114 :
TF2(name, f, xmin, xmax, ymin, ymax, npar, ndim),
215 Double_t rsign = (findmax) ? -1. : 1.;
216 TF3 &
function =
const_cast<TF3&
>(*this);
217 Double_t xxmin = 0, yymin = 0, zzmin = 0, ttmin = 0;
227 xx[0]=
fXmin + (i+0.5)*dx;
229 xx[1]=
fYmin+(j+0.5)*dy;
231 xx[2] =
fZmin+(k+0.5)*dz;
233 if (rsign*tt < rsign*ttmin) {xxmin = xx[0], yymin = xx[1]; zzmin = xx[2]; ttmin=
tt;}
246 zzmin =
function(xx);
253 if (rsign*fmin < rsign*zzmin) {
254 if (x) {x[0] = xx[0]; x[1] = xx[1]; x[2] = xx[2];}
258 if (x) { x[0] = xxmin; x[1] = yymin; x[2] = zzmin; }
280 double xx[3] = { 0,0,0 };
283 x = xx[0]; y = xx[1]; z = xx[2];
296 double xx[3] = { 0,0,0 };
299 x = xx[0]; y = xx[1]; z = xx[2];
338 Int_t intNegative = 0;
340 for (k=0;k<
fNpz;k++) {
341 xx[2] =
fZmin+(k+0.5)*dz;
342 for (j=0;j<
fNpy;j++) {
343 xx[1] =
fYmin+(j+0.5)*dy;
344 for (i=0;i<
fNpx;i++) {
345 xx[0] =
fXmin+(i+0.5)*dx;
346 integ =
EvalPar(xx,parameters);
347 if (integ < 0) {intNegative++; integ = -integ;}
353 if (intNegative > 0) {
354 Warning(
"GetRandom3",
"function:%s has %d negative values: abs assumed",
GetName(),intNegative);
357 Error(
"GetRandom3",
"Integral of function is zero");
360 for (i=1;i<=ncells;i++) {
397 if (
fSave.empty())
return 0;
410 if (x < xmin || x > xmax)
return 0;
411 if (dx <= 0)
return 0;
414 if (y < ymin || y > ymax)
return 0;
415 if (dy <= 0)
return 0;
418 if (z < zmin || z > zmax)
return 0;
419 if (dz <= 0)
return 0;
431 Int_t k1 = (ibin ) + (npx+1)*((jbin ) + (npy+1)*(kbin ));
432 Int_t k2 = (ibin+1) + (npx+1)*((jbin ) + (npy+1)*(kbin ));
433 Int_t k3 = (ibin+1) + (npx+1)*((jbin+1) + (npy+1)*(kbin ));
434 Int_t k4 = (ibin ) + (npx+1)*((jbin+1) + (npy+1)*(kbin ));
435 Int_t k5 = (ibin ) + (npx+1)*((jbin ) + (npy+1)*(kbin+1));
436 Int_t k6 = (ibin+1) + (npx+1)*((jbin ) + (npy+1)*(kbin+1));
437 Int_t k7 = (ibin+1) + (npx+1)*((jbin+1) + (npy+1)*(kbin+1));
438 Int_t k8 = (ibin ) + (npx+1)*((jbin+1) + (npy+1)*(kbin+1));
463 Warning(
"Integral",
"failed code=%d, maxpts=%d, epsrel=%g, nfnevl=%d, relerr=%g ",ifail,maxpts,epsrel,nfnevl,relerr);
538 Int_t fNsave = nsave+9;
541 fSave.resize(fNsave);
548 xmin =
fXmin +0.5*dx;
549 xmax =
fXmax -0.5*dx;
553 ymin =
fYmin +0.5*dy;
554 ymax =
fYmax -0.5*dy;
558 zmin =
fZmin +0.5*dz;
559 zmax =
fZmax -0.5*dz;
564 for (k=0;k<=
fNpz;k++) {
566 for (j=0;j<=
fNpy;j++) {
568 for (i=0;i<=
fNpx;i++) {
579 fSave[nsave+4] = zmin;
580 fSave[nsave+5] = zmax;
607 out<<
" "<<
GetName()<<
"->SetFillColor(ci);" << std::endl;
614 out<<
" "<<
GetName()<<
"->SetLineColor(ci);" << std::endl;
629 out<<
" "<<
GetName()<<
"->SetParLimits("<<i<<
","<<parmin<<
","<<parmax<<
");"<<std::endl;
632 <<quote<<option<<quote<<
");"<<std::endl;
667 Warning(
"SetNpz",
"Number of points must be >=4 && <= 10000, fNpz set to 4");
669 }
else if(npz > 10000) {
670 Warning(
"SetNpz",
"Number of points must be >=4 && <= 10000, fNpz set to 10000");
695 void TF3::Streamer(
TBuffer &R__b)
711 if (saved) {
fSave.clear(); }
723 Error(
"Moment3",
"Integral zero over range");
727 TF3 fnc(
"TF3_ExpValHelper",
Form(
"%s*pow(x,%f)*pow(y,%f)*pow(z,%f)",
GetName(),nx,ny,nz));
739 Error(
"CentralMoment3",
"Integral zero over range");
758 TF3 fnc(
"TF3_ExpValHelper",
Form(
"%s*pow(x-%f,%f)*pow(y-%f,%f)*pow(z-%f,%f)",
GetName(),xbar,nx,ybar,ny,zbar,nz));
virtual const char * GetTitle() const
Returns title of object.
virtual void Paint(Option_t *option="")
Control routine to paint any kind of histograms.
virtual Int_t WriteClassBuffer(const TClass *cl, void *pointer)=0
virtual Double_t Moment3(Double_t nx, Double_t ax, Double_t bx, Double_t ny, Double_t ay, Double_t by, Double_t nz, Double_t az, Double_t bz, Double_t epsilon=0.000001)
Return x^nx * y^ny * z^nz moment of a 3d function in range [ax,bx],[ay,by],[az,bz].
virtual Bool_t IsInside(const Double_t *x) const
Return kTRUE is the point is inside the function range.
virtual void Draw(Option_t *option="")
Draw this function with its current attributes.
TVirtualHistPainter * GetPainter(Option_t *option="")
Return pointer to painter.
virtual void SetDirectory(TDirectory *dir)
By default when an histogram is created, it is added to the list of histogram objects in the current ...
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 TH1 * CreateHistogram()
Create a histogram for axis range.
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 ~TF3()
F3 default destructor.
virtual void ExecuteEvent(Int_t event, Int_t px, Int_t py)
Execute action corresponding to one event.
Buffer base class used for serializing objects.
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.
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.
virtual void GetRange(Double_t &xmin, Double_t &xmax) const
Return range of a 1-D function.
virtual void SavePrimitive(std::ostream &out, Option_t *option="")
Save primitive as a C++ statement(s) on output stream out.
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.
virtual Double_t GetMinimumXYZ(Double_t &x, Double_t &y, Double_t &z)
Compute the X, Y and Z values corresponding to the minimum value of the function on its range...
virtual void AppendPad(Option_t *option="")
Append graphics object to current pad.
virtual void SetClippingBoxOn(Double_t xclip=0, Double_t yclip=0, Double_t zclip=0)
Set the function clipping box (for drawing) "on" and define the clipping box.
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 FindMinMax(Double_t *x, bool findmax) const
Return minimum/maximum value of the function.
virtual void ExecuteEvent(Int_t event, Int_t px, Int_t py)
Execute action corresponding to one event.
Double_t GetChisquare() const
virtual Double_t GetMaximumXYZ(Double_t &x, Double_t &y, Double_t &z)
Compute the X, Y and Z values corresponding to the maximum value of the function on its range...
TF3()
F3 default constructor.
virtual void Paint(Option_t *option="")
Paint this 3-D function with its current attributes.
virtual void ProcessMessage(const char *mess, const TObject *obj)=0
virtual void Error(const char *method, const char *msgfmt,...) const
Issue error message.
std::vector< Double_t > fIntegral
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.
TH1 * fHistogram
Parent object hooking this function (if one)
virtual Double_t Integral(Double_t ax, Double_t bx, Double_t ay, Double_t by, Double_t az, Double_t bz, Double_t epsrel=1.e-6)
Return Integral of a 3d function in range [ax,bx],[ay,by],[az,bz] with a desired relative accuracy...
virtual Double_t CentralMoment3(Double_t nx, Double_t ax, Double_t bx, Double_t ny, Double_t ay, Double_t by, Double_t nz, Double_t az, Double_t bz, Double_t epsilon=0.000001)
Return x^nx * y^ny * z^nz central moment of a 3d function in range [ax,bx],[ay,by],[az,bz].
virtual void Copy(TObject &f3) const
Copy this F3 to a new F3.
A 3-Dim function with parameters.
virtual Color_t GetFillColor() const
Return the fill area color.
virtual Int_t GetNdim() const
char * Form(const char *fmt,...)
virtual void GetRandom3(Double_t &xrandom, Double_t &yrandom, Double_t &zrandom)
Return 3 random numbers following this function shape.
virtual void Copy(TObject &f2) const
Copy this F2 to a new F2.
virtual const char * GetName() const
Returns name of object.
virtual void Update()
Called by functions such as SetRange, SetNpx, SetParameters to force the deletion of the associated h...
A 2-Dim function with parameters.
R__EXTERN TRandom * gRandom
virtual Color_t GetLineColor() const
Return the line color.
TF3 & operator=(const TF3 &rhs)
Operator =.
virtual Int_t ReadClassBuffer(const TClass *cl, void *pointer, const TClass *onfile_class=0)=0
std::vector< Double_t > fSave
virtual void SetClippingBoxOff()
Set the function clipping box (for drawing) "off".
virtual void InitArgs(const Double_t *x, const Double_t *params)
Initialize parameters addresses.
THist< 3, float, THistStatContent, THistStatUncertainty > TH3F
virtual Double_t GetParameter(Int_t ipar) const
Mother of all ROOT objects.
virtual Double_t * GetParameters() const
you should not use this method at all Int_t Int_t z
TMethodCall * fMethodCall
Pointer to histogram used for visualisation.
virtual void SetRange(Double_t xmin, Double_t xmax)
Initialize the upper and lower bounds to draw the function.
Bool_t Contains(const char *pat, ECaseCompare cmp=kExact) const
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
virtual Double_t EvalPar(const Double_t *x, const Double_t *params=0)
Evaluate function with given coordinates and parameters.
double norm(double *x, double *p)
virtual Int_t GetNpar() const
Long64_t BinarySearch(Long64_t n, const T *array, T value)
virtual Int_t DistancetoPrimitive(Int_t px, Int_t py)
Compute distance from point px,py to a function.
virtual Double_t GetSave(const Double_t *x)
Get value corresponding to X in array of fSave values.
virtual Version_t ReadVersion(UInt_t *start=0, UInt_t *bcnt=0, const TClass *cl=0)=0
virtual void SetNpz(Int_t npz=30)
Set the number of points used to draw the function.
virtual void Warning(const char *method, const char *msgfmt,...) const
Issue warning message.