source: trunk/MagicSoft/Mars/mpedestal/MPedCalcPedRun.cc@ 5200

Last change on this file since 5200 was 4629, checked in by tbretz, 20 years ago
*** empty log message ***
File size: 20.4 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): Josep Flix 04/2001 <mailto:jflix@ifae.es>
19! Author(s): Thomas Bretz 05/2001 <mailto:tbretz@astro.uni-wuerzburg.de>
20! Author(s): Sebastian Commichau 12/2003
21! Author(s): Javier Rico 01/2004 <mailto:jrico@ifae.es>
22! Author(s): Markus Gaug 01/2004 <mailto:markus@ifae.es>
23!
24! Copyright: MAGIC Software Development, 2000-2004
25!
26!
27\* ======================================================================== */
28/////////////////////////////////////////////////////////////////////////////
29//
30// MPedCalcPedRun
31//
32// This task takes a pedestal run file and fills MPedestalCam during
33// the Process() with the pedestal and rms computed in an event basis.
34// In the PostProcess() MPedestalCam is finally filled with the pedestal
35// mean and rms computed in a run basis.
36// More than one run (file) can be merged
37//
38// MPedCalcPedRun applies the following formula (1):
39//
40// Pedestal per slice = sum(x_i) / n / slices
41// PedRMS per slice = Sqrt( ( sum(x_i^2) - sum(x_i)^2/n ) / n-1 / slices )
42//
43// where x_i is the sum of "slices" FADC slices and sum means the sum over all
44// events. "n" is the number of events, "slices" is the number of summed FADC samples.
45//
46// Note that the slice-to-slice fluctuations are not Gaussian, but Poissonian, thus
47// asymmetric and they are correlated.
48//
49// It is important to know that the Pedestal per slice and PedRMS per slice depend
50// on the number of used FADC slices, as seen in the following plots:
51//
52//Begin_Html
53/*
54<img src="images/PedestalStudyInner.gif">
55*/
56//End_Html
57//
58//Begin_Html
59/*
60<img src="images/PedestalStudyOuter.gif">
61*/
62//
63// The plots show the inner and outer pixels, respectivly and have the following meaning:
64//
65// 1) The calculated mean pedestal per slice (from MPedCalcPedRun)
66// 2) The fitted mean pedestal per slice (from MHPedestalCam)
67// 3) The calculated pedestal RMS per slice (from MPedCalcPedRun)
68// 4) The fitted sigma of the pedestal distribution per slice
69// (from MHPedestalCam)
70// 5) The relative difference between calculation and histogram fit
71// for the mean
72// 6) The relative difference between calculation and histogram fit
73// for the sigma or RMS, respectively.
74//
75// The calculated means do not change significantly except for the case of 2 slices,
76// however the RMS changes from 5.7 per slice in the case of 2 extracted slices
77// to 8.3 per slice in the case of 26 extracted slices. This change is very significant.
78//
79// The plots have been produced on run 20123. You can reproduce them using
80// the macro pedestalstudies.C
81//
82// Usage of this class:
83// ====================
84//
85// Call: SetRange(higainfirst, higainlast, logainfirst, logainlast)
86// to modify the ranges in which the window is allowed to move.
87// Defaults are:
88//
89// fHiGainFirst = fgHiGainFirst = 0
90// fHiGainLast = fgHiGainLast = 29
91// fLoGainFirst = fgLoGainFirst = 0
92// fLoGainLast = fgLoGainLast = 14
93//
94// Call: SetWindowSize(windowhigain, windowlogain)
95// to modify the sliding window widths. Windows have to be an even number.
96// In case of odd numbers, the window will be modified.
97//
98// Defaults are:
99//
100// fHiGainWindowSize = fgHiGainWindowSize = 14
101// fLoGainWindowSize = fgLoGainWindowSize = 0
102//
103//
104// ToDo:
105// - Take a setup file (ReadEnv-implementation) as input
106//
107//
108// Input Containers:
109// MRawEvtData
110// MRawRunHeader
111// MGeomCam
112//
113// Output Containers:
114// MPedestalCam
115//
116// See also: MPedestalCam, MPedestalPix, MHPedestalCam, MExtractor
117//
118/////////////////////////////////////////////////////////////////////////////
119#include "MPedCalcPedRun.h"
120#include "MExtractor.h"
121
122#include "MParList.h"
123
124#include "MLog.h"
125#include "MLogManip.h"
126
127#include "MRawRunHeader.h"
128#include "MRawEvtPixelIter.h"
129#include "MRawEvtData.h"
130
131#include "MPedestalPix.h"
132#include "MPedestalCam.h"
133
134#include "MGeomPix.h"
135#include "MGeomCam.h"
136
137ClassImp(MPedCalcPedRun);
138
139using namespace std;
140
141const Byte_t MPedCalcPedRun::fgHiGainFirst = 0;
142const Byte_t MPedCalcPedRun::fgHiGainLast = 29;
143const Byte_t MPedCalcPedRun::fgLoGainFirst = 0;
144const Byte_t MPedCalcPedRun::fgLoGainLast = 14;
145const Byte_t MPedCalcPedRun::fgHiGainWindowSize = 14;
146const Byte_t MPedCalcPedRun::fgLoGainWindowSize = 0;
147// --------------------------------------------------------------------------
148//
149// Default constructor:
150//
151// Sets:
152// - all pointers to NULL
153// - fWindowSizeHiGain to fgHiGainWindowSize
154// - fWindowSizeLoGain to fgLoGainWindowSize
155//
156// Calls:
157// - AddToBranchList("fHiGainPixId");
158// - AddToBranchList("fHiGainFadcSamples");
159// - SetRange(fgHiGainFirst, fgHiGainLast, fgLoGainFirst, fgLoGainLast)
160// - Clear()
161//
162MPedCalcPedRun::MPedCalcPedRun(const char *name, const char *title)
163 : fWindowSizeHiGain(fgHiGainWindowSize),
164 fWindowSizeLoGain(fgLoGainWindowSize),
165 fGeom(NULL)
166{
167 fName = name ? name : "MPedCalcPedRun";
168 fTitle = title ? title : "Task to calculate pedestals from pedestal runs raw data";
169
170 AddToBranchList("fHiGainPixId");
171 AddToBranchList("fHiGainFadcSamples");
172
173 SetRange(fgHiGainFirst, fgHiGainLast, fgLoGainFirst, fgLoGainLast);
174 Clear();
175}
176
177// --------------------------------------------------------------------------
178//
179// Sets:
180// - fNumSamplesTot to 0
181// - fRawEvt to NULL
182// - fRunHeader to NULL
183// - fPedestals to NULL
184//
185void MPedCalcPedRun::Clear(const Option_t *o)
186{
187 fNumSamplesTot = 0;
188
189 fRawEvt = NULL;
190 fRunHeader = NULL;
191 fPedestals = NULL;
192}
193
194// --------------------------------------------------------------------------
195//
196// SetRange:
197//
198// Calls:
199// - MExtractor::SetRange(hifirst,hilast,lofirst,lolast);
200// - SetWindowSize(fWindowSizeHiGain,fWindowSizeLoGain);
201//
202void MPedCalcPedRun::SetRange(Byte_t hifirst, Byte_t hilast, Byte_t lofirst, Byte_t lolast)
203{
204 MExtractor::SetRange(hifirst,hilast,lofirst,lolast);
205
206 //
207 // Redo the checks if the window is still inside the ranges
208 //
209 SetWindowSize(fWindowSizeHiGain,fWindowSizeLoGain);
210}
211
212// --------------------------------------------------------------------------
213//
214// Checks:
215// - if a window is odd, subtract one
216// - if a window is bigger than the one defined by the ranges, set it to the available range
217// - if a window is smaller than 2, set it to 2
218//
219// Sets:
220// - fNumHiGainSamples to: (Float_t)fWindowSizeHiGain
221// - fNumLoGainSamples to: (Float_t)fWindowSizeLoGain
222// - fSqrtHiGainSamples to: TMath::Sqrt(fNumHiGainSamples)
223// - fSqrtLoGainSamples to: TMath::Sqrt(fNumLoGainSamples)
224//
225void MPedCalcPedRun::SetWindowSize(Byte_t windowh, Byte_t windowl)
226{
227
228 fWindowSizeHiGain = windowh & ~1;
229 fWindowSizeLoGain = windowl & ~1;
230
231 if (fWindowSizeHiGain != windowh)
232 *fLog << warn << GetDescriptor() << ": Hi Gain window size has to be even, set to: "
233 << int(fWindowSizeHiGain) << " samples " << endl;
234
235 if (fWindowSizeLoGain != windowl)
236 *fLog << warn << GetDescriptor() << ": Lo Gain window size has to be even, set to: "
237 << int(fWindowSizeLoGain) << " samples " << endl;
238
239 const Byte_t availhirange = (fHiGainLast-fHiGainFirst+1) & ~1;
240 const Byte_t availlorange = (fLoGainLast-fLoGainFirst+1) & ~1;
241
242 if (fWindowSizeHiGain > availhirange)
243 {
244 *fLog << warn << GetDescriptor()
245 << Form("%s%2i%s%2i%s%2i%s",": Hi Gain window size: ",(int)fWindowSizeHiGain,
246 " is bigger than available range: [",(int)fHiGainFirst,",",(int)fHiGainLast,"]") << endl;
247 *fLog << warn << GetDescriptor()
248 << ": Will set window size to: " << (int)availhirange << endl;
249 fWindowSizeHiGain = availhirange;
250 }
251
252 if (fWindowSizeLoGain > availlorange)
253 {
254 *fLog << warn << GetDescriptor()
255 << Form("%s%2i%s%2i%s%2i%s",": Lo Gain window size: ",(int)fWindowSizeLoGain,
256 " is bigger than available range: [",(int)fLoGainFirst,",",(int)fLoGainLast,"]") << endl;
257 *fLog << warn << GetDescriptor()
258 << ": Will set window size to: " << (int)availlorange << endl;
259 fWindowSizeLoGain = availlorange;
260 }
261
262
263 fNumHiGainSamples = (Float_t)fWindowSizeHiGain;
264 fNumLoGainSamples = (Float_t)fWindowSizeLoGain;
265
266 fSqrtHiGainSamples = TMath::Sqrt(fNumHiGainSamples);
267 fSqrtLoGainSamples = TMath::Sqrt(fNumLoGainSamples);
268}
269
270// --------------------------------------------------------------------------
271//
272// Look for the following input containers:
273//
274// - MRawEvtData
275// - MRawRunHeader
276// - MGeomCam
277//
278// The following output containers are also searched and created if
279// they were not found:
280//
281// - MPedestalCam
282//
283Int_t MPedCalcPedRun::PreProcess( MParList *pList )
284{
285 Clear();
286
287 fRawEvt = (MRawEvtData*)pList->FindObject("MRawEvtData");
288 if (!fRawEvt)
289 {
290 *fLog << err << "MRawEvtData not found... aborting." << endl;
291 return kFALSE;
292 }
293
294 fRunHeader = (MRawRunHeader*)pList->FindObject(AddSerialNumber("MRawRunHeader"));
295 if (!fRunHeader)
296 {
297 *fLog << err << AddSerialNumber("MRawRunHeader") << " not found... aborting." << endl;
298 return kFALSE;
299 }
300
301 fGeom = (MGeomCam*)pList->FindObject("MGeomCam");
302 if (!fGeom)
303 {
304 *fLog << err << "MGeomCam not found... aborting." << endl;
305 return kFALSE;
306 }
307
308 fPedestals = (MPedestalCam*)pList->FindCreateObj("MPedestalCam");
309 if (!fPedestals)
310 return kFALSE;
311
312 return kTRUE;
313}
314
315// --------------------------------------------------------------------------
316//
317// The ReInit searches for:
318// - MRawRunHeader::GetNumSamplesHiGain()
319// - MRawRunHeader::GetNumSamplesLoGain()
320//
321// In case that the variables fHiGainLast and fLoGainLast are smaller than
322// the even part of the number of samples obtained from the run header, a
323// warning is given an the range is set back accordingly. A call to:
324// - SetRange(fHiGainFirst, fHiGainLast-diff, fLoGainFirst, fLoGainLast) or
325// - SetRange(fHiGainFirst, fHiGainLast, fLoGainFirst, fLoGainLast-diff)
326// is performed in that case. The variable diff means here the difference
327// between the requested range (fHiGainLast) and the available one. Note that
328// the functions SetRange() are mostly overloaded and perform more checks,
329// modifying the ranges again, if necessary.
330//
331Bool_t MPedCalcPedRun::ReInit(MParList *pList)
332{
333 Int_t lastdesired = (Int_t)fLoGainLast;
334 Int_t lastavailable = (Int_t)fRunHeader->GetNumSamplesLoGain()-1;
335
336 if (lastdesired > lastavailable)
337 {
338 const Int_t diff = lastdesired - lastavailable;
339 *fLog << endl;
340 *fLog << warn << GetDescriptor()
341 << Form("%s%2i%s%2i%s%2i%s",": Selected Lo Gain FADC Window [",
342 (int)fLoGainFirst,",",lastdesired,
343 "] ranges out of the available limits: [0,",lastavailable,"].") << endl;
344 *fLog << GetDescriptor() << ": Will reduce the upper edge to " << (int)(fLoGainLast - diff) << endl;
345 SetRange(fHiGainFirst, fHiGainLast, fLoGainFirst, fLoGainLast-diff);
346 }
347
348 lastdesired = (Int_t)fHiGainLast;
349 lastavailable = (Int_t)fRunHeader->GetNumSamplesHiGain()-1;
350
351 if (lastdesired > lastavailable)
352 {
353 const Int_t diff = lastdesired - lastavailable;
354 *fLog << endl;
355 *fLog << warn << GetDescriptor()
356 << Form("%s%2i%s%2i%s%2i%s",": Selected Hi Gain Range [",
357 (int)fHiGainFirst,",",lastdesired,
358 "] ranges out of the available limits: [0,",lastavailable,"].") << endl;
359 *fLog << warn << GetDescriptor()
360 << Form("%s%2i%s",": Will possibly use ",diff," samples from the Low-Gain for the High-Gain range")
361 << endl;
362 fHiGainLast -= diff;
363 fHiLoLast = diff;
364 }
365
366 lastdesired = (Int_t)fHiGainFirst+fWindowSizeHiGain-1;
367 lastavailable = (Int_t)fRunHeader->GetNumSamplesHiGain()-1;
368
369 if (lastdesired > lastavailable)
370 {
371 const Int_t diff = lastdesired - lastavailable;
372 *fLog << endl;
373 *fLog << warn << GetDescriptor()
374 << Form("%s%2i%s%2i%s",": Selected Hi Gain FADC Window size ",
375 (int)fWindowSizeHiGain,
376 " ranges out of the available limits: [0,",lastavailable,"].") << endl;
377 *fLog << warn << GetDescriptor()
378 << Form("%s%2i%s",": Will use ",diff," samples from the Low-Gain for the High-Gain extraction")
379 << endl;
380
381 if ((Int_t)fWindowSizeHiGain > diff)
382 {
383 fWindowSizeHiGain -= diff;
384 fWindowSizeLoGain += diff;
385 }
386 else
387 {
388 fWindowSizeLoGain += fWindowSizeHiGain;
389 fLoGainFirst = diff-fWindowSizeHiGain;
390 fWindowSizeHiGain = 0;
391 }
392 }
393
394
395 const Int_t npixels = fPedestals->GetSize();
396 const Int_t areas = fPedestals->GetAverageAreas();
397 const Int_t sectors = fPedestals->GetAverageSectors();
398
399 if (fSumx.GetSize()==0)
400 {
401 fSumx. Set(npixels);
402 fSumx2.Set(npixels);
403
404 fAreaSumx. Set(areas);
405 fAreaSumx2.Set(areas);
406 fAreaValid.Set(areas);
407
408 fSectorSumx. Set(sectors);
409 fSectorSumx2.Set(sectors);
410 fSectorValid.Set(sectors);
411
412 fSumx.Reset();
413 fSumx2.Reset();
414 }
415
416 if (fWindowSizeHiGain == 0 && fWindowSizeLoGain == 0)
417 {
418 *fLog << err << GetDescriptor()
419 << ": Number of extracted Slices is 0, abort ... " << endl;
420 return kFALSE;
421 }
422
423
424 *fLog << endl;
425 *fLog << inf << GetDescriptor() << ": Taking " << (int)fWindowSizeHiGain
426 << " HiGain FADC samples starting with slice: " << (int)fHiGainFirst << endl;
427 *fLog << inf << GetDescriptor() << ": Taking " << (int)fWindowSizeLoGain
428 << " LoGain FADC samples starting with slice: " << (int)fLoGainFirst << endl;
429
430 return kTRUE;
431}
432
433// --------------------------------------------------------------------------
434//
435// Fill the MPedestalCam container with the signal mean and rms for the event.
436// Store the measured signal in arrays fSumx and fSumx2 so that we can
437// calculate the overall mean and rms in the PostProcess()
438//
439Int_t MPedCalcPedRun::Process()
440{
441 MRawEvtPixelIter pixel(fRawEvt);
442
443 while (pixel.Next())
444 {
445 const UInt_t idx = pixel.GetPixelId();
446 const UInt_t aidx = (*fGeom)[idx].GetAidx();
447 const UInt_t sector = (*fGeom)[idx].GetSector();
448
449 Byte_t *ptr = pixel.GetHiGainSamples() + fHiGainFirst;
450 Byte_t *end = ptr + fWindowSizeHiGain;
451
452 UInt_t sum = 0;
453 UInt_t sqr = 0;
454
455 if (fWindowSizeHiGain != 0)
456 {
457 do
458 {
459 sum += *ptr;
460 sqr += *ptr * *ptr;
461 }
462 while (++ptr != end);
463 }
464
465 if (fWindowSizeLoGain != 0)
466 {
467
468 ptr = pixel.GetLoGainSamples() + fLoGainFirst;
469 end = ptr + fWindowSizeLoGain;
470
471 do
472 {
473 sum += *ptr;
474 sqr += *ptr * *ptr;
475 }
476 while (++ptr != end);
477
478 }
479
480 const Float_t msum = (Float_t)sum;
481
482 //
483 // These three lines have been uncommented by Markus Gaug
484 // If anybody needs them, please contact me!!
485 //
486 // const Float_t higainped = msum/fNumHiGainSlices;
487 // const Float_t higainrms = TMath::Sqrt((msqr-msum*msum/fNumHiGainSlices)/(fNumHiGainSlices-1.));
488 // (*fPedestals)[idx].Set(higainped, higainrms);
489
490 fSumx[idx] += msum;
491 fAreaSumx[aidx] += msum;
492 fSectorSumx[sector] += msum;
493 //
494 // The old version:
495 //
496 // const Float_t msqr = (Float_t)sqr;
497 // fSumx2[idx] += msqr;
498 //
499 // The new version:
500 //
501 const Float_t sqrsum = msum*msum;
502 fSumx2[idx] += sqrsum;
503 fAreaSumx2[aidx] += sqrsum;
504 fSectorSumx2[sector] += sqrsum;
505 }
506
507 fPedestals->SetReadyToSave();
508 fNumSamplesTot += fWindowSizeHiGain + fWindowSizeLoGain;
509
510 return kTRUE;
511}
512
513// --------------------------------------------------------------------------
514//
515// Compute signal mean and rms in the whole run and store it in MPedestalCam
516//
517Int_t MPedCalcPedRun::PostProcess()
518{
519 // Compute pedestals and rms from the whole run
520 const ULong_t n = fNumSamplesTot;
521 const ULong_t nevts = GetNumExecutions();
522
523 MRawEvtPixelIter pixel(fRawEvt);
524
525 while (pixel.Next())
526 {
527
528 const Int_t pixid = pixel.GetPixelId();
529 const UInt_t aidx = (*fGeom)[pixid].GetAidx();
530 const UInt_t sector = (*fGeom)[pixid].GetSector();
531
532 fAreaValid [aidx]++;
533 fSectorValid[sector]++;
534
535 const Float_t sum = fSumx.At(pixid);
536 const Float_t sum2 = fSumx2.At(pixid);
537 const Float_t higainped = sum/n;
538 //
539 // The old version:
540 //
541 // const Float_t higainrms = TMath::Sqrt((sum2-sum*sum/n)/(n-1.));
542 //
543 // The new version:
544 //
545 // 1. Calculate the Variance of the sums:
546 Float_t higainVar = (sum2-sum*sum/nevts)/(nevts-1.);
547 // 2. Scale the variance to the number of slices:
548 higainVar /= (Float_t)(fWindowSizeHiGain+fWindowSizeLoGain);
549 // 3. Calculate the RMS from the Variance:
550 (*fPedestals)[pixid].Set(higainped, higainVar < 0 ? 0. : TMath::Sqrt(higainVar));
551
552 }
553
554 //
555 // Loop over the (two) area indices to get the averaged pedestal per aidx
556 //
557 for (Int_t aidx=0; aidx<fAreaValid.GetSize(); aidx++)
558 {
559
560 const Int_t napix = fAreaValid.At(aidx);
561
562 if (napix == 0)
563 continue;
564
565 const Float_t sum = fAreaSumx.At(aidx);
566 const Float_t sum2 = fAreaSumx2.At(aidx);
567 const ULong_t an = napix * n;
568 const ULong_t aevts = napix * nevts;
569
570 const Float_t higainped = sum/an;
571
572 // 1. Calculate the Variance of the sums:
573 Float_t higainVar = (sum2-sum*sum/aevts)/(aevts-1.);
574 // 2. Scale the variance to the number of slices:
575 higainVar /= fWindowSizeHiGain+fWindowSizeLoGain;
576 // 3. Calculate the RMS from the Variance:
577 Float_t higainrms = TMath::Sqrt(higainVar);
578 // 4. Re-scale it with the square root of the number of involved pixels
579 // in order to be comparable to the mean of pedRMS of that area
580 higainrms *= TMath::Sqrt((Float_t)napix);
581
582 fPedestals->GetAverageArea(aidx).Set(higainped, higainrms);
583 }
584
585 //
586 // Loop over the (six) sector indices to get the averaged pedestal per sector
587 //
588 for (Int_t sector=0; sector<fSectorValid.GetSize(); sector++)
589 {
590
591 const Int_t nspix = fSectorValid.At(sector);
592
593 if (nspix == 0)
594 continue;
595
596 const Float_t sum = fSectorSumx.At(sector);
597 const Float_t sum2 = fSectorSumx2.At(sector);
598 const ULong_t sn = nspix * n;
599 const ULong_t sevts = nspix * nevts;
600
601 const Float_t higainped = sum/sn;
602
603 // 1. Calculate the Variance of the sums:
604 Float_t higainVar = (sum2-sum*sum/sevts)/(sevts-1.);
605 // 2. Scale the variance to the number of slices:
606 higainVar /= fWindowSizeHiGain+fWindowSizeLoGain;
607 // 3. Calculate the RMS from the Variance:
608 Float_t higainrms = TMath::Sqrt(higainVar);
609 // 4. Re-scale it with the square root of the number of involved pixels
610 // in order to be comparable to the mean of pedRMS of that sector
611 higainrms *= TMath::Sqrt((Float_t)nspix);
612
613 fPedestals->GetAverageSector(sector).Set(higainped, higainrms);
614 }
615
616 fPedestals->SetTotalEntries(fNumSamplesTot);
617 fPedestals->SetReadyToSave();
618
619 return kTRUE;
620}
621
622Int_t MPedCalcPedRun::ReadEnv(const TEnv &env, TString prefix, Bool_t print)
623{
624 if (MExtractor::ReadEnv(env, prefix, print)==kERROR)
625 return kERROR;
626
627 Byte_t hw = fWindowSizeHiGain;
628 Byte_t lw = fWindowSizeLoGain;
629
630 Bool_t rc = kFALSE;
631
632 if (IsEnvDefined(env, prefix, "WindowSizeHiGain", print))
633 {
634 hw = GetEnvValue(env, prefix, "WindowSizeHiGain", hw);
635 rc=kTRUE;
636 }
637
638 if (IsEnvDefined(env, prefix, "WindowSizeLoGain", print))
639 {
640 lw = GetEnvValue(env, prefix, "WindowSizeLoGain", lw);
641 rc=kTRUE;
642 }
643
644 if (rc)
645 SetWindowSize(hw, lw);
646
647 return rc;
648}
Note: See TracBrowser for help on using the repository browser.