//
//-----------------------------------------------------------------------------
+#include <TClonesArray.h>
+#include <TMath.h>
+#include <TRandom.h>
+
#include "AliH2F.h"
-#include "TClonesArray.h"
-#include "TRandom.h"
ClassImp(AliH2F)
for (Int_t i = 0 ;i<dimx;i++)
for (Int_t j = 0 ;j<dimy;j++)
{
- SetCellContent(i,j,0);
- SetCellError(i,j,0);
+ SetBinContent(GetBin(i,j),0);
+ SetBinError(GetBin(i,j),0);
}
}
for (Int_t j = 0 ;j<dimy;j++)
{
Float_t noise = gRandom->Gaus(0,sn);
- Float_t oldv =GetCellContent(i,j);
- Float_t olds =GetCellError(i,j);
+ Float_t oldv =GetBinContent(GetBin(i,j));
+ Float_t olds =GetBinError(GetBin(i,j));
if (noise >0)
{
- SetCellContent(i,j,noise+oldv);
- SetCellError(i,j,TMath::Sqrt((noise*noise+olds*olds)));
+ SetBinContent(GetBin(i,j),noise+oldv);
+ SetBinError(GetBin(i,j),TMath::Sqrt((noise*noise+olds*olds)));
}
}
}
Float_t dx2 = (x2-x)*(x2-x);
Float_t dy2 = (y2-y)*(y2-y);
Float_t amp =max*exp(-(dx2/(2*sx*sx)+dy2/(2*sy*sy)));
- //Float_t oldv =GetCellContent(i+1,j+1);
- // SetCellContent(i+1,j+1,amp+oldv);
+ //Float_t oldv =GetBinContent(GetBin(i+1,j+1));
+ // SetBinContent(GetBin(i+1,j+1),amp+oldv);
Fill(x2,y2,amp);
}
}
for (Int_t i = 0 ;i<=dimx;i++)
for (Int_t j = 0 ;j<=dimy;j++)
{
- Float_t oldv =GetCellContent(i,j);
+ Float_t oldv =GetBinContent(GetBin(i,j));
if (oldv <threshold)
- SetCellContent(i,j,0);
+ SetBinContent(GetBin(i,j),0);
}
}
for (Int_t i = 0 ;i<=dimx;i++)
for (Int_t j = 0 ;j<=dimy;j++)
{
- Float_t oldv =GetCellContent(i,j);
+ Float_t oldv =GetBinContent(GetBin(i,j));
oldv=(Int_t)oldv;
- SetCellContent(i,j,oldv);
+ SetBinContent(GetBin(i,j),oldv);
}
}
// sub->SetBinContent(index2,val);
// Float_t err = GetBinError(index1);
//sub->SetBinError(index2,GetBinError(index1));
- sub->SetCellContent(i,j,val);
+ sub->SetBinContent(GetBin(i,j),val);
}
return sub;
}