source: tags/Mars-V0.8.6/mmc/MMcRunHeader.cxx

Last change on this file was 3894, checked in by moralejo, 21 years ago
*** empty log message ***
  • Property svn:executable set to *
File size: 7.3 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// 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
76ClassImp(MMcRunHeader);
77
78using namespace std;
79
80// --------------------------------------------------------------------------
81//
82// Default constructor.
83//
84//
85MMcRunHeader::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//
142MMcRunHeader::~MMcRunHeader()
143{
144
145}
146
147// -------------------------------------------------------------------------
148//
149// Fill. Put data in the container
150//
151void 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//
243void 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//
253void MMcRunHeader::GetStarFieldDec(Int_t *degree, Int_t *minute, Int_t *second) const
254{
255 *degree = fStarFieldDeD;
256 *minute = fStarFieldDeM;
257 *second = fStarFieldDeS;
258}
Note: See TracBrowser for help on using the repository browser.