1 | #include "MCamDisplay.h"
|
---|
2 |
|
---|
3 | #include <math.h>
|
---|
4 | #include <fstream.h>
|
---|
5 |
|
---|
6 | #include <TClonesArray.h>
|
---|
7 | #include <TCanvas.h>
|
---|
8 | #include <TStyle.h>
|
---|
9 | #include <TBox.h>
|
---|
10 | #include <TText.h>
|
---|
11 |
|
---|
12 | #include "MHexagon.h"
|
---|
13 | #include "MGeomCam.h"
|
---|
14 |
|
---|
15 | #include "MCerPhotPix.h"
|
---|
16 | #include "MCerPhotEvt.h"
|
---|
17 |
|
---|
18 | #define kItemsLegend 25 // see SetPalette(1,0)
|
---|
19 |
|
---|
20 | ClassImp(MCamDisplay);
|
---|
21 |
|
---|
22 | // ------------------------------------------------------------------------
|
---|
23 | //
|
---|
24 | // default constructor
|
---|
25 | //
|
---|
26 | MCamDisplay::MCamDisplay(MGeomCam *geom)
|
---|
27 | : fAutoScale(kTRUE), fMinPhe(-2), fMaxPhe(50), fW(0), fH(0), fDrawingPad(NULL)
|
---|
28 | {
|
---|
29 | //
|
---|
30 | // create the hexagons of the display
|
---|
31 | //
|
---|
32 | fNumPixels = geom->GetNumPixels();
|
---|
33 | fRange = geom->GetMaxRadius();
|
---|
34 |
|
---|
35 | //
|
---|
36 | // Construct all hexagons. Use new-operator with placement
|
---|
37 | //
|
---|
38 |
|
---|
39 | // root 3.02
|
---|
40 | // * base/inc/TObject.h:
|
---|
41 | // register BIT(8) as kNoContextMenu. If an object has this bit set it will
|
---|
42 | // not get an automatic context menu when clicked with the right mouse button.
|
---|
43 |
|
---|
44 | fPixels = new TClonesArray("MHexagon", fNumPixels);
|
---|
45 | for (UInt_t i=0; i<fNumPixels; i++)
|
---|
46 | new ((*fPixels)[i]) MHexagon((*geom)[i]);
|
---|
47 |
|
---|
48 | //
|
---|
49 | // set the color palette for the TBox elements
|
---|
50 | //
|
---|
51 | gStyle->SetPalette(1, 0);
|
---|
52 |
|
---|
53 | //
|
---|
54 | // set up the Legend
|
---|
55 | //
|
---|
56 | fLegend = new TClonesArray("TBox", kItemsLegend);
|
---|
57 | fLegText = new TClonesArray("TText", kItemsLegend);
|
---|
58 |
|
---|
59 | for (Int_t i = 0; i<kItemsLegend; i++)
|
---|
60 | {
|
---|
61 | TBox *newbox = new ((*fLegend)[i]) TBox;
|
---|
62 | TText *newtxt = new ((*fLegText)[i]) TText;
|
---|
63 |
|
---|
64 | const Float_t lvl = 50. / kItemsLegend * i;
|
---|
65 |
|
---|
66 | newbox->SetFillColor(GetColor(lvl));
|
---|
67 |
|
---|
68 | newtxt->SetTextSize(0.025);
|
---|
69 | newtxt->SetTextAlign(12);
|
---|
70 | }
|
---|
71 | }
|
---|
72 |
|
---|
73 | // ------------------------------------------------------------------------
|
---|
74 | //
|
---|
75 | // Destructor. Deletes TClonesArrays for hexagons and legend elements.
|
---|
76 | //
|
---|
77 | MCamDisplay::~MCamDisplay()
|
---|
78 | {
|
---|
79 | delete fPixels;
|
---|
80 | delete fLegend;
|
---|
81 | delete fLegText;
|
---|
82 | }
|
---|
83 |
|
---|
84 | inline void MCamDisplay::SetPixColor(const MCerPhotPix &pix)
|
---|
85 | {
|
---|
86 | (*this)[pix.GetPixId()].SetFillColor(GetColor(pix.GetNumPhotons()));
|
---|
87 | }
|
---|
88 |
|
---|
89 | // ------------------------------------------------------------------------
|
---|
90 | //
|
---|
91 | // This is called at any time the canvas should get repainted.
|
---|
92 | // Here we maintain an aspect ratio of 5/4=1.15. This makes sure,
|
---|
93 | // that the camera image doesn't get distorted by resizing the canvas.
|
---|
94 | //
|
---|
95 | void MCamDisplay::Paint(Option_t *opt)
|
---|
96 | {
|
---|
97 | const UInt_t w = (UInt_t)(gPad->GetWw()*gPad->GetAbsWNDC());
|
---|
98 | const UInt_t h = (UInt_t)(gPad->GetWh()*gPad->GetAbsHNDC());
|
---|
99 |
|
---|
100 | //
|
---|
101 | // Check for a change in width or height, and make sure, that the
|
---|
102 | // first call also sets the range
|
---|
103 | //
|
---|
104 | if (w*fH == h*fW && fW && fH)
|
---|
105 | return;
|
---|
106 |
|
---|
107 | //
|
---|
108 | // Calculate aspect ratio (5/4=1.25 recommended)
|
---|
109 | //
|
---|
110 | const Double_t ratio = (Double_t)w/h;
|
---|
111 |
|
---|
112 | Float_t x;
|
---|
113 | Float_t y;
|
---|
114 |
|
---|
115 | if (ratio>1.25)
|
---|
116 | {
|
---|
117 | x = (ratio*2-1)*fRange;
|
---|
118 | y = fRange;
|
---|
119 | }
|
---|
120 | else
|
---|
121 | {
|
---|
122 | x = fRange*1.5;
|
---|
123 | y = fRange*1.25/ratio;
|
---|
124 | }
|
---|
125 |
|
---|
126 | fH = h;
|
---|
127 | fW = w;
|
---|
128 |
|
---|
129 | //
|
---|
130 | // Set new range
|
---|
131 | //
|
---|
132 | gPad->Range(-fRange, -y, x, y);
|
---|
133 | }
|
---|
134 |
|
---|
135 | // ------------------------------------------------------------------------
|
---|
136 | //
|
---|
137 | // Call this function to draw the camera layout into your canvas.
|
---|
138 | // Setup a drawing canvas. Add this object and all child objects
|
---|
139 | // (hexagons, etc) to the current pad. If no pad exists a new one is
|
---|
140 | // created.
|
---|
141 | //
|
---|
142 | void MCamDisplay::Draw(Option_t *option)
|
---|
143 | {
|
---|
144 | //
|
---|
145 | // if no canvas is yet existing to draw into, create a new one
|
---|
146 | //
|
---|
147 | if (!gPad)
|
---|
148 | fDrawingPad = new TCanvas("CamDisplay", "Magic Camera Display", 0, 0, 750, 600);
|
---|
149 | else
|
---|
150 | fDrawingPad = gPad;
|
---|
151 |
|
---|
152 | fDrawingPad->SetBorderMode(0);
|
---|
153 |
|
---|
154 | //
|
---|
155 | // Append this object, so that the aspect ratio is maintained
|
---|
156 | // (Paint-function is called)
|
---|
157 | //
|
---|
158 | AppendPad(option);
|
---|
159 |
|
---|
160 | //
|
---|
161 | // Setup the correct environment
|
---|
162 | //
|
---|
163 | gStyle->SetPalette(1, 0);
|
---|
164 |
|
---|
165 | gPad->SetFillColor(22);
|
---|
166 |
|
---|
167 | //
|
---|
168 | // Draw all pixels of the camera
|
---|
169 | // (means apend all pixelobjects to the current pad)
|
---|
170 | //
|
---|
171 | for (UInt_t i=0; i<fNumPixels; i++)
|
---|
172 | (*this)[i].Draw();
|
---|
173 |
|
---|
174 | //
|
---|
175 | // draw legend
|
---|
176 | //
|
---|
177 | const Float_t H = 0.9*fRange;
|
---|
178 | const Float_t h = 2./kItemsLegend;
|
---|
179 |
|
---|
180 | const Float_t w = fRange/sqrt(fNumPixels);
|
---|
181 |
|
---|
182 | for (Int_t i=0; i<kItemsLegend; i++)
|
---|
183 | {
|
---|
184 | TBox *box = GetBox(i);
|
---|
185 | box->SetX1(fRange);
|
---|
186 | box->SetX2(fRange+w);
|
---|
187 | box->SetY1(H*( i *h - 1.));
|
---|
188 | box->SetY2(H*((i+1)*h - 1.));
|
---|
189 | box->Draw();
|
---|
190 |
|
---|
191 | TText *txt = GetText(i);
|
---|
192 | txt->SetX(fRange+1.5*w);
|
---|
193 | txt->SetY(H*((i+0.5)*h - 1.));
|
---|
194 | txt->Draw();
|
---|
195 | }
|
---|
196 |
|
---|
197 | //fDrawingPad->SetEditable(kFALSE);
|
---|
198 | }
|
---|
199 |
|
---|
200 | // ------------------------------------------------------------------------
|
---|
201 | //
|
---|
202 | // Call this function to draw the number of photo electron into the
|
---|
203 | // camera.
|
---|
204 | //
|
---|
205 | void MCamDisplay::DrawPhotNum(const MCerPhotEvt *event)
|
---|
206 | {
|
---|
207 | if (!fDrawingPad)
|
---|
208 | Draw();
|
---|
209 |
|
---|
210 | fDrawingPad->cd();
|
---|
211 |
|
---|
212 | //
|
---|
213 | // Reset pixel colors to default value
|
---|
214 | //
|
---|
215 | Reset();
|
---|
216 |
|
---|
217 | //
|
---|
218 | // if the autoscale is true, set the values for the range for
|
---|
219 | // each event
|
---|
220 | //
|
---|
221 | if (fAutoScale)
|
---|
222 | {
|
---|
223 | fMinPhe = event->GetNumPhotonsMin();
|
---|
224 | fMaxPhe = event->GetNumPhotonsMax();
|
---|
225 |
|
---|
226 | if (fMaxPhe < 20.)
|
---|
227 | fMaxPhe = 20.;
|
---|
228 |
|
---|
229 | UpdateLegend();
|
---|
230 | }
|
---|
231 |
|
---|
232 | //
|
---|
233 | // update the colors in the picture
|
---|
234 | //
|
---|
235 | const Int_t entries = event->GetNumPixels();
|
---|
236 |
|
---|
237 | for (Int_t i=0; i<entries; i++)
|
---|
238 | {
|
---|
239 | const MCerPhotPix &pix = (*event)[i];
|
---|
240 |
|
---|
241 | if (!pix.IsPixelUsed())
|
---|
242 | continue;
|
---|
243 |
|
---|
244 | SetPixColor(pix);
|
---|
245 | }
|
---|
246 |
|
---|
247 | //
|
---|
248 | // Update display physically
|
---|
249 | //
|
---|
250 | gPad->Modified();
|
---|
251 | gPad->Update();
|
---|
252 | }
|
---|
253 |
|
---|
254 | // ------------------------------------------------------------------------
|
---|
255 | //
|
---|
256 | // reset the all pixel colors to a default value
|
---|
257 | //
|
---|
258 | void MCamDisplay::Reset()
|
---|
259 | {
|
---|
260 | for (UInt_t i=0; i<fNumPixels; i++)
|
---|
261 | (*this)[i].SetFillColor(10);
|
---|
262 | }
|
---|
263 |
|
---|
264 | // ------------------------------------------------------------------------
|
---|
265 | //
|
---|
266 | // Here we calculate the color index for the current value.
|
---|
267 | // The color index is defined with the class TStyle and the
|
---|
268 | // Color palette inside. We use the command gStyle->SetPalette(1,0)
|
---|
269 | // for the display. So we have to convert the value "wert" into
|
---|
270 | // a color index that fits the color palette.
|
---|
271 | // The range of the color palette is defined by the values fMinPhe
|
---|
272 | // and fMaxRange. Between this values we have 50 color index, starting
|
---|
273 | // with 0 up to 49.
|
---|
274 | //
|
---|
275 | Int_t MCamDisplay::GetColor(Float_t val)
|
---|
276 | {
|
---|
277 | //
|
---|
278 | // first treat the over- and under-flows
|
---|
279 | //
|
---|
280 | const Float_t maxcolidx = 49.0;
|
---|
281 |
|
---|
282 | if (val >= fMaxPhe)
|
---|
283 | return gStyle->GetColorPalette(maxcolidx);
|
---|
284 |
|
---|
285 | if (val <= fMinPhe)
|
---|
286 | return gStyle->GetColorPalette(0);
|
---|
287 |
|
---|
288 | //
|
---|
289 | // calculate the color index
|
---|
290 | //
|
---|
291 | const Float_t ratio = (val-fMinPhe) / (fMaxPhe-fMinPhe);
|
---|
292 | const Int_t colidx = (Int_t)(maxcolidx*ratio + .5);
|
---|
293 |
|
---|
294 | return gStyle->GetColorPalette(colidx);
|
---|
295 | }
|
---|
296 |
|
---|
297 | // ------------------------------------------------------------------------
|
---|
298 | //
|
---|
299 | // change the text on the legend according to the range of the
|
---|
300 | // Display
|
---|
301 | //
|
---|
302 | void MCamDisplay::UpdateLegend()
|
---|
303 | {
|
---|
304 | char text[10];
|
---|
305 |
|
---|
306 | for (Int_t i=0; i<kItemsLegend; i++)
|
---|
307 | {
|
---|
308 | const Float_t val = fMinPhe + (Float_t)i/kItemsLegend * (fMaxPhe-fMinPhe) ;
|
---|
309 |
|
---|
310 | sprintf(text, "%5.1f", val);
|
---|
311 |
|
---|
312 | TText &txt = *GetText(i);
|
---|
313 |
|
---|
314 | txt.SetText(txt.GetX(), txt.GetY(), text);
|
---|
315 | }
|
---|
316 | }
|
---|
317 |
|
---|
318 | // ------------------------------------------------------------------------
|
---|
319 | //
|
---|
320 | // Save primitive as a C++ statement(s) on output stream out
|
---|
321 | //
|
---|
322 | void MCamDisplay::SavePrimitive(ofstream &out, Option_t *opt)
|
---|
323 | {
|
---|
324 | if (!gROOT->ClassSaved(TCanvas::Class()))
|
---|
325 | fDrawingPad->SavePrimitive(out, opt);
|
---|
326 |
|
---|
327 | out << " " << fDrawingPad->GetName() << "->SetWindowSize(";
|
---|
328 | out << fDrawingPad->GetWw() << "," << fDrawingPad->GetWh() << ");" << endl;
|
---|
329 | }
|
---|