94 :
TH1(
name,title,nbinsx,xlow,xup),
99 Warning(
"TH3",
"nbinsy is <=0 - set to nbinsy = 1");
103 Warning(
"TH3",
"nbinsz is <=0 - set to nbinsz = 1");
108 fNcells = (nbinsx+2)*(nbinsy+2)*(nbinsz+2);
140 if (nbinsy <= 0) {
Warning(
"TH3",
"nbinsy is <=0 - set to nbinsy = 1"); nbinsy = 1; }
141 if (nbinsz <= 0) nbinsz = 1;
146 fNcells = (nbinsx+2)*(nbinsy+2)*(nbinsz+2);
178 if (nbinsy <= 0) {
Warning(
"TH3",
"nbinsy is <=0 - set to nbinsy = 1"); nbinsy = 1; }
179 if (nbinsz <= 0) nbinsz = 1;
184 fNcells = (nbinsx+2)*(nbinsy+2)*(nbinsz+2);
197 ((
TH3&)
h).Copy(*
this);
238 if (!nbentries)
return 0;
241 if (action == 0)
return 0;
242 nbentries = -nbentries;
257 for (
Int_t i=1;i<nbentries;i++) {
265 if (z < zmin) zmin = z;
266 if (z > zmax) zmax = z;
285 for (
Int_t i=0;i<nbentries;i++) {
286 Fill(buffer[4*i+2],buffer[4*i+3],buffer[4*i+4],buffer[4*i+1]);
313 nbentries = -nbentries;
339 Error(
"Fill",
"Invalid signature - do nothing");
354 Int_t binx, biny, binz, bin;
359 if (binx <0 || biny <0 || binz<0)
return -1;
402 Int_t binx, biny, binz, bin;
407 if (binx <0 || biny <0 || binz<0)
return -1;
447 Int_t binx, biny, binz, bin;
452 if (binx <0 || biny <0 || binz<0)
return -1;
495 Int_t binx, biny, binz, bin;
500 if (binx <0 || biny <0 || binz<0)
return -1;
543 Int_t binx, biny, binz, bin;
548 if (binx <0 || biny <0 || binz<0)
return -1;
591 Int_t binx, biny, binz, bin;
596 if (binx <0 || biny <0 || binz<0)
return -1;
637 Int_t binx, biny, binz, bin;
642 if (binx < 0 || biny < 0 || binz < 0)
689 Int_t binx, biny, binz, bin;
694 if (binx <0 || biny <0 || binz<0)
return -1;
735 Int_t binx, biny, binz, bin;
740 if (binx <0 || biny <0 || binz<0)
return -1;
793 Int_t bin, binx, biny, binz, ibin, loop;
797 if (!fobj) {
Error(
"FillRandom",
"Unknown function: %s",fname);
return; }
798 TF3 *
f1 =
dynamic_cast<TF3*
>( fobj );
799 if (!
f1) {
Error(
"FillRandom",
"Function: %s is not a TF3, is a %s",fname,fobj->
IsA()->
GetName());
return; }
809 Info(
"FillRandom",
"Using function axis and range ([%g,%g],[%g,%g],[%g,%g])",
xmin,
xmax,
ymin,
ymax,zmin,zmax);
819 Int_t nxy = nbinsx*nbinsy;
820 Int_t nbins = nbinsx*nbinsy*nbinsz;
826 for (binz=1;binz<=nbinsz;binz++) {
828 for (biny=1;biny<=nbinsy;biny++) {
830 for (binx=1;binx<=nbinsx;binx++) {
838 integral[ibin] = integral[ibin-1] + fint;
844 if (integral[nbins] == 0 ) {
846 Error(
"FillRandom",
"Integral = zero");
return;
848 for (bin=1;bin<=nbins;bin++) integral[bin] /= integral[nbins];
853 for (loop=0;loop<ntimes;loop++) {
857 biny = (ibin - nxy*binz)/nbinsx;
858 binx = 1 + ibin - nbinsx*(biny + nbinsy*binz);
888 if (!
h) {
Error(
"FillRandom",
"Null histogram");
return; }
890 Error(
"FillRandom",
"Histograms with different dimensions");
return;
893 if (
h->ComputeIntegral() == 0)
return;
898 for (loop=0;loop<ntimes;loop++) {
943 auto computeFirstAndLastBin = [](
const TAxis & outerAxis,
Int_t &firstbin,
Int_t &lastbin) {
951 if (firstbin == 0 && lastbin == 0) {
956 if (firstbin < 0) firstbin = 0;
957 if (lastbin < 0 || lastbin > nbins + 1) lastbin = nbins + 1;
958 if (lastbin < firstbin) {firstbin = 0; lastbin = nbins + 1;}
961 computeFirstAndLastBin(
fXaxis, binminx, binmaxx);
962 computeFirstAndLastBin(
fYaxis, binminy, binmaxy);
965 auto computeAxisLimits = [](
const TAxis & outerAxis,
Int_t firstbin,
Int_t lastbin,
967 Int_t firstOutBin = std::max(firstbin,1);
969 nBins = lastOutBin-firstOutBin+1;
977 Int_t firstBinXaxis = computeAxisLimits(
fXaxis, binminx, binmaxx, nbinsX, xMin, xMax);
980 Int_t firstBinYaxis = computeAxisLimits(
fYaxis, binminy, binmaxy, nbinsY, yMin, yMax);
997 std::vector<TH1*> hlist(npar+1);
1000 for (ipar=0;ipar<= npar;ipar++) {
1008 title =
"chisquare";
1010 if (
xbins->fN == 0 && ybins->
fN == 0) {
1011 hlist[ipar] =
new TH2D(
name, title,
1013 nbinsY, yMin, yMax);
1014 }
else if (
xbins->fN > 0 && ybins->
fN > 0 ) {
1015 hlist[ipar] =
new TH2D(
name, title,
1016 nbinsX, &
xbins->fArray[firstBinXaxis],
1017 nbinsY, &ybins->
fArray[firstBinYaxis]);
1025 TH1 * hchi2 = hlist.back();
1028 TH1D *hpz =
nullptr;
1033 for (
Int_t biny=binminy; biny<=binmaxy; biny++) {
1034 for (
Int_t binx=binminx; binx<=binmaxx; binx++) {
1041 Info(
"FitSlicesZ",
"Slice (%d,%d) skipped, the number of entries is zero or smaller than the given cut value, n=%f",binx,biny,
nentries);
1047 int ibx,iby,ibz = 0;
1048 hlist[0]->GetBinXYZ(bin,ibx,iby,ibz);
1054 if (npfits > npar && npfits >= cut) {
1055 for (ipar=0;ipar<npar;ipar++) {
1063 Info(
"FitSlicesZ",
"Fitted slice (%d,%d) skipped, the number of fitted points is too small, n=%d",binx,biny,npfits);
1078 if (biny < 0) biny = 0;
1079 if (biny > ofy) biny = ofy;
1082 if (binz < 0) binz = 0;
1083 if (binz > ofz) binz = ofz;
1119 Error(
"GetBinWithContent3",
"function is only valid for 3-D histograms");
1122 if (firstx <= 0) firstx = 1;
1124 if (firsty <= 0) firsty = 1;
1126 if (firstz <= 0) firstz = 1;
1128 Int_t binminx = 0, binminy=0, binminz=0;
1130 for (
Int_t k=firstz;k<=lastz;k++) {
1131 for (
Int_t j=firsty;j<=lasty;j++) {
1132 for (
Int_t i=firstx;i<=lastx;i++) {
1134 if (diff <= 0) {binx = i; biny=j; binz=k;
return diff;}
1135 if (diff < curmax && diff <= maxdiff) {curmax = diff, binminx=i; binminy=j;binminz=k;}
1151 if (axis1 < 1 || axis2 < 1 || axis1 > 3 || axis2 > 3) {
1152 Error(
"GetCorrelationFactor",
"Wrong parameters");
1155 if (axis1 == axis2)
return 1;
1157 if (stddev1 == 0)
return 0;
1159 if (stddev2 == 0)
return 0;
1169 if (axis1 < 1 || axis2 < 1 || axis1 > 3 || axis2 > 3) {
1170 Error(
"GetCovariance",
"Wrong parameters");
1186 if (sumw == 0)
return 0;
1187 if (axis1 == 1 && axis2 == 1) {
1188 return TMath::Abs(sumwx2/sumw - sumwx*sumwx/(sumw*sumw));
1190 if (axis1 == 2 && axis2 == 2) {
1191 return TMath::Abs(sumwy2/sumw - sumwy*sumwy/(sumw*sumw));
1193 if (axis1 == 3 && axis2 == 3) {
1194 return TMath::Abs(sumwz2/sumw - sumwz*sumwz/(sumw*sumw));
1196 if ((axis1 == 1 && axis2 == 2) || (axis1 == 2 && axis2 == 1)) {
1197 return sumwxy/sumw - sumwx*sumwy/(sumw*sumw);
1199 if ((axis1 == 1 && axis2 == 3) || (axis1 == 3 && axis2 == 1)) {
1200 return sumwxz/sumw - sumwx*sumwz/(sumw*sumw);
1202 if ((axis1 == 2 && axis2 == 3) || (axis1 == 3 && axis2 == 2)) {
1203 return sumwyz/sumw - sumwy*sumwz/(sumw*sumw);
1222 Int_t nxy = nbinsx*nbinsy;
1223 Int_t nbins = nxy*nbinsz;
1232 if (integral == 0 ) {
x = 0;
y = 0; z = 0;
return;}
1239 Int_t binz = ibin/nxy;
1240 Int_t biny = (ibin - nxy*binz)/nbinsx;
1241 Int_t binx = ibin - nbinsx*(biny + nbinsy*binz);
1269 Int_t bin, binx, biny, binz;
1273 for (bin=0;bin<11;bin++) stats[bin] = 0;
1284 if (firstBinX == 1) firstBinX = 0;
1288 if (firstBinY == 1) firstBinY = 0;
1292 if (firstBinZ == 1) firstBinZ = 0;
1302 for (binz = firstBinZ; binz <= lastBinZ; binz++) {
1304 for (biny = firstBinY; biny <= lastBinY; biny++) {
1306 for (binx = firstBinX; binx <= lastBinX; binx++) {
1307 bin =
GetBin(binx,biny,binz);
1313 stats[1] += err*err;
1391 Error(
"Interpolate",
"This function must be called with 3 arguments for a TH3");
1401 Error(
"Interpolate",
"This function must be called with 3 arguments for a TH3");
1420 Int_t obx = ubx + 1;
1424 Int_t oby = uby + 1;
1428 Int_t obz = ubz + 1;
1433 if (ubx <=0 || uby <=0 || ubz <= 0 ||
1435 Error(
"Interpolate",
"Cannot interpolate outside histogram domain.");
1460 Double_t w1 = i1 * (1 - yd) + i2 * yd;
1461 Double_t w2 = j1 * (1 - yd) + j2 * yd;
1497 if (h2 == 0)
return 0;
1513 Error(
"KolmogorovTest",
"Histograms must be 3-D\n");
1519 Error(
"KolmogorovTest",
"Number of channels in X is different, %d and %d\n",ncx1,ncx2);
1523 Error(
"KolmogorovTest",
"Number of channels in Y is different, %d and %d\n",ncy1,ncy2);
1527 Error(
"KolmogorovTest",
"Number of channels in Z is different, %d and %d\n",ncz1,ncz2);
1537 if (diff1 > difprec || diff2 > difprec) {
1538 Error(
"KolmogorovTest",
"histograms with different binning along X");
1543 if (diff1 > difprec || diff2 > difprec) {
1544 Error(
"KolmogorovTest",
"histograms with different binning along Y");
1549 if (diff1 > difprec || diff2 > difprec) {
1550 Error(
"KolmogorovTest",
"histograms with different binning along Z");
1555 Int_t ibeg = 1, jbeg = 1, kbeg = 1;
1556 Int_t iend = ncx1, jend = ncy1, kend = ncz1;
1557 if (opt.
Contains(
"U")) {ibeg = 0; jbeg = 0; kbeg = 0;}
1558 if (opt.
Contains(
"O")) {iend = ncx1+1; jend = ncy1+1; kend = ncz1+1;}
1565 for (i = ibeg; i <= iend; i++) {
1566 for (j = jbeg; j <= jend; j++) {
1567 for (k = kbeg; k <= kend; k++) {
1582 Error(
"KolmogorovTest",
"Integral is zero for h1=%s\n",
h1->
GetName());
1586 Error(
"KolmogorovTest",
"Integral is zero for h2=%s\n",h2->
GetName());
1594 esum1 = sum1 * sum1 / w1;
1599 esum2 = sum2 * sum2 / w2;
1603 if (afunc2 && afunc1) {
1604 Error(
"KolmogorovTest",
"Errors are zero for both histograms\n");
1610 int order[3] = {0,1,2};
1614 binbeg[0] = ibeg; binbeg[1] = jbeg; binbeg[2] = kbeg;
1615 binend[0] = iend; binend[1] = jend; binend[2] = kend;
1624 for (i = binbeg[order[0] ]; i <= binend[order[0] ]; i++) {
1625 for ( j = binbeg[order[1] ]; j <= binend[order[1] ]; j++) {
1626 for ( k = binbeg[order[2] ]; k <= binend[order[2] ]; k++) {
1627 ibin[ order[0] ] = i;
1628 ibin[ order[1] ] = j;
1629 ibin[ order[2] ] = k;
1630 bin =
h1->
GetBin(ibin[0],ibin[1],ibin[2]);
1637 vdfmax[icomb] = dmax;
1656 if (opt.
Contains(
"N") && !(afunc1 || afunc2 ) ) {
1660 Double_t chi2 = d12*d12/(esum1+esum2);
1663 if (prb > 0 && prb2 > 0) prb = prb*prb2*(1-
TMath::Log(prb*prb2));
1669 printf(
" Kolmo Prob h1 = %s, sum1=%g\n",
h1->
GetName(),sum1);
1670 printf(
" Kolmo Prob h2 = %s, sum2=%g\n",h2->
GetName(),sum2);
1671 printf(
" Kolmo Probabil = %f, Max Dist = %g\n",prb,dfmax);
1673 printf(
" Kolmo Probabil = %f for shape alone, =%f for normalisation alone\n",prb1,prb2);
1679 if (opt.
Contains(
"M"))
return dfmax;
1812 computeErrors =
kTRUE;
1817 originalRange =
kTRUE;
1821 TH1D *
h1 =
DoProject1D(
name, title, projAxis, &out1, &out2, computeErrors, originalRange,
true,
true);
1840 if (padsav) padsav->
cd();
1852 bool computeErrors,
bool originalRange,
1853 bool useUF,
bool useOF)
const
1862 Int_t nx = ixmax-ixmin+1;
1868 Error(
"DoProject1D",
"Histogram with name %s must be a TH1D and is a %s",
name,h1obj->
ClassName());
1875 if ( originalRange )
1877 if (bins->
fN == 0) {
1883 if (bins->
fN == 0) {
1893 if ( originalRange )
1895 if (bins->
fN == 0) {
1901 if (bins->
fN == 0) {
1931 if (out1 ==
nullptr && out2 ==
nullptr) {
1943 R__ASSERT(out1 !=
nullptr && out2 !=
nullptr);
1945 Int_t *refX = 0, *refY = 0, *refZ = 0;
1946 Int_t ixbin, out1bin, out2bin;
1962 R__ASSERT (refX != 0 && refY != 0 && refZ != 0);
1981 for (ixbin=0;ixbin<=1+projX->
GetNbins();ixbin++) {
1988 for (out1bin = out1min; out1bin <= out1max; out1bin++) {
1989 for (out2bin = out2min; out2bin <= out2max; out2bin++) {
1995 if (computeErrors) {
2013 bool resetStats =
true;
2014 double eps = 1.E-12;
2018 bool resetEntries = resetStats;
2020 resetEntries |= !useUF || !useOF;
2027 stats[2] = stats[4];
2028 stats[3] = stats[5];
2031 stats[2] = stats[7];
2032 stats[3] = stats[8];
2061 bool computeErrors,
bool originalRange,
2062 bool useUF,
bool useOF)
const
2072 Int_t nx = ixmax-ixmin+1;
2073 Int_t ny = iymax-iymin+1;
2081 Error(
"DoProject2D",
"Histogram with name %s must be a TH2D and is a %s",
name,h2obj->
ClassName());
2089 if ( originalRange ) {
2111 if ( originalRange )
2113 if (
xbins->fN == 0 && ybins->
fN == 0) {
2116 }
else if (ybins->
fN == 0) {
2119 }
else if (
xbins->fN == 0) {
2126 if (
xbins->fN == 0 && ybins->
fN == 0) {
2129 }
else if (ybins->
fN == 0) {
2131 ,nx,&
xbins->fArray[ixmin-1]);
2132 }
else if (
xbins->fN == 0) {
2176 const TAxis* out = 0;
2185 Int_t *refX = 0, *refY = 0, *refZ = 0;
2186 Int_t ixbin, iybin, outbin;
2187 if ( projX ==
GetXaxis() && projY ==
GetYaxis() ) { refX = &ixbin; refY = &iybin; refZ = &outbin; }
2188 if ( projX ==
GetYaxis() && projY ==
GetXaxis() ) { refX = &iybin; refY = &ixbin; refZ = &outbin; }
2189 if ( projX ==
GetXaxis() && projY ==
GetZaxis() ) { refX = &ixbin; refY = &outbin; refZ = &iybin; }
2190 if ( projX ==
GetZaxis() && projY ==
GetXaxis() ) { refX = &iybin; refY = &outbin; refZ = &ixbin; }
2191 if ( projX ==
GetYaxis() && projY ==
GetZaxis() ) { refX = &outbin; refY = &ixbin; refZ = &iybin; }
2192 if ( projX ==
GetZaxis() && projY ==
GetYaxis() ) { refX = &outbin; refY = &iybin; refZ = &ixbin; }
2193 R__ASSERT (refX != 0 && refY != 0 && refZ != 0);
2202 if (outmin == 0 && outmax == 0) { outmin = 1; outmax = out->
GetNbins(); }
2207 for (ixbin=0;ixbin<=1+projX->
GetNbins();ixbin++) {
2211 for (iybin=0;iybin<=1+projY->
GetNbins();iybin++) {
2219 for (outbin = outmin; outbin <= outmax; outbin++) {
2225 if (computeErrors) {
2243 bool resetStats =
true;
2244 double eps = 1.E-12;
2248 bool resetEntries = resetStats;
2250 resetEntries |= !useUF || !useOF;
2255 for (
Int_t i = 0; i <
kNstat; ++i) { oldst[i] = 0; }
2257 std::copy(oldst,oldst+
kNstat,stats);
2261 stats[4] = oldst[7];
2262 stats[5] = oldst[8];
2263 stats[6] = oldst[9];
2266 stats[2] = oldst[4];
2267 stats[3] = oldst[5];
2269 stats[4] = oldst[2];
2270 stats[5] = oldst[3];
2273 stats[4] = oldst[7];
2274 stats[5] = oldst[8];
2275 stats[6] = oldst[10];
2279 stats[2] = oldst[7];
2280 stats[3] = oldst[8];
2282 stats[4] = oldst[2];
2283 stats[5] = oldst[3];
2284 stats[6] = oldst[9];
2287 stats[4] = oldst[4];
2288 stats[5] = oldst[5];
2289 stats[6] = oldst[10];
2304 if (!computeErrors) entries =
TMath::Floor( entries + 0.5);
2371 if (opt.
Contains(
"x")) { pcase = 1; ptype =
"x"; }
2372 if (opt.
Contains(
"y")) { pcase = 2; ptype =
"y"; }
2373 if (opt.
Contains(
"z")) { pcase = 3; ptype =
"z"; }
2374 if (opt.
Contains(
"xy")) { pcase = 4; ptype =
"xy"; }
2375 if (opt.
Contains(
"yx")) { pcase = 5; ptype =
"yx"; }
2376 if (opt.
Contains(
"xz")) { pcase = 6; ptype =
"xz"; }
2377 if (opt.
Contains(
"zx")) { pcase = 7; ptype =
"zx"; }
2378 if (opt.
Contains(
"yz")) { pcase = 8; ptype =
"yz"; }
2379 if (opt.
Contains(
"zy")) { pcase = 9; ptype =
"zy"; }
2382 Error(
"Project3D",
"No projection axis specified - return a NULL pointer");
2389 computeErrors =
kTRUE;
2406 originalRange =
kTRUE;
2417 title +=
" "; title += ptype; title +=
" projection";
2423 computeErrors, originalRange, useUF, useOF);
2429 computeErrors, originalRange, useUF, useOF);
2435 computeErrors, originalRange, useUF, useOF);
2441 computeErrors, originalRange, useUF, useOF);
2447 computeErrors, originalRange, useUF, useOF);
2453 computeErrors, originalRange, useUF, useOF);
2459 computeErrors, originalRange, useUF, useOF);
2465 computeErrors, originalRange, useUF, useOF);
2471 computeErrors, originalRange, useUF, useOF);
2487 if (padsav) padsav->
cd();
2505 if (useWeights && binSumw2.
fN <= 0) useWeights =
false;
2512 if (outBin <0)
return;
2514 if ( useWeights ) tmp = binSumw2.
fArray[outBin];
2515 p2->
Fill( u ,
v,
w, cont);
2525 bool originalRange,
bool useUF,
bool useOF)
const
2533 Int_t nx = ixmax-ixmin+1;
2534 Int_t ny = iymax-iymin+1;
2544 Error(
"DoProjectProfile2D",
"Histogram with name %s must be a TProfile2D and is a %s",
name,p2obj->
ClassName());
2552 if ( originalRange ) {
2573 if ( originalRange ) {
2574 if (
xbins->fN == 0 && ybins->
fN == 0) {
2577 }
else if (ybins->
fN == 0) {
2580 }
else if (
xbins->fN == 0) {
2587 if (
xbins->fN == 0 && ybins->
fN == 0) {
2590 }
else if (ybins->
fN == 0) {
2592 ,nx,&
xbins->fArray[ixmin-1]);
2593 }
else if (
xbins->fN == 0) {
2603 const TAxis* outAxis = 0;
2618 Int_t *refX = 0, *refY = 0, *refZ = 0;
2619 Int_t ixbin, iybin, outbin;
2620 if ( projX ==
GetXaxis() && projY ==
GetYaxis() ) { refX = &ixbin; refY = &iybin; refZ = &outbin; }
2621 if ( projX ==
GetYaxis() && projY ==
GetXaxis() ) { refX = &iybin; refY = &ixbin; refZ = &outbin; }
2622 if ( projX ==
GetXaxis() && projY ==
GetZaxis() ) { refX = &ixbin; refY = &outbin; refZ = &iybin; }
2623 if ( projX ==
GetZaxis() && projY ==
GetXaxis() ) { refX = &iybin; refY = &outbin; refZ = &ixbin; }
2624 if ( projX ==
GetYaxis() && projY ==
GetZaxis() ) { refX = &outbin; refY = &ixbin; refZ = &iybin; }
2625 if ( projX ==
GetZaxis() && projY ==
GetYaxis() ) { refX = &outbin; refY = &iybin; refZ = &ixbin; }
2626 R__ASSERT (refX != 0 && refY != 0 && refZ != 0);
2636 if (useWeights && binSumw2.
fN <= 0) useWeights =
false;
2640 for (ixbin=0;ixbin<=1+projX->
GetNbins();ixbin++) {
2642 for ( iybin=0;iybin<=1+projY->
GetNbins();iybin++) {
2647 if (poutBin <0)
continue;
2649 for (outbin = outmin; outbin <= outmax; outbin++) {
2656 if (!cont)
continue;
2660 if ( useWeights ) tmp = binSumw2.
fArray[poutBin];
2670 bool resetStats =
true;
2681 if (!useWeights) entries =
TMath::Floor( entries + 0.5);
2735 if (opt.
Contains(
"xy")) { pcase = 4; ptype =
"xy"; }
2736 if (opt.
Contains(
"yx")) { pcase = 5; ptype =
"yx"; }
2737 if (opt.
Contains(
"xz")) { pcase = 6; ptype =
"xz"; }
2738 if (opt.
Contains(
"zx")) { pcase = 7; ptype =
"zx"; }
2739 if (opt.
Contains(
"yz")) { pcase = 8; ptype =
"yz"; }
2740 if (opt.
Contains(
"zy")) { pcase = 9; ptype =
"zy"; }
2743 Error(
"Project3D",
"No projection axis specified - return a NULL pointer");
2761 originalRange =
kTRUE;
2770 title +=
" profile "; title += ptype; title +=
" projection";
2832 return Rebin3D(ngroup, 1, 1, newname);
2842 return Rebin3D(1, ngroup, 1, newname);
2852 return Rebin3D(1, 1, ngroup, newname);
2883 Int_t i,j,k,xbin,ybin,zbin;
2893 if ((nxgroup <= 0) || (nxgroup > nxbins)) {
2894 Error(
"Rebin",
"Illegal value of nxgroup=%d",nxgroup);
2897 if ((nygroup <= 0) || (nygroup > nybins)) {
2898 Error(
"Rebin",
"Illegal value of nygroup=%d",nygroup);
2901 if ((nzgroup <= 0) || (nzgroup > nzbins)) {
2902 Error(
"Rebin",
"Illegal value of nzgroup=%d",nzgroup);
2906 Int_t newxbins = nxbins/nxgroup;
2907 Int_t newybins = nybins/nygroup;
2908 Int_t newzbins = nzbins/nzgroup;
2926 if (newname && strlen(newname)) {
2934 bool resetStat =
false;
2938 if (newxbins*nxgroup != nxbins) {
2942 if (newybins*nygroup != nybins) {
2946 if (newzbins*nzgroup != nzbins) {
2988 if (nxgroup != 1 || nygroup != 1 || nzgroup != 1) {
2997 hnew->
SetBins(newxbins,
xbins, newybins, ybins, newzbins, zbins);
3010 for (xbin = 1; xbin <= newxbins; xbin++) {
3013 for (ybin = 1; ybin <= newybins; ybin++) {
3015 for (zbin = 1; zbin <= newzbins; zbin++) {
3018 for (i = 0; i < nxgroup; i++) {
3019 if (oldxbin+i > nxbins)
break;
3020 for (j =0; j < nygroup; j++) {
3021 if (oldybin+j > nybins)
break;
3022 for (k =0; k < nzgroup; k++) {
3023 if (oldzbin+k > nzbins)
break;
3025 bin = oldxbin + i + (oldybin + j)*(nxbins + 2) + (oldzbin + k)*(nxbins + 2)*(nybins + 2);
3026 binContent += oldBins[bin];
3027 if (oldSumw2) binSumw2 += oldSumw2[bin];
3042 for (
Int_t xover = 0; xover <= 1; xover++) {
3043 for (
Int_t yover = 0; yover <= 1; yover++) {
3044 for (
Int_t zover = 0; zover <= 1; zover++) {
3048 for (xbin = xover*oldxbin; xbin <= xover*(nxbins+1); xbin++) {
3049 for (ybin = yover*oldybin; ybin <= yover*(nybins+1); ybin++) {
3050 for (zbin = zover*oldzbin; zbin <= zover*(nzbins+1); zbin++) {
3051 bin =
GetBin(xbin,ybin,zbin);
3052 binContent += oldBins[bin];
3053 if (oldSumw2) binSumw2 += oldSumw2[bin];
3058 yover*(newybins+1), zover*(newzbins+1) );
3065 Double_t binContent0, binContent2, binContent3, binContent4;
3066 Double_t binError0, binError2, binError3, binError4;
3067 Int_t oldxbin2, oldybin2, oldzbin2;
3068 Int_t ufbin, ofbin, ofbin2, ofbin3, ofbin4;
3074 for (xbin = 1; xbin<=newxbins; xbin++) {
3076 for (zbin = 1; zbin<=newzbins; zbin++) {
3077 binContent0 = binContent2 = 0;
3078 binError0 = binError2 = 0;
3079 for (i=0; i<nxgroup; i++) {
3080 if (oldxbin2+i > nxbins)
break;
3081 for (k=0; k<nzgroup; k++) {
3082 if (oldzbin2+k > nzbins)
break;
3084 ufbin = oldxbin2 + i + (nxbins+2)*(nybins+2)*(oldzbin2+k);
3085 binContent0 += oldBins[ufbin];
3086 if (oldSumw2) binError0 += oldSumw2[ufbin];
3087 for (ybin = oldybin; ybin <= nybins + 1; ybin++) {
3089 ofbin = ufbin + ybin*(nxbins+2);
3090 binContent2 += oldBins[ofbin];
3091 if (oldSumw2) binError2 += oldSumw2[ofbin];
3101 oldzbin2 += nzgroup;
3103 oldxbin2 += nxgroup;
3110 for (ybin = 1; ybin<=newybins; ybin++) {
3112 for (zbin = 1; zbin<=newzbins; zbin++) {
3113 binContent0 = binContent2 = 0;
3114 binError0 = binError2 = 0;
3115 for (j=0; j<nygroup; j++) {
3116 if (oldybin2+j > nybins)
break;
3117 for (k=0; k<nzgroup; k++) {
3118 if (oldzbin2+k > nzbins)
break;
3120 ufbin = (oldybin2 + j)*(nxbins+2) + (nxbins+2)*(nybins+2)*(oldzbin2+k);
3121 binContent0 += oldBins[ufbin];
3122 if (oldSumw2) binError0 += oldSumw2[ufbin];
3123 for (xbin = oldxbin; xbin <= nxbins + 1; xbin++) {
3125 ofbin = ufbin + xbin;
3126 binContent2 += oldBins[ofbin];
3127 if (oldSumw2) binError2 += oldSumw2[ofbin];
3137 oldzbin2 += nzgroup;
3139 oldybin2 += nygroup;
3146 for (xbin = 1; xbin<=newxbins; xbin++) {
3148 for (ybin = 1; ybin<=newybins; ybin++) {
3149 binContent0 = binContent2 = 0;
3150 binError0 = binError2 = 0;
3151 for (i=0; i<nxgroup; i++) {
3152 if (oldxbin2+i > nxbins)
break;
3153 for (j=0; j<nygroup; j++) {
3154 if (oldybin2+j > nybins)
break;
3156 ufbin = oldxbin2 + i + (nxbins+2)*(oldybin2+j);
3157 binContent0 += oldBins[ufbin];
3158 if (oldSumw2) binError0 += oldSumw2[ufbin];
3159 for (zbin = oldzbin; zbin <= nzbins + 1; zbin++) {
3161 ofbin = ufbin + (nxbins+2)*(nybins+2)*zbin;
3162 binContent2 += oldBins[ofbin];
3163 if (oldSumw2) binError2 += oldSumw2[ofbin];
3173 oldybin2 += nygroup;
3175 oldxbin2 += nxgroup;
3182 for (xbin = 1; xbin<=newxbins; xbin++) {
3191 for (i=0; i<nxgroup; i++) {
3192 if (oldxbin2+i > nxbins)
break;
3193 ufbin = oldxbin2 + i;
3194 binContent0 += oldBins[ufbin];
3195 if (oldSumw2) binError0 += oldSumw2[ufbin];
3196 for (ybin = oldybin; ybin <= nybins + 1; ybin++) {
3197 ofbin3 = ufbin+ybin*(nxbins+2);
3198 binContent3 += oldBins[ ofbin3 ];
3199 if (oldSumw2) binError3 += oldSumw2[ofbin3];
3200 for (zbin = oldzbin; zbin <= nzbins + 1; zbin++) {
3202 ofbin4 = oldxbin2 + i + ybin*(nxbins+2) + (nxbins+2)*(nybins+2)*zbin;
3203 binContent4 += oldBins[ofbin4];
3204 if (oldSumw2) binError4 += oldSumw2[ofbin4];
3207 for (zbin = oldzbin; zbin <= nzbins + 1; zbin++) {
3208 ofbin2 = ufbin+zbin*(nxbins+2)*(nybins+2);
3209 binContent2 += oldBins[ ofbin2 ];
3210 if (oldSumw2) binError2 += oldSumw2[ofbin2];
3216 hnew->
SetBinContent(xbin,newybins+1,newzbins+1,binContent4);
3223 oldxbin2 += nxgroup;
3230 for (zbin = 1; zbin<=newzbins; zbin++) {
3239 for (i=0; i<nzgroup; i++) {
3240 if (oldzbin2+i > nzbins)
break;
3241 ufbin = (oldzbin2 + i)*(nxbins+2)*(nybins+2);
3242 binContent0 += oldBins[ufbin];
3243 if (oldSumw2) binError0 += oldSumw2[ufbin];
3244 for (ybin = oldybin; ybin <= nybins + 1; ybin++) {
3245 ofbin3 = ufbin+ybin*(nxbins+2);
3246 binContent3 += oldBins[ ofbin3 ];
3247 if (oldSumw2) binError3 += oldSumw2[ofbin3];
3248 for (xbin = oldxbin; xbin <= nxbins + 1; xbin++) {
3250 ofbin4 = ufbin + xbin + ybin*(nxbins+2);
3251 binContent4 += oldBins[ofbin4];
3252 if (oldSumw2) binError4 += oldSumw2[ofbin4];
3255 for (xbin = oldxbin; xbin <= nxbins + 1; xbin++) {
3256 ofbin2 = xbin +(oldzbin2+i)*(nxbins+2)*(nybins+2);
3257 binContent2 += oldBins[ ofbin2 ];
3258 if (oldSumw2) binError2 += oldSumw2[ofbin2];
3264 hnew->
SetBinContent(newxbins+1,newybins+1,zbin,binContent4);
3271 oldzbin2 += nzgroup;
3278 for (ybin = 1; ybin<=newybins; ybin++) {
3287 for (i=0; i<nygroup; i++) {
3288 if (oldybin2+i > nybins)
break;
3289 ufbin = (oldybin2 + i)*(nxbins+2);
3290 binContent0 += oldBins[ufbin];
3291 if (oldSumw2) binError0 += oldSumw2[ufbin];
3292 for (xbin = oldxbin; xbin <= nxbins + 1; xbin++) {
3293 ofbin3 = ufbin+xbin;
3294 binContent3 += oldBins[ ofbin3 ];
3295 if (oldSumw2) binError3 += oldSumw2[ofbin3];
3296 for (zbin = oldzbin; zbin <= nzbins + 1; zbin++) {
3298 ofbin4 = xbin + (nxbins+2)*(nybins+2)*zbin+(oldybin2+i)*(nxbins+2);
3299 binContent4 += oldBins[ofbin4];
3300 if (oldSumw2) binError4 += oldSumw2[ofbin4];
3303 for (zbin = oldzbin; zbin <= nzbins + 1; zbin++) {
3304 ofbin2 = (oldybin2+i)*(nxbins+2)+zbin*(nxbins+2)*(nybins+2);
3305 binContent2 += oldBins[ ofbin2 ];
3306 if (oldSumw2) binError2 += oldSumw2[ofbin2];
3312 hnew->
SetBinContent(newxbins+1,ybin,newzbins+1,binContent4);
3319 oldybin2 += nygroup;
3362 if (!resetStat) hnew->
PutStats(stat);
3365 if (oldSumw2)
delete [] oldSumw2;
3396 if (bin < 0)
return;
3459 :
TH3(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
3475 :
TH3(
name,title,nbinsx,
xbins,nbinsy,ybins,nbinsz,zbins)
3489 :
TH3(
name,title,nbinsx,
xbins,nbinsy,ybins,nbinsz,zbins)
3521 if (newval > -128 && newval < 128) {
fArray[bin] =
Char_t(newval);
return;}
3522 if (newval < -127)
fArray[bin] = -127;
3523 if (newval > 127)
fArray[bin] = 127;
3724 :
TH3(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
3740 :
TH3(
name,title,nbinsx,
xbins,nbinsy,ybins,nbinsz,zbins)
3754 :
TH3(
name,title,nbinsx,
xbins,nbinsy,ybins,nbinsz,zbins)
3786 if (newval > -32768 && newval < 32768) {
fArray[bin] =
Short_t(newval);
return;}
3787 if (newval < -32767)
fArray[bin] = -32767;
3788 if (newval > 32767)
fArray[bin] = 32767;
3960 :
TH3(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
3976 :
TH3(
name,title,nbinsx,
xbins,nbinsy,ybins,nbinsz,zbins)
3990 :
TH3(
name,title,nbinsx,
xbins,nbinsy,ybins,nbinsz,zbins)
4022 if (newval > -INT_MAX && newval < INT_MAX) {
fArray[bin] =
Int_t(newval);
return;}
4023 if (newval < -INT_MAX)
fArray[bin] = -INT_MAX;
4024 if (newval > INT_MAX)
fArray[bin] = INT_MAX;
4163 :
TH3(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
4179 :
TH3(
name,title,nbinsx,
xbins,nbinsy,ybins,nbinsz,zbins)
4193 :
TH3(
name,title,nbinsx,
xbins,nbinsy,ybins,nbinsz,zbins)
4378 :
TH3(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
4394 :
TH3(
name,title,nbinsx,
xbins,nbinsy,ybins,nbinsz,zbins)
4408 :
TH3(
name,title,nbinsx,
xbins,nbinsy,ybins,nbinsz,zbins)
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 result
TH3C operator/(TH3C &h1, TH3C &h2)
Operator /.
TH3C operator*(Float_t c1, TH3C &h1)
Operator *.
TH3C operator-(TH3C &h1, TH3C &h2)
Operator -.
TH3C operator+(TH3C &h1, TH3C &h2)
Operator +.
R__EXTERN TRandom * gRandom
Array of chars or bytes (8 bits per element).
void Streamer(TBuffer &) override
Stream a TArrayC object.
void Copy(TArrayC &array) const
void Set(Int_t n) override
Set size of this array to n chars.
Array of doubles (64 bits per element).
void Streamer(TBuffer &) override
Stream a TArrayD object.
void Copy(TArrayD &array) const
void Set(Int_t n) override
Set size of this array to n doubles.
Array of floats (32 bits per element).
void Copy(TArrayF &array) const
void Set(Int_t n) override
Set size of this array to n floats.
void Streamer(TBuffer &) override
Stream a TArrayF object.
Array of integers (32 bits per element).
void Set(Int_t n) override
Set size of this array to n ints.
void Copy(TArrayI &array) const
Array of shorts (16 bits per element).
void Set(Int_t n) override
Set size of this array to n shorts.
void Streamer(TBuffer &) override
Stream a TArrayS object.
void Copy(TArrayS &array) const
Use this attribute class when an object should have 3D capabilities.
virtual void Streamer(TBuffer &)
virtual Color_t GetTitleColor() const
virtual Color_t GetLabelColor() const
virtual Int_t GetNdivisions() const
virtual Color_t GetAxisColor() const
virtual void SetTitleOffset(Float_t offset=1)
Set distance between the axis and the axis title.
virtual Style_t GetTitleFont() const
virtual Float_t GetLabelOffset() const