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 | // MMcRunHeader
|
---|
28 | //
|
---|
29 | // Root storage container for the RUN MONTE CARLO HEADER information
|
---|
30 | //
|
---|
31 | // This the second version of this output class. Old root files, which have
|
---|
32 | // a previous version of this class, are still compatibles and can be used.
|
---|
33 | // But of course, you can no try to get infromatino in these old files about
|
---|
34 | // the new data members.
|
---|
35 | //
|
---|
36 | // The following data member have been introduced in this second version
|
---|
37 | // and they do not exist in the previous one:
|
---|
38 | //
|
---|
39 | // Float_t fMcRunNumber; Run Number
|
---|
40 | // UInt_t fProductionSite; code to know where the run was generated
|
---|
41 | // Float_t fDateRunMMCs; Date of the MMCs production
|
---|
42 | // Float_t fDateRunCamera; Date, when the Camera program is run.
|
---|
43 | // Byte_t fRawEvt; RawEvt[Data,Hedaer] stored or not
|
---|
44 | // Byte_t fElecNoise; Electronic Noise simulated or not
|
---|
45 | // Byte_t fStarFieldRotate; Is the starfield rotation switched on (1) or
|
---|
46 | // off (0)
|
---|
47 | // Float_t fCWaveLower; Wavelength limits for the Cerenkov photons
|
---|
48 | // Float_t fCWaveUpper;
|
---|
49 | // UInt_t fNumObsLev; Observation levels
|
---|
50 | // Float_t fHeightLev[10];
|
---|
51 | // Float_t fSlopeSpec; Spectral index
|
---|
52 | //
|
---|
53 | // Third version:
|
---|
54 | //
|
---|
55 | // Byte_t fOpticLinksNoise; Flag to state if the optic noise is simualted or not
|
---|
56 | //
|
---|
57 | //
|
---|
58 | // Note: All the azimuth Phi angles in this and other MC classes follow
|
---|
59 | // the convention in the Corsika program (see Corsika manual). There, phi
|
---|
60 | // is the azimuth of the momentum vector of particles, and is measured
|
---|
61 | // from the north direction, anticlockwise (i.e, west is phi=90 degrees).
|
---|
62 | // When it refers to the telescope orientation, it is the azimuth of a
|
---|
63 | // vector along the telescope axis, going from the camera to the mirror.
|
---|
64 | // So, fTelesTheta=90, fTelesPhi = 0 means the telescope is pointing
|
---|
65 | // horizontally towards South. For an explanation, see also TDAS 02-11.
|
---|
66 | //
|
---|
67 | ////////////////////////////////////////////////////////////////////////////
|
---|
68 |
|
---|
69 | #include "MMcRunHeader.hxx"
|
---|
70 |
|
---|
71 | #include <fstream>
|
---|
72 | #include <iomanip>
|
---|
73 |
|
---|
74 | #include "MLog.h"
|
---|
75 |
|
---|
76 | ClassImp(MMcRunHeader);
|
---|
77 |
|
---|
78 | using namespace std;
|
---|
79 |
|
---|
80 | // --------------------------------------------------------------------------
|
---|
81 | //
|
---|
82 | // Default constructor.
|
---|
83 | //
|
---|
84 | //
|
---|
85 | MMcRunHeader::MMcRunHeader(const char *name, const char *title)
|
---|
86 | {
|
---|
87 | fName = name ? name : "MMcRunHeader";
|
---|
88 | fTitle = title ? title : "Raw Run Header Information";
|
---|
89 |
|
---|
90 | fMcRunNumber =0;
|
---|
91 | fProductionSite = 0;
|
---|
92 | fDateRunMMCs = 0;
|
---|
93 | fDateRunCamera = 0;
|
---|
94 | fNumTrigCond = 0;
|
---|
95 | fAllEvtsTriggered = 0 ;
|
---|
96 | fMcEvt = 0;
|
---|
97 | fMcTrig = 0;
|
---|
98 | fMcFadc = 0;
|
---|
99 | fRawEvt = 0;
|
---|
100 | fElecNoise = 0;
|
---|
101 | fStarFieldRotate = 0;
|
---|
102 | fNumAnalisedPixels = 0;
|
---|
103 | fNumSimulatedShowers = 0;
|
---|
104 | fNumStoredShowers = 0;
|
---|
105 | fNumEvents = 0;
|
---|
106 |
|
---|
107 | fStarFieldRaH = 0;
|
---|
108 | fStarFieldRaM = 0;
|
---|
109 | fStarFieldRaS = 0;
|
---|
110 | fStarFieldDeD = 0;
|
---|
111 | fStarFieldDeM = 0;
|
---|
112 | fStarFieldDeS = 0;
|
---|
113 |
|
---|
114 | fNumPheFromDNSB = 0.0;
|
---|
115 | fTelesTheta = 0.0;
|
---|
116 | fTelesPhi = 0.0;
|
---|
117 | fShowerThetaMax = 0.0;
|
---|
118 | fShowerThetaMin = 0.0;
|
---|
119 | fShowerPhiMax = 0.0;
|
---|
120 | fShowerPhiMin = 0.0;
|
---|
121 |
|
---|
122 | fCWaveLower = 0.0;
|
---|
123 | fCWaveUpper = 0.0;
|
---|
124 |
|
---|
125 | fNumObsLev = 0;
|
---|
126 | for (int i=0; i<10; i++){
|
---|
127 | fHeightLev[i]=0.0;
|
---|
128 | }
|
---|
129 | fSlopeSpec = 0.0;
|
---|
130 |
|
---|
131 | fCorsikaVersion = 0;
|
---|
132 | fReflVersion = 0;
|
---|
133 | fCamVersion = 0;
|
---|
134 | fOpticLinksNoise= 0;
|
---|
135 |
|
---|
136 | }
|
---|
137 |
|
---|
138 | // --------------------------------------------------------------------------
|
---|
139 | //
|
---|
140 | // Destructor.
|
---|
141 | //
|
---|
142 | MMcRunHeader::~MMcRunHeader()
|
---|
143 | {
|
---|
144 |
|
---|
145 | }
|
---|
146 |
|
---|
147 | // -------------------------------------------------------------------------
|
---|
148 | //
|
---|
149 | // Fill. Put data in the container
|
---|
150 | //
|
---|
151 | void MMcRunHeader::Fill(const Float_t runnumber,
|
---|
152 | const UInt_t productionsite,
|
---|
153 | const Float_t daterunMMCs,
|
---|
154 | const Float_t daterunCamera,
|
---|
155 | const UInt_t numtrigcond,
|
---|
156 | const Byte_t allevts,
|
---|
157 | const Byte_t mcevt,
|
---|
158 | const Byte_t mctrig,
|
---|
159 | const Byte_t mcfadc,
|
---|
160 | const Byte_t rawevt,
|
---|
161 | const Byte_t elecnoise,
|
---|
162 | const Int_t numanalpixels,
|
---|
163 | const UInt_t numsim,
|
---|
164 | const UInt_t numsto,
|
---|
165 | const Byte_t starfieldrotate,
|
---|
166 | const Int_t sfRaH,
|
---|
167 | const Int_t sfRaM,
|
---|
168 | const Int_t sfRaS,
|
---|
169 | const Int_t sfDeD,
|
---|
170 | const Int_t sfDeM,
|
---|
171 | const Int_t sfDeS,
|
---|
172 | const Float_t numdnsb,
|
---|
173 | const Float_t telestheta,
|
---|
174 | const Float_t telesphi,
|
---|
175 | const Float_t shthetamax,
|
---|
176 | const Float_t shthetamin,
|
---|
177 | const Float_t shphimax,
|
---|
178 | const Float_t shphimin,
|
---|
179 | const Float_t impactmax,
|
---|
180 | const Float_t cwavelower,
|
---|
181 | const Float_t cwaveupper,
|
---|
182 | const Float_t slopespec,
|
---|
183 | const UInt_t numObslev,
|
---|
184 | const Float_t heightlev[10],
|
---|
185 | const UInt_t corsika,
|
---|
186 | const UInt_t refl,
|
---|
187 | const UInt_t cam,
|
---|
188 | const Byte_t opticnoise)
|
---|
189 | {
|
---|
190 | fMcRunNumber =runnumber;
|
---|
191 | fProductionSite = productionsite;
|
---|
192 | fDateRunMMCs = daterunMMCs;
|
---|
193 | fDateRunCamera = daterunCamera;
|
---|
194 | fNumTrigCond = numtrigcond;
|
---|
195 | fAllEvtsTriggered = allevts;
|
---|
196 | fMcEvt = mcevt;
|
---|
197 | fMcTrig = mctrig;
|
---|
198 | fMcFadc = mcfadc;
|
---|
199 | fRawEvt = rawevt;
|
---|
200 | fElecNoise = elecnoise;
|
---|
201 | fStarFieldRotate = starfieldrotate;
|
---|
202 | fNumAnalisedPixels = numanalpixels;
|
---|
203 | fNumSimulatedShowers = numsim;
|
---|
204 | fNumStoredShowers = numsto;
|
---|
205 | fNumEvents = numsto;
|
---|
206 |
|
---|
207 | fStarFieldRaH = sfRaH;
|
---|
208 | fStarFieldRaM = sfRaM;
|
---|
209 | fStarFieldRaS = sfRaS;
|
---|
210 | fStarFieldDeD = sfDeD;
|
---|
211 | fStarFieldDeM = sfDeM;
|
---|
212 | fStarFieldDeS = sfDeS;
|
---|
213 |
|
---|
214 | fNumPheFromDNSB = numdnsb;
|
---|
215 | fTelesTheta = telestheta;
|
---|
216 | fTelesPhi = telesphi;
|
---|
217 | fShowerThetaMax = shthetamax;
|
---|
218 | fShowerThetaMin = shthetamin;
|
---|
219 | fShowerPhiMax = shphimax;
|
---|
220 | fShowerPhiMin = shphimin;
|
---|
221 |
|
---|
222 | fImpactMax=impactmax;
|
---|
223 |
|
---|
224 | fCWaveLower = cwavelower;
|
---|
225 | fCWaveUpper = cwaveupper;
|
---|
226 |
|
---|
227 | fNumObsLev = numObslev;
|
---|
228 | for (UInt_t i=0; i<numObslev; i++){
|
---|
229 | fHeightLev[i]=heightlev[i];
|
---|
230 | }
|
---|
231 | fSlopeSpec = slopespec;
|
---|
232 |
|
---|
233 | fCorsikaVersion = corsika;
|
---|
234 | fReflVersion = refl;
|
---|
235 | fCamVersion = cam;
|
---|
236 | fOpticLinksNoise= opticnoise;
|
---|
237 | }
|
---|
238 |
|
---|
239 | // -------------------------------------------------------------------------
|
---|
240 | //
|
---|
241 | // GetStarFieldRa. Get RA coordinates of the starfield
|
---|
242 | //
|
---|
243 | void MMcRunHeader::GetStarFieldRa(Int_t *hour, Int_t *minute, Int_t *second) const
|
---|
244 | {
|
---|
245 | *hour = fStarFieldRaH;
|
---|
246 | *minute = fStarFieldRaM;
|
---|
247 | *second = fStarFieldRaS;
|
---|
248 | }
|
---|
249 | // -------------------------------------------------------------------------
|
---|
250 | //
|
---|
251 | // GetStarFieldDec. Get DE coordinates of the starfield
|
---|
252 | //
|
---|
253 | void MMcRunHeader::GetStarFieldDec(Int_t *degree, Int_t *minute, Int_t *second) const
|
---|
254 | {
|
---|
255 | *degree = fStarFieldDeD;
|
---|
256 | *minute = fStarFieldDeM;
|
---|
257 | *second = fStarFieldDeS;
|
---|
258 | }
|
---|