584class DifferentDimension:
public std::exception {};
585class DifferentNumberOfBins:
public std::exception {};
586class DifferentAxisLimits:
public std::exception {};
587class DifferentBinLimits:
public std::exception {};
588class DifferentLabels:
public std::exception {};
688 if (nbins <= 0) {
Warning(
"TH1",
"nbins is <=0 - set to nbins = 1"); nbins = 1; }
710 if (nbins <= 0) {
Warning(
"TH1",
"nbins is <=0 - set to nbins = 1"); nbins = 1; }
732 if (nbins <= 0) {
Warning(
"TH1",
"nbins is <=0 - set to nbins = 1"); nbins = 1; }
760 Draw(
b ?
b->GetDrawOption() :
"");
825 Error(
"Add",
"Attempt to add a non-existing function");
845 for (
Int_t i = 0; i < 10; ++i)
s1[i] = 0;
851 Int_t bin, binx, biny, binz;
856 for (binz = 0; binz < ncellsz; ++binz) {
858 for (biny = 0; biny < ncellsy; ++biny) {
860 for (binx = 0; binx < ncellsx; ++binx) {
864 bin = binx + ncellsx * (biny + ncellsy * binz);
916 Error(
"Add",
"Attempt to add a non-existing histogram");
927 }
catch(DifferentNumberOfBins&) {
929 Info(
"Add",
"Attempt to add histograms with different number of bins - trying to use TH1::Merge");
931 Error(
"Add",
"Attempt to add histograms with different number of bins : nbins h1 = %d , nbins h2 = %d",
GetNbinsX(),
h1->
GetNbinsX());
934 }
catch(DifferentAxisLimits&) {
936 Info(
"Add",
"Attempt to add histograms with different axis limits - trying to use TH1::Merge");
938 Warning(
"Add",
"Attempt to add histograms with different axis limits");
939 }
catch(DifferentBinLimits&) {
941 Info(
"Add",
"Attempt to add histograms with different bin limits - trying to use TH1::Merge");
943 Warning(
"Add",
"Attempt to add histograms with different bin limits");
944 }
catch(DifferentLabels&) {
947 Info(
"Add",
"Attempt to add histograms with different labels - trying to use TH1::Merge");
949 Info(
"Warning",
"Attempt to add histograms with different labels");
997 if (e1sq) w1 = 1. / e1sq;
1002 double sf = (s2[0] != 0) ? s2[1]/s2[0] : 1;
1006 if (e2sq) w2 = 1. / e2sq;
1011 double sf = (
s1[0] != 0) ?
s1[1]/
s1[0] : 1;
1016 double y = (w1*y1 + w2*y2)/(w1 + w2);
1019 double err2 = 1./(w1 + w2);
1020 if (err2 < 1.E-200) err2 = 0;
1036 if (i == 1)
s1[i] +=
c1*
c1*s2[i];
1037 else s1[i] +=
c1*s2[i];
1084 Error(
"Add",
"Attempt to add a non-existing histogram");
1092 if (
h1 == h2 &&
c2 < 0) {
c2 = 0; normWidth =
kTRUE;}
1101 }
catch(DifferentNumberOfBins&) {
1103 Info(
"Add",
"Attempt to add histograms with different number of bins - trying to use TH1::Merge");
1105 Error(
"Add",
"Attempt to add histograms with different number of bins : nbins h1 = %d , nbins h2 = %d",
GetNbinsX(),
h1->
GetNbinsX());
1108 }
catch(DifferentAxisLimits&) {
1110 Info(
"Add",
"Attempt to add histograms with different axis limits - trying to use TH1::Merge");
1112 Warning(
"Add",
"Attempt to add histograms with different axis limits");
1113 }
catch(DifferentBinLimits&) {
1115 Info(
"Add",
"Attempt to add histograms with different bin limits - trying to use TH1::Merge");
1117 Warning(
"Add",
"Attempt to add histograms with different bin limits");
1118 }
catch(DifferentLabels&) {
1121 Info(
"Add",
"Attempt to add histograms with different labels - trying to use TH1::Merge");
1123 Info(
"Warning",
"Attempt to add histograms with different labels");
1130 l.Add(
const_cast<TH1*
>(h2));
1152 Bool_t resetStats = (
c1*
c2 < 0) || normWidth;
1159 if (i == 1) s3[i] =
c1*
c1*
s1[i] +
c2*
c2*s2[i];
1161 else s3[i] =
c1*
s1[i] +
c2*s2[i];
1177 Int_t bin, binx, biny, binz;
1178 for (binz = 0; binz < nbinsz; ++binz) {
1180 for (biny = 0; biny < nbinsy; ++biny) {
1182 for (binx = 0; binx < nbinsx; ++binx) {
1184 bin =
GetBin(binx, biny, binz);
1205 if (e1sq) w1 = 1./ e1sq;
1209 double sf = (
s1[0] != 0) ?
s1[1]/
s1[0] : 1;
1213 if (e2sq) w2 = 1./ e2sq;
1217 double sf = (s2[0] != 0) ? s2[1]/s2[0] : 1;
1222 double y = (w1*y1 + w2*y2)/(w1 + w2);
1225 double err2 = 1./(w1 + w2);
1226 if (err2 < 1.E-200) err2 = 0;
1300 return ((next &&
x > 0.) || (!next &&
x <= 0.)) ? std::ldexp(std::copysign(1., f2), nn)
1301 : std::ldexp(std::copysign(1., f2), --nn);
1314 return (
Int_t)std::ldexp(1., nn);
1354 Double_t rr = (xhma - xhmi) / (xma - xmi);
1364 Int_t nbup = (xhma - xma) / bw;
1367 if (nbup != nbside) {
1369 xhma -= bw * (nbup - nbside);
1370 nb -= (nbup - nbside);
1374 Int_t nblw = (xmi - xhmi) / bw;
1377 if (nblw != nbside) {
1379 xhmi += bw * (nblw - nbside);
1380 nb -= (nblw - nbside);
1410 if (nbentries == 0) {
1420 if (nbentries < 0 && action == 0)
return 0;
1423 if (nbentries < 0) {
1424 nbentries = -nbentries;
1436 for (
Int_t i=1;i<nbentries;i++) {
1446 "inconsistency found by power-of-2 autobin algorithm: fallback to standard method");
1464 DoFillN(nbentries,&buffer[2],&buffer[1],2);
1498 if (nbentries < 0) {
1501 nbentries = -nbentries;
1533 if ( h2Array->
fN != fN ) {
1534 throw DifferentBinLimits();
1538 for (
int i = 0; i < fN; ++i ) {
1543 throw DifferentBinLimits();
1564 throw DifferentLabels();
1569 throw DifferentLabels();
1572 for (
int i = 1; i <= a1->
GetNbins(); ++i) {
1575 if (label1 != label2) {
1576 throw DifferentLabels();
1594 throw DifferentAxisLimits();
1607 ::Info(
"CheckEqualAxes",
"Axes have different number of bins : nbin1 = %d nbin2 = %d",a1->
GetNbins(),a2->
GetNbins() );
1612 }
catch (DifferentAxisLimits&) {
1613 ::Info(
"CheckEqualAxes",
"Axes have different limits");
1618 }
catch (DifferentBinLimits&) {
1619 ::Info(
"CheckEqualAxes",
"Axes have different bin limits");
1626 }
catch (DifferentLabels&) {
1627 ::Info(
"CheckEqualAxes",
"Axes have different labels");
1642 Int_t nbins1 = lastBin1-firstBin1 + 1;
1650 if (firstBin2 < lastBin2) {
1652 nbins2 = lastBin1-firstBin1 + 1;
1657 if (nbins1 != nbins2 ) {
1658 ::Info(
"CheckConsistentSubAxes",
"Axes have different number of bins");
1666 ::Info(
"CheckConsistentSubAxes",
"Axes have different limits");
1678 if (
h1 == h2)
return true;
1681 throw DifferentDimension();
1693 (dim > 1 && nbinsy != h2->
GetNbinsY()) ||
1694 (dim > 2 && nbinsz != h2->
GetNbinsZ()) ) {
1695 throw DifferentNumberOfBins();
2010 Int_t ndf = 0, igood = 0;
2018 printf(
"Chi2 = %f, Prob = %g, NDF = %d, igood = %d\n", chi2,prob,ndf,igood);
2021 if (ndf == 0)
return 0;
2069 Int_t i_start, i_end;
2070 Int_t j_start, j_end;
2071 Int_t k_start, k_end;
2100 Error(
"Chi2TestX",
"Histograms have different dimensions.");
2105 if (nbinx1 != nbinx2) {
2106 Error(
"Chi2TestX",
"different number of x channels");
2108 if (nbiny1 != nbiny2) {
2109 Error(
"Chi2TestX",
"different number of y channels");
2111 if (nbinz1 != nbinz2) {
2112 Error(
"Chi2TestX",
"different number of z channels");
2116 i_start = j_start = k_start = 1;
2147 ndf = (i_end - i_start + 1) * (j_end - j_start + 1) * (k_end - k_start + 1) - 1;
2154 if (scaledHistogram && !comparisonUU) {
2155 Info(
"Chi2TestX",
"NORM option should be used together with UU option. It is ignored");
2162 Double_t effEntries1 = (s[1] ? s[0] * s[0] / s[1] : 0.0);
2166 Double_t effEntries2 = (s[1] ? s[0] * s[0] / s[1] : 0.0);
2168 if (!comparisonUU && !comparisonUW && !comparisonWW ) {
2170 if (
TMath::Abs(sumBinContent1 - effEntries1) < 1) {
2171 if (
TMath::Abs(sumBinContent2 - effEntries2) < 1) comparisonUU =
true;
2172 else comparisonUW =
true;
2174 else comparisonWW =
true;
2178 if (
TMath::Abs(sumBinContent1 - effEntries1) >= 1) {
2179 Warning(
"Chi2TestX",
"First histogram is not unweighted and option UW has been requested");
2182 if ( (!scaledHistogram && comparisonUU) ) {
2183 if ( (
TMath::Abs(sumBinContent1 - effEntries1) >= 1) || (
TMath::Abs(sumBinContent2 - effEntries2) >= 1) ) {
2184 Warning(
"Chi2TestX",
"Both histograms are not unweighted and option UU has been requested");
2190 if (comparisonUU && scaledHistogram) {
2191 for (
Int_t i = i_start; i <= i_end; ++i) {
2192 for (
Int_t j = j_start; j <= j_end; ++j) {
2193 for (
Int_t k = k_start; k <= k_end; ++k) {
2202 if (e1sq > 0.0) cnt1 =
TMath::Floor(cnt1 * cnt1 / e1sq + 0.5);
2205 if (e2sq > 0.0) cnt2 =
TMath::Floor(cnt2 * cnt2 / e2sq + 0.5);
2216 if (sumw1 <= 0.0 || sumw2 <= 0.0) {
2217 Error(
"Chi2TestX",
"Cannot use option NORM when one histogram has all zero errors");
2222 for (
Int_t i = i_start; i <= i_end; ++i) {
2223 for (
Int_t j = j_start; j <= j_end; ++j) {
2224 for (
Int_t k = k_start; k <= k_end; ++k) {
2238 if (sum1 == 0.0 || sum2 == 0.0) {
2239 Error(
"Chi2TestX",
"one histogram is empty");
2243 if ( comparisonWW && ( sumw1 <= 0.0 && sumw2 <= 0.0 ) ){
2244 Error(
"Chi2TestX",
"Hist1 and Hist2 have both all zero errors\n");
2254 for (
Int_t i = i_start; i <= i_end; ++i) {
2255 for (
Int_t j = j_start; j <= j_end; ++j) {
2256 for (
Int_t k = k_start; k <= k_end; ++k) {
2263 if (scaledHistogram) {
2268 if (e1sq > 0) cnt1 =
TMath::Floor(cnt1 * cnt1 / e1sq + 0.5);
2271 if (e2sq > 0) cnt2 =
TMath::Floor(cnt2 * cnt2 / e2sq + 0.5);
2275 if (
Int_t(cnt1) == 0 &&
Int_t(cnt2) == 0) --ndf;
2282 if (res) res[i - i_start] = (cnt1 - nexp1) /
TMath::Sqrt(nexp1);
2291 Double_t delta = sum2 * cnt1 - sum1 * cnt2;
2292 chi2 += delta * delta / cntsum;
2297 chi2 /= sum1 * sum2;
2302 Info(
"Chi2TestX",
"There is a bin in h1 with less than 1 event.\n");
2306 Info(
"Chi2TestX",
"There is a bin in h2 with less than 1 event.\n");
2317 if ( comparisonUW ) {
2318 for (
Int_t i = i_start; i <= i_end; ++i) {
2319 for (
Int_t j = j_start; j <= j_end; ++j) {
2320 for (
Int_t k = k_start; k <= k_end; ++k) {
2329 if (cnt1 * cnt1 == 0 && cnt2 * cnt2 == 0) {
2335 if (cnt2 * cnt2 == 0 && e2sq == 0) {
2339 e2sq = sumw2 / sum2;
2344 Error(
"Chi2TestX",
"Hist2 has in bin (%d,%d,%d) zero content and zero errors\n", i, j, k);
2350 if (e2sq > 0 && cnt2 * cnt2 / e2sq < 10)
n++;
2352 Double_t var1 = sum2 * cnt2 - sum1 * e2sq;
2353 Double_t var2 = var1 * var1 + 4. * sum2 * sum2 * cnt1 * e2sq;
2358 while (var1 * var1 + cnt1 == 0 || var1 + var2 == 0) {
2361 var1 = sum2 * cnt2 - sum1 * e2sq;
2362 var2 = var1 * var1 + 4. * sum2 * sum2 * cnt1 * e2sq;
2366 while (var1 + var2 == 0) {
2369 var1 = sum2 * cnt2 - sum1 * e2sq;
2370 var2 = var1 * var1 + 4. * sum2 * sum2 * cnt1 * e2sq;
2371 while (var1 * var1 + cnt1 == 0 || var1 + var2 == 0) {
2374 var1 = sum2 * cnt2 - sum1 * e2sq;
2375 var2 = var1 * var1 + 4. * sum2 * sum2 * cnt1 * e2sq;
2380 Double_t probb = (var1 + var2) / (2. * sum2 * sum2);
2388 chi2 += delta1 * delta1 / nexp1;
2391 chi2 += delta2 * delta2 / e2sq;
2396 Double_t temp1 = sum2 * e2sq / var2;
2397 Double_t temp2 = 1.0 + (sum1 * e2sq - sum2 * cnt2) / var2;
2398 temp2 = temp1 * temp1 * sum1 * probb * (1.0 - probb) + temp2 * temp2 * e2sq / 4.0;
2411 Info(
"Chi2TestX",
"There is a bin in h1 with less than 1 event.\n");
2415 Info(
"Chi2TestX",
"There is a bin in h2 with less than 10 effective events.\n");
2425 for (
Int_t i = i_start; i <= i_end; ++i) {
2426 for (
Int_t j = j_start; j <= j_end; ++j) {
2427 for (
Int_t k = k_start; k <= k_end; ++k) {
2437 if (cnt1 * cnt1 == 0 && cnt2 * cnt2 == 0) {
2442 if (e1sq == 0 && e2sq == 0) {
2444 Error(
"Chi2TestX",
"h1 and h2 both have bin %d,%d,%d with all zero errors\n", i,j,k);
2449 Double_t delta = sum2 * cnt1 - sum1 * cnt2;
2450 chi2 += delta * delta /
sigma;
2453 Double_t temp = cnt1 * sum1 * e2sq + cnt2 * sum2 * e1sq;
2466 res[i - i_start] = z;
2469 if (e1sq > 0 && cnt1 * cnt1 / e1sq < 10)
m++;
2470 if (e2sq > 0 && cnt2 * cnt2 / e2sq < 10)
n++;
2476 Info(
"Chi2TestX",
"There is a bin in h1 with less than 10 effective events.\n");
2480 Info(
"Chi2TestX",
"There is a bin in h2 with less than 10 effective events.\n");
2497 Error(
"Chisquare",
"Function pointer is Null - return -1");
2543 Int_t nbins = nbinsx * nbinsy * nbinsz;
2548 for (
Int_t binz=1; binz <= nbinsz; ++binz) {
2549 for (
Int_t biny=1; biny <= nbinsy; ++biny) {
2550 for (
Int_t binx=1; binx <= nbinsx; ++binx) {
2553 if (onlyPositive &&
y < 0) {
2554 Error(
"ComputeIntegral",
"Bin content is negative - return a NaN value");
2565 Error(
"ComputeIntegral",
"Integral = zero");
return 0;
2620 hintegrated->
Reset();
2624 for (
Int_t binz = firstZ; binz <= lastZ; ++binz) {
2625 for (
Int_t biny = firstY; biny <= lastY; ++biny) {
2626 for (
Int_t binx = firstX; binx <= lastX; ++binx) {
2627 const Int_t bin = hintegrated->
GetBin(binx, biny, binz);
2638 for (
Int_t binz = lastZ; binz >= firstZ; --binz) {
2639 for (
Int_t biny = lastY; biny >= firstY; --biny) {
2640 for (
Int_t binx = lastX; binx >= firstX; --binx) {
2641 const Int_t bin = hintegrated->
GetBin(binx, biny, binz);
2672 ((
TH1&)obj).fDirectory = 0;
2686 delete [] ((
TH1&)obj).fBuffer;
2687 ((
TH1&)obj).fBuffer = 0;
2693 ((
TH1&)obj).fBuffer = buf;
2719 ((
TH1&)obj).fXaxis.SetParent(&obj);
2720 ((
TH1&)obj).fYaxis.SetParent(&obj);
2721 ((
TH1&)obj).fZaxis.SetParent(&obj);
2731 ((
TH1&)obj).fFunctions->UseRWLock();
2744 TH1* obj = (
TH1*)IsA()->GetNew()(0);
2759 oldparent = oldstats->GetParent();
2760 oldstats->SetParent(
nullptr);
2766 oldstats->SetParent(oldparent);
2778 if(newname && strlen(newname) ) {
2833 Error(
"Divide",
"Attempt to divide by a non-existing function");
2851 Int_t bin, binx, biny, binz;
2856 for (binz = 0; binz < nz; ++binz) {
2858 for (biny = 0; biny < ny; ++biny) {
2860 for (binx = 0; binx < nx; ++binx) {
2864 bin = binx + nx * (biny + ny * binz);
2901 Error(
"Divide",
"Input histogram passed does not exist (NULL).");
2910 }
catch(DifferentNumberOfBins&) {
2911 Error(
"Divide",
"Cannot divide histograms with different number of bins");
2913 }
catch(DifferentAxisLimits&) {
2914 Warning(
"Divide",
"Dividing histograms with different axis limits");
2915 }
catch(DifferentBinLimits&) {
2916 Warning(
"Divide",
"Dividing histograms with different bin limits");
2917 }
catch(DifferentLabels&) {
2918 Warning(
"Divide",
"Dividing histograms with different labels");
2973 Error(
"Divide",
"At least one of the input histograms passed does not exist (NULL).");
2983 }
catch(DifferentNumberOfBins&) {
2984 Error(
"Divide",
"Cannot divide histograms with different number of bins");
2986 }
catch(DifferentAxisLimits&) {
2987 Warning(
"Divide",
"Dividing histograms with different axis limits");
2988 }
catch(DifferentBinLimits&) {
2989 Warning(
"Divide",
"Dividing histograms with different bin limits");
2990 }
catch(DifferentLabels&) {
2991 Warning(
"Divide",
"Dividing histograms with different labels");
2996 Error(
"Divide",
"Coefficient of dividing histogram cannot be zero");
3033 fSumw2.
fArray[i] = c1sq * c2sq * (e1sq * b2sq + e2sq * b1sq) / (c2sq * c2sq * b2sq * b2sq);
3086 if (index>indb && index<indk) index = -1;
3092 if (!
gPad->IsEditable())
gROOT->MakeDefCanvas();
3094 if (
gPad->GetX1() == 0 &&
gPad->GetX2() == 1 &&
3095 gPad->GetY1() == 0 &&
gPad->GetY2() == 1 &&
3096 gPad->GetListOfPrimitives()->GetSize()==0) opt2.
Remove(index,4);
3103 gPad->IncrementPaletteColor(1, opt1);
3105 if (index>=0) opt2.
Remove(index,4);
3131 if (
gPad)
gPad->IncrementPaletteColor(1, opt);
3156 Error(
"DrawNormalized",
"Sum of weights is null. Cannot normalize histogram: %s",
GetName());
3168 if (opt.
IsNull() || opt ==
"SAME") opt +=
"HIST";
3203 Int_t range, stat, add;
3223 for (
Int_t binz = 1; binz <= nbinsz; ++binz) {
3225 for (
Int_t biny = 1; biny <= nbinsy; ++biny) {
3227 for (
Int_t binx = 1; binx <= nbinsx; ++binx) {
3327 for (
Int_t binx = 1; binx<=ndim[0]; binx++) {
3328 for (
Int_t biny=1; biny<=ndim[1]; biny++) {
3329 for (
Int_t binz=1; binz<=ndim[2]; binz++) {
3359 if (bin <0)
return -1;
3392 if (bin <0)
return -1;
3425 if (bin <0)
return -1;
3461 for (i=0;i<ntimes;i+=stride) {
3467 if (i < ntimes &&
fBuffer==0) {
3468 auto weights = w ? &w[i] :
nullptr;
3469 DoFillN((ntimes-i)/stride,&
x[i],weights,stride);
3489 for (i=0;i<ntimes;i+=stride) {
3491 if (bin <0)
continue;
3496 if (bin == 0 || bin > nbins) {
3529 Int_t bin, binx, ibin, loop;
3533 if (!
f1) {
Error(
"FillRandom",
"Unknown function: %s",fname);
return; }
3543 Info(
"FillRandom",
"Using function axis and range [%g,%g]",
xmin,
xmax);
3553 for (binx=1;binx<=nbinsx;binx++) {
3555 integral[binx] = integral[binx-1] + fint;
3559 if (integral[nbinsx] == 0 ) {
3561 Error(
"FillRandom",
"Integral = zero");
return;
3563 for (bin=1;bin<=nbinsx;bin++) integral[bin] /= integral[nbinsx];
3566 for (loop=0;loop<ntimes;loop++) {
3572 +xAxis->
GetBinWidth(ibin+
first)*(r1-integral[ibin])/(integral[ibin+1] - integral[ibin]);
3600 if (!
h) {
Error(
"FillRandom",
"Null histogram");
return; }
3602 Error(
"FillRandom",
"Histograms with different dimensions");
return;
3604 if (std::isnan(
h->ComputeIntegral(
true))) {
3605 Error(
"FillRandom",
"Histograms contains negative bins, does not represent probabilities");
3614 if (ntimes > 10*nbins) {
3618 if (sumw == 0)
return;
3621 Double_t mean =
h->RetrieveBinContent(bin)*ntimes/sumw;
3632 if (sumgen < ntimes) {
3634 for (i =
Int_t(sumgen+0.5); i < ntimes; ++i)
3640 else if (sumgen > ntimes) {
3642 i =
Int_t(sumgen+0.5);
3643 while( i > ntimes) {
3658 catch(std::exception&) {}
3662 if (
h->ComputeIntegral() ==0)
return;
3665 for (loop=0;loop<ntimes;loop++) {
3691 return binx + nx*biny;
3699 return binx + nx*(biny +ny*binz);
3724 return binx + nx*biny;
3732 return binx + nx*(biny +ny*binz);
3750 Warning(
"FindFirstBinAbove",
"Invalid axis number : %d, axis x assumed\n",axis);
3764 for (
Int_t binx = firstBin; binx <= lastBin; binx++) {
3765 for (
Int_t biny = 1; biny <= nbinsy; biny++) {
3766 for (
Int_t binz = 1; binz <= nbinsz; binz++) {
3772 else if (axis == 2) {
3776 for (
Int_t biny = firstBin; biny <= lastBin; biny++) {
3777 for (
Int_t binx = 1; binx <= nbinsx; binx++) {
3778 for (
Int_t binz = 1; binz <= nbinsz; binz++) {
3784 else if (axis == 3) {
3788 for (
Int_t binz = firstBin; binz <= lastBin; binz++) {
3789 for (
Int_t binx = 1; binx <= nbinsx; binx++) {
3790 for (
Int_t biny = 1; biny <= nbinsy; biny++) {
3814 Warning(
"FindFirstBinAbove",
"Invalid axis number : %d, axis x assumed\n",axis);
3828 for (
Int_t binx = lastBin; binx >= firstBin; binx--) {
3829 for (
Int_t biny = 1; biny <= nbinsy; biny++) {
3830 for (
Int_t binz = 1; binz <= nbinsz; binz++) {
3836 else if (axis == 2) {
3840 for (
Int_t biny = lastBin; biny >= firstBin; biny--) {
3841 for (
Int_t binx = 1; binx <= nbinsx; binx++) {
3842 for (
Int_t binz = 1; binz <= nbinsz; binz++) {
3848 else if (axis == 3) {
3852 for (
Int_t binz = lastBin; binz >= firstBin; binz--) {
3853 for (
Int_t binx = 1; binx <= nbinsx; binx++) {
3854 for (
Int_t biny = 1; biny <= nbinsy; biny++) {
3897 linear= (
char*)strstr(fname,
"++");
3901 TF1 f1(fname, fname, xxmin, xxmax);
3902 return Fit(&
f1,option,goption,xxmin,xxmax);
3905 TF2 f2(fname, fname);
3906 return Fit(&f2,option,goption,xxmin,xxmax);
3909 TF3 f3(fname, fname);
3910 return Fit(&f3,option,goption,xxmin,xxmax);
3915 if (!
f1) {
Printf(
"Unknown function: %s",fname);
return -1; }
3916 return Fit(
f1,option,goption,xxmin,xxmax);
4252 gROOT->MakeDefCanvas();
4255 Error(
"FitPanel",
"Unable to create a default canvas");
4262 if (handler && handler->
LoadPlugin() != -1) {
4264 Error(
"FitPanel",
"Unable to create the FitPanel");
4267 Error(
"FitPanel",
"Unable to find the FitPanel plug-in");
4324 asym->
Divide(top,bottom);
4338 for(
Int_t k=1; k<= zmax; k++){
4417 return (s[1] ? s[0]*s[0]/s[1] :
TMath::Abs(s[0]) );
4428 Info(
"SetHighlight",
"Supported only 1-D or 2-D histograms");
4433 Info(
"SetHighlight",
"Need to draw histogram first");
4447 return ((
TH1*)
this)->GetPainter()->GetObjectInfo(px,py);
4548 Error(
"GetQuantiles",
"Only available for 1-d histograms");
4558 Int_t nq = nprobSum;
4563 for (i=1;i<nq;i++) {
4568 for (i = 0; i < nq; i++) {
4570 while (ibin < nbins-1 &&
fIntegral[ibin+1] == prob[i]) {
4571 if (
fIntegral[ibin+2] == prob[i]) ibin++;
4579 if (!probSum)
delete [] prob;
4597 Double_t allcha, sumx, sumx2,
x, val, stddev, mean;
4608 allcha = sumx = sumx2 = 0;
4609 for (bin=hxfirst;bin<=hxlast;bin++) {
4612 if (val > valmax) valmax = val;
4617 if (allcha == 0)
return;
4619 stddev = sumx2/allcha - mean*mean;
4622 if (stddev == 0) stddev = binwidx*(hxlast-hxfirst+1)/4;
4629 Double_t constant = 0.5*(valmax+binwidx*allcha/(sqrtpi*stddev));
4637 if ((mean < xmin || mean >
xmax) && stddev > (
xmax-
xmin)) {
4658 Int_t nchanx = hxlast - hxfirst + 1;
4679 Int_t nchanx = hxlast - hxfirst + 1;
4682 if (nchanx <=1 || npar == 1) {
4705 const Int_t idim = 20;
4716 if (
m > idim ||
m >
n)
return;
4719 for (
l = 2;
l <=
m; ++
l) {
4721 b[
m +
l*20 - 21] = zero;
4728 for (k = hxfirst; k <= hxlast; ++k) {
4733 for (
l = 2;
l <=
m; ++
l) {
4736 da[
l-1] += power*yk;
4738 for (
l = 2;
l <=
m; ++
l) {
4740 b[
m +
l*20 - 21] += power;
4743 for (i = 3; i <=
m; ++i) {
4744 for (k = i; k <=
m; ++k) {
4745 b[k - 1 + (i-1)*20 - 21] =
b[k + (i-2)*20 - 21];
4750 for (i=0; i<
m; ++i)
a[i] = da[i];
4770 xbar = ybar = x2bar = xybar = 0;
4775 for (i = hxfirst; i <= hxlast; ++i) {
4779 if (yk <= 0) yk = 1
e-9;
4788 det = fn*x2bar - xbar*xbar;
4796 a0 = (x2bar*ybar - xbar*xybar) / det;
4797 a1 = (fn*xybar - xbar*ybar) / det;
4808 Int_t a_dim1, a_offset, b_dim1, b_offset;
4810 Int_t im1, jp1, nm1, nmi;
4816 b_offset = b_dim1 + 1;
4819 a_offset = a_dim1 + 1;
4822 if (idim <
n)
return;
4825 for (j = 1; j <=
n; ++j) {
4826 if (
a[j + j*a_dim1] <= 0) { ifail = -1;
return; }
4827 a[j + j*a_dim1] = one /
a[j + j*a_dim1];
4828 if (j ==
n)
continue;
4830 for (
l = jp1;
l <=
n; ++
l) {
4831 a[j +
l*a_dim1] =
a[j + j*a_dim1] *
a[
l + j*a_dim1];
4832 s1 = -
a[
l + (j+1)*a_dim1];
4833 for (i = 1; i <= j; ++i) {
s1 =
a[
l + i*a_dim1] *
a[i + (j+1)*a_dim1] +
s1; }
4834 a[
l + (j+1)*a_dim1] = -
s1;
4839 for (
l = 1;
l <= k; ++
l) {
4840 b[
l*b_dim1 + 1] =
a[a_dim1 + 1]*
b[
l*b_dim1 + 1];
4843 for (
l = 1;
l <= k; ++
l) {
4844 for (i = 2; i <=
n; ++i) {
4846 s21 = -
b[i +
l*b_dim1];
4847 for (j = 1; j <= im1; ++j) {
4848 s21 =
a[i + j*a_dim1]*
b[j +
l*b_dim1] + s21;
4850 b[i +
l*b_dim1] = -
a[i + i*a_dim1]*s21;
4853 for (i = 1; i <= nm1; ++i) {
4855 s22 = -
b[nmi +
l*b_dim1];
4856 for (j = 1; j <= i; ++j) {
4858 s22 =
a[nmi + nmjp1*a_dim1]*
b[nmjp1 +
l*b_dim1] + s22;
4860 b[nmi +
l*b_dim1] = -s22;
4898 if (binx < 0) binx = 0;
4899 if (binx > ofx) binx = ofx;
4914 binx = binglobal%nx;
4920 binx = binglobal%nx;
4921 biny = ((binglobal-binx)/nx)%ny;
4926 binx = binglobal%nx;
4927 biny = ((binglobal-binx)/nx)%ny;
4928 binz = ((binglobal-binx)/nx -biny)/ny;
4947 Error(
"GetRandom",
"Function only valid for 1-d histograms");
4957 integral = ((
TH1*)
this)->ComputeIntegral(
true);
4959 if (integral == 0)
return 0;
4998 if (bin < 0) bin = 0;
5024 Error(
"GetBinWithContent",
"function is only valid for 1-D histograms");
5030 if (firstx <= 0) firstx = 1;
5034 for (
Int_t i=firstx;i<=lastx;i++) {
5036 if (diff <= 0) {binx = i;
return diff;}
5037 if (diff < curmax && diff <= maxdiff) {curmax = diff, binminx=i;}
5072 return y0 + (
x-x0)*((y1-y0)/(
x1-x0));
5081 Error(
"Interpolate",
"This function must be called with 1 argument for a TH1");
5090 Error(
"Interpolate",
"This function must be called with 1 argument for a TH1");
5118 Int_t binx, biny, binz;
5140 Error(
"IsBinOverflow",
"Invalid axis value");
5150 Int_t binx, biny, binz;
5157 return (binx <= 0 || biny <= 0);
5159 return (binx <= 0 || biny <= 0 || binz <= 0);
5170 Error(
"IsBinUnderflow",
"Invalid axis value");