1 | #ifndef __RFL_HEADER__
|
---|
2 | #define __RFL_HEADER__
|
---|
3 |
|
---|
4 | /* RUNH from cerfile. See CORSIKA manual for explanations */
|
---|
5 | typedef struct
|
---|
6 | { char RUNH[4];
|
---|
7 | float RunNumber;
|
---|
8 | float date;
|
---|
9 | float Corsika_version;
|
---|
10 | float NumObsLev;
|
---|
11 | float HeightLev[10];
|
---|
12 | float SlopeSpec; /* Slope of primaries' energy spectrum */
|
---|
13 | float ELowLim;
|
---|
14 | float EUppLim; /* Limits of energy range for generation */
|
---|
15 | float EGS4_flag;
|
---|
16 | float NKG_flag;
|
---|
17 | float Ecutoffh;
|
---|
18 | float Ecutoffm;
|
---|
19 | float Ecutoffe;
|
---|
20 | float Ecutoffg;
|
---|
21 | /* Physical constants and interaction flags (see CORSIKA manual): */
|
---|
22 | float C[50];
|
---|
23 | float dummy1[20]; /* not used */
|
---|
24 | float CKA[40];
|
---|
25 | float CETA[5];
|
---|
26 | float CSTRBA[11];
|
---|
27 | float dummy2[104]; /* not used */
|
---|
28 | float AATM[5];
|
---|
29 | float BATM[5];
|
---|
30 | float CATM[5];
|
---|
31 | float NFL[4];
|
---|
32 | } CerRunHeader;
|
---|
33 |
|
---|
34 |
|
---|
35 | /* EVTH from cerfile. See CORSIKA manual for explanations */
|
---|
36 | typedef struct
|
---|
37 | { char EVTH[4];
|
---|
38 | float EvtNumber;
|
---|
39 | float PrimaryID;
|
---|
40 | float Etotal;
|
---|
41 | float Thick0;
|
---|
42 | float FirstTarget;
|
---|
43 | float zFirstInt;
|
---|
44 | float p[3];
|
---|
45 | float Theta;
|
---|
46 | float Phi;
|
---|
47 |
|
---|
48 | float NumRndSeq;
|
---|
49 | float RndData[10][3];
|
---|
50 |
|
---|
51 | float RunNumber;
|
---|
52 | float DateRun;
|
---|
53 | float Corsika_version;
|
---|
54 |
|
---|
55 | float NumObsLev;
|
---|
56 | float HeightLev[10];
|
---|
57 |
|
---|
58 | float SlopeSpec;
|
---|
59 | float ELowLim;
|
---|
60 | float EUppLim;
|
---|
61 |
|
---|
62 | float Ecutoffh;
|
---|
63 | float Ecutoffm;
|
---|
64 | float Ecutoffe;
|
---|
65 | float Ecutoffg;
|
---|
66 |
|
---|
67 | float NFLAIN;
|
---|
68 | float NFLDIF;
|
---|
69 | float NFLPI0;
|
---|
70 | float NFLPIF;
|
---|
71 | float NFLCHE;
|
---|
72 | float NFRAGM;
|
---|
73 |
|
---|
74 | float Bx;
|
---|
75 | float By;
|
---|
76 |
|
---|
77 | float EGS4yn;
|
---|
78 | float NKGyn;
|
---|
79 | float GHEISHAyn;
|
---|
80 | float VENUSyn;
|
---|
81 | float CERENKOVyn;
|
---|
82 | float NEUTRINOyn;
|
---|
83 | float HORIZONTyn;
|
---|
84 | float COMPUTER;
|
---|
85 |
|
---|
86 | float ThetaMin;
|
---|
87 | float ThetaMax;
|
---|
88 | float PhiMin;
|
---|
89 | float PhiMax;
|
---|
90 |
|
---|
91 | float CBunchSize;
|
---|
92 | float CDetInX,CDetInY;
|
---|
93 | float CSpacInX,CSpacInY;
|
---|
94 | float CLenInX,CLenInY;
|
---|
95 | float COutput;
|
---|
96 |
|
---|
97 | float AngleNorthX;
|
---|
98 | float MuonInfo;
|
---|
99 |
|
---|
100 | float StepLength;
|
---|
101 | float CWaveLower;
|
---|
102 | float CWaveUpper;
|
---|
103 | float Multipl;
|
---|
104 | float CorePos[2][20];
|
---|
105 | float SIBYLL[2];
|
---|
106 | float QGSJET[2];
|
---|
107 | float DPMJET[2];
|
---|
108 | float VENUS_cross;
|
---|
109 | float mu_mult_scat;
|
---|
110 | float NKG_range;
|
---|
111 | float EFRCTHN[2];
|
---|
112 | float WMAX[2];
|
---|
113 | float rthin_rmax;
|
---|
114 |
|
---|
115 | float viewcone_angles[2]; /* (degrees) Inner and outer angles in
|
---|
116 | * Corsika's VIEWCONE option
|
---|
117 | */
|
---|
118 | float dmmy2[119];
|
---|
119 | } CerEventHeader;
|
---|
120 |
|
---|
121 | /* EVTH from rflfile. We include first all the information from the
|
---|
122 | Corsika event header (except empty bytes) */
|
---|
123 |
|
---|
124 | typedef struct
|
---|
125 | { char EVTH[4];
|
---|
126 | float EvtNumber;
|
---|
127 | float PrimaryID;
|
---|
128 | float Etotal;
|
---|
129 | float Thick0;
|
---|
130 | float FirstTarget;
|
---|
131 | float zFirstInt;
|
---|
132 | float p[3];
|
---|
133 | float Theta;
|
---|
134 | float Phi;
|
---|
135 |
|
---|
136 | float NumRndSeq;
|
---|
137 | float RndData[10][3];
|
---|
138 |
|
---|
139 | float RunNumber;
|
---|
140 | float DateRun;
|
---|
141 | float Corsika_version;
|
---|
142 |
|
---|
143 | float NumObsLev;
|
---|
144 | float HeightLev[10];
|
---|
145 |
|
---|
146 | float SlopeSpec;
|
---|
147 | float ELowLim;
|
---|
148 | float EUppLim;
|
---|
149 |
|
---|
150 | float Ecutoffh;
|
---|
151 | float Ecutoffm;
|
---|
152 | float Ecutoffe;
|
---|
153 | float Ecutoffg;
|
---|
154 |
|
---|
155 | float NFLAIN;
|
---|
156 | float NFLDIF;
|
---|
157 | float NFLPI0;
|
---|
158 | float NFLPIF;
|
---|
159 | float NFLCHE;
|
---|
160 | float NFRAGM;
|
---|
161 |
|
---|
162 | float Bx;
|
---|
163 | float By;
|
---|
164 |
|
---|
165 | float EGS4yn;
|
---|
166 | float NKGyn;
|
---|
167 | float GHEISHAyn;
|
---|
168 | float VENUSyn;
|
---|
169 | float CERENKOVyn;
|
---|
170 | float NEUTRINOyn;
|
---|
171 | float HORIZONTyn;
|
---|
172 | float COMPUTER;
|
---|
173 |
|
---|
174 | float ThetaMin;
|
---|
175 | float ThetaMax;
|
---|
176 | float PhiMin;
|
---|
177 | float PhiMax;
|
---|
178 |
|
---|
179 | float CBunchSize;
|
---|
180 | float CDetInX,CDetInY;
|
---|
181 | float CSpacInX,CSpacInY;
|
---|
182 | float CLenInX,CLenInY;
|
---|
183 | float COutput;
|
---|
184 |
|
---|
185 | float AngleNorthX;
|
---|
186 | float MuonInfo;
|
---|
187 |
|
---|
188 | float StepLength;
|
---|
189 | float CWaveLower;
|
---|
190 | float CWaveUpper;
|
---|
191 | float Multipl;
|
---|
192 | float CorePos[2][20];
|
---|
193 | float SIBYLL[2];
|
---|
194 | float QGSJET[2];
|
---|
195 | float DPMJET[2];
|
---|
196 | float VENUS_cross;
|
---|
197 | float mu_mult_scat;
|
---|
198 | float NKG_range;
|
---|
199 | float EFRCTHN[2];
|
---|
200 | float WMAX[2];
|
---|
201 | float rthin_rmax;
|
---|
202 |
|
---|
203 | float viewcone_angles[2];
|
---|
204 | /* (degrees) Inner and outer angles in
|
---|
205 | * Corsika's VIEWCONE option. This is
|
---|
206 | * only possible with Corsika>6 versions. In
|
---|
207 | * that case, PhiMin=PhiMax and
|
---|
208 | * ThetaMin=ThetaMax (also in this header)
|
---|
209 | * indicate the axis of this cone.
|
---|
210 | */
|
---|
211 |
|
---|
212 | /* ^^^ Up to here, the info from the CORSIKA event header. */
|
---|
213 |
|
---|
214 | /* Telescope orientation: */
|
---|
215 | float telescopePhi; /* rad */
|
---|
216 | float telescopeTheta; /* rad */
|
---|
217 |
|
---|
218 | /* Time of first and last photon: */
|
---|
219 | float TimeFirst;
|
---|
220 | float TimeLast;
|
---|
221 |
|
---|
222 | /* 6 parameters and chi2 of the NKG fit to the longitudinal
|
---|
223 | * particle distribution (see CORSIKA manual for explanation):
|
---|
224 | */
|
---|
225 | float longi_Nmax;
|
---|
226 | float longi_t0;
|
---|
227 | float longi_tmax;
|
---|
228 | float longi_a;
|
---|
229 | float longi_b;
|
---|
230 | float longi_c;
|
---|
231 | float longi_chi2;
|
---|
232 |
|
---|
233 | /* Now the photon statistics of the event: */
|
---|
234 | float CORSIKAPhs; /* Original photons written by Corsika */
|
---|
235 | float AtmAbsPhs; /* Photons absorbed by the atmosphere */
|
---|
236 | float MirrAbsPhs; /* Photons absorbed by the mirror */
|
---|
237 | float OutOfMirrPhs; /* Photons outside the mirror */
|
---|
238 | float BlackSpotPhs; /* Photons lost in the "black spot" */
|
---|
239 | float OutOfChamPhs; /* Photons outside the camera */
|
---|
240 | float CPhotons; /* Photons reaching the camera */
|
---|
241 |
|
---|
242 | float dummy[10]; /* not used */
|
---|
243 | } RflEventHeader;
|
---|
244 |
|
---|
245 |
|
---|
246 | typedef struct /* EVTE from cerfile */
|
---|
247 | {
|
---|
248 | float dummy1[255];
|
---|
249 | float longi_Nmax;
|
---|
250 | float longi_t0;
|
---|
251 | float longi_tmax;
|
---|
252 | float longi_a;
|
---|
253 | float longi_b;
|
---|
254 | float longi_c;
|
---|
255 | float longi_chi2;
|
---|
256 | float dummy2[11];
|
---|
257 | } Event_end;
|
---|
258 |
|
---|
259 | void TranslateHeader(RflEventHeader *r, CerEventHeader *c);
|
---|
260 |
|
---|
261 | #endif
|
---|