source: tags/Mars-V0.5/manalysis/MMcPedestalCopy.cc

Last change on this file was 1035, checked in by tbretz, 23 years ago
*** empty log message ***
File size: 3.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 12/2000 (tbretz@uni-sw.gwdg.de)
19!
20! Copyright: MAGIC Software Development, 2000-2001
21!
22!
23\* ======================================================================== */
24
25/////////////////////////////////////////////////////////////////////////////
26// //
27// MMcPedestalCopy //
28// //
29// Input Containers: //
30// MMcFadcHeader //
31// //
32// Output Containers: //
33// MPedestalCam //
34// //
35/////////////////////////////////////////////////////////////////////////////
36
37#include "MMcPedestalCopy.h"
38
39#include "MParList.h"
40
41#include "MLog.h"
42#include "MLogManip.h"
43
44#include "MPedestalCam.h"
45#include "MRawRunHeader.h"
46#include "MMcFadcHeader.hxx"
47
48ClassImp(MMcPedestalCopy);
49
50MMcPedestalCopy::MMcPedestalCopy(const char *name, const char *title)
51{
52 fName = name ? name : "MMcPedestalCopy";
53 fTitle = title ? title : "Task to copy monte carlo pedestals into MPedestal Container";
54
55 AddToBranchList("fPedestalMean");
56 AddToBranchList("fElecNoise");
57}
58
59Bool_t MMcPedestalCopy::PreProcess(MParList *pList)
60{
61 MRawRunHeader *run = (MRawRunHeader*)pList->FindObject("MRawRunHeader");
62 if (run)
63 {
64 if (run->GetRunType() != kRTMonteCarlo)
65 return kTRUE;
66 }
67
68 const MMcFadcHeader *mcped = (MMcFadcHeader*)pList->FindObject("MMcFadcHeader");
69 if (!mcped)
70 {
71 *fLog << dbginf << "MMcFadcHeader not found... aborting." << endl;
72 return kFALSE;
73 }
74
75 MPedestalCam *pedcam = (MPedestalCam*)pList->FindCreateObj("MPedestalCam");
76 if (!pedcam)
77 return kFALSE;
78
79 const int num = mcped->GetNumPixel();
80
81 pedcam->InitSize(num);
82
83 for (int i=0; i<num; i++)
84 {
85 MPedestalPix &pix = (*pedcam)[i];
86
87 const Float_t pedest = mcped->GetPedestal(i);
88 const Float_t pedrms = mcped->GetPedestalRms(i);
89
90 const Float_t sigma = pedest*sqrt(num);
91 const Float_t sigrms = sigma/sqrt(num*2);
92
93 pix.SetPedestal(pedest, sigma);
94 pix.SetPedestalRms(pedrms, sigrms);
95 }
96
97 return kTRUE;
98}
99
Note: See TracBrowser for help on using the repository browser.