source: trunk/MagicSoft/Mars/mpedestal/MMcPedestalCopy.cc@ 8260

Last change on this file since 8260 was 7188, checked in by tbretz, 19 years ago
*** empty log message ***
File size: 4.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): Thomas Bretz, 12/2000 <mailto:tbretz@astro.uni-wuerzburg.de>
19!
20! Copyright: MAGIC Software Development, 2000-2003
21!
22!
23\* ======================================================================== */
24
25/////////////////////////////////////////////////////////////////////////////
26//
27// MMcPedestalCopy
28//
29// This task looks for the ìnformation about FADC pedestals in
30// MMcFadcHeader and translates it to the pedestal values in
31// MPedestalCam
32//
33// Input Containers:
34// MMcFadcHeader
35// [MMcRunHeader]
36// [MRawRunHeader]
37//
38// Output Containers:
39// MPedestalCam
40//
41/////////////////////////////////////////////////////////////////////////////
42#include "MMcPedestalCopy.h"
43
44#include "MParList.h"
45
46#include "MLog.h"
47#include "MLogManip.h"
48
49#include "MPedestalPix.h"
50#include "MPedestalCam.h"
51
52#include "MRawRunHeader.h"
53#include "MMcRunHeader.hxx"
54#include "MMcFadcHeader.hxx"
55
56ClassImp(MMcPedestalCopy);
57
58using namespace std;
59
60MMcPedestalCopy::MMcPedestalCopy(const char *name, const char *title)
61{
62 fName = name ? name : "MMcPedestalCopy";
63 fTitle = title ? title : "Copy MC pedestals into MPedestal Container";
64
65 //
66 // This is not needed here using MReadMarsFile because for the
67 // RunHeader tree the auto scheme is disabled by default
68 //
69 AddToBranchList("MMcFadcHeader.fPedesMean");
70 AddToBranchList("MMcFadcHeader.fElecNoise");
71}
72
73// --------------------------------------------------------------------------
74//
75// Check for the run type. Return kTRUE if it is a MC run or if there
76// is no MC run header (old camera files) kFALSE in case of a different
77// run type
78//
79Bool_t MMcPedestalCopy::CheckRunType(MParList *pList) const
80{
81 const MRawRunHeader *run = (MRawRunHeader*)pList->FindObject(AddSerialNumber("MRawRunHeader"));
82 if (!run)
83 {
84 *fLog << warn << dbginf << "Warning - cannot check file type, "
85 << AddSerialNumber("MRawRunHeader") << " not found." << endl;
86 return kTRUE;
87 }
88
89 return run->IsMonteCarloRun();
90}
91
92
93
94// --------------------------------------------------------------------------
95//
96// Make sure that there is a MPedestalCam object in the parameter list.
97//
98Int_t MMcPedestalCopy::PreProcess(MParList *pList)
99{
100 if ( ! pList->FindObject(AddSerialNumber("MPedestalCam")) )
101 pList->FindCreateObj(AddSerialNumber("MPedestalCam"));
102
103 return kTRUE;
104}
105
106// --------------------------------------------------------------------------
107//
108// Check for the runtype.
109// Search for MPedestalCam and MMcFadcHeader.
110//
111Bool_t MMcPedestalCopy::ReInit(MParList *pList)
112{
113 //
114 // If it is no MC file skip this function...
115 //
116 if (!CheckRunType(pList))
117 {
118 *fLog << inf << "This is no MC file... skipping." << endl;
119 return kTRUE;
120 }
121
122 //
123 // Now check the existance of all necessary containers. This has
124 // to be done only if this is a MC file.
125 //
126 MMcFadcHeader *fadc = (MMcFadcHeader*)pList->FindObject(AddSerialNumber("MMcFadcHeader"));
127 if (!fadc)
128 {
129 *fLog << err << dbginf << "MMcFadcHeader not found... aborting." << endl;
130 return kFALSE;
131 }
132
133 MPedestalCam *pedcam = (MPedestalCam*)pList->FindObject(AddSerialNumber("MPedestalCam"));
134 if (!pedcam)
135 {
136 *fLog << err << dbginf << "Cannot create " << AddSerialNumber("MPedestalCam") <<"... Exiting." << endl;
137
138 return kFALSE;
139 }
140
141 MMcRunHeader *mcrun = (MMcRunHeader*)pList->FindObject(AddSerialNumber("MMcRunHeader"));
142 if (!mcrun)
143 *fLog << warn << dbginf << AddSerialNumber("MMcRunHeader")
144 << " not found... assuming camera<0.7" << endl;
145
146 const int num = pedcam->GetSize();
147
148 const Bool_t camver70 = mcrun && mcrun->GetCamVersion()>=70;
149
150 for (int i=0; i<num; i++)
151 {
152 MPedestalPix &pix = (*pedcam)[i];
153
154 // Here one should compute the Pedestal taking into account how
155 // the MC makes the transformation analogic-digital for the FADC.
156 // This is done only once per file -> not time critical.
157
158 const Float_t pedest = fadc->GetPedestal(i);
159 const Float_t sigma = camver70 ? fadc->GetPedestalRmsHigh(i) : fadc->GetElecNoise(i);
160
161 pix.Set(pedest, sigma);
162
163 }
164
165 if (camver70)
166 pedcam->SetReadyToSave();
167
168 return kTRUE;
169}
Note: See TracBrowser for help on using the repository browser.