#include "TROOT.h"
#include "TClass.h"
#include "TH3.h"
#include "TH2.h"
#include "TF1.h"
#include "TVirtualPad.h"
#include "THLimitsFinder.h"
#include "TRandom.h"
#include "TFile.h"
#include "TError.h"
#include "TObjString.h"
ClassImp(TH3)
TH3::TH3()
{
fDimension = 3;
fTsumwy = fTsumwy2 = fTsumwxy = 0;
fTsumwz = fTsumwz2 = fTsumwxz = fTsumwyz = 0;
SetBinsLength(27);
if (fgDefaultSumw2) Sumw2();
}
TH3::TH3(const char *name,const char *title,Int_t nbinsx,Double_t xlow,Double_t xup
,Int_t nbinsy,Double_t ylow,Double_t yup
,Int_t nbinsz,Double_t zlow,Double_t zup)
:TH1(name,title,nbinsx,xlow,xup),
TAtt3D()
{
fDimension = 3;
if (nbinsy <= 0) nbinsy = 1;
if (nbinsz <= 0) nbinsz = 1;
fYaxis.Set(nbinsy,ylow,yup);
fZaxis.Set(nbinsz,zlow,zup);
fNcells = (nbinsx+2)*(nbinsy+2)*(nbinsz+2);
fTsumwy = fTsumwy2 = fTsumwxy = 0;
fTsumwz = fTsumwz2 = fTsumwxz = fTsumwyz = 0;
if (fgDefaultSumw2) Sumw2();
}
TH3::TH3(const char *name,const char *title,Int_t nbinsx,const Float_t *xbins
,Int_t nbinsy,const Float_t *ybins
,Int_t nbinsz,const Float_t *zbins)
:TH1(name,title,nbinsx,xbins),
TAtt3D()
{
fDimension = 3;
if (nbinsy <= 0) nbinsy = 1;
if (nbinsz <= 0) nbinsz = 1;
if (ybins) fYaxis.Set(nbinsy,ybins);
else fYaxis.Set(nbinsy,0,1);
if (zbins) fZaxis.Set(nbinsz,zbins);
else fZaxis.Set(nbinsz,0,1);
fNcells = (nbinsx+2)*(nbinsy+2)*(nbinsz+2);
fTsumwy = fTsumwy2 = fTsumwxy = 0;
fTsumwz = fTsumwz2 = fTsumwxz = fTsumwyz = 0;
if (fgDefaultSumw2) Sumw2();
}
TH3::TH3(const char *name,const char *title,Int_t nbinsx,const Double_t *xbins
,Int_t nbinsy,const Double_t *ybins
,Int_t nbinsz,const Double_t *zbins)
:TH1(name,title,nbinsx,xbins),
TAtt3D()
{
fDimension = 3;
if (nbinsy <= 0) nbinsy = 1;
if (nbinsz <= 0) nbinsz = 1;
if (ybins) fYaxis.Set(nbinsy,ybins);
else fYaxis.Set(nbinsy,0,1);
if (zbins) fZaxis.Set(nbinsz,zbins);
else fZaxis.Set(nbinsz,0,1);
fNcells = (nbinsx+2)*(nbinsy+2)*(nbinsz+2);
fTsumwy = fTsumwy2 = fTsumwxy = 0;
fTsumwz = fTsumwz2 = fTsumwxz = fTsumwyz = 0;
if (fgDefaultSumw2) Sumw2();
}
TH3::TH3(const TH3 &h) : TH1(), TAtt3D()
{
Copy((TObject&)h);
}
TH3::~TH3()
{
}
void TH3::Copy(TObject &obj) const
{
TH1::Copy(obj);
((TH3&)obj).fTsumwy = fTsumwy;
((TH3&)obj).fTsumwy2 = fTsumwy2;
((TH3&)obj).fTsumwxy = fTsumwxy;
((TH3&)obj).fTsumwz = fTsumwz;
((TH3&)obj).fTsumwz2 = fTsumwz2;
((TH3&)obj).fTsumwxz = fTsumwxz;
((TH3&)obj).fTsumwyz = fTsumwyz;
}
Int_t TH3::BufferEmpty(Int_t action)
{
if (!fBuffer) return 0;
Int_t nbentries = (Int_t)fBuffer[0];
if (!nbentries) return 0;
Double_t *buffer = fBuffer;
if (nbentries < 0) {
if (action == 0) return 0;
nbentries = -nbentries;
fBuffer=0;
Reset();
fBuffer = buffer;
}
if (TestBit(kCanRebin) || fXaxis.GetXmax() <= fXaxis.GetXmin() ||
fYaxis.GetXmax() <= fYaxis.GetXmin() ||
fZaxis.GetXmax() <= fZaxis.GetXmin()) {
Double_t xmin = fBuffer[2];
Double_t xmax = xmin;
Double_t ymin = fBuffer[3];
Double_t ymax = ymin;
Double_t zmin = fBuffer[4];
Double_t zmax = zmin;
for (Int_t i=1;i<nbentries;i++) {
Double_t x = fBuffer[4*i+2];
if (x < xmin) xmin = x;
if (x > xmax) xmax = x;
Double_t y = fBuffer[4*i+3];
if (y < ymin) ymin = y;
if (y > ymax) ymax = y;
Double_t z = fBuffer[4*i+4];
if (z < zmin) zmin = z;
if (z > zmax) zmax = z;
}
if (fXaxis.GetXmax() <= fXaxis.GetXmin() || fYaxis.GetXmax() <= fYaxis.GetXmin() || fZaxis.GetXmax() <= fZaxis.GetXmin()) {
THLimitsFinder::GetLimitsFinder()->FindGoodLimits(this,xmin,xmax,ymin,ymax,zmin,zmax);
} else {
fBuffer = 0;
Int_t keep = fBufferSize; fBufferSize = 0;
if (xmin < fXaxis.GetXmin()) RebinAxis(xmin,"X");
if (xmax >= fXaxis.GetXmax()) RebinAxis(xmax,"X");
if (ymin < fYaxis.GetXmin()) RebinAxis(ymin,"Y");
if (ymax >= fYaxis.GetXmax()) RebinAxis(ymax,"Y");
if (zmin < fZaxis.GetXmin()) RebinAxis(zmin,"Z");
if (zmax >= fZaxis.GetXmax()) RebinAxis(zmax,"Z");
fBuffer = buffer;
fBufferSize = keep;
}
}
fBuffer = 0;
for (Int_t i=0;i<nbentries;i++) {
Fill(buffer[4*i+2],buffer[4*i+3],buffer[4*i+4],buffer[4*i+1]);
}
fBuffer = buffer;
if (action > 0) { delete [] fBuffer; fBuffer = 0; fBufferSize = 0;}
else {
if (nbentries == (Int_t)fEntries) fBuffer[0] = -nbentries;
else fBuffer[0] = 0;
}
return nbentries;
}
Int_t TH3::BufferFill(Double_t x, Double_t y, Double_t z, Double_t w)
{
if (!fBuffer) return -3;
Int_t nbentries = (Int_t)fBuffer[0];
if (nbentries < 0) {
nbentries = -nbentries;
fBuffer[0] = nbentries;
if (fEntries > 0) {
Double_t *buffer = fBuffer; fBuffer=0;
Reset();
fBuffer = buffer;
}
}
if (4*nbentries+4 >= fBufferSize) {
BufferEmpty(1);
return Fill(x,y,z,w);
}
fBuffer[4*nbentries+1] = w;
fBuffer[4*nbentries+2] = x;
fBuffer[4*nbentries+3] = y;
fBuffer[4*nbentries+4] = z;
fBuffer[0] += 1;
return -3;
}
Int_t TH3::Fill(Double_t x, Double_t y, Double_t z)
{
if (fBuffer) return BufferFill(x,y,z,1);
Int_t binx, biny, binz, bin;
fEntries++;
binx = fXaxis.FindBin(x);
biny = fYaxis.FindBin(y);
binz = fZaxis.FindBin(z);
bin = binx + (fXaxis.GetNbins()+2)*(biny + (fYaxis.GetNbins()+2)*binz);
AddBinContent(bin);
if (fSumw2.fN) ++fSumw2.fArray[bin];
if (binx == 0 || binx > fXaxis.GetNbins()) {
if (!fgStatOverflows) return -1;
}
if (biny == 0 || biny > fYaxis.GetNbins()) {
if (!fgStatOverflows) return -1;
}
if (binz == 0 || binz > fZaxis.GetNbins()) {
if (!fgStatOverflows) return -1;
}
++fTsumw;
++fTsumw2;
fTsumwx += x;
fTsumwx2 += x*x;
fTsumwy += y;
fTsumwy2 += y*y;
fTsumwxy += x*y;
fTsumwz += z;
fTsumwz2 += z*z;
fTsumwxz += x*z;
fTsumwyz += y*z;
return bin;
}
Int_t TH3::Fill(Double_t x, Double_t y, Double_t z, Double_t w)
{
if (fBuffer) return BufferFill(x,y,z,w);
Int_t binx, biny, binz, bin;
fEntries++;
binx = fXaxis.FindBin(x);
biny = fYaxis.FindBin(y);
binz = fZaxis.FindBin(z);
bin = binx + (fXaxis.GetNbins()+2)*(biny + (fYaxis.GetNbins()+2)*binz);
AddBinContent(bin,w);
if (fSumw2.fN) fSumw2.fArray[bin] += w*w;
if (binx == 0 || binx > fXaxis.GetNbins()) {
if (!fgStatOverflows) return -1;
}
if (biny == 0 || biny > fYaxis.GetNbins()) {
if (!fgStatOverflows) return -1;
}
if (binz == 0 || binz > fZaxis.GetNbins()) {
if (!fgStatOverflows) return -1;
}
fTsumw += w;
fTsumw2 += w*w;
fTsumwx += w*x;
fTsumwx2 += w*x*x;
fTsumwy += w*y;
fTsumwy2 += w*y*y;
fTsumwxy += w*x*y;
fTsumwz += w*z;
fTsumwz2 += w*z*z;
fTsumwxz += w*x*z;
fTsumwyz += w*y*z;
return bin;
}
Int_t TH3::Fill(const char *namex, const char *namey, const char *namez, Double_t w)
{
Int_t binx, biny, binz, bin;
fEntries++;
binx = fXaxis.FindBin(namex);
biny = fYaxis.FindBin(namey);
binz = fZaxis.FindBin(namez);
bin = binx + (fXaxis.GetNbins()+2)*(biny + (fYaxis.GetNbins()+2)*binz);
AddBinContent(bin,w);
if (fSumw2.fN) fSumw2.fArray[bin] += w*w;
if (binx == 0 || binx > fXaxis.GetNbins()) return -1;
if (biny == 0 || biny > fYaxis.GetNbins()) return -1;
if (binz == 0 || binz > fZaxis.GetNbins()) return -1;
Double_t x = fXaxis.GetBinCenter(binx);
Double_t y = fYaxis.GetBinCenter(biny);
Double_t z = fZaxis.GetBinCenter(binz);
Double_t v = (w > 0 ? w : -w);
fTsumw += v;
fTsumw2 += v*v;
fTsumwx += v*x;
fTsumwx2 += v*x*x;
fTsumwy += v*y;
fTsumwy2 += v*y*y;
fTsumwxy += v*x*y;
fTsumwz += v*z;
fTsumwz2 += v*z*z;
fTsumwxz += v*x*z;
fTsumwyz += v*y*z;
return bin;
}
Int_t TH3::Fill(const char *namex, Double_t y, const char *namez, Double_t w)
{
Int_t binx, biny, binz, bin;
fEntries++;
binx = fXaxis.FindBin(namex);
biny = fYaxis.FindBin(y);
binz = fZaxis.FindBin(namez);
bin = binx + (fXaxis.GetNbins()+2)*(biny + (fYaxis.GetNbins()+2)*binz);
AddBinContent(bin,w);
if (fSumw2.fN) fSumw2.fArray[bin] += w*w;
if (binx == 0 || binx > fXaxis.GetNbins()) return -1;
if (biny == 0 || biny > fYaxis.GetNbins()) {
if (!fgStatOverflows) return -1;
}
if (binz == 0 || binz > fZaxis.GetNbins()) return -1;
Double_t x = fXaxis.GetBinCenter(binx);
Double_t z = fZaxis.GetBinCenter(binz);
Double_t v = (w > 0 ? w : -w);
fTsumw += v;
fTsumw2 += v*v;
fTsumwx += v*x;
fTsumwx2 += v*x*x;
fTsumwy += v*y;
fTsumwy2 += v*y*y;
fTsumwxy += v*x*y;
fTsumwz += v*z;
fTsumwz2 += v*z*z;
fTsumwxz += v*x*z;
fTsumwyz += v*y*z;
return bin;
}
Int_t TH3::Fill(const char *namex, const char *namey, Double_t z, Double_t w)
{
Int_t binx, biny, binz, bin;
fEntries++;
binx = fXaxis.FindBin(namex);
biny = fYaxis.FindBin(namey);
binz = fZaxis.FindBin(z);
bin = binx + (fXaxis.GetNbins()+2)*(biny + (fYaxis.GetNbins()+2)*binz);
AddBinContent(bin,w);
if (fSumw2.fN) fSumw2.fArray[bin] += w*w;
if (binx == 0 || binx > fXaxis.GetNbins()) return -1;
if (biny == 0 || biny > fYaxis.GetNbins()) return -1;
if (binz == 0 || binz > fZaxis.GetNbins()) {
if (!fgStatOverflows) return -1;
}
Double_t x = fXaxis.GetBinCenter(binx);
Double_t y = fYaxis.GetBinCenter(biny);
Double_t v = (w > 0 ? w : -w);
fTsumw += v;
fTsumw2 += v*v;
fTsumwx += v*x;
fTsumwx2 += v*x*x;
fTsumwy += v*y;
fTsumwy2 += v*y*y;
fTsumwxy += v*x*y;
fTsumwz += v*z;
fTsumwz2 += v*z*z;
fTsumwxz += v*x*z;
fTsumwyz += v*y*z;
return bin;
}
Int_t TH3::Fill(Double_t x, const char *namey, const char *namez, Double_t w)
{
Int_t binx, biny, binz, bin;
fEntries++;
binx = fXaxis.FindBin(x);
biny = fYaxis.FindBin(namey);
binz = fZaxis.FindBin(namez);
bin = binx + (fXaxis.GetNbins()+2)*(biny + (fYaxis.GetNbins()+2)*binz);
AddBinContent(bin,w);
if (fSumw2.fN) fSumw2.fArray[bin] += w*w;
if (binx == 0 || binx > fXaxis.GetNbins()) {
if (!fgStatOverflows) return -1;
}
if (biny == 0 || biny > fYaxis.GetNbins()) return -1;
if (binz == 0 || binz > fZaxis.GetNbins()) return -1;
Double_t y = fYaxis.GetBinCenter(biny);
Double_t z = fZaxis.GetBinCenter(binz);
Double_t v = (w > 0 ? w : -w);
fTsumw += v;
fTsumw2 += v*v;
fTsumwx += v*x;
fTsumwx2 += v*x*x;
fTsumwy += v*y;
fTsumwy2 += v*y*y;
fTsumwxy += v*x*y;
fTsumwz += v*z;
fTsumwz2 += v*z*z;
fTsumwxz += v*x*z;
fTsumwyz += v*y*z;
return bin;
}
Int_t TH3::Fill(Double_t x, const char *namey, Double_t z, Double_t w)
{
Int_t binx, biny, binz, bin;
fEntries++;
binx = fXaxis.FindBin(x);
biny = fYaxis.FindBin(namey);
binz = fZaxis.FindBin(z);
bin = binx + (fXaxis.GetNbins()+2)*(biny + (fYaxis.GetNbins()+2)*binz);
AddBinContent(bin,w);
if (fSumw2.fN) fSumw2.fArray[bin] += w*w;
if (binx == 0 || binx > fXaxis.GetNbins()) {
if (!fgStatOverflows) return -1;
}
if (biny == 0 || biny > fYaxis.GetNbins()) return -1;
if (binz == 0 || binz > fZaxis.GetNbins()) {
if (!fgStatOverflows) return -1;
}
Double_t y = fYaxis.GetBinCenter(biny);
Double_t v = (w > 0 ? w : -w);
fTsumw += v;
fTsumw2 += v*v;
fTsumwx += v*x;
fTsumwx2 += v*x*x;
fTsumwy += v*y;
fTsumwy2 += v*y*y;
fTsumwxy += v*x*y;
fTsumwz += v*z;
fTsumwz2 += v*z*z;
fTsumwxz += v*x*z;
fTsumwyz += v*y*z;
return bin;
}
Int_t TH3::Fill(Double_t x, Double_t y, const char *namez, Double_t w)
{
Int_t binx, biny, binz, bin;
fEntries++;
binx = fXaxis.FindBin(x);
biny = fYaxis.FindBin(y);
binz = fZaxis.FindBin(namez);
bin = binx + (fXaxis.GetNbins()+2)*(biny + (fYaxis.GetNbins()+2)*binz);
AddBinContent(bin,w);
if (fSumw2.fN) fSumw2.fArray[bin] += w*w;
if (binx == 0 || binx > fXaxis.GetNbins()) {
if (!fgStatOverflows) return -1;
}
if (biny == 0 || biny > fYaxis.GetNbins()) {
if (!fgStatOverflows) return -1;
}
if (binz == 0 || binz > fZaxis.GetNbins()) return -1;
Double_t z = fZaxis.GetBinCenter(binz);
Double_t v = (w > 0 ? w : -w);
fTsumw += v;
fTsumw2 += v*v;
fTsumwx += v*x;
fTsumwx2 += v*x*x;
fTsumwy += v*y;
fTsumwy2 += v*y*y;
fTsumwxy += v*x*y;
fTsumwz += v*z;
fTsumwz2 += v*z*z;
fTsumwxz += v*x*z;
fTsumwyz += v*y*z;
return bin;
}
void TH3::FillRandom(const char *fname, Int_t ntimes)
{
Int_t bin, binx, biny, binz, ibin, loop;
Double_t r1, x, y,z, xv[3];
TF1 *f1 = (TF1*)gROOT->GetFunction(fname);
if (!f1) { Error("FillRandom", "Unknown function: %s",fname); return; }
Int_t nbinsx = GetNbinsX();
Int_t nbinsy = GetNbinsY();
Int_t nbinsz = GetNbinsZ();
Int_t nxy = nbinsx*nbinsy;
Int_t nbins = nxy*nbinsz;
Double_t *integral = new Double_t[nbins+1];
ibin = 0;
integral[ibin] = 0;
for (binz=1;binz<=nbinsz;binz++) {
xv[2] = fZaxis.GetBinCenter(binz);
for (biny=1;biny<=nbinsy;biny++) {
xv[1] = fYaxis.GetBinCenter(biny);
for (binx=1;binx<=nbinsx;binx++) {
xv[0] = fXaxis.GetBinCenter(binx);
ibin++;
integral[ibin] = integral[ibin-1] + f1->Eval(xv[0],xv[1],xv[2]);
}
}
}
if (integral[nbins] == 0 ) {
Error("FillRandom", "Integral = zero"); return;
}
for (bin=1;bin<=nbins;bin++) integral[bin] /= integral[nbins];
if (fDimension < 2) nbinsy = -1;
if (fDimension < 3) nbinsz = -1;
for (loop=0;loop<ntimes;loop++) {
r1 = gRandom->Rndm(loop);
ibin = TMath::BinarySearch(nbins,&integral[0],r1);
binz = ibin/nxy;
biny = (ibin - nxy*binz)/nbinsx;
binx = 1 + ibin - nbinsx*(biny + nbinsy*binz);
if (nbinsz) binz++;
if (nbinsy) biny++;
x = fXaxis.GetBinCenter(binx);
y = fYaxis.GetBinCenter(biny);
z = fZaxis.GetBinCenter(binz);
Fill(x,y,z, 1.);
}
delete [] integral;
}
void TH3::FillRandom(TH1 *h, Int_t ntimes)
{
if (!h) { Error("FillRandom", "Null histogram"); return; }
if (fDimension != h->GetDimension()) {
Error("FillRandom", "Histograms with different dimensions"); return;
}
if (h->ComputeIntegral() == 0) return;
TH3 *h3 = (TH3*)h;
Int_t loop;
Double_t x,y,z;
for (loop=0;loop<ntimes;loop++) {
h3->GetRandom3(x,y,z);
Fill(x,y,z,1.);
}
}
void TH3::FitSlicesZ(TF1 *f1, Int_t binminx, Int_t binmaxx, Int_t binminy, Int_t binmaxy, Int_t cut, Option_t *option)
{
Int_t nbinsx = fXaxis.GetNbins();
Int_t nbinsy = fYaxis.GetNbins();
Int_t nbinsz = fZaxis.GetNbins();
if (binminx < 1) binminx = 1;
if (binmaxx > nbinsx) binmaxx = nbinsx;
if (binmaxx < binminx) {binminx = 1; binmaxx = nbinsx;}
if (binminy < 1) binminy = 1;
if (binmaxy > nbinsy) binmaxy = nbinsy;
if (binmaxy < binminy) {binminy = 1; binmaxy = nbinsy;}
if (f1 == 0) {
f1 = (TF1*)gROOT->GetFunction("gaus");
if (f1 == 0) f1 = new TF1("gaus","gaus",fZaxis.GetXmin(),fZaxis.GetXmax());
else f1->SetRange(fZaxis.GetXmin(),fZaxis.GetXmax());
}
const char *fname = f1->GetName();
Int_t npar = f1->GetNpar();
Double_t *parsave = new Double_t[npar];
f1->GetParameters(parsave);
Int_t ipar;
char name[80], title[80];
TH2D *hlist[25];
const TArrayD *xbins = fXaxis.GetXbins();
const TArrayD *ybins = fYaxis.GetXbins();
for (ipar=0;ipar<npar;ipar++) {
sprintf(name,"%s_%d",GetName(),ipar);
sprintf(title,"Fitted value of par[%d]=%s",ipar,f1->GetParName(ipar));
if (xbins->fN == 0) {
hlist[ipar] = new TH2D(name, title,
nbinsx, fXaxis.GetXmin(), fXaxis.GetXmax(),
nbinsy, fYaxis.GetXmin(), fYaxis.GetXmax());
} else {
hlist[ipar] = new TH2D(name, title,
nbinsx, xbins->fArray,
nbinsy, ybins->fArray);
}
hlist[ipar]->GetXaxis()->SetTitle(fXaxis.GetTitle());
hlist[ipar]->GetYaxis()->SetTitle(fYaxis.GetTitle());
}
sprintf(name,"%s_chi2",GetName());
TH2D *hchi2 = new TH2D(name,"chisquare", nbinsx, fXaxis.GetXmin(), fXaxis.GetXmax()
, nbinsy, fYaxis.GetXmin(), fYaxis.GetXmax());
TH1D *hpz = new TH1D("R_temp","_temp",nbinsz, fZaxis.GetXmin(), fZaxis.GetXmax());
Int_t bin,binx,biny,binz;
for (biny=binminy;biny<=binmaxy;biny++) {
Float_t y = fYaxis.GetBinCenter(biny);
for (binx=binminx;binx<=binmaxx;binx++) {
Float_t x = fXaxis.GetBinCenter(binx);
hpz->Reset();
Int_t nfill = 0;
for (binz=1;binz<=nbinsz;binz++) {
bin = GetBin(binx,biny,binz);
Float_t w = GetBinContent(bin);
if (w == 0) continue;
hpz->Fill(fZaxis.GetBinCenter(binz),w);
hpz->SetBinError(binz,GetBinError(bin));
nfill++;
}
if (nfill < cut) continue;
f1->SetParameters(parsave);
hpz->Fit(fname,option);
Int_t npfits = f1->GetNumberFitPoints();
if (npfits > npar && npfits >= cut) {
for (ipar=0;ipar<npar;ipar++) {
hlist[ipar]->Fill(x,y,f1->GetParameter(ipar));
hlist[ipar]->SetCellError(binx,biny,f1->GetParError(ipar));
}
hchi2->Fill(x,y,f1->GetChisquare()/(npfits-npar));
}
}
}
delete [] parsave;
delete hpz;
}
Double_t TH3::GetBinWithContent3(Double_t c, Int_t &binx, Int_t &biny, Int_t &binz, Int_t firstx, Int_t lastx, Int_t firsty, Int_t lasty, Int_t firstz, Int_t lastz, Double_t maxdiff) const
{
if (fDimension != 3) {
binx = 0;
biny = 0;
binz = 0;
Error("GetBinWithContent3","function is only valid for 3-D histograms");
return 0;
}
if (firstx <= 0) firstx = 1;
if (lastx < firstx) lastx = fXaxis.GetNbins();
if (firsty <= 0) firsty = 1;
if (lasty < firsty) lasty = fYaxis.GetNbins();
if (firstz <= 0) firstz = 1;
if (lastz < firstz) lastz = fZaxis.GetNbins();
Int_t binminx = 0, binminy=0, binminz=0;
Double_t diff, curmax = 1.e240;
for (Int_t k=firstz;k<=lastz;k++) {
for (Int_t j=firsty;j<=lasty;j++) {
for (Int_t i=firstx;i<=lastx;i++) {
diff = TMath::Abs(GetBinContent(i,j,k)-c);
if (diff <= 0) {binx = i; biny=j; binz=k; return diff;}
if (diff < curmax && diff <= maxdiff) {curmax = diff, binminx=i; binminy=j;binminz=k;}
}
}
}
binx = binminx;
biny = binminy;
binz = binminz;
return curmax;
}
Double_t TH3::GetCorrelationFactor(Int_t axis1, Int_t axis2) const
{
if (axis1 < 1 || axis2 < 1 || axis1 > 3 || axis2 > 3) {
Error("GetCorrelationFactor","Wrong parameters");
return 0;
}
if (axis1 == axis2) return 1;
Double_t rms1 = GetRMS(axis1);
if (rms1 == 0) return 0;
Double_t rms2 = GetRMS(axis2);
if (rms2 == 0) return 0;
return GetCovariance(axis1,axis2)/rms1/rms2;
}
Double_t TH3::GetCovariance(Int_t axis1, Int_t axis2) const
{
if (axis1 < 1 || axis2 < 1 || axis1 > 3 || axis2 > 3) {
Error("GetCovariance","Wrong parameters");
return 0;
}
Double_t stats[11];
GetStats(stats);
Double_t sumw = stats[0];
Double_t sumw2 = stats[1];
Double_t sumwx = stats[2];
Double_t sumwx2 = stats[3];
Double_t sumwy = stats[4];
Double_t sumwy2 = stats[5];
Double_t sumwxy = stats[6];
Double_t sumwz = stats[7];
Double_t sumwz2 = stats[8];
Double_t sumwxz = stats[9];
Double_t sumwyz = stats[10];
if (sumw == 0) return 0;
if (axis1 == 1 && axis2 == 1) {
return TMath::Abs(sumwx2/sumw - sumwx*sumwx/sumw2);
}
if (axis1 == 2 && axis2 == 2) {
return TMath::Abs(sumwy2/sumw - sumwy*sumwy/sumw2);
}
if (axis1 == 3 && axis2 == 3) {
return TMath::Abs(sumwz2/sumw - sumwz*sumwz/sumw2);
}
if ((axis1 == 1 && axis2 == 2) || (axis1 == 2 && axis1 == 1)) {
return sumwxy/sumw - sumwx/sumw*sumwy/sumw;
}
if ((axis1 == 1 && axis2 == 3) || (axis1 == 3 && axis2 == 1)) {
return sumwxz/sumw - sumwx/sumw*sumwz/sumw;
}
if ((axis1 == 2 && axis2 == 3) || (axis1 == 3 && axis2 == 2)) {
return sumwyz/sumw - sumwy/sumw*sumwz/sumw;
}
return 0;
}
void TH3::GetRandom3(Double_t &x, Double_t &y, Double_t &z)
{
Int_t nbinsx = GetNbinsX();
Int_t nbinsy = GetNbinsY();
Int_t nbinsz = GetNbinsZ();
Int_t nxy = nbinsx*nbinsy;
Int_t nbins = nxy*nbinsz;
Double_t integral;
if (fIntegral) {
if (fIntegral[nbins+1] != fEntries) integral = ComputeIntegral();
} else {
integral = ComputeIntegral();
if (integral == 0 || fIntegral == 0) return;
}
Float_t r1 = gRandom->Rndm();
Int_t ibin = TMath::BinarySearch(nbins,fIntegral,r1);
Int_t binz = ibin/nxy;
Int_t biny = (ibin - nxy*binz)/nbinsx;
Int_t binx = ibin - nbinsx*(biny + nbinsy*binz);
x = fXaxis.GetBinLowEdge(binx+1);
if (r1 > fIntegral[ibin]) x +=
fXaxis.GetBinWidth(binx+1)*(r1-fIntegral[ibin])/(fIntegral[ibin+1] - fIntegral[ibin]);
y = fYaxis.GetBinLowEdge(biny+1) + fYaxis.GetBinWidth(biny+1)*gRandom->Rndm();
z = fZaxis.GetBinLowEdge(binz+1) + fZaxis.GetBinWidth(binz+1)*gRandom->Rndm();
}
void TH3::GetStats(Double_t *stats) const
{
if (fBuffer) ((TH3*)this)->BufferEmpty();
Int_t bin, binx, biny, binz;
Double_t w,err;
Double_t x,y,z;
if (fTsumw == 0 || fXaxis.TestBit(TAxis::kAxisRange) || fYaxis.TestBit(TAxis::kAxisRange) || fZaxis.TestBit(TAxis::kAxisRange)) {
for (bin=0;bin<9;bin++) stats[bin] = 0;
Int_t firstBinX = fXaxis.GetFirst();
Int_t lastBinX = fXaxis.GetLast();
Int_t firstBinY = fYaxis.GetFirst();
Int_t lastBinY = fYaxis.GetLast();
Int_t firstBinZ = fZaxis.GetFirst();
Int_t lastBinZ = fZaxis.GetLast();
if (fgStatOverflows) {
if ( !fXaxis.TestBit(TAxis::kAxisRange) ) {
if (firstBinX == 1) firstBinX = 0;
if (lastBinX == fXaxis.GetNbins() ) lastBinX += 1;
}
if ( !fYaxis.TestBit(TAxis::kAxisRange) ) {
if (firstBinY == 1) firstBinY = 0;
if (lastBinY == fYaxis.GetNbins() ) lastBinY += 1;
}
if ( !fZaxis.TestBit(TAxis::kAxisRange) ) {
if (firstBinZ == 1) firstBinZ = 0;
if (lastBinZ == fZaxis.GetNbins() ) lastBinZ += 1;
}
}
for (binz = firstBinZ; binz <= lastBinZ; binz++) {
z = fZaxis.GetBinCenter(binz);
for (biny = firstBinY; biny <= lastBinY; biny++) {
y = fYaxis.GetBinCenter(biny);
for (binx = firstBinX; binx <= lastBinX; binx++) {
bin = GetBin(binx,biny,binz);
x = fXaxis.GetBinCenter(binx);
w = TMath::Abs(GetBinContent(bin));
err = TMath::Abs(GetBinError(bin));
stats[0] += w;
stats[1] += err*err;
stats[2] += w*x;
stats[3] += w*x*x;
stats[4] += w*y;
stats[5] += w*y*y;
stats[6] += w*x*y;
stats[7] += w*z;
stats[8] += w*z*z;
stats[9] += w*x*z;
stats[10]+= w*y*z;
}
}
}
} else {
stats[0] = fTsumw;
stats[1] = fTsumw2;
stats[2] = fTsumwx;
stats[3] = fTsumwx2;
stats[4] = fTsumwy;
stats[5] = fTsumwy2;
stats[6] = fTsumwxy;
stats[7] = fTsumwz;
stats[8] = fTsumwz2;
stats[9] = fTsumwxz;
stats[10]= fTsumwyz;
}
}
Double_t TH3::Integral(Option_t *option) const
{
return Integral(fXaxis.GetFirst(),fXaxis.GetLast(),
fYaxis.GetFirst(),fYaxis.GetLast(),
fZaxis.GetFirst(),fZaxis.GetLast(),option);
}
Double_t TH3::Integral(Int_t binx1, Int_t binx2, Int_t biny1, Int_t biny2, Int_t binz1, Int_t binz2, Option_t *option) const
{
Int_t nbinsx = GetNbinsX();
Int_t nbinsy = GetNbinsY();
Int_t nbinsz = GetNbinsZ();
if (binx1 < 0) binx1 = 0;
if (binx2 > nbinsx+1) binx2 = nbinsx+1;
if (binx2 < binx1) binx2 = nbinsx;
if (biny1 < 0) biny1 = 0;
if (biny2 > nbinsy+1) biny2 = nbinsy+1;
if (biny2 < biny1) biny2 = nbinsy;
if (binz1 < 0) binz1 = 0;
if (binz2 > nbinsz+1) binz2 = nbinsz+1;
if (binz2 < binz1) binz2 = nbinsz;
Double_t integral = 0;
TString opt = option;
opt.ToLower();
Bool_t width = kFALSE;
if (opt.Contains("width")) width = kTRUE;
Int_t bin, binx, biny, binz;
for (binz=binz1;binz<=binz2;binz++) {
for (biny=biny1;biny<=biny2;biny++) {
for (binx=binx1;binx<=binx2;binx++) {
bin = binx +(nbinsx+2)*(biny + (nbinsy+2)*binz);
if (width) integral += GetBinContent(bin)*fXaxis.GetBinWidth(binx)*fYaxis.GetBinWidth(biny)*fZaxis.GetBinWidth(binz);
else integral += GetBinContent(bin);
}
}
}
return integral;
}
Double_t TH3::KolmogorovTest(const TH1 *h2, Option_t *option) const
{
TString opt = option;
opt.ToUpper();
Double_t prb = 0;
TH1 *h1 = (TH1*)this;
if (h2 == 0) return 0;
TAxis *xaxis1 = h1->GetXaxis();
TAxis *xaxis2 = h2->GetXaxis();
TAxis *yaxis1 = h1->GetYaxis();
TAxis *yaxis2 = h2->GetYaxis();
TAxis *zaxis1 = h1->GetZaxis();
TAxis *zaxis2 = h2->GetZaxis();
Int_t ncx1 = xaxis1->GetNbins();
Int_t ncx2 = xaxis2->GetNbins();
Int_t ncy1 = yaxis1->GetNbins();
Int_t ncy2 = yaxis2->GetNbins();
Int_t ncz1 = zaxis1->GetNbins();
Int_t ncz2 = zaxis2->GetNbins();
if (h1->GetDimension() != 3 || h2->GetDimension() != 3) {
Error("KolmogorovTest","Histograms must be 3-D\n");
return 0;
}
if (ncx1 != ncx2) {
Error("KolmogorovTest","Number of channels in X is different, %d and %d\n",ncx1,ncx2);
return 0;
}
if (ncy1 != ncy2) {
Error("KolmogorovTest","Number of channels in Y is different, %d and %d\n",ncy1,ncy2);
return 0;
}
if (ncz1 != ncz2) {
Error("KolmogorovTest","Number of channels in Z is different, %d and %d\n",ncz1,ncz2);
return 0;
}
Bool_t afunc1 = kFALSE;
Bool_t afunc2 = kFALSE;
Double_t difprec = 1e-5;
Double_t diff1 = TMath::Abs(xaxis1->GetXmin() - xaxis2->GetXmin());
Double_t diff2 = TMath::Abs(xaxis1->GetXmax() - xaxis2->GetXmax());
if (diff1 > difprec || diff2 > difprec) {
Error("KolmogorovTest","histograms with different binning along X");
return 0;
}
diff1 = TMath::Abs(yaxis1->GetXmin() - yaxis2->GetXmin());
diff2 = TMath::Abs(yaxis1->GetXmax() - yaxis2->GetXmax());
if (diff1 > difprec || diff2 > difprec) {
Error("KolmogorovTest","histograms with different binning along Y");
return 0;
}
diff1 = TMath::Abs(zaxis1->GetXmin() - zaxis2->GetXmin());
diff2 = TMath::Abs(zaxis1->GetXmax() - zaxis2->GetXmax());
if (diff1 > difprec || diff2 > difprec) {
Error("KolmogorovTest","histograms with different binning along Z");
return 0;
}
Int_t ibeg = 1, jbeg = 1, kbeg = 1;
Int_t iend = ncx1, jend = ncy1, kend = ncz1;
if (opt.Contains("U")) {ibeg = 0; jbeg = 0; kbeg = 0;}
if (opt.Contains("O")) {iend = ncx1+1; jend = ncy1+1; kend = ncz1+1;}
Int_t i,j,k,bin;
Double_t hsav;
Double_t sum1 = 0;
Double_t tsum1 = 0;
for (i=0;i<=ncx1+1;i++) {
for (j=0;j<=ncy1+1;j++) {
for (k=0;k<=ncz1+1;k++) {
bin = h1->GetBin(i,j,k);
hsav = h1->GetBinContent(bin);
tsum1 += hsav;
if (i >= ibeg && i <= iend && j >= jbeg && j <= jend && k >= kbeg && k <= kend) sum1 += hsav;
}
}
}
Double_t sum2 = 0;
Double_t tsum2 = 0;
for (i=0;i<=ncx1+1;i++) {
for (j=0;j<=ncy1+1;j++) {
for (k=0;k<=ncz1+1;k++) {
bin = h1->GetBin(i,j,k);
hsav = h1->GetBinContent(bin);
tsum2 += hsav;
if (i >= ibeg && i <= iend && j >= jbeg && j <= jend&& k >= kbeg && k <= kend) sum2 += hsav;
}
}
}
if (sum1 == 0) {
Error("KolmogorovTest","Integral is zero for h1=%s\n",h1->GetName());
return 0;
}
if (sum2 == 0) {
Error("KolmogorovTest","Integral is zero for h2=%s\n",h2->GetName());
return 0;
}
Double_t num1 = h1->GetEntries();
Double_t num2 = h2->GetEntries();
if (num1 != tsum1) {
Warning("KolmogorovTest","Saturation or weighted events for h1=%s, num1=%g, tsum1=%g\n",h1->GetName(),num1,tsum1);
}
if (num2 != tsum2) {
Warning("KolmogorovTest","Saturation or weighted events for h2=%s, num2=%g, tsum2=%g\n",h2->GetName(),num2,tsum2);
}
Double_t s1 = 1/sum1;
Double_t s2 = 1/sum2;
Double_t dfmax = 0;
Double_t rsum1=0, rsum2=0;
for (i=ibeg;i<=iend;i++) {
for (j=jbeg;j<=jend;j++) {
for (k=kbeg;k<=kend;k++) {
bin = h1->GetBin(i,j,k);
rsum1 += s1*h1->GetBinContent(bin);
rsum2 += s2*h2->GetBinContent(bin);
dfmax = TMath::Max(dfmax, TMath::Abs(rsum1-rsum2));
}
}
}
Double_t dfmax2 = 0;
rsum1=0, rsum2=0;
for (k=kbeg;k<=kend;k++) {
for (j=jbeg;j<=jend;j++) {
for (i=ibeg;i<=iend;i++) {
bin = h1->GetBin(i,j,k);
rsum1 += s1*h1->GetBinContent(bin);
rsum2 += s2*h2->GetBinContent(bin);
dfmax2 = TMath::Max(dfmax2, TMath::Abs(rsum1-rsum2));
}
}
}
Double_t factnm;
if (afunc1) factnm = TMath::Sqrt(sum2);
else if (afunc2) factnm = TMath::Sqrt(sum1);
else factnm = TMath::Sqrt(sum1*sum2/(sum1+sum2));
Double_t z = dfmax*factnm;
Double_t z2 = dfmax2*factnm;
prb = TMath::KolmogorovProb(0.5*(z+z2));
Double_t prb1=0, prb2=0;
Double_t resum1, resum2, chi2, d12;
if (opt.Contains("N")) {
prb1 = prb;
resum1 = sum1; if (afunc1) resum1 = 0;
resum2 = sum2; if (afunc2) resum2 = 0;
d12 = sum1-sum2;
chi2 = d12*d12/(resum1+resum2);
prb2 = TMath::Prob(chi2,1);
if (prb > 0 && prb2 > 0) prb = prb*prb2*(1-TMath::Log(prb*prb2));
else prb = 0;
}
if (opt.Contains("D")) {
printf(" Kolmo Prob h1 = %s, sum1=%g\n",h1->GetName(),sum1);
printf(" Kolmo Prob h2 = %s, sum2=%g\n",h2->GetName(),sum2);
printf(" Kolmo Probabil = %f, Max Dist = %g\n",prb,dfmax);
if (opt.Contains("N"))
printf(" Kolmo Probabil = %f for shape alone, =%f for normalisation alone\n",prb1,prb2);
}
if (TMath::Abs(rsum1-1) > 0.002) Warning("KolmogorovTest","Numerical problems with h1=%s\n",h1->GetName());
if (TMath::Abs(rsum2-1) > 0.002) Warning("KolmogorovTest","Numerical problems with h2=%s\n",h2->GetName());
return prb;
}
Long64_t TH3::Merge(TCollection *list)
{
if (!list) return 0;
if (list->IsEmpty()) return (Int_t) GetEntries();
TList inlist;
TH1* hclone = (TH1*)Clone("FirstClone");
R__ASSERT(hclone);
BufferEmpty(1);
Reset();
SetEntries(0);
inlist.Add(hclone);
inlist.AddAll(list);
TAxis newXAxis;
TAxis newYAxis;
TAxis newZAxis;
Bool_t initialLimitsFound = kFALSE;
Bool_t same = kTRUE;
Bool_t allHaveLimits = kTRUE;
TIter next(&inlist);
while (TObject *o = next()) {
TH3* h = dynamic_cast<TH3*> (o);
if (!h) {
Error("Add","Attempt to add object of class: %s to a %s",
o->ClassName(),this->ClassName());
return -1;
}
Bool_t hasLimits = h->GetXaxis()->GetXmin() < h->GetXaxis()->GetXmax();
allHaveLimits = allHaveLimits && hasLimits;
if (hasLimits) {
h->BufferEmpty();
if (!initialLimitsFound) {
initialLimitsFound = kTRUE;
newXAxis.Set(h->GetXaxis()->GetNbins(), h->GetXaxis()->GetXmin(),
h->GetXaxis()->GetXmax());
newYAxis.Set(h->GetYaxis()->GetNbins(), h->GetYaxis()->GetXmin(),
h->GetYaxis()->GetXmax());
newZAxis.Set(h->GetZaxis()->GetNbins(), h->GetZaxis()->GetXmin(),
h->GetZaxis()->GetXmax());
}
else {
if (!RecomputeAxisLimits(newXAxis, *(h->GetXaxis()))) {
Error("Merge", "Cannot merge histograms - limits are inconsistent:\n "
"first: (%d, %f, %f), second: (%d, %f, %f)",
newXAxis.GetNbins(), newXAxis.GetXmin(), newXAxis.GetXmax(),
h->GetXaxis()->GetNbins(), h->GetXaxis()->GetXmin(),
h->GetXaxis()->GetXmax());
}
if (!RecomputeAxisLimits(newYAxis, *(h->GetYaxis()))) {
Error("Merge", "Cannot merge histograms - limits are inconsistent:\n "
"first: (%d, %f, %f), second: (%d, %f, %f)",
newYAxis.GetNbins(), newYAxis.GetXmin(), newYAxis.GetXmax(),
h->GetYaxis()->GetNbins(), h->GetYaxis()->GetXmin(),
h->GetYaxis()->GetXmax());
}
if (!RecomputeAxisLimits(newZAxis, *(h->GetZaxis()))) {
Error("Merge", "Cannot merge histograms - limits are inconsistent:\n "
"first: (%d, %f, %f), second: (%d, %f, %f)",
newZAxis.GetNbins(), newZAxis.GetXmin(), newZAxis.GetXmax(),
h->GetZaxis()->GetNbins(), h->GetZaxis()->GetXmin(),
h->GetZaxis()->GetXmax());
}
}
}
}
next.Reset();
same = same && SameLimitsAndNBins(newXAxis, *GetXaxis())
&& SameLimitsAndNBins(newYAxis, *GetYaxis());
if (!same && initialLimitsFound)
SetBins(newXAxis.GetNbins(), newXAxis.GetXmin(), newXAxis.GetXmax(),
newYAxis.GetNbins(), newYAxis.GetXmin(), newYAxis.GetXmax(),
newZAxis.GetNbins(), newZAxis.GetXmin(), newZAxis.GetXmax());
if (!allHaveLimits) {
while (TH3* h = dynamic_cast<TH3*> (next())) {
if (h->GetXaxis()->GetXmin() >= h->GetXaxis()->GetXmax() && h->fBuffer) {
Int_t nbentries = (Int_t)h->fBuffer[0];
for (Int_t i = 0; i < nbentries; i++)
Fill(h->fBuffer[4*i + 2], h->fBuffer[4*i + 3],
h->fBuffer[4*i + 4], h->fBuffer[4*i + 1]);
}
}
if (!initialLimitsFound)
return (Int_t) GetEntries();
next.Reset();
}
const Int_t kNstat = 11;
Double_t stats[kNstat], totstats[kNstat];
for (Int_t i=0;i<kNstat;i++) {totstats[i] = stats[i] = 0;}
GetStats(totstats);
Double_t nentries = GetEntries();
Int_t binx, biny, binz, ix, iy, iz, nx, ny, nz, bin, ibin;
Double_t cu;
Int_t nbix = fXaxis.GetNbins();
Int_t nbiy = fYaxis.GetNbins();
Bool_t canRebin=TestBit(kCanRebin);
ResetBit(kCanRebin);
while (TH1* h=(TH1*)next()) {
if (h->GetXaxis()->GetXmin() < h->GetXaxis()->GetXmax()) {
h->GetStats(stats);
for (Int_t i = 0; i < kNstat; i++)
totstats[i] += stats[i];
nentries += h->GetEntries();
nx = h->GetXaxis()->GetNbins();
ny = h->GetYaxis()->GetNbins();
nz = h->GetZaxis()->GetNbins();
for (binz = 0; binz <= nz + 1; binz++) {
iz = fZaxis.FindBin(h->GetZaxis()->GetBinCenter(binz));
for (biny = 0; biny <= ny + 1; biny++) {
iy = fYaxis.FindBin(h->GetYaxis()->GetBinCenter(biny));
for (binx = 0; binx <= nx + 1; binx++) {
ix = fXaxis.FindBin(h->GetXaxis()->GetBinCenter(binx));
bin = binx +(nx+2)*(biny + (ny+2)*binz);
ibin = ix +(nbix+2)*(iy + (nbiy+2)*iz);
cu = h->GetBinContent(bin);
if ((!same) && (binx == 0 || binx == nx + 1
|| biny == 0 || biny == ny + 1
|| binz == 0 || binz == nz + 1)) {
if (cu != 0) {
Error("Merge", "Cannot merge histograms - the histograms have"
" different limits and undeflows/overflows are present."
" The initial histogram is now broken!");
return -1;
}
}
AddBinContent(ibin,cu);
if (fSumw2.fN) {
Double_t error1 = h->GetBinError(bin);
fSumw2.fArray[ibin] += error1*error1;
}
}
}
}
}
}
if (canRebin) SetBit(kCanRebin);
PutStats(totstats);
SetEntries(nentries);
inlist.Remove(hclone);
delete hclone;
return (Long64_t)nentries;
}
TH1D *TH3::ProjectionZ(const char *name, Int_t ixmin, Int_t ixmax, Int_t iymin, Int_t iymax, Option_t *option) const
{
TString opt = option;
opt.ToLower();
Int_t nx = GetNbinsX();
Int_t ny = GetNbinsY();
Int_t nz = GetNbinsZ();
if (ixmin < 0) ixmin = 1;
if (ixmax < 0) ixmax = nx;
if (ixmax > nx+1) ixmax = nx;
if (iymin < 0) iymin = 1;
if (iymax < 0) iymax = ny;
if (iymax > ny+1) iymax = ny;
char *pname = (char*)name;
if (strcmp(name,"_pz") == 0) {
Int_t nch = strlen(GetName()) + 4;
pname = new char[nch];
sprintf(pname,"%s%s",GetName(),name);
}
TH1D *h1=0;
TObject *h1obj = gROOT->FindObject(pname);
if (h1obj && h1obj->InheritsFrom("TH1D")) {
h1 = (TH1D*)h1obj;
h1->Reset();
}
if (!h1) {
const TArrayD *bins = fZaxis.GetXbins();
if (bins->fN == 0) {
h1 = new TH1D(pname,GetTitle(),nz,fZaxis.GetXmin(),fZaxis.GetXmax());
} else {
h1 = new TH1D(pname,GetTitle(),nz,bins->fArray);
}
}
if (opt.Contains("e")) h1->Sumw2();
if (pname != name) delete [] pname;
h1->GetXaxis()->ImportAttributes(this->GetZaxis());
THashList* labels=GetZaxis()->GetLabels();
if (labels) {
TIter iL(labels);
TObjString* lb;
Int_t i = 1;
while ((lb=(TObjString*)iL())) {
h1->GetXaxis()->SetBinLabel(i,lb->String().Data());
i++;
}
}
h1->SetLineColor(this->GetLineColor());
h1->SetFillColor(this->GetFillColor());
h1->SetMarkerColor(this->GetMarkerColor());
h1->SetMarkerStyle(this->GetMarkerStyle());
Float_t cont,e,e1;
Double_t entries = 0;
Double_t newerror = 0;
for (Int_t ixbin=ixmin;ixbin<=ixmax;ixbin++){
for (Int_t iybin=iymin;iybin<=iymax;iybin++){
for (Int_t binz=0;binz<=(nz+1);binz++){
Int_t bin = GetBin(ixbin,iybin,binz);
cont = GetBinContent(bin);
if (h1->GetSumw2N()) {
e = GetBinError(bin);
e1 = h1->GetBinError(binz);
newerror = TMath::Sqrt(e*e + e1*e1);
}
if (cont) {
h1->Fill(fZaxis.GetBinCenter(binz), cont);
entries += cont;
}
if (h1->GetSumw2N()) h1->SetBinError(binz,newerror);
}
}
}
if (iymin <=1 && iymax >= ny && ixmin <=1 && ixmax >= nx) h1->SetEntries(fEntries);
else h1->SetEntries(entries);
if (opt.Contains("d")) {
TVirtualPad *padsav = gPad;
TVirtualPad *pad = gROOT->GetSelectedPad();
if (pad) pad->cd();
char optin[100];
strcpy(optin,opt.Data());
char *d = (char*)strstr(optin,"d"); if (d) {*d = ' '; if (*(d+1) == 0) *d=0;}
char *e = (char*)strstr(optin,"e"); if (e) {*e = ' '; if (*(e+1) == 0) *e=0;}
if (!gPad->FindObject(h1)) {
h1->Draw(optin);
} else {
h1->Paint(optin);
}
if (padsav) padsav->cd();
}
return h1;
}
TH1 *TH3::Project3D(Option_t *option) const
{
TString opt = option; opt.ToLower();
Int_t ixmin = fXaxis.GetFirst();
Int_t ixmax = fXaxis.GetLast();
Int_t iymin = fYaxis.GetFirst();
Int_t iymax = fYaxis.GetLast();
Int_t izmin = fZaxis.GetFirst();
Int_t izmax = fZaxis.GetLast();
Int_t nx = ixmax-ixmin+1;
Int_t ny = iymax-iymin+1;
Int_t nz = izmax-izmin+1;
Int_t pcase = 0;
if (opt.Contains("x")) pcase = 1;
if (opt.Contains("y")) pcase = 2;
if (opt.Contains("z")) pcase = 3;
if (opt.Contains("xy")) pcase = 4;
if (opt.Contains("yx")) pcase = 5;
if (opt.Contains("xz")) pcase = 6;
if (opt.Contains("zx")) pcase = 7;
if (opt.Contains("yz")) pcase = 8;
if (opt.Contains("zy")) pcase = 9;
TH1D *h1 = 0;
TH2D *h2 = 0;
Int_t nch = strlen(GetName()) +opt.Length() +2;
char *name = new char[nch];
sprintf(name,"%s_%s",GetName(),option);
nch = strlen(GetTitle()) +opt.Length() +2;
char *title = new char[nch];
sprintf(title,"%s_%s",GetTitle(),option);
TObject *h1obj = gROOT->FindObject(name);
if (h1obj && !h1obj->InheritsFrom("TH1")) h1obj = 0;
const TArrayD *bins;
const TArrayD *xbins;
const TArrayD *ybins;
const TArrayD *zbins;
switch (pcase) {
case 1:
if (h1obj && h1obj->InheritsFrom("TH1D")) {
h1 = (TH1D*)h1obj;
h1->Reset();
break;
}
bins = fXaxis.GetXbins();
if (bins->fN == 0) {
h1 = new TH1D(name,title,nx,fXaxis.GetBinLowEdge(ixmin),fXaxis.GetBinUpEdge(ixmax));
} else {
h1 = new TH1D(name,title,nx,&bins->fArray[ixmin-1]);
}
break;
case 2:
if (h1obj && h1obj->InheritsFrom("TH1D")) {
h1 = (TH1D*)h1obj;
h1->Reset();
break;
}
bins = fYaxis.GetXbins();
if (bins->fN == 0) {
h1 = new TH1D(name,title,ny,fYaxis.GetBinLowEdge(iymin),fYaxis.GetBinUpEdge(iymax));
} else {
h1 = new TH1D(name,title,ny,&bins->fArray[iymin-1]);
}
break;
case 3:
if (h1obj && h1obj->InheritsFrom("TH1D")) {
h1 = (TH1D*)h1obj;
h1->Reset();
break;
}
bins = fZaxis.GetXbins();
if (bins->fN == 0) {
h1 = new TH1D(name,title,nz,fZaxis.GetBinLowEdge(izmin),fZaxis.GetBinUpEdge(izmax));
} else {
h1 = new TH1D(name,title,nz,&bins->fArray[izmin-1]);
}
break;
case 4:
if (h1obj && h1obj->InheritsFrom("TH2D")) {
h2 = (TH2D*)h1obj;
h2->Reset();
break;
}
xbins = fXaxis.GetXbins();
ybins = fYaxis.GetXbins();
if (xbins->fN == 0 && ybins->fN == 0) {
h2 = new TH2D(name,title,ny,fYaxis.GetBinLowEdge(iymin),fYaxis.GetBinUpEdge(iymax)
,nx,fXaxis.GetBinLowEdge(ixmin),fXaxis.GetBinUpEdge(ixmax));
} else if (ybins->fN == 0) {
h2 = new TH2D(name,title,ny,fYaxis.GetBinLowEdge(iymin),fYaxis.GetBinUpEdge(iymax)
,nx,&xbins->fArray[ixmin-1]);
} else if (xbins->fN == 0) {
h2 = new TH2D(name,title,ny,&ybins->fArray[iymin-1]
,nx,fXaxis.GetBinLowEdge(ixmin),fXaxis.GetBinUpEdge(ixmax));
} else {
h2 = new TH2D(name,title,ny,&ybins->fArray[iymin-1],nx,&xbins->fArray[ixmin-1]);
}
break;
case 5:
if (h1obj && h1obj->InheritsFrom("TH2D")) {
h2 = (TH2D*)h1obj;
h2->Reset();
break;
}
xbins = fXaxis.GetXbins();
ybins = fYaxis.GetXbins();
if (xbins->fN == 0 && ybins->fN == 0) {
h2 = new TH2D(name,title,nx,fXaxis.GetBinLowEdge(ixmin),fXaxis.GetBinUpEdge(ixmax)
,ny,fYaxis.GetBinLowEdge(iymin),fYaxis.GetBinUpEdge(iymax));
} else if (xbins->fN == 0) {
h2 = new TH2D(name,title,nx,fXaxis.GetBinLowEdge(ixmin),fXaxis.GetBinUpEdge(ixmax)
,ny,&ybins->fArray[iymin-1]);
} else if (ybins->fN == 0) {
h2 = new TH2D(name,title,nx,&xbins->fArray[ixmin-1]
,ny,fYaxis.GetBinLowEdge(iymin),fYaxis.GetBinUpEdge(iymax));
} else {
h2 = new TH2D(name,title,nx,&xbins->fArray[ixmin-1],ny,&ybins->fArray[iymin-1]);
}
break;
case 6:
if (h1obj && h1obj->InheritsFrom("TH2D")) {
h2 = (TH2D*)h1obj;
h2->Reset();
break;
}
xbins = fXaxis.GetXbins();
zbins = fZaxis.GetXbins();
if (xbins->fN == 0 && zbins->fN == 0) {
h2 = new TH2D(name,title,nz,fZaxis.GetBinLowEdge(izmin),fZaxis.GetBinUpEdge(izmax)
,nx,fXaxis.GetBinLowEdge(ixmin),fXaxis.GetBinUpEdge(ixmax));
} else if (zbins->fN == 0) {
h2 = new TH2D(name,title,nz,fZaxis.GetBinLowEdge(izmin),fZaxis.GetBinUpEdge(izmax)
,nx,&xbins->fArray[ixmin-1]);
} else if (xbins->fN == 0) {
h2 = new TH2D(name,title,nz,&zbins->fArray[izmin-1]
,nx,fXaxis.GetBinLowEdge(ixmin),fXaxis.GetBinUpEdge(ixmax));
} else {
h2 = new TH2D(name,title,nz,&zbins->fArray[izmin-1],nx,&xbins->fArray[ixmin-1]);
}
break;
case 7:
if (h1obj && h1obj->InheritsFrom("TH2D")) {
h2 = (TH2D*)h1obj;
h2->Reset();
break;
}
xbins = fXaxis.GetXbins();
zbins = fZaxis.GetXbins();
if (xbins->fN == 0 && zbins->fN == 0) {
h2 = new TH2D(name,title,nx,fXaxis.GetBinLowEdge(ixmin),fXaxis.GetBinUpEdge(ixmax)
,nz,fZaxis.GetBinLowEdge(izmin),fZaxis.GetBinUpEdge(izmax));
} else if (xbins->fN == 0) {
h2 = new TH2D(name,title,nx,fXaxis.GetBinLowEdge(ixmin),fXaxis.GetBinUpEdge(ixmax)
,nz,&zbins->fArray[izmin-1]);
} else if (zbins->fN == 0) {
h2 = new TH2D(name,title,nx,&xbins->fArray[ixmin-1]
,nz,fZaxis.GetBinLowEdge(izmin),fZaxis.GetBinUpEdge(izmax));
} else {
h2 = new TH2D(name,title,nx,&xbins->fArray[ixmin-1],nz,&zbins->fArray[izmin-1]);
}
break;
case 8:
if (h1obj && h1obj->InheritsFrom("TH2D")) {
h2 = (TH2D*)h1obj;
h2->Reset();
break;
}
ybins = fYaxis.GetXbins();
zbins = fZaxis.GetXbins();
if (ybins->fN == 0 && zbins->fN == 0) {
h2 = new TH2D(name,title,nz,fZaxis.GetBinLowEdge(izmin),fZaxis.GetBinUpEdge(izmax)
,ny,fYaxis.GetBinLowEdge(iymin),fYaxis.GetBinUpEdge(iymax));
} else if (zbins->fN == 0) {
h2 = new TH2D(name,title,nz,fZaxis.GetBinLowEdge(izmin),fZaxis.GetBinUpEdge(izmax)
,ny,&ybins->fArray[iymin-1]);
} else if (ybins->fN == 0) {
h2 = new TH2D(name,title,nz,&zbins->fArray[izmin-1]
,ny,fYaxis.GetBinLowEdge(iymin),fYaxis.GetBinUpEdge(iymax));
} else {
h2 = new TH2D(name,title,nz,&zbins->fArray[izmin-1],ny,&ybins->fArray[iymin-1]);
}
break;
case 9:
if (h1obj && h1obj->InheritsFrom("TH2D")) {
h2 = (TH2D*)h1obj;
h2->Reset();
break;
}
ybins = fYaxis.GetXbins();
zbins = fZaxis.GetXbins();
if (ybins->fN == 0 && zbins->fN == 0) {
h2 = new TH2D(name,title,ny,fYaxis.GetBinLowEdge(iymin),fYaxis.GetBinUpEdge(iymax)
,nz,fZaxis.GetBinLowEdge(izmin),fZaxis.GetBinUpEdge(izmax));
} else if (ybins->fN == 0) {
h2 = new TH2D(name,title,ny,fYaxis.GetBinLowEdge(iymin),fYaxis.GetBinUpEdge(iymax)
,nz,&zbins->fArray[izmin-1]);
} else if (zbins->fN == 0) {
h2 = new TH2D(name,title,ny,&ybins->fArray[iymin-1]
,nz,fZaxis.GetBinLowEdge(izmin),fZaxis.GetBinUpEdge(izmax));
} else {
h2 = new TH2D(name,title,ny,&ybins->fArray[iymin-1],nz,&zbins->fArray[izmin-1]);
}
break;
}
if (pcase >=1 && pcase <=3) {
THashList* labels = 0;
switch (pcase) {
case 1:
h1->GetXaxis()->ImportAttributes(this->GetXaxis());
labels = GetXaxis()->GetLabels();
break;
case 2:
h1->GetXaxis()->ImportAttributes(this->GetYaxis());
labels = GetYaxis()->GetLabels();
break;
case 3:
h1->GetXaxis()->ImportAttributes(this->GetZaxis());
labels = GetZaxis()->GetLabels();
break;
}
if (labels) {
TIter iL(labels);
TObjString* lb;
Int_t i = 1;
while ((lb=(TObjString*)iL())) {
h1->GetXaxis()->SetBinLabel(i,lb->String().Data());
i++;
}
}
h1->SetLineColor(this->GetLineColor());
h1->SetFillColor(this->GetFillColor());
h1->SetMarkerColor(this->GetMarkerColor());
h1->SetMarkerStyle(this->GetMarkerStyle());
}
if (pcase >=4 && pcase <=9) {
THashList* labels1 = 0;
THashList* labels2 = 0;
switch (pcase) {
case 4:
h2->GetXaxis()->ImportAttributes(this->GetYaxis());
h2->GetYaxis()->ImportAttributes(this->GetXaxis());
labels1 = GetYaxis()->GetLabels();
labels2 = GetXaxis()->GetLabels();
break;
case 5:
h2->GetXaxis()->ImportAttributes(this->GetXaxis());
h2->GetYaxis()->ImportAttributes(this->GetYaxis());
labels1 = GetXaxis()->GetLabels();
labels2 = GetYaxis()->GetLabels();
break;
case 6:
h2->GetXaxis()->ImportAttributes(this->GetZaxis());
h2->GetYaxis()->ImportAttributes(this->GetXaxis());
labels1 = GetZaxis()->GetLabels();
labels2 = GetXaxis()->GetLabels();
break;
case 7:
h2->GetXaxis()->ImportAttributes(this->GetXaxis());
h2->GetYaxis()->ImportAttributes(this->GetZaxis());
labels1 = GetXaxis()->GetLabels();
labels2 = GetZaxis()->GetLabels();
break;
case 8:
h2->GetXaxis()->ImportAttributes(this->GetZaxis());
h2->GetYaxis()->ImportAttributes(this->GetYaxis());
labels1 = GetZaxis()->GetLabels();
labels2 = GetYaxis()->GetLabels();
break;
case 9:
h2->GetXaxis()->ImportAttributes(this->GetYaxis());
h2->GetYaxis()->ImportAttributes(this->GetZaxis());
labels1 = GetYaxis()->GetLabels();
labels2 = GetZaxis()->GetLabels();
break;
}
if (labels1) {
TIter iL(labels1);
TObjString* lb;
Int_t i = 1;
while ((lb=(TObjString*)iL())) {
h2->GetXaxis()->SetBinLabel(i,lb->String().Data());
i++;
}
}
if (labels2) {
TIter iL(labels2);
TObjString* lb;
Int_t i = 1;
while ((lb=(TObjString*)iL())) {
h2->GetYaxis()->SetBinLabel(i,lb->String().Data());
i++;
}
}
h2->SetLineColor(this->GetLineColor());
h2->SetFillColor(this->GetFillColor());
h2->SetMarkerColor(this->GetMarkerColor());
h2->SetMarkerStyle(this->GetMarkerStyle());
}
delete [] name;
delete [] title;
TH1 *h = h1;
if (h2) h = h2;
if (h == 0) return 0;
Bool_t computeErrors = kFALSE;
if (opt.Contains("e")) {h->Sumw2(); computeErrors = kTRUE;}
if (!fXaxis.TestBit(TAxis::kAxisRange)) {ixmin--; ixmax++;}
if (!fYaxis.TestBit(TAxis::kAxisRange)) {iymin--; iymax++;}
if (!fZaxis.TestBit(TAxis::kAxisRange)) {izmin--; izmax++;}
Float_t cont,e,e1;
Double_t entries = 0;
Double_t newerror = 0;
for (Int_t ixbin=0;ixbin<=1+fXaxis.GetNbins();ixbin++){
Int_t ix = ixbin-ixmin;
if (ix < 0) ix=0; if (ix > nx+1) ix = nx+1;
for (Int_t iybin=0;iybin<=1+fYaxis.GetNbins();iybin++){
Int_t iy = iybin-iymin;
if (iy < 0) iy=0; if (iy > ny+1) iy = ny+1;
for (Int_t izbin=0;izbin<=1+fZaxis.GetNbins();izbin++){
Int_t iz = izbin-izmin;
if (iz < 0) iz=0; if (iz > nz+1) iz = nz+1;
Int_t bin = GetBin(ixbin,iybin,izbin);
cont = GetBinContent(bin);
switch (pcase) {
case 1:
if (iybin < iymin || iybin > iymax) continue;
if (izbin < izmin || izbin > izmax) continue;
e1 = h1->GetBinError(ix);
if (cont) h1->Fill(fXaxis.GetBinCenter(ixbin), cont);
if (computeErrors) {
e = GetBinError(bin);
newerror = TMath::Sqrt(e*e + e1*e1);
h1->SetBinError(ix,newerror);
}
break;
case 2:
if (ixbin < ixmin || ixbin > ixmax) continue;
if (izbin < izmin || izbin > izmax) continue;
e1 = h1->GetBinError(iy);
if (cont) h1->Fill(fYaxis.GetBinCenter(iybin), cont);
if (computeErrors) {
e = GetBinError(bin);
newerror = TMath::Sqrt(e*e + e1*e1);
h1->SetBinError(iy,newerror);
}
break;
case 3:
if (ixbin < ixmin || ixbin > ixmax) continue;
if (iybin < iymin || iybin > iymax) continue;
e1 = h1->GetBinError(iz);
if (cont) h1->Fill(fZaxis.GetBinCenter(izbin), cont);
if (computeErrors) {
e = GetBinError(bin);
newerror = TMath::Sqrt(e*e + e1*e1);
h1->SetBinError(iz,newerror);
}
break;
case 4:
if (izbin < izmin || izbin > izmax) continue;
e1 = h2->GetCellError(iy,ix);
if (cont) h2->Fill(fYaxis.GetBinCenter(iybin),fXaxis.GetBinCenter(ixbin), cont);
if (computeErrors) {
e = GetBinError(bin);
newerror = TMath::Sqrt(e*e + e1*e1);
h2->SetCellError(iy,ix,newerror);
}
break;
case 5:
if (izbin < izmin || izbin > izmax) continue;
e1 = h2->GetCellError(ix,iy);
if (cont) h2->Fill(fXaxis.GetBinCenter(ixbin),fYaxis.GetBinCenter(iybin), cont);
if (computeErrors) {
e = GetBinError(bin);
newerror = TMath::Sqrt(e*e + e1*e1);
h2->SetCellError(ix,iy,newerror);
}
break;
case 6:
if (iybin < iymin || iybin > iymax) continue;
e1 = h2->GetCellError(iz,ix);
if (cont) h2->Fill(fZaxis.GetBinCenter(izbin),fXaxis.GetBinCenter(ixbin), cont);
if (computeErrors) {
e = GetBinError(bin);
newerror = TMath::Sqrt(e*e + e1*e1);
h2->SetCellError(iz,ix,newerror);
}
break;
case 7:
if (iybin < iymin || iybin > iymax) continue;
e1 = h2->GetCellError(ix,iz);
if (cont) h2->Fill(fXaxis.GetBinCenter(ixbin),fZaxis.GetBinCenter(izbin), cont);
if (computeErrors) {
e = GetBinError(bin);
newerror = TMath::Sqrt(e*e + e1*e1);
h2->SetCellError(ix,iz,newerror);
}
break;
case 8:
if (ixbin < ixmin || ixbin > ixmax) continue;
e1 = h2->GetCellError(iz,iy);
if (cont) h2->Fill(fZaxis.GetBinCenter(izbin),fYaxis.GetBinCenter(iybin), cont);
if (computeErrors) {
e = GetBinError(bin);
newerror = TMath::Sqrt(e*e + e1*e1);
h2->SetCellError(iz,iy,newerror);
}
break;
case 9:
if (ixbin < ixmin || ixbin > ixmax) continue;
e1 = h2->GetCellError(iy,iz);
if (cont) h2->Fill(fYaxis.GetBinCenter(iybin),fZaxis.GetBinCenter(izbin), cont);
if (computeErrors) {
e = GetBinError(bin);
newerror = TMath::Sqrt(e*e + e1*e1);
h2->SetCellError(iy,iz,newerror);
}
break;
}
if (cont) {
entries += cont;
}
}
}
}
if (izmin <=1 && izmax >= nz && iymin <=1 && iymax >= ny && ixmin <=1 && ixmax >= nx) h->SetEntries(fEntries);
else h->SetEntries(entries);
return h;
}
TProfile2D *TH3::Project3DProfile(Option_t *option) const
{
TString opt = option; opt.ToLower();
Int_t ixmin = fXaxis.GetFirst();
Int_t ixmax = fXaxis.GetLast();
Int_t iymin = fYaxis.GetFirst();
Int_t iymax = fYaxis.GetLast();
Int_t izmin = fZaxis.GetFirst();
Int_t izmax = fZaxis.GetLast();
Int_t nx = ixmax-ixmin+1;
Int_t ny = iymax-iymin+1;
Int_t nz = izmax-izmin+1;
Int_t pcase = 0;
if (opt.Contains("xy")) pcase = 4;
if (opt.Contains("yx")) pcase = 5;
if (opt.Contains("xz")) pcase = 6;
if (opt.Contains("zx")) pcase = 7;
if (opt.Contains("yz")) pcase = 8;
if (opt.Contains("zy")) pcase = 9;
TProfile2D *p2 = 0;
Int_t nch = strlen(GetName()) +opt.Length() +3;
char *name = new char[nch];
sprintf(name,"%s_p%s",GetName(),option);
nch = strlen(GetTitle()) +opt.Length() +3;
char *title = new char[nch];
sprintf(title,"%s_p%s",GetTitle(),option);
delete gROOT->FindObject(name);
const TArrayD *xbins;
const TArrayD *ybins;
const TArrayD *zbins;
switch (pcase) {
case 4:
xbins = fXaxis.GetXbins();
ybins = fYaxis.GetXbins();
if (xbins->fN == 0 && ybins->fN == 0) {
p2 = new TProfile2D(name,title,ny,fYaxis.GetBinLowEdge(iymin),fYaxis.GetBinUpEdge(iymax)
,nx,fXaxis.GetBinLowEdge(ixmin),fXaxis.GetBinUpEdge(ixmax));
} else if (ybins->fN == 0) {
p2 = new TProfile2D(name,title,ny,fYaxis.GetBinLowEdge(iymin),fYaxis.GetBinUpEdge(iymax)
,nx,&xbins->fArray[ixmin-1]);
} else if (xbins->fN == 0) {
p2 = new TProfile2D(name,title,ny,&ybins->fArray[iymin-1]
,nx,fXaxis.GetBinLowEdge(ixmin),fXaxis.GetBinUpEdge(ixmax));
} else {
p2 = new TProfile2D(name,title,ny,&ybins->fArray[iymin-1],nx,&xbins->fArray[ixmin-1]);
}
break;
case 5:
xbins = fXaxis.GetXbins();
ybins = fYaxis.GetXbins();
if (xbins->fN == 0 && ybins->fN == 0) {
p2 = new TProfile2D(name,title,nx,fXaxis.GetBinLowEdge(ixmin),fXaxis.GetBinUpEdge(ixmax)
,ny,fYaxis.GetBinLowEdge(iymin),fYaxis.GetBinUpEdge(iymax));
} else if (xbins->fN == 0) {
p2 = new TProfile2D(name,title,nx,fXaxis.GetBinLowEdge(ixmin),fXaxis.GetBinUpEdge(ixmax)
,ny,&ybins->fArray[iymin-1]);
} else if (ybins->fN == 0) {
p2 = new TProfile2D(name,title,nx,&xbins->fArray[ixmin-1]
,ny,fYaxis.GetBinLowEdge(iymin),fYaxis.GetBinUpEdge(iymax));
} else {
p2 = new TProfile2D(name,title,nx,&xbins->fArray[ixmin-1],ny,&ybins->fArray[iymin-1]);
}
break;
case 6:
xbins = fXaxis.GetXbins();
zbins = fZaxis.GetXbins();
if (xbins->fN == 0 && zbins->fN == 0) {
p2 = new TProfile2D(name,title,nz,fZaxis.GetBinLowEdge(izmin),fZaxis.GetBinUpEdge(izmax)
,nx,fXaxis.GetBinLowEdge(ixmin),fXaxis.GetBinUpEdge(ixmax));
} else if (zbins->fN == 0) {
p2 = new TProfile2D(name,title,nz,fZaxis.GetBinLowEdge(izmin),fZaxis.GetBinUpEdge(izmax)
,nx,&xbins->fArray[ixmin-1]);
} else if (xbins->fN == 0) {
p2 = new TProfile2D(name,title,nz,&zbins->fArray[izmin-1]
,nx,fXaxis.GetBinLowEdge(ixmin),fXaxis.GetBinUpEdge(ixmax));
} else {
p2 = new TProfile2D(name,title,nz,&zbins->fArray[izmin-1],nx,&xbins->fArray[ixmin-1]);
}
break;
case 7:
xbins = fXaxis.GetXbins();
zbins = fZaxis.GetXbins();
if (xbins->fN == 0 && zbins->fN == 0) {
p2 = new TProfile2D(name,title,nx,fXaxis.GetBinLowEdge(ixmin),fXaxis.GetBinUpEdge(ixmax)
,nz,fZaxis.GetBinLowEdge(izmin),fZaxis.GetBinUpEdge(izmax));
} else if (xbins->fN == 0) {
p2 = new TProfile2D(name,title,nx,fXaxis.GetBinLowEdge(ixmin),fXaxis.GetBinUpEdge(ixmax)
,nz,&zbins->fArray[izmin-1]);
} else if (zbins->fN == 0) {
p2 = new TProfile2D(name,title,nx,&xbins->fArray[ixmin-1]
,nz,fZaxis.GetBinLowEdge(izmin),fZaxis.GetBinUpEdge(izmax));
} else {
p2 = new TProfile2D(name,title,nx,&xbins->fArray[ixmin-1],nz,&zbins->fArray[izmin-1]);
}
break;
case 8:
ybins = fYaxis.GetXbins();
zbins = fZaxis.GetXbins();
if (ybins->fN == 0 && zbins->fN == 0) {
p2 = new TProfile2D(name,title,nz,fZaxis.GetBinLowEdge(izmin),fZaxis.GetBinUpEdge(izmax)
,ny,fYaxis.GetBinLowEdge(iymin),fYaxis.GetBinUpEdge(iymax));
} else if (zbins->fN == 0) {
p2 = new TProfile2D(name,title,nz,fZaxis.GetBinLowEdge(izmin),fZaxis.GetBinUpEdge(izmax)
,ny,&ybins->fArray[iymin-1]);
} else if (ybins->fN == 0) {
p2 = new TProfile2D(name,title,nz,&zbins->fArray[izmin-1]
,ny,fYaxis.GetBinLowEdge(iymin),fYaxis.GetBinUpEdge(iymax));
} else {
p2 = new TProfile2D(name,title,nz,&zbins->fArray[izmin-1],ny,&ybins->fArray[iymin-1]);
}
break;
case 9:
ybins = fYaxis.GetXbins();
zbins = fZaxis.GetXbins();
if (ybins->fN == 0 && zbins->fN == 0) {
p2 = new TProfile2D(name,title,ny,fYaxis.GetBinLowEdge(iymin),fYaxis.GetBinUpEdge(iymax)
,nz,fZaxis.GetBinLowEdge(izmin),fZaxis.GetBinUpEdge(izmax));
} else if (ybins->fN == 0) {
p2 = new TProfile2D(name,title,ny,fYaxis.GetBinLowEdge(iymin),fYaxis.GetBinUpEdge(iymax)
,nz,&zbins->fArray[izmin-1]);
} else if (zbins->fN == 0) {
p2 = new TProfile2D(name,title,ny,&ybins->fArray[iymin-1]
,nz,fZaxis.GetBinLowEdge(izmin),fZaxis.GetBinUpEdge(izmax));
} else {
p2 = new TProfile2D(name,title,ny,&ybins->fArray[iymin-1],nz,&zbins->fArray[izmin-1]);
}
break;
}
delete [] name;
delete [] title;
if (p2 == 0) return 0;
if (!fXaxis.TestBit(TAxis::kAxisRange)) {ixmin--; ixmax++;}
if (!fYaxis.TestBit(TAxis::kAxisRange)) {iymin--; izmax++;}
if (!fZaxis.TestBit(TAxis::kAxisRange)) {izmin--; izmax++;}
Double_t cont,e,e2;
Double_t entries = 0;
for (Int_t ixbin=0;ixbin<=1+fXaxis.GetNbins();ixbin++){
Int_t ix = ixbin-ixmin;
if (ix < 0) ix=0; if (ix > nx+1) ix = nx+1;
for (Int_t iybin=0;iybin<=1+fYaxis.GetNbins();iybin++){
Int_t iy = iybin-iymin;
if (iy < 0) iy=0; if (iy > ny+1) iy = ny+1;
for (Int_t izbin=0;izbin<=1+fZaxis.GetNbins();izbin++){
Int_t iz = izbin-izmin;
if (iz < 0) iz=0; if (iz > nz+1) iz = nz+1;
Int_t bin = GetBin(ixbin,iybin,izbin);
cont = GetBinContent(bin);
e = GetBinError(bin);
e2 = e*e;
switch (pcase) {
case 4:
if (izbin < izmin || izbin > izmax) continue;
if (cont) p2->Fill(fYaxis.GetBinCenter(iybin),fXaxis.GetBinCenter(ixbin),fZaxis.GetBinCenter(izbin), e2);
break;
case 5:
if (izbin < izmin || izbin > izmax) continue;
if (cont) p2->Fill(fXaxis.GetBinCenter(ixbin),fYaxis.GetBinCenter(iybin),fZaxis.GetBinCenter(izbin), e2);
break;
case 6:
if (iybin < iymin || iybin > iymax) continue;
if (cont) p2->Fill(fZaxis.GetBinCenter(izbin),fXaxis.GetBinCenter(ixbin),fYaxis.GetBinCenter(izbin), e2);
break;
case 7:
if (iybin < iymin || iybin > iymax) continue;
if (cont) p2->Fill(fXaxis.GetBinCenter(ixbin),fZaxis.GetBinCenter(izbin),fYaxis.GetBinCenter(izbin), e2);
break;
case 8:
if (ixbin < ixmin || ixbin > ixmax) continue;
if (cont) p2->Fill(fZaxis.GetBinCenter(izbin),fYaxis.GetBinCenter(iybin),fXaxis.GetBinCenter(izbin), e2);
break;
case 9:
if (ixbin < ixmin || ixbin > ixmax) continue;
if (cont) p2->Fill(fYaxis.GetBinCenter(iybin),fZaxis.GetBinCenter(izbin),fXaxis.GetBinCenter(izbin), e2);
break;
}
if (cont) {
entries += cont;
}
}
}
}
if (iymin <=1 && iymax >= ny && ixmin <=1 && ixmax >= nx) p2->SetEntries(fEntries);
else p2->SetEntries(entries);
return p2;
}
void TH3::PutStats(Double_t *stats)
{
TH1::PutStats(stats);
fTsumwy = stats[4];
fTsumwy2 = stats[5];
fTsumwxy = stats[6];
fTsumwz = stats[7];
fTsumwz2 = stats[8];
fTsumwxz = stats[9];
fTsumwyz = stats[10];
}
void TH3::Reset(Option_t *option)
{
TH1::Reset(option);
TString opt = option;
opt.ToUpper();
if (opt.Contains("ICE")) return;
fTsumwy = 0;
fTsumwy2 = 0;
fTsumwxy = 0;
fTsumwz = 0;
fTsumwz2 = 0;
fTsumwxz = 0;
fTsumwyz = 0;
}
void TH3::Streamer(TBuffer &R__b)
{
if (R__b.IsReading()) {
UInt_t R__s, R__c;
Version_t R__v = R__b.ReadVersion(&R__s, &R__c);
if (R__v > 2) {
TH3::Class()->ReadBuffer(R__b, this, R__v, R__s, R__c);
return;
}
TH1::Streamer(R__b);
TAtt3D::Streamer(R__b);
R__b.CheckByteCount(R__s, R__c, TH3::IsA());
} else {
TH3::Class()->WriteBuffer(R__b,this);
}
}
ClassImp(TH3C)
TH3C::TH3C(): TH3()
{
}
TH3C::~TH3C()
{
}
TH3C::TH3C(const char *name,const char *title,Int_t nbinsx,Double_t xlow,Double_t xup
,Int_t nbinsy,Double_t ylow,Double_t yup
,Int_t nbinsz,Double_t zlow,Double_t zup)
:TH3(name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
{
TArrayC::Set(fNcells);
}
TH3C::TH3C(const char *name,const char *title,Int_t nbinsx,const Float_t *xbins
,Int_t nbinsy,const Float_t *ybins
,Int_t nbinsz,const Float_t *zbins)
:TH3(name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
{
TArrayC::Set(fNcells);
}
TH3C::TH3C(const char *name,const char *title,Int_t nbinsx,const Double_t *xbins
,Int_t nbinsy,const Double_t *ybins
,Int_t nbinsz,const Double_t *zbins)
:TH3(name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
{
TArrayC::Set(fNcells);
}
TH3C::TH3C(const TH3C &h3c) : TH3(), TArrayC()
{
((TH3C&)h3c).Copy(*this);
}
void TH3C::AddBinContent(Int_t bin)
{
if (fArray[bin] < 127) fArray[bin]++;
}
void TH3C::AddBinContent(Int_t bin, Double_t w)
{
Int_t newval = fArray[bin] + Int_t(w);
if (newval > -128 && newval < 128) {fArray[bin] = Char_t(newval); return;}
if (newval < -127) fArray[bin] = -127;
if (newval > 127) fArray[bin] = 127;
}
void TH3C::Copy(TObject &newth3) const
{
TH3::Copy((TH3C&)newth3);
TArrayC::Copy((TH3C&)newth3);
}
TH1 *TH3C::DrawCopy(Option_t *option) const
{
TString opt = option;
opt.ToLower();
if (gPad && !opt.Contains("same")) gPad->Clear();
TH3C *newth3 = (TH3C*)Clone();
newth3->SetDirectory(0);
newth3->SetBit(kCanDelete);
newth3->AppendPad(option);
return newth3;
}
Double_t TH3C::GetBinContent(Int_t bin) const
{
if (fBuffer) ((TH3C*)this)->BufferEmpty();
if (bin < 0) bin = 0;
if (bin >= fNcells) bin = fNcells-1;
if (!fArray) return 0;
return Double_t (fArray[bin]);
}
void TH3C::Reset(Option_t *option)
{
TH3::Reset(option);
TArrayC::Reset();
}
void TH3C::SetBinsLength(Int_t n)
{
if (n < 0) n = (fXaxis.GetNbins()+2)*(fYaxis.GetNbins()+2)*(fZaxis.GetNbins()+2);
fNcells = n;
TArrayC::Set(n);
}
void TH3C::SetBinContent(Int_t bin, Double_t content)
{
if (bin < 0) return;
if (bin >= fNcells) return;
fArray[bin] = Char_t (content);
fEntries++;
}
void TH3::SetShowProjection(const char *option,Int_t nbins)
{
GetPainter();
if (fPainter) fPainter->SetShowProjection(option,nbins);
}
void TH3C::Streamer(TBuffer &R__b)
{
if (R__b.IsReading()) {
UInt_t R__s, R__c;
TFile *file = (TFile*)R__b.GetParent();
if (file && file->GetVersion() < 22300) return;
Version_t R__v = R__b.ReadVersion(&R__s, &R__c);
if (R__v > 2) {
TH3C::Class()->ReadBuffer(R__b, this, R__v, R__s, R__c);
return;
}
if (R__v < 2) {
R__b.ReadVersion();
TH1::Streamer(R__b);
TArrayC::Streamer(R__b);
R__b.ReadVersion(&R__s, &R__c);
TAtt3D::Streamer(R__b);
} else {
TH3::Streamer(R__b);
TArrayC::Streamer(R__b);
R__b.CheckByteCount(R__s, R__c, TH3C::IsA());
}
} else {
TH3C::Class()->WriteBuffer(R__b,this);
}
}
TH3C& TH3C::operator=(const TH3C &h1)
{
if (this != &h1) ((TH3C&)h1).Copy(*this);
return *this;
}
TH3C operator*(Float_t c1, TH3C &h1)
{
TH3C hnew = h1;
hnew.Scale(c1);
hnew.SetDirectory(0);
return hnew;
}
TH3C operator+(TH3C &h1, TH3C &h2)
{
TH3C hnew = h1;
hnew.Add(&h2,1);
hnew.SetDirectory(0);
return hnew;
}
TH3C operator-(TH3C &h1, TH3C &h2)
{
TH3C hnew = h1;
hnew.Add(&h2,-1);
hnew.SetDirectory(0);
return hnew;
}
TH3C operator*(TH3C &h1, TH3C &h2)
{
TH3C hnew = h1;
hnew.Multiply(&h2);
hnew.SetDirectory(0);
return hnew;
}
TH3C operator/(TH3C &h1, TH3C &h2)
{
TH3C hnew = h1;
hnew.Divide(&h2);
hnew.SetDirectory(0);
return hnew;
}
ClassImp(TH3S)
TH3S::TH3S(): TH3()
{
}
TH3S::~TH3S()
{
}
TH3S::TH3S(const char *name,const char *title,Int_t nbinsx,Double_t xlow,Double_t xup
,Int_t nbinsy,Double_t ylow,Double_t yup
,Int_t nbinsz,Double_t zlow,Double_t zup)
:TH3(name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
{
TH3S::Set(fNcells);
}
TH3S::TH3S(const char *name,const char *title,Int_t nbinsx,const Float_t *xbins
,Int_t nbinsy,const Float_t *ybins
,Int_t nbinsz,const Float_t *zbins)
:TH3(name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
{
TH3S::Set(fNcells);
}
TH3S::TH3S(const char *name,const char *title,Int_t nbinsx,const Double_t *xbins
,Int_t nbinsy,const Double_t *ybins
,Int_t nbinsz,const Double_t *zbins)
:TH3(name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
{
TH3S::Set(fNcells);
}
TH3S::TH3S(const TH3S &h3s) : TH3(), TArrayS()
{
((TH3S&)h3s).Copy(*this);
}
void TH3S::AddBinContent(Int_t bin)
{
if (fArray[bin] < 32767) fArray[bin]++;
}
void TH3S::AddBinContent(Int_t bin, Double_t w)
{
Int_t newval = fArray[bin] + Int_t(w);
if (newval > -32768 && newval < 32768) {fArray[bin] = Short_t(newval); return;}
if (newval < -32767) fArray[bin] = -32767;
if (newval > 32767) fArray[bin] = 32767;
}
void TH3S::Copy(TObject &newth3) const
{
TH3::Copy((TH3S&)newth3);
TArrayS::Copy((TH3S&)newth3);
}
TH1 *TH3S::DrawCopy(Option_t *option) const
{
TString opt = option;
opt.ToLower();
if (gPad && !opt.Contains("same")) gPad->Clear();
TH3S *newth3 = (TH3S*)Clone();
newth3->SetDirectory(0);
newth3->SetBit(kCanDelete);
newth3->AppendPad(option);
return newth3;
}
Double_t TH3S::GetBinContent(Int_t bin) const
{
if (fBuffer) ((TH3S*)this)->BufferEmpty();
if (bin < 0) bin = 0;
if (bin >= fNcells) bin = fNcells-1;
if (!fArray) return 0;
return Double_t (fArray[bin]);
}
void TH3S::Reset(Option_t *option)
{
TH3::Reset(option);
TArrayS::Reset();
}
void TH3S::SetBinContent(Int_t bin, Double_t content)
{
if (bin < 0) return;
if (bin >= fNcells) return;
fArray[bin] = Short_t (content);
fEntries++;
}
void TH3S::SetBinsLength(Int_t n)
{
if (n < 0) n = (fXaxis.GetNbins()+2)*(fYaxis.GetNbins()+2)*(fZaxis.GetNbins()+2);
fNcells = n;
TArrayS::Set(n);
}
void TH3S::Streamer(TBuffer &R__b)
{
if (R__b.IsReading()) {
UInt_t R__s, R__c;
TFile *file = (TFile*)R__b.GetParent();
if (file && file->GetVersion() < 22300) return;
Version_t R__v = R__b.ReadVersion(&R__s, &R__c);
if (R__v > 2) {
TH3S::Class()->ReadBuffer(R__b, this, R__v, R__s, R__c);
return;
}
if (R__v < 2) {
R__b.ReadVersion();
TH1::Streamer(R__b);
TArrayS::Streamer(R__b);
R__b.ReadVersion(&R__s, &R__c);
TAtt3D::Streamer(R__b);
} else {
TH3::Streamer(R__b);
TArrayS::Streamer(R__b);
R__b.CheckByteCount(R__s, R__c, TH3S::IsA());
}
} else {
TH3S::Class()->WriteBuffer(R__b,this);
}
}
TH3S& TH3S::operator=(const TH3S &h1)
{
if (this != &h1) ((TH3S&)h1).Copy(*this);
return *this;
}
TH3S operator*(Float_t c1, TH3S &h1)
{
TH3S hnew = h1;
hnew.Scale(c1);
hnew.SetDirectory(0);
return hnew;
}
TH3S operator+(TH3S &h1, TH3S &h2)
{
TH3S hnew = h1;
hnew.Add(&h2,1);
hnew.SetDirectory(0);
return hnew;
}
TH3S operator-(TH3S &h1, TH3S &h2)
{
TH3S hnew = h1;
hnew.Add(&h2,-1);
hnew.SetDirectory(0);
return hnew;
}
TH3S operator*(TH3S &h1, TH3S &h2)
{
TH3S hnew = h1;
hnew.Multiply(&h2);
hnew.SetDirectory(0);
return hnew;
}
TH3S operator/(TH3S &h1, TH3S &h2)
{
TH3S hnew = h1;
hnew.Divide(&h2);
hnew.SetDirectory(0);
return hnew;
}
ClassImp(TH3I)
TH3I::TH3I(): TH3()
{
}
TH3I::~TH3I()
{
}
TH3I::TH3I(const char *name,const char *title,Int_t nbinsx,Double_t xlow,Double_t xup
,Int_t nbinsy,Double_t ylow,Double_t yup
,Int_t nbinsz,Double_t zlow,Double_t zup)
:TH3(name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
{
TH3I::Set(fNcells);
}
TH3I::TH3I(const char *name,const char *title,Int_t nbinsx,const Float_t *xbins
,Int_t nbinsy,const Float_t *ybins
,Int_t nbinsz,const Float_t *zbins)
:TH3(name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
{
TArrayI::Set(fNcells);
}
TH3I::TH3I(const char *name,const char *title,Int_t nbinsx,const Double_t *xbins
,Int_t nbinsy,const Double_t *ybins
,Int_t nbinsz,const Double_t *zbins)
:TH3(name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
{
TArrayI::Set(fNcells);
}
TH3I::TH3I(const TH3I &h3i) : TH3(), TArrayI()
{
((TH3I&)h3i).Copy(*this);
}
void TH3I::AddBinContent(Int_t bin)
{
if (fArray[bin] < 2147483647) fArray[bin]++;
}
void TH3I::AddBinContent(Int_t bin, Double_t w)
{
Int_t newval = fArray[bin] + Int_t(w);
if (newval > -2147483647 && newval < 2147483647) {fArray[bin] = Int_t(newval); return;}
if (newval < -2147483647) fArray[bin] = -2147483647;
if (newval > 2147483647) fArray[bin] = 2147483647;
}
void TH3I::Copy(TObject &newth3) const
{
TH3::Copy((TH3I&)newth3);
TArrayI::Copy((TH3I&)newth3);
}
TH1 *TH3I::DrawCopy(Option_t *option) const
{
TString opt = option;
opt.ToLower();
if (gPad && !opt.Contains("same")) gPad->Clear();
TH3I *newth3 = (TH3I*)Clone();
newth3->SetDirectory(0);
newth3->SetBit(kCanDelete);
newth3->AppendPad(option);
return newth3;
}
Double_t TH3I::GetBinContent(Int_t bin) const
{
if (fBuffer) ((TH3I*)this)->BufferEmpty();
if (bin < 0) bin = 0;
if (bin >= fNcells) bin = fNcells-1;
if (!fArray) return 0;
return Double_t (fArray[bin]);
}
void TH3I::Reset(Option_t *option)
{
TH3::Reset(option);
TArrayI::Reset();
}
void TH3I::SetBinContent(Int_t bin, Double_t content)
{
if (bin < 0) return;
if (bin >= fNcells) return;
fArray[bin] = Int_t (content);
fEntries++;
}
void TH3I::SetBinsLength(Int_t n)
{
if (n < 0) n = (fXaxis.GetNbins()+2)*(fYaxis.GetNbins()+2)*(fZaxis.GetNbins()+2);
fNcells = n;
TArrayI::Set(n);
}
TH3I& TH3I::operator=(const TH3I &h1)
{
if (this != &h1) ((TH3I&)h1).Copy(*this);
return *this;
}
TH3I operator*(Float_t c1, TH3I &h1)
{
TH3I hnew = h1;
hnew.Scale(c1);
hnew.SetDirectory(0);
return hnew;
}
TH3I operator+(TH3I &h1, TH3I &h2)
{
TH3I hnew = h1;
hnew.Add(&h2,1);
hnew.SetDirectory(0);
return hnew;
}
TH3I operator-(TH3I &h1, TH3I &h2)
{
TH3I hnew = h1;
hnew.Add(&h2,-1);
hnew.SetDirectory(0);
return hnew;
}
TH3I operator*(TH3I &h1, TH3I &h2)
{
TH3I hnew = h1;
hnew.Multiply(&h2);
hnew.SetDirectory(0);
return hnew;
}
TH3I operator/(TH3I &h1, TH3I &h2)
{
TH3I hnew = h1;
hnew.Divide(&h2);
hnew.SetDirectory(0);
return hnew;
}
ClassImp(TH3F)
TH3F::TH3F(): TH3()
{
}
TH3F::~TH3F()
{
}
TH3F::TH3F(const char *name,const char *title,Int_t nbinsx,Double_t xlow,Double_t xup
,Int_t nbinsy,Double_t ylow,Double_t yup
,Int_t nbinsz,Double_t zlow,Double_t zup)
:TH3(name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
{
TArrayF::Set(fNcells);
}
TH3F::TH3F(const char *name,const char *title,Int_t nbinsx,const Float_t *xbins
,Int_t nbinsy,const Float_t *ybins
,Int_t nbinsz,const Float_t *zbins)
:TH3(name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
{
TArrayF::Set(fNcells);
}
TH3F::TH3F(const char *name,const char *title,Int_t nbinsx,const Double_t *xbins
,Int_t nbinsy,const Double_t *ybins
,Int_t nbinsz,const Double_t *zbins)
:TH3(name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
{
TArrayF::Set(fNcells);
}
TH3F::TH3F(const TH3F &h3f) : TH3(), TArrayF()
{
((TH3F&)h3f).Copy(*this);
}
void TH3F::Copy(TObject &newth3) const
{
TH3::Copy((TH3F&)newth3);
TArrayF::Copy((TH3F&)newth3);
}
TH1 *TH3F::DrawCopy(Option_t *option) const
{
TString opt = option;
opt.ToLower();
if (gPad && !opt.Contains("same")) gPad->Clear();
TH3F *newth3 = (TH3F*)Clone();
newth3->SetDirectory(0);
newth3->SetBit(kCanDelete);
newth3->AppendPad(option);
return newth3;
}
Double_t TH3F::GetBinContent(Int_t bin) const
{
if (fBuffer) ((TH3F*)this)->BufferEmpty();
if (bin < 0) bin = 0;
if (bin >= fNcells) bin = fNcells-1;
if (!fArray) return 0;
return Double_t (fArray[bin]);
}
void TH3F::Reset(Option_t *option)
{
TH3::Reset(option);
TArrayF::Reset();
}
void TH3F::SetBinContent(Int_t bin, Double_t content)
{
if (bin < 0) return;
if (bin >= fNcells) return;
fArray[bin] = Float_t (content);
fEntries++;
}
void TH3F::SetBinsLength(Int_t n)
{
if (n < 0) n = (fXaxis.GetNbins()+2)*(fYaxis.GetNbins()+2)*(fZaxis.GetNbins()+2);
fNcells = n;
TArrayF::Set(n);
}
void TH3F::Streamer(TBuffer &R__b)
{
if (R__b.IsReading()) {
UInt_t R__s, R__c;
TFile *file = (TFile*)R__b.GetParent();
if (file && file->GetVersion() < 22300) return;
Version_t R__v = R__b.ReadVersion(&R__s, &R__c);
if (R__v > 2) {
TH3F::Class()->ReadBuffer(R__b, this, R__v, R__s, R__c);
return;
}
if (R__v < 2) {
R__b.ReadVersion();
TH1::Streamer(R__b);
TArrayF::Streamer(R__b);
R__b.ReadVersion(&R__s, &R__c);
TAtt3D::Streamer(R__b);
} else {
TH3::Streamer(R__b);
TArrayF::Streamer(R__b);
R__b.CheckByteCount(R__s, R__c, TH3F::IsA());
}
} else {
TH3F::Class()->WriteBuffer(R__b,this);
}
}
TH3F& TH3F::operator=(const TH3F &h1)
{
if (this != &h1) ((TH3F&)h1).Copy(*this);
return *this;
}
TH3F operator*(Float_t c1, TH3F &h1)
{
TH3F hnew = h1;
hnew.Scale(c1);
hnew.SetDirectory(0);
return hnew;
}
TH3F operator+(TH3F &h1, TH3F &h2)
{
TH3F hnew = h1;
hnew.Add(&h2,1);
hnew.SetDirectory(0);
return hnew;
}
TH3F operator-(TH3F &h1, TH3F &h2)
{
TH3F hnew = h1;
hnew.Add(&h2,-1);
hnew.SetDirectory(0);
return hnew;
}
TH3F operator*(TH3F &h1, TH3F &h2)
{
TH3F hnew = h1;
hnew.Multiply(&h2);
hnew.SetDirectory(0);
return hnew;
}
TH3F operator/(TH3F &h1, TH3F &h2)
{
TH3F hnew = h1;
hnew.Divide(&h2);
hnew.SetDirectory(0);
return hnew;
}
ClassImp(TH3D)
TH3D::TH3D(): TH3()
{
}
TH3D::~TH3D()
{
}
TH3D::TH3D(const char *name,const char *title,Int_t nbinsx,Double_t xlow,Double_t xup
,Int_t nbinsy,Double_t ylow,Double_t yup
,Int_t nbinsz,Double_t zlow,Double_t zup)
:TH3(name,title,nbinsx,xlow,xup,nbinsy,ylow,yup,nbinsz,zlow,zup)
{
TArrayD::Set(fNcells);
}
TH3D::TH3D(const char *name,const char *title,Int_t nbinsx,const Float_t *xbins
,Int_t nbinsy,const Float_t *ybins
,Int_t nbinsz,const Float_t *zbins)
:TH3(name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
{
TArrayD::Set(fNcells);
}
TH3D::TH3D(const char *name,const char *title,Int_t nbinsx,const Double_t *xbins
,Int_t nbinsy,const Double_t *ybins
,Int_t nbinsz,const Double_t *zbins)
:TH3(name,title,nbinsx,xbins,nbinsy,ybins,nbinsz,zbins)
{
TArrayD::Set(fNcells);
}
TH3D::TH3D(const TH3D &h3d) : TH3(), TArrayD()
{
((TH3D&)h3d).Copy(*this);
}
void TH3D::Copy(TObject &newth3) const
{
TH3::Copy((TH3D&)newth3);
TArrayD::Copy((TH3D&)newth3);
}
TH1 *TH3D::DrawCopy(Option_t *option) const
{
TString opt = option;
opt.ToLower();
if (gPad && !opt.Contains("same")) gPad->Clear();
TH3D *newth3 = (TH3D*)Clone();
newth3->SetDirectory(0);
newth3->SetBit(kCanDelete);
newth3->AppendPad(option);
return newth3;
}
Double_t TH3D::GetBinContent(Int_t bin) const
{
if (fBuffer) ((TH3D*)this)->BufferEmpty();
if (bin < 0) bin = 0;
if (bin >= fNcells) bin = fNcells-1;
if (!fArray) return 0;
return Double_t (fArray[bin]);
}
void TH3D::Reset(Option_t *option)
{
TH3::Reset(option);
TArrayD::Reset();
}
void TH3D::SetBinContent(Int_t bin, Double_t content)
{
if (bin < 0) return;
if (bin >= fNcells) return;
fArray[bin] = Double_t (content);
fEntries++;
}
void TH3D::SetBinsLength(Int_t n)
{
if (n < 0) n = (fXaxis.GetNbins()+2)*(fYaxis.GetNbins()+2)*(fZaxis.GetNbins()+2);
fNcells = n;
TArrayD::Set(n);
}
void TH3D::Streamer(TBuffer &R__b)
{
if (R__b.IsReading()) {
UInt_t R__s, R__c;
TFile *file = (TFile*)R__b.GetParent();
if (file && file->GetVersion() < 22300) return;
Version_t R__v = R__b.ReadVersion(&R__s, &R__c);
if (R__v > 2) {
TH3D::Class()->ReadBuffer(R__b, this, R__v, R__s, R__c);
return;
}
if (R__v < 2) {
R__b.ReadVersion();
TH1::Streamer(R__b);
TArrayD::Streamer(R__b);
R__b.ReadVersion(&R__s, &R__c);
TAtt3D::Streamer(R__b);
} else {
TH3::Streamer(R__b);
TArrayD::Streamer(R__b);
R__b.CheckByteCount(R__s, R__c, TH3D::IsA());
}
} else {
TH3D::Class()->WriteBuffer(R__b,this);
}
}
TH3D& TH3D::operator=(const TH3D &h1)
{
if (this != &h1) ((TH3D&)h1).Copy(*this);
return *this;
}
TH3D operator*(Float_t c1, TH3D &h1)
{
TH3D hnew = h1;
hnew.Scale(c1);
hnew.SetDirectory(0);
return hnew;
}
TH3D operator+(TH3D &h1, TH3D &h2)
{
TH3D hnew = h1;
hnew.Add(&h2,1);
hnew.SetDirectory(0);
return hnew;
}
TH3D operator-(TH3D &h1, TH3D &h2)
{
TH3D hnew = h1;
hnew.Add(&h2,-1);
hnew.SetDirectory(0);
return hnew;
}
TH3D operator*(TH3D &h1, TH3D &h2)
{
TH3D hnew = h1;
hnew.Multiply(&h2);
hnew.SetDirectory(0);
return hnew;
}
TH3D operator/(TH3D &h1, TH3D &h2)
{
TH3D hnew = h1;
hnew.Divide(&h2);
hnew.SetDirectory(0);
return hnew;
}
ROOT page - Class index - Class Hierarchy - Top of the page
This page has been automatically generated. If you have any comments or suggestions about the page layout send a mail to ROOT support, or contact the developers with any questions or problems regarding ROOT.