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

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