97 :
TH1(
name,title,nbinsx,xlow,xup),
102 Warning(
"TH3",
"nbinsy is <=0 - set to nbinsy = 1");
106 Warning(
"TH3",
"nbinsz is <=0 - set to nbinsz = 1");
111 fNcells = (nbinsx+2)*(nbinsy+2)*(nbinsz+2);
143 if (nbinsy <= 0) {
Warning(
"TH3",
"nbinsy is <=0 - set to nbinsy = 1"); nbinsy = 1; }
144 if (nbinsz <= 0) nbinsz = 1;
149 fNcells = (nbinsx+2)*(nbinsy+2)*(nbinsz+2);
181 if (nbinsy <= 0) {
Warning(
"TH3",
"nbinsy is <=0 - set to nbinsy = 1"); nbinsy = 1; }
182 if (nbinsz <= 0) nbinsz = 1;
187 fNcells = (nbinsx+2)*(nbinsy+2)*(nbinsz+2);
230 if (!nbentries)
return 0;
233 if (action == 0)
return 0;
234 nbentries = -nbentries;
249 for (
Int_t i=1;i<nbentries;i++) {
257 if (z < zmin) zmin = z;
258 if (z > zmax) zmax = z;
277 for (
Int_t i=0;i<nbentries;i++) {
278 Fill(buffer[4*i+2],buffer[4*i+3],buffer[4*i+4],buffer[4*i+1]);
305 nbentries = -nbentries;
331 Error(
"Fill",
"Invalid signature - do nothing");
346 Int_t binx, biny, binz, bin;
351 if (binx <0 || biny <0 || binz<0)
return -1;
394 Int_t binx, biny, binz, bin;
399 if (binx <0 || biny <0 || binz<0)
return -1;
439 Int_t binx, biny, binz, bin;
444 if (binx <0 || biny <0 || binz<0)
return -1;
487 Int_t binx, biny, binz, bin;
492 if (binx <0 || biny <0 || binz<0)
return -1;
535 Int_t binx, biny, binz, bin;
540 if (binx <0 || biny <0 || binz<0)
return -1;
583 Int_t binx, biny, binz, bin;
588 if (binx <0 || biny <0 || binz<0)
return -1;
631 Int_t binx, biny, binz, bin;
636 if (binx < 0 || biny < 0 || binz < 0)
685 Int_t binx, biny, binz, bin;
690 if (binx <0 || biny <0 || binz<0)
return -1;
734 Int_t binx, biny, binz, bin;
739 if (binx <0 || biny <0 || binz<0)
return -1;
796 Int_t bin, binx, biny, binz, ibin, loop;
800 if (!fobj) {
Error(
"FillRandom",
"Unknown function: %s",fname);
return; }
801 TF3 *
f1 =
dynamic_cast<TF3*
>( fobj );
802 if (!
f1) {
Error(
"FillRandom",
"Function: %s is not a TF3, is a %s",fname,fobj->
IsA()->
GetName());
return; }
812 Info(
"FillRandom",
"Using function axis and range ([%g,%g],[%g,%g],[%g,%g])",
xmin,
xmax,
ymin,
ymax,zmin,zmax);
822 Int_t nxy = nbinsx*nbinsy;
823 Int_t nbins = nbinsx*nbinsy*nbinsz;
829 for (binz=1;binz<=nbinsz;binz++) {
831 for (biny=1;biny<=nbinsy;biny++) {
833 for (binx=1;binx<=nbinsx;binx++) {
841 integral[ibin] = integral[ibin-1] + fint;
847 if (integral[nbins] == 0 ) {
849 Error(
"FillRandom",
"Integral = zero");
return;
851 for (bin=1;bin<=nbins;bin++) integral[bin] /= integral[nbins];
856 for (loop=0;loop<ntimes;loop++) {
860 biny = (ibin - nxy*binz)/nbinsx;
861 binx = 1 + ibin - nbinsx*(biny + nbinsy*binz);
891 if (!
h) {
Error(
"FillRandom",
"Null histogram");
return; }
893 Error(
"FillRandom",
"Histograms with different dimensions");
return;
896 if (
h->ComputeIntegral() == 0)
return;
901 for (loop=0;loop<ntimes;loop++) {
902 h3->GetRandom3(
x,
y,z,rng);
946 auto computeFirstAndLastBin = [](
const TAxis & outerAxis,
Int_t &firstbin,
Int_t &lastbin) {
954 if (firstbin == 0 && lastbin == 0) {
959 if (firstbin < 0) firstbin = 0;
960 if (lastbin < 0 || lastbin > nbins + 1) lastbin = nbins + 1;
961 if (lastbin < firstbin) {firstbin = 0; lastbin = nbins + 1;}
964 computeFirstAndLastBin(
fXaxis, binminx, binmaxx);
965 computeFirstAndLastBin(
fYaxis, binminy, binmaxy);
968 auto computeAxisLimits = [](
const TAxis & outerAxis,
Int_t firstbin,
Int_t lastbin,
970 Int_t firstOutBin = std::max(firstbin,1);
972 nBins = lastOutBin-firstOutBin+1;
980 Int_t firstBinXaxis = computeAxisLimits(
fXaxis, binminx, binmaxx, nbinsX, xMin, xMax);
983 Int_t firstBinYaxis = computeAxisLimits(
fYaxis, binminy, binmaxy, nbinsY, yMin, yMax);
1000 std::vector<TH1*> hlist(npar+1);
1003 for (ipar=0;ipar<= npar;ipar++) {
1011 title =
"chisquare";
1013 if (xbins->fN == 0 && ybins->
fN == 0) {
1014 hlist[ipar] =
new TH2D(
name, title,
1016 nbinsY, yMin, yMax);
1017 }
else if (xbins->fN > 0 && ybins->
fN > 0 ) {
1018 hlist[ipar] =
new TH2D(
name, title,
1019 nbinsX, &xbins->fArray[firstBinXaxis],
1020 nbinsY, &ybins->
fArray[firstBinYaxis]);
1028 TH1 * hchi2 = hlist.back();
1031 TH1D *hpz =
nullptr;
1036 for (
Int_t biny=binminy; biny<=binmaxy; biny++) {
1037 for (
Int_t binx=binminx; binx<=binmaxx; binx++) {
1044 Info(
"FitSlicesZ",
"Slice (%d,%d) skipped, the number of entries is zero or smaller than the given cut value, n=%f",binx,biny,
nentries);
1050 int ibx,iby,ibz = 0;
1051 hlist[0]->GetBinXYZ(bin,ibx,iby,ibz);
1057 if (npfits > npar && npfits >= cut) {
1058 for (ipar=0;ipar<npar;ipar++) {
1066 Info(
"FitSlicesZ",
"Fitted slice (%d,%d) skipped, the number of fitted points is too small, n=%d",binx,biny,npfits);
1081 if (biny < 0) biny = 0;
1082 if (biny > ofy) biny = ofy;
1085 if (binz < 0) binz = 0;
1086 if (binz > ofz) binz = ofz;
1122 Error(
"GetBinWithContent3",
"function is only valid for 3-D histograms");
1125 if (firstx <= 0) firstx = 1;
1127 if (firsty <= 0) firsty = 1;
1129 if (firstz <= 0) firstz = 1;
1131 Int_t binminx = 0, binminy=0, binminz=0;
1133 for (
Int_t k=firstz;k<=lastz;k++) {
1134 for (
Int_t j=firsty;j<=lasty;j++) {
1135 for (
Int_t i=firstx;i<=lastx;i++) {
1137 if (diff <= 0) {binx = i; biny=j; binz=k;
return diff;}
1138 if (diff < curmax && diff <= maxdiff) {curmax = diff, binminx=i; binminy=j;binminz=k;}
1154 if (axis1 < 1 || axis2 < 1 || axis1 > 3 || axis2 > 3) {
1155 Error(
"GetCorrelationFactor",
"Wrong parameters");
1158 if (axis1 == axis2)
return 1;
1160 if (stddev1 == 0)
return 0;
1162 if (stddev2 == 0)
return 0;
1172 if (axis1 < 1 || axis2 < 1 || axis1 > 3 || axis2 > 3) {
1173 Error(
"GetCovariance",
"Wrong parameters");
1189 if (sumw == 0)
return 0;
1190 if (axis1 == 1 && axis2 == 1) {
1191 return TMath::Abs(sumwx2/sumw - sumwx*sumwx/(sumw*sumw));
1193 if (axis1 == 2 && axis2 == 2) {
1194 return TMath::Abs(sumwy2/sumw - sumwy*sumwy/(sumw*sumw));
1196 if (axis1 == 3 && axis2 == 3) {
1197 return TMath::Abs(sumwz2/sumw - sumwz*sumwz/(sumw*sumw));
1199 if ((axis1 == 1 && axis2 == 2) || (axis1 == 2 && axis2 == 1)) {
1200 return sumwxy/sumw - sumwx*sumwy/(sumw*sumw);
1202 if ((axis1 == 1 && axis2 == 3) || (axis1 == 3 && axis2 == 1)) {
1203 return sumwxz/sumw - sumwx*sumwz/(sumw*sumw);
1205 if ((axis1 == 2 && axis2 == 3) || (axis1 == 3 && axis2 == 2)) {
1206 return sumwyz/sumw - sumwy*sumwz/(sumw*sumw);
1225 Int_t nxy = nbinsx*nbinsy;
1226 Int_t nbins = nxy*nbinsz;
1235 if (integral == 0 ) {
x = 0;
y = 0; z = 0;
return;}
1242 Int_t binz = ibin/nxy;
1243 Int_t biny = (ibin - nxy*binz)/nbinsx;
1244 Int_t binx = ibin - nbinsx*(biny + nbinsy*binz);
1272 Int_t bin, binx, biny, binz;
1276 for (bin=0;bin<11;bin++) stats[bin] = 0;
1287 if (firstBinX == 1) firstBinX = 0;
1291 if (firstBinY == 1) firstBinY = 0;
1295 if (firstBinZ == 1) firstBinZ = 0;
1305 for (binz = firstBinZ; binz <= lastBinZ; binz++) {
1307 for (biny = firstBinY; biny <= lastBinY; biny++) {
1309 for (binx = firstBinX; binx <= lastBinX; binx++) {
1310 bin =
GetBin(binx,biny,binz);
1316 stats[1] += err*err;
1394 Error(
"Interpolate",
"This function must be called with 3 arguments for a TH3");
1404 Error(
"Interpolate",
"This function must be called with 3 arguments for a TH3");
1423 Int_t obx = ubx + 1;
1427 Int_t oby = uby + 1;
1431 Int_t obz = ubz + 1;
1436 if (ubx <=0 || uby <=0 || ubz <= 0 ||
1438 Error(
"Interpolate",
"Cannot interpolate outside histogram domain.");
1463 Double_t w1 = i1 * (1 - yd) + i2 * yd;
1464 Double_t w2 = j1 * (1 - yd) + j2 * yd;
1500 if (h2 ==
nullptr)
return 0;
1502 const TAxis *xaxis2 = h2->GetXaxis();
1504 const TAxis *yaxis2 = h2->GetYaxis();
1506 const TAxis *zaxis2 = h2->GetZaxis();
1516 Error(
"KolmogorovTest",
"Histograms must be 3-D\n");
1522 Error(
"KolmogorovTest",
"Number of channels in X is different, %d and %d\n",ncx1,ncx2);
1526 Error(
"KolmogorovTest",
"Number of channels in Y is different, %d and %d\n",ncy1,ncy2);
1530 Error(
"KolmogorovTest",
"Number of channels in Z is different, %d and %d\n",ncz1,ncz2);
1540 if (diff1 > difprec || diff2 > difprec) {
1541 Error(
"KolmogorovTest",
"histograms with different binning along X");
1546 if (diff1 > difprec || diff2 > difprec) {
1547 Error(
"KolmogorovTest",
"histograms with different binning along Y");
1552 if (diff1 > difprec || diff2 > difprec) {
1553 Error(
"KolmogorovTest",
"histograms with different binning along Z");
1558 Int_t ibeg = 1, jbeg = 1, kbeg = 1;
1559 Int_t iend = ncx1, jend = ncy1, kend = ncz1;
1560 if (opt.
Contains(
"U")) {ibeg = 0; jbeg = 0; kbeg = 0;}
1561 if (opt.
Contains(
"O")) {iend = ncx1+1; jend = ncy1+1; kend = ncz1+1;}
1568 for (i = ibeg; i <= iend; i++) {
1569 for (j = jbeg; j <= jend; j++) {
1570 for (k = kbeg; k <= kend; k++) {
1573 sum2 += h2->GetBinContent(bin);
1575 Double_t ew2 = h2->GetBinError(bin);
1585 Error(
"KolmogorovTest",
"Integral is zero for h1=%s\n",
h1->
GetName());
1589 Error(
"KolmogorovTest",
"Integral is zero for h2=%s\n",h2->GetName());
1597 esum1 = sum1 * sum1 / w1;
1602 esum2 = sum2 * sum2 / w2;
1606 if (afunc2 && afunc1) {
1607 Error(
"KolmogorovTest",
"Errors are zero for both histograms\n");
1613 int order[3] = {0,1,2};
1617 binbeg[0] = ibeg; binbeg[1] = jbeg; binbeg[2] = kbeg;
1618 binend[0] = iend; binend[1] = jend; binend[2] = kend;
1627 for (i = binbeg[order[0] ]; i <= binend[order[0] ]; i++) {
1628 for ( j = binbeg[order[1] ]; j <= binend[order[1] ]; j++) {
1629 for ( k = binbeg[order[2] ]; k <= binend[order[2] ]; k++) {
1630 ibin[ order[0] ] = i;
1631 ibin[ order[1] ] = j;
1632 ibin[ order[2] ] = k;
1633 bin =
h1->
GetBin(ibin[0],ibin[1],ibin[2]);
1635 rsum2 += s2*h2->GetBinContent(bin);
1640 vdfmax[icomb] = dmax;
1659 if (opt.
Contains(
"N") && !(afunc1 || afunc2 ) ) {
1663 Double_t chi2 = d12*d12/(esum1+esum2);
1666 if (prb > 0 && prb2 > 0) prb = prb*prb2*(1-
TMath::Log(prb*prb2));
1672 printf(
" Kolmo Prob h1 = %s, sum1=%g\n",
h1->
GetName(),sum1);
1673 printf(
" Kolmo Prob h2 = %s, sum2=%g\n",h2->GetName(),sum2);
1674 printf(
" Kolmo Probabil = %f, Max Dist = %g\n",prb,dfmax);
1676 printf(
" Kolmo Probabil = %f for shape alone, =%f for normalisation alone\n",prb1,prb2);
1680 if (
TMath::Abs(rsum2-1) > 0.002)
Warning(
"KolmogorovTest",
"Numerical problems with h2=%s\n",h2->GetName());
1682 if (opt.
Contains(
"M"))
return dfmax;
1815 computeErrors =
kTRUE;
1820 originalRange =
kTRUE;
1824 TH1D *
h1 =
DoProject1D(
name, title, projAxis, &out1, &out2, computeErrors, originalRange,
true,
true);
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 =
nullptr, *refY =
nullptr, *refZ =
nullptr;
1946 Int_t ixbin, out1bin, out2bin;
1962 R__ASSERT (refX !=
nullptr && refY !=
nullptr && refZ !=
nullptr);
1984 for (ixbin=0;ixbin<=1+projX->
GetNbins();ixbin++) {
1991 for (out1bin = out1min; out1bin <= out1max; out1bin++) {
1992 for (out2bin = out2min; out2bin <= out2max; out2bin++) {
1998 if (computeErrors) {
2017 bool resetStats =
true;
2018 double eps = 1.E-12;
2022 bool resetEntries = resetStats;
2024 resetEntries |= !useUF || !useOF;
2031 stats[2] = stats[4];
2032 stats[3] = stats[5];
2035 stats[2] = stats[7];
2036 stats[3] = stats[8];
2065 bool computeErrors,
bool originalRange,
2066 bool useUF,
bool useOF)
const
2076 Int_t nx = ixmax-ixmin+1;
2077 Int_t ny = iymax-iymin+1;
2085 Error(
"DoProject2D",
"Histogram with name %s must be a TH2D and is a %s",
name,h2obj->
ClassName());
2093 if ( originalRange ) {
2100 h2->GetYaxis()->Set(projX->
GetNbins(),&xbins->fArray[ixmin-1]);
2105 h2->GetXaxis()->Set(ny,&ybins->
fArray[iymin-1]);
2107 h2->GetYaxis()->Set(nx,&xbins->fArray[ixmin-1]);
2115 if ( originalRange )
2117 if (xbins->fN == 0 && ybins->
fN == 0) {
2120 }
else if (ybins->
fN == 0) {
2122 ,projX->
GetNbins(),&xbins->fArray[ixmin-1]);
2123 }
else if (xbins->fN == 0) {
2130 if (xbins->fN == 0 && ybins->
fN == 0) {
2133 }
else if (ybins->
fN == 0) {
2135 ,nx,&xbins->fArray[ixmin-1]);
2136 }
else if (xbins->fN == 0) {
2140 h2 =
new TH2D(
name,title,ny,&ybins->
fArray[iymin-1],nx,&xbins->fArray[ixmin-1]);
2149 h2->GetXaxis()->ImportAttributes(projY);
2150 h2->GetYaxis()->ImportAttributes(projX);
2158 h2->GetXaxis()->SetBinLabel(i,lb->
String().
Data());
2167 h2->GetYaxis()->SetBinLabel(i,lb->
String().
Data());
2177 if ( computeErrors && (h2->GetSumw2N() != h2->GetNcells()) ) h2->Sumw2();
2180 const TAxis* out =
nullptr;
2189 Int_t *refX =
nullptr, *refY =
nullptr, *refZ =
nullptr;
2190 Int_t ixbin, iybin, outbin;
2191 if ( projX ==
GetXaxis() && projY ==
GetYaxis() ) { refX = &ixbin; refY = &iybin; refZ = &outbin; }
2192 if ( projX ==
GetYaxis() && projY ==
GetXaxis() ) { refX = &iybin; refY = &ixbin; refZ = &outbin; }
2193 if ( projX ==
GetXaxis() && projY ==
GetZaxis() ) { refX = &ixbin; refY = &outbin; refZ = &iybin; }
2194 if ( projX ==
GetZaxis() && projY ==
GetXaxis() ) { refX = &iybin; refY = &outbin; refZ = &ixbin; }
2195 if ( projX ==
GetYaxis() && projY ==
GetZaxis() ) { refX = &outbin; refY = &ixbin; refZ = &iybin; }
2196 if ( projX ==
GetZaxis() && projY ==
GetYaxis() ) { refX = &outbin; refY = &iybin; refZ = &ixbin; }
2197 R__ASSERT (refX !=
nullptr && refY !=
nullptr && refZ !=
nullptr);
2203 Int_t outmin = out->GetFirst();
2204 Int_t outmax = out->GetLast();
2206 if (outmin == 0 && outmax == 0) { outmin = 1; outmax = out->GetNbins(); }
2211 for (ixbin=0;ixbin<=1+projX->
GetNbins();ixbin++) {
2215 for (iybin=0;iybin<=1+projY->
GetNbins();iybin++) {
2223 for (outbin = outmin; outbin <= outmax; outbin++) {
2229 if (computeErrors) {
2237 h2->SetBinContent(iy , ix, cont);
2238 if (computeErrors) h2->SetBinError(iy, ix,
TMath::Sqrt(err2) );
2247 bool resetStats =
true;
2248 double eps = 1.E-12;
2252 bool resetEntries = resetStats;
2254 resetEntries |= !useUF || !useOF;
2259 for (
Int_t i = 0; i <
kNstat; ++i) { oldst[i] = 0; }
2261 std::copy(oldst,oldst+
kNstat,stats);
2265 stats[4] = oldst[7];
2266 stats[5] = oldst[8];
2267 stats[6] = oldst[9];
2270 stats[2] = oldst[4];
2271 stats[3] = oldst[5];
2273 stats[4] = oldst[2];
2274 stats[5] = oldst[3];
2277 stats[4] = oldst[7];
2278 stats[5] = oldst[8];
2279 stats[6] = oldst[10];
2283 stats[2] = oldst[7];
2284 stats[3] = oldst[8];
2286 stats[4] = oldst[2];
2287 stats[5] = oldst[3];
2288 stats[6] = oldst[9];
2291 stats[4] = oldst[4];
2292 stats[5] = oldst[5];
2293 stats[6] = oldst[10];
2297 h2->PutStats(stats);
2307 Double_t entries = h2->GetEffectiveEntries();
2308 if (!computeErrors) entries =
TMath::Floor( entries + 0.5);
2309 h2->SetEntries( entries );
2374 Int_t underscore = extra_name.
Last(
'_');
2375 if (underscore > 0) {
2376 extra_name.
Remove(underscore,extra_name.
Length()-underscore);
2377 opt.
Remove(0,underscore+1);
2383 if (opt.
Contains(
"x")) { pcase = 1; ptype =
"x"; }
2384 if (opt.
Contains(
"y")) { pcase = 2; ptype =
"y"; }
2385 if (opt.
Contains(
"z")) { pcase = 3; ptype =
"z"; }
2386 if (opt.
Contains(
"xy")) { pcase = 4; ptype =
"xy"; }
2387 if (opt.
Contains(
"yx")) { pcase = 5; ptype =
"yx"; }
2388 if (opt.
Contains(
"xz")) { pcase = 6; ptype =
"xz"; }
2389 if (opt.
Contains(
"zx")) { pcase = 7; ptype =
"zx"; }
2390 if (opt.
Contains(
"yz")) { pcase = 8; ptype =
"yz"; }
2391 if (opt.
Contains(
"zy")) { pcase = 9; ptype =
"zy"; }
2394 Error(
"Project3D",
"No projection axis specified - return a NULL pointer");
2401 computeErrors =
kTRUE;
2418 originalRange =
kTRUE;
2428 if (underscore > 0) {
2433 title +=
" "; title += ptype; title +=
" projection";
2439 computeErrors, originalRange, useUF, useOF);
2445 computeErrors, originalRange, useUF, useOF);
2451 computeErrors, originalRange, useUF, useOF);
2457 computeErrors, originalRange, useUF, useOF);
2463 computeErrors, originalRange, useUF, useOF);
2469 computeErrors, originalRange, useUF, useOF);
2475 computeErrors, originalRange, useUF, useOF);
2481 computeErrors, originalRange, useUF, useOF);
2487 computeErrors, originalRange, useUF, useOF);
2518 if (useWeights && binSumw2.
fN <= 0) useWeights =
false;
2525 if (outBin <0)
return;
2527 if ( useWeights ) tmp = binSumw2.
fArray[outBin];
2528 p2->
Fill( u ,
v,
w, cont);
2538 bool originalRange,
bool useUF,
bool useOF)
const
2546 Int_t nx = ixmax-ixmin+1;
2547 Int_t ny = iymax-iymin+1;
2557 Error(
"DoProjectProfile2D",
"Histogram with name %s must be a TProfile2D and is a %s",
name,p2obj->
ClassName());
2565 if ( originalRange ) {
2586 if ( originalRange ) {
2587 if (xbins->fN == 0 && ybins->
fN == 0) {
2590 }
else if (ybins->
fN == 0) {
2592 ,projX->
GetNbins(),&xbins->fArray[ixmin-1]);
2593 }
else if (xbins->fN == 0) {
2600 if (xbins->fN == 0 && ybins->
fN == 0) {
2603 }
else if (ybins->
fN == 0) {
2605 ,nx,&xbins->fArray[ixmin-1]);
2606 }
else if (xbins->fN == 0) {
2640 const TAxis* outAxis =
nullptr;
2655 Int_t *refX =
nullptr, *refY =
nullptr, *refZ =
nullptr;
2656 Int_t ixbin, iybin, outbin;
2657 if ( projX ==
GetXaxis() && projY ==
GetYaxis() ) { refX = &ixbin; refY = &iybin; refZ = &outbin; }
2658 if ( projX ==
GetYaxis() && projY ==
GetXaxis() ) { refX = &iybin; refY = &ixbin; refZ = &outbin; }
2659 if ( projX ==
GetXaxis() && projY ==
GetZaxis() ) { refX = &ixbin; refY = &outbin; refZ = &iybin; }
2660 if ( projX ==
GetZaxis() && projY ==
GetXaxis() ) { refX = &iybin; refY = &outbin; refZ = &ixbin; }
2661 if ( projX ==
GetYaxis() && projY ==
GetZaxis() ) { refX = &outbin; refY = &ixbin; refZ = &iybin; }
2662 if ( projX ==
GetZaxis() && projY ==
GetYaxis() ) { refX = &outbin; refY = &iybin; refZ = &ixbin; }
2663 R__ASSERT (refX !=
nullptr && refY !=
nullptr && refZ !=
nullptr);
2673 if (useWeights && binSumw2.
fN <= 0) useWeights =
false;
2677 for (ixbin=0;ixbin<=1+projX->
GetNbins();ixbin++) {
2679 for ( iybin=0;iybin<=1+projY->
GetNbins();iybin++) {
2684 if (poutBin <0)
continue;
2686 for (outbin = outmin; outbin <= outmax; outbin++) {
2693 if (!cont)
continue;
2697 if ( useWeights ) tmp = binSumw2.
fArray[poutBin];
2707 bool resetStats =
true;
2718 if (!useWeights) entries =
TMath::Floor( entries + 0.5);
2772 if (opt.
Contains(
"xy")) { pcase = 4; ptype =
"xy"; }
2773 if (opt.
Contains(
"yx")) { pcase = 5; ptype =
"yx"; }
2774 if (opt.
Contains(
"xz")) { pcase = 6; ptype =
"xz"; }
2775 if (opt.
Contains(
"zx")) { pcase = 7; ptype =
"zx"; }
2776 if (opt.
Contains(
"yz")) { pcase = 8; ptype =
"yz"; }
2777 if (opt.
Contains(
"zy")) { pcase = 9; ptype =
"zy"; }
2780 Error(
"Project3D",
"No projection axis specified - return a NULL pointer");
2798 originalRange =
kTRUE;
2807 title +=
" profile "; title += ptype; title +=
" projection";
2869 return Rebin3D(ngroup, 1, 1, newname);
2879 return Rebin3D(1, ngroup, 1, newname);
2889 return Rebin3D(1, 1, ngroup, newname);
2920 Int_t i,j,k,xbin,ybin,zbin;
2930 if ((nxgroup <= 0) || (nxgroup > nxbins)) {
2931 Error(
"Rebin",
"Illegal value of nxgroup=%d",nxgroup);
2934 if ((nygroup <= 0) || (nygroup > nybins)) {
2935 Error(
"Rebin",
"Illegal value of nygroup=%d",nygroup);
2938 if ((nzgroup <= 0) || (nzgroup > nzbins)) {
2939 Error(
"Rebin",
"Illegal value of nzgroup=%d",nzgroup);
2943 Int_t newxbins = nxbins/nxgroup;
2944 Int_t newybins = nybins/nygroup;
2945 Int_t newzbins = nzbins/nzgroup;
2963 if (newname && strlen(newname)) {
2971 bool resetStat =
false;
2975 if (newxbins*nxgroup != nxbins) {
2979 if (newybins*nygroup != nybins) {
2983 if (newzbins*nzgroup != nzbins) {
3025 if (nxgroup != 1 || nygroup != 1 || nzgroup != 1) {
3034 hnew->
SetBins(newxbins,xbins, newybins, ybins, newzbins, zbins);
3047 for (xbin = 1; xbin <= newxbins; xbin++) {
3050 for (ybin = 1; ybin <= newybins; ybin++) {
3052 for (zbin = 1; zbin <= newzbins; zbin++) {
3055 for (i = 0; i < nxgroup; i++) {
3056 if (oldxbin+i > nxbins)
break;
3057 for (j =0; j < nygroup; j++) {
3058 if (oldybin+j > nybins)
break;
3059 for (k =0; k < nzgroup; k++) {
3060 if (oldzbin+k > nzbins)
break;
3062 bin = oldxbin + i + (oldybin + j)*(nxbins + 2) + (oldzbin + k)*(nxbins + 2)*(nybins + 2);
3063 binContent += oldBins[bin];
3064 if (oldSumw2) binSumw2 += oldSumw2[bin];
3079 for (
Int_t xover = 0; xover <= 1; xover++) {
3080 for (
Int_t yover = 0; yover <= 1; yover++) {
3081 for (
Int_t zover = 0; zover <= 1; zover++) {
3085 for (xbin = xover*oldxbin; xbin <= xover*(nxbins+1); xbin++) {
3086 for (ybin = yover*oldybin; ybin <= yover*(nybins+1); ybin++) {
3087 for (zbin = zover*oldzbin; zbin <= zover*(nzbins+1); zbin++) {
3088 bin =
GetBin(xbin,ybin,zbin);
3089 binContent += oldBins[bin];
3090 if (oldSumw2) binSumw2 += oldSumw2[bin];
3095 yover*(newybins+1), zover*(newzbins+1) );
3102 Double_t binContent0, binContent2, binContent3, binContent4;
3103 Double_t binError0, binError2, binError3, binError4;
3104 Int_t oldxbin2, oldybin2, oldzbin2;
3105 Int_t ufbin, ofbin, ofbin2, ofbin3, ofbin4;
3111 for (xbin = 1; xbin<=newxbins; xbin++) {
3113 for (zbin = 1; zbin<=newzbins; zbin++) {
3114 binContent0 = binContent2 = 0;
3115 binError0 = binError2 = 0;
3116 for (i=0; i<nxgroup; i++) {
3117 if (oldxbin2+i > nxbins)
break;
3118 for (k=0; k<nzgroup; k++) {
3119 if (oldzbin2+k > nzbins)
break;
3121 ufbin = oldxbin2 + i + (nxbins+2)*(nybins+2)*(oldzbin2+k);
3122 binContent0 += oldBins[ufbin];
3123 if (oldSumw2) binError0 += oldSumw2[ufbin];
3124 for (ybin = oldybin; ybin <= nybins + 1; ybin++) {
3126 ofbin = ufbin + ybin*(nxbins+2);
3127 binContent2 += oldBins[ofbin];
3128 if (oldSumw2) binError2 += oldSumw2[ofbin];
3138 oldzbin2 += nzgroup;
3140 oldxbin2 += nxgroup;
3147 for (ybin = 1; ybin<=newybins; ybin++) {
3149 for (zbin = 1; zbin<=newzbins; zbin++) {
3150 binContent0 = binContent2 = 0;
3151 binError0 = binError2 = 0;
3152 for (j=0; j<nygroup; j++) {
3153 if (oldybin2+j > nybins)
break;
3154 for (k=0; k<nzgroup; k++) {
3155 if (oldzbin2+k > nzbins)
break;
3157 ufbin = (oldybin2 + j)*(nxbins+2) + (nxbins+2)*(nybins+2)*(oldzbin2+k);
3158 binContent0 += oldBins[ufbin];
3159 if (oldSumw2) binError0 += oldSumw2[ufbin];
3160 for (xbin = oldxbin; xbin <= nxbins + 1; xbin++) {
3162 ofbin = ufbin + xbin;
3163 binContent2 += oldBins[ofbin];
3164 if (oldSumw2) binError2 += oldSumw2[ofbin];
3174 oldzbin2 += nzgroup;
3176 oldybin2 += nygroup;
3183 for (xbin = 1; xbin<=newxbins; xbin++) {
3185 for (ybin = 1; ybin<=newybins; ybin++) {
3186 binContent0 = binContent2 = 0;
3187 binError0 = binError2 = 0;
3188 for (i=0; i<nxgroup; i++) {
3189 if (oldxbin2+i > nxbins)
break;
3190 for (j=0; j<nygroup; j++) {
3191 if (oldybin2+j > nybins)
break;
3193 ufbin = oldxbin2 + i + (nxbins+2)*(oldybin2+j);
3194 binContent0 += oldBins[ufbin];
3195 if (oldSumw2) binError0 += oldSumw2[ufbin];
3196 for (zbin = oldzbin; zbin <= nzbins + 1; zbin++) {
3198 ofbin = ufbin + (nxbins+2)*(nybins+2)*zbin;
3199 binContent2 += oldBins[ofbin];
3200 if (oldSumw2) binError2 += oldSumw2[ofbin];
3210 oldybin2 += nygroup;
3212 oldxbin2 += nxgroup;
3219 for (xbin = 1; xbin<=newxbins; xbin++) {
3228 for (i=0; i<nxgroup; i++) {
3229 if (oldxbin2+i > nxbins)
break;
3230 ufbin = oldxbin2 + i;
3231 binContent0 += oldBins[ufbin];
3232 if (oldSumw2) binError0 += oldSumw2[ufbin];
3233 for (ybin = oldybin; ybin <= nybins + 1; ybin++) {
3234 ofbin3 = ufbin+ybin*(nxbins+2);
3235 binContent3 += oldBins[ ofbin3 ];
3236 if (oldSumw2) binError3 += oldSumw2[ofbin3];
3237 for (zbin = oldzbin; zbin <= nzbins + 1; zbin++) {
3239 ofbin4 = oldxbin2 + i + ybin*(nxbins+2) + (nxbins+2)*(nybins+2)*zbin;
3240 binContent4 += oldBins[ofbin4];
3241 if (oldSumw2) binError4 += oldSumw2[ofbin4];
3244 for (zbin = oldzbin; zbin <= nzbins + 1; zbin++) {
3245 ofbin2 = ufbin+zbin*(nxbins+2)*(nybins+2);
3246 binContent2 += oldBins[ ofbin2 ];
3247 if (oldSumw2) binError2 += oldSumw2[ofbin2];
3253 hnew->
SetBinContent(xbin,newybins+1,newzbins+1,binContent4);
3260 oldxbin2 += nxgroup;
3267 for (zbin = 1; zbin<=newzbins; zbin++) {
3276 for (i=0; i<nzgroup; i++) {
3277 if (oldzbin2+i > nzbins)
break;
3278 ufbin = (oldzbin2 + i)*(nxbins+2)*(nybins+2);
3279 binContent0 += oldBins[ufbin];
3280 if (oldSumw2) binError0 += oldSumw2[ufbin];
3281 for (ybin = oldybin; ybin <= nybins + 1; ybin++) {
3282 ofbin3 = ufbin+ybin*(nxbins+2);
3283 binContent3 += oldBins[ ofbin3 ];
3284 if (oldSumw2) binError3 += oldSumw2[ofbin3];
3285 for (xbin = oldxbin; xbin <= nxbins + 1; xbin++) {
3287 ofbin4 = ufbin + xbin + ybin*(nxbins+2);
3288 binContent4 += oldBins[ofbin4];
3289 if (oldSumw2) binError4 += oldSumw2[ofbin4];
3292 for (xbin = oldxbin; xbin <= nxbins + 1; xbin++) {
3293 ofbin2 = xbin +(oldzbin2+i)*(nxbins+2)*(nybins+2);
3294 binContent2 += oldBins[ ofbin2 ];
3295 if (oldSumw2) binError2 += oldSumw2[ofbin2];
3301 hnew->
SetBinContent(newxbins+1,newybins+1,zbin,binContent4);
3308 oldzbin2 += nzgroup;
3315 for (ybin = 1; ybin<=newybins; ybin++) {
3324 for (i=0; i<nygroup; i++) {
3325 if (oldybin2+i > nybins)
break;
3326 ufbin = (oldybin2 + i)*(nxbins+2);
3327 binContent0 += oldBins[ufbin];
3328 if (oldSumw2) binError0 += oldSumw2[ufbin];
3329 for (xbin = oldxbin; xbin <= nxbins + 1; xbin++) {
3330 ofbin3 = ufbin+xbin;
3331 binContent3 += oldBins[ ofbin3 ];
3332 if (oldSumw2) binError3 += oldSumw2[ofbin3];
3333 for (zbin = oldzbin; zbin <= nzbins + 1; zbin++) {
3335 ofbin4 = xbin + (nxbins+2)*(nybins+2)*zbin+(oldybin2+i)*(nxbins+2);
3336 binContent4 += oldBins[ofbin4];
3337 if (oldSumw2) binError4 += oldSumw2[ofbin4];
3340 for (zbin = oldzbin; zbin <= nzbins + 1; zbin++) {
3341 ofbin2 = (oldybin2+i)*(nxbins+2)+zbin*(nxbins+2)*(nybins+2);
3342 binContent2 += oldBins[ ofbin2 ];
3343 if (oldSumw2) binError2 += oldSumw2[ofbin2];
3349 hnew->
SetBinContent(newxbins+1,ybin,newzbins+1,binContent4);
3356 oldybin2 += nygroup;
3399 if (!resetStat) hnew->
PutStats(stat);
3402 if (oldSumw2)
delete [] oldSumw2;
3433 if (bin < 0)
return;
3466 bool originalRange,
bool useUF,
bool useOF)
3468 return h.DoProject1D(
name, title, projX,
nullptr,
nullptr, computeErrors, originalRange, useUF, useOF);
3475 bool computeErrors,
bool originalRange,
bool useUF,
bool useOF)
3477 return h.DoProject2D(
name, title, projX, projY, computeErrors, originalRange, useUF, useOF);
3513 :
TH3(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
3529 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
3543 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
3556 h3c.TH3C::Copy(*
this);
3577 if (newval > -128 && newval < 128) {
fArray[bin] =
Char_t(newval);
return;}
3578 if (newval < -127)
fArray[bin] = -127;
3579 if (newval > 127)
fArray[bin] = 127;
3683 h3c.TH3C::Copy(*
this);
3781 :
TH3(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
3797 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
3811 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
3824 h3s.TH3S::Copy(*
this);
3845 if (newval > -32768 && newval < 32768) {
fArray[bin] =
Short_t(newval);
return;}
3846 if (newval < -32767)
fArray[bin] = -32767;
3847 if (newval > 32767)
fArray[bin] = 32767;
3922 h3s.TH3S::Copy(*
this);
4020 :
TH3(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
4036 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
4050 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
4063 h3i.TH3I::Copy(*
this);
4084 if (newval > -INT_MAX && newval < INT_MAX) {
fArray[bin] =
Int_t(newval);
return;}
4085 if (newval < -INT_MAX)
fArray[bin] = -INT_MAX;
4086 if (newval > INT_MAX)
fArray[bin] = INT_MAX;
4128 h3i.TH3I::Copy(*
this);
4226 :
TH3(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
4242 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
4256 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
4269 h3l.TH3L::Copy(*
this);
4290 if (newval > -LLONG_MAX && newval < LLONG_MAX) {
fArray[bin] =
Int_t(newval);
return;}
4291 if (newval < -LLONG_MAX)
fArray[bin] = -LLONG_MAX;
4292 if (newval > LLONG_MAX)
fArray[bin] = LLONG_MAX;
4334 h3l.TH3L::Copy(*
this);
4432 :
TH3(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
4448 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
4462 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
4475 h3f.TH3F::Copy(*
this);
4550 h3f.TH3F::Copy(*
this);
4648 :
TH3(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
4664 :
TH3(
name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
4678 :
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 long64s (64 bits per element).
void Set(Int_t n) override
Set size of this array to n long64s.
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.
Class to manage histogram axis.
virtual void SetBinLabel(Int_t bin, const char *label)
Set label for bin.
Bool_t IsAlphanumeric() const
const char * GetTitle() const override
Returns title of object.
virtual Double_t GetBinCenter(Int_t bin) const
Return center of bin.
const TArrayD * GetXbins() const
void SetCanExtend(Bool_t canExtend)
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.
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 Version_t ReadVersion(UInt_t *start=nullptr, UInt_t *bcnt=nullptr, const TClass *cl=nullptr)=0
TObject * GetParent() const
Return pointer to parent of this buffer.
virtual Int_t CheckByteCount(UInt_t startpos, UInt_t bcnt, const TClass *clss)=0
virtual Int_t ReadClassBuffer(const TClass *cl, void *pointer, const TClass *onfile_class=nullptr)=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
Return the Chisquare after fitting. See ROOT::Fit::FitResult::Chi2()
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 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 Double_t EvalPar(const Double_t *x, const Double_t *params=nullptr)
Evaluate function with given coordinates and parameters.
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)
void Reset(Option_t *option="") override
Reset.
TH1 is the base class of all histogram classes in ROOT.
virtual void SetDirectory(TDirectory *dir)
By default, when a 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.
void Copy(TObject &hnew) const override
Copy this histogram structure to newth1.
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
virtual void AddBinContent(Int_t bin)
Increment bin content by 1.
virtual Double_t GetBinError(Int_t bin) const
Return value of error associated to bin number bin.
virtual Int_t GetNbinsZ() const
virtual Int_t GetDimension() const
void Streamer(TBuffer &) override
Stream a class object.
@ kIsNotW
Histogram is forced to be not weighted even when the histogram is filled with weighted.
virtual Bool_t CanExtendAllAxes() const
Returns true if all axes are extendable.
virtual void Reset(Option_t *option="")
Reset this histogram: contents, errors, etc.
virtual Int_t GetNcells() const
virtual void PutStats(Double_t *stats)
Replace current statistics with the values in array stats.
TVirtualHistPainter * GetPainter(Option_t *option="")
Return pointer to painter.
virtual TFitResultPtr Fit(const char *formula, Option_t *option="", Option_t *goption="", Double_t xmin=0, Double_t xmax=0)
Fit histogram with function fname.
virtual Int_t GetBin(Int_t binx, Int_t biny=0, Int_t binz=0) const
Return Global bin number corresponding to binx,y,z.
virtual Int_t GetNbinsX() const
virtual Bool_t Add(TF1 *h1, Double_t c1=1, Option_t *option="")
Performs the operation: this = this + c1*f1 if errors are defined (see TH1::Sumw2),...
Int_t fBufferSize
fBuffer size
virtual Double_t RetrieveBinContent(Int_t bin) const
Raw retrieval of bin content on internal data structure see convention for numbering bins in TH1::Get...
Int_t fDimension
! Histogram dimension (1, 2 or 3 dim)
virtual void SetBinError(Int_t bin, Double_t error)
Set the bin Error Note that this resets the bin eror option to be of Normal Type and for the non-empt...
static Int_t fgBufferSize
! Default buffer size for automatic histograms
void Draw(Option_t *option="") override
Draw this histogram with options.
UInt_t GetAxisLabelStatus() const
Internal function used in TH1::Fill to see which axis is full alphanumeric, i.e.
Double_t * fIntegral
! Integral of bins used by GetRandom
virtual void SetBinContent(Int_t bin, Double_t content)
Set bin content see convention for numbering bins in TH1::GetBin In case the bin number is greater th...
virtual Double_t GetBinLowEdge(Int_t bin) const
Return bin lower edge for 1D histogram.
virtual Double_t GetEntries() const
Return the current number of entries.
void SetName(const char *name) override
Change the name of this histogram.
void Paint(Option_t *option="") override
Control routine to paint any kind of histograms.
virtual void ResetStats()
Reset the statistics including the number of entries and replace with values calculated from bin cont...
virtual void SetBuffer(Int_t buffersize, Option_t *option="")
Set the maximum number of entries to be kept in the buffer.
@ kNstat
Size of statistics data (up to TProfile3D)
Double_t fEntries
Number of entries.
TAxis fZaxis
Z axis descriptor.
virtual void UpdateBinContent(Int_t bin, Double_t content)
Raw update of bin content on internal data structure see convention for numbering bins in TH1::GetBin...
virtual Double_t GetBinContent(Int_t bin) const
Return content of bin number bin.
TAxis fXaxis
X axis descriptor.
virtual void ExtendAxis(Double_t x, TAxis *axis)
Histogram is resized along axis such that x is in the axis range.
TArrayD fSumw2
Array of sum of squares of weights.
virtual void Scale(Double_t c1=1, Option_t *option="")
Multiply this histogram by a constant c1.
virtual Int_t GetSumw2N() const
virtual Int_t FindBin(Double_t x, Double_t y=0, Double_t z=0)
Return Global bin number corresponding to x,y,z.
Bool_t GetStatOverflowsBehaviour() const
TObject * Clone(const char *newname="") const override
Make a complete copy of the underlying object.
virtual Bool_t Divide(TF1 *f1, Double_t c1=1)
Performs the operation: this = this/(c1*f1) if errors are defined (see TH1::Sumw2),...
TAxis fYaxis
Y axis descriptor.
TVirtualHistPainter * fPainter
! Pointer to histogram painter
virtual void SetBins(Int_t nx, Double_t xmin, Double_t xmax)
Redefine x axis parameters.
virtual void Sumw2(Bool_t flag=kTRUE)
Create structure to store sum of squares of weights.
virtual void SetEntries(Double_t n)
static Bool_t fgDefaultSumw2
! Flag to call TH1::Sumw2 automatically at histogram creation time
Double_t fTsumwx
Total Sum of weight*X.
virtual Double_t ComputeIntegral(Bool_t onlyPositive=false)
Compute integral (cumulative sum of bins) The result stored in fIntegral is used by the GetRandom fun...
2-D histogram with a double per channel (see TH1 documentation)
3-D histogram with a byte per channel (see TH1 documentation)
void SetBinsLength(Int_t n=-1) override
Set total number of bins including under/overflow Reallocate bin contents array.
TClass * IsA() const override
~TH3C() override
Destructor.
void Reset(Option_t *option="") override
Reset this histogram: contents, errors, etc.
void AddBinContent(Int_t bin) override
Increment bin content by 1.
TH3C & operator=(const TH3C &h1)
Operator =.
void Streamer(TBuffer &) override
Stream an object of class TH3C.
void Copy(TObject &hnew) const override
Copy this 3-D histogram structure to newth3.
3-D histogram with a double per channel (see TH1 documentation)
TClass * IsA() const override
void Streamer(TBuffer &) override
Stream an object of class TH3D.
void SetBinsLength(Int_t n=-1) override
Set total number of bins including under/overflow Reallocate bin contents array.
~TH3D() override
Destructor.
void Copy(TObject &hnew) const override
Copy this 3-D histogram structure to newth3.
TH3D & operator=(const TH3D &h1)
Operator =.
3-D histogram with a float per channel (see TH1 documentation)
TH3F & operator=(const TH3F &h1)
Operator =.
void SetBinsLength(Int_t n=-1) override
Set total number of bins including under/overflow Reallocate bin contents array.
~TH3F() override
Destructor.
void Streamer(TBuffer &) override
Stream an object of class TH3F.
void Copy(TObject &hnew) const override
Copy this 3-D histogram structure to newth3.
TClass * IsA() const override
3-D histogram with an int per channel (see TH1 documentation)
TH3I & operator=(const TH3I &h1)
Operator =.
void AddBinContent(Int_t bin) override
Increment bin content by 1.
void Copy(TObject &hnew) const override
Copy this 3-D histogram structure to newth3.
void SetBinsLength(Int_t n=-1) override
Set total number of bins including under/overflow Reallocate bin contents array.
~TH3I() override
Destructor.
3-D histogram with a long64 per channel (see TH1 documentation)
void Copy(TObject &hnew) const override
Copy this 3-D histogram structure to newth3.
TH3L & operator=(const TH3L &h1)
Operator =.
void AddBinContent(Int_t bin) override
Increment bin content by 1.
void SetBinsLength(Int_t n=-1) override
Set total number of bins including under/overflow Reallocate bin contents array.
~TH3L() override
Destructor.
3-D histogram with a short per channel (see TH1 documentation)
void Streamer(TBuffer &) override
Stream an object of class TH3S.
void SetBinsLength(Int_t n=-1) override
Set total number of bins including under/overflow Reallocate bin contents array.
void AddBinContent(Int_t bin) override
Increment bin content by 1.
~TH3S() override
Destructor.
void Copy(TObject &hnew) const override
Copy this 3-D histogram structure to newth3.
TClass * IsA() const override
TH3S & operator=(const TH3S &h1)
Operator =.
The 3-D histogram classes derived from the 1-D histogram classes.
virtual TH3 * Rebin3D(Int_t nxgroup=2, Int_t nygroup=2, Int_t nzgroup=2, const char *newname="")
Rebin this histogram grouping nxgroup/nygroup/nzgroup bins along the xaxis/yaxis/zaxis together.
Int_t BufferEmpty(Int_t action=0) override
Fill histogram with all entries in the buffer.
Double_t fTsumwy
Total Sum of weight*Y.
Double_t fTsumwy2
Total Sum of weight*Y*Y.
virtual Double_t GetCovariance(Int_t axis1=1, Int_t axis2=2) const
Return covariance between axis1 and axis2.
void GetStats(Double_t *stats) const override
Fill the array stats from the contents of this histogram The array stats must be correctly dimensione...
void Copy(TObject &hnew) const override
Copy.
virtual TH2D * DoProject2D(const char *name, const char *title, const TAxis *projX, const TAxis *projY, bool computeErrors, bool originalRange, bool useUF, bool useOF) const
internal method performing the projection to a 2D histogram called from TH3::Project3D
Double_t fTsumwxz
Total Sum of weight*X*Z.
void FillRandom(const char *fname, Int_t ntimes=5000, TRandom *rng=nullptr) override
Fill histogram following distribution in function fname.
Double_t KolmogorovTest(const TH1 *h2, Option_t *option="") const override
Statistical test of compatibility in shape between THIS histogram and h2, using Kolmogorov test.
virtual TH1D * ProjectionY(const char *name="_py", Int_t ixmin=0, Int_t ixmax=-1, Int_t izmin=0, Int_t izmax=-1, Option_t *option="") const
Project a 3-D histogram into a 1-D histogram along Y.
Double_t Interpolate(Double_t x, Double_t y) const override
Not yet implemented.
virtual void GetRandom3(Double_t &x, Double_t &y, Double_t &, TRandom *rng=nullptr)
Return 3 random numbers along axis x , y and z distributed according to the cell-contents of this 3-d...
void Reset(Option_t *option="") override
Reset this histogram: contents, errors, etc.
~TH3() override
Destructor.
Int_t Fill(Double_t) override
Invalid Fill method.
virtual TH3 * RebinY(Int_t ngroup=2, const char *newname="")
Rebin only the Y axis see Rebin3D.
virtual Double_t IntegralAndError(Int_t binx1, Int_t binx2, Int_t biny1, Int_t biny2, Int_t binz1, Int_t binz2, Double_t &err, Option_t *option="") const
Return integral of bin contents in range [binx1,binx2],[biny1,biny2],[binz1,binz2] for a 3-D histogra...
virtual TH1D * ProjectionZ(const char *name="_pz", Int_t ixmin=0, Int_t ixmax=-1, Int_t iymin=0, Int_t iymax=-1, Option_t *option="") const
Project a 3-D histogram into a 1-D histogram along Z.
virtual TH1D * ProjectionX(const char *name="_px", Int_t iymin=0, Int_t iymax=-1, Int_t izmin=0, Int_t izmax=-1, Option_t *option="") const
Project a 3-D histogram into a 1-D histogram along X.
virtual TProfile2D * Project3DProfile(Option_t *option="xy") const
Project a 3-d histogram into a 2-d profile histograms depending on the option parameter option may co...
Double_t fTsumwz2
Total Sum of weight*Z*Z.
Double_t fTsumwxy
Total Sum of weight*X*Y.
virtual TH1 * Project3D(Option_t *option="x") const
Project a 3-d histogram into 1 or 2-d histograms depending on the option parameter,...
virtual Double_t GetBinWithContent3(Double_t c, Int_t &binx, Int_t &biny, Int_t &binz, Int_t firstx=0, Int_t lastx=0, Int_t firsty=0, Int_t lasty=0, Int_t firstz=0, Int_t lastz=0, Double_t maxdiff=0) const
Compute first cell (binx,biny,binz) in the range [firstx,lastx](firsty,lasty][firstz,...
void DoFillProfileProjection(TProfile2D *p2, const TAxis &a1, const TAxis &a2, const TAxis &a3, Int_t bin1, Int_t bin2, Int_t bin3, Int_t inBin, Bool_t useWeights) const
internal function to fill the bins of the projected profile 2D histogram called from DoProjectProfile...
virtual TH3 * RebinZ(Int_t ngroup=2, const char *newname="")
Rebin only the Z axis see Rebin3D.
void Streamer(TBuffer &) override
Stream an object of class TH3.
Double_t Integral(Option_t *option="") const override
Return integral of bin contents.
virtual Int_t BufferFill(Double_t x, Double_t y, Double_t z, Double_t w)
Accumulate arguments in buffer.
TClass * IsA() const override
virtual void SetShowProjection(const char *option="xy", Int_t nbins=1)
When the mouse is moved in a pad containing a 3-d view of this histogram a second canvas shows a proj...
TH3 * RebinX(Int_t ngroup=2, const char *newname="") override
Rebin only the X axis see Rebin3D.
virtual Double_t GetCorrelationFactor(Int_t axis1=1, Int_t axis2=2) const
Return correlation factor between axis1 and axis2.
virtual TH1D * DoProject1D(const char *name, const char *title, int imin1, int imax1, int imin2, int imax2, const TAxis *projAxis, const TAxis *axis1, const TAxis *axis2, Option_t *option) const
internal method performing the projection to 1D histogram called from TH3::Project3D
Double_t fTsumwz
Total Sum of weight*Z.
Double_t GetBinContent(Int_t binx, Int_t biny, Int_t binz) const override
void SetBinContent(Int_t bin, Double_t content) override
Set bin content.
Double_t fTsumwyz
Total Sum of weight*Y*Z.
TH3()
Default constructor.
Int_t GetBin(Int_t binx, Int_t biny, Int_t binz) const override
See comments in TH1::GetBin.
virtual void FitSlicesZ(TF1 *f1=nullptr, Int_t binminx=1, Int_t binmaxx=0, Int_t binminy=1, Int_t binmaxy=0, Int_t cut=0, Option_t *option="QNR")
Project slices along Z in case of a 3-D histogram, then fit each slice with function f1 and make a 2-...
virtual TProfile2D * DoProjectProfile2D(const char *name, const char *title, const TAxis *projX, const TAxis *projY, bool originalRange, bool useUF, bool useOF) const
internal method to project to a 2D Profile called from TH3::Project3DProfile
void PutStats(Double_t *stats) override
Replace current statistics with the values in array stats.
static THLimitsFinder * GetLimitsFinder()
Return pointer to the current finder.
virtual Int_t FindGoodLimits(TH1 *h, Double_t xmin, Double_t xmax)
Compute the best axis limits for the X axis.
THashList implements a hybrid collection class consisting of a hash table and a list to store TObject...
const char * GetName() const override
Returns name of object.
const char * GetTitle() const override
Returns title of object.
Collectable string class.
Mother of all ROOT objects.
R__ALWAYS_INLINE Bool_t TestBit(UInt_t f) const
virtual const char * ClassName() const
Returns name of class to which the object belongs.
virtual void Warning(const char *method, const char *msgfmt,...) const
Issue warning message.
virtual TObject * FindObject(const char *name) const
Must be redefined in derived classes.
void SetBit(UInt_t f, Bool_t set)
Set or unset the user status bits as specified in f.
virtual Bool_t InheritsFrom(const char *classname) const
Returns kTRUE if object inherits from class "classname".
virtual void Error(const char *method, const char *msgfmt,...) const
Issue error message.
virtual TClass * IsA() const
virtual void Info(const char *method, const char *msgfmt,...) const
Issue info message.
Profile2D histograms are used to display the mean value of Z and its error for each cell in X,...
void PutStats(Double_t *stats) override
Replace current statistics with the values in array stats.
Int_t Fill(const Double_t *v)
void Sumw2(Bool_t flag=kTRUE) override
Create/Delete structure to store sum of squares of weights per bin.
void SetBins(const Int_t *nbins, const Double_t *range)
virtual TArrayD * GetBinSumw2()
void Reset(Option_t *option="") override
Reset contents of a Profile2D histogram.
This is the base class for the ROOT Random number generators.
Double_t Rndm() override
Machine independent random number generator.
void ToLower()
Change string to lower-case.
Ssiz_t First(char c) const
Find first occurrence of a character c.
const char * Data() const
Ssiz_t Last(char c) const
Find last occurrence of a character c.
void ToUpper()
Change string to upper case.
TString & Remove(Ssiz_t pos)
static TString Format(const char *fmt,...)
Static method which formats a string using a printf style format descriptor and return a TString.
Bool_t Contains(const char *pat, ECaseCompare cmp=kExact) const
Ssiz_t Index(const char *pat, Ssiz_t i=0, ECaseCompare cmp=kExact) const
virtual void SetShowProjection(const char *option, Int_t nbins)=0
small helper class to store/restore gPad context in TPad methods
Short_t Max(Short_t a, Short_t b)
Returns the largest of a and b.
Double_t Prob(Double_t chi2, Int_t ndf)
Computation of the probability for a certain Chi-squared (chi2) and number of degrees of freedom (ndf...
Bool_t Permute(Int_t n, Int_t *a)
Simple recursive algorithm to find the permutations of n natural numbers, not necessarily all distinc...
Double_t QuietNaN()
Returns a quiet NaN as defined by IEEE 754.
Double_t Floor(Double_t x)
Rounds x downward, returning the largest integral value that is not greater than x.
Double_t Log(Double_t x)
Returns the natural logarithm of x.
Double_t Sqrt(Double_t x)
Returns the square root of x.
Double_t Mean(Long64_t n, const T *a, const Double_t *w=nullptr)
Returns the weighted mean of an array a with length n.
Double_t KolmogorovProb(Double_t z)
Calculates the Kolmogorov distribution function,.
Long64_t BinarySearch(Long64_t n, const T *array, T value)
Binary search in an array of n values to locate value.
Short_t Abs(Short_t d)
Returns the absolute value of parameter Short_t d.