source: trunk/MagicSoft/Mars/mhist/MHCamera.h@ 3378

Last change on this file since 3378 was 3237, checked in by tbretz, 21 years ago
*** empty log message ***
File size: 16.5 KB
Line 
1#ifndef MARS_MHCamera
2#define MARS_MHCamera
3
4#ifndef MARS_MAGIC
5#include "MAGIC.h"
6#endif
7#ifndef ROOT_TArrayC
8#include <TArrayC.h>
9#endif
10#ifndef ROOT_TArrayI
11#include <TArrayI.h>
12#endif
13#ifndef ROOT_TClonesArray
14#include <TClonesArray.h>
15#endif
16#ifndef ROOT_TH1
17#include <TH1.h>
18#endif
19
20class TPaveStats;
21
22class MGeomCam;
23class MCamEvent;
24class MRflEvtData;
25class MCerPhotEvt;
26class MImgCleanStd;
27
28class MHCamera : public TH1D
29{
30public:
31 enum {
32 kProfile = BIT(18), // FIXME: When changing change max/min!
33 kFreezed = BIT(19),
34 kNoLegend = BIT(20),
35 kVariance = BIT(21)
36 };
37private:
38 MGeomCam *fGeomCam; // pointer to camera geometry (y-axis)
39 TArrayC fUsed; // array containing flags
40
41 TArrayI fColors; //! Color conversion table
42 TList *fNotify; //!
43
44 Bool_t fFreezed; //! Just a dummy!!!! ([Set,Is]Freezed)
45
46 //Int_t fOptStat;
47 //TGStatusBar *fStatusBar;
48
49 void Init();
50
51 Stat_t Profile(Stat_t val) const
52 {
53 if (!TestBit(kProfile))
54 return val;
55
56 const Stat_t n = TH1D::GetEntries();
57 return n>0 ? val/n : val;
58 }
59
60 Int_t GetColor(Float_t val, Float_t min, Float_t max, Bool_t islog);
61
62 void PaintIndices(Int_t type);
63 void Update(Bool_t islog, Bool_t isbox, Bool_t iscol);
64 void UpdateLegend(Float_t min, Float_t max, Bool_t islog);
65 void SetRange();
66
67 TPaveStats *GetStatisticBox();
68
69 Int_t GetPixelIndex(Int_t px, Int_t py) const;
70
71 void PaintAxisTitle();
72
73 enum {
74 kIsUsed = BIT(1)
75 };
76
77 void SetUsed(Int_t idx) { SETBIT(fUsed[idx], kIsUsed); }
78 void ResetUsed(Int_t idx) { CLRBIT(fUsed[idx], kIsUsed); }
79
80 Bool_t MatchSector(Int_t idx, Int_t sector) const;
81
82 // This is a trick to remove TH1 entries from the context menu
83 TH1 *Rebin(Int_t ngroup=2, const char*newname="") { return this; }
84 void DrawPanel() {}
85
86 Int_t Fill(Axis_t x);
87 Int_t Fill(Axis_t x, Stat_t w);
88 Int_t Fill(const char *name, Stat_t w) { return -1; }
89 void FillN(Int_t ntimes, const Axis_t *x, const Double_t *w, Int_t stride=1) {}
90 void FillN(Int_t, const Axis_t *, const Axis_t *, const Double_t *, Int_t) {}
91
92public:
93 MHCamera();
94 MHCamera(const MGeomCam &geom, const char *name="", const char *title="");
95 ~MHCamera();
96
97 void SetGeometry(const MGeomCam &geom, const char *name="", const char *title="");
98 const MGeomCam* GetGeometry() const { return fGeomCam; }
99
100 Bool_t IsUsed(Int_t idx) const { return TESTBIT(const_cast<TArrayC&>(fUsed)[idx], kIsUsed); }
101
102 Int_t Fill(Axis_t x, Axis_t y, Stat_t w);
103
104 //void AddPixContent(Int_t idx) const { AddBinContent(idx+1); }
105 //void AddPixContent(Int_t idx, Stat_t w) const { AddBinContent(idx+1, w); }
106
107 // This is a trick to remove TH1 entries from the context menu
108 /*
109 void Add(const TH1 *h1, const TH1 *h2, Double_t c1=1, Double_t c2=1) { TH1::Add(h1, h2, c1, c2); }
110 void Add(TF1 *h1, Double_t c1=1) { TH1::Add(h1, c1); }
111 void Add(const TH1 *h1, Double_t c1=1) { TH1::Add(h1, c1); }
112 void Divide(TF1 *f1, Double_t c1=1) { TH1::Divide(f1, c1); }
113 void Divide(const TH1 *h1) { TH1::Divide(h1); }
114 void Divide(const TH1 *h1, const TH1 *h2, Double_t c1=1, Double_t c2=1, Option_t *option="") { TH1::Divide(h1, h2, c1, c2, option); }
115 void Multiply(TF1 *h1, Double_t c1=1) { TH1::Multiply(h1, c1); }
116 void Multiply(const TH1 *h1) { TH1::Multiply(h1); }
117 void Multiply(const TH1 *h1, const TH1 *h2, Double_t c1=1, Double_t c2=1, Option_t *option="") { TH1::Multiply(h1, h2, c1, c2, option); }
118 */
119
120 void FitPanel() { TH1::FitPanel(); }
121
122 virtual Double_t GetPixContent(Int_t idx) const { return GetBinContent(idx+1); }
123 virtual void AddCamContent(const MCamEvent &evt, Int_t type=0);
124 virtual void AddCamContent(const MHCamera &evt, Int_t type=0);
125 virtual void AddCamContent(const TArrayD &arr, const TArrayC *used=NULL);
126 virtual void SetCamContent(const MCamEvent &evt, Int_t type=0) { Reset(); AddCamContent(evt, type); }
127 virtual void SetCamContent(const TArrayD &evt, const TArrayC *used=NULL) { Reset(); AddCamContent(evt, used); }
128 virtual void SetCamContent(const MHCamera &d, Int_t type=0) { Reset(); AddCamContent(d, type); fEntries=d.fEntries; }
129
130 virtual void SetCamError(const MCamEvent &evt, Int_t type=0);
131
132 virtual void CntCamContent(const MCamEvent &evt, Double_t threshold, Int_t type=0);
133 virtual void CntCamContent(const TArrayD &evt, Double_t threshold, Bool_t ispos=kTRUE);
134
135 Stat_t GetBinContent(Int_t bin) const { return Profile(TH1D::GetBinContent(bin)); }
136 Stat_t GetBinContent(Int_t binx, Int_t biny) const { return GetBinContent(binx); }
137 Stat_t GetBinContent(Int_t binx, Int_t biny, Int_t binz) const { return GetBinContent(binx); }
138 Stat_t GetBinError(Int_t bin) const;
139 Stat_t GetBinError(Int_t binx, Int_t biny) const { return GetBinError(binx); }
140 Stat_t GetBinError(Int_t binx, Int_t biny, Int_t binz) const { return GetBinError(binx); }
141
142 Double_t GetMinimum(Bool_t all) const { return GetMinimumSector(-1, all); }
143 Double_t GetMaximum(Bool_t all) const { return GetMaximumSector(-1, all); }
144
145 Double_t GetMinimum() const { return GetMinimumSector(-1, kFALSE); }
146 Double_t GetMaximum() const { return GetMaximumSector(-1, kFALSE); }
147
148 Double_t GetMinimumSector(Int_t sector, Bool_t all=kFALSE) const;
149 Double_t GetMaximumSector(Int_t sector, Bool_t all=kFALSE) const;
150
151 void SetLevels(const TArrayF &arr);
152
153 void FillRandom(const char *fname, Int_t ntimes=5000) { TH1::FillRandom(fname, ntimes); }
154 void FillRandom(TH1 *h, Int_t ntimes=5000) { TH1::FillRandom(h, ntimes); }
155 void FillRandom();
156
157 void PrintInfo() const { Print(""); } // *MENU*
158 void Reset(Option_t *);
159 void Reset() { Reset(""); } // *MENU*
160 TH1 *DrawCopy() const/* { gPad=NULL; return TH1D::DrawCopy(); }*/; // *MENU*
161 TH1 *DrawCopy(Option_t *o) const { return TH1D::DrawCopy(o); }
162
163 void Print(Option_t *) const;
164 void Paint(Option_t *option="");
165 void Draw(Option_t *option="");
166 TObject *DrawClone(Option_t *option="") const;
167 void SavePrimitive(ofstream &out, Option_t *);
168 Int_t DistancetoPrimitive(Int_t px, Int_t py);
169 char *GetObjectInfo(Int_t px, Int_t py) const;
170 void ExecuteEvent(Int_t event, Int_t px, Int_t py);
171
172 void SetPalette(Int_t ncolors, Int_t *colors);
173
174 void SetPrettyPalette(); // *MENU*
175 void SetDeepBlueSeaPalette(); // *MENU*
176 void SetInvDeepBlueSeaPalette(); // *MENU*
177
178 void SetAutoScale() { fMinimum = fMaximum = -1111; } // *MENU*
179 void DisplayAsHistogram() { SetDrawOption("histEP"); } // *MENU*
180 void DisplayAsCamera() { SetDrawOption(""); } // *MENU*
181
182 void SetFreezed(Bool_t f=kTRUE) { f ? SetBit(kFreezed) : ResetBit(kFreezed); } // *TOGGLE*
183 Bool_t IsFreezed() const { return TestBit(kFreezed); }
184 //void SetOptStat(Int_t os=-1) { fOptStat = os; } // *MENU*
185
186 void AddNotify(TObject *event);
187
188 Stat_t GetMean(Bool_t all) const { return GetMeanSector(-1, all); }
189 Stat_t GetRMS(Bool_t all) const { return GetRmsSector(-1, all); }
190
191 Stat_t GetMean(Int_t sector=-1) const { return GetMeanSector(-1, kFALSE); }
192 Stat_t GetRMS(Int_t sector=-1) const { return GetRmsSector(-1, kFALSE); }
193
194 Stat_t GetMeanSector(Int_t sector, Bool_t all=kFALSE) const;
195 Stat_t GetRmsSector(Int_t sector, Bool_t all=kFALSE) const;
196
197 UInt_t GetNumPixels() const;
198
199 TH1D *Projection(const char *name="_py") const { return ProjectionS(-1, name); }
200 TH1D *ProjectionS(Int_t sector, const char *name="_py") const;
201
202 const MGeomCam &GetGeomCam() const { return *fGeomCam; }
203
204 ClassDef(MHCamera, 1) // Displays the magic camera
205};
206
207#endif
208
209/* ------------ OK ---------------
210 virtual void Browse(TBrowser *b);
211 virtual void FillRandom(const char *fname, Int_t ntimes=5000);
212 virtual void FillRandom(TH1 *h, Int_t ntimes=5000);
213
214 virtual Int_t GetQuantiles(Int_t nprobSum, Double_t *q, const Double_t *probSum=0);
215 virtual Axis_t GetRandom();
216 virtual void GetStats(Stat_t *stats) const;
217 virtual Stat_t GetSumOfWeights() const;
218 virtual Int_t GetSumw2N() const {return fSumw2.fN;}
219 virtual Stat_t GetRMS(Int_t axis=1) const;
220
221 virtual Int_t GetNbinsX() const {return fXaxis.GetNbins();}
222 virtual Int_t GetNbinsY() const {return fYaxis.GetNbins();}
223 virtual Int_t GetNbinsZ() const {return fZaxis.GetNbins();}
224
225 // ------------- to check -------------------
226
227 virtual Double_t ComputeIntegral();
228 virtual Int_t DistancetoPrimitive(Int_t px, Int_t py);
229 virtual void Draw(Option_t *option="");
230 virtual TH1 *DrawCopy(Option_t *option="") const;
231 virtual TH1 *DrawNormalized(Option_t *option="", Double_t norm=1) const;
232 virtual Int_t BufferEmpty(Bool_t deleteBuffer=kFALSE);
233 virtual void Eval(TF1 *f1, Option_t *option="");
234 virtual void ExecuteEvent(Int_t event, Int_t px, Int_t py);
235 virtual void FillN(Int_t ntimes, const Axis_t *x, const Double_t *w, Int_t stride=1);
236 virtual void FillN(Int_t, const Axis_t *, const Axis_t *, const Double_t *, Int_t) {;}
237 virtual Int_t FindBin(Axis_t x, Axis_t y=0, Axis_t z=0);
238 virtual TObject *FindObject(const char *name) const;
239 virtual TObject *FindObject(const TObject *obj) const;
240 virtual Int_t Fit(const char *formula ,Option_t *option="" ,Option_t *goption="", Axis_t xmin=0, Axis_t xmax=0); // *MENU*
241 virtual Int_t Fit(TF1 *f1 ,Option_t *option="" ,Option_t *goption="", Axis_t xmin=0, Axis_t xmax=0);
242 virtual void FitPanel(); // *MENU*
243 TH1 *GetAsymmetry(TH1* h2, Double_t c2=1, Double_t dc2=0);
244 Int_t GetBufferLength() const {return (Int_t)fBuffer[0];}
245 Int_t GetBufferSize () const {return fBufferSize;}
246 const Double_t *GetBuffer() const {return fBuffer;}
247 static Int_t GetDefaultBufferSize();
248 virtual Double_t *GetIntegral() {return fIntegral;}
249
250 TList *GetListOfFunctions() const { return fFunctions; }
251
252 virtual Int_t GetNdivisions(Option_t *axis="X") const;
253 virtual Color_t GetAxisColor(Option_t *axis="X") const;
254 virtual Color_t GetLabelColor(Option_t *axis="X") const;
255 virtual Style_t GetLabelFont(Option_t *axis="X") const;
256 virtual Float_t GetLabelOffset(Option_t *axis="X") const;
257 virtual Float_t GetLabelSize(Option_t *axis="X") const;
258 virtual Float_t GetTitleOffset(Option_t *axis="X") const;
259 virtual Float_t GetTitleSize(Option_t *axis="X") const;
260 virtual Float_t GetTickLength(Option_t *axis="X") const;
261 virtual Float_t GetBarOffset() const {return Float_t(0.001*Float_t(fBarOffset));}
262 virtual Float_t GetBarWidth() const {return Float_t(0.001*Float_t(fBarWidth));}
263 virtual Int_t GetContour(Double_t *levels=0);
264 virtual Double_t GetContourLevel(Int_t level) const;
265 virtual Double_t GetContourLevelPad(Int_t level) const;
266
267 virtual void GetCenter(Axis_t *center) const {fXaxis.GetCenter(center);}
268 TDirectory *GetDirectory() const {return fDirectory;}
269 virtual Stat_t GetEntries() const;
270 virtual TF1 *GetFunction(const char *name) const;
271 virtual Int_t GetDimension() const { return fDimension; }
272 virtual void GetLowEdge(Axis_t *edge) const {fXaxis.GetLowEdge(edge);}
273 virtual Double_t GetMaximum() const;
274 virtual Int_t GetMaximumBin() const;
275 virtual Int_t GetMaximumBin(Int_t &locmax, Int_t &locmay, Int_t &locmaz) const;
276 virtual Double_t GetMaximumStored() const {return fMaximum;}
277 virtual Double_t GetMinimum() const;
278 virtual Int_t GetMinimumBin() const;
279 virtual Int_t GetMinimumBin(Int_t &locmix, Int_t &locmiy, Int_t &locmiz) const;
280 virtual Double_t GetMinimumStored() const {return fMinimum;}
281 virtual Stat_t GetMean(Int_t axis=1) const;
282 virtual Double_t GetNormFactor() const {return fNormFactor;}
283 virtual char *GetObjectInfo(Int_t px, Int_t py) const;
284 Option_t *GetOption() const {return fOption.Data();}
285
286 TVirtualHistPainter *GetPainter();
287
288 TAxis *GetXaxis() const;
289 TAxis *GetYaxis() const;
290 TAxis *GetZaxis() const;
291 virtual Stat_t Integral(Option_t *option="") const;
292 virtual Stat_t Integral(Int_t binx1, Int_t binx2, Option_t *option="") const;
293 virtual Stat_t Integral(Int_t, Int_t, Int_t, Int_t, Option_t * ="") const {return 0;}
294 virtual Stat_t Integral(Int_t, Int_t, Int_t, Int_t, Int_t, Int_t, Option_t * ="" ) const {return 0;}
295 virtual Double_t KolmogorovTest(TH1 *h2, Option_t *option="") const;
296 virtual void LabelsDeflate(Option_t *axis="X");
297 virtual void LabelsInflate(Option_t *axis="X");
298 virtual void LabelsOption(Option_t *option="h", Option_t *axis="X");
299 virtual Int_t Merge(TCollection *list);
300 virtual void Multiply(TF1 *h1, Double_t c1=1);
301 virtual void Multiply(const TH1 *h1);
302 virtual void Multiply(const TH1 *h1, const TH1 *h2, Double_t c1=1, Double_t c2=1, Option_t *option=""); // *MENU*
303 virtual void Paint(Option_t *option="");
304 virtual void Print(Option_t *option="") const;
305 virtual void PutStats(Stat_t *stats);
306 virtual TH1 *Rebin(Int_t ngroup=2, const char*newname=""); // *MENU*
307 virtual void RebinAxis(Axis_t x, Option_t *axis="X");
308 virtual void Rebuild(Option_t *option="");
309 virtual void RecursiveRemove(TObject *obj);
310 virtual void Reset(Option_t *option="");
311 virtual void SavePrimitive(ofstream &out, Option_t *option);
312 virtual void Scale(Double_t c1=1);
313 virtual void SetAxisColor(Color_t color=1, Option_t *axis="X");
314 virtual void SetAxisRange(Axis_t xmin, Axis_t xmax, Option_t *axis="X");
315 virtual void SetBarOffset(Float_t offset=0.25) {fBarOffset = Short_t(1000*offset);}
316 virtual void SetBarWidth(Float_t width=0.5) {fBarWidth = Short_t(1000*width);}
317 virtual void SetBinContent(Int_t bin, Stat_t content);
318 virtual void SetBinContent(Int_t binx, Int_t biny, Stat_t content);
319 virtual void SetBinContent(Int_t binx, Int_t biny, Int_t binz, Stat_t content);
320 virtual void SetBinError(Int_t bin, Stat_t error);
321 virtual void SetBinError(Int_t binx, Int_t biny, Stat_t error);
322 virtual void SetBinError(Int_t binx, Int_t biny, Int_t binz, Stat_t error);
323 virtual void SetBins(Int_t nx, Axis_t xmin, Axis_t xmax);
324 virtual void SetBins(Int_t nx, Axis_t xmin, Axis_t xmax, Int_t ny, Axis_t ymin, Axis_t ymax);
325 virtual void SetBins(Int_t nx, Axis_t xmin, Axis_t xmax, Int_t ny, Axis_t ymin, Axis_t ymax,
326 Int_t nz, Axis_t zmin, Axis_t zmax);
327 virtual void SetBinsLength(Int_t = -1) { } //refefined in derived classes
328 virtual void SetBuffer(Int_t buffersize, Option_t *option="");
329 virtual void SetCellContent(Int_t binx, Int_t biny, Stat_t content);
330 virtual void SetCellError(Int_t binx, Int_t biny, Stat_t content);
331 virtual void SetContent(const Stat_t *content);
332 virtual void SetContour(Int_t nlevels, const Double_t *levels=0);
333 virtual void SetContourLevel(Int_t level, Double_t value);
334 static void SetDefaultBufferSize(Int_t buffersize=1000);
335 virtual void SetDirectory(TDirectory *dir);
336 virtual void SetEntries(Stat_t n) {fEntries = n;};
337 virtual void SetError(const Stat_t *error);
338 virtual void SetLabelColor(Color_t color=1, Option_t *axis="X");
339 virtual void SetLabelFont(Style_t font=62, Option_t *axis="X");
340 virtual void SetLabelOffset(Float_t offset=0.005, Option_t *axis="X");
341 virtual void SetLabelSize(Float_t size=0.02, Option_t *axis="X");
342
343 virtual void SetMaximum(Double_t maximum=-1111); // *MENU*
344 virtual void SetMinimum(Double_t minimum=-1111); // *MENU*
345 virtual void SetName(const char *name); // *MENU*
346 virtual void SetNameTitle(const char *name, const char *title);
347 virtual void SetNdivisions(Int_t n=510, Option_t *axis="X");
348 virtual void SetNormFactor(Double_t factor=1) {fNormFactor = factor;}
349 virtual void SetStats(Bool_t stats=kTRUE);
350 virtual void SetOption(Option_t *option=" ") {fOption = option;}
351 virtual void SetTickLength(Float_t length=0.02, Option_t *axis="X");
352 virtual void SetTitleOffset(Float_t offset=1, Option_t *axis="X");
353 virtual void SetTitleSize(Float_t size=0.02, Option_t *axis="X");
354 virtual void SetTitle(const char *title);
355 virtual void SetXTitle(const char *title) {fXaxis.SetTitle(title);}
356 virtual void SetYTitle(const char *title) {fYaxis.SetTitle(title);}
357 virtual void SetZTitle(const char *title) {fZaxis.SetTitle(title);}
358 virtual void Smooth(Int_t ntimes=1); // *MENU*
359 static void SmoothArray(Int_t NN, Double_t *XX, Int_t ntimes=1);
360 static Double_t SmoothMedian(Int_t n, Double_t *a);
361 virtual void Sumw2();
362 void UseCurrentStyle();
363
364 ClassDef(TH1,4) //1-Dim histogram base class
365};
366
367*/
Note: See TracBrowser for help on using the repository browser.