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");
1137 if (sumw == 0)
return 0;
1138 if (axis1 == 1 && axis2 == 1) {
1139 return TMath::Abs(sumwx2/sumw - sumwx*sumwx/(sumw*sumw));
1141 if (axis1 == 2 && axis2 == 2) {
1142 return TMath::Abs(sumwy2/sumw - sumwy*sumwy/(sumw*sumw));
1144 if (axis1 == 3 && axis2 == 3) {
1145 return TMath::Abs(sumwz2/sumw - sumwz*sumwz/(sumw*sumw));
1147 if ((axis1 == 1 && axis2 == 2) || (axis1 == 2 && axis2 == 1)) {
1148 return sumwxy/sumw - sumwx*sumwy/(sumw*sumw);
1150 if ((axis1 == 1 && axis2 == 3) || (axis1 == 3 && axis2 == 1)) {
1151 return sumwxz/sumw - sumwx*sumwz/(sumw*sumw);
1153 if ((axis1 == 2 && axis2 == 3) || (axis1 == 3 && axis2 == 2)) {
1154 return sumwyz/sumw - sumwy*sumwz/(sumw*sumw);
1173 Int_t nxy = nbinsx*nbinsy;
1174 Int_t nbins = nxy*nbinsz;
1183 if (integral == 0 ) {
x = 0;
y = 0; z = 0;
return;}
1190 Int_t binz = ibin/nxy;
1191 Int_t biny = (ibin - nxy*binz)/nbinsx;
1192 Int_t binx = ibin - nbinsx*(biny + nbinsy*binz);
1220 Int_t bin, binx, biny, binz;
1224 for (bin=0;bin<11;bin++) stats[bin] = 0;
1235 if (firstBinX == 1) firstBinX = 0;
1239 if (firstBinY == 1) firstBinY = 0;
1243 if (firstBinZ == 1) firstBinZ = 0;
1253 for (binz = firstBinZ; binz <= lastBinZ; binz++) {
1255 for (biny = firstBinY; biny <= lastBinY; biny++) {
1257 for (binx = firstBinX; binx <= lastBinX; binx++) {
1258 bin =
GetBin(binx,biny,binz);
1264 stats[1] += err*err;
1318 return DoIntegral(binx1,binx2,biny1,biny2,binz1,binz2,err,option);
1334 return DoIntegral(binx1,binx2,biny1,biny2,binz1,binz2,error,option,
kTRUE);
1342 Error(
"Interpolate",
"This function must be called with 3 arguments for a TH3");
1352 Error(
"Interpolate",
"This function must be called with 3 arguments for a TH3");
1371 Int_t obx = ubx + 1;
1375 Int_t oby = uby + 1;
1379 Int_t obz = ubz + 1;
1384 if (ubx <=0 || uby <=0 || ubz <= 0 ||
1386 Error(
"Interpolate",
"Cannot interpolate outside histogram domain.");
1411 Double_t w1 = i1 * (1 - yd) + i2 * yd;
1412 Double_t w2 = j1 * (1 - yd) + j2 * yd;
1415 Double_t result = w1 * (1 - xd) + w2 * xd;
1448 if (h2 == 0)
return 0;
1464 Error(
"KolmogorovTest",
"Histograms must be 3-D\n");
1470 Error(
"KolmogorovTest",
"Number of channels in X is different, %d and %d\n",ncx1,ncx2);
1474 Error(
"KolmogorovTest",
"Number of channels in Y is different, %d and %d\n",ncy1,ncy2);
1478 Error(
"KolmogorovTest",
"Number of channels in Z is different, %d and %d\n",ncz1,ncz2);
1488 if (diff1 > difprec || diff2 > difprec) {
1489 Error(
"KolmogorovTest",
"histograms with different binning along X");
1494 if (diff1 > difprec || diff2 > difprec) {
1495 Error(
"KolmogorovTest",
"histograms with different binning along Y");
1500 if (diff1 > difprec || diff2 > difprec) {
1501 Error(
"KolmogorovTest",
"histograms with different binning along Z");
1506 Int_t ibeg = 1, jbeg = 1, kbeg = 1;
1507 Int_t iend = ncx1, jend = ncy1, kend = ncz1;
1508 if (opt.
Contains(
"U")) {ibeg = 0; jbeg = 0; kbeg = 0;}
1509 if (opt.
Contains(
"O")) {iend = ncx1+1; jend = ncy1+1; kend = ncz1+1;}
1516 for (i = ibeg; i <= iend; i++) {
1517 for (j = jbeg; j <= jend; j++) {
1518 for (k = kbeg; k <= kend; k++) {
1533 Error(
"KolmogorovTest",
"Integral is zero for h1=%s\n",
h1->
GetName());
1537 Error(
"KolmogorovTest",
"Integral is zero for h2=%s\n",h2->
GetName());
1545 esum1 = sum1 * sum1 / w1;
1550 esum2 = sum2 * sum2 / w2;
1554 if (afunc2 && afunc1) {
1555 Error(
"KolmogorovTest",
"Errors are zero for both histograms\n");
1561 int order[3] = {0,1,2};
1565 binbeg[0] = ibeg; binbeg[1] = jbeg; binbeg[2] = kbeg;
1566 binend[0] = iend; binend[1] = jend; binend[2] = kend;
1575 for (i = binbeg[order[0] ]; i <= binend[order[0] ]; i++) {
1576 for ( j = binbeg[order[1] ]; j <= binend[order[1] ]; j++) {
1577 for ( k = binbeg[order[2] ]; k <= binend[order[2] ]; k++) {
1578 ibin[ order[0] ] = i;
1579 ibin[ order[1] ] = j;
1580 ibin[ order[2] ] = k;
1581 bin =
h1->
GetBin(ibin[0],ibin[1],ibin[2]);
1588 vdfmax[icomb] = dmax;
1607 if (opt.
Contains(
"N") && !(afunc1 || afunc2 ) ) {
1611 Double_t chi2 = d12*d12/(esum1+esum2);
1614 if (prb > 0 && prb2 > 0) prb = prb*prb2*(1-
TMath::Log(prb*prb2));
1620 printf(
" Kolmo Prob h1 = %s, sum1=%g\n",
h1->
GetName(),sum1);
1621 printf(
" Kolmo Prob h2 = %s, sum2=%g\n",h2->
GetName(),sum2);
1622 printf(
" Kolmo Probabil = %f, Max Dist = %g\n",prb,dfmax);
1624 printf(
" Kolmo Probabil = %f for shape alone, =%f for normalisation alone\n",prb1,prb2);
1630 if (opt.
Contains(
"M"))
return dfmax;
1763 computeErrors =
kTRUE;
1768 originalRange =
kTRUE;
1772 TH1D *
h1 =
DoProject1D(
name, title, projAxis, &out1, &out2, computeErrors, originalRange,
true,
true);
1791 if (padsav) padsav->
cd();
1803 bool computeErrors,
bool originalRange,
1804 bool useUF,
bool useOF)
const
1813 Int_t nx = ixmax-ixmin+1;
1818 if (h1obj->IsA() != TH1D::Class() ) {
1819 Error(
"DoProject1D",
"Histogram with name %s must be a TH1D and is a %s",
name,h1obj->
ClassName());
1826 if ( originalRange )
1828 if (bins->
fN == 0) {
1834 if (bins->
fN == 0) {
1844 if ( originalRange )
1846 if (bins->
fN == 0) {
1852 if (bins->
fN == 0) {
1882 if (out1 ==
nullptr && out2 ==
nullptr) {
1894 R__ASSERT(out1 !=
nullptr && out2 !=
nullptr);
1896 Int_t *refX = 0, *refY = 0, *refZ = 0;
1897 Int_t ixbin, out1bin, out2bin;
1913 R__ASSERT (refX != 0 && refY != 0 && refZ != 0);
1932 for (ixbin=0;ixbin<=1+projX->
GetNbins();ixbin++) {
1939 for (out1bin = out1min; out1bin <= out1max; out1bin++) {
1940 for (out2bin = out2min; out2bin <= out2max; out2bin++) {
1946 if (computeErrors) {
1964 bool resetStats =
true;
1965 double eps = 1.E-12;
1966 if (IsA() == TH3F::Class() ) eps = 1.E-6;
1969 bool resetEntries = resetStats;
1971 resetEntries |= !useUF || !useOF;
1978 stats[2] = stats[4];
1979 stats[3] = stats[5];
1982 stats[2] = stats[7];
1983 stats[3] = stats[8];
2012 bool computeErrors,
bool originalRange,
2013 bool useUF,
bool useOF)
const
2023 Int_t nx = ixmax-ixmin+1;
2024 Int_t ny = iymax-iymin+1;
2031 if ( h2obj->IsA() != TH2D::Class() ) {
2032 Error(
"DoProject2D",
"Histogram with name %s must be a TH2D and is a %s",
name,h2obj->
ClassName());
2040 if ( originalRange ) {
2062 if ( originalRange )
2064 if (xbins->fN == 0 && ybins->
fN == 0) {
2067 }
else if (ybins->
fN == 0) {
2069 ,projX->
GetNbins(),&xbins->fArray[ixmin-1]);
2070 }
else if (xbins->fN == 0) {
2077 if (xbins->fN == 0 && ybins->
fN == 0) {
2080 }
else if (ybins->
fN == 0) {
2082 ,nx,&xbins->fArray[ixmin-1]);
2083 }
else if (xbins->fN == 0) {
2087 h2 =
new TH2D(
name,title,ny,&ybins->
fArray[iymin-1],nx,&xbins->fArray[ixmin-1]);
2127 const TAxis* out = 0;
2136 Int_t *refX = 0, *refY = 0, *refZ = 0;
2137 Int_t ixbin, iybin, outbin;
2138 if ( projX ==
GetXaxis() && projY ==
GetYaxis() ) { refX = &ixbin; refY = &iybin; refZ = &outbin; }
2139 if ( projX ==
GetYaxis() && projY ==
GetXaxis() ) { refX = &iybin; refY = &ixbin; refZ = &outbin; }
2140 if ( projX ==
GetXaxis() && projY ==
GetZaxis() ) { refX = &ixbin; refY = &outbin; refZ = &iybin; }
2141 if ( projX ==
GetZaxis() && projY ==
GetXaxis() ) { refX = &iybin; refY = &outbin; refZ = &ixbin; }
2142 if ( projX ==
GetYaxis() && projY ==
GetZaxis() ) { refX = &outbin; refY = &ixbin; refZ = &iybin; }
2143 if ( projX ==
GetZaxis() && projY ==
GetYaxis() ) { refX = &outbin; refY = &iybin; refZ = &ixbin; }
2144 R__ASSERT (refX != 0 && refY != 0 && refZ != 0);
2153 if (outmin == 0 && outmax == 0) { outmin = 1; outmax = out->
GetNbins(); }
2158 for (ixbin=0;ixbin<=1+projX->
GetNbins();ixbin++) {
2162 for (iybin=0;iybin<=1+projY->
GetNbins();iybin++) {
2170 for (outbin = outmin; outbin <= outmax; outbin++) {
2176 if (computeErrors) {
2194 bool resetStats =
true;
2195 double eps = 1.E-12;
2196 if (IsA() == TH3F::Class() ) eps = 1.E-6;
2199 bool resetEntries = resetStats;
2201 resetEntries |= !useUF || !useOF;
2206 for (
Int_t i = 0; i <
kNstat; ++i) { oldst[i] = 0; }
2208 std::copy(oldst,oldst+
kNstat,stats);
2212 stats[4] = oldst[7];
2213 stats[5] = oldst[8];
2214 stats[6] = oldst[9];
2217 stats[2] = oldst[4];
2218 stats[3] = oldst[5];
2220 stats[4] = oldst[2];
2221 stats[5] = oldst[3];
2224 stats[4] = oldst[7];
2225 stats[5] = oldst[8];
2226 stats[6] = oldst[10];
2230 stats[2] = oldst[7];
2231 stats[3] = oldst[8];
2233 stats[4] = oldst[2];
2234 stats[5] = oldst[3];
2235 stats[6] = oldst[9];
2238 stats[4] = oldst[4];
2239 stats[5] = oldst[5];
2240 stats[6] = oldst[10];
2255 if (!computeErrors) entries =
TMath::Floor( entries + 0.5);
2318 if (opt.
Contains(
"x")) { pcase = 1; ptype =
"x"; }
2319 if (opt.
Contains(
"y")) { pcase = 2; ptype =
"y"; }
2320 if (opt.
Contains(
"z")) { pcase = 3; ptype =
"z"; }
2321 if (opt.
Contains(
"xy")) { pcase = 4; ptype =
"xy"; }
2322 if (opt.
Contains(
"yx")) { pcase = 5; ptype =
"yx"; }
2323 if (opt.
Contains(
"xz")) { pcase = 6; ptype =
"xz"; }
2324 if (opt.
Contains(
"zx")) { pcase = 7; ptype =
"zx"; }
2325 if (opt.
Contains(
"yz")) { pcase = 8; ptype =
"yz"; }
2326 if (opt.
Contains(
"zy")) { pcase = 9; ptype =
"zy"; }
2329 Error(
"Project3D",
"No projection axis specified - return a NULL pointer");
2336 computeErrors =
kTRUE;
2353 originalRange =
kTRUE;
2364 title +=
" "; title += ptype; title +=
" projection";
2370 computeErrors, originalRange, useUF, useOF);
2376 computeErrors, originalRange, useUF, useOF);
2382 computeErrors, originalRange, useUF, useOF);
2388 computeErrors, originalRange, useUF, useOF);
2394 computeErrors, originalRange, useUF, useOF);
2400 computeErrors, originalRange, useUF, useOF);
2406 computeErrors, originalRange, useUF, useOF);
2412 computeErrors, originalRange, useUF, useOF);
2418 computeErrors, originalRange, useUF, useOF);
2434 if (padsav) padsav->
cd();
2452 if (useWeights && binSumw2.
fN <= 0) useWeights =
false;
2459 if (outBin <0)
return;
2461 if ( useWeights ) tmp = binSumw2.
fArray[outBin];
2462 p2->
Fill( u ,
v, w, cont);
2472 bool originalRange,
bool useUF,
bool useOF)
const
2480 Int_t nx = ixmax-ixmin+1;
2481 Int_t ny = iymax-iymin+1;
2490 if (p2obj->IsA() != TProfile2D::Class() ) {
2491 Error(
"DoProjectProfile2D",
"Histogram with name %s must be a TProfile2D and is a %s",
name,p2obj->
ClassName());
2499 if ( originalRange ) {
2520 if ( originalRange ) {
2521 if (xbins->fN == 0 && ybins->
fN == 0) {
2524 }
else if (ybins->
fN == 0) {
2526 ,projX->
GetNbins(),&xbins->fArray[ixmin-1]);
2527 }
else if (xbins->fN == 0) {
2534 if (xbins->fN == 0 && ybins->
fN == 0) {
2537 }
else if (ybins->
fN == 0) {
2539 ,nx,&xbins->fArray[ixmin-1]);
2540 }
else if (xbins->fN == 0) {
2550 const TAxis* outAxis = 0;
2565 Int_t *refX = 0, *refY = 0, *refZ = 0;
2566 Int_t ixbin, iybin, outbin;
2567 if ( projX ==
GetXaxis() && projY ==
GetYaxis() ) { refX = &ixbin; refY = &iybin; refZ = &outbin; }
2568 if ( projX ==
GetYaxis() && projY ==
GetXaxis() ) { refX = &iybin; refY = &ixbin; refZ = &outbin; }
2569 if ( projX ==
GetXaxis() && projY ==
GetZaxis() ) { refX = &ixbin; refY = &outbin; refZ = &iybin; }
2570 if ( projX ==
GetZaxis() && projY ==
GetXaxis() ) { refX = &iybin; refY = &outbin; refZ = &ixbin; }
2571 if ( projX ==
GetYaxis() && projY ==
GetZaxis() ) { refX = &outbin; refY = &ixbin; refZ = &iybin; }
2572 if ( projX ==
GetZaxis() && projY ==
GetYaxis() ) { refX = &outbin; refY = &iybin; refZ = &ixbin; }
2573 R__ASSERT (refX != 0 && refY != 0 && refZ != 0);
2583 if (useWeights && binSumw2.
fN <= 0) useWeights =
false;
2587 for (ixbin=0;ixbin<=1+projX->
GetNbins();ixbin++) {
2589 for ( iybin=0;iybin<=1+projY->
GetNbins();iybin++) {
2594 if (poutBin <0)
continue;
2596 for (outbin = outmin; outbin <= outmax; outbin++) {
2603 if (!cont)
continue;
2607 if ( useWeights ) tmp = binSumw2.
fArray[poutBin];
2617 bool resetStats =
true;
2628 if (!useWeights) entries =
TMath::Floor( entries + 0.5);
2682 if (opt.
Contains(
"xy")) { pcase = 4; ptype =
"xy"; }
2683 if (opt.
Contains(
"yx")) { pcase = 5; ptype =
"yx"; }
2684 if (opt.
Contains(
"xz")) { pcase = 6; ptype =
"xz"; }
2685 if (opt.
Contains(
"zx")) { pcase = 7; ptype =
"zx"; }
2686 if (opt.
Contains(
"yz")) { pcase = 8; ptype =
"yz"; }
2687 if (opt.
Contains(
"zy")) { pcase = 9; ptype =
"zy"; }
2690 Error(
"Project3D",
"No projection axis specified - return a NULL pointer");
2708 originalRange =
kTRUE;
2717 title +=
" profile "; title += ptype; title +=
" projection";
2779 return Rebin3D(ngroup, 1, 1, newname);
2789 return Rebin3D(1, ngroup, 1, newname);
2799 return Rebin3D(1, 1, ngroup, newname);
2830 Int_t i,j,k,xbin,ybin,zbin;
2840 if ((nxgroup <= 0) || (nxgroup > nxbins)) {
2841 Error(
"Rebin",
"Illegal value of nxgroup=%d",nxgroup);
2844 if ((nygroup <= 0) || (nygroup > nybins)) {
2845 Error(
"Rebin",
"Illegal value of nygroup=%d",nygroup);
2848 if ((nzgroup <= 0) || (nzgroup > nzbins)) {
2849 Error(
"Rebin",
"Illegal value of nzgroup=%d",nzgroup);
2853 Int_t newxbins = nxbins/nxgroup;
2854 Int_t newybins = nybins/nygroup;
2855 Int_t newzbins = nzbins/nzgroup;
2873 if (newname && strlen(newname)) {
2881 bool resetStat =
false;
2885 if (newxbins*nxgroup != nxbins) {
2889 if (newybins*nygroup != nybins) {
2893 if (newzbins*nzgroup != nzbins) {
2935 if (nxgroup != 1 || nygroup != 1 || nzgroup != 1) {
2944 hnew->
SetBins(newxbins,xbins, newybins, ybins, newzbins, zbins);
2957 for (xbin = 1; xbin <= newxbins; xbin++) {
2960 for (ybin = 1; ybin <= newybins; ybin++) {
2962 for (zbin = 1; zbin <= newzbins; zbin++) {
2965 for (i = 0; i < nxgroup; i++) {
2966 if (oldxbin+i > nxbins)
break;
2967 for (j =0; j < nygroup; j++) {
2968 if (oldybin+j > nybins)
break;
2969 for (k =0; k < nzgroup; k++) {
2970 if (oldzbin+k > nzbins)
break;
2972 bin = oldxbin + i + (oldybin + j)*(nxbins + 2) + (oldzbin + k)*(nxbins + 2)*(nybins + 2);
2973 binContent += oldBins[bin];
2974 if (oldSumw2) binSumw2 += oldSumw2[bin];
2989 for (
Int_t xover = 0; xover <= 1; xover++) {
2990 for (
Int_t yover = 0; yover <= 1; yover++) {
2991 for (
Int_t zover = 0; zover <= 1; zover++) {
2995 for (xbin = xover*oldxbin; xbin <= xover*(nxbins+1); xbin++) {
2996 for (ybin = yover*oldybin; ybin <= yover*(nybins+1); ybin++) {
2997 for (zbin = zover*oldzbin; zbin <= zover*(nzbins+1); zbin++) {
2998 bin =
GetBin(xbin,ybin,zbin);
2999 binContent += oldBins[bin];
3000 if (oldSumw2) binSumw2 += oldSumw2[bin];
3005 yover*(newybins+1), zover*(newzbins+1) );
3012 Double_t binContent0, binContent2, binContent3, binContent4;
3013 Double_t binError0, binError2, binError3, binError4;
3014 Int_t oldxbin2, oldybin2, oldzbin2;
3015 Int_t ufbin, ofbin, ofbin2, ofbin3, ofbin4;
3021 for (xbin = 1; xbin<=newxbins; xbin++) {
3023 for (zbin = 1; zbin<=newzbins; zbin++) {
3024 binContent0 = binContent2 = 0;
3025 binError0 = binError2 = 0;
3026 for (i=0; i<nxgroup; i++) {
3027 if (oldxbin2+i > nxbins)
break;
3028 for (k=0; k<nzgroup; k++) {
3029 if (oldzbin2+k > nzbins)
break;
3031 ufbin = oldxbin2 + i + (nxbins+2)*(nybins+2)*(oldzbin2+k);
3032 binContent0 += oldBins[ufbin];
3033 if (oldSumw2) binError0 += oldSumw2[ufbin];
3034 for (ybin = oldybin; ybin <= nybins + 1; ybin++) {
3036 ofbin = ufbin + ybin*(nxbins+2);
3037 binContent2 += oldBins[ofbin];
3038 if (oldSumw2) binError2 += oldSumw2[ofbin];
3048 oldzbin2 += nzgroup;
3050 oldxbin2 += nxgroup;
3057 for (ybin = 1; ybin<=newybins; ybin++) {
3059 for (zbin = 1; zbin<=newzbins; zbin++) {
3060 binContent0 = binContent2 = 0;
3061 binError0 = binError2 = 0;
3062 for (j=0; j<nygroup; j++) {
3063 if (oldybin2+j > nybins)
break;
3064 for (k=0; k<nzgroup; k++) {
3065 if (oldzbin2+k > nzbins)
break;
3067 ufbin = (oldybin2 + j)*(nxbins+2) + (nxbins+2)*(nybins+2)*(oldzbin2+k);
3068 binContent0 += oldBins[ufbin];
3069 if (oldSumw2) binError0 += oldSumw2[ufbin];
3070 for (xbin = oldxbin; xbin <= nxbins + 1; xbin++) {
3072 ofbin = ufbin + xbin;
3073 binContent2 += oldBins[ofbin];
3074 if (oldSumw2) binError2 += oldSumw2[ofbin];
3084 oldzbin2 += nzgroup;
3086 oldybin2 += nygroup;
3093 for (xbin = 1; xbin<=newxbins; xbin++) {
3095 for (ybin = 1; ybin<=newybins; ybin++) {
3096 binContent0 = binContent2 = 0;
3097 binError0 = binError2 = 0;
3098 for (i=0; i<nxgroup; i++) {
3099 if (oldxbin2+i > nxbins)
break;
3100 for (j=0; j<nygroup; j++) {
3101 if (oldybin2+j > nybins)
break;
3103 ufbin = oldxbin2 + i + (nxbins+2)*(oldybin2+j);
3104 binContent0 += oldBins[ufbin];
3105 if (oldSumw2) binError0 += oldSumw2[ufbin];
3106 for (zbin = oldzbin; zbin <= nzbins + 1; zbin++) {
3108 ofbin = ufbin + (nxbins+2)*(nybins+2)*zbin;
3109 binContent2 += oldBins[ofbin];
3110 if (oldSumw2) binError2 += oldSumw2[ofbin];
3120 oldybin2 += nygroup;
3122 oldxbin2 += nxgroup;
3129 for (xbin = 1; xbin<=newxbins; xbin++) {
3138 for (i=0; i<nxgroup; i++) {
3139 if (oldxbin2+i > nxbins)
break;
3140 ufbin = oldxbin2 + i;
3141 binContent0 += oldBins[ufbin];
3142 if (oldSumw2) binError0 += oldSumw2[ufbin];
3143 for (ybin = oldybin; ybin <= nybins + 1; ybin++) {
3144 ofbin3 = ufbin+ybin*(nxbins+2);
3145 binContent3 += oldBins[ ofbin3 ];
3146 if (oldSumw2) binError3 += oldSumw2[ofbin3];
3147 for (zbin = oldzbin; zbin <= nzbins + 1; zbin++) {
3149 ofbin4 = oldxbin2 + i + ybin*(nxbins+2) + (nxbins+2)*(nybins+2)*zbin;
3150 binContent4 += oldBins[ofbin4];
3151 if (oldSumw2) binError4 += oldSumw2[ofbin4];
3154 for (zbin = oldzbin; zbin <= nzbins + 1; zbin++) {
3155 ofbin2 = ufbin+zbin*(nxbins+2)*(nybins+2);
3156 binContent2 += oldBins[ ofbin2 ];
3157 if (oldSumw2) binError2 += oldSumw2[ofbin2];
3163 hnew->
SetBinContent(xbin,newybins+1,newzbins+1,binContent4);
3170 oldxbin2 += nxgroup;
3177 for (zbin = 1; zbin<=newzbins; zbin++) {
3186 for (i=0; i<nzgroup; i++) {
3187 if (oldzbin2+i > nzbins)
break;
3188 ufbin = (oldzbin2 + i)*(nxbins+2)*(nybins+2);
3189 binContent0 += oldBins[ufbin];
3190 if (oldSumw2) binError0 += oldSumw2[ufbin];
3191 for (ybin = oldybin; ybin <= nybins + 1; ybin++) {
3192 ofbin3 = ufbin+ybin*(nxbins+2);
3193 binContent3 += oldBins[ ofbin3 ];
3194 if (oldSumw2) binError3 += oldSumw2[ofbin3];
3195 for (xbin = oldxbin; xbin <= nxbins + 1; xbin++) {
3197 ofbin4 = ufbin + xbin + ybin*(nxbins+2);
3198 binContent4 += oldBins[ofbin4];
3199 if (oldSumw2) binError4 += oldSumw2[ofbin4];
3202 for (xbin = oldxbin; xbin <= nxbins + 1; xbin++) {
3203 ofbin2 = xbin +(oldzbin2+i)*(nxbins+2)*(nybins+2);
3204 binContent2 += oldBins[ ofbin2 ];
3205 if (oldSumw2) binError2 += oldSumw2[ofbin2];
3211 hnew->
SetBinContent(newxbins+1,newybins+1,zbin,binContent4);
3218 oldzbin2 += nzgroup;
3225 for (ybin = 1; ybin<=newybins; ybin++) {
3234 for (i=0; i<nygroup; i++) {
3235 if (oldybin2+i > nybins)
break;
3236 ufbin = (oldybin2 + i)*(nxbins+2);
3237 binContent0 += oldBins[ufbin];
3238 if (oldSumw2) binError0 += oldSumw2[ufbin];
3239 for (xbin = oldxbin; xbin <= nxbins + 1; xbin++) {
3240 ofbin3 = ufbin+xbin;
3241 binContent3 += oldBins[ ofbin3 ];
3242 if (oldSumw2) binError3 += oldSumw2[ofbin3];
3243 for (zbin = oldzbin; zbin <= nzbins + 1; zbin++) {
3245 ofbin4 = xbin + (nxbins+2)*(nybins+2)*zbin+(oldybin2+i)*(nxbins+2);
3246 binContent4 += oldBins[ofbin4];
3247 if (oldSumw2) binError4 += oldSumw2[ofbin4];
3250 for (zbin = oldzbin; zbin <= nzbins + 1; zbin++) {
3251 ofbin2 = (oldybin2+i)*(nxbins+2)+zbin*(nxbins+2)*(nybins+2);
3252 binContent2 += oldBins[ ofbin2 ];
3253 if (oldSumw2) binError2 += oldSumw2[ofbin2];
3259 hnew->
SetBinContent(newxbins+1,ybin,newzbins+1,binContent4);
3266 oldybin2 += nygroup;
3309 if (!resetStat) hnew->
PutStats(stat);
3312 if (oldSumw2)
delete [] oldSumw2;
3343 if (bin < 0)
return;
3352void TH3::Streamer(
TBuffer &R__b)
3362 TH1::Streamer(R__b);
3363 TAtt3D::Streamer(R__b);
3405 :
TH3(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
3420 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
3433 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
3464 if (newval > -128 && newval < 128) {
fArray[bin] =
Char_t(newval);
return;}
3465 if (newval < -127)
fArray[bin] = -127;
3466 if (newval > 127)
fArray[bin] = 127;
3534void TH3C::Streamer(
TBuffer &R__b)
3547 TH1::Streamer(R__b);
3548 TArrayC::Streamer(R__b);
3550 TAtt3D::Streamer(R__b);
3552 TH3::Streamer(R__b);
3553 TArrayC::Streamer(R__b);
3666 :
TH3(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
3681 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
3694 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
3725 if (newval > -32768 && newval < 32768) {
fArray[bin] =
Short_t(newval);
return;}
3726 if (newval < -32767)
fArray[bin] = -32767;
3727 if (newval > 32767)
fArray[bin] = 32767;
3766void TH3S::Streamer(
TBuffer &R__b)
3779 TH1::Streamer(R__b);
3780 TArrayS::Streamer(R__b);
3782 TAtt3D::Streamer(R__b);
3784 TH3::Streamer(R__b);
3785 TArrayS::Streamer(R__b);
3898 :
TH3(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
3913 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
3926 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
3957 if (newval > -INT_MAX && newval < INT_MAX) {
fArray[bin] =
Int_t(newval);
return;}
3958 if (newval < -INT_MAX)
fArray[bin] = -INT_MAX;
3959 if (newval > INT_MAX)
fArray[bin] = INT_MAX;
4097 :
TH3(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
4112 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
4125 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
4176void TH3F::Streamer(
TBuffer &R__b)
4189 TH1::Streamer(R__b);
4190 TArrayF::Streamer(R__b);
4192 TAtt3D::Streamer(R__b);
4194 TH3::Streamer(R__b);
4195 TArrayF::Streamer(R__b);
4308 :
TH3(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
4323 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
4336 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
4387void TH3D::Streamer(
TBuffer &R__b)
4400 TH1::Streamer(R__b);
4401 TArrayD::Streamer(R__b);
4403 TAtt3D::Streamer(R__b);
4405 TH3::Streamer(R__b);
4406 TArrayD::Streamer(R__b);
void Warning(const char *location, const char *msgfmt,...)
Use this function in warning situations.
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 Set(Int_t n)
Set size of this array to n chars.
void Copy(TArrayC &array) const
Array of doubles (64 bits per element).
void Copy(TArrayD &array) const
void Set(Int_t n)
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)
Set size of this array to n floats.
Array of integers (32 bits per element).
void Set(Int_t n)
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)
Set size of this array to n shorts.
void Copy(TArrayS &array) const
virtual void Set(Int_t n)=0
Use this attribute class when an object should have 3D capabilities.
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.
Class to manage histogram axis.
virtual void SetBinLabel(Int_t bin, const char *label)
Set label for bin.
virtual Double_t GetBinCenter(Int_t bin) const
Return center of bin.
const TArrayD * GetXbins() const
virtual Int_t FindBin(Double_t x)
Find bin number corresponding to abscissa x.
virtual Double_t GetBinLowEdge(Int_t bin) const
Return low edge of bin.
virtual void Set(Int_t nbins, Double_t xmin, Double_t xmax)
Initialize axis with fix bins.
virtual Int_t FindFixBin(Double_t x) const
Find bin number corresponding to abscissa x.
Int_t GetLast() const
Return last bin on the axis i.e.
virtual void ImportAttributes(const TAxis *axis)
Copy axis attributes to this.
const char * GetTitle() const
Returns title of object.
virtual void SetRange(Int_t first=0, Int_t last=0)
Set the viewing range for the axis using bin numbers.
virtual Double_t GetBinWidth(Int_t bin) const
Return bin width.
virtual Double_t GetBinUpEdge(Int_t bin) const
Return up edge of bin.
Int_t GetFirst() const
Return first bin on the axis i.e.
THashList * GetLabels() const
Buffer base class used for serializing objects.
virtual Int_t ReadClassBuffer(const TClass *cl, void *pointer, const TClass *onfile_class=0)=0
TObject * GetParent() const
Return pointer to parent of this buffer.
virtual Version_t ReadVersion(UInt_t *start=0, UInt_t *bcnt=0, const TClass *cl=0)=0
virtual Int_t CheckByteCount(UInt_t startpos, UInt_t bcnt, const TClass *clss)=0
virtual Int_t GetVersionOwner() const =0
virtual Int_t WriteClassBuffer(const TClass *cl, void *pointer)=0
virtual TH1 * GetHistogram() const
Return a pointer to the histogram used to visualise the function Note that this histogram is managed ...
virtual Double_t GetParError(Int_t ipar) const
Return value of parameter number ipar.
Double_t GetChisquare() const
virtual void SetRange(Double_t xmin, Double_t xmax)
Initialize the upper and lower bounds to draw the function.
virtual Int_t GetNpar() const
virtual Int_t GetNumberFitPoints() const
virtual Double_t * GetParameters() const
virtual Double_t EvalPar(const Double_t *x, const Double_t *params=0)
Evaluate function with given coordinates and parameters.
virtual void GetRange(Double_t *xmin, Double_t *xmax) const
Return range of a generic N-D function.
virtual const char * GetParName(Int_t ipar) const
virtual void SetParameters(const Double_t *params)
virtual Double_t GetParameter(Int_t ipar) const
A 3-Dim function with parameters.
1-D histogram with a double per channel (see TH1 documentation)}
virtual void Reset(Option_t *option="")
Reset.
TH1 is the base class of all histogram classes in ROOT.
virtual void SetDirectory(TDirectory *dir)
By default when an histogram is created, it is added to the list of histogram objects in the current ...
Double_t * fBuffer
[fBufferSize] entry buffer
virtual Double_t GetEffectiveEntries() const
Number of effective entries of the histogram.
virtual Bool_t Multiply(TF1 *f1, Double_t c1=1)
Performs the operation:
Int_t fNcells
number of bins(1D), cells (2D) +U/Overflows
Double_t fTsumw
Total Sum of weights.
Double_t fTsumw2
Total Sum of squares of weights.
virtual Double_t DoIntegral(Int_t ix1, Int_t ix2, Int_t iy1, Int_t iy2, Int_t iz1, Int_t iz2, Double_t &err, Option_t *opt, Bool_t doerr=kFALSE) const
Internal function compute integral and optionally the error between the limits specified by the bin n...
Double_t fTsumwx2
Total Sum of weight*X*X.
virtual Double_t GetStdDev(Int_t axis=1) const
Returns the Standard Deviation (Sigma).
virtual Int_t GetNbinsY() const