Logo ROOT   6.30.04
Reference Guide
 All Namespaces Files Pages
TProfile2D.h
Go to the documentation of this file.
1 // @(#)root/hist:$Id$
2 // Author: Rene Brun 16/04/2000
3 
4 /*************************************************************************
5  * Copyright (C) 1995-2000, Rene Brun and Fons Rademakers. *
6  * All rights reserved. *
7  * *
8  * For the licensing terms see $ROOTSYS/LICENSE. *
9  * For the list of contributors see $ROOTSYS/README/CREDITS. *
10  *************************************************************************/
11 
12 #ifndef ROOT_TProfile2D
13 #define ROOT_TProfile2D
14 
15 
16 //////////////////////////////////////////////////////////////////////////
17 // //
18 // TProfile2D //
19 // //
20 // Profile2D histogram class. //
21 // //
22 //////////////////////////////////////////////////////////////////////////
23 
24 #include "TH2.h"
25 #include "TProfile.h"
26 
27 class TProfile2D : public TH2D {
28 
29 public:
30  friend class TProfileHelper;
31 
32 protected:
33  TArrayD fBinEntries; //number of entries per bin
34  EErrorType fErrorMode; //Option to compute errors
35  Double_t fZmin; //Lower limit in Z (if set)
36  Double_t fZmax; //Upper limit in Z (if set)
37  Bool_t fScaling; //!True when TProfile2D::Scale is called
38  Double_t fTsumwz; //Total Sum of weight*Z
39  Double_t fTsumwz2; //Total Sum of weight*Z*Z
40  TArrayD fBinSumw2; //Array of sum of squares of weights per bin
41  static Bool_t fgApproximate; //bin error approximation option
42 
43  virtual Int_t BufferFill(Double_t, Double_t) {return -2;} //may not use
44  virtual Int_t BufferFill(Double_t, Double_t, Double_t) {return -2;} //may not use
45  virtual Int_t BufferFill(Double_t x, Double_t y, Double_t z, Double_t w);
46 
47  // helper methods for the Merge unification in TProfileHelper
48  void SetBins(const Int_t* nbins, const Double_t* range) { SetBins(nbins[0], range[0], range[1],
49  nbins[1], range[2], range[3]); };
50  Int_t Fill(const Double_t* v) { return Fill(v[0], v[1], v[2], v[3]); };
51 
52  virtual TProfile *DoProfile(bool onX, const char *name, Int_t firstbin, Int_t lastbin, Option_t *option) const;
53 
54  using TH2::Fill;
55  Int_t Fill(Double_t, Double_t) {return TH2::Fill(0); } //MayNotUse
56 
57  virtual Double_t RetrieveBinContent(Int_t bin) const { return (fBinEntries.fArray[bin] > 0) ? fArray[bin]/fBinEntries.fArray[bin] : 0; }
58  //virtual void UpdateBinContent(Int_t bin, Double_t content);
59  virtual Double_t GetBinErrorSqUnchecked(Int_t bin) const { Double_t err = GetBinError(bin); return err*err; }
60 
61 private:
62  Double_t *GetB() {return &fBinEntries.fArray[0];}
63  Double_t *GetB2() {return (fBinSumw2.fN ? &fBinSumw2.fArray[0] : 0 ); }
64  Double_t *GetW() {return &fArray[0];}
65  Double_t *GetW2() {return &fSumw2.fArray[0];}
66  void SetBins(Int_t, Double_t, Double_t)
67  { MayNotUse("SetBins(Int_t, Double_t, Double_t"); }
68  void SetBins(Int_t, const Double_t*)
69  { MayNotUse("SetBins(Int_t, const Double_t*"); }
70  void SetBins(Int_t, Double_t, Double_t, Int_t, Double_t, Double_t, Int_t, Double_t, Double_t)
71  { MayNotUse("SetBins(Int_t, Double_t, Double_t, Int_t, Double_t, Double_t, Int_t, Double_t, Double_t"); }
72  void SetBins(Int_t, const Double_t *, Int_t, const Double_t *, Int_t, const Double_t *)
73  { MayNotUse("SetBins(Int_t, const Double_t*, Int_t, const Double_t*, Int_t, const Double_t*"); }
74 
75 public:
76  TProfile2D();
77  TProfile2D(const char *name,const char *title,Int_t nbinsx,Double_t xlow,Double_t xup
78  ,Int_t nbinsy,Double_t ylow,Double_t yup
79  ,Double_t zlow, Double_t zup,Option_t *option="");
80  TProfile2D(const char *name,const char *title,Int_t nbinsx,Double_t xlow,Double_t xup
81  ,Int_t nbinsy,Double_t ylow,Double_t yup,Option_t *option="");
82  TProfile2D(const char *name,const char *title,Int_t nbinsx,const Double_t *xbins
83  ,Int_t nbinsy,Double_t ylow,Double_t yup,Option_t *option="");
84  TProfile2D(const char *name,const char *title,Int_t nbinsx,Double_t xlow,Double_t xup
85  ,Int_t nbinsy,const Double_t *ybins,Option_t *option="");
86  TProfile2D(const char *name,const char *title,Int_t nbinsx,const Double_t *xbins
87  ,Int_t nbinsy,const Double_t *ybins,Option_t *option="");
88  TProfile2D(const TProfile2D &profile);
89  TProfile2D &operator=(const TProfile2D &profile);
90  virtual ~TProfile2D();
91  virtual Bool_t Add(TF1 *h1, Double_t c1=1, Option_t *option="");
92  virtual Bool_t Add(const TH1 *h1, Double_t c1=1);
93  virtual Bool_t Add(const TH1 *h1, const TH1 *h2, Double_t c1=1, Double_t c2=1); // *MENU*
94  static void Approximate(Bool_t approx=kTRUE);
95  void BuildOptions(Double_t zmin, Double_t zmax, Option_t *option);
96  virtual Int_t BufferEmpty(Int_t action=0);
97  virtual void Copy(TObject &hnew) const;
98  virtual Bool_t Divide(TF1 *h1, Double_t c1=1);
99  virtual Bool_t Divide(const TH1 *h1);
100  virtual Bool_t Divide(const TH1 *h1, const TH1 *h2, Double_t c1=1, Double_t c2=1, Option_t *option=""); // *MENU*
101  virtual void ExtendAxis(Double_t x, TAxis *axis);
102  Int_t Fill(Double_t x, Double_t y, Double_t z);
103  virtual Int_t Fill(Double_t x, const char *namey, Double_t z);
104  virtual Int_t Fill(const char *namex, Double_t y, Double_t z);
105  virtual Int_t Fill(const char *namex, const char *namey, Double_t z);
106  virtual Int_t Fill(Double_t x, Double_t y, Double_t z, Double_t w);
107  virtual Double_t GetBinContent(Int_t bin) const;
108  virtual Double_t GetBinContent(Int_t binx, Int_t biny) const {return GetBinContent(GetBin(binx,biny));}
109  virtual Double_t GetBinContent(Int_t binx, Int_t biny, Int_t) const {return GetBinContent(GetBin(binx,biny));}
110  virtual Double_t GetBinError(Int_t bin) const;
111  virtual Double_t GetBinError(Int_t binx, Int_t biny) const {return GetBinError(GetBin(binx,biny));}
112  virtual Double_t GetBinError(Int_t binx, Int_t biny, Int_t) const {return GetBinError(GetBin(binx,biny));}
113  virtual Double_t GetBinEntries(Int_t bin) const;
114  virtual Double_t GetBinEffectiveEntries(Int_t bin);
115  virtual TArrayD *GetBinSumw2() {return &fBinSumw2;}
116  virtual const TArrayD *GetBinSumw2() const {return &fBinSumw2;}
117  Option_t *GetErrorOption() const;
118  virtual void GetStats(Double_t *stats) const;
119  virtual Double_t GetZmin() const {return fZmin;}
120  virtual Double_t GetZmax() const {return fZmax;}
121  virtual void LabelsDeflate(Option_t *axis="X");
122  virtual void LabelsInflate(Option_t *axis="X");
123  virtual void LabelsOption(Option_t *option="h", Option_t *axis="X");
124  virtual Long64_t Merge(TCollection *list);
125  virtual Bool_t Multiply(TF1 *h1, Double_t c1=1);
126  virtual Bool_t Multiply(const TH1 *h1);
127  virtual Bool_t Multiply(const TH1 *h1, const TH1 *h2, Double_t c1=1, Double_t c2=1, Option_t *option=""); // *MENU*
128  TH2D *ProjectionXY(const char *name="_pxy", Option_t *option="e") const;
129  TProfile *ProfileX(const char *name="_pfx", Int_t firstybin=0, Int_t lastybin=-1, Option_t *option="") const; // *MENU*
130  TProfile *ProfileY(const char *name="_pfy", Int_t firstxbin=0, Int_t lastxbin=-1, Option_t *option="") const; // *MENU*
131  virtual void PutStats(Double_t *stats);
132  virtual void Reset(Option_t *option="");
133  virtual TProfile2D *Rebin2D(Int_t nxgroup=2, Int_t nygroup=2, const char *newname="");
134  virtual TProfile2D *RebinX(Int_t ngroup=2, const char *newname="");
135  virtual TProfile2D *RebinY(Int_t ngroup=2, const char *newname="");
136  virtual void SavePrimitive(std::ostream &out, Option_t *option = "");
137  virtual void Scale(Double_t c1=1, Option_t *option="");
138  virtual void SetBinEntries(Int_t bin, Double_t w);
139  virtual void SetBins(Int_t nbinsx, Double_t xmin, Double_t xmax, Int_t nbinsy, Double_t ymin, Double_t ymax);
140  virtual void SetBins(Int_t nx, const Double_t *xBins, Int_t ny, const Double_t *yBins);
141  virtual void SetBinsLength(Int_t n=-1);
142  virtual void SetBuffer(Int_t buffersize, Option_t *option="");
143  virtual void SetErrorOption(Option_t *option=""); // *MENU*
144  virtual void Sumw2(Bool_t flag = kTRUE);
145  Double_t GetNumberOfBins() { return fBinEntries.GetSize(); }
146 
147  ClassDef(TProfile2D,8) //Profile2D histogram class
148 };
149 
150 #endif