73 :
TH1(
name,title,nbinsx,xlow,xup)
78 if (nbinsy <= 0) {
Warning(
"TH2",
"nbinsy is <=0 - set to nbinsy = 1"); nbinsy = 1; }
94 if (nbinsy <= 0) {
Warning(
"TH2",
"nbinsy is <=0 - set to nbinsy = 1"); nbinsy = 1; }
109 fTsumwy = fTsumwy2 = fTsumwxy = 0;
110 if (nbinsy <= 0) {
Warning(
"TH2",
"nbinsy is <=0 - set to nbinsy = 1"); nbinsy = 1; }
111 if (ybins) fYaxis.Set(nbinsy,ybins);
112 else fYaxis.Set(nbinsy,0,1);
113 fNcells = fNcells*(nbinsy+2);
127 if (nbinsy <= 0) {
Warning(
"TH2",
"nbinsy is <=0 - set to nbinsy = 1"); nbinsy = 1; }
144 if (nbinsy <= 0) {
Warning(
"TH2",
"nbinsy is <=0 - set to nbinsy = 1"); nbinsy = 1; }
185 if (nbentries == 0)
return 0;
186 if (nbentries < 0 && action == 0)
return 0;
190 nbentries = -nbentries;
205 for (
Int_t i=1;i<nbentries;i++) {
228 for (
Int_t i=0;i<nbentries;i++) {
229 Fill(buffer[3*i+2],buffer[3*i+3],buffer[3*i+1]);
256 nbentries = -nbentries;
294 Error(
"Fill",
"Invalid signature - do nothing");
318 Int_t binx, biny, bin;
322 if (binx <0 || biny <0)
return -1;
362 Int_t binx, biny, bin;
366 if (binx <0 || biny <0)
return -1;
406 Int_t binx, biny, bin;
410 if (binx <0 || biny <0)
return -1;
448 Int_t binx, biny, bin;
452 if (binx <0 || biny <0)
return -1;
491 Int_t binx, biny, bin;
495 if (binx <0 || biny <0)
return -1;
535 Int_t binx, biny, bin, i;
542 for (i=0;i<ntimes;i+=stride) {
555 for (i=ifirst;i<ntimes;i+=stride) {
559 if (binx <0 || biny <0)
continue;
599 Int_t bin, binx, biny, ibin, loop;
603 if (!fobj) {
Error(
"FillRandom",
"Unknown function: %s",fname);
return; }
604 TF2 *
f1 =
dynamic_cast<TF2*
>(fobj);
605 if (!
f1) {
Error(
"FillRandom",
"Function: %s is not a TF2, is a %s",fname,fobj->IsA()->
GetName());
return; }
624 Int_t nbins = nbinsx*nbinsy;
630 for (biny=1;biny<=nbinsy;biny++) {
631 for (binx=1;binx<=nbinsx;binx++) {
634 integral[ibin] = integral[ibin-1] + fint;
639 if (integral[nbins] == 0 ) {
641 Error(
"FillRandom",
"Integral = zero");
return;
643 for (bin=1;bin<=nbins;bin++) integral[bin] /= integral[nbins];
646 for (loop=0;loop<ntimes;loop++) {
650 binx = 1 + ibin - nbinsx*biny;
674 if (!
h) {
Error(
"FillRandom",
"Null histogram");
return; }
676 Error(
"FillRandom",
"Histograms with different dimensions");
return;
679 if (
h->ComputeIntegral() == 0)
return;
684 for (loop=0;loop<ntimes;loop++) {
700 if (firstbin < 0) firstbin = 0;
701 if (lastbin < 0 || lastbin > nbins + 1) lastbin = nbins + 1;
702 if (lastbin < firstbin) {firstbin = 0; lastbin = nbins + 1;}
707 if (i1>=0 && i2>i1) {
708 proj_opt += opt(i1,i2-i1+1);
719 Int_t nstep = ngroup;
729 if (npar <= 0)
return;
741 char *
name =
new char[2000];
742 char *title =
new char[2000];
744 for (ipar=0;ipar<npar;ipar++) {
755 (*arr)[ipar] = hlist[ipar];
767 (*arr)[npar] = hchi2;
774 for (bin=firstbin;bin+ngroup-1<=lastbin;bin += nstep) {
777 hp=
ProjectionX(
"_temp",bin,bin+ngroup-1,proj_opt);
779 hp=
ProjectionY(
"_temp",bin,bin+ngroup-1,proj_opt);
780 if (hp == 0)
continue;
786 if (npfits > npar && npfits >= cut) {
787 Int_t binOn = bin + ngroup/2;
788 for (ipar=0;ipar<npar;ipar++) {
860 DoFitSlices(
true,
f1, firstybin, lastybin, cut, option, arr);
925 DoFitSlices(
false,
f1, firstxbin, lastxbin, cut, option, arr);
932 if (biny < 0) biny = 0;
933 if (biny > ofy) biny = ofy;
965 Error(
"GetBinWithContent2",
"function is only valid for 2-D histograms");
968 if (firstxbin < 0) firstxbin = 1;
970 if (firstybin < 0) firstybin = 1;
973 for (
Int_t j = firstybin; j <= lastybin; j++) {
974 for (
Int_t i = firstxbin; i <= lastxbin; i++) {
976 if (diff <= 0) {binx = i; biny=j;
return diff;}
977 if (diff < curmax && diff <= maxdiff) {curmax = diff, binx=i; biny=j;}
989 if (axis1 < 1 || axis2 < 1 || axis1 > 2 || axis2 > 2) {
990 Error(
"GetCorrelationFactor",
"Wrong parameters");
993 if (axis1 == axis2)
return 1;
995 if (stddev1 == 0)
return 0;
997 if (stddev2 == 0)
return 0;
1007 if (axis1 < 1 || axis2 < 1 || axis1 > 2 || axis2 > 2) {
1008 Error(
"GetCovariance",
"Wrong parameters");
1021 if (sumw == 0)
return 0;
1022 if (axis1 == 1 && axis2 == 1) {
1023 return TMath::Abs(sumwx2/sumw - sumwx/sumw*sumwx/sumw);
1025 if (axis1 == 2 && axis2 == 2) {
1026 return TMath::Abs(sumwy2/sumw - sumwy/sumw*sumwy/sumw);
1028 return sumwxy/sumw - sumwx/sumw*sumwy/sumw;
1041 Int_t nbins = nbinsx*nbinsy;
1050 if (integral == 0 ) {
x = 0;
y = 0;
return;}
1056 Int_t biny = ibin/nbinsx;
1057 Int_t binx = ibin - nbinsx*biny;
1104 if (firstBinX == 1) firstBinX = 0;
1108 if (firstBinY == 1) firstBinY = 0;
1112 for (
Int_t biny = firstBinY; biny <= lastBinY; ++biny) {
1114 for (
Int_t binx = firstBinX; binx <= lastBinX; ++binx) {
1166 return DoIntegral(firstxbin,lastxbin,firstybin,lastybin,-1,0,err,option);
1179 return DoIntegral(firstxbin,lastxbin,firstybin,lastybin,-1,0,error,option,
kTRUE);
1187 Error(
"Interpolate",
"This function must be called with 2 arguments for a TH2");
1206 Error(
"Interpolate",
"Cannot interpolate outside histogram domain.");
1248 if(bin_x1<1) bin_x1=1;
1252 if(bin_y1<1) bin_y1=1;
1264 f = 1.0*q11/
d*(
x2-
x)*(y2-
y)+1.0*q21/
d*(
x-
x1)*(y2-
y)+1.0*q12/
d*(
x2-
x)*(
y-y1)+1.0*q22/
d*(
x-
x1)*(
y-y1);
1274 Error(
"Interpolate",
"This function must be called with 2 arguments for a TH2");
1309 if (h2 == 0)
return 0;
1321 Error(
"KolmogorovTest",
"Histograms must be 2-D\n");
1327 Error(
"KolmogorovTest",
"Number of channels in X is different, %d and %d\n",ncx1,ncx2);
1331 Error(
"KolmogorovTest",
"Number of channels in Y is different, %d and %d\n",ncy1,ncy2);
1341 if (diff1 > difprec || diff2 > difprec) {
1342 Error(
"KolmogorovTest",
"histograms with different binning along X");
1347 if (diff1 > difprec || diff2 > difprec) {
1348 Error(
"KolmogorovTest",
"histograms with different binning along Y");
1353 Int_t ibeg = 1, jbeg = 1;
1354 Int_t iend = ncx1, jend = ncy1;
1355 if (opt.
Contains(
"U")) {ibeg = 0; jbeg = 0;}
1356 if (opt.
Contains(
"O")) {iend = ncx1+1; jend = ncy1+1;}
1363 for (i = ibeg; i <= iend; i++) {
1364 for (j = jbeg; j <= jend; j++) {
1377 Error(
"KolmogorovTest",
"Integral is zero for h1=%s\n",
h1->
GetName());
1381 Error(
"KolmogorovTest",
"Integral is zero for h2=%s\n",h2->
GetName());
1389 esum1 = sum1 * sum1 / w1;
1394 esum2 = sum2 * sum2 / w2;
1398 if (afunc2 && afunc1) {
1399 Error(
"KolmogorovTest",
"Errors are zero for both histograms\n");
1408 for (i=ibeg;i<=iend;i++) {
1409 for (j=jbeg;j<=jend;j++) {
1419 for (j=jbeg;j<=jend;j++) {
1420 for (i=ibeg;i<=iend;i++) {
1431 else factnm =
TMath::Sqrt(esum1*sum2/(esum1+esum2));
1434 Double_t dfmax = 0.5*(dfmax1+dfmax2);
1441 if (opt.
Contains(
"N") && !(afunc1 || afunc2 ) ) {
1445 Double_t chi2 = d12*d12/(esum1+esum2);
1448 if (prb > 0 && prb2 > 0) prb = prb*prb2*(1-
TMath::Log(prb*prb2));
1454 printf(
" Kolmo Prob h1 = %s, sum1=%g\n",
h1->
GetName(),sum1);
1455 printf(
" Kolmo Prob h2 = %s, sum2=%g\n",h2->
GetName(),sum2);
1456 printf(
" Kolmo Probabil = %f, Max Dist = %g\n",prb,dfmax);
1458 printf(
" Kolmo Probabil = %f for shape alone, =%f for normalisation alone\n",prb1,prb2);
1464 if(opt.
Contains(
"M"))
return dfmax;
1476 return Rebin2D(ngroup, 1, newname);
1486 return Rebin2D(1, ngroup, newname);
1497 if (xbins !=
nullptr) {
1498 Error(
"Rebin",
"Rebinning a 2-d histogram into variable bins is not supported (it is possible only for 1-d histograms). Return a nullptr");
1501 Info(
"Rebin",
"Rebinning only the x-axis. Use Rebin2D for rebinning both axes");
1502 return RebinX(ngroup, newname);
1534 Int_t nx = nxbins + 2;
1535 Int_t ny = nybins + 2;
1542 Error(
"Rebin2D",
"Histogram must be TH2. This histogram has %d dimensions.",
GetDimension());
1545 if ((nxgroup <= 0) || (nxgroup > nxbins)) {
1546 Error(
"Rebin2D",
"Illegal value of nxgroup=%d",nxgroup);
1549 if ((nygroup <= 0) || (nygroup > nybins)) {
1550 Error(
"Rebin2D",
"Illegal value of nygroup=%d",nygroup);
1554 Int_t newxbins = nxbins / nxgroup;
1555 Int_t newybins = nybins / nygroup;
1556 Int_t newnx = newxbins + 2;
1557 Int_t newny = newybins + 2;
1571 if (newname && strlen(newname)) {
1576 bool resetStat =
false;
1579 if(newxbins * nxgroup != nxbins) {
1583 if(newybins * nygroup != nybins) {
1615 if (nxgroup != 1 || nygroup != 1) {
1622 hnew->
SetBins(newxbins, xbins, newybins, ybins);
1631 if (oldErrors) hnew->
fSumw2[0] = 0;
1634 for(
Int_t binx = 1, oldbinx = 1; binx < newnx; ++binx, oldbinx += nxgroup){
1635 Double_t binContent = 0.0, binErrorSq = 0.0;
1636 for (
Int_t i = 0; i < nxgroup && (oldbinx + i) < nx; ++i) {
1637 Int_t bin = oldbinx + i;
1638 binContent += oldBins[bin];
1639 if(oldErrors) binErrorSq += oldErrors[bin];
1641 Int_t newbin = binx;
1643 if (oldErrors) hnew->
fSumw2[newbin] = binErrorSq;
1647 for(
Int_t biny = 1, oldbiny = 1; biny < newny; ++biny, oldbiny += nygroup){
1648 Double_t binContent = 0.0, binErrorSq = 0.0;
1649 for (
Int_t j = 0; j < nygroup && (oldbiny + j) < ny; ++j) {
1650 Int_t bin = (oldbiny + j) * nx;
1651 binContent += oldBins[bin];
1652 if(oldErrors) binErrorSq += oldErrors[bin];
1654 Int_t newbin = biny * newnx;
1656 if (oldErrors) hnew->
fSumw2[newbin] = binErrorSq;
1660 for (
Int_t binx = 1, oldbinx = 1; binx < newnx; ++binx, oldbinx += nxgroup) {
1661 for (
Int_t biny = 1, oldbiny = 1; biny < newny; ++biny, oldbiny += nygroup) {
1662 Double_t binContent = 0.0, binErrorSq = 0.0;
1663 for (
Int_t i = 0; i < nxgroup && (oldbinx + i) < nx; ++i) {
1664 for (
Int_t j = 0; j < nygroup && (oldbiny + j) < ny; ++j) {
1665 Int_t bin = oldbinx + i + (oldbiny + j) * nx;
1666 binContent += oldBins[bin];
1667 if (oldErrors) binErrorSq += oldErrors[bin];
1670 Int_t newbin = binx + biny * newnx;
1672 if (oldErrors) hnew->
fSumw2[newbin] = binErrorSq;
1705 if (oldErrors)
delete [] oldErrors;
1720 cut = opt(i1,i2-i1+1);
1723 bool originalRange = opt.
Contains(
"o");
1728 const char *expectedName = ( onX ?
"_pfx" :
"_pfy" );
1741 if (firstbin == 0 && lastbin == 0)
1747 if (firstbin < 0) firstbin = 1;
1748 if (lastbin < 0) lastbin = inN;
1749 if (lastbin > inN+1) lastbin = inN;
1752 char *pname = (
char*)
name;
1753 if (
name && strcmp(
name, expectedName) == 0) {
1755 pname =
new char[nch];
1764 Error(
"DoProfile",
"Histogram with name %s must be a TProfile and is a %s",
name,h1obj->
ClassName());
1773 if (xbins->
fN == 0) {
1774 if ( originalRange )
1789 ((
TH2 *)
this)->GetPainter();
1795 if (bins->
fN == 0) {
1796 if ( originalRange )
1810 if (pname !=
name)
delete [] pname;
1821 TArrayD & binSumw2 = *(
h1->GetBinSumw2());
1835 for (
Int_t outbin = 0; outbin <= outAxis.
GetNbins() + 1; ++outbin) {
1841 if (binOut <0)
continue;
1843 for (
Int_t inbin = firstbin ; inbin <= lastbin ; ++inbin) {
1845 if (onX) { binx = outbin; biny=inbin; }
1846 else { binx = inbin; biny=outbin; }
1858 if ( useWeights ) tmp = binSumw2.
fArray[binOut];
1886 if (padsav) padsav->
cd();
1997 const char *expectedName = 0;
1999 const TAxis* outAxis;
2000 const TAxis* inAxis;
2007 cut = opt(i1,i2-i1+1);
2010 bool originalRange = opt.
Contains(
"o");
2014 expectedName =
"_px";
2021 expectedName =
"_py";
2038 if (firstbin == 0 && lastbin == 0)
2044 if (firstbin < 0) firstbin = 0;
2045 if (lastbin < 0) lastbin = inNbin + 1;
2046 if (lastbin > inNbin+1) lastbin = inNbin + 1;
2049 char *pname = (
char*)
name;
2050 if (
name && strcmp(
name,expectedName) == 0) {
2052 pname =
new char[nch];
2062 Error(
"DoProjection",
"Histogram with name %s must be a TH1D and is a %s",
name,h1obj->
ClassName());
2071 if (xbins->
fN == 0) {
2072 if ( originalRange )
2087 ((
TH2 *)
this)->GetPainter();
2093 if (bins->
fN == 0) {
2094 if ( originalRange )
2108 if (pname !=
name)
delete [] pname;
2136 for (
Int_t outbin = 0; outbin <= outAxis->
GetNbins() + 1; ++outbin) {
2141 for (
Int_t inbin = firstbin ; inbin <= lastbin ; ++inbin) {
2143 if (onX) { binx = outbin; biny=inbin; }
2144 else { binx = inbin; biny=outbin; }
2151 if (computeErrors) {
2165 bool reuseStats =
false;
2171 double eps = 1.E-12;
2176 if (ncuts) reuseStats =
false;
2178 bool reuseEntries = reuseStats;
2180 reuseEntries &= (firstbin==0 && lastbin == inNbin+1);
2185 stats[2] = stats[4];
2186 stats[3] = stats[5];
2221 if (padsav) padsav->
cd();
2354 const TAxis *outAxis = 0;
2363 if (qname.
IsNull() || qname ==
"_qx" || qname ==
"_qy") {
2364 const char * qtype = (onX) ?
"qx" :
"qy";
2372 h1 =
dynamic_cast<TH1D*
>(h1obj);
2374 Error(
"DoQuantiles",
"Histogram with name %s must be a TH1D and is a %s",qname.
Data(),h1obj->
ClassName());
2398 for (
int ibin = outAxis->
GetFirst() ; ibin <= outAxis->GetLast() ; ++ibin) {
2403 if (slice->
GetSum() == 0)
continue;
2417 if (slice)
delete slice;
2445 if (bin < 0)
return;
2489 return (
TH1*)
gROOT->ProcessLineFast(
Form(
"TSpectrum2::StaticBackground((TH1*)0x%lx,%d,\"%s\")",
2490 (
ULong_t)
this, niter, option));
2505 return (
Int_t)
gROOT->ProcessLineFast(
Form(
"TSpectrum2::StaticSearch((TH1*)0x%lx,%g,\"%s\",%g)",
2533 Double_t k5a[5][5] = { { 0, 0, 1, 0, 0 },
2537 { 0, 0, 1, 0, 0 } };
2538 Double_t k5b[5][5] = { { 0, 1, 2, 1, 0 },
2542 { 0, 1, 2, 1, 0 } };
2543 Double_t k3a[3][3] = { { 0, 1, 0 },
2548 Warning(
"Smooth",
"Currently only ntimes=1 is supported");
2555 if (opt.
Contains(
"k5b")) kernel = &k5b[0][0];
2557 kernel = &k3a[0][0];
2572 Int_t bufSize = (nx+2)*(ny+2);
2579 for (i=ifirst; i<=ilast; i++){
2580 for (j=jfirst; j<=jlast; j++){
2588 Int_t x_push = (ksize_x-1)/2;
2589 Int_t y_push = (ksize_y-1)/2;
2592 for (i=ifirst; i<=ilast; i++){
2593 for (j=jfirst; j<=jlast; j++) {
2602 if ( (xb >= 1) && (xb <= nx) && (yb >= 1) && (yb <= ny) ) {
2608 content += k*buf[bin];
2609 if (ebuf) error += k*k*ebuf[bin]*ebuf[bin];
2615 if ( norm != 0.0 ) {
2618 error /= (norm*norm);
2634void TH2::Streamer(
TBuffer &R__b)
2644 TH1::Streamer(R__b);
2688 :
TH2(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup)
2702 :
TH2(
name,title,nbinsx,xbins,nbinsy,ylow,yup)
2714 :
TH2(
name,title,nbinsx,xlow,xup,nbinsy,ybins)
2726 :
TH2(
name,title,nbinsx,xbins,nbinsy,ybins)
2738 :
TH2(
name,title,nbinsx,xbins,nbinsy,ybins)
2769 if (newval > -128 && newval < 128) {
fArray[bin] =
Char_t(newval);
return;}
2770 if (newval < -127)
fArray[bin] = -127;
2771 if (newval > 127)
fArray[bin] = 127;
2809void TH2C::Streamer(
TBuffer &R__b)
2821 TH1::Streamer(R__b);
2822 TArrayC::Streamer(R__b);
2829 TH2::Streamer(R__b);
2830 TArrayC::Streamer(R__b);
2942 :
TH2(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup)
2956 :
TH2(
name,title,nbinsx,xbins,nbinsy,ylow,yup)
2968 :
TH2(
name,title,nbinsx,xlow,xup,nbinsy,ybins)
2980 :
TH2(
name,title,nbinsx,xbins,nbinsy,ybins)
2992 :
TH2(
name,title,nbinsx,xbins,nbinsy,ybins)
3023 if (newval > -32768 && newval < 32768) {
fArray[bin] =
Short_t(newval);
return;}
3024 if (newval < -32767)
fArray[bin] = -32767;
3025 if (newval > 32767)
fArray[bin] = 32767;
3063void TH2S::Streamer(
TBuffer &R__b)
3075 TH1::Streamer(R__b);
3076 TArrayS::Streamer(R__b);
3083 TH2::Streamer(R__b);
3084 TArrayS::Streamer(R__b);
3196 :
TH2(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup)
3210 :
TH2(
name,title,nbinsx,xbins,nbinsy,ylow,yup)
3222 :
TH2(
name,title,nbinsx,xlow,xup,nbinsy,ybins)
3234 :
TH2(
name,title,nbinsx,xbins,nbinsy,ybins)
3246 :
TH2(
name,title,nbinsx,xbins,nbinsy,ybins)
3277 if (newval > -2147483647 && newval < 2147483647) {
fArray[bin] =
Int_t(newval);
return;}
3278 if (newval < -2147483647)
fArray[bin] = -2147483647;
3279 if (newval > 2147483647)
fArray[bin] = 2147483647;
3415 :
TH2(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup)
3429 :
TH2(
name,title,nbinsx,xbins,nbinsy,ylow,yup)
3441 :
TH2(
name,title,nbinsx,xlow,xup,nbinsy,ybins)
3453 :
TH2(
name,title,nbinsx,xbins,nbinsy,ybins)
3465 :
TH2(
name,title,nbinsx,xbins,nbinsy,ybins)
3476:
TH2(
"TMatrixFBase",
"",
m.GetNcols(),
m.GetColLwb(),1+
m.GetColUpb(),
m.GetNrows(),
m.GetRowLwb(),1+
m.GetRowUpb())
3479 Int_t ilow =
m.GetRowLwb();
3480 Int_t iup =
m.GetRowUpb();
3481 Int_t jlow =
m.GetColLwb();
3482 Int_t jup =
m.GetColUpb();
3483 for (
Int_t i=ilow;i<=iup;i++) {
3484 for (
Int_t j=jlow;j<=jup;j++) {
3534void TH2F::Streamer(
TBuffer &R__b)
3546 TH1::Streamer(R__b);
3547 TArrayF::Streamer(R__b);
3554 TH2::Streamer(R__b);
3555 TArrayF::Streamer(R__b);
3679 :
TH2(
name,title,nbinsx,xlow,xup,nbinsy,ylow,yup)
3693 :
TH2(
name,title,nbinsx,xbins,nbinsy,ylow,yup)
3705 :
TH2(
name,title,nbinsx,xlow,xup,nbinsy,ybins)
3717 :
TH2(
name,title,nbinsx,xbins,nbinsy,ybins)
3729 :
TH2(
name,title,nbinsx,xbins,nbinsy,ybins)
3740:
TH2(
"TMatrixDBase",
"",
m.GetNcols(),
m.GetColLwb(),1+
m.GetColUpb(),
m.GetNrows(),
m.GetRowLwb(),1+
m.GetRowUpb())
3743 Int_t ilow =
m.GetRowLwb();
3744 Int_t iup =
m.GetRowUpb();
3745 Int_t jlow =
m.GetColLwb();
3746 Int_t jup =
m.GetColUpb();
3747 for (
Int_t i=ilow;i<=iup;i++) {
3748 for (
Int_t j=jlow;j<=jup;j++) {
3799void TH2D::Streamer(
TBuffer &R__b)
3811 TH1::Streamer(R__b);
3812 TArrayD::Streamer(R__b);
3819 TH2::Streamer(R__b);
3820 TArrayD::Streamer(R__b);
static const double x2[5]
static const double x1[5]
void Warning(const char *location, const char *msgfmt,...)
TH2C operator-(TH2C &h1, TH2C &h2)
Operator -.
TH2C operator+(TH2C &h1, TH2C &h2)
Operator +.
TH2C operator/(TH2C &h1, TH2C &h2)
Operator /.
TH2C operator*(Float_t c1, TH2C &h1)
Operator *.
R__EXTERN TRandom * gRandom
char * Form(const char *fmt,...)
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 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 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
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 WriteClassBuffer(const TClass *cl, void *pointer)=0
virtual void SetOwner(Bool_t enable=kTRUE)
Set whether this collection is the owner (enable==true) of its content.
virtual TH1 * GetHistogram() const
Return a pointer to the histogram used to visualise the function.
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 Double_t Integral(Double_t a, Double_t b, Double_t epsrel=1.e-12)
IntegralOneDim or analytical integral.
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 void SetParameters(const Double_t *params)
virtual Double_t GetParameter(Int_t ipar) const
A 2-Dim function with parameters.
1-D histogram with a double per channel (see TH1 documentation)}
virtual void Reset(Option_t *option="")
Reset.
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 Int_t GetQuantiles(Int_t nprobSum, Double_t *q, const Double_t *probSum=0)
Compute Quantiles for this histogram Quantile x_q of a probability distribution Function F is defined...
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 Double_t GetMean(Int_t axis=1) const
For axis = 1,2 or 3 returns the mean value of the histogram along X,Y or Z axis.
virtual Int_t GetDimension() const
@ kIsNotW
Histogram is forced to be not weighted even when the histogram is filled with weighted different than...
virtual Bool_t CanExtendAllAxes() const
Returns true if all axes are extendable.
virtual void Reset(Option_t *option="")
Reset this histogram: contents, errors, etc.
TAxis * GetXaxis()
Get the behaviour adopted by the object about the statoverflows. See EStatOverflows for more informat...
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.
TObject * Clone(const char *newname=0) const
Make a complete copy of the underlying object.
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
virtual Int_t Fill(Double_t x)
Increment bin with abscissa X by 1.
virtual Double_t GetBinErrorSqUnchecked(Int_t bin) const
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 GetEntries() const
Return the current number of entries.
virtual void Copy(TObject &hnew) const
Copy this histogram structure to newth1.
virtual void Draw(Option_t *option="")
Draw this histogram with options.
virtual void ResetStats()
Reset the statistics including the number of entries and replace with values calculates from bin cont...
virtual void SetBuffer(Int_t buffersize, Option_t *option="")
Set the maximum number of entries to be kept in the buffer.
Double_t fEntries
Number of entries.
virtual void SetName(const char *name)
Change the name of this histogram.
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 void Paint(Option_t *option="")
Control routine to paint any kind of histograms.
virtual Int_t GetSumw2N() const
Bool_t GetStatOverflowsBehaviour() const
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 byte per channel (see TH1 documentation)
virtual void Reset(Option_t *option="")
Reset this histogram: contents, errors, etc.
virtual void AddBinContent(Int_t bin)
Increment bin content by 1.
virtual ~TH2C()
Destructor.
virtual void SetBinsLength(Int_t n=-1)
Set total number of bins including under/overflow Reallocate bin contents array.
TH2C & operator=(const TH2C &h1)
Operator =.
virtual void Copy(TObject &hnew) const
Copy.
2-D histogram with a double per channel (see TH1 documentation)}
virtual void SetBinsLength(Int_t n=-1)
Set total number of bins including under/overflow Reallocate bin contents array.
virtual ~TH2D()
Destructor.
virtual void Copy(TObject &hnew) const
Copy.
TH2D & operator=(const TH2D &h1)
Operator =.
2-D histogram with a float per channel (see TH1 documentation)}
TH2F & operator=(const TH2F &h1)
Operator =.
virtual ~TH2F()
Destructor.
virtual void SetBinsLength(Int_t n=-1)
Set total number of bins including under/overflow Reallocate bin contents array.
virtual void Copy(TObject &hnew) const
Copy.
2-D histogram with an int per channel (see TH1 documentation)}
virtual void Copy(TObject &hnew) const
Copy.
virtual void SetBinsLength(Int_t n=-1)
Set total number of bins including under/overflow Reallocate bin contents array.
virtual void AddBinContent(Int_t bin)
Increment bin content by 1.
virtual ~TH2I()
Destructor.
TH2I & operator=(const TH2I &h1)
Operator =.
2-D histogram with a short per channel (see TH1 documentation)
virtual void SetBinsLength(Int_t n=-1)
Set total number of bins including under/overflow Reallocate bin contents array.
TH2S & operator=(const TH2S &h1)
Operator =.
virtual ~TH2S()
Destructor.
virtual void Copy(TObject &hnew) const
Copy.
virtual void AddBinContent(Int_t bin)
Increment bin content by 1.
Service class for 2-Dim histogram classes.
virtual void PutStats(Double_t *stats)
Replace current statistics with the values in array stats.
TH1D * ProjectionY(const char *name="_py", Int_t firstxbin=0, Int_t lastxbin=-1, Option_t *option="") const
Project a 2-D histogram into a 1-D histogram along Y.
virtual Int_t BufferEmpty(Int_t action=0)
Fill histogram with all entries in the buffer.
virtual Double_t GetCorrelationFactor(Int_t axis1=1, Int_t axis2=2) const
Return correlation factor between axis1 and axis2.
virtual TProfile * DoProfile(bool onX, const char *name, Int_t firstbin, Int_t lastbin, Option_t *option) const
virtual Double_t GetBinWithContent2(Double_t c, Int_t &binx, Int_t &biny, Int_t firstxbin=1, Int_t lastxbin=-1, Int_t firstybin=1, Int_t lastybin=-1, Double_t maxdiff=0) const
compute first cell (binx,biny) in the range [firstxbin,lastxbin][firstybin,lastybin] for which diff =...
TProfile * ProfileX(const char *name="_pfx", Int_t firstybin=1, Int_t lastybin=-1, Option_t *option="") const
Project a 2-D histogram into a profile histogram along X.
TH1D * QuantilesY(Double_t prob=0.5, const char *name="_qy") const
Compute the Y distribution of quantiles in the other variable X name is the name of the returned hist...
TProfile * ProfileY(const char *name="_pfy", Int_t firstxbin=1, Int_t lastxbin=-1, Option_t *option="") const
Project a 2-D histogram into a profile histogram along Y.
virtual void Reset(Option_t *option="")
Reset this histogram: contents, errors, etc.
virtual TH1D * DoQuantiles(bool onX, const char *name, Double_t prob) const
Implementation of quantiles for x or y.
Int_t Fill(Double_t)
Invalid Fill method.
virtual TH1 * ShowBackground(Int_t niter=20, Option_t *option="same")
This function calculates the background spectrum in this histogram.
virtual Int_t ShowPeaks(Double_t sigma=2, Option_t *option="", Double_t threshold=0.05)
Interface to TSpectrum2::Search the function finds peaks in this histogram where the width is > sigma...
virtual void DoFitSlices(bool onX, TF1 *f1, Int_t firstbin, Int_t lastbin, Int_t cut, Option_t *option, TObjArray *arr)
TH1D * QuantilesX(Double_t prob=0.5, const char *name="_qx") const
Compute the X distribution of quantiles in the other variable Y name is the name of the returned hist...
virtual TH2 * RebinX(Int_t ngroup=2, const char *newname="")
Rebin only the X axis see Rebin2D.
virtual void SetShowProjectionY(Int_t nbins=1)
When the mouse is moved in a pad containing a 2-d view of this histogram a second canvas shows the pr...
virtual void FillRandom(const char *fname, Int_t ntimes=5000)
Fill histogram following distribution in function fname.
virtual TH2 * Rebin(Int_t ngroup=2, const char *newname="", const Double_t *xbins=0)
Override TH1::Rebin as TH2::RebinX Rebinning in variable binning as for TH1 is not allowed If a non-n...
virtual void GetStats(Double_t *stats) const
Fill the array stats from the contents of this histogram The array stats must be correctly dimensione...
virtual TH1D * DoProjection(bool onX, const char *name, Int_t firstbin, Int_t lastbin, Option_t *option) const
Internal (protected) method for performing projection on the X or Y axis called by ProjectionX or Pro...
virtual Double_t GetCovariance(Int_t axis1=1, Int_t axis2=2) const
Return covariance between axis1 and axis2.
virtual void GetRandom2(Double_t &x, Double_t &y)
Return 2 random numbers along axis x and y distributed according the cell-contents of a 2-dim histogr...
virtual void Smooth(Int_t ntimes=1, Option_t *option="")
Smooth bin contents of this 2-d histogram using kernel algorithms similar to the ones used in the ras...
virtual void FillN(Int_t, const Double_t *, const Double_t *, Int_t)
Fill this histogram with an array x and weights w.
TH1D * ProjectionX(const char *name="_px", Int_t firstybin=0, Int_t lastybin=-1, Option_t *option="") const
Project a 2-D histogram into a 1-D histogram along X.
virtual void FitSlicesX(TF1 *f1=0, Int_t firstybin=0, Int_t lastybin=-1, Int_t cut=0, Option_t *option="QNR", TObjArray *arr=0)
Project slices along X in case of a 2-D histogram, then fit each slice with function f1 and make a hi...
virtual Int_t GetBin(Int_t binx, Int_t biny, Int_t binz=0) const
Return Global bin number corresponding to binx,y,z.
virtual Double_t Integral(Option_t *option="") const
Return integral of bin contents.
virtual Double_t IntegralAndError(Int_t binx1, Int_t binx2, Int_t biny1, Int_t biny2, Double_t &err, Option_t *option="") const
Return integral of bin contents in range [firstxbin,lastxbin],[firstybin,lastybin] for a 2-D histogra...
virtual Double_t KolmogorovTest(const TH1 *h2, Option_t *option="") const
Statistical test of compatibility in shape between THIS histogram and h2, using Kolmogorov test.
virtual void SetShowProjectionX(Int_t nbins=1)
When the mouse is moved in a pad containing a 2-d view of this histogram a second canvas shows the pr...
virtual Double_t Interpolate(Double_t x) const
illegal for a TH2
virtual void FitSlicesY(TF1 *f1=0, Int_t firstxbin=0, Int_t lastxbin=-1, Int_t cut=0, Option_t *option="QNR", TObjArray *arr=0)
Project slices along Y in case of a 2-D histogram, then fit each slice with function f1 and make a hi...
virtual Double_t GetBinContent(Int_t bin) const
Return content of bin number bin.
virtual TH2 * Rebin2D(Int_t nxgroup=2, Int_t nygroup=2, const char *newname="")
Rebin this histogram grouping nxgroup/nygroup bins along the xaxis/yaxis together.
virtual Int_t BufferFill(Double_t x, Double_t y, Double_t w)
accumulate arguments in buffer.
virtual void SetBinContent(Int_t bin, Double_t content)
Set bin content.
virtual ~TH2()
Destructor.
virtual TH2 * RebinY(Int_t ngroup=2, const char *newname="")
Rebin only the Y axis see Rebin2D.
virtual void Copy(TObject &hnew) const
Copy.
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...
virtual void SetTitle(const char *title="")
Set the title of the TNamed.
virtual const char * GetTitle() const
Returns title of object.
virtual const char * GetName() const
Returns name of object.
virtual void Expand(Int_t newSize)
Expand or shrink the array to newSize elements.
Collectable string class.
Mother of all ROOT objects.
virtual const char * GetName() const
Returns name of object.
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.
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 void Info(const char *method, const char *msgfmt,...) const
Issue info message.
virtual Double_t Rndm()
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
TString & ReplaceAll(const TString &s1, const TString &s2)
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
virtual Int_t MakeCuts(char *cutsopt)=0
virtual Bool_t IsInside(Int_t x, Int_t y)=0
TVirtualPad is an abstract base class for the Pad and Canvas classes.
virtual TVirtualPad * cd(Int_t subpadnumber=0)=0
Double_t Gaus(Double_t x, Double_t mean=0, Double_t sigma=1, Bool_t norm=kFALSE)
Calculate a gaussian function with mean and sigma.
Short_t Max(Short_t a, Short_t 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...
Double_t QuietNaN()
Returns a quiet NaN as defined by IEEE 754
Double_t Floor(Double_t x)
Double_t Sqrt(Double_t x)
Double_t KolmogorovProb(Double_t z)
Calculates the Kolmogorov distribution function,.
Long64_t BinarySearch(Long64_t n, const T *array, T value)