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

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