53 TNamed(name, title), fNdimensions(dim), fAxes(dim), fBrowsables(dim),
54 fEntries(0), fTsumw(0), fTsumw2(-1.), fTsumwx(dim), fTsumwx2(dim),
55 fIntegral(0), fIntegralStatus(kNoInt)
58 TAxis* axis =
new TAxis(nbins[i], xmin ? xmin[i] : 0., xmax ? xmax[i] : 1.);
84 Int_t chunkSize = 1024 * 16;
86 chunkSize = ((
const THnSparse*)
this)->GetChunkSize();
88 ret->
Init(name, title, axes, keepTargetAxis, chunkSize);
103 const TAxis* axis = 0;
106 while ((axis = (
TAxis*)iAxis())) {
117 Int_t nBins = binLast - binFirst + 1;
146 Bool_t keepTargetAxis )
const {
147 const int ndim = axes->
GetSize();
152 hist =
new TH1D(name, title, 1, 0., 1.);
154 hist =
new TH2D(name, title, 1, 0., 1., 1, 0., 1.);
156 hist =
new TH3D(name, title, 1, 0., 1., 1, 0., 1., 1, 0., 1.);
158 Error(
"CreateHist",
"Cannot create histogram %s with %d dimensions!", name, ndim);
163 for (
Int_t d = 0; d < ndim; ++d) {
170 Int_t nBins = binLast - binFirst + 1;
204 int nbins[3] = {0,0,0};
208 for (
int i = 0; i < ndim; ++i) {
210 minRange[i] = axis[i]->
GetXmin();
211 maxRange[i] = axis[i]->
GetXmax();
220 if (cname[0] ==
'T' && cname[1] ==
'H' 221 && cname[2] >=
'1' && cname[2] <=
'3' && cname[4] == 0) {
223 #define R__THNBCASE(TAG) \ 225 s = new _NAME2_(THnSparse,TAG)(name, title, ndim, nbins, \ 226 minRange, maxRange, chunkSize); \ 228 s = new _NAME2_(THn,TAG)(name, title, ndim, nbins, \ 229 minRange, maxRange); \ 243 ::Warning(
"THnSparse::CreateHnAny",
"Unknown Type of Histogram");
247 for (
int i = 0; i < ndim; ++i) {
254 ::Warning(
"THnSparse::CreateHnAny",
"Unknown Type of Histogram");
273 if (sparse) type = hn->IsA();
283 hn->
Error(
"CreateHnAny",
"Type %s not implemented; please inform the ROOT team!",
290 if (!sparse) type = hn->IsA();
300 hn->
Error(
"CreateHnAny",
"Type THnSparse with Long64_t bins is not available!");
308 hn->
Error(
"CreateHnAny",
"Unhandled type %s, not deriving from THn nor THnSparse!",
313 hn->
Error(
"CreateHnAny",
"Unhandled type %s, please inform the ROOT team!",
335 for (
int i = 0; i <
nbins; ++i) {
338 if (!value && !error)
continue;
376 if (!fitOption.
Chi2) fitOption.
Like =
true;
403 const Int_t nStaticBins = 40;
404 Int_t bin[nStaticBins];
438 if (coord[i] < min || coord[i] > max)
461 for (
Int_t d = 0; d < ndim; ++d) {
468 if (posInsert ==
kNPOS) {
469 title +=
" projection ";
470 for (
Int_t d = 0; d < ndim; ++d)
473 for (
Int_t d = ndim - 1; d >= 0; --d) {
476 title.
Insert(posInsert,
", ");
478 title.
Insert(posInsert,
" projection ");
482 for (
Int_t d = 0; d < ndim; ++d) {
491 Bool_t ignoreTargetRange = (option && (strchr(option,
'A') || strchr(option,
'a')));
492 Bool_t keepTargetAxis = ignoreTargetRange || (option && (strchr(option,
'O') || strchr(option,
'o')));
493 if (ignoreTargetRange) {
494 hadRange =
new Bool_t[ndim];
495 for (
Int_t d = 0; d < ndim; ++d){
503 ret = hn =
CloneEmpty(name, title, &newaxes, keepTargetAxis);
505 ret = hist =
CreateHist(name, title, &newaxes, keepTargetAxis);
507 if (keepTargetAxis) {
510 for (
Int_t d = 0; d < ndim; ++d) {
521 Bool_t wantErrors = haveErrors || (option && (strchr(option,
'E') || strchr(option,
'e')));
528 while ((myLinBin = iter.
Next()) >= 0) {
531 for (
Int_t d = 0; d < ndim; ++d) {
536 if (binOffset > 0) --binOffset;
537 bins[d] -= binOffset;
543 if (ndim == 1) targetLinBin = bins[0];
544 else if (ndim == 2) targetLinBin = hist->
GetBin(bins[0], bins[1]);
545 else if (ndim == 3) targetLinBin = hist->
GetBin(bins[0], bins[1], bins[2]);
596 for (
Int_t d = 0; d < ndim; ++d)
618 while ((i = iter.
Next()) >= 0) {
637 Warning(
"RebinnedAdd",
"Different number of dimensions, cannot carry out operation on the histograms");
659 while ((i = iter.
Next(coord)) >= 0) {
729 while ((addMeObj = iter())) {
738 while ((addMeObj = iter())) {
741 Error(
"Merge",
"Object named %s is not THnBase! Skipping it.",
766 if (wantErrors)
Sumw2();
774 while ((i = iter.
Next(coord)) >= 0) {
810 if (wantErrors)
Sumw2();
815 while ((i = iter.
Next(coord)) >= 0) {
860 if (wantErrors)
Sumw2();
868 while ((i = iter.
Next(coord)) >= 0) {
878 Warning(
"Divide(h)",
"Histogram h has empty bins - division by zero! Setting bin to 0.");
887 Double_t err = (err1 * err1 + err2 * err2) / (b22 * b22);
915 Error(
"Divide",
"Coefficient of dividing histogram cannot be zero");
937 while ((i = iter.
Next(coord)) >= 0) {
947 Warning(
"Divide(h1, h2)",
"Histogram h2 has empty bins - division by zero! Setting bin to 0.");
962 errSq =
TMath::Abs( ( (1. - 2.*w) * err1 * err1 + err2 * err2 ) / (v2 * v2) );
970 errSq = c1 * c2 * (err1 * err1 + err2 * err2) / (b22 * b22);
989 Warning(tag,
"Different number of dimensions, cannot carry out operation on the histograms");
994 Warning(tag,
"Different number of bins on axis %i, cannot carry out operation on the histograms", dim);
1026 if (endHistTitle >= 0) {
1028 Int_t posTitle = endHistTitle + 1;
1031 while (posTitle > 0 && posTitle < lenTitle && dim <
fNdimensions){
1034 axisTitle.
ReplaceAll(
"#semicolon", 10,
";", 1);
1038 posTitle = endTitle + 1;
1076 for (
Int_t d = 0; d < ndim; ++d)
1077 name +=
Form(
"_%d", group[d]);
1082 if (posInsert ==
kNPOS) {
1084 for (
Int_t d = 0; d < ndim; ++d)
1085 title +=
Form(
"{%d}", group[d]);
1087 for (
Int_t d = ndim - 1; d >= 0; --d)
1088 title.
Insert(posInsert,
Form(
"{%d}", group[d]));
1089 title.
Insert(posInsert,
" rebin ");
1094 for (
Int_t d = 0; d < ndim; ++d) {
1098 Int_t newbins = (newaxis->
GetNbins() + group[d] - 1) / group[d];
1102 for (
Int_t i = 0; i < newbins + 1; ++i)
1103 if (group[d] * i <= newaxis->
GetNbins())
1104 edges[i] = newaxis->
GetXbins()->
At(group[d] * i);
1105 else edges[i] = newaxis->
GetXmax();
1106 newaxis->
Set(newbins, edges);
1116 Bool_t wantErrors = haveErrors;
1123 while ((i = iter.
Next(coord)) >= 0) {
1125 for (
Int_t d = 0; d < ndim; ++d) {
1175 Error(
"ComputeIntegral",
"The histogram must have at least one bin.");
1187 while ((i = iter.
Next(coord)) >= 0) {
1191 bool regularBin =
true;
1200 if (!regularBin) v = 0.;
1208 Error(
"ComputeIntegral",
"No hits in regular bins (non over/underflow).");
1257 if (v == 0. && err == 0. && options && strchr(options,
'0')) {
1270 Printf(
"Bin at (%s) = %g (+/- %g)", coord.
Data(),
v, err);
1287 if (from < 0) from = 0;
1288 if (howmany == -1) howmany =
GetNbins();
1292 if (options && (strchr(options,
'x') || strchr(options,
'X'))) {
1296 bin[dim] = from % nbins[dim];
1300 for (
Long64_t i = 0; i < howmany; ++i) {
1306 if (bin[dim] >= nbins[dim]) {
1318 for (
Long64_t i = from; i < from + howmany; ++i) {
1336 Bool_t optAxis = options && (strchr(options,
'A') || (strchr(options,
'a')));
1337 Bool_t optMem = options && (strchr(options,
'M') || (strchr(options,
'm')));
1338 Bool_t optStat = options && (strchr(options,
'S') || (strchr(options,
's')));
1339 Bool_t optContent = options && (strchr(options,
'C') || (strchr(options,
'c')));
1347 Printf(
" axis %d \"%s\": %d bins (%g..%g), %s bin sizes", dim,
1349 (axis->
GetXbins() ?
"variable" :
"fixed"));
1365 Printf(
" coordinates stored in %d chunks of %d entries\n %g of bins filled using %g of memory compared to an array",
1433 fHist(hist), fAxis(axis), fProj(0)
Abstract array base class.
void Add(TObject *obj, const char *name=0, Int_t check=-1)
Add object with name to browser.
Double_t ComputeIntegral()
Calculate the integral of the histogram.
Double_t GetBinError(const Int_t *idx) const
virtual const char * GetName() const
Returns name of object.
Double_t GetSparseFractionMem() const
Return the amount of used memory over memory that would be used by a non-sparse n-dimensional histogr...
virtual Int_t GetNcells() const
virtual Double_t GetEffectiveEntries() const
Number of effective entries of the histogram.
virtual void Reset(Option_t *option="")=0
Bool_t GetCalculateErrors() const
void Divide(const THnBase *h)
Divide this histogram by h this = this/(h) Note that if h has Sumw2 set, Sumw2 is automatically calle...
Double_t Floor(Double_t x)
Int_t GetFirst() const
Return first bin on the axis i.e.
void SetEntries(Double_t entries)
TFitResultPtr Fit(TF1 *f1, Option_t *option="", Option_t *goption="")
Fit a THnSparse with function f.
THnBase * RebinBase(Int_t group) const
Combine the content of "group" neighboring bins into a new bin and return the resulting THnBase...
virtual Long64_t GetNbins() const =0
void Multiply(const THnBase *h)
Multiply this histogram by histogram h this = this * h Note that if h has Sumw2 set, Sumw2 is automatically called for this if not already set.
virtual Double_t GetBinLowEdge(Int_t bin) const
Return low edge of bin.
TString & ReplaceAll(const TString &s1, const TString &s2)
const Double_t * GetArray() const
Bool_t TestBit(UInt_t f) const
virtual void SetName(const char *name)
Set the name of the TNamed.
virtual void SetOwner(Bool_t enable=kTRUE)
Set whether this collection is the owner (enable==true) of its content.
virtual Double_t GetBinContent(Int_t bin) const
Return content of bin number bin.
Double_t GetBinContent(const Int_t *idx) const
Ssiz_t Index(const char *pat, Ssiz_t i=0, ECaseCompare cmp=kExact) const
void Scale(Double_t c)
Scale contents and errors of this histogram by c: this = this * c It does not modify the histogram's ...
virtual void GetBinXYZ(Int_t binglobal, Int_t &binx, Int_t &biny, Int_t &binz) const
Return binx, biny, binz corresponding to the global bin number globalbin see TH1::GetBin function abo...
static THnBase * CreateHnAny(const char *name, const char *title, const TH1 *h1, Bool_t sparse, Int_t chunkSize=1024 *16)
Create a THn / THnSparse object from a histogram deriving from TH1.
void GetRandom(Double_t *rand, Bool_t subBinRandom=kTRUE)
Generate an n-dimensional random tuple based on the histogrammed distribution.
void ToLower()
Change string to lower-case.
TObject * At(Int_t idx) const
virtual void SetNameTitle(const char *name, const char *title)
Set all the TNamed parameters (name and title).
TString & Insert(Ssiz_t pos, const char *s)
void AddBinContent(const Int_t *x, Double_t v=1.)
void SetBit(UInt_t f, Bool_t set)
Set or unset the user status bits as specified in f.
virtual Double_t GetBinUpEdge(Int_t bin) const
Return up edge of bin.
virtual Int_t GetDimension() const
Double_t fEntries
browser-helpers for each axis
THnBaseBrowsable(THnBase *hist, Int_t axis)
Construct a THnBaseBrowsable.
virtual const char * ClassName() const
Returns name of class to which the object belongs.
TBrowser helper for THnBase.
static Int_t FitOptionsMake(Option_t *option, Foption_t &Foption)
Decode string choptin and fill fitOption structure.
Efficient multidimensional histogram.
static TString Format(const char *fmt,...)
Static method which formats a string using a printf style format descriptor and return a TString...
The TNamed class is the base class for all named ROOT classes.
TH1D * Projection(Int_t xDim, Option_t *option="") const
THnBase * CloneEmpty(const char *name, const char *title, const TObjArray *axes, Bool_t keepTargetAxis) const
Create a new THnBase object that is of the same type as *this, but with dimensions and bins given by ...
Int_t GetCoord(Int_t dim) const
virtual void Reserve(Long64_t)
virtual Double_t GetBinCenter(Int_t bin) const
Return center of bin.
void PrintEntries(Long64_t from=0, Long64_t howmany=-1, Option_t *options=0) const
Print "howmany" entries starting at "from".
void SetBinContent(const Int_t *idx, Double_t v)
std::vector< std::vector< double > > Data
Double_t GetSumw2() const
virtual ~THnBaseBinIter()
Destruct a bin iterator.
Ssiz_t First(char c) const
Find first occurrence of a character c.
Bool_t CheckConsistency(const THnBase *h, const char *tag) const
Consistency check on (some of) the parameters of two histograms (for operations). ...
void Browse(TBrowser *b)
Browse an axis of a THnBase, i.e. draw its projection.
virtual void AddBinContent(Int_t bin)
Increment bin content by 1.
virtual void ResetStats()
Reset the statistics including the number of entries and replace with values calculates from bin cont...
virtual void SetBinError(Int_t bin, Double_t error)
See convention for numbering bins in TH1::GetBin.
virtual Double_t GetBinError2(Long64_t linidx) const =0
Bool_t PrintBin(Long64_t idx, Int_t *coord, Option_t *options) const
Print one bin.
virtual Double_t Rndm()
Machine independent random number generator.
THist< 3, double, THistStatContent, THistStatUncertainty > TH3D
Using a TBrowser one can browse all ROOT objects.
RooArgSet S(const RooAbsArg &v1)
void AddInternal(const THnBase *h, Double_t c, Bool_t rebinned)
Add() implementation for both rebinned histograms and those with identical binning.
virtual void SetRange(Int_t first=0, Int_t last=0)
Set the viewing range for the axis from bin first to last.
virtual void InitStorage(Int_t *nbins, Int_t chunkSize)=0
virtual void SetFilledBins(Long64_t)
TAxis * GetAxis(Int_t dim) const
virtual void AddAtAndExpand(TObject *obj, Int_t idx)
Add object at position idx.
void Print(Option_t *option="") const
Print a THnBase.
const char * GetTitle() const
Returns title of object.
Int_t GetLast() const
Return last bin on the axis i.e.
Class to manage histogram axis.
virtual void Draw(Option_t *option="")
Draw this histogram with options.
Long64_t Next(Int_t *coord=0)
virtual Long64_t GetBin(const Int_t *idx) const =0
void Init(const char *name, const char *title, const TObjArray *axes, Bool_t keepTargetAxis, Int_t chunkSize=1024 *16)
Initialize axes and name.
Provides an indirection to the TFitResult class and with a semantics identical to a TFitResult pointe...
virtual Bool_t InheritsFrom(const char *classname) const
Returns kTRUE if object inherits from class "classname".
Collection abstract base class.
virtual ~THnBase()
Destruct a THnBase.
Int_t GetEntriesFast() const
Double_t GetSumwx2(Int_t dim) const
virtual void Error(const char *method, const char *msgfmt,...) const
Issue error message.
char * Form(const char *fmt,...)
static void RejectPoint(Bool_t reject=kTRUE)
Static function to set the global flag to reject points the fgRejectPoint global flag is tested by al...
The ROOT global object gROOT contains a list of all defined classes.
void SetBinEdges(Int_t idim, const Double_t *bins)
Set the axis # of bins and bin limits on dimension idim.
R__EXTERN TRandom * gRandom
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 void Rebuild(Option_t *option="")
Using the current bin info, recompute the arrays for contents and errors.
virtual void AddAt(TObject *obj, Int_t idx)
Add object at position ids.
Double_t GetEntries() const
Double_t At(Int_t i) const
Double_t GetSumwx(Int_t dim) const
TString & Remove(Ssiz_t pos)
void AddRange(unsigned int icoord, double xmin, double xmax)
add a range [xmin,xmax] for the new coordinate icoord Adding a range does not delete existing one...
void SetBinError(const Int_t *idx, Double_t e)
virtual Bool_t IsEmpty() const
virtual void SetBinError2(Long64_t bin, Double_t e2)=0
class describing the range in the coordinates it supports multiple range in a coordinate.
Bool_t IsInRange(Int_t *coord) const
Check whether bin coord is in range, as defined by TAxis::SetRange().
void Add(const THnBase *h, Double_t c=1.)
Add contents of h scaled by c to this histogram: this = this + c * h Note that if h has Sumw2 set...
Long64_t Merge(TCollection *list)
Merge this with a list of THnBase's.
Bool_t Contains(const char *pat, ECaseCompare cmp=kExact) const
~THnBaseBrowsable()
Destruct a THnBaseBrowsable.
enum THnBase::@97 fIntegralStatus
array with bin weight sums
THist< 2, double, THistStatContent, THistStatUncertainty > TH2D
virtual void AddBinError2(Long64_t bin, Double_t e2)=0
Double_t GetSparseFractionBins() const
Return the amount of filled bins over all bins.
static TClass * GetClass(const char *name, Bool_t load=kTRUE, Bool_t silent=kFALSE)
Static method returning pointer to TClass of the specified class name.
TObject * ProjectionAny(Int_t ndim, const Int_t *dim, Bool_t wantNDim, Option_t *option="") const
Project all bins into a ndim-dimensional THn / THnSparse (whatever *this is) or if (ndim < 4 and !wan...
Int_t GetNdimensions() const
void Browse(TBrowser *b)
Browse a THnSparse: create an entry (ROOT::THnSparseBrowsable) for each dimension.
Mother of all ROOT objects.
virtual Bool_t IsInside(const Double_t *x) const
return kTRUE if the point is inside the function range
virtual Double_t GetBinWidth(Int_t bin) const
Return bin width.
Bool_t HaveSkippedBin() const
Iterator over THnBase bins.
you should not use this method at all Int_t Int_t Double_t Double_t Double_t Int_t Double_t Double_t Double_t Double_t b
THist< 1, double, THistStatContent, THistStatUncertainty > TH1D
TObjArray * GetListOfAxes()
virtual void SetEntries(Double_t n)
Multidimensional histogram base.
TFitResultPtr FitObject(TH1 *h1, TF1 *f1, Foption_t &option, const ROOT::Math::MinimizerOptions &moption, const char *goption, ROOT::Fit::DataRange &range)
fitting function for a TH1 (called from TH1::Fit)
Option_t * GetDrawOption() const
Get option used by the graphics system to draw this object.
Double_t Sqrt(Double_t x)
TH1 * CreateHist(const char *name, const char *title, const TObjArray *axes, Bool_t keepTargetAxis) const
Create an empty histogram with name and title with a given set of axes.
virtual const char * GetName() const
Returns name of object.
virtual Int_t GetSize() const
virtual void Set(Int_t nbins, Double_t xmin, Double_t xmax)
Initialize axis with fix bins.
Int_t GetChunkSize() const
virtual Double_t EvalPar(const Double_t *x, const Double_t *params=0)
Evaluate function with given coordinates and parameters.
virtual void SetTitle(const char *title="")
Set the title of the TNamed.
void RebinnedAdd(const THnBase *h, Double_t c=1.)
Add contents of h scaled by c to this histogram: this = this + c * h Note that if h has Sumw2 set...
Int_t CeilNint(Double_t x)
Long64_t BinarySearch(Long64_t n, const T *array, T value)
const TArrayD * GetXbins() const
void ResetBase(Option_t *option="")
Clear the histogram.
virtual void Warning(const char *method, const char *msgfmt,...) const
Issue warning message.
void SetTitle(const char *title)
Change (i.e.
virtual const char * GetTitle() const
Returns title of object.
virtual Double_t GetBinError(Int_t bin) const
Return value of error associated to bin number bin.
void * New(ENewType defConstructor=kClassNew, Bool_t quiet=kFALSE) const
Return a pointer to a newly allocated object of this class.
const char * Data() const