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);
227 if (!nbentries)
return 0;
230 if (action == 0)
return 0;
231 nbentries = -nbentries;
246 for (
Int_t i=1;i<nbentries;i++) {
254 if (z < zmin) zmin = z;
255 if (z > zmax) zmax = z;
274 for (
Int_t i=0;i<nbentries;i++) {
275 Fill(buffer[4*i+2],buffer[4*i+3],buffer[4*i+4],buffer[4*i+1]);
302 nbentries = -nbentries;
328 Error(
"Fill",
"Invalid signature - do nothing");
343 Int_t binx, biny, binz, bin;
348 if (binx <0 || biny <0 || binz<0)
return -1;
391 Int_t binx, biny, binz, bin;
396 if (binx <0 || biny <0 || binz<0)
return -1;
436 Int_t binx, biny, binz, bin;
441 if (binx <0 || biny <0 || binz<0)
return -1;
484 Int_t binx, biny, binz, bin;
489 if (binx <0 || biny <0 || binz<0)
return -1;
532 Int_t binx, biny, binz, bin;
537 if (binx <0 || biny <0 || binz<0)
return -1;
580 Int_t binx, biny, binz, bin;
585 if (binx <0 || biny <0 || binz<0)
return -1;
628 Int_t binx, biny, binz, bin;
633 if (binx < 0 || biny < 0 || binz < 0)
682 Int_t binx, biny, binz, bin;
687 if (binx <0 || biny <0 || binz<0)
return -1;
731 Int_t binx, biny, binz, bin;
736 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++) {
899 h3->GetRandom3(
x,
y,z,rng);
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 ==
nullptr)
return 0;
1499 const TAxis *xaxis2 = h2->GetXaxis();
1501 const TAxis *yaxis2 = h2->GetYaxis();
1503 const TAxis *zaxis2 = h2->GetZaxis();
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++) {
1570 sum2 += h2->GetBinContent(bin);
1572 Double_t ew2 = h2->GetBinError(bin);
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]);
1632 rsum2 += s2*h2->GetBinContent(bin);
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);
1677 if (
TMath::Abs(rsum2-1) > 0.002)
Warning(
"KolmogorovTest",
"Numerical problems with h2=%s\n",h2->GetName());
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);
1849 bool computeErrors,
bool originalRange,
1850 bool useUF,
bool useOF)
const
1859 Int_t nx = ixmax-ixmin+1;
1865 Error(
"DoProject1D",
"Histogram with name %s must be a TH1D and is a %s",
name,h1obj->
ClassName());
1872 if ( originalRange )
1874 if (bins->
fN == 0) {
1880 if (bins->
fN == 0) {
1890 if ( originalRange )
1892 if (bins->
fN == 0) {
1898 if (bins->
fN == 0) {
1928 if (out1 ==
nullptr && out2 ==
nullptr) {
1940 R__ASSERT(out1 !=
nullptr && out2 !=
nullptr);
1942 Int_t *refX =
nullptr, *refY =
nullptr, *refZ =
nullptr;
1943 Int_t ixbin, out1bin, out2bin;
1959 R__ASSERT (refX !=
nullptr && refY !=
nullptr && refZ !=
nullptr);
1978 for (ixbin=0;ixbin<=1+projX->
GetNbins();ixbin++) {
1985 for (out1bin = out1min; out1bin <= out1max; out1bin++) {
1986 for (out2bin = out2min; out2bin <= out2max; out2bin++) {
1992 if (computeErrors) {
2010 bool resetStats =
true;
2011 double eps = 1.E-12;
2015 bool resetEntries = resetStats;
2017 resetEntries |= !useUF || !useOF;
2024 stats[2] = stats[4];
2025 stats[3] = stats[5];
2028 stats[2] = stats[7];
2029 stats[3] = stats[8];
2058 bool computeErrors,
bool originalRange,
2059 bool useUF,
bool useOF)
const
2069 Int_t nx = ixmax-ixmin+1;
2070 Int_t ny = iymax-iymin+1;
2078 Error(
"DoProject2D",
"Histogram with name %s must be a TH2D and is a %s",
name,h2obj->
ClassName());
2086 if ( originalRange ) {
2093 h2->GetYaxis()->Set(projX->
GetNbins(),&xbins->fArray[ixmin-1]);
2098 h2->GetXaxis()->Set(ny,&ybins->
fArray[iymin-1]);
2100 h2->GetYaxis()->Set(nx,&xbins->fArray[ixmin-1]);
2108 if ( originalRange )
2110 if (xbins->fN == 0 && ybins->
fN == 0) {
2113 }
else if (ybins->
fN == 0) {
2115 ,projX->
GetNbins(),&xbins->fArray[ixmin-1]);
2116 }
else if (xbins->fN == 0) {
2123 if (xbins->fN == 0 && ybins->
fN == 0) {
2126 }
else if (ybins->
fN == 0) {
2128 ,nx,&xbins->fArray[ixmin-1]);
2129 }
else if (xbins->fN == 0) {
2133 h2 =
new TH2D(
name,title,ny,&ybins->
fArray[iymin-1],nx,&xbins->fArray[ixmin-1]);
2142 h2->GetXaxis()->ImportAttributes(projY);
2143 h2->GetYaxis()->ImportAttributes(projX);
2151 h2->GetXaxis()->SetBinLabel(i,lb->
String().
Data());
2160 h2->GetYaxis()->SetBinLabel(i,lb->
String().
Data());
2170 if ( computeErrors && (h2->GetSumw2N() != h2->GetNcells()) ) h2->Sumw2();
2173 const TAxis* out =
nullptr;
2182 Int_t *refX =
nullptr, *refY =
nullptr, *refZ =
nullptr;
2183 Int_t ixbin, iybin, outbin;
2184 if ( projX ==
GetXaxis() && projY ==
GetYaxis() ) { refX = &ixbin; refY = &iybin; refZ = &outbin; }
2185 if ( projX ==
GetYaxis() && projY ==
GetXaxis() ) { refX = &iybin; refY = &ixbin; refZ = &outbin; }
2186 if ( projX ==
GetXaxis() && projY ==
GetZaxis() ) { refX = &ixbin; refY = &outbin; refZ = &iybin; }
2187 if ( projX ==
GetZaxis() && projY ==
GetXaxis() ) { refX = &iybin; refY = &outbin; refZ = &ixbin; }
2188 if ( projX ==
GetYaxis() && projY ==
GetZaxis() ) { refX = &outbin; refY = &ixbin; refZ = &iybin; }
2189 if ( projX ==
GetZaxis() && projY ==
GetYaxis() ) { refX = &outbin; refY = &iybin; refZ = &ixbin; }
2190 R__ASSERT (refX !=
nullptr && refY !=
nullptr && refZ !=
nullptr);
2196 Int_t outmin = out->GetFirst();
2197 Int_t outmax = out->GetLast();
2199 if (outmin == 0 && outmax == 0) { outmin = 1; outmax = out->GetNbins(); }
2204 for (ixbin=0;ixbin<=1+projX->
GetNbins();ixbin++) {
2208 for (iybin=0;iybin<=1+projY->
GetNbins();iybin++) {
2216 for (outbin = outmin; outbin <= outmax; outbin++) {
2222 if (computeErrors) {
2230 h2->SetBinContent(iy , ix, cont);
2231 if (computeErrors) h2->SetBinError(iy, ix,
TMath::Sqrt(err2) );
2240 bool resetStats =
true;
2241 double eps = 1.E-12;
2245 bool resetEntries = resetStats;
2247 resetEntries |= !useUF || !useOF;
2252 for (
Int_t i = 0; i <
kNstat; ++i) { oldst[i] = 0; }
2254 std::copy(oldst,oldst+
kNstat,stats);
2258 stats[4] = oldst[7];
2259 stats[5] = oldst[8];
2260 stats[6] = oldst[9];
2263 stats[2] = oldst[4];
2264 stats[3] = oldst[5];
2266 stats[4] = oldst[2];
2267 stats[5] = oldst[3];
2270 stats[4] = oldst[7];
2271 stats[5] = oldst[8];
2272 stats[6] = oldst[10];
2276 stats[2] = oldst[7];
2277 stats[3] = oldst[8];
2279 stats[4] = oldst[2];
2280 stats[5] = oldst[3];
2281 stats[6] = oldst[9];
2284 stats[4] = oldst[4];
2285 stats[5] = oldst[5];
2286 stats[6] = oldst[10];
2290 h2->PutStats(stats);
2300 Double_t entries = h2->GetEffectiveEntries();
2301 if (!computeErrors) entries =
TMath::Floor( entries + 0.5);
2302 h2->SetEntries( entries );
2368 if (opt.
Contains(
"x")) { pcase = 1; ptype =
"x"; }
2369 if (opt.
Contains(
"y")) { pcase = 2; ptype =
"y"; }
2370 if (opt.
Contains(
"z")) { pcase = 3; ptype =
"z"; }
2371 if (opt.
Contains(
"xy")) { pcase = 4; ptype =
"xy"; }
2372 if (opt.
Contains(
"yx")) { pcase = 5; ptype =
"yx"; }
2373 if (opt.
Contains(
"xz")) { pcase = 6; ptype =
"xz"; }
2374 if (opt.
Contains(
"zx")) { pcase = 7; ptype =
"zx"; }
2375 if (opt.
Contains(
"yz")) { pcase = 8; ptype =
"yz"; }
2376 if (opt.
Contains(
"zy")) { pcase = 9; ptype =
"zy"; }
2379 Error(
"Project3D",
"No projection axis specified - return a NULL pointer");
2386 computeErrors =
kTRUE;
2403 originalRange =
kTRUE;
2414 title +=
" "; title += ptype; title +=
" projection";
2420 computeErrors, originalRange, useUF, useOF);
2426 computeErrors, originalRange, useUF, useOF);
2432 computeErrors, originalRange, useUF, useOF);
2438 computeErrors, originalRange, useUF, useOF);
2444 computeErrors, originalRange, useUF, useOF);
2450 computeErrors, originalRange, useUF, useOF);
2456 computeErrors, originalRange, useUF, useOF);
2462 computeErrors, originalRange, useUF, useOF);
2468 computeErrors, originalRange, useUF, useOF);
2499 if (useWeights && binSumw2.
fN <= 0) useWeights =
false;
2506 if (outBin <0)
return;
2508 if ( useWeights ) tmp = binSumw2.
fArray[outBin];
2509 p2->
Fill( u ,
v,
w, cont);
2519 bool originalRange,
bool useUF,
bool useOF)
const
2527 Int_t nx = ixmax-ixmin+1;
2528 Int_t ny = iymax-iymin+1;
2538 Error(
"DoProjectProfile2D",
"Histogram with name %s must be a TProfile2D and is a %s",
name,p2obj->
ClassName());
2546 if ( originalRange ) {
2567 if ( originalRange ) {
2568 if (xbins->fN == 0 && ybins->
fN == 0) {
2571 }
else if (ybins->
fN == 0) {
2573 ,projX->
GetNbins(),&xbins->fArray[ixmin-1]);
2574 }
else if (xbins->fN == 0) {
2581 if (xbins->fN == 0 && ybins->
fN == 0) {
2584 }
else if (ybins->
fN == 0) {
2586 ,nx,&xbins->fArray[ixmin-1]);
2587 }
else if (xbins->fN == 0) {
2597 const TAxis* outAxis =
nullptr;
2612 Int_t *refX =
nullptr, *refY =
nullptr, *refZ =
nullptr;
2613 Int_t ixbin, iybin, outbin;
2614 if ( projX ==
GetXaxis() && projY ==
GetYaxis() ) { refX = &ixbin; refY = &iybin; refZ = &outbin; }
2615 if ( projX ==
GetYaxis() && projY ==
GetXaxis() ) { refX = &iybin; refY = &ixbin; refZ = &outbin; }
2616 if ( projX ==
GetXaxis() && projY ==
GetZaxis() ) { refX = &ixbin; refY = &outbin; refZ = &iybin; }
2617 if ( projX ==
GetZaxis() && projY ==
GetXaxis() ) { refX = &iybin; refY = &outbin; refZ = &ixbin; }
2618 if ( projX ==
GetYaxis() && projY ==
GetZaxis() ) { refX = &outbin; refY = &ixbin; refZ = &iybin; }
2619 if ( projX ==
GetZaxis() && projY ==
GetYaxis() ) { refX = &outbin; refY = &iybin; refZ = &ixbin; }
2620 R__ASSERT (refX !=
nullptr && refY !=
nullptr && refZ !=
nullptr);
2630 if (useWeights && binSumw2.
fN <= 0) useWeights =
false;
2634 for (ixbin=0;ixbin<=1+projX->
GetNbins();ixbin++) {
2636 for ( iybin=0;iybin<=1+projY->
GetNbins();iybin++) {
2641 if (poutBin <0)
continue;
2643 for (outbin = outmin; outbin <= outmax; outbin++) {
2650 if (!cont)
continue;
2654 if ( useWeights ) tmp = binSumw2.
fArray[poutBin];
2664 bool resetStats =
true;
2675 if (!useWeights) entries =
TMath::Floor( entries + 0.5);
2729 if (opt.
Contains(
"xy")) { pcase = 4; ptype =
"xy"; }
2730 if (opt.
Contains(
"yx")) { pcase = 5; ptype =
"yx"; }
2731 if (opt.
Contains(
"xz")) { pcase = 6; ptype =
"xz"; }
2732 if (opt.
Contains(
"zx")) { pcase = 7; ptype =
"zx"; }
2733 if (opt.
Contains(
"yz")) { pcase = 8; ptype =
"yz"; }
2734 if (opt.
Contains(
"zy")) { pcase = 9; ptype =
"zy"; }
2737 Error(
"Project3D",
"No projection axis specified - return a NULL pointer");
2755 originalRange =
kTRUE;
2764 title +=
" profile "; title += ptype; title +=
" projection";
2826 return Rebin3D(ngroup, 1, 1, newname);
2836 return Rebin3D(1, ngroup, 1, newname);
2846 return Rebin3D(1, 1, ngroup, newname);
2877 Int_t i,j,k,xbin,ybin,zbin;
2887 if ((nxgroup <= 0) || (nxgroup > nxbins)) {
2888 Error(
"Rebin",
"Illegal value of nxgroup=%d",nxgroup);
2891 if ((nygroup <= 0) || (nygroup > nybins)) {
2892 Error(
"Rebin",
"Illegal value of nygroup=%d",nygroup);
2895 if ((nzgroup <= 0) || (nzgroup > nzbins)) {
2896 Error(
"Rebin",
"Illegal value of nzgroup=%d",nzgroup);
2900 Int_t newxbins = nxbins/nxgroup;
2901 Int_t newybins = nybins/nygroup;
2902 Int_t newzbins = nzbins/nzgroup;
2920 if (newname && strlen(newname)) {
2928 bool resetStat =
false;
2932 if (newxbins*nxgroup != nxbins) {
2936 if (newybins*nygroup != nybins) {
2940 if (newzbins*nzgroup != nzbins) {
2982 if (nxgroup != 1 || nygroup != 1 || nzgroup != 1) {
2991 hnew->
SetBins(newxbins,xbins, newybins, ybins, newzbins, zbins);
3004 for (xbin = 1; xbin <= newxbins; xbin++) {
3007 for (ybin = 1; ybin <= newybins; ybin++) {
3009 for (zbin = 1; zbin <= newzbins; zbin++) {
3012 for (i = 0; i < nxgroup; i++) {
3013 if (oldxbin+i > nxbins)
break;
3014 for (j =0; j < nygroup; j++) {
3015 if (oldybin+j > nybins)
break;
3016 for (k =0; k < nzgroup; k++) {
3017 if (oldzbin+k > nzbins)
break;
3019 bin = oldxbin + i + (oldybin + j)*(nxbins + 2) + (oldzbin + k)*(nxbins + 2)*(nybins + 2);
3020 binContent += oldBins[bin];
3021 if (oldSumw2) binSumw2 += oldSumw2[bin];
3036 for (
Int_t xover = 0; xover <= 1; xover++) {
3037 for (
Int_t yover = 0; yover <= 1; yover++) {
3038 for (
Int_t zover = 0; zover <= 1; zover++) {
3042 for (xbin = xover*oldxbin; xbin <= xover*(nxbins+1); xbin++) {
3043 for (ybin = yover*oldybin; ybin <= yover*(nybins+1); ybin++) {
3044 for (zbin = zover*oldzbin; zbin <= zover*(nzbins+1); zbin++) {
3045 bin =
GetBin(xbin,ybin,zbin);
3046 binContent += oldBins[bin];
3047 if (oldSumw2) binSumw2 += oldSumw2[bin];
3052 yover*(newybins+1), zover*(newzbins+1) );
3059 Double_t binContent0, binContent2, binContent3, binContent4;
3060 Double_t binError0, binError2, binError3, binError4;
3061 Int_t oldxbin2, oldybin2, oldzbin2;
3062 Int_t ufbin, ofbin, ofbin2, ofbin3, ofbin4;
3068 for (xbin = 1; xbin<=newxbins; xbin++) {
3070 for (zbin = 1; zbin<=newzbins; zbin++) {
3071 binContent0 = binContent2 = 0;
3072 binError0 = binError2 = 0;
3073 for (i=0; i<nxgroup; i++) {
3074 if (oldxbin2+i > nxbins)
break;
3075 for (k=0; k<nzgroup; k++) {
3076 if (oldzbin2+k > nzbins)
break;
3078 ufbin = oldxbin2 + i + (nxbins+2)*(nybins+2)*(oldzbin2+k);
3079 binContent0 += oldBins[ufbin];
3080 if (oldSumw2) binError0 += oldSumw2[ufbin];
3081 for (ybin = oldybin; ybin <= nybins + 1; ybin++) {
3083 ofbin = ufbin + ybin*(nxbins+2);
3084 binContent2 += oldBins[ofbin];
3085 if (oldSumw2) binError2 += oldSumw2[ofbin];
3095 oldzbin2 += nzgroup;
3097 oldxbin2 += nxgroup;
3104 for (ybin = 1; ybin<=newybins; ybin++) {
3106 for (zbin = 1; zbin<=newzbins; zbin++) {
3107 binContent0 = binContent2 = 0;
3108 binError0 = binError2 = 0;
3109 for (j=0; j<nygroup; j++) {
3110 if (oldybin2+j > nybins)
break;
3111 for (k=0; k<nzgroup; k++) {
3112 if (oldzbin2+k > nzbins)
break;
3114 ufbin = (oldybin2 + j)*(nxbins+2) + (nxbins+2)*(nybins+2)*(oldzbin2+k);
3115 binContent0 += oldBins[ufbin];
3116 if (oldSumw2) binError0 += oldSumw2[ufbin];
3117 for (xbin = oldxbin; xbin <= nxbins + 1; xbin++) {
3119 ofbin = ufbin + xbin;
3120 binContent2 += oldBins[ofbin];
3121 if (oldSumw2) binError2 += oldSumw2[ofbin];
3131 oldzbin2 += nzgroup;
3133 oldybin2 += nygroup;
3140 for (xbin = 1; xbin<=newxbins; xbin++) {
3142 for (ybin = 1; ybin<=newybins; ybin++) {
3143 binContent0 = binContent2 = 0;
3144 binError0 = binError2 = 0;
3145 for (i=0; i<nxgroup; i++) {
3146 if (oldxbin2+i > nxbins)
break;
3147 for (j=0; j<nygroup; j++) {
3148 if (oldybin2+j > nybins)
break;
3150 ufbin = oldxbin2 + i + (nxbins+2)*(oldybin2+j);
3151 binContent0 += oldBins[ufbin];
3152 if (oldSumw2) binError0 += oldSumw2[ufbin];
3153 for (zbin = oldzbin; zbin <= nzbins + 1; zbin++) {
3155 ofbin = ufbin + (nxbins+2)*(nybins+2)*zbin;
3156 binContent2 += oldBins[ofbin];
3157 if (oldSumw2) binError2 += oldSumw2[ofbin];
3167 oldybin2 += nygroup;
3169 oldxbin2 += nxgroup;
3176 for (xbin = 1; xbin<=newxbins; xbin++) {
3185 for (i=0; i<nxgroup; i++) {
3186 if (oldxbin2+i > nxbins)
break;
3187 ufbin = oldxbin2 + i;
3188 binContent0 += oldBins[ufbin];
3189 if (oldSumw2) binError0 += oldSumw2[ufbin];
3190 for (ybin = oldybin; ybin <= nybins + 1; ybin++) {
3191 ofbin3 = ufbin+ybin*(nxbins+2);
3192 binContent3 += oldBins[ ofbin3 ];
3193 if (oldSumw2) binError3 += oldSumw2[ofbin3];
3194 for (zbin = oldzbin; zbin <= nzbins + 1; zbin++) {
3196 ofbin4 = oldxbin2 + i + ybin*(nxbins+2) + (nxbins+2)*(nybins+2)*zbin;
3197 binContent4 += oldBins[ofbin4];
3198 if (oldSumw2) binError4 += oldSumw2[ofbin4];
3201 for (zbin = oldzbin; zbin <= nzbins + 1; zbin++) {
3202 ofbin2 = ufbin+zbin*(nxbins+2)*(nybins+2);
3203 binContent2 += oldBins[ ofbin2 ];
3204 if (oldSumw2) binError2 += oldSumw2[ofbin2];
3210 hnew->
SetBinContent(xbin,newybins+1,newzbins+1,binContent4);
3217 oldxbin2 += nxgroup;
3224 for (zbin = 1; zbin<=newzbins; zbin++) {
3233 for (i=0; i<nzgroup; i++) {
3234 if (oldzbin2+i > nzbins)
break;
3235 ufbin = (oldzbin2 + i)*(nxbins+2)*(nybins+2);
3236 binContent0 += oldBins[ufbin];
3237 if (oldSumw2) binError0 += oldSumw2[ufbin];
3238 for (ybin = oldybin; ybin <= nybins + 1; ybin++) {
3239 ofbin3 = ufbin+ybin*(nxbins+2);
3240 binContent3 += oldBins[ ofbin3 ];
3241 if (oldSumw2) binError3 += oldSumw2[ofbin3];
3242 for (xbin = oldxbin; xbin <= nxbins + 1; xbin++) {
3244 ofbin4 = ufbin + xbin + ybin*(nxbins+2);
3245 binContent4 += oldBins[ofbin4];
3246 if (oldSumw2) binError4 += oldSumw2[ofbin4];
3249 for (xbin = oldxbin; xbin <= nxbins + 1; xbin++) {
3250 ofbin2 = xbin +(oldzbin2+i)*(nxbins+2)*(nybins+2);
3251 binContent2 += oldBins[ ofbin2 ];
3252 if (oldSumw2) binError2 += oldSumw2[ofbin2];
3258 hnew->
SetBinContent(newxbins+1,newybins+1,zbin,binContent4);
3265 oldzbin2 += nzgroup;
3272 for (ybin = 1; ybin<=newybins; ybin++) {
3281 for (i=0; i<nygroup; i++) {
3282 if (oldybin2+i > nybins)
break;
3283 ufbin = (oldybin2 + i)*(nxbins+2);
3284 binContent0 += oldBins[ufbin];
3285 if (oldSumw2) binError0 += oldSumw2[ufbin];
3286 for (xbin = oldxbin; xbin <= nxbins + 1; xbin++) {
3287 ofbin3 = ufbin+xbin;
3288 binContent3 += oldBins[ ofbin3 ];
3289 if (oldSumw2) binError3 += oldSumw2[ofbin3];
3290 for (zbin = oldzbin; zbin <= nzbins + 1; zbin++) {
3292 ofbin4 = xbin + (nxbins+2)*(nybins+2)*zbin+(oldybin2+i)*(nxbins+2);
3293 binContent4 += oldBins[ofbin4];
3294 if (oldSumw2) binError4 += oldSumw2[ofbin4];
3297 for (zbin = oldzbin; zbin <= nzbins + 1; zbin++) {
3298 ofbin2 = (oldybin2+i)*(nxbins+2)+zbin*(nxbins+2)*(nybins+2);
3299 binContent2 += oldBins[ ofbin2 ];
3300 if (oldSumw2) binError2 += oldSumw2[ofbin2];
3306 hnew->
SetBinContent(newxbins+1,ybin,newzbins+1,binContent4);
3313 oldybin2 += nygroup;
3356 if (!resetStat) hnew->
PutStats(stat);
3359 if (oldSumw2)
delete [] oldSumw2;
3390 if (bin < 0)
return;
3423 bool originalRange,
bool useUF,
bool useOF)
3425 return h.DoProject1D(
name, title, projX,
nullptr,
nullptr, computeErrors, originalRange, useUF, useOF);
3432 bool computeErrors,
bool originalRange,
bool useUF,
bool useOF)
3434 return h.DoProject2D(
name, title, projX, projY, computeErrors, originalRange, useUF, useOF);
3470 :
TH3(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
3486 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
3500 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
3513 h3c.TH3C::Copy(*
this);
3532 if (newval > -128 && newval < 128) {
fArray[bin] =
Char_t(newval);
return;}
3533 if (newval < -127)
fArray[bin] = -127;
3534 if (newval > 127)
fArray[bin] = 127;
3638 h3c.TH3C::Copy(*
this);
3736 :
TH3(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
3752 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
3766 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
3779 h3s.TH3S::Copy(*
this);
3798 if (newval > -32768 && newval < 32768) {
fArray[bin] =
Short_t(newval);
return;}
3799 if (newval < -32767)
fArray[bin] = -32767;
3800 if (newval > 32767)
fArray[bin] = 32767;
3875 h3s.TH3S::Copy(*
this);
3973 :
TH3(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
3989 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
4003 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
4016 h3i.TH3I::Copy(*
this);
4035 if (newval > -INT_MAX && newval < INT_MAX) {
fArray[bin] =
Int_t(newval);
return;}
4036 if (newval < -INT_MAX)
fArray[bin] = -INT_MAX;
4037 if (newval > INT_MAX)
fArray[bin] = INT_MAX;
4079 h3i.TH3I::Copy(*
this);
4177 :
TH3(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
4193 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
4207 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
4220 h3f.TH3F::Copy(*
this);
4295 h3f.TH3F::Copy(*
this);
4393 :
TH3(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
4409 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
4423 :
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 &h3c)
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 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 Set(Int_t n) override
Set size of this array to n doubles.
Array of floats (32 bits per element).
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.
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.
virtual void Set(Int_t n)=0
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
virtual void SetAxisColor(Color_t color=1, Float_t alpha=1.)
Set color of the line axis and tick marks.
virtual void SetLabelSize(Float_t size=0.04)
Set size of axis labels.
virtual Style_t GetLabelFont() const
virtual void SetTitleFont(Style_t font=62)
Set the title font.
virtual void SetLabelOffset(Float_t offset=0.005)
Set distance between the axis and the labels.
virtual void SetLabelFont(Style_t font=62)
Set labels' font.
virtual void SetTitleSize(Float_t size=0.04)
Set size of axis title.
virtual void SetTitleColor(Color_t color=1)
Set color of axis title.
virtual Float_t GetTitleSize() const
virtual Float_t GetLabelSize() const
virtual Float_t GetTickLength() const
virtual Float_t GetTitleOffset() const
virtual void SetTickLength(Float_t length=0.03)
Set tick mark length.
virtual void SetNdivisions(Int_t n=510, Bool_t optim=kTRUE)
Set the number of divisions for this axis.
virtual void SetLabelColor(Color_t color=1, Float_t alpha=1.)
Set color of labels.
virtual Color_t GetFillColor() const
Return the fill area color.
virtual void SetFillColor(Color_t fcolor)
Set the fill area color.
virtual Color_t GetLineColor() const
Return the line color.
virtual void SetLineColor(Color_t lcolor)
Set the line color.
virtual Style_t GetMarkerStyle() const
Return the marker style.
virtual void SetMarkerColor(Color_t mcolor=1)
Set the marker color.
virtual Color_t GetMarkerColor() const
Return the marker color.
virtual void SetMarkerStyle(Style_t mstyle=1)
Set the marker style.