source: trunk/MagicSoft/Mars/mcalib/MCalibrationCam.cc@ 4136

Last change on this file since 4136 was 3936, checked in by gaug, 21 years ago
*** empty log message ***
File size: 11.3 KB
Line 
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
62ClassImp(MCalibrationCam);
63
64using namespace std;
65
66const 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//
89MCalibrationCam::MCalibrationCam(const char *name, const char *title)
90 : fNumUnsuitable(),
91 fNumUnreliable(),
92 fNumHiGainFADCSlices(0.), fNumLoGainFADCSlices(0.), fPulserColor(kNONE),
93 fPixels(NULL), fAverageAreas(NULL), fAverageSectors(NULL)
94{
95 fName = name ? name : "MCalibrationCam";
96 fTitle = title ? title : "Base class Storage container for Camera Calibration";
97
98 fAverageBadAreas = new TClonesArray("MBadPixelsPix",1);
99 fAverageBadSectors = new TClonesArray("MBadPixelsPix",1);
100
101}
102
103// --------------------------------------------------------------------------
104//
105// Deletes the following TClonesArray's of MCalibrationPix containers (if exist):
106// - fPixels
107// - fAverageAreas
108// - fAverageSectors
109//
110// Deletes the following TClonesArray's of MBadPixelsPix containers (if exist):
111// - fAverageBadAreas
112// - fAverageBadSectors
113//
114MCalibrationCam::~MCalibrationCam()
115{
116
117 //
118 // delete fPixels should delete all Objects stored inside
119 //
120 if (fPixels)
121 delete fPixels;
122
123 if (fAverageAreas)
124 delete fAverageAreas;
125
126 if (fAverageSectors)
127 delete fAverageSectors;
128
129 delete fAverageBadAreas;
130 delete fAverageBadSectors;
131
132}
133
134// --------------------------------------
135//
136// Calls the ForEach macro for the TClonesArray fPixels with the argument Clear()
137//
138// Loops over the fAverageAreas, calling the function Clear() for
139// every entry in:
140// - fAverageAreas
141// - fAverageBadAreas
142//
143// Loops over the fAverageSectors, calling the function Clear() for
144// every entry in:
145// - fAverageSectors
146// - fAverageBadSectors
147//
148void MCalibrationCam::Clear(Option_t *o)
149{
150
151 fPixels->ForEach(TObject, Clear)();
152
153 //
154 // another ForEach does not compile, thus have to do the loop ourselves:
155 //
156 for (Int_t i=0;i<GetAverageAreas();i++)
157 {
158 fAverageAreas[i].Clear();
159 fAverageBadAreas[i].Clear();
160 }
161
162 //
163 // another ForEach does not compile, thus have to do the loop ourselves:
164 //
165 for (Int_t i=0;i<GetAverageSectors();i++)
166 {
167 fAverageSectors[i].Clear();
168 fAverageBadSectors[i].Clear();
169 }
170
171 return;
172}
173
174// -------------------------------------------------------------------
175//
176// Calls TClonesArray::ExpandCreate() for fPixels
177//
178void MCalibrationCam::InitSize(const UInt_t i)
179{
180 fPixels->ExpandCreate(i);
181}
182
183// -------------------------------------------------------------------
184//
185// Calls TClonesArray::ExpandCreate() for:
186// - fAverageAreas
187// - fAverageBadAreas
188//
189void MCalibrationCam::InitAverageAreas(const UInt_t i)
190{
191 fAverageAreas->ExpandCreate(i);
192 fAverageBadAreas->ExpandCreate(i);
193
194 fNumUnsuitable.Set(i);
195 fNumUnreliable.Set(i);
196}
197
198// -------------------------------------------------------------------
199//
200// Calls TClonesArray::ExpandCreate() for:
201// - fAverageSectors
202// - fAverageBadSectors
203//
204void MCalibrationCam::InitAverageSectors(const UInt_t i)
205{
206 fAverageSectors->ExpandCreate(i);
207 fAverageBadSectors->ExpandCreate(i);
208}
209
210// -------------------------------------------------------------------
211//
212// Calls:
213// - InitSize()
214// - InitAverageAreas()
215// - InitAverageSectors()
216//
217void MCalibrationCam::Init(const MGeomCam &geom)
218{
219 InitSize (geom.GetNumPixels() );
220 InitAverageAreas (geom.GetNumAreas() );
221 InitAverageSectors(geom.GetNumSectors());
222}
223
224// --------------------------------------------------------------------------
225//
226// Returns the number of un-suitable pixels per area index and -1 if
227// the area index exceeds the initialized array.
228//
229const Int_t MCalibrationCam::GetNumUnsuitable( Int_t aidx ) const
230{
231 if (aidx < 0)
232 return -1;
233
234 return aidx > fNumUnsuitable.GetSize() ? -1 : fNumUnsuitable[aidx];
235}
236
237// --------------------------------------------------------------------------
238//
239// Returns the number of un-reliable pixels per area index and -1 if
240// the area index exceeds the initialized array.
241//
242const Int_t MCalibrationCam::GetNumUnreliable( Int_t aidx ) const
243{
244 if (aidx < 0)
245 return -1;
246
247 return aidx > fNumUnreliable.GetSize() ? -1 : fNumUnreliable[aidx];
248}
249
250
251// --------------------------------------------------------------------------
252//
253// Returns the current size of the TClonesArray fAverageAreas
254// independently if the MCalibrationPix is filled with values or not.
255//
256const Int_t MCalibrationCam::GetAverageAreas() const
257{
258 return fAverageAreas->GetEntriesFast();
259}
260
261// --------------------------------------------------------------------------
262//
263// Returns the current size of the TClonesArray fAverageSectors
264// independently if the MCalibrationPix is filled with values or not.
265//
266const Int_t MCalibrationCam::GetAverageSectors() const
267{
268 return fAverageSectors->GetEntriesFast();
269}
270
271
272// --------------------------------------------------------------------------
273//
274// Get i-th pixel (pixel number)
275//
276MCalibrationPix &MCalibrationCam::operator[](UInt_t i)
277{
278 return *static_cast<MCalibrationPix*>(fPixels->UncheckedAt(i));
279}
280
281// --------------------------------------------------------------------------
282//
283// Get i-th pixel (pixel number)
284//
285const MCalibrationPix &MCalibrationCam::operator[](UInt_t i) const
286{
287 return *static_cast<MCalibrationPix*>(fPixels->UncheckedAt(i));
288}
289
290// Returns the current size of the TClonesArray fPixels
291// independently if the MCalibrationPix is filled with values or not.
292const Int_t MCalibrationCam::GetSize() const
293{
294 return fPixels->GetEntriesFast();
295}
296
297// --------------------------------------------------------------------------
298//
299// Get i-th average pixel (area number)
300//
301MCalibrationPix &MCalibrationCam::GetAverageArea(UInt_t i)
302{
303 return *static_cast<MCalibrationPix*>(fAverageAreas->UncheckedAt(i));
304}
305
306// --------------------------------------------------------------------------
307//
308// Get i-th average pixel (area number)
309//
310const MCalibrationPix &MCalibrationCam::GetAverageArea(UInt_t i) const
311{
312 return *static_cast<MCalibrationPix*>(fAverageAreas->UncheckedAt(i));
313}
314
315// --------------------------------------------------------------------------
316//
317// Get i-th average pixel (sector number)
318//
319MCalibrationPix &MCalibrationCam::GetAverageSector(UInt_t i)
320{
321 return *static_cast<MCalibrationPix*>(fAverageSectors->UncheckedAt(i));
322}
323
324// --------------------------------------------------------------------------
325//
326// Get i-th average pixel (sector number)
327//
328const MCalibrationPix &MCalibrationCam::GetAverageSector(UInt_t i) const
329{
330 return *static_cast<MCalibrationPix*>(fAverageSectors->UncheckedAt(i));
331}
332
333// --------------------------------------------------------------------------
334//
335// Get i-th average pixel (area number)
336//
337MBadPixelsPix &MCalibrationCam::GetAverageBadArea(UInt_t i)
338{
339 return *static_cast<MBadPixelsPix*>(fAverageBadAreas->UncheckedAt(i));
340}
341
342// --------------------------------------------------------------------------
343//
344// Get i-th average pixel (area number)
345//
346const MBadPixelsPix &MCalibrationCam::GetAverageBadArea(UInt_t i) const
347{
348 return *static_cast<MBadPixelsPix*>(fAverageBadAreas->UncheckedAt(i));
349}
350
351// --------------------------------------------------------------------------
352//
353// Get i-th average pixel (sector number)
354//
355MBadPixelsPix &MCalibrationCam::GetAverageBadSector(UInt_t i)
356{
357 return *static_cast<MBadPixelsPix*>(fAverageBadSectors->UncheckedAt(i));
358}
359
360// --------------------------------------------------------------------------
361//
362// Get i-th average pixel (sector number)
363//
364const MBadPixelsPix &MCalibrationCam::GetAverageBadSector(UInt_t i) const
365{
366 return *static_cast<MBadPixelsPix*>(fAverageBadSectors->UncheckedAt(i));
367}
368
369
370
371// --------------------------------------------------------------------------
372//
373// Dummy needed for compilation with MCamEvent
374//
375Bool_t MCalibrationCam::GetPixelContent(Double_t &val, Int_t idx, const MGeomCam &cam, Int_t type) const
376{
377 return kTRUE;
378}
379
380
381
382// --------------------------------------------------------------------------
383//
384// Calls MCalibrationPix::DrawClone()
385//
386void MCalibrationCam::DrawPixelContent(Int_t idx) const
387{
388 (*this)[idx].DrawClone();
389}
390
391void MCalibrationCam::SetNumUnsuitable( const UInt_t i, const Int_t aidx)
392{
393 if (aidx < 0)
394 return;
395
396 if (aidx < fNumUnsuitable.GetSize())
397 fNumUnsuitable[aidx] = i;
398}
399
400void MCalibrationCam::SetNumUnreliable( const UInt_t i, const Int_t aidx)
401{
402 if (aidx < 0)
403 return;
404 if (aidx < fNumUnreliable.GetSize())
405 fNumUnreliable[aidx] = i;
406}
Note: See TracBrowser for help on using the repository browser.