source: trunk/MagicSoft/Mars/manalysis/MCerPhotCalc.cc@ 2247

Last change on this file since 2247 was 2237, checked in by tbretz, 22 years ago
*** empty log message ***
File size: 6.9 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): Abelardo Moralejo 7/2002 <mailto:moralejo@pd.infn.it>
19! Author(s): Thomas Bretz 2002 <mailto:tbretz@astro.uni-wuerzburg.de>
20!
21! Copyright: MAGIC Software Development, 2002-2003
22!
23!
24\* ======================================================================== */
25
26//////////////////////////////////////////////////////////////////////////////
27//
28// MCerPhotCalc
29//
30// This is a task which calculates the number of photons from the FADC
31// time slices. It weights the each slice according to the numbers in
32// the array fWeight (default: all slices added up with weight 1).
33//
34// The weights are rescaled, such that sum(weigths)=num slices
35//
36// Input Containers:
37// MRawRunHeader, MRawEvtData, MPedestalCam
38//
39// Output Containers:
40// MCerPhotEvt
41//
42//////////////////////////////////////////////////////////////////////////////
43#include "MCerPhotCalc.h"
44
45#include "MParList.h"
46
47#include "MLog.h"
48#include "MLogManip.h"
49
50#include "MMcRunHeader.hxx"
51
52#include "MRawRunHeader.h"
53#include "MRawEvtData.h" // MRawEvtData::GetNumPixels
54#include "MCerPhotEvt.h"
55#include "MPedestalPix.h"
56#include "MPedestalCam.h"
57#include "MRawEvtPixelIter.h"
58
59ClassImp(MCerPhotCalc);
60
61using namespace std;
62
63// --------------------------------------------------------------------------
64//
65// Default constructor.
66//
67MCerPhotCalc::MCerPhotCalc(const char *name, const char *title)
68{
69 fName = name ? name : "MCerPhotCalc";
70 fTitle = title ? title : "Calculate pixel signal from FADC data";
71
72 AddToBranchList("MRawEvtData.fHiGainPixId");
73 AddToBranchList("MRawEvtData.fLoGainPixId");
74 AddToBranchList("MRawEvtData.fHiGainFadcSamples");
75 AddToBranchList("MRawEvtData.fLoGainFadcSamples");
76
77 SetDefaultWeights();
78}
79
80// --------------------------------------------------------------------------
81//
82// The PreProcess searches for the following input containers:
83// - MRawRunHeader
84// - MRawEvtData
85// - MPedestalCam
86//
87// The following output containers are also searched and created if
88// they were not found:
89// - MCerPhotEvt
90//
91Int_t MCerPhotCalc::PreProcess(MParList *pList)
92{
93 fRunHeader = (MRawRunHeader*)pList->FindObject("MRawRunHeader");
94 if (!fRunHeader)
95 {
96 *fLog << dbginf << "MRawRunHeader not found... aborting." << endl;
97 return kFALSE;
98 }
99
100 fRawEvt = (MRawEvtData*)pList->FindObject("MRawEvtData");
101 if (!fRawEvt)
102 {
103 *fLog << dbginf << "MRawEvtData not found... aborting." << endl;
104 return kFALSE;
105 }
106
107 fPedestals = (MPedestalCam*)pList->FindCreateObj("MPedestalCam");
108 if (!fPedestals)
109 {
110 *fLog << dbginf << "MPedestalCam not found... aborting." << endl;
111 return kFALSE;
112 }
113
114 fCerPhotEvt = (MCerPhotEvt*)pList->FindCreateObj("MCerPhotEvt");
115 if (!fCerPhotEvt)
116 return kFALSE;
117
118 // Calculate quadratic sum of weights:
119 fSumWeights = 0;
120 fSumQuadWeights = 0;
121 for (Int_t i=0; i<fWeight.GetSize(); i++)
122 {
123 fSumWeights += fWeight[i];
124 fSumQuadWeights += fWeight[i]*fWeight[i];
125 }
126
127 fSumQuadWeights = sqrt(fSumQuadWeights);
128
129 return kTRUE;
130}
131
132// --------------------------------------------------------------------------
133//
134// Check for the run type and camera version.
135// If the file is a MC file and the used camera version is <= 40
136// we enable a fix for truncated pedestal means in this version.
137//
138Bool_t MCerPhotCalc::ReInit(MParList *pList)
139{
140 const MRawRunHeader *runheader = (MRawRunHeader*)pList->FindObject("MRawRunHeader");
141 if (!runheader)
142 {
143 *fLog << warn << dbginf << "Warning - cannot check file type, MRawRunHeader not found." << endl;
144 return kTRUE;
145 }
146
147 if (fRunHeader->GetNumSamplesHiGain() != fWeight.GetSize())
148 {
149 *fLog << dbginf << "Number of FADC slices (" << fRunHeader->GetNumSamplesHiGain() <<") is different from assumed one (" << fWeight.GetSize() << ")... aborting." << endl;
150 return kFALSE;
151 }
152
153 Bool_t fIsMcFile = runheader->GetRunType() == kRTMonteCarlo;
154 if (!fIsMcFile)
155 return kTRUE;
156
157 ScalePedestals();
158
159 MMcRunHeader *mcrunheader = (MMcRunHeader*)pList->FindObject("MMcRunHeader");
160 if (!mcrunheader)
161 {
162 *fLog << warn << dbginf << "Warning - cannot check for camera version, MC run header not found." << endl;
163 return kTRUE;
164 }
165
166 fEnableFix = kFALSE;
167 if (mcrunheader->GetCamVersion() <= 40)
168 fEnableFix = kTRUE;
169
170 return kTRUE;
171}
172
173void MCerPhotCalc::ScalePedestals()
174{
175 Int_t n = 577; //FIXME: fPedestals->GetNumPixel();
176
177 for (int idx=0; idx<n; idx++)
178 {
179 MPedestalPix &ped = (*fPedestals)[idx];
180
181 const Double_t offset = fEnableFix ? ped.GetPedestal()-0.5 : ped.GetPedestal();
182 ped.Set(offset*fSumWeights, ped.GetPedestalRms()*fSumQuadWeights);
183 }
184
185 fPedestals->SetReadyToSave();
186}
187
188// --------------------------------------------------------------------------
189//
190// Calculate the integral of the FADC time slices and store them as a new
191// pixel in the MCerPhotEvt container.
192//
193Int_t MCerPhotCalc::Process()
194{
195 //fCerPhotEvt->InitSize(fRawEvt->GetNumPixels());
196
197 if (fIsMcFile)
198 ScalePedestals();
199
200 MRawEvtPixelIter pixel(fRawEvt);
201
202 while (pixel.Next())
203 {
204 const UInt_t idx = pixel.GetPixelId();
205 MPedestalPix &ped = (*fPedestals)[idx];
206
207 //
208 // sanity check (old MC files sometimes have pixids>577)
209 //
210 if (!fPedestals->CheckBounds(idx))
211 {
212 *fLog << inf << "Pixel Index larger than camera... skipping event." << endl;
213 return kCONTINUE;
214 }
215
216 //
217 // Calculate pixel signal unless it has all FADC slices empty:
218 //
219 Byte_t *ptr = pixel.GetHiGainSamples();
220
221 Float_t nphot = 0;
222 for(Int_t i=0; i<fWeight.GetSize(); i++)
223 nphot += ptr[i]*fWeight[i];
224
225 nphot -= ped.GetPedestal();
226
227 fCerPhotEvt->AddPixel(idx, nphot, 0);
228
229 // FIXME! Handling of Lo Gains is missing!
230 }
231
232 fCerPhotEvt->FixSize();
233 fCerPhotEvt->SetReadyToSave();
234
235 return kTRUE;
236}
237
238// --------------------------------------------------------------------------
239//
240// Set default values for the number of slices and weights:
241//
242void MCerPhotCalc::SetDefaultWeights()
243{
244 const Float_t dummy[15] = { 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 };
245 fWeight.Set(15, dummy);
246}
Note: See TracBrowser for help on using the repository browser.