1 | /* ======================================================================== *\
|
---|
2 | !
|
---|
3 | ! *
|
---|
4 | ! * This file is part of MARS, the MAGIC Analysis and Reconstruction
|
---|
5 | ! * Software. It is distributed to you in the hope that it can be a useful
|
---|
6 | ! * and timesaving tool in analysing Data of imaging Cerenkov telescopes.
|
---|
7 | ! * It is distributed WITHOUT ANY WARRANTY.
|
---|
8 | ! *
|
---|
9 | ! * Permission to use, copy, modify and distribute this software and its
|
---|
10 | ! * documentation for any purpose is hereby granted without fee,
|
---|
11 | ! * provided that the above copyright notice appear in all copies and
|
---|
12 | ! * that both that copyright notice and this permission notice appear
|
---|
13 | ! * in supporting documentation. It is provided "as is" without express
|
---|
14 | ! * or implied warranty.
|
---|
15 | ! *
|
---|
16 | !
|
---|
17 | !
|
---|
18 | ! Author(s): Markus Gaug 11/2003 <mailto:markus@ifae.es>
|
---|
19 | !
|
---|
20 | ! Copyright: MAGIC Software Development, 2000-2004
|
---|
21 | !
|
---|
22 | !
|
---|
23 | \* ======================================================================== */
|
---|
24 |
|
---|
25 | /////////////////////////////////////////////////////////////////////////////
|
---|
26 | //
|
---|
27 | // MCalibrationCam
|
---|
28 | //
|
---|
29 | // Base class for Camera Calibration results.
|
---|
30 | //
|
---|
31 | // Contains TClonesArrays for the following objects:
|
---|
32 | // - fPixels: Array of classes derived from MCalibrationPix, one entry
|
---|
33 | // per pixel. Has to be created
|
---|
34 | // - fAverageAreas: Array of classes derived from MCalibrationPix, one entry
|
---|
35 | // per pixel AREA. Has to be created
|
---|
36 | // - fAverageSectors: Array of classes derived from MCalibrationPix, one entry
|
---|
37 | // per camera SECTOR. Has to be created
|
---|
38 | //
|
---|
39 | // - fAverageBadAreas: Array of classes derived from MBadPixelsPix, one entry
|
---|
40 | // per pixel AREA. Is created automatically.
|
---|
41 | // - fAverageBadSectors: Array of classes derived from MBadPixelsPix, one entry
|
---|
42 | // per camera SECTOR. Is created automatically.
|
---|
43 | //
|
---|
44 | // All TClonesArrays have to enlarged by the corresponding calls to (e.g. in MGeomApply):
|
---|
45 | // - InitSize()
|
---|
46 | // - InitAverageAreas()
|
---|
47 | // - InitAverageSectors()
|
---|
48 | //
|
---|
49 | /////////////////////////////////////////////////////////////////////////////
|
---|
50 | #include "MCalibrationCam.h"
|
---|
51 |
|
---|
52 | #include <TClonesArray.h>
|
---|
53 |
|
---|
54 | #include "MGeomCam.h"
|
---|
55 | #include "MGeomPix.h"
|
---|
56 |
|
---|
57 | #include "MBadPixelsCam.h"
|
---|
58 | #include "MBadPixelsPix.h"
|
---|
59 |
|
---|
60 | #include "MCalibrationPix.h"
|
---|
61 |
|
---|
62 | ClassImp(MCalibrationCam);
|
---|
63 |
|
---|
64 | using namespace std;
|
---|
65 |
|
---|
66 | const Int_t MCalibrationCam::gkNumPulserColors = 4;
|
---|
67 | // --------------------------------------------------------------------------
|
---|
68 | //
|
---|
69 | // Default constructor.
|
---|
70 | //
|
---|
71 | // Set the following pointer to NULL:
|
---|
72 | // - fPixels
|
---|
73 | // - fAverageAreas
|
---|
74 | // - fAverageSectors
|
---|
75 | //
|
---|
76 | // Initializes:
|
---|
77 | // - fPulserColor to kNONE
|
---|
78 | // - fNumHiGainFADCSlices to 0.
|
---|
79 | // - fNumLoGainFADCSlices to 0.
|
---|
80 | //
|
---|
81 | // Creates a TClonesArray of MBadPixelsPix containers for the TClonesArray's:
|
---|
82 | // - fAverageBadAreas
|
---|
83 | // - fAverageBadSectors
|
---|
84 | // all initialized to 1 entry
|
---|
85 | //
|
---|
86 | // Later, a call to InitAverageAreas() and InitAverageSectors() (or Init())
|
---|
87 | // has to be performed in order to get the dimension correctly.
|
---|
88 | //
|
---|
89 | MCalibrationCam::MCalibrationCam(const char *name, const char *title)
|
---|
90 | : fNumHiGainFADCSlices(0.), fNumLoGainFADCSlices(0.), fPulserColor(kNONE),
|
---|
91 | fPixels(NULL), fAverageAreas(NULL), fAverageSectors(NULL)
|
---|
92 | {
|
---|
93 | fName = name ? name : "MCalibrationCam";
|
---|
94 | fTitle = title ? title : "Base class Storage container for Camera Calibration";
|
---|
95 |
|
---|
96 | fAverageBadAreas = new TClonesArray("MBadPixelsPix",1);
|
---|
97 | fAverageBadSectors = new TClonesArray("MBadPixelsPix",1);
|
---|
98 |
|
---|
99 | }
|
---|
100 |
|
---|
101 | // --------------------------------------------------------------------------
|
---|
102 | //
|
---|
103 | // Deletes the following TClonesArray's of MCalibrationPix containers (if exist):
|
---|
104 | // - fPixels
|
---|
105 | // - fAverageAreas
|
---|
106 | // - fAverageSectors
|
---|
107 | //
|
---|
108 | // Deletes the following TClonesArray's of MBadPixelsPix containers (if exist):
|
---|
109 | // - fAverageBadAreas
|
---|
110 | // - fAverageBadSectors
|
---|
111 | //
|
---|
112 | MCalibrationCam::~MCalibrationCam()
|
---|
113 | {
|
---|
114 |
|
---|
115 | //
|
---|
116 | // delete fPixels should delete all Objects stored inside
|
---|
117 | //
|
---|
118 | if (fPixels)
|
---|
119 | delete fPixels;
|
---|
120 |
|
---|
121 | if (fAverageAreas)
|
---|
122 | delete fAverageAreas;
|
---|
123 |
|
---|
124 | if (fAverageSectors)
|
---|
125 | delete fAverageSectors;
|
---|
126 |
|
---|
127 | delete fAverageBadAreas;
|
---|
128 | delete fAverageBadSectors;
|
---|
129 |
|
---|
130 | }
|
---|
131 |
|
---|
132 | // --------------------------------------
|
---|
133 | //
|
---|
134 | // Calls the ForEach macro for the TClonesArray fPixels with the argument Clear()
|
---|
135 | //
|
---|
136 | // Loops over the fAverageAreas, calling the function Clear() for
|
---|
137 | // every entry in:
|
---|
138 | // - fAverageAreas
|
---|
139 | // - fAverageBadAreas
|
---|
140 | //
|
---|
141 | // Loops over the fAverageSectors, calling the function Clear() for
|
---|
142 | // every entry in:
|
---|
143 | // - fAverageSectors
|
---|
144 | // - fAverageBadSectors
|
---|
145 | //
|
---|
146 | void MCalibrationCam::Clear(Option_t *o)
|
---|
147 | {
|
---|
148 |
|
---|
149 | fPixels->ForEach(TObject, Clear)();
|
---|
150 |
|
---|
151 | //
|
---|
152 | // another ForEach does not compile, thus have to do the loop ourselves:
|
---|
153 | //
|
---|
154 | for (Int_t i=0;i<GetAverageAreas();i++)
|
---|
155 | {
|
---|
156 | fAverageAreas[i].Clear();
|
---|
157 | fAverageBadAreas[i].Clear();
|
---|
158 | }
|
---|
159 |
|
---|
160 | //
|
---|
161 | // another ForEach does not compile, thus have to do the loop ourselves:
|
---|
162 | //
|
---|
163 | for (Int_t i=0;i<GetAverageSectors();i++)
|
---|
164 | {
|
---|
165 | fAverageSectors[i].Clear();
|
---|
166 | fAverageBadSectors[i].Clear();
|
---|
167 | }
|
---|
168 |
|
---|
169 | return;
|
---|
170 | }
|
---|
171 |
|
---|
172 | void MCalibrationCam::Copy(TObject& object) const
|
---|
173 | {
|
---|
174 |
|
---|
175 | MCalibrationCam &calib = (MCalibrationCam&)object;
|
---|
176 |
|
---|
177 | MParContainer::Copy(object);
|
---|
178 |
|
---|
179 | calib.fNumHiGainFADCSlices = fNumHiGainFADCSlices;
|
---|
180 | calib.fNumLoGainFADCSlices = fNumLoGainFADCSlices;
|
---|
181 | calib.fPulserColor = fPulserColor;
|
---|
182 |
|
---|
183 | const Int_t n3 = GetSize();
|
---|
184 | if (n3 != 0)
|
---|
185 | {
|
---|
186 | calib.InitSize(n3);
|
---|
187 | for (int i=0; i<n3; i++)
|
---|
188 | (*this)[i].Copy(calib[i]);
|
---|
189 | }
|
---|
190 |
|
---|
191 | const Int_t n4 = GetAverageAreas();
|
---|
192 | if (n4 != 0)
|
---|
193 | {
|
---|
194 | calib.InitAverageAreas(n4);
|
---|
195 | for (int i=0; i<n4; i++)
|
---|
196 | GetAverageArea(i).Copy(calib.GetAverageArea(i));
|
---|
197 | for (int i=0; i<n4; i++)
|
---|
198 | GetAverageBadArea(i).Copy(calib.GetAverageBadArea(i));
|
---|
199 | for (int i=0; i<n4; i++)
|
---|
200 | calib.fNumUnsuitable[i] = fNumUnsuitable[i];
|
---|
201 | for (int i=0; i<n4; i++)
|
---|
202 | calib.fNumUnreliable[i] = fNumUnreliable[i];
|
---|
203 | }
|
---|
204 |
|
---|
205 | const Int_t n5 = GetAverageSectors();
|
---|
206 | if (n5 != 0)
|
---|
207 | {
|
---|
208 | calib.InitAverageSectors(n5);
|
---|
209 | for (int i=0; i<n5; i++)
|
---|
210 | GetAverageSector(i).Copy(calib.GetAverageSector(i));
|
---|
211 | for (int i=0; i<n5; i++)
|
---|
212 | GetAverageBadSector(i).Copy(calib.GetAverageBadSector(i));
|
---|
213 | }
|
---|
214 |
|
---|
215 | }
|
---|
216 |
|
---|
217 |
|
---|
218 | // -------------------------------------------------------------------
|
---|
219 | //
|
---|
220 | // Calls TClonesArray::ExpandCreate() for fPixels
|
---|
221 | //
|
---|
222 | void MCalibrationCam::InitSize(const UInt_t i)
|
---|
223 | {
|
---|
224 | fPixels->ExpandCreate(i);
|
---|
225 | }
|
---|
226 |
|
---|
227 | // -------------------------------------------------------------------
|
---|
228 | //
|
---|
229 | // Calls TClonesArray::ExpandCreate() for:
|
---|
230 | // - fAverageAreas
|
---|
231 | // - fAverageBadAreas
|
---|
232 | //
|
---|
233 | void MCalibrationCam::InitAverageAreas(const UInt_t i)
|
---|
234 | {
|
---|
235 | fAverageAreas->ExpandCreate(i);
|
---|
236 | fAverageBadAreas->ExpandCreate(i);
|
---|
237 |
|
---|
238 | for (UInt_t j=0; j<i; j++)
|
---|
239 | GetAverageArea(j).SetPixId(j);
|
---|
240 |
|
---|
241 | fNumUnsuitable.Set(i);
|
---|
242 | fNumUnreliable.Set(i);
|
---|
243 | }
|
---|
244 |
|
---|
245 | // -------------------------------------------------------------------
|
---|
246 | //
|
---|
247 | // Calls TClonesArray::ExpandCreate() for:
|
---|
248 | // - fAverageSectors
|
---|
249 | // - fAverageBadSectors
|
---|
250 | //
|
---|
251 | void MCalibrationCam::InitAverageSectors(const UInt_t i)
|
---|
252 | {
|
---|
253 | fAverageSectors->ExpandCreate(i);
|
---|
254 | fAverageBadSectors->ExpandCreate(i);
|
---|
255 |
|
---|
256 | for (UInt_t j=0; j<i; j++)
|
---|
257 | GetAverageSector(j).SetPixId(j);
|
---|
258 |
|
---|
259 | }
|
---|
260 |
|
---|
261 | // -------------------------------------------------------------------
|
---|
262 | //
|
---|
263 | // Calls:
|
---|
264 | // - InitSize()
|
---|
265 | // - InitAverageAreas()
|
---|
266 | // - InitAverageSectors()
|
---|
267 | //
|
---|
268 | void MCalibrationCam::Init(const MGeomCam &geom)
|
---|
269 | {
|
---|
270 | InitSize (geom.GetNumPixels() );
|
---|
271 | InitAverageAreas (geom.GetNumAreas() );
|
---|
272 | InitAverageSectors(geom.GetNumSectors());
|
---|
273 | }
|
---|
274 |
|
---|
275 | // --------------------------------------------------------------------------
|
---|
276 | //
|
---|
277 | // Returns the number of un-suitable pixels per area index and -1 if
|
---|
278 | // the area index exceeds the initialized array.
|
---|
279 | //
|
---|
280 | const Int_t MCalibrationCam::GetNumUnsuitable( Int_t aidx ) const
|
---|
281 | {
|
---|
282 | if (aidx < 0)
|
---|
283 | return -1;
|
---|
284 |
|
---|
285 | return aidx > fNumUnsuitable.GetSize() ? -1 : fNumUnsuitable[aidx];
|
---|
286 | }
|
---|
287 |
|
---|
288 | // --------------------------------------------------------------------------
|
---|
289 | //
|
---|
290 | // Returns the number of un-reliable pixels per area index and -1 if
|
---|
291 | // the area index exceeds the initialized array.
|
---|
292 | //
|
---|
293 | const Int_t MCalibrationCam::GetNumUnreliable( Int_t aidx ) const
|
---|
294 | {
|
---|
295 | if (aidx < 0)
|
---|
296 | return -1;
|
---|
297 |
|
---|
298 | return aidx > fNumUnreliable.GetSize() ? -1 : fNumUnreliable[aidx];
|
---|
299 | }
|
---|
300 |
|
---|
301 |
|
---|
302 | // --------------------------------------------------------------------------
|
---|
303 | //
|
---|
304 | // Returns the current size of the TClonesArray fAverageAreas
|
---|
305 | // independently if the MCalibrationPix is filled with values or not.
|
---|
306 | //
|
---|
307 | const Int_t MCalibrationCam::GetAverageAreas() const
|
---|
308 | {
|
---|
309 | return fAverageAreas->GetEntriesFast();
|
---|
310 | }
|
---|
311 |
|
---|
312 | // --------------------------------------------------------------------------
|
---|
313 | //
|
---|
314 | // Returns the current size of the TClonesArray fAverageSectors
|
---|
315 | // independently if the MCalibrationPix is filled with values or not.
|
---|
316 | //
|
---|
317 | const Int_t MCalibrationCam::GetAverageSectors() const
|
---|
318 | {
|
---|
319 | return fAverageSectors->GetEntriesFast();
|
---|
320 | }
|
---|
321 |
|
---|
322 |
|
---|
323 | // --------------------------------------------------------------------------
|
---|
324 | //
|
---|
325 | // Get i-th pixel (pixel number)
|
---|
326 | //
|
---|
327 | MCalibrationPix &MCalibrationCam::operator[](UInt_t i)
|
---|
328 | {
|
---|
329 | return *static_cast<MCalibrationPix*>(fPixels->UncheckedAt(i));
|
---|
330 | }
|
---|
331 |
|
---|
332 | // --------------------------------------------------------------------------
|
---|
333 | //
|
---|
334 | // Get i-th pixel (pixel number)
|
---|
335 | //
|
---|
336 | const MCalibrationPix &MCalibrationCam::operator[](UInt_t i) const
|
---|
337 | {
|
---|
338 | return *static_cast<MCalibrationPix*>(fPixels->UncheckedAt(i));
|
---|
339 | }
|
---|
340 |
|
---|
341 | // Returns the current size of the TClonesArray fPixels
|
---|
342 | // independently if the MCalibrationPix is filled with values or not.
|
---|
343 | const Int_t MCalibrationCam::GetSize() const
|
---|
344 | {
|
---|
345 | return fPixels->GetEntriesFast();
|
---|
346 | }
|
---|
347 |
|
---|
348 | // --------------------------------------------------------------------------
|
---|
349 | //
|
---|
350 | // Get i-th average pixel (area number)
|
---|
351 | //
|
---|
352 | MCalibrationPix &MCalibrationCam::GetAverageArea(UInt_t i)
|
---|
353 | {
|
---|
354 | return *static_cast<MCalibrationPix*>(fAverageAreas->UncheckedAt(i));
|
---|
355 | }
|
---|
356 |
|
---|
357 | // --------------------------------------------------------------------------
|
---|
358 | //
|
---|
359 | // Get i-th average pixel (area number)
|
---|
360 | //
|
---|
361 | const MCalibrationPix &MCalibrationCam::GetAverageArea(UInt_t i) const
|
---|
362 | {
|
---|
363 | return *static_cast<MCalibrationPix*>(fAverageAreas->UncheckedAt(i));
|
---|
364 | }
|
---|
365 |
|
---|
366 | // --------------------------------------------------------------------------
|
---|
367 | //
|
---|
368 | // Get i-th average pixel (sector number)
|
---|
369 | //
|
---|
370 | MCalibrationPix &MCalibrationCam::GetAverageSector(UInt_t i)
|
---|
371 | {
|
---|
372 | return *static_cast<MCalibrationPix*>(fAverageSectors->UncheckedAt(i));
|
---|
373 | }
|
---|
374 |
|
---|
375 | // --------------------------------------------------------------------------
|
---|
376 | //
|
---|
377 | // Get i-th average pixel (sector number)
|
---|
378 | //
|
---|
379 | const MCalibrationPix &MCalibrationCam::GetAverageSector(UInt_t i) const
|
---|
380 | {
|
---|
381 | return *static_cast<MCalibrationPix*>(fAverageSectors->UncheckedAt(i));
|
---|
382 | }
|
---|
383 |
|
---|
384 | // --------------------------------------------------------------------------
|
---|
385 | //
|
---|
386 | // Get i-th average pixel (area number)
|
---|
387 | //
|
---|
388 | MBadPixelsPix &MCalibrationCam::GetAverageBadArea(UInt_t i)
|
---|
389 | {
|
---|
390 | return *static_cast<MBadPixelsPix*>(fAverageBadAreas->UncheckedAt(i));
|
---|
391 | }
|
---|
392 |
|
---|
393 | // --------------------------------------------------------------------------
|
---|
394 | //
|
---|
395 | // Get i-th average pixel (area number)
|
---|
396 | //
|
---|
397 | const MBadPixelsPix &MCalibrationCam::GetAverageBadArea(UInt_t i) const
|
---|
398 | {
|
---|
399 | return *static_cast<MBadPixelsPix*>(fAverageBadAreas->UncheckedAt(i));
|
---|
400 | }
|
---|
401 |
|
---|
402 | // --------------------------------------------------------------------------
|
---|
403 | //
|
---|
404 | // Get i-th average pixel (sector number)
|
---|
405 | //
|
---|
406 | MBadPixelsPix &MCalibrationCam::GetAverageBadSector(UInt_t i)
|
---|
407 | {
|
---|
408 | return *static_cast<MBadPixelsPix*>(fAverageBadSectors->UncheckedAt(i));
|
---|
409 | }
|
---|
410 |
|
---|
411 | // --------------------------------------------------------------------------
|
---|
412 | //
|
---|
413 | // Get i-th average pixel (sector number)
|
---|
414 | //
|
---|
415 | const MBadPixelsPix &MCalibrationCam::GetAverageBadSector(UInt_t i) const
|
---|
416 | {
|
---|
417 | return *static_cast<MBadPixelsPix*>(fAverageBadSectors->UncheckedAt(i));
|
---|
418 | }
|
---|
419 |
|
---|
420 |
|
---|
421 |
|
---|
422 | // --------------------------------------------------------------------------
|
---|
423 | //
|
---|
424 | // Dummy needed for compilation with MCamEvent
|
---|
425 | //
|
---|
426 | Bool_t MCalibrationCam::GetPixelContent(Double_t &val, Int_t idx, const MGeomCam &cam, Int_t type) const
|
---|
427 | {
|
---|
428 | return kTRUE;
|
---|
429 | }
|
---|
430 |
|
---|
431 |
|
---|
432 |
|
---|
433 | // --------------------------------------------------------------------------
|
---|
434 | //
|
---|
435 | // Calls MCalibrationPix::DrawClone()
|
---|
436 | //
|
---|
437 | void MCalibrationCam::DrawPixelContent(Int_t idx) const
|
---|
438 | {
|
---|
439 | (*this)[idx].DrawClone();
|
---|
440 | }
|
---|
441 |
|
---|
442 | void MCalibrationCam::SetNumUnsuitable( const UInt_t i, const Int_t aidx)
|
---|
443 | {
|
---|
444 | if (aidx < 0)
|
---|
445 | return;
|
---|
446 |
|
---|
447 | if (aidx < fNumUnsuitable.GetSize())
|
---|
448 | fNumUnsuitable[aidx] = i;
|
---|
449 | }
|
---|
450 |
|
---|
451 | void MCalibrationCam::SetNumUnreliable( const UInt_t i, const Int_t aidx)
|
---|
452 | {
|
---|
453 | if (aidx < 0)
|
---|
454 | return;
|
---|
455 | if (aidx < fNumUnreliable.GetSize())
|
---|
456 | fNumUnreliable[aidx] = i;
|
---|
457 | }
|
---|