source: tags/Mars-V0.8.6/mbadpixels/MBadPixelsCalc.cc

Last change on this file was 4635, checked in by tbretz, 20 years ago
*** empty log message ***
File size: 7.2 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): Thomas Bretz, 02/2004 <mailto:tbretz@astro.uni.wuerzburg.de>
19!
20! Copyright: MAGIC Software Development, 2000-2004
21!
22!
23\* ======================================================================== */
24
25/////////////////////////////////////////////////////////////////////////////
26//
27// MBadPixelsCalc
28//
29//
30// The job of the task is to determin bad pixels event-wise. This must be
31// redone for each event. This particular task is for what is explained
32// below.
33// New algorithms may enter new tasks.
34//
35//
36// Check the pedestal RMS of every pixel with respect to the mean
37// pedestal RMS of the camera.
38//
39// The pixels can be set as blind if the pedestalRMS is too big or 0.
40//
41// If you don't want to use this option set the PedestalLevel<=0;
42//
43// MBadPixelsCalc calc;
44// calc.SetPedestalLevel(-1);
45//
46//
47// Input Containers:
48// [MPedPhotCam]
49// [MGeomCam]
50//
51// Output Containers:
52// MBadPixels
53//
54/////////////////////////////////////////////////////////////////////////////
55#include "MBadPixelsCalc.h"
56
57#include <TEnv.h>
58#include <TArrayD.h>
59
60#include "MLog.h"
61#include "MLogManip.h"
62
63#include "MParList.h"
64
65#include "MGeomCam.h"
66#include "MGeomPix.h"
67
68//#include "MSigmabar.h"
69
70#include "MPedPhotCam.h"
71#include "MPedPhotPix.h"
72
73#include "MBadPixelsCam.h"
74#include "MBadPixelsPix.h"
75
76ClassImp(MBadPixelsCalc);
77
78using namespace std;
79
80static const TString gsDefName = "MBadPixelsCalc";
81static const TString gsDefTitle = "Find hot spots (star, broken pixels, etc)";
82
83// --------------------------------------------------------------------------
84//
85// Default constructor.
86//
87MBadPixelsCalc::MBadPixelsCalc(const char *name, const char *title)
88 : fPedestalLevel(3), fNamePedPhotContainer("MPedPhotCam")
89{
90 fName = name ? name : gsDefName.Data();
91 fTitle = title ? title : gsDefTitle.Data();
92}
93
94// --------------------------------------------------------------------------
95//
96//
97Int_t MBadPixelsCalc::PreProcess (MParList *pList)
98{
99 fBadPixels = (MBadPixelsCam*)pList->FindCreateObj(AddSerialNumber("MBadPixelsCam"));
100 if (!fBadPixels)
101 return kFALSE;
102
103 if (fPedestalLevel>0)
104 {
105 fPedPhotCam = (MPedPhotCam*)pList->FindObject(AddSerialNumber(fNamePedPhotContainer), "MPedPhotCam");
106 if (!fPedPhotCam)
107 {
108 *fLog << err << "MPedPhotCam not found... aborting." << endl;
109 return kFALSE;
110 }
111
112 fGeomCam = (MGeomCam*)pList->FindObject(AddSerialNumber("MGeomCam"));
113 if (!fGeomCam)
114 {
115 *fLog << err << "MGeomCam not found... aborting." << endl;
116 return kFALSE;
117 }
118 }
119 return kTRUE;
120}
121
122// --------------------------------------------------------------------------
123//
124// Check the pedestal RMS of every pixel with respect to the mean pedestal RMS
125// of the camera (Sigmabar).
126//
127// The pixels can be set as blind if the pedestalRMS is too big or 0.
128//
129// If you don't want to use this option set the PedestalLevel<=0;
130//
131// MBadPixelsCalc calc;
132// calc.SetPedestalLevel(-1);
133/*
134void MBadPixelsCalc::CheckPedestalRMS() const
135{
136 const Int_t entries = fPedPhotCam->GetSize();
137
138 const Float_t meanPedRMS = fSigmabar->GetSigmabar();
139
140 for (Int_t i=0; i<entries; i++)
141 {
142 //
143 // get pixel as entry from list
144 //
145 const Double_t nratio = fGeomCam->GetPixRatio(i);
146 const Double_t pixPedRms = (*fPedPhotCam)[i].GetRms();
147
148 if (pixPedRms*nratio > fPedestalLevel * meanPedRMS || pixPedRms == 0)
149 (*fBadPixels)[i].SetUnsuitable(MBadPixelsPix::kUnsuitableEvt);
150 }
151}
152*/
153
154// --------------------------------------------------------------------------
155//
156// Check the pedestal Rms of the pixels: compute with 2 iterations the mean
157// for inner and outer pixels. Set as blind the pixels with too small or
158// too high pedestal Rms with respect to the mean.
159//
160Bool_t MBadPixelsCalc::CheckPedestalRms() const
161{
162 const Int_t entries = fPedPhotCam->GetSize();
163
164 const Int_t na = fGeomCam->GetNumAreas();
165
166 TArrayD meanrms(na);
167 TArrayI npix(na);
168
169 for (Int_t i=0; i<entries; i++)
170 {
171 const Double_t rms = (*fPedPhotCam)[i].GetRms();
172
173 if (rms<=0 || rms>=200*fGeomCam->GetPixRatioSqrt(i))
174 continue;
175
176 const Byte_t aidx = (*fGeomCam)[i].GetAidx();
177
178 meanrms[aidx] += rms;
179 npix[aidx]++;
180 }
181
182 //if no pixel has a minimum signal, return
183 for (int i=0; i<na; i++)
184 {
185 if (npix[i]==0 || meanrms[i]==0)
186 {
187 //fErrors[1]++; //no valid Pedestals Rms
188 return kFALSE;
189 }
190
191 meanrms[i] /= npix[i];
192 npix[i]=0;
193 }
194
195 TArrayD meanrms2(na);
196 for (Int_t i=0; i<entries; i++)
197 {
198 const Double_t rms = (*fPedPhotCam)[i].GetRms();
199 const Byte_t aidx = (*fGeomCam)[i].GetAidx();
200
201 //Calculate the corrected means:
202
203 if (rms<=0.5*meanrms[aidx] || rms>=1.5*meanrms[aidx])
204 continue;
205
206 meanrms2[aidx] += rms;
207 npix[aidx]++;
208 }
209
210 //if no pixel has a minimum signal, return
211 for (int i=0; i<na; i++)
212 {
213 if (npix[i]==0 || meanrms2[i]==0)
214 {
215 //fErrors[1]++; //no valid Pedestals Rms
216 return kFALSE;
217 }
218
219 meanrms2[i] /= npix[i];
220 }
221
222 Int_t bads = 0;
223
224 //Blind the Bad Pixels
225 for (Int_t i=0; i<entries; i++)
226 {
227 const Double_t rms = (*fPedPhotCam)[i].GetRms();
228 const Byte_t aidx = (*fGeomCam)[i].GetAidx();
229
230 if (rms>meanrms2[aidx]/fPedestalLevel && rms<=meanrms2[aidx]*fPedestalLevel)
231 continue;
232
233 (*fBadPixels)[i].SetUnsuitable(MBadPixelsPix::kUnsuitableEvt);
234 bads++;
235 }
236
237 // Check if the number of pixels to blind is > 60% of total number of pixels
238 //
239 // if (bads>0.6*entries)
240 // {
241 // fErrors[2]++;
242 // return kFALSE;
243 // }
244
245 return kTRUE;
246}
247
248
249// --------------------------------------------------------------------------
250//
251//
252Int_t MBadPixelsCalc::Process()
253{
254 if (fPedestalLevel>0)
255 {
256 CheckPedestalRms();
257 fPedPhotCam->ReCalc(*fGeomCam, fBadPixels);
258 }
259 return kTRUE;
260}
261
262// --------------------------------------------------------------------------
263//
264// Read the setup from a TEnv, eg:
265// MBadPixelsCalc.PedestalLevel: 3.0
266//
267Int_t MBadPixelsCalc::ReadEnv(const TEnv &env, TString prefix, Bool_t print)
268{
269 Bool_t rc = kFALSE;
270 if (IsEnvDefined(env, prefix, "PedestalLevel", print))
271 {
272 rc = kTRUE;
273 SetPedestalLevel(GetEnvValue(env, prefix, "PedestalLevel", fPedestalLevel));
274 }
275 return rc;
276}
Note: See TracBrowser for help on using the repository browser.