682 if (nbins <= 0) {
Warning(
"TH1",
"nbins is <=0 - set to nbins = 1"); nbins = 1; }
704 if (nbins <= 0) {
Warning(
"TH1",
"nbins is <=0 - set to nbins = 1"); nbins = 1; }
726 if (nbins <= 0) {
Warning(
"TH1",
"nbins is <=0 - set to nbins = 1"); nbins = 1; }
755 Draw(
b ?
b->GetDrawOption() :
"");
820 Error(
"Add",
"Attempt to add a non-existing function");
840 for (
Int_t i = 0; i < 10; ++i)
s1[i] = 0;
846 Int_t bin, binx, biny, binz;
851 for (binz = 0; binz < ncellsz; ++binz) {
853 for (biny = 0; biny < ncellsy; ++biny) {
855 for (binx = 0; binx < ncellsx; ++binx) {
859 bin = binx + ncellsx * (biny + ncellsy * binz);
911 Error(
"Add",
"Attempt to add a non-existing histogram");
924 Info(
"Add",
"Attempt to add histograms with different number of bins - trying to use TH1::Merge");
926 Error(
"Add",
"Attempt to add histograms with different number of bins : nbins h1 = %d , nbins h2 = %d",
GetNbinsX(),
h1->
GetNbinsX());
931 Info(
"Add",
"Attempt to add histograms with different axis limits - trying to use TH1::Merge");
933 Warning(
"Add",
"Attempt to add histograms with different axis limits");
936 Info(
"Add",
"Attempt to add histograms with different bin limits - trying to use TH1::Merge");
938 Warning(
"Add",
"Attempt to add histograms with different bin limits");
942 Info(
"Add",
"Attempt to add histograms with different labels - trying to use TH1::Merge");
944 Info(
"Warning",
"Attempt to add histograms with different labels");
949 l.Add(
const_cast<TH1*
>(
h1));
992 if (e1sq) w1 = 1. / e1sq;
997 double sf = (s2[0] != 0) ? s2[1]/s2[0] : 1;
1001 if (e2sq) w2 = 1. / e2sq;
1006 double sf = (
s1[0] != 0) ?
s1[1]/
s1[0] : 1;
1011 double y = (w1*
y1 + w2*
y2)/(w1 + w2);
1014 double err2 = 1./(w1 + w2);
1015 if (err2 < 1.E-200) err2 = 0;
1031 if (i == 1)
s1[i] +=
c1*
c1*s2[i];
1032 else s1[i] +=
c1*s2[i];
1079 Error(
"Add",
"Attempt to add a non-existing histogram");
1087 if (
h1 == h2 &&
c2 < 0) {
c2 = 0; normWidth =
kTRUE;}
1098 Info(
"Add",
"Attempt to add histograms with different number of bins - trying to use TH1::Merge");
1100 Error(
"Add",
"Attempt to add histograms with different number of bins : nbins h1 = %d , nbins h2 = %d",
GetNbinsX(),
h1->
GetNbinsX());
1105 Info(
"Add",
"Attempt to add histograms with different axis limits - trying to use TH1::Merge");
1107 Warning(
"Add",
"Attempt to add histograms with different axis limits");
1110 Info(
"Add",
"Attempt to add histograms with different bin limits - trying to use TH1::Merge");
1112 Warning(
"Add",
"Attempt to add histograms with different bin limits");
1116 Info(
"Add",
"Attempt to add histograms with different labels - trying to use TH1::Merge");
1118 Info(
"Warning",
"Attempt to add histograms with different labels");
1124 l.Add(
const_cast<TH1*
>(
h1));
1125 l.Add(
const_cast<TH1*
>(h2));
1147 Bool_t resetStats = (
c1*
c2 < 0) || normWidth;
1154 if (i == 1) s3[i] =
c1*
c1*
s1[i] +
c2*
c2*s2[i];
1156 else s3[i] =
c1*
s1[i] +
c2*s2[i];
1172 Int_t bin, binx, biny, binz;
1173 for (binz = 0; binz < nbinsz; ++binz) {
1175 for (biny = 0; biny < nbinsy; ++biny) {
1177 for (binx = 0; binx < nbinsx; ++binx) {
1179 bin =
GetBin(binx, biny, binz);
1200 if (e1sq) w1 = 1./ e1sq;
1204 double sf = (
s1[0] != 0) ?
s1[1]/
s1[0] : 1;
1208 if (e2sq) w2 = 1./ e2sq;
1212 double sf = (s2[0] != 0) ? s2[1]/s2[0] : 1;
1217 double y = (w1*
y1 + w2*
y2)/(w1 + w2);
1220 double err2 = 1./(w1 + w2);
1221 if (err2 < 1.E-200) err2 = 0;
1295 return ((next &&
x > 0.) || (!next &&
x <= 0.)) ? std::ldexp(std::copysign(1., f2), nn)
1296 : std::ldexp(std::copysign(1., f2), --nn);
1309 return (
Int_t)std::ldexp(1., nn);
1349 Double_t rr = (xhma - xhmi) / (xma - xmi);
1359 Int_t nbup = (xhma - xma) / bw;
1362 if (nbup != nbside) {
1364 xhma -= bw * (nbup - nbside);
1365 nb -= (nbup - nbside);
1369 Int_t nblw = (xmi - xhmi) / bw;
1372 if (nblw != nbside) {
1374 xhmi += bw * (nblw - nbside);
1375 nb -= (nblw - nbside);
1405 if (nbentries == 0) {
1415 if (nbentries < 0 && action == 0)
return 0;
1418 if (nbentries < 0) {
1419 nbentries = -nbentries;
1431 for (
Int_t i=0;i<nbentries;i++) {
1434 if (!std::isfinite(
x))
continue;
1443 "inconsistency found by power-of-2 autobin algorithm: fallback to standard method");
1461 DoFillN(nbentries,&buffer[2],&buffer[1],2);
1495 if (nbentries < 0) {
1498 nbentries = -nbentries;
1530 if ( h2Array->
fN != fN ) {
1535 for (
int i = 0; i < fN; ++i ) {
1569 for (
int i = 1; i <= a1->
GetNbins(); ++i) {
1572 if (label1 != label2) {
1604 ::Info(
"CheckEqualAxes",
"Axes have different number of bins : nbin1 = %d nbin2 = %d",a1->
GetNbins(),a2->
GetNbins() );
1610 ::Info(
"CheckEqualAxes",
"Axes have different limits");
1616 ::Info(
"CheckEqualAxes",
"Axes have different bin limits");
1624 ::Info(
"CheckEqualAxes",
"Axes have different labels");
1639 Int_t nbins1 = lastBin1-firstBin1 + 1;
1647 if (firstBin2 < lastBin2) {
1649 nbins2 = lastBin1-firstBin1 + 1;
1654 if (nbins1 != nbins2 ) {
1655 ::Info(
"CheckConsistentSubAxes",
"Axes have different number of bins");
1663 ::Info(
"CheckConsistentSubAxes",
"Axes have different limits");
1675 if (
h1 == h2)
return true;
1690 (dim > 1 && nbinsy != h2->
GetNbinsY()) ||
1691 (dim > 2 && nbinsz != h2->
GetNbinsZ()) ) {
2007 Int_t ndf = 0, igood = 0;
2015 printf(
"Chi2 = %f, Prob = %g, NDF = %d, igood = %d\n", chi2,prob,ndf,igood);
2018 if (ndf == 0)
return 0;
2066 Int_t i_start, i_end;
2067 Int_t j_start, j_end;
2068 Int_t k_start, k_end;
2097 Error(
"Chi2TestX",
"Histograms have different dimensions.");
2102 if (nbinx1 != nbinx2) {
2103 Error(
"Chi2TestX",
"different number of x channels");
2105 if (nbiny1 != nbiny2) {
2106 Error(
"Chi2TestX",
"different number of y channels");
2108 if (nbinz1 != nbinz2) {
2109 Error(
"Chi2TestX",
"different number of z channels");
2113 i_start = j_start = k_start = 1;
2144 ndf = (i_end - i_start + 1) * (j_end - j_start + 1) * (k_end - k_start + 1) - 1;
2151 if (scaledHistogram && !comparisonUU) {
2152 Info(
"Chi2TestX",
"NORM option should be used together with UU option. It is ignored");
2159 Double_t effEntries1 = (
s[1] ?
s[0] *
s[0] /
s[1] : 0.0);
2163 Double_t effEntries2 = (
s[1] ?
s[0] *
s[0] /
s[1] : 0.0);
2165 if (!comparisonUU && !comparisonUW && !comparisonWW ) {
2167 if (
TMath::Abs(sumBinContent1 - effEntries1) < 1) {
2168 if (
TMath::Abs(sumBinContent2 - effEntries2) < 1) comparisonUU =
true;
2169 else comparisonUW =
true;
2171 else comparisonWW =
true;
2175 if (
TMath::Abs(sumBinContent1 - effEntries1) >= 1) {
2176 Warning(
"Chi2TestX",
"First histogram is not unweighted and option UW has been requested");
2179 if ( (!scaledHistogram && comparisonUU) ) {
2180 if ( (
TMath::Abs(sumBinContent1 - effEntries1) >= 1) || (
TMath::Abs(sumBinContent2 - effEntries2) >= 1) ) {
2181 Warning(
"Chi2TestX",
"Both histograms are not unweighted and option UU has been requested");
2187 if (comparisonUU && scaledHistogram) {
2188 for (
Int_t i = i_start; i <= i_end; ++i) {
2189 for (
Int_t j = j_start; j <= j_end; ++j) {
2190 for (
Int_t k = k_start; k <= k_end; ++k) {
2199 if (e1sq > 0.0) cnt1 =
TMath::Floor(cnt1 * cnt1 / e1sq + 0.5);
2202 if (e2sq > 0.0) cnt2 =
TMath::Floor(cnt2 * cnt2 / e2sq + 0.5);
2213 if (sumw1 <= 0.0 || sumw2 <= 0.0) {
2214 Error(
"Chi2TestX",
"Cannot use option NORM when one histogram has all zero errors");
2219 for (
Int_t i = i_start; i <= i_end; ++i) {
2220 for (
Int_t j = j_start; j <= j_end; ++j) {
2221 for (
Int_t k = k_start; k <= k_end; ++k) {
2235 if (sum1 == 0.0 || sum2 == 0.0) {
2236 Error(
"Chi2TestX",
"one histogram is empty");
2240 if ( comparisonWW && ( sumw1 <= 0.0 && sumw2 <= 0.0 ) ){
2241 Error(
"Chi2TestX",
"Hist1 and Hist2 have both all zero errors\n");
2251 for (
Int_t i = i_start; i <= i_end; ++i) {
2252 for (
Int_t j = j_start; j <= j_end; ++j) {
2253 for (
Int_t k = k_start; k <= k_end; ++k) {
2260 if (scaledHistogram) {
2265 if (e1sq > 0) cnt1 =
TMath::Floor(cnt1 * cnt1 / e1sq + 0.5);
2268 if (e2sq > 0) cnt2 =
TMath::Floor(cnt2 * cnt2 / e2sq + 0.5);
2272 if (
Int_t(cnt1) == 0 &&
Int_t(cnt2) == 0) --ndf;
2279 if (res) res[i - i_start] = (cnt1 - nexp1) /
TMath::Sqrt(nexp1);
2288 Double_t delta = sum2 * cnt1 - sum1 * cnt2;
2289 chi2 += delta * delta / cntsum;
2294 chi2 /= sum1 * sum2;
2299 Info(
"Chi2TestX",
"There is a bin in h1 with less than 1 event.\n");
2303 Info(
"Chi2TestX",
"There is a bin in h2 with less than 1 event.\n");
2314 if ( comparisonUW ) {
2315 for (
Int_t i = i_start; i <= i_end; ++i) {
2316 for (
Int_t j = j_start; j <= j_end; ++j) {
2317 for (
Int_t k = k_start; k <= k_end; ++k) {
2326 if (cnt1 * cnt1 == 0 && cnt2 * cnt2 == 0) {
2332 if (cnt2 * cnt2 == 0 && e2sq == 0) {
2336 e2sq = sumw2 / sum2;
2341 Error(
"Chi2TestX",
"Hist2 has in bin (%d,%d,%d) zero content and zero errors\n", i, j, k);
2347 if (e2sq > 0 && cnt2 * cnt2 / e2sq < 10)
n++;
2349 Double_t var1 = sum2 * cnt2 - sum1 * e2sq;
2350 Double_t var2 = var1 * var1 + 4. * sum2 * sum2 * cnt1 * e2sq;
2355 while (var1 * var1 + cnt1 == 0 || var1 + var2 == 0) {
2358 var1 = sum2 * cnt2 - sum1 * e2sq;
2359 var2 = var1 * var1 + 4. * sum2 * sum2 * cnt1 * e2sq;
2363 while (var1 + var2 == 0) {
2366 var1 = sum2 * cnt2 - sum1 * e2sq;
2367 var2 = var1 * var1 + 4. * sum2 * sum2 * cnt1 * e2sq;
2368 while (var1 * var1 + cnt1 == 0 || var1 + var2 == 0) {
2371 var1 = sum2 * cnt2 - sum1 * e2sq;
2372 var2 = var1 * var1 + 4. * sum2 * sum2 * cnt1 * e2sq;
2377 Double_t probb = (var1 + var2) / (2. * sum2 * sum2);
2385 chi2 += delta1 * delta1 / nexp1;
2388 chi2 += delta2 * delta2 / e2sq;
2393 Double_t temp1 = sum2 * e2sq / var2;
2394 Double_t temp2 = 1.0 + (sum1 * e2sq - sum2 * cnt2) / var2;
2395 temp2 = temp1 * temp1 * sum1 * probb * (1.0 - probb) + temp2 * temp2 * e2sq / 4.0;
2408 Info(
"Chi2TestX",
"There is a bin in h1 with less than 1 event.\n");
2412 Info(
"Chi2TestX",
"There is a bin in h2 with less than 10 effective events.\n");
2422 for (
Int_t i = i_start; i <= i_end; ++i) {
2423 for (
Int_t j = j_start; j <= j_end; ++j) {
2424 for (
Int_t k = k_start; k <= k_end; ++k) {
2434 if (cnt1 * cnt1 == 0 && cnt2 * cnt2 == 0) {
2439 if (e1sq == 0 && e2sq == 0) {
2441 Error(
"Chi2TestX",
"h1 and h2 both have bin %d,%d,%d with all zero errors\n", i,j,k);
2446 Double_t delta = sum2 * cnt1 - sum1 * cnt2;
2447 chi2 += delta * delta /
sigma;
2450 Double_t temp = cnt1 * sum1 * e2sq + cnt2 * sum2 * e1sq;
2463 res[i - i_start] = z;
2466 if (e1sq > 0 && cnt1 * cnt1 / e1sq < 10)
m++;
2467 if (e2sq > 0 && cnt2 * cnt2 / e2sq < 10)
n++;
2473 Info(
"Chi2TestX",
"There is a bin in h1 with less than 10 effective events.\n");
2477 Info(
"Chi2TestX",
"There is a bin in h2 with less than 10 effective events.\n");
2494 Error(
"Chisquare",
"Function pointer is Null - return -1");
2540 Int_t nbins = nbinsx * nbinsy * nbinsz;
2545 for (
Int_t binz=1; binz <= nbinsz; ++binz) {
2546 for (
Int_t biny=1; biny <= nbinsy; ++biny) {
2547 for (
Int_t binx=1; binx <= nbinsx; ++binx) {
2550 if (onlyPositive &&
y < 0) {
2551 Error(
"ComputeIntegral",
"Bin content is negative - return a NaN value");
2562 Error(
"ComputeIntegral",
"Integral = zero");
return 0;
2617 hintegrated->
Reset();
2621 for (
Int_t binz = firstZ; binz <= lastZ; ++binz) {
2622 for (
Int_t biny = firstY; biny <= lastY; ++biny) {
2623 for (
Int_t binx = firstX; binx <= lastX; ++binx) {
2624 const Int_t bin = hintegrated->
GetBin(binx, biny, binz);
2635 for (
Int_t binz = lastZ; binz >= firstZ; --binz) {
2636 for (
Int_t biny = lastY; biny >= firstY; --biny) {
2637 for (
Int_t binx = lastX; binx >= firstX; --binx) {
2638 const Int_t bin = hintegrated->
GetBin(binx, biny, binz);
2668 ((
TH1&)obj).fDirectory->Remove(&obj);
2669 ((
TH1&)obj).fDirectory =
nullptr;
2683 delete [] ((
TH1&)obj).fBuffer;
2684 ((
TH1&)obj).fBuffer =
nullptr;
2690 ((
TH1&)obj).fBuffer = buf;
2716 ((
TH1&)obj).fXaxis.SetParent(&obj);
2717 ((
TH1&)obj).fYaxis.SetParent(&obj);
2718 ((
TH1&)obj).fZaxis.SetParent(&obj);
2728 ((
TH1&)obj).fFunctions->UseRWLock();
2756 oldparent = oldstats->GetParent();
2757 oldstats->SetParent(
nullptr);
2763 oldstats->SetParent(oldparent);
2775 if(newname && strlen(newname) ) {
2830 Error(
"Divide",
"Attempt to divide by a non-existing function");
2848 Int_t bin, binx, biny, binz;
2853 for (binz = 0; binz < nz; ++binz) {
2855 for (biny = 0; biny < ny; ++biny) {
2857 for (binx = 0; binx < nx; ++binx) {
2861 bin = binx + nx * (biny + ny * binz);
2898 Error(
"Divide",
"Input histogram passed does not exist (NULL).");
2908 Error(
"Divide",
"Cannot divide histograms with different number of bins");
2911 Warning(
"Divide",
"Dividing histograms with different axis limits");
2913 Warning(
"Divide",
"Dividing histograms with different bin limits");
2915 Warning(
"Divide",
"Dividing histograms with different labels");
2970 Error(
"Divide",
"At least one of the input histograms passed does not exist (NULL).");
2981 Error(
"Divide",
"Cannot divide histograms with different number of bins");
2984 Warning(
"Divide",
"Dividing histograms with different axis limits");
2986 Warning(
"Divide",
"Dividing histograms with different bin limits");
2988 Warning(
"Divide",
"Dividing histograms with different labels");
2993 Error(
"Divide",
"Coefficient of dividing histogram cannot be zero");
3030 fSumw2.
fArray[i] = c1sq * c2sq * (e1sq * b2sq + e2sq * b1sq) / (c2sq * c2sq * b2sq * b2sq);
3089 if (!
gPad->IsEditable())
gROOT->MakeDefCanvas();
3091 if (
gPad->GetX1() == 0 &&
gPad->GetX2() == 1 &&
3092 gPad->GetY1() == 0 &&
gPad->GetY2() == 1 &&
3100 gPad->IncrementPaletteColor(1, opt1);
3125 if (name_postfix) newName.
Form(
"%s%s",
GetName(), name_postfix);
3129 if (
gPad)
gPad->IncrementPaletteColor(1, opt);
3154 Error(
"DrawNormalized",
"Sum of weights is null. Cannot normalize histogram: %s",
GetName());
3166 if (opt.
IsNull() || opt ==
"SAME") opt +=
"HIST";
3201 Int_t range, stat, add;
3221 for (
Int_t binz = 1; binz <= nbinsz; ++binz) {
3223 for (
Int_t biny = 1; biny <= nbinsy; ++biny) {
3225 for (
Int_t binx = 1; binx <= nbinsx; ++binx) {
3325 for (
Int_t binx = 1; binx<=ndim[0]; binx++) {
3326 for (
Int_t biny=1; biny<=ndim[1]; biny++) {
3327 for (
Int_t binz=1; binz<=ndim[2]; binz++) {
3357 if (bin <0)
return -1;
3390 if (bin <0)
return -1;
3423 if (bin <0)
return -1;
3459 for (i=0;i<ntimes;i+=stride) {
3466 auto weights =
w ? &
w[i] :
nullptr;
3467 DoFillN((ntimes-i)/stride,&
x[i],weights,stride);
3487 for (i=0;i<ntimes;i+=stride) {
3489 if (bin <0)
continue;
3494 if (bin == 0 || bin > nbins) {
3527 Int_t bin, binx, ibin, loop;
3531 if (!
f1) {
Error(
"FillRandom",
"Unknown function: %s",fname);
return; }
3541 Info(
"FillRandom",
"Using function axis and range [%g,%g]",
xmin,
xmax);
3551 for (binx=1;binx<=nbinsx;binx++) {
3553 integral[binx] = integral[binx-1] + fint;
3557 if (integral[nbinsx] == 0 ) {
3559 Error(
"FillRandom",
"Integral = zero");
return;
3561 for (bin=1;bin<=nbinsx;bin++) integral[bin] /= integral[nbinsx];
3564 for (loop=0;loop<ntimes;loop++) {
3570 +xAxis->
GetBinWidth(ibin+
first)*(r1-integral[ibin])/(integral[ibin+1] - integral[ibin]);
3598 if (!
h) {
Error(
"FillRandom",
"Null histogram");
return; }
3600 Error(
"FillRandom",
"Histograms with different dimensions");
return;
3602 if (std::isnan(
h->ComputeIntegral(
true))) {
3603 Error(
"FillRandom",
"Histograms contains negative bins, does not represent probabilities");
3612 if (ntimes > 10*nbins) {
3616 if (sumw == 0)
return;
3619 Double_t mean =
h->RetrieveBinContent(bin)*ntimes/sumw;
3630 if (sumgen < ntimes) {
3632 for (i =
Int_t(sumgen+0.5); i < ntimes; ++i)
3638 else if (sumgen > ntimes) {
3640 i =
Int_t(sumgen+0.5);
3641 while( i > ntimes) {
3656 catch(std::exception&) {}
3660 if (
h->ComputeIntegral() ==0)
return;
3663 for (loop=0;loop<ntimes;loop++) {
3689 return binx + nx*biny;
3697 return binx + nx*(biny +ny*binz);
3722 return binx + nx*biny;
3730 return binx + nx*(biny +ny*binz);
3748 Warning(
"FindFirstBinAbove",
"Invalid axis number : %d, axis x assumed\n",axis);
3762 for (
Int_t binx = firstBin; binx <= lastBin; binx++) {
3763 for (
Int_t biny = 1; biny <= nbinsy; biny++) {
3764 for (
Int_t binz = 1; binz <= nbinsz; binz++) {
3770 else if (axis == 2) {
3774 for (
Int_t biny = firstBin; biny <= lastBin; biny++) {
3775 for (
Int_t binx = 1; binx <= nbinsx; binx++) {
3776 for (
Int_t binz = 1; binz <= nbinsz; binz++) {
3782 else if (axis == 3) {
3786 for (
Int_t binz = firstBin; binz <= lastBin; binz++) {
3787 for (
Int_t binx = 1; binx <= nbinsx; binx++) {
3788 for (
Int_t biny = 1; biny <= nbinsy; biny++) {
3812 Warning(
"FindFirstBinAbove",
"Invalid axis number : %d, axis x assumed\n",axis);
3826 for (
Int_t binx = lastBin; binx >= firstBin; binx--) {
3827 for (
Int_t biny = 1; biny <= nbinsy; biny++) {
3828 for (
Int_t binz = 1; binz <= nbinsz; binz++) {
3834 else if (axis == 2) {
3838 for (
Int_t biny = lastBin; biny >= firstBin; biny--) {
3839 for (
Int_t binx = 1; binx <= nbinsx; binx++) {
3840 for (
Int_t binz = 1; binz <= nbinsz; binz++) {
3846 else if (axis == 3) {
3850 for (
Int_t binz = lastBin; binz >= firstBin; binz--) {
3851 for (
Int_t binx = 1; binx <= nbinsx; binx++) {
3852 for (
Int_t biny = 1; biny <= nbinsy; biny++) {
3909 linear= (
char*)strstr(fname,
"++");
3913 TF1 f1(fname, fname, xxmin, xxmax);
3917 TF2 f2(fname, fname);
3918 return Fit(&f2,
option,goption,xxmin,xxmax);
3921 TF3 f3(fname, fname);
3922 return Fit(&f3,
option,goption,xxmin,xxmax);
3927 if (!
f1) {
Printf(
"Unknown function: %s",fname);
return -1; }
4295 gROOT->MakeDefCanvas();
4298 Error(
"FitPanel",
"Unable to create a default canvas");
4305 if (handler && handler->
LoadPlugin() != -1) {
4307 Error(
"FitPanel",
"Unable to create the FitPanel");
4310 Error(
"FitPanel",
"Unable to find the FitPanel plug-in");
4367 asym->
Divide(top,bottom);
4381 for(
Int_t k=1; k<= zmax; k++){
4472 Info(
"SetHighlight",
"Supported only 1-D or 2-D histograms");
4489 return ((
TH1*)
this)->GetPainter()->GetObjectInfo(px,py);
4590 Error(
"GetQuantiles",
"Only available for 1-d histograms");
4600 Int_t nq = nprobSum;
4605 for (i=1;i<nq;i++) {
4610 for (i = 0; i < nq; i++) {
4612 while (ibin < nbins-1 &&
fIntegral[ibin+1] == prob[i]) {
4613 if (
fIntegral[ibin+2] == prob[i]) ibin++;
4621 if (!probSum)
delete [] prob;
4639 Double_t allcha, sumx, sumx2,
x, val, stddev, mean;
4650 allcha = sumx = sumx2 = 0;
4651 for (bin=hxfirst;bin<=hxlast;bin++) {
4654 if (val > valmax) valmax = val;
4659 if (allcha == 0)
return;
4661 stddev = sumx2/allcha - mean*mean;
4664 if (stddev == 0) stddev = binwidx*(hxlast-hxfirst+1)/4;
4671 Double_t constant = 0.5*(valmax+binwidx*allcha/(sqrtpi*stddev));
4679 if ((mean < xmin || mean >
xmax) && stddev > (
xmax-
xmin)) {
4700 Int_t nchanx = hxlast - hxfirst + 1;
4721 Int_t nchanx = hxlast - hxfirst + 1;
4724 if (nchanx <=1 || npar == 1) {
4747 const Int_t idim = 20;
4758 if (
m > idim ||
m >
n)
return;
4761 for (
l = 2;
l <=
m; ++
l) {
4763 b[
m +
l*20 - 21] = zero;
4770 for (k = hxfirst; k <= hxlast; ++k) {
4775 for (
l = 2;
l <=
m; ++
l) {
4778 da[
l-1] += power*yk;
4780 for (
l = 2;
l <=
m; ++
l) {
4782 b[
m +
l*20 - 21] += power;
4785 for (i = 3; i <=
m; ++i) {
4786 for (k = i; k <=
m; ++k) {
4787 b[k - 1 + (i-1)*20 - 21] =
b[k + (i-2)*20 - 21];
4792 for (i=0; i<
m; ++i)
a[i] = da[i];
4812 xbar = ybar = x2bar = xybar = 0;
4817 for (i = hxfirst; i <= hxlast; ++i) {
4821 if (yk <= 0) yk = 1
e-9;
4830 det = fn*x2bar - xbar*xbar;
4838 a0 = (x2bar*ybar - xbar*xybar) / det;
4839 a1 = (fn*xybar - xbar*ybar) / det;
4850 Int_t a_dim1, a_offset, b_dim1, b_offset;
4852 Int_t im1, jp1, nm1, nmi;
4858 b_offset = b_dim1 + 1;
4861 a_offset = a_dim1 + 1;
4864 if (idim <
n)
return;
4867 for (j = 1; j <=
n; ++j) {
4868 if (
a[j + j*a_dim1] <= 0) { ifail = -1;
return; }
4869 a[j + j*a_dim1] = one /
a[j + j*a_dim1];
4870 if (j ==
n)
continue;
4872 for (
l = jp1;
l <=
n; ++
l) {
4873 a[j +
l*a_dim1] =
a[j + j*a_dim1] *
a[
l + j*a_dim1];
4874 s1 = -
a[
l + (j+1)*a_dim1];
4875 for (i = 1; i <= j; ++i) {
s1 =
a[
l + i*a_dim1] *
a[i + (j+1)*a_dim1] +
s1; }
4876 a[
l + (j+1)*a_dim1] = -
s1;
4881 for (
l = 1;
l <= k; ++
l) {
4882 b[
l*b_dim1 + 1] =
a[a_dim1 + 1]*
b[
l*b_dim1 + 1];
4885 for (
l = 1;
l <= k; ++
l) {
4886 for (i = 2; i <=
n; ++i) {
4888 s21 = -
b[i +
l*b_dim1];
4889 for (j = 1; j <= im1; ++j) {
4890 s21 =
a[i + j*a_dim1]*
b[j +
l*b_dim1] + s21;
4892 b[i +
l*b_dim1] = -
a[i + i*a_dim1]*s21;
4895 for (i = 1; i <= nm1; ++i) {
4897 s22 = -
b[nmi +
l*b_dim1];
4898 for (j = 1; j <= i; ++j) {
4900 s22 =
a[nmi + nmjp1*a_dim1]*
b[nmjp1 +
l*b_dim1] + s22;
4902 b[nmi +
l*b_dim1] = -s22;
4940 if (binx < 0) binx = 0;
4941 if (binx > ofx) binx = ofx;
4956 binx = binglobal%nx;
4962 binx = binglobal%nx;
4963 biny = ((binglobal-binx)/nx)%ny;
4968 binx = binglobal%nx;
4969 biny = ((binglobal-binx)/nx)%ny;
4970 binz = ((binglobal-binx)/nx -biny)/ny;
4989 Error(
"GetRandom",
"Function only valid for 1-d histograms");
4999 integral = ((
TH1*)
this)->ComputeIntegral(
true);
5001 if (integral == 0)
return 0;
5040 if (bin < 0) bin = 0;
5066 Error(
"GetBinWithContent",
"function is only valid for 1-D histograms");
5072 if (firstx <= 0) firstx = 1;
5076 for (
Int_t i=firstx;i<=lastx;i++) {
5078 if (diff <= 0) {binx = i;
return diff;}
5079 if (diff < curmax && diff <= maxdiff) {curmax = diff, binminx=i;}
5114 return y0 + (
x-x0)*((
y1-y0)/(
x1-x0));
5123 Error(
"Interpolate",
"This function must be called with 1 argument for a TH1");
5132 Error(
"Interpolate",
"This function must be called with 1 argument for a TH1");