67 TAttLine::operator=(
sp);
68 TAttFill::operator=(
sp);
69 TAttMarker::operator=(
sp);
163 for (i=0;i<=
fNpx;i++) {
175 for (i=1;i<=
fNpx;i++) {
193 char *o = (
char *) opt.
Data();
198 if(o[i]==
'p') graph=
kTRUE ;
else o[
j++]=o[i];
227 if (
R__b.IsReading()) {
356 fValBeg(
valbeg), fValEnd(
valend), fBegCond(0), fEndCond(0)
366 for (
Int_t i=0; i<
n; ++i) {
386 fBegCond(0), fEndCond(0)
396 for (
Int_t i=0; i<
n; ++i) {
415 fBegCond(0), fEndCond(0)
425 for (
Int_t i=0; i<
n; ++i) {
441 const TF1 *func,
Int_t n,
const char *opt,
445 fBegCond(0), fEndCond(0)
458 for (
Int_t i=0; i<
n; ++i) {
474 const TGraph *
g,
const char *opt,
478 fBegCond(0), fEndCond(0)
490 g->GetPoint(i,
xx,
yy);
508 fBegCond(0), fEndCond(0)
519 fPoly[i].
X()=
h->GetXaxis()->GetBinCenter(i+1);
520 fPoly[i].
Y()=
h->GetBinContent(i+1);
533 :
TSpline(
sp3), fValBeg(
sp3.fValBeg), fValEnd(
sp3.fValEnd), fBegCond(
sp3.fBegCond), fEndCond(
sp3.fEndCond)
570 const char *b1 =
strstr(opt,
"b1");
572 const char *b2 =
strstr(opt,
"b2");
575 Error(
"SetCond",
"Cannot specify first and second derivative at first point");
577 Error(
"SetCond",
"Cannot specify first and second derivative at last point");
617 printf(
"1 TEST OF TSpline3 WITH NONEQUIDISTANT KNOTS\n");
634 for (i = 0; i <
n; ++i)
635 spline->GetCoeff(i,
hx,
a[i],
a[i+200],
a[i+400],
a[i+600]);
637 for (i = 0; i <
mm1; ++i)
diff[i] =
com[i] = 0;
638 for (k = 0; k <
n; ++k) {
639 for (i = 0; i < mm; ++i)
c[i] =
a[k+i*200];
640 printf(
" ---------------------------------------%3d --------------------------------------------\n",k+1);
645 for (i = 0; i < mm; ++i)
printf(
"%16.8f",
c[i]);
647 for (i = 0; i <
mm1; ++i)
650 for (i = 1; i < mm; ++i)
651 for (
jj = i;
jj < mm; ++
jj) {
655 for (i = 0; i < mm; ++i)
printf(
"%16.8f",
c[i]);
657 for (i = 0; i <
mm1; ++i)
658 if (!(k >=
n-2 && i != 0))
663 printf(
" MAXIMUM ABSOLUTE VALUES OF DIFFERENCES \n");
666 printf(
" MAXIMUM ABSOLUTE VALUES OF COEFFICIENTS \n");
672 for (
n = 10;
n <= 100;
n += 10) {
676 for (i = 0; i <
nm1; i += 2) {
688 for (i = 0; i <
n; ++i)
691 for (i = 0; i <
mm1; ++i)
693 for (k = 0; k <
n; ++k) {
694 for (i = 0; i < mm; ++i)
697 printf(
" ---------------------------------------%3d --------------------------------------------\n",k+1);
699 if (k ==
n-1)
printf(
"%16.8f\n",
c[0]);
703 for (i = 0; i < mm; ++i)
printf(
"%16.8f",
c[i]);
706 for (i = 0; i <
mm1; ++i)
710 for (i = 1; i < mm; ++i)
711 for (
jj = i;
jj < mm; ++
jj) {
716 for (i = 0; i < mm; ++i)
printf(
"%16.8f",
c[i]);
719 for (i = 0; i <
mm1; ++i)
720 if (!(k >=
n-2 && i != 0))
724 printf(
" MAXIMUM ABSOLUTE VALUES OF DIFFERENCES \n");
727 printf(
" MAXIMUM ABSOLUTE VALUES OF COEFFICIENTS \n");
770 "Binary search failed x(%d) = %f < x= %f < x(%d) = %f\n",
804 std::ofstream *
f =
new std::ofstream(
filename,std::ios::out);
814 char *dot =
strstr(buffer,
".");
816 strlcat(buffer,
"(double x) {\n",512);
825 snprintf(buffer,512,
" const double fX[%d] = {",
fNp);
830 for (i=0;i<
fNp;i++) {
835 if (i%5 == 4 || i ==
fNp-1) {
843 snprintf(buffer,512,
" const double fY[%d] = {",
fNp);
846 for (i=0;i<
fNp;i++) {
851 if (i%5 == 4 || i ==
fNp-1) {
859 snprintf(buffer,512,
" const double fB[%d] = {",
fNp);
862 for (i=0;i<
fNp;i++) {
867 if (i%5 == 4 || i ==
fNp-1) {
875 snprintf(buffer,512,
" const double fC[%d] = {",
fNp);
878 for (i=0;i<
fNp;i++) {
883 if (i%5 == 4 || i ==
fNp-1) {
891 snprintf(buffer,512,
" const double fD[%d] = {",
fNp);
894 for (i=0;i<
fNp;i++) {
899 if (i%5 == 4 || i ==
fNp-1) {
908 snprintf(buffer,512,
" int klow=0;\n");
911 snprintf(buffer,512,
" // If out of boundaries, extrapolate. It may be badly wrong\n");
912 snprintf(buffer,512,
" if(x<=fXmin) klow=0;\n");
914 snprintf(buffer,512,
" else if(x>=fXmax) klow=fNp-1;\n");
918 snprintf(buffer,512,
" if(fKstep) {\n");
921 snprintf(buffer,512,
" // Equidistant knots, use histogramming\n");
923 snprintf(buffer,512,
" klow = int((x-fXmin)/fDelta);\n");
925 snprintf(buffer,512,
" if (klow > fNp-1) klow = fNp-1;\n");
929 snprintf(buffer,512,
" int khig=fNp-1, khalf;\n");
932 snprintf(buffer,512,
" // Non equidistant knots, binary search\n");
934 snprintf(buffer,512,
" while(khig-klow>1)\n");
936 snprintf(buffer,512,
" if(x>fX[khalf=(klow+khig)/2]) klow=khalf;\n");
938 snprintf(buffer,512,
" else khig=khalf;\n");
944 snprintf(buffer,512,
" // Evaluate now\n");
946 snprintf(buffer,512,
" double dx=x-fX[klow];\n");
948 snprintf(buffer,512,
" return (fY[klow]+dx*(fB[klow]+dx*(fC[klow]+dx*fD[klow])));\n");
954 if (
f) {
f->close();
delete f;}
973 out <<
" spline3->SetNpx(" <<
fNpx <<
");\n";
976 out <<
" spline3->SetPoint(" << i <<
"," <<
fPoly[i].
X() <<
"," <<
fPoly[i].
Y() <<
");\n";
977 out <<
" spline3->SetPointCoeff(" << i <<
"," <<
fPoly[i].
B() <<
"," <<
fPoly[i].
C() <<
"," <<
fPoly[i].
D()
1084 for (
m=1;
m<
l; ++
m) {
1153 for (i=1; i<
fNp; ++i) {
1167 if (
R__b.IsReading()) {
1221 for (
Int_t i=0; i<
n; ++i) {
1255 for (
Int_t i=0; i<
n; ++i) {
1288 for (
Int_t i=0; i<
n; i++) {
1322 for (
Int_t i=0; i<
n; ++i) {
1359 g->GetPoint(i,
xx,
yy);
1392 fPoly[i+
beg].
X()=
h->GetXaxis()->GetBinCenter(i+1);
1440 const char *&
cb1,
const char *&
ce1,
1441 const char *&
cb2,
const char *&
ce2)
1471 const char *
cb1,
const char *
ce1,
const char *
cb2,
1546 "Binary search failed x(%d) = %f < x(%d) = %f\n",
1577 std::ofstream *
f =
new std::ofstream(
filename,std::ios::out);
1587 char *dot =
strstr(buffer,
".");
1589 strlcat(buffer,
"(double x) {\n",512);
1598 snprintf(buffer,512,
" const double fX[%d] = {",
fNp);
1603 for (i=0;i<
fNp;i++) {
1608 if (i%5 == 4 || i ==
fNp-1) {
1616 snprintf(buffer,512,
" const double fY[%d] = {",
fNp);
1619 for (i=0;i<
fNp;i++) {
1624 if (i%5 == 4 || i ==
fNp-1) {
1632 snprintf(buffer,512,
" const double fB[%d] = {",
fNp);
1635 for (i=0;i<
fNp;i++) {
1640 if (i%5 == 4 || i ==
fNp-1) {
1648 snprintf(buffer,512,
" const double fC[%d] = {",
fNp);
1651 for (i=0;i<
fNp;i++) {
1656 if (i%5 == 4 || i ==
fNp-1) {
1664 snprintf(buffer,512,
" const double fD[%d] = {",
fNp);
1667 for (i=0;i<
fNp;i++) {
1672 if (i%5 == 4 || i ==
fNp-1) {
1680 snprintf(buffer,512,
" const double fE[%d] = {",
fNp);
1683 for (i=0;i<
fNp;i++) {
1688 if (i%5 == 4 || i ==
fNp-1) {
1696 snprintf(buffer,512,
" const double fF[%d] = {",
fNp);
1699 for (i=0;i<
fNp;i++) {
1704 if (i%5 == 4 || i ==
fNp-1) {
1713 snprintf(buffer,512,
" int klow=0;\n");
1716 snprintf(buffer,512,
" // If out of boundaries, extrapolate. It may be badly wrong\n");
1717 snprintf(buffer,512,
" if(x<=fXmin) klow=0;\n");
1719 snprintf(buffer,512,
" else if(x>=fXmax) klow=fNp-1;\n");
1723 snprintf(buffer,512,
" if(fKstep) {\n");
1726 snprintf(buffer,512,
" // Equidistant knots, use histogramming\n");
1728 snprintf(buffer,512,
" klow = int((x-fXmin)/fDelta);\n");
1730 snprintf(buffer,512,
" if (klow > fNp-1) klow = fNp-1;\n");
1732 snprintf(buffer,512,
" } else {\n");
1734 snprintf(buffer,512,
" int khig=fNp-1, khalf;\n");
1737 snprintf(buffer,512,
" // Non equidistant knots, binary search\n");
1739 snprintf(buffer,512,
" while(khig-klow>1)\n");
1741 snprintf(buffer,512,
" if(x>fX[khalf=(klow+khig)/2]) klow=khalf;\n");
1743 snprintf(buffer,512,
" else khig=khalf;\n");
1749 snprintf(buffer,512,
" // Evaluate now\n");
1751 snprintf(buffer,512,
" double dx=x-fX[klow];\n");
1753 snprintf(buffer,512,
" return (fY[klow]+dx*(fB[klow]+dx*(fC[klow]+dx*(fD[klow]+dx*(fE[klow]+dx*fF[klow])))));\n");
1759 if (
f) {
f->close();
delete f;}
1773 TString::Format(
"\"%s\", %g, %g, (TF1 *)nullptr, %d, \"\", %g, %g, %g, %g",
1783 out <<
" spline5->SetNpx(" <<
fNpx <<
");\n";
1786 out <<
" spline5->SetPoint(" << i <<
"," <<
fPoly[i].
X() <<
"," <<
fPoly[i].
Y() <<
");\n";
1787 out <<
" spline5->SetPointCoeff(" << i <<
"," <<
fPoly[i].
B() <<
"," <<
fPoly[i].
C() <<
"," <<
fPoly[i].
D()
1788 <<
"," <<
fPoly[i].
E() <<
"," <<
fPoly[i].
F() <<
");\n";
1919 for (i = 1; i <
m; ++i) {
1949 for (i = 1; i <
fNp; ++i) {
1958 for (i = 2; i <
fNp; ++i) {
1976 for (i = 2; i <
m; ++i) {
1989 for (i=
fNp-3; i > 0; --i)
2007 for (i = 1; i <
m; ++i) {
2085 printf(
"1 TEST OF TSpline5 WITH NONEQUIDISTANT KNOTS\n");
2102 for (i = 0; i <
n; ++i)
2104 a[i+600],
a[i+800],
a[i+1000]);
2106 for (i = 0; i <
mm1; ++i)
diff[i] =
com[i] = 0;
2107 for (k = 0; k <
n; ++k) {
2108 for (i = 0; i < mm; ++i)
c[i] =
a[k+i*200];
2109 printf(
" ---------------------------------------%3d --------------------------------------------\n",k+1);
2114 for (i = 0; i < mm; ++i)
printf(
"%16.8f",
c[i]);
2116 for (i = 0; i <
mm1; ++i)
2119 for (i = 1; i < mm; ++i)
2120 for (
jj = i;
jj < mm; ++
jj) {
2122 c[
j-2] =
c[
j-1]*z+
c[
j-2];
2124 for (i = 0; i < mm; ++i)
printf(
"%16.8f",
c[i]);
2126 for (i = 0; i <
mm1; ++i)
2127 if (!(k >=
n-2 && i != 0))
2132 printf(
" MAXIMUM ABSOLUTE VALUES OF DIFFERENCES \n");
2135 printf(
" MAXIMUM ABSOLUTE VALUES OF COEFFICIENTS \n");
2141 for (
n = 10;
n <= 100;
n += 10) {
2145 for (i = 0; i <
nm1; i += 2) {
2157 for (i = 0; i <
n; ++i)
2159 a[i+600],
a[i+800],
a[i+1000]);
2161 for (i = 0; i <
mm1; ++i)
2163 for (k = 0; k <
n; ++k) {
2164 for (i = 0; i < mm; ++i)
2167 printf(
" ---------------------------------------%3d --------------------------------------------\n",k+1);
2169 if (k ==
n-1)
printf(
"%16.8f\n",
c[0]);
2171 if (k ==
n-1)
break;
2173 for (i = 0; i < mm; ++i)
printf(
"%16.8f",
c[i]);
2176 for (i = 0; i <
mm1; ++i)
2180 for (i = 1; i < mm; ++i)
2181 for (
jj = i;
jj < mm; ++
jj) {
2183 c[
j-2] =
c[
j-1]*z+
c[
j-2];
2186 for (i = 0; i < mm; ++i)
printf(
"%16.8f",
c[i]);
2189 for (i = 0; i <
mm1; ++i)
2190 if (!(k >=
n-2 && i != 0))
2194 printf(
" MAXIMUM ABSOLUTE VALUES OF DIFFERENCES \n");
2197 printf(
" MAXIMUM ABSOLUTE VALUES OF COEFFICIENTS \n");
2205 printf(
"1 TEST OF TSpline5 WITH NONEQUIDISTANT DOUBLE KNOTS\n");
2233 for (i = 0; i <
nn; ++i)
2235 a[i+600],
a[i+800],
a[i+1000]);
2237 for (i = 0; i <
mm1; ++i)
2239 for (k = 0; k <
nn; ++k) {
2240 for (i = 0; i < mm; ++i)
2242 printf(
" ---------------------------------------%3d --------------------------------------------\n",k+1);
2248 for (i = 0; i < mm; ++i)
printf(
"%16.8f",
c[i]);
2250 for (i = 0; i <
mm1; ++i)
2253 for (i = 1; i < mm; ++i)
2254 for (
jj = i;
jj < mm; ++
jj) {
2256 c[
j-2] =
c[
j-1]*z+
c[
j-2];
2258 for (i = 0; i < mm; ++i)
printf(
"%16.8f",
c[i]);
2260 for (i = 0; i <
mm1; ++i)
2261 if (!(k >=
nn-2 && i != 0))
2265 printf(
" MAXIMUM ABSOLUTE VALUES OF DIFFERENCES \n");
2266 for (i = 1; i <=
mm1; ++i) {
2272 printf(
" MAXIMUM ABSOLUTE VALUES OF COEFFICIENTS \n");
2276 for (
n = 10;
n <= 100;
n += 10) {
2281 for (i = 0; i <
n; ++i) {
2290 for (i = 0; i <
nn; ++i)
2292 a[i+600],
a[i+800],
a[i+1000]);
2294 for (i = 0; i <
mm1; ++i)
2296 for (k = 0; k <
nn; ++k) {
2297 for (i = 0; i < mm; ++i)
2300 printf(
" ---------------------------------------%3d --------------------------------------------\n",k+1);
2304 if (k ==
nn-1)
break;
2306 for (i = 0; i < mm; ++i)
printf(
"%16.8f",
c[i]);
2309 for (i = 0; i <
mm1; ++i)
2312 for (i = 1; i < mm; ++i) {
2313 for (
jj = i;
jj < mm; ++
jj) {
2315 c[
j-2] =
c[
j-1]*z+
c[
j-2];
2319 for (i = 0; i < mm; ++i)
printf(
"%16.8f",
c[i]);
2322 for (i = 0; i <
mm1; ++i)
2323 if (!(k >=
nn-2 && i != 0))
2327 printf(
" MAXIMUM ABSOLUTE VALUES OF DIFFERENCES \n");
2330 printf(
" MAXIMUM ABSOLUTE VALUES OF COEFFICIENTS \n");
2339 printf(
"1 TEST OF TSpline5 WITH NONEQUIDISTANT KNOTS,\n");
2340 printf(
" ONE DOUBLE, ONE TRIPLE KNOT\n");
2363 for (i = 0; i <
n; ++i)
2365 a[i+600],
a[i+800],
a[i+1000]);
2367 for (i = 0; i <
mm1; ++i)
2369 for (k = 0; k <
n; ++k) {
2370 for (i = 0; i < mm; ++i)
2372 printf(
" ---------------------------------------%3d --------------------------------------------\n",k+1);
2378 for (i = 0; i < mm; ++i)
printf(
"%16.8f",
c[i]);
2380 for (i = 0; i <
mm1; ++i)
2383 for (i = 1; i < mm; ++i)
2384 for (
jj = i;
jj < mm; ++
jj) {
2386 c[
j-2] =
c[
j-1]*z+
c[
j-2];
2388 for (i = 0; i < mm; ++i)
printf(
"%16.8f",
c[i]);
2390 for (i = 0; i <
mm1; ++i)
2391 if (!(k >=
n-2 && i != 0))
2395 printf(
" MAXIMUM ABSOLUTE VALUES OF DIFFERENCES \n");
2398 printf(
" MAXIMUM ABSOLUTE VALUES OF COEFFICIENTS \n");
2406 printf(
"1 TEST OF TSpline5 WITH NONEQUIDISTANT KNOTS,\n");
2407 printf(
" TWO DOUBLE, ONE TRIPLE KNOT\n");
2434 for (i = 0; i <
n; ++i)
2436 a[i+600],
a[i+800],
a[i+1000]);
2438 for (i = 0; i <
mm1; ++i)
2440 for (k = 0; k <
n; ++k) {
2441 for (i = 0; i < mm; ++i)
2443 printf(
" ---------------------------------------%3d --------------------------------------------\n",k+1);
2449 for (i = 0; i < mm; ++i)
printf(
"%16.8f",
c[i]);
2451 for (i = 0; i <
mm1; ++i)
2454 for (i = 1; i < mm; ++i)
2455 for (
jj = i;
jj < mm; ++
jj) {
2457 c[
j-2] =
c[
j-1]*z+
c[
j-2];
2459 for (i = 0; i < mm; ++i)
printf(
"%16.8f",
c[i]);
2461 for (i = 0; i <
mm1; ++i)
2462 if (!(k >=
n-2 && i != 0))
2466 printf(
" MAXIMUM ABSOLUTE VALUES OF DIFFERENCES \n");
2469 printf(
" MAXIMUM ABSOLUTE VALUES OF COEFFICIENTS \n");
2481 if (
R__b.IsReading()) {
int Int_t
Signed integer 4 bytes (int)
short Version_t
Class version identifier (short)
const char Option_t
Option string (const char)
ROOT::Detail::TRangeCast< T, true > TRangeDynCast
TRangeDynCast is an adapter class that allows the typed iteration through a TCollection.
winID h TVirtualViewer3D TVirtualGLPainter p
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void char Point_t Rectangle_t WindowAttributes_t Float_t Float_t Float_t Int_t Int_t UInt_t UInt_t Rectangle_t Int_t Int_t Window_t TString Int_t GCValues_t GetPrimarySelectionOwner GetDisplay GetScreen GetColormap GetNativeEvent const char const char dpyName wid window const char font_name cursor keysym reg const char only_if_exist regb h Point_t winding char text const char depth char const char Int_t count const char ColorStruct_t color const char filename
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void char Point_t Rectangle_t WindowAttributes_t Float_t r
R__EXTERN TSystem * gSystem
Fill Area Attributes class.
virtual void Streamer(TBuffer &)
virtual Color_t GetFillColor() const
Return the fill area color.
virtual Style_t GetFillStyle() const
Return the fill area style.
virtual void SetFillColor(Color_t fcolor)
Set the fill area color.
virtual void SetFillStyle(Style_t fstyle)
Set the fill area style.
virtual void SaveFillAttributes(std::ostream &out, const char *name, Int_t coldef=1, Int_t stydef=1001)
Save fill attributes as C++ statement(s) on output stream out.
virtual void Streamer(TBuffer &)
virtual Color_t GetLineColor() const
Return the line color.
virtual void SetLineStyle(Style_t lstyle)
Set the line style.
virtual Width_t GetLineWidth() const
Return the line width.
virtual void SetLineWidth(Width_t lwidth)
Set the line width.
virtual void SetLineColor(Color_t lcolor)
Set the line color.
virtual Style_t GetLineStyle() const
Return the line style.
virtual void SaveLineAttributes(std::ostream &out, const char *name, Int_t coldef=1, Int_t stydef=1, Int_t widdef=1)
Save line attributes as C++ statement(s) on output stream out.
virtual void SaveMarkerAttributes(std::ostream &out, const char *name, Int_t coldef=1, Int_t stydef=1, Int_t sizdef=1)
Save line attributes as C++ statement(s) on output stream out.
virtual Style_t GetMarkerStyle() const
Return the marker style.
virtual void SetMarkerColor(Color_t mcolor=1)
Set the marker color.
virtual Color_t GetMarkerColor() const
Return the marker color.
virtual Size_t GetMarkerSize() const
Return the marker size.
virtual void SetMarkerStyle(Style_t mstyle=1)
Set the marker style.
virtual void Streamer(TBuffer &)
virtual void SetMarkerSize(Size_t msize=1)
Set the marker size.
virtual void SetLimits(Double_t xmin, Double_t xmax)
Buffer base class used for serializing objects.
A TGraph is an object made of two arrays X and Y with npoints each.
void Paint(Option_t *chopt="") override
Draw this graph with its current attributes.
1-D histogram with a float per channel (see TH1 documentation)
TH1 is the base class of all histogram classes in ROOT.
virtual void SetDirectory(TDirectory *dir)
By default, when a histogram is created, it is added to the list of histogram objects in the current ...
virtual Double_t GetBinCenter(Int_t bin) const
Return bin center for 1D histogram.
Int_t DistancetoPrimitive(Int_t px, Int_t py) override
Compute distance from point px,py to a line.
@ kLogX
X-axis in log scale.
@ kNoStats
Don't draw stats box.
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...
void Paint(Option_t *option="") override
Control routine to paint any kind of histograms.
void ExecuteEvent(Int_t event, Int_t px, Int_t py) override
Execute action corresponding to one event.
The TNamed class is the base class for all named ROOT classes.
const char * GetName() const override
Returns name of object.
void Streamer(TBuffer &) override
Stream an object of class TObject.
const char * GetTitle() const override
Returns title of object.
TNamed & operator=(const TNamed &rhs)
TNamed assignment operator.
TObject & operator=(const TObject &rhs) noexcept
TObject assignment operator.
R__ALWAYS_INLINE Bool_t TestBit(UInt_t f) const
virtual void AppendPad(Option_t *option="")
Append graphics object to current pad.
void SetBit(UInt_t f, Bool_t set)
Set or unset the user status bits as specified in f.
virtual void Error(const char *method, const char *msgfmt,...) const
Issue error message.
static void SavePrimitiveDraw(std::ostream &out, const char *variable_name, Option_t *option=nullptr)
Save invocation of primitive Draw() method Skipped if option contains "nodraw" string.
static void SavePrimitiveConstructor(std::ostream &out, TClass *cl, const char *variable_name, const char *constructor_agrs="", Bool_t empty_line=kTRUE)
Save object constructor in the output stream "out".
Class to create third splines to interpolate knots Arbitrary conditions can be introduced for first a...
void SaveAs(const char *filename="", Option_t *option="") const override
Write this spline as a C++ function that can be executed without ROOT the name of the function is the...
Int_t fEndCond
0=no end cond, 1=first derivative, 2=second derivative
Int_t fBegCond
0=no beg cond, 1=first derivative, 2=second derivative
Int_t FindX(Double_t x) const
Find X.
void SavePrimitive(std::ostream &out, Option_t *option="") override
Save primitive as a C++ statement(s) on output stream out.
static void Test()
Test method for TSpline5.
Double_t fValBeg
Initial value of first or second derivative.
void BuildCoeff() override
Build coefficients.
void Streamer(TBuffer &) override
Stream an object of class TSpline3.
Double_t Eval(Double_t x) const override
Eval this spline at x.
Double_t fValEnd
End value of first or second derivative.
void SetCond(const char *opt)
Check the boundary conditions.
Double_t Derivative(Double_t x) const
Derivative.
TSplinePoly3 * fPoly
[fNp] Array of polynomial terms
TSpline3 & operator=(const TSpline3 &)
Assignment operator.
virtual void SetPoint(Int_t i, Double_t x, Double_t y)
Set point number i.
virtual void SetPointCoeff(Int_t i, Double_t b, Double_t c, Double_t d)
Set point coefficient number i.
Class to create quintic natural splines to interpolate knots Arbitrary conditions can be introduced f...
static void Test()
Test method for TSpline5.
Double_t Eval(Double_t x) const override
Eval this spline at x.
void Streamer(TBuffer &) override
Stream an object of class TSpline5.
void BuildCoeff() override
Algorithm 600, collected algorithms from acm.
TSplinePoly5 * fPoly
[fNp] Array of polynomial terms
virtual void SetPointCoeff(Int_t i, Double_t b, Double_t c, Double_t d, Double_t e, Double_t f)
Set point coefficient number i.
void SaveAs(const char *filename="", Option_t *option="") const override
Write this spline as a C++ function that can be executed without ROOT the name of the function is the...
Double_t Derivative(Double_t x) const
Derivative.
void BoundaryConditions(const char *opt, Int_t &beg, Int_t &end, const char *&cb1, const char *&ce1, const char *&cb2, const char *&ce2)
Check the boundary conditions and the amount of extra double knots needed.
void SavePrimitive(std::ostream &out, Option_t *option="") override
Save primitive as a C++ statement(s) on output stream out.
void SetBoundaries(Double_t b1, Double_t e1, Double_t b2, Double_t e2, const char *cb1, const char *ce1, const char *cb2, const char *ce2)
Set the boundary conditions at double/triple knots.
Int_t FindX(Double_t x) const
Find X.
virtual void SetPoint(Int_t i, Double_t x, Double_t y)
Set point number i.
TSpline5 & operator=(const TSpline5 &)
Assignment operator.
Double_t fC
Second order expansion coefficient : fC*2! is the second derivative at x.
Double_t Eval(Double_t x) const override
Double_t fD
Third order expansion coefficient : fD*3! is the third derivative at x.
Double_t Derivative(Double_t x) const
Double_t fB
First order expansion coefficient : fB*1! is the first derivative at x.
void CopyPoly(TSplinePoly3 const &other)
Utility called by the copy constructors and = operator.
TSplinePoly3 & operator=(TSplinePoly3 const &other)
Assignment operator.
void Streamer(TBuffer &) override
Stream an object of class TObject.
void CopyPoly(TSplinePoly5 const &other)
Utility called by the copy constructors and = operator.
TSplinePoly5 & operator=(TSplinePoly5 const &other)
Assignment operator.
Double_t fF
Fifth order expansion coefficient : fF*5! is the fifth derivative at x.
Double_t Derivative(Double_t x) const
void Streamer(TBuffer &) override
Stream an object of class TObject.
Double_t fB
First order expansion coefficient : fB*1! is the first derivative at x.
Double_t fC
Second order expansion coefficient : fC*2! is the second derivative at x.
Double_t fD
Third order expansion coefficient : fD*3! is the third derivative at x.
Double_t Eval(Double_t x) const override
Double_t fE
Fourth order expansion coefficient : fE*4! is the fourth derivative at x.
Base class for TSpline knot.
void CopyPoly(TSplinePoly const &other)
Utility called by the copy constructors and = operator.
TSplinePoly & operator=(TSplinePoly const &other)
Assignment operator.
Double_t fY
Constant term.
Base class for spline implementation containing the Draw/Paint methods.
void Streamer(TBuffer &) override
Stream an object of class TSpline.
virtual Double_t Eval(Double_t x) const =0
void ExecuteEvent(Int_t event, Int_t px, Int_t py) override
Execute action corresponding to one event.
TGraph * fGraph
Graph for drawing the knots.
Int_t DistancetoPrimitive(Int_t px, Int_t py) override
Compute distance from point px,py to a spline.
Double_t fXmin
Minimum value of abscissa.
TClass * IsA() const override
TH1F * fHistogram
Temporary histogram.
void Draw(Option_t *option="") override
Draw this function with its current attributes.
Double_t fDelta
Distance between equidistant knots.
virtual void GetKnot(Int_t i, Double_t &x, Double_t &y) const =0
TSpline & operator=(const TSpline &)
Assignment operator.
void Paint(Option_t *option="") override
Paint this function with its current attributes.
Bool_t fKstep
True of equidistant knots.
Int_t fNp
Number of knots.
~TSpline() override
Destructor.
Double_t fXmax
Maximum value of abscissa.
Int_t fNpx
Number of points used for graphical representation.
void ToLower()
Change string to lower-case.
TString & ReplaceSpecialCppChars()
Find special characters which are typically used in printf() calls and replace them by appropriate es...
const char * Data() const
static TString Format(const char *fmt,...)
Static method which formats a string using a printf style format descriptor and return a TString.
Bool_t Contains(const char *pat, ECaseCompare cmp=kExact) const
virtual Bool_t AccessPathName(const char *path, EAccessMode mode=kFileExists)
Returns FALSE if one can access a file using the specified access mode.
Short_t Max(Short_t a, Short_t b)
Returns the largest of a and b.
Int_t FloorNint(Double_t x)
Returns the nearest integer of TMath::Floor(x).
Short_t Min(Short_t a, Short_t b)
Returns the smallest of a and b.
Double_t Cos(Double_t)
Returns the cosine of an angle of x radians.
Double_t Sin(Double_t)
Returns the sine of an angle of x radians.
Double_t Log10(Double_t x)
Returns the common (base-10) logarithm of x.
Short_t Abs(Short_t d)
Returns the absolute value of parameter Short_t d.