78 :
TH1(
name,title,nbinsx,xlow,xup),
83 Warning(
"TH3",
"nbinsy is <=0 - set to nbinsy = 1");
87 Warning(
"TH3",
"nbinsz is <=0 - set to nbinsz = 1");
90 fYaxis.Set(nbinsy,ylow,yup);
91 fZaxis.Set(nbinsz,zlow,zup);
92 fNcells = (nbinsx+2)*(nbinsy+2)*(nbinsz+2);
93 fTsumwy = fTsumwy2 = fTsumwxy = 0;
94 fTsumwz = fTsumwz2 = fTsumwxz = fTsumwyz = 0;
108 if (nbinsy <= 0) {
Warning(
"TH3",
"nbinsy is <=0 - set to nbinsy = 1"); nbinsy = 1; }
109 if (nbinsz <= 0) nbinsz = 1;
114 fNcells = (nbinsx+2)*(nbinsy+2)*(nbinsz+2);
130 if (nbinsy <= 0) {
Warning(
"TH3",
"nbinsy is <=0 - set to nbinsy = 1"); nbinsy = 1; }
131 if (nbinsz <= 0) nbinsz = 1;
136 fNcells = (nbinsx+2)*(nbinsy+2)*(nbinsz+2);
148 ((
TH3&)
h).Copy(*
this);
189 if (!nbentries)
return 0;
192 if (action == 0)
return 0;
193 nbentries = -nbentries;
208 for (
Int_t i=1;i<nbentries;i++) {
216 if (z < zmin) zmin = z;
217 if (z > zmax) zmax = z;
236 for (
Int_t i=0;i<nbentries;i++) {
237 Fill(buffer[4*i+2],buffer[4*i+3],buffer[4*i+4],buffer[4*i+1]);
263 nbentries = -nbentries;
289 Error(
"Fill",
"Invalid signature - do nothing");
304 Int_t binx, biny, binz, bin;
309 if (binx <0 || biny <0 || binz<0)
return -1;
352 Int_t binx, biny, binz, bin;
357 if (binx <0 || biny <0 || binz<0)
return -1;
397 Int_t binx, biny, binz, bin;
402 if (binx <0 || biny <0 || binz<0)
return -1;
445 Int_t binx, biny, binz, bin;
450 if (binx <0 || biny <0 || binz<0)
return -1;
493 Int_t binx, biny, binz, bin;
498 if (binx <0 || biny <0 || binz<0)
return -1;
541 Int_t binx, biny, binz, bin;
546 if (binx <0 || biny <0 || binz<0)
return -1;
587 Int_t binx, biny, binz, bin;
592 if (binx < 0 || biny < 0 || binz < 0)
639 Int_t binx, biny, binz, bin;
644 if (binx <0 || biny <0 || binz<0)
return -1;
685 Int_t binx, biny, binz, bin;
690 if (binx <0 || biny <0 || binz<0)
return -1;
744 Int_t bin, binx, biny, binz, ibin, loop;
748 if (!fobj) {
Error(
"FillRandom",
"Unknown function: %s",fname);
return; }
749 TF3 *
f1 =
dynamic_cast<TF3*
>( fobj );
750 if (!
f1) {
Error(
"FillRandom",
"Function: %s is not a TF3, is a %s",fname,fobj->IsA()->
GetName());
return; }
760 Info(
"FillRandom",
"Using function axis and range ([%g,%g],[%g,%g],[%g,%g])",
xmin,
xmax,
ymin,
ymax,zmin,zmax);
770 Int_t nxy = nbinsx*nbinsy;
771 Int_t nbins = nbinsx*nbinsy*nbinsz;
777 for (binz=1;binz<=nbinsz;binz++) {
779 for (biny=1;biny<=nbinsy;biny++) {
781 for (binx=1;binx<=nbinsx;binx++) {
789 integral[ibin] = integral[ibin-1] + fint;
795 if (integral[nbins] == 0 ) {
797 Error(
"FillRandom",
"Integral = zero");
return;
799 for (bin=1;bin<=nbins;bin++) integral[bin] /= integral[nbins];
804 for (loop=0;loop<ntimes;loop++) {
808 biny = (ibin - nxy*binz)/nbinsx;
809 binx = 1 + ibin - nbinsx*(biny + nbinsy*binz);
839 if (!
h) {
Error(
"FillRandom",
"Null histogram");
return; }
841 Error(
"FillRandom",
"Histograms with different dimensions");
return;
844 if (
h->ComputeIntegral() == 0)
return;
849 for (loop=0;loop<ntimes;loop++) {
894 auto computeFirstAndLastBin = [](
const TAxis & outerAxis,
Int_t &firstbin,
Int_t &lastbin) {
902 if (firstbin == 0 && lastbin == 0) {
907 if (firstbin < 0) firstbin = 0;
908 if (lastbin < 0 || lastbin > nbins + 1) lastbin = nbins + 1;
909 if (lastbin < firstbin) {firstbin = 0; lastbin = nbins + 1;}
912 computeFirstAndLastBin(
fXaxis, binminx, binmaxx);
913 computeFirstAndLastBin(
fYaxis, binminy, binmaxy);
916 auto computeAxisLimits = [](
const TAxis & outerAxis,
Int_t firstbin,
Int_t lastbin,
918 Int_t firstOutBin = std::max(firstbin,1);
920 nBins = lastOutBin-firstOutBin+1;
928 Int_t firstBinXaxis = computeAxisLimits(
fXaxis, binminx, binmaxx, nbinsX, xMin, xMax);
931 Int_t firstBinYaxis = computeAxisLimits(
fYaxis, binminy, binmaxy, nbinsY, yMin, yMax);
948 std::vector<TH1*> hlist(npar+1);
951 for (ipar=0;ipar<= npar;ipar++) {
961 if (
xbins->fN == 0 && ybins->
fN == 0) {
965 }
else if (
xbins->fN > 0 && ybins->
fN > 0 ) {
967 nbinsX, &
xbins->fArray[firstBinXaxis],
968 nbinsY, &ybins->
fArray[firstBinYaxis]);
976 TH1 * hchi2 = hlist.back();
984 for (
Int_t biny=binminy; biny<=binmaxy; biny++) {
985 for (
Int_t binx=binminx; binx<=binmaxx; binx++) {
992 Info(
"FitSlicesZ",
"Slice (%d,%d) skipped, the number of entries is zero or smaller than the given cut value, n=%f",binx,biny,
nentries);
999 hlist[0]->GetBinXYZ(bin,ibx,iby,ibz);
1005 if (npfits > npar && npfits >= cut) {
1006 for (ipar=0;ipar<npar;ipar++) {
1014 Info(
"FitSlicesZ",
"Fitted slice (%d,%d) skipped, the number of fitted points is too small, n=%d",binx,biny,npfits);
1029 if (biny < 0) biny = 0;
1030 if (biny > ofy) biny = ofy;
1033 if (binz < 0) binz = 0;
1034 if (binz > ofz) binz = ofz;
1070 Error(
"GetBinWithContent3",
"function is only valid for 3-D histograms");
1073 if (firstx <= 0) firstx = 1;
1075 if (firsty <= 0) firsty = 1;
1077 if (firstz <= 0) firstz = 1;
1079 Int_t binminx = 0, binminy=0, binminz=0;
1081 for (
Int_t k=firstz;k<=lastz;k++) {
1082 for (
Int_t j=firsty;j<=lasty;j++) {
1083 for (
Int_t i=firstx;i<=lastx;i++) {
1085 if (diff <= 0) {binx = i; biny=j; binz=k;
return diff;}
1086 if (diff < curmax && diff <= maxdiff) {curmax = diff, binminx=i; binminy=j;binminz=k;}
1102 if (axis1 < 1 || axis2 < 1 || axis1 > 3 || axis2 > 3) {
1103 Error(
"GetCorrelationFactor",
"Wrong parameters");
1106 if (axis1 == axis2)
return 1;
1108 if (stddev1 == 0)
return 0;
1110 if (stddev2 == 0)
return 0;
1120 if (axis1 < 1 || axis2 < 1 || axis1 > 3 || axis2 > 3) {
1121 Error(
"GetCovariance",
"Wrong parameters");
1138 if (sumw == 0)
return 0;
1139 if (axis1 == 1 && axis2 == 1) {
1140 return TMath::Abs(sumwx2/sumw - sumwx*sumwx/sumw2);
1142 if (axis1 == 2 && axis2 == 2) {
1143 return TMath::Abs(sumwy2/sumw - sumwy*sumwy/sumw2);
1145 if (axis1 == 3 && axis2 == 3) {
1146 return TMath::Abs(sumwz2/sumw - sumwz*sumwz/sumw2);
1148 if ((axis1 == 1 && axis2 == 2) || (axis1 == 2 && axis2 == 1)) {
1149 return sumwxy/sumw - sumwx/sumw*sumwy/sumw;
1151 if ((axis1 == 1 && axis2 == 3) || (axis1 == 3 && axis2 == 1)) {
1152 return sumwxz/sumw - sumwx/sumw*sumwz/sumw;
1154 if ((axis1 == 2 && axis2 == 3) || (axis1 == 3 && axis2 == 2)) {
1155 return sumwyz/sumw - sumwy/sumw*sumwz/sumw;
1174 Int_t nxy = nbinsx*nbinsy;
1175 Int_t nbins = nxy*nbinsz;
1184 if (integral == 0 ) {
x = 0;
y = 0; z = 0;
return;}
1191 Int_t binz = ibin/nxy;
1192 Int_t biny = (ibin - nxy*binz)/nbinsx;
1193 Int_t binx = ibin - nbinsx*(biny + nbinsy*binz);
1221 Int_t bin, binx, biny, binz;
1225 for (bin=0;bin<11;bin++) stats[bin] = 0;
1236 if (firstBinX == 1) firstBinX = 0;
1240 if (firstBinY == 1) firstBinY = 0;
1244 if (firstBinZ == 1) firstBinZ = 0;
1254 for (binz = firstBinZ; binz <= lastBinZ; binz++) {
1256 for (biny = firstBinY; biny <= lastBinY; biny++) {
1258 for (binx = firstBinX; binx <= lastBinX; binx++) {
1259 bin =
GetBin(binx,biny,binz);
1265 stats[1] += err*err;
1319 return DoIntegral(binx1,binx2,biny1,biny2,binz1,binz2,err,option);
1335 return DoIntegral(binx1,binx2,biny1,biny2,binz1,binz2,error,option,
kTRUE);
1343 Error(
"Interpolate",
"This function must be called with 3 arguments for a TH3");
1353 Error(
"Interpolate",
"This function must be called with 3 arguments for a TH3");
1372 Int_t obx = ubx + 1;
1376 Int_t oby = uby + 1;
1380 Int_t obz = ubz + 1;
1385 if (ubx <=0 || uby <=0 || ubz <= 0 ||
1387 Error(
"Interpolate",
"Cannot interpolate outside histogram domain.");
1412 Double_t w1 = i1 * (1 - yd) + i2 * yd;
1413 Double_t w2 = j1 * (1 - yd) + j2 * yd;
1416 Double_t result = w1 * (1 - xd) + w2 * xd;
1449 if (h2 == 0)
return 0;
1465 Error(
"KolmogorovTest",
"Histograms must be 3-D\n");
1471 Error(
"KolmogorovTest",
"Number of channels in X is different, %d and %d\n",ncx1,ncx2);
1475 Error(
"KolmogorovTest",
"Number of channels in Y is different, %d and %d\n",ncy1,ncy2);
1479 Error(
"KolmogorovTest",
"Number of channels in Z is different, %d and %d\n",ncz1,ncz2);
1489 if (diff1 > difprec || diff2 > difprec) {
1490 Error(
"KolmogorovTest",
"histograms with different binning along X");
1495 if (diff1 > difprec || diff2 > difprec) {
1496 Error(
"KolmogorovTest",
"histograms with different binning along Y");
1501 if (diff1 > difprec || diff2 > difprec) {
1502 Error(
"KolmogorovTest",
"histograms with different binning along Z");
1507 Int_t ibeg = 1, jbeg = 1, kbeg = 1;
1508 Int_t iend = ncx1, jend = ncy1, kend = ncz1;
1509 if (opt.
Contains(
"U")) {ibeg = 0; jbeg = 0; kbeg = 0;}
1510 if (opt.
Contains(
"O")) {iend = ncx1+1; jend = ncy1+1; kend = ncz1+1;}
1517 for (i = ibeg; i <= iend; i++) {
1518 for (j = jbeg; j <= jend; j++) {
1519 for (k = kbeg; k <= kend; k++) {
1534 Error(
"KolmogorovTest",
"Integral is zero for h1=%s\n",
h1->
GetName());
1538 Error(
"KolmogorovTest",
"Integral is zero for h2=%s\n",h2->
GetName());
1546 esum1 = sum1 * sum1 / w1;
1551 esum2 = sum2 * sum2 / w2;
1555 if (afunc2 && afunc1) {
1556 Error(
"KolmogorovTest",
"Errors are zero for both histograms\n");
1562 int order[3] = {0,1,2};
1566 binbeg[0] = ibeg; binbeg[1] = jbeg; binbeg[2] = kbeg;
1567 binend[0] = iend; binend[1] = jend; binend[2] = kend;
1576 for (i = binbeg[order[0] ]; i <= binend[order[0] ]; i++) {
1577 for ( j = binbeg[order[1] ]; j <= binend[order[1] ]; j++) {
1578 for ( k = binbeg[order[2] ]; k <= binend[order[2] ]; k++) {
1579 ibin[ order[0] ] = i;
1580 ibin[ order[1] ] = j;
1581 ibin[ order[2] ] = k;
1582 bin =
h1->
GetBin(ibin[0],ibin[1],ibin[2]);
1589 vdfmax[icomb] = dmax;
1608 if (opt.
Contains(
"N") && !(afunc1 || afunc2 ) ) {
1612 Double_t chi2 = d12*d12/(esum1+esum2);
1615 if (prb > 0 && prb2 > 0) prb = prb*prb2*(1-
TMath::Log(prb*prb2));
1621 printf(
" Kolmo Prob h1 = %s, sum1=%g\n",
h1->
GetName(),sum1);
1622 printf(
" Kolmo Prob h2 = %s, sum2=%g\n",h2->
GetName(),sum2);
1623 printf(
" Kolmo Probabil = %f, Max Dist = %g\n",prb,dfmax);
1625 printf(
" Kolmo Probabil = %f for shape alone, =%f for normalisation alone\n",prb1,prb2);
1631 if (opt.
Contains(
"M"))
return dfmax;
1764 computeErrors =
kTRUE;
1769 originalRange =
kTRUE;
1773 TH1D *
h1 =
DoProject1D(
name, title, projAxis, &out1, &out2, computeErrors, originalRange,
true,
true);
1792 if (padsav) padsav->
cd();
1804 bool computeErrors,
bool originalRange,
1805 bool useUF,
bool useOF)
const
1814 Int_t nx = ixmax-ixmin+1;
1820 Error(
"DoProject1D",
"Histogram with name %s must be a TH1D and is a %s",
name,h1obj->
ClassName());
1827 if ( originalRange )
1829 if (bins->
fN == 0) {
1835 if (bins->
fN == 0) {
1845 if ( originalRange )
1847 if (bins->
fN == 0) {
1853 if (bins->
fN == 0) {
1883 if (out1 ==
nullptr && out2 ==
nullptr) {
1895 R__ASSERT(out1 !=
nullptr && out2 !=
nullptr);
1897 Int_t *refX = 0, *refY = 0, *refZ = 0;
1898 Int_t ixbin, out1bin, out2bin;
1914 R__ASSERT (refX != 0 && refY != 0 && refZ != 0);
1933 for (ixbin=0;ixbin<=1+projX->
GetNbins();ixbin++) {
1940 for (out1bin = out1min; out1bin <= out1max; out1bin++) {
1941 for (out2bin = out2min; out2bin <= out2max; out2bin++) {
1947 if (computeErrors) {
1965 bool resetStats =
true;
1966 double eps = 1.E-12;
1970 bool resetEntries = resetStats;
1972 resetEntries |= !useUF || !useOF;
1979 stats[2] = stats[4];
1980 stats[3] = stats[5];
1983 stats[2] = stats[7];
1984 stats[3] = stats[8];
2013 bool computeErrors,
bool originalRange,
2014 bool useUF,
bool useOF)
const
2024 Int_t nx = ixmax-ixmin+1;
2025 Int_t ny = iymax-iymin+1;
2033 Error(
"DoProject2D",
"Histogram with name %s must be a TH2D and is a %s",
name,h2obj->
ClassName());
2041 if ( originalRange ) {
2063 if ( originalRange )
2065 if (
xbins->fN == 0 && ybins->
fN == 0) {
2068 }
else if (ybins->
fN == 0) {
2071 }
else if (
xbins->fN == 0) {
2078 if (
xbins->fN == 0 && ybins->
fN == 0) {
2081 }
else if (ybins->
fN == 0) {
2083 ,nx,&
xbins->fArray[ixmin-1]);
2084 }
else if (
xbins->fN == 0) {
2128 const TAxis* out = 0;
2137 Int_t *refX = 0, *refY = 0, *refZ = 0;
2138 Int_t ixbin, iybin, outbin;
2139 if ( projX ==
GetXaxis() && projY ==
GetYaxis() ) { refX = &ixbin; refY = &iybin; refZ = &outbin; }
2140 if ( projX ==
GetYaxis() && projY ==
GetXaxis() ) { refX = &iybin; refY = &ixbin; refZ = &outbin; }
2141 if ( projX ==
GetXaxis() && projY ==
GetZaxis() ) { refX = &ixbin; refY = &outbin; refZ = &iybin; }
2142 if ( projX ==
GetZaxis() && projY ==
GetXaxis() ) { refX = &iybin; refY = &outbin; refZ = &ixbin; }
2143 if ( projX ==
GetYaxis() && projY ==
GetZaxis() ) { refX = &outbin; refY = &ixbin; refZ = &iybin; }
2144 if ( projX ==
GetZaxis() && projY ==
GetYaxis() ) { refX = &outbin; refY = &iybin; refZ = &ixbin; }
2145 R__ASSERT (refX != 0 && refY != 0 && refZ != 0);
2154 if (outmin == 0 && outmax == 0) { outmin = 1; outmax = out->
GetNbins(); }
2159 for (ixbin=0;ixbin<=1+projX->
GetNbins();ixbin++) {
2163 for (iybin=0;iybin<=1+projY->
GetNbins();iybin++) {
2171 for (outbin = outmin; outbin <= outmax; outbin++) {
2177 if (computeErrors) {
2195 bool resetStats =
true;
2196 double eps = 1.E-12;
2200 bool resetEntries = resetStats;
2202 resetEntries |= !useUF || !useOF;
2207 for (
Int_t i = 0; i <
kNstat; ++i) { oldst[i] = 0; }
2209 std::copy(oldst,oldst+
kNstat,stats);
2213 stats[4] = oldst[7];
2214 stats[5] = oldst[8];
2215 stats[6] = oldst[9];
2218 stats[2] = oldst[4];
2219 stats[3] = oldst[5];
2221 stats[4] = oldst[2];
2222 stats[5] = oldst[3];
2225 stats[4] = oldst[7];
2226 stats[5] = oldst[8];
2227 stats[6] = oldst[10];
2231 stats[2] = oldst[7];
2232 stats[3] = oldst[8];
2234 stats[4] = oldst[2];
2235 stats[5] = oldst[3];
2236 stats[6] = oldst[9];
2239 stats[4] = oldst[4];
2240 stats[5] = oldst[5];
2241 stats[6] = oldst[10];
2256 if (!computeErrors) entries =
TMath::Floor( entries + 0.5);
2319 if (opt.
Contains(
"x")) { pcase = 1; ptype =
"x"; }
2320 if (opt.
Contains(
"y")) { pcase = 2; ptype =
"y"; }
2321 if (opt.
Contains(
"z")) { pcase = 3; ptype =
"z"; }
2322 if (opt.
Contains(
"xy")) { pcase = 4; ptype =
"xy"; }
2323 if (opt.
Contains(
"yx")) { pcase = 5; ptype =
"yx"; }
2324 if (opt.
Contains(
"xz")) { pcase = 6; ptype =
"xz"; }
2325 if (opt.
Contains(
"zx")) { pcase = 7; ptype =
"zx"; }
2326 if (opt.
Contains(
"yz")) { pcase = 8; ptype =
"yz"; }
2327 if (opt.
Contains(
"zy")) { pcase = 9; ptype =
"zy"; }
2330 Error(
"Project3D",
"No projection axis specified - return a NULL pointer");
2337 computeErrors =
kTRUE;
2354 originalRange =
kTRUE;
2365 title +=
" "; title += ptype; title +=
" projection";
2371 computeErrors, originalRange, useUF, useOF);
2377 computeErrors, originalRange, useUF, useOF);
2383 computeErrors, originalRange, useUF, useOF);
2389 computeErrors, originalRange, useUF, useOF);
2395 computeErrors, originalRange, useUF, useOF);
2401 computeErrors, originalRange, useUF, useOF);
2407 computeErrors, originalRange, useUF, useOF);
2413 computeErrors, originalRange, useUF, useOF);
2419 computeErrors, originalRange, useUF, useOF);
2435 if (padsav) padsav->
cd();
2453 if (useWeights && binSumw2.
fN <= 0) useWeights =
false;
2460 if (outBin <0)
return;
2462 if ( useWeights ) tmp = binSumw2.
fArray[outBin];
2463 p2->
Fill( u ,
v, w, cont);
2473 bool originalRange,
bool useUF,
bool useOF)
const
2481 Int_t nx = ixmax-ixmin+1;
2482 Int_t ny = iymax-iymin+1;
2492 Error(
"DoProjectProfile2D",
"Histogram with name %s must be a TProfile2D and is a %s",
name,p2obj->
ClassName());
2500 if ( originalRange ) {
2521 if ( originalRange ) {
2522 if (
xbins->fN == 0 && ybins->
fN == 0) {
2525 }
else if (ybins->
fN == 0) {
2528 }
else if (
xbins->fN == 0) {
2535 if (
xbins->fN == 0 && ybins->
fN == 0) {
2538 }
else if (ybins->
fN == 0) {
2540 ,nx,&
xbins->fArray[ixmin-1]);
2541 }
else if (
xbins->fN == 0) {
2551 const TAxis* outAxis = 0;
2566 Int_t *refX = 0, *refY = 0, *refZ = 0;
2567 Int_t ixbin, iybin, outbin;
2568 if ( projX ==
GetXaxis() && projY ==
GetYaxis() ) { refX = &ixbin; refY = &iybin; refZ = &outbin; }
2569 if ( projX ==
GetYaxis() && projY ==
GetXaxis() ) { refX = &iybin; refY = &ixbin; refZ = &outbin; }
2570 if ( projX ==
GetXaxis() && projY ==
GetZaxis() ) { refX = &ixbin; refY = &outbin; refZ = &iybin; }
2571 if ( projX ==
GetZaxis() && projY ==
GetXaxis() ) { refX = &iybin; refY = &outbin; refZ = &ixbin; }
2572 if ( projX ==
GetYaxis() && projY ==
GetZaxis() ) { refX = &outbin; refY = &ixbin; refZ = &iybin; }
2573 if ( projX ==
GetZaxis() && projY ==
GetYaxis() ) { refX = &outbin; refY = &iybin; refZ = &ixbin; }
2574 R__ASSERT (refX != 0 && refY != 0 && refZ != 0);
2584 if (useWeights && binSumw2.
fN <= 0) useWeights =
false;
2588 for (ixbin=0;ixbin<=1+projX->
GetNbins();ixbin++) {
2590 for ( iybin=0;iybin<=1+projY->
GetNbins();iybin++) {
2595 if (poutBin <0)
continue;
2597 for (outbin = outmin; outbin <= outmax; outbin++) {
2604 if (!cont)
continue;
2608 if ( useWeights ) tmp = binSumw2.
fArray[poutBin];
2618 bool resetStats =
true;
2629 if (!useWeights) entries =
TMath::Floor( entries + 0.5);
2683 if (opt.
Contains(
"xy")) { pcase = 4; ptype =
"xy"; }
2684 if (opt.
Contains(
"yx")) { pcase = 5; ptype =
"yx"; }
2685 if (opt.
Contains(
"xz")) { pcase = 6; ptype =
"xz"; }
2686 if (opt.
Contains(
"zx")) { pcase = 7; ptype =
"zx"; }
2687 if (opt.
Contains(
"yz")) { pcase = 8; ptype =
"yz"; }
2688 if (opt.
Contains(
"zy")) { pcase = 9; ptype =
"zy"; }
2691 Error(
"Project3D",
"No projection axis specified - return a NULL pointer");
2709 originalRange =
kTRUE;
2718 title +=
" profile "; title += ptype; title +=
" projection";
2780 return Rebin3D(ngroup, 1, 1, newname);
2790 return Rebin3D(1, ngroup, 1, newname);
2800 return Rebin3D(1, 1, ngroup, newname);
2831 Int_t i,j,k,xbin,ybin,zbin;
2841 if ((nxgroup <= 0) || (nxgroup > nxbins)) {
2842 Error(
"Rebin",
"Illegal value of nxgroup=%d",nxgroup);
2845 if ((nygroup <= 0) || (nygroup > nybins)) {
2846 Error(
"Rebin",
"Illegal value of nygroup=%d",nygroup);
2849 if ((nzgroup <= 0) || (nzgroup > nzbins)) {
2850 Error(
"Rebin",
"Illegal value of nzgroup=%d",nzgroup);
2854 Int_t newxbins = nxbins/nxgroup;
2855 Int_t newybins = nybins/nygroup;
2856 Int_t newzbins = nzbins/nzgroup;
2874 if (newname && strlen(newname)) {
2882 bool resetStat =
false;
2886 if (newxbins*nxgroup != nxbins) {
2890 if (newybins*nygroup != nybins) {
2894 if (newzbins*nzgroup != nzbins) {
2936 if (nxgroup != 1 || nygroup != 1 || nzgroup != 1) {
2945 hnew->
SetBins(newxbins,
xbins, newybins, ybins, newzbins, zbins);
2958 for (xbin = 1; xbin <= newxbins; xbin++) {
2961 for (ybin = 1; ybin <= newybins; ybin++) {
2963 for (zbin = 1; zbin <= newzbins; zbin++) {
2966 for (i = 0; i < nxgroup; i++) {
2967 if (oldxbin+i > nxbins)
break;
2968 for (j =0; j < nygroup; j++) {
2969 if (oldybin+j > nybins)
break;
2970 for (k =0; k < nzgroup; k++) {
2971 if (oldzbin+k > nzbins)
break;
2973 bin = oldxbin + i + (oldybin + j)*(nxbins + 2) + (oldzbin + k)*(nxbins + 2)*(nybins + 2);
2974 binContent += oldBins[bin];
2975 if (oldSumw2) binSumw2 += oldSumw2[bin];
2990 for (
Int_t xover = 0; xover <= 1; xover++) {
2991 for (
Int_t yover = 0; yover <= 1; yover++) {
2992 for (
Int_t zover = 0; zover <= 1; zover++) {
2996 for (xbin = xover*oldxbin; xbin <= xover*(nxbins+1); xbin++) {
2997 for (ybin = yover*oldybin; ybin <= yover*(nybins+1); ybin++) {
2998 for (zbin = zover*oldzbin; zbin <= zover*(nzbins+1); zbin++) {
2999 bin =
GetBin(xbin,ybin,zbin);
3000 binContent += oldBins[bin];
3001 if (oldSumw2) binSumw2 += oldSumw2[bin];
3006 yover*(newybins+1), zover*(newzbins+1) );
3013 Double_t binContent0, binContent2, binContent3, binContent4;
3014 Double_t binError0, binError2, binError3, binError4;
3015 Int_t oldxbin2, oldybin2, oldzbin2;
3016 Int_t ufbin, ofbin, ofbin2, ofbin3, ofbin4;
3022 for (xbin = 1; xbin<=newxbins; xbin++) {
3024 for (zbin = 1; zbin<=newzbins; zbin++) {
3025 binContent0 = binContent2 = 0;
3026 binError0 = binError2 = 0;
3027 for (i=0; i<nxgroup; i++) {
3028 if (oldxbin2+i > nxbins)
break;
3029 for (k=0; k<nzgroup; k++) {
3030 if (oldzbin2+k > nzbins)
break;
3032 ufbin = oldxbin2 + i + (nxbins+2)*(nybins+2)*(oldzbin2+k);
3033 binContent0 += oldBins[ufbin];
3034 if (oldSumw2) binError0 += oldSumw2[ufbin];
3035 for (ybin = oldybin; ybin <= nybins + 1; ybin++) {
3037 ofbin = ufbin + ybin*(nxbins+2);
3038 binContent2 += oldBins[ofbin];
3039 if (oldSumw2) binError2 += oldSumw2[ofbin];
3049 oldzbin2 += nzgroup;
3051 oldxbin2 += nxgroup;
3058 for (ybin = 1; ybin<=newybins; ybin++) {
3060 for (zbin = 1; zbin<=newzbins; zbin++) {
3061 binContent0 = binContent2 = 0;
3062 binError0 = binError2 = 0;
3063 for (j=0; j<nygroup; j++) {
3064 if (oldybin2+j > nybins)
break;
3065 for (k=0; k<nzgroup; k++) {
3066 if (oldzbin2+k > nzbins)
break;
3068 ufbin = (oldybin2 + j)*(nxbins+2) + (nxbins+2)*(nybins+2)*(oldzbin2+k);
3069 binContent0 += oldBins[ufbin];
3070 if (oldSumw2) binError0 += oldSumw2[ufbin];
3071 for (xbin = oldxbin; xbin <= nxbins + 1; xbin++) {
3073 ofbin = ufbin + xbin;
3074 binContent2 += oldBins[ofbin];
3075 if (oldSumw2) binError2 += oldSumw2[ofbin];
3085 oldzbin2 += nzgroup;
3087 oldybin2 += nygroup;
3094 for (xbin = 1; xbin<=newxbins; xbin++) {
3096 for (ybin = 1; ybin<=newybins; ybin++) {
3097 binContent0 = binContent2 = 0;
3098 binError0 = binError2 = 0;
3099 for (i=0; i<nxgroup; i++) {
3100 if (oldxbin2+i > nxbins)
break;
3101 for (j=0; j<nygroup; j++) {
3102 if (oldybin2+j > nybins)
break;
3104 ufbin = oldxbin2 + i + (nxbins+2)*(oldybin2+j);
3105 binContent0 += oldBins[ufbin];
3106 if (oldSumw2) binError0 += oldSumw2[ufbin];
3107 for (zbin = oldzbin; zbin <= nzbins + 1; zbin++) {
3109 ofbin = ufbin + (nxbins+2)*(nybins+2)*zbin;
3110 binContent2 += oldBins[ofbin];
3111 if (oldSumw2) binError2 += oldSumw2[ofbin];
3121 oldybin2 += nygroup;
3123 oldxbin2 += nxgroup;
3130 for (xbin = 1; xbin<=newxbins; xbin++) {
3139 for (i=0; i<nxgroup; i++) {
3140 if (oldxbin2+i > nxbins)
break;
3141 ufbin = oldxbin2 + i;
3142 binContent0 += oldBins[ufbin];
3143 if (oldSumw2) binError0 += oldSumw2[ufbin];
3144 for (ybin = oldybin; ybin <= nybins + 1; ybin++) {
3145 ofbin3 = ufbin+ybin*(nxbins+2);
3146 binContent3 += oldBins[ ofbin3 ];
3147 if (oldSumw2) binError3 += oldSumw2[ofbin3];
3148 for (zbin = oldzbin; zbin <= nzbins + 1; zbin++) {
3150 ofbin4 = oldxbin2 + i + ybin*(nxbins+2) + (nxbins+2)*(nybins+2)*zbin;
3151 binContent4 += oldBins[ofbin4];
3152 if (oldSumw2) binError4 += oldSumw2[ofbin4];
3155 for (zbin = oldzbin; zbin <= nzbins + 1; zbin++) {
3156 ofbin2 = ufbin+zbin*(nxbins+2)*(nybins+2);
3157 binContent2 += oldBins[ ofbin2 ];
3158 if (oldSumw2) binError2 += oldSumw2[ofbin2];
3164 hnew->
SetBinContent(xbin,newybins+1,newzbins+1,binContent4);
3171 oldxbin2 += nxgroup;
3178 for (zbin = 1; zbin<=newzbins; zbin++) {
3187 for (i=0; i<nzgroup; i++) {
3188 if (oldzbin2+i > nzbins)
break;
3189 ufbin = (oldzbin2 + i)*(nxbins+2)*(nybins+2);
3190 binContent0 += oldBins[ufbin];
3191 if (oldSumw2) binError0 += oldSumw2[ufbin];
3192 for (ybin = oldybin; ybin <= nybins + 1; ybin++) {
3193 ofbin3 = ufbin+ybin*(nxbins+2);
3194 binContent3 += oldBins[ ofbin3 ];
3195 if (oldSumw2) binError3 += oldSumw2[ofbin3];
3196 for (xbin = oldxbin; xbin <= nxbins + 1; xbin++) {
3198 ofbin4 = ufbin + xbin + ybin*(nxbins+2);
3199 binContent4 += oldBins[ofbin4];
3200 if (oldSumw2) binError4 += oldSumw2[ofbin4];
3203 for (xbin = oldxbin; xbin <= nxbins + 1; xbin++) {
3204 ofbin2 = xbin +(oldzbin2+i)*(nxbins+2)*(nybins+2);
3205 binContent2 += oldBins[ ofbin2 ];
3206 if (oldSumw2) binError2 += oldSumw2[ofbin2];
3212 hnew->
SetBinContent(newxbins+1,newybins+1,zbin,binContent4);
3219 oldzbin2 += nzgroup;
3226 for (ybin = 1; ybin<=newybins; ybin++) {
3235 for (i=0; i<nygroup; i++) {
3236 if (oldybin2+i > nybins)
break;
3237 ufbin = (oldybin2 + i)*(nxbins+2);
3238 binContent0 += oldBins[ufbin];
3239 if (oldSumw2) binError0 += oldSumw2[ufbin];
3240 for (xbin = oldxbin; xbin <= nxbins + 1; xbin++) {
3241 ofbin3 = ufbin+xbin;
3242 binContent3 += oldBins[ ofbin3 ];
3243 if (oldSumw2) binError3 += oldSumw2[ofbin3];
3244 for (zbin = oldzbin; zbin <= nzbins + 1; zbin++) {
3246 ofbin4 = xbin + (nxbins+2)*(nybins+2)*zbin+(oldybin2+i)*(nxbins+2);
3247 binContent4 += oldBins[ofbin4];
3248 if (oldSumw2) binError4 += oldSumw2[ofbin4];
3251 for (zbin = oldzbin; zbin <= nzbins + 1; zbin++) {
3252 ofbin2 = (oldybin2+i)*(nxbins+2)+zbin*(nxbins+2)*(nybins+2);
3253 binContent2 += oldBins[ ofbin2 ];
3254 if (oldSumw2) binError2 += oldSumw2[ofbin2];
3260 hnew->
SetBinContent(newxbins+1,ybin,newzbins+1,binContent4);
3267 oldybin2 += nygroup;
3310 if (!resetStat) hnew->
PutStats(stat);
3313 if (oldSumw2)
delete [] oldSumw2;
3344 if (bin < 0)
return;
3353 void TH3::Streamer(
TBuffer &R__b)
3363 TH1::Streamer(R__b);
3364 TAtt3D::Streamer(R__b);
3406 :
TH3(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
3421 :
TH3(
name,title,nbinsx,
xbins,nbinsy,ybins,nbinsz,zbins)
3434 :
TH3(
name,title,nbinsx,
xbins,nbinsy,ybins,nbinsz,zbins)
3465 if (newval > -128 && newval < 128) {
fArray[bin] =
Char_t(newval);
return;}
3466 if (newval < -127)
fArray[bin] = -127;
3467 if (newval > 127)
fArray[bin] = 127;
3535 void TH3C::Streamer(
TBuffer &R__b)
3548 TH1::Streamer(R__b);
3549 TArrayC::Streamer(R__b);
3551 TAtt3D::Streamer(R__b);
3553 TH3::Streamer(R__b);
3554 TArrayC::Streamer(R__b);
3667 :
TH3(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
3682 :
TH3(
name,title,nbinsx,
xbins,nbinsy,ybins,nbinsz,zbins)
3695 :
TH3(
name,title,nbinsx,
xbins,nbinsy,ybins,nbinsz,zbins)
3726 if (newval > -32768 && newval < 32768) {
fArray[bin] =
Short_t(newval);
return;}
3727 if (newval < -32767)
fArray[bin] = -32767;
3728 if (newval > 32767)
fArray[bin] = 32767;
3767 void TH3S::Streamer(
TBuffer &R__b)
3780 TH1::Streamer(R__b);
3781 TArrayS::Streamer(R__b);
3783 TAtt3D::Streamer(R__b);
3785 TH3::Streamer(R__b);
3786 TArrayS::Streamer(R__b);
3899 :
TH3(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
3914 :
TH3(
name,title,nbinsx,
xbins,nbinsy,ybins,nbinsz,zbins)
3927 :
TH3(
name,title,nbinsx,
xbins,nbinsy,ybins,nbinsz,zbins)
3958 if (newval > -INT_MAX && newval < INT_MAX) {
fArray[bin] =
Int_t(newval);
return;}
3959 if (newval < -INT_MAX)
fArray[bin] = -INT_MAX;
3960 if (newval > INT_MAX)
fArray[bin] = INT_MAX;
4098 :
TH3(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
4113 :
TH3(
name,title,nbinsx,
xbins,nbinsy,ybins,nbinsz,zbins)
4126 :
TH3(
name,title,nbinsx,
xbins,nbinsy,ybins,nbinsz,zbins)
4177 void TH3F::Streamer(
TBuffer &R__b)
4190 TH1::Streamer(R__b);
4191 TArrayF::Streamer(R__b);
4193 TAtt3D::Streamer(R__b);
4195 TH3::Streamer(R__b);
4196 TArrayF::Streamer(R__b);
4309 :
TH3(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
4324 :
TH3(
name,title,nbinsx,
xbins,nbinsy,ybins,nbinsz,zbins)
4337 :
TH3(
name,title,nbinsx,
xbins,nbinsy,ybins,nbinsz,zbins)
4388 void TH3D::Streamer(
TBuffer &R__b)
4401 TH1::Streamer(R__b);
4402 TArrayD::Streamer(R__b);
4404 TAtt3D::Streamer(R__b);
4406 TH3::Streamer(R__b);
4407 TArrayD::Streamer(R__b);