source: trunk/MagicSoft/Simulation/Detector/Camera/creadparam.h@ 6347

Last change on this file since 6347 was 5600, checked in by MAGIC, 20 years ago
Added input card option pmt_jitter_ns to simulate the time jitter of PMTs. The input parameter is the sigma of a gaussian, which by default is sigma=0.25 ns. This jitter is applied to each phe- independently. We have not applied this to the NSB noise, since the arrival time of NSB photons is random and nothing would change.
File size: 10.9 KB
Line 
1//=//////////////////////////////////////////////////////////////////////
2//=
3//= creadparam
4//=
5//= @file creadparam.h
6//= @desc Header file
7//= @author J C Gonzalez
8//= @email gonzalez@mppmu.mpg.de
9//= @date Thu May 7 16:24:22 1998
10//=
11//=----------------------------------------------------------------------
12//=
13//= Created: Thu May 7 16:24:22 1998
14//= Author: Jose Carlos Gonzalez
15//= Purpose: Program for reflector simulation
16//= Notes: See files README for details
17//=
18//=----------------------------------------------------------------------
19//=
20//= $RCSfile: creadparam.h,v $
21//= $Revision: 1.31 $
22//= $Author: MAGIC $
23//= $Date: 2004-12-15 01:56:39 $
24//=
25//=//////////////////////////////////////////////////////////////////////
26
27#ifndef _creadparam_
28#define _creadparam_
29
30#ifndef _this_
31#define _this_ creadparam
32#endif
33
34//!@}
35
36//!@subsection Include files.
37
38//!@{
39
40#include <iostream>
41#include <fstream>
42#include <string>
43
44#include <stdlib.h>
45#include <stdio.h>
46#include <math.h>
47#include <float.h>
48
49#include "jcmacros.h"
50#include "jcdebug.h"
51
52#include "camera-v.h"
53
54#include "MTriggerDefine.h"
55#include "MFadcDefine.h"
56#include "Mdefine.h"
57
58//!@}
59
60//!@subsection Macro-definitions, and constants.
61
62//!@{
63
64// now we define the list PM_ITEM_LIST of possible items in
65// the parameters file. note that they are functions of
66// T(x). we will change T(x) to be the real item or the
67// string of this item when needed
68
69#define PM_ITEM_LIST /* LIST OF ITEMS IN THE PARAMETERS FILE */ \
70T(input_file), /* input file */ \
71T(starfield_file), /* starfield file */ \
72T(starfield_center),/* center of thestarfield FOV */ \
73T(data_file), /* data file */ \
74T(root_file), /* ROOT file */ \
75T(ct_num), /* number of CT, which is number of reflector files */ \
76T(ct_geom), /* class with the geometry of the CT camera */ \
77T(qe_file), /* file with the characteristics of the CT */ \
78T(elec_noise_off), /* de-activates ElecNoise for FADC and Trigger channels */\
79T(fadc_pedestal), /* Pedestal value for FADC */\
80T(fadc_noise), /* ElecNoise for FADC */\
81T(trig_noise), /* ElecNoise for trigger */\
82T(sfr_on), /* starfield rotation */ \
83T(nsb_on), /* activates NSB simulation */ \
84T(nsb_off), /* de-activates NSB simulation */ \
85T(nsb_mean), /* mean value of NSB contribution per pixel */ \
86T(nsb_directory), /* database for the NSB */ \
87T(nsb_dir_outer), /* database for the NSB */ \
88T(pixel_thres), /* value of q0 for trigger for a given pixel */ \
89T(secure_disc), /* Ac value to rise discrimnator threshold */ \
90T(seeds), /* seeds for random number generation */ \
91T(data_from_stdin), /* to read data from STDIN */ \
92T(skip), /* skip pathological showers */ \
93T(write_all_events),/* write to file .phe ALL images (even w.o. trigger)*/ \
94T(nowrite_McEvt), /* nowrite to file .root McEvt infos */ \
95T(write_McTrig), /* write to file .root McTrig infos */ \
96T(write_McFadc), /* write to file .root McFadc infos */ \
97T(nowrite_RawEvt), /* write to file .root RawData infos */ \
98T(select_energy), /* energy range to read: only for .phe files */ \
99T(trigger_scan), /* show signal in the diskriminator */ \
100T(fadc_scan), /* show fadc signal */ \
101T(trigger_loop), /* make loop over trigger conditions */ \
102T(trigger_prop), /* values of trigger properties */ \
103T(fadc_prop), /* values of fadc properties */ \
104T(fadc_outer), /* values of fadc properties for outer pixels */ \
105T(fadc_GHz), /* FADC sampling frequency (GHz) */ \
106T(trigger_single), /* trigger conditions */ \
107T(Trigger_Loop_Output_Only_Specialists), /* special steercard */ \
108T(high_to_low), /* High to Low gain ratio */ \
109T(sigma_xy_cm_spot), /* Additional gaussian spread of mirror spot on camera */ \
110T(sigma_x_cm_spot), /* Additional gaussian spread of mirror in x axis spot on camera */ \
111T(sigma_y_cm_spot), /* Additional gaussian spread of mirror in y axis spot on camera */ \
112T(misspoint_deg), /* Additional misspointing of Telescope direction*/ \
113T(trigger_delay), /* Delay in ns between beginning of FADC
114 time window and the trigger instant.*/ \
115T(telescope_location_cm), /* Telescope coordinates in Corsika system */ \
116T(mirror_fraction), /* Fraction of dish mirror really operative */ \
117T(calibration_run), /* option to produce a calibration file */ \
118T(gain_fluctuations_off), /* Switches off PMT gain fluctuations for signal */ \
119T(noise_gain_fluctuations_off), /* Switches off PMT gain fluctuations for NSB noise */ \
120T(pmt_jitter_ns), /* Sets the time jitter of each phe in the PMT */ \
121T(end_file) /* end of the parameters file */
122
123#define T(x) x // define T() as the name as it is
124
125enum ITEM_TYPE {
126 PM_ITEM_LIST
127};
128
129#undef T
130
131#define T(x) #x // define T() as the string of x
132
133const char *const ITEM_NAMES[] = {
134 PM_ITEM_LIST
135};
136
137#undef T
138
139#define LINE_MAX_LENGTH 400
140#define ITEM_MAX_LENGTH 40
141#define PATH_MAX_LENGTH 256
142
143//!@}
144
145//!@subsection Prototypes of functions.
146
147//!@{
148
149//++
150// prototypes
151//--
152
153void readparam(char * filename);
154char *get_input_filename(int i);
155char *get_starfield_filename(void);
156void get_starfield_center(int *rh,int *rm,int *rs,int *dd,int *dm,int *ds);
157int get_starfield_rotate(void);
158char *get_data_filename(void);
159char *get_root_filename(void);
160char *get_loop_filename(void);
161int get_ct_number(void);
162int get_ct_geometry(int ict);
163char *get_qe_filename(int ict=0);
164char *get_nsb_directory(void);
165char *get_nsb_directory_outer(void);
166int add_elec_noise(float *fadcinner, float* fadcouter,
167 float *digi, float *trig);
168int get_nsb(float *n, int *m);
169long int get_seeds(int n);
170void get_skip_showers( int *s );
171int get_nskip_showers( void );
172int get_data_from_stdin(void);
173int get_write_all_events(void);
174int get_write_McEvt(void);
175int get_write_McTrig(void);
176int get_write_McFadc(void);
177int get_write_RawEvt(void);
178int get_select_energy(float *le, float *ue);
179int get_FADC_Scan(void);
180int get_Trigger_Scan(void);
181void get_Trigger_properties(float *gl, float *ot, float *ra, float *rf);
182void get_FADC_properties(int *shape, float *ra, float *rf,int *shapeo, float *rao, float *rfo,
183 float *fadc_spns, int *fadc_slices);
184float get_FADC_pedestal(void);
185int get_Trigger_Loop(float *lt, float *ut, float *st, int *lm, int *um, int *lg, int *ug);
186void get_Trigger_Single(float **t, int *m, int *g);
187void get_threshold(float *t);
188void get_secure_threhold(float *ac, float *disc);
189int get_indi_thres_pixel(void);
190float get_High_to_Low(void);
191float get_sigma_xy_cm_spot(float *x, float *y);
192int get_misspointing(float *x, float *y);
193float get_trig_delay(void);
194float get_telescope_location_cm(int j, int icoor);
195float get_mirror_fraction(int j);
196float get_pmt_jitter_ns();
197
198int is_calibration_run();
199void get_calibration_properties(float *lambda, float *sigma_lambda, float *phot_per_pix,
200 float *sigma_time, int *nevents, int *selected_pixel);
201int apply_gain_fluctuations();
202int apply_noise_gain_fluctuations();
203//!@}
204
205//!@{
206
207#endif // ! _creadparam_
208
209//!@}
210
211//=------------------------------------------------------------
212//!@subsection Log of this file.
213
214//!@{
215
216/*
217 * $Log: not supported by cvs2svn $
218 * Revision 1.30 2004/11/17 11:36:07 moralejo
219 * *** empty log message ***
220 *
221 * Revision 1.29 2004/10/26 13:58:59 moralejo
222 *
223 * Added option gain_fluctuations_off
224 *
225 * Revision 1.28 2004/10/13 17:05:05 moralejo
226 * *** empty log message ***
227 *
228 * Revision 1.27 2004/10/12 13:32:02 moralejo
229 * *** empty log message ***
230 *
231 * Revision 1.26 2004/09/17 14:48:21 moralejo
232 *
233 * Changes to adapt headers to c++ style.
234 *
235 * Revision 1.25 2004/09/16 15:27:08 moralejo
236 * Updated in CVS after some time (see changes below)
237 *
238 *
239 * Revision 1.23 2004/01/30 09:55:33 blanch
240 * New "trigger_delay" and "sigma_xy_cm_spot" commands.
241 *
242 * Revision 1.22 2003/10/26 19:51:30 blanch
243 * Header file for creadparam.cxx version 1.25
244 *
245 * Revision 1.21 2003/10/17 19:40:47 blanch
246 * *** empty log message ***
247 *
248 * Revision 1.20 2003/09/23 17:38:59 blanch
249 * *** empty log message ***
250 *
251 * Revision 1.19 2003/09/23 16:52:43 blanch
252 * Header file for creadparam.cxx 1.22
253 *
254 * Revision 1.18 2003/09/15 10:23:31 blanch
255 * Hedaer file for creadparam.cxx 1.21
256 *
257 * Revision 1.17 2003/07/17 18:05:13 blanch
258 * Header file for creadparam.cxx 1.20
259 *
260 * Revision 1.16 2003/01/14 13:37:32 blanch
261 * Option to set a dc value to rise the discriminator threshold has been added.
262 *
263 * Revision 1.15 2002/07/16 16:20:59 blanch
264 * Modifications done for the camera.cxx version, where a first implementation
265 * of the Star Field Rotation has been introduced.
266 *
267 * Revision 1.14 2002/03/15 15:17:57 blanch
268 * Header file for creadparam.cxx v 1.16
269 *
270 * Revision 1.13 2002/03/04 17:20:50 blanch
271 * Header file for creadparam.cxx 1.15
272 *
273 * Revision 1.12 2002/02/28 15:08:01 blanch
274 * Header file for creadparam.cxx 1.14.
275 *
276 * Revision 1.11 2002/01/18 17:45:37 blanch
277 * Header file for creadparam.cxx 1.13.
278 *
279 * Revision 1.10 2001/11/13 17:06:01 blanch
280 * Header file for creadparam.cxx version 1.12
281 *
282 * Revision 1.9 2001/07/19 09:29:02 blanch
283 * Header file for creadparam.cxx vrsion 1.10
284 *
285 * Revision 1.8 2001/03/05 10:44:44 blanch
286 * Header file of creadparam.cxx version 1.9
287 *
288 * Revision 1.7 2001/02/23 10:56:10 magicsol
289 * Header file for creadparam.cxx version 1.8
290 *
291 * Revision 1.6 2001/01/15 12:38:00 magicsol
292 * *** empty log message ***
293 *
294 * Revision 1.5 2000/07/04 14:13:34 MagicSol
295 * *** empty log message ***
296 *
297 * Revision 1.4 2000/05/11 14:23:28 blanch
298 * Thi version owns to the version 1.4 of creadparam.cxx.
299 *
300 * Revision 1.3 2000/03/24 18:15:46 blanch
301 * *** empty log message ***
302 *
303 * Revision 1.2 2000/02/18 17:50:07 petry
304 * This version belongs to camera.cxx 1.5.
305 * It version has been put in the repository in order to be
306 * able to share the further development with others.
307 *
308 * If you need something working, wait or take an earlier one.
309 * See file README
310 *
311 * Revision 1.1.1.1 1999/11/05 11:59:34 harald
312 * This the starting point for CVS controlled further developments of the
313 * camera program. The program was originally written by Jose Carlos.
314 * But here you can find a "rootified" version to the program. This means
315 * that there is no hbook stuff in it now. Also the output of the
316 * program changed to the MagicRawDataFormat.
317 *
318 * The "rootification" was done by Dirk Petry and Harald Kornmayer.
319 *
320 * Revision 1.2 1999/10/22 15:01:29 petry
321 * version sent to H.K. and N.M. on Fri Oct 22 1999
322 *
323 * Revision 1.1.1.1 1999/10/21 16:35:11 petry
324 * first synthesised version
325 *
326 * Revision 1.7 1999/03/15 14:59:09 gonzalez
327 * camera-1_1
328 *
329 * Revision 1.6 1999/03/02 09:56:13 gonzalez
330 * *** empty log message ***
331 *
332 * Revision 1.5 1999/01/14 17:32:43 gonzalez
333 * Added to camera the STDIN input option (data_from_input)
334 *
335 */
336
337//!@}
338//=EOF
Note: See TracBrowser for help on using the repository browser.