Changeset 6710
- Timestamp:
- 03/02/05 20:49:19 (20 years ago)
- Location:
- trunk/MagicSoft/Simulation/Detector
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Simulation/Detector/Camera/creadparam.cxx
r6560 r6710 19 19 //= 20 20 //= $RCSfile: creadparam.cxx,v $ 21 //= $Revision: 1.3 7$21 //= $Revision: 1.38 $ 22 22 //= $Author: moralejo $ 23 //= $Date: 2005-0 2-17 09:15:28$23 //= $Date: 2005-03-02 20:49:19 $ 24 24 //= 25 25 //=////////////////////////////////////////////////////////////////////// … … 103 103 static float misp_y = 0.; // Mispointing in y 104 104 105 static float trig_delay = 25.; // Delay in ns between beginning of FADC105 static float trig_delay = 19.; // Delay in ns between beginning of FADC 106 106 // time window and the trigger instant. 107 107 … … 1445 1445 // 1446 1446 // $Log: not supported by cvs2svn $ 1447 // Revision 1.37 2005/02/17 09:15:28 moralejo 1448 // 1449 // Set as default option that of writing all event headers to output file, 1450 // not only those of the triggered events. To disable it, set the input card 1451 // flag "no_write_all_event_headers". 1452 // 1453 // Changed such that output images for events below the minimum number of 1454 // photoelectrons nphe2NSB required to simulate the noise (NSB & electronic) 1455 // will be empty. This will avoid the problem of these events being processed, 1456 // without any noise, later in the chain. Although those images are not in the 1457 // output, one can still check in the headers (MMcTrig) how many such events 1458 // with less than nphe2NSB photoelectrons would have triggered. 1459 // 1447 1460 // Revision 1.36 2005/02/10 19:28:10 moralejo 1448 1461 // -
trunk/MagicSoft/Simulation/Detector/Camera/input.card
r6588 r6710 5 5 ct_geom 1 6 6 # Quantum efficiency file: 7 qe_file 0 / users/emc/moralejo/MagicSoft/Simulation/Detector/Data/qe-emi-coat.RFL.dat7 qe_file 0 /home/magic/MagicSoft/Simulation/Detector/Data/qe-emi-coat.RFL.dat 8 8 # Input file (one per telescope): 9 9 # input_file 0 /data1/magic/reflex/Gamma_zbin0_0_7_1000to1009_w0.rfl 10 10 # Perform calibration run: lambda sigma_lambda phot_per_pixel time_fwhm n_events [selected_pixel] 11 11 # The values below correspond roughly to 10 LED UV 12 #calibration_run 375. 12. 120. 2.5 300013 calibration_run 375. 12. 120. 2.5 1012 #calibration_run 375. 12. 120. 2.5 5000 13 calibration_run 375. 12. 120. 0.01 1 14 14 # line below shows how to create a pedestal run (= cal. with 0-photon pulses): 15 15 #calibration_run 0. 0. 0. 0. 1000 … … 26 26 # L1 Trigger condition: CT number, threshold (mV), multiplicity and topology: 27 27 trigger_single 0 4 4 2 28 # To shift the pulses in the FADC window, modify the trigger delay (ns): 29 #trigger_delay 18. 28 30 # Correction to overall light collection efficiency: CT# fraction 29 31 mirror_fraction 0 0.73 30 32 # Switch on NSB: 31 nsb_on32 #nsb_off33 #nsb_on 34 nsb_off 33 35 # Number of photons from the diffuse NSB (nphe / ns 0.1*0.1 deg^2 239 m^2) and 34 36 # minimum number of phe from shower required to simulate NSB: 35 37 nsb_mean 0.183 10 36 38 # Starfield (see Starfieldadder program) 37 # starfield_file / users/emc/moralejo/MagicSoft/Simulation/Detector/Starfield/starfield.rfl39 # starfield_file /home/magic/MagicSoft/Simulation/Detector/Starfield/starfield.rfl 38 40 # Electronic noise in FADC (sigma in ADC counts): Inner pixels, outer pixels, digital noise: 39 fadc_noise 1.3 2.4 1.40 #elec_noise_off41 #fadc_noise 1.3 2.4 1. 42 elec_noise_off 41 43 # Mean pedestal per slice (ADC counts): 42 44 fadc_pedestal 10. … … 46 48 seeds 66767 45069 47 49 # Directory where NSB database can be found for inner and outer pixels: 48 nsb_directory / users/emc/moralejo/MagicSoft/Simulation/Detector/StarLight/inner/49 nsb_dir_outer / users/emc/moralejo/MagicSoft/Simulation/Detector/StarLight/outer/50 nsb_directory /home/magic/MagicSoft/Simulation/Detector/StarLight/inner/ 51 nsb_dir_outer /home/magic/MagicSoft/Simulation/Detector/StarLight/outer/ 50 52 # 51 53 # FADC properties: shape of single phe response (1 means realistic one, from -
trunk/MagicSoft/Simulation/Detector/include-MFadc/MFadc.cxx
r6695 r6710 227 227 fadc_time_offset = trigger_delay - p2 / FADC_SLICES_PER_NSEC; // ns 228 228 229 229 230 for (i=0; i< fResponseSlicesFadc ; i++ ) 230 231 { … … 241 242 p3+p4*exp(-p1*(exp(-p1*zed_slices)+ 242 243 p5*zed_slices))+p6*d); 244 response_sum_inner += sing_resp[i]; 245 246 247 // Now the low gain: 243 248 244 249 zed_slices = x * FADC_SLICES_PER_NSEC - p2_LG; 245 250 d = (zed_slices>0)? 0.5 : -0.5; 251 246 252 sing_resp_lowgain[i] = (Float_t) (p0_LG*exp(-p1_LG*(exp(-p1_LG*zed_slices)+zed_slices))+ 247 p3+p4*exp(-p1_LG*(exp(-p1_LG*zed_slices)+ 248 p5*zed_slices))+p6*d); 249 250 response_sum_inner += sing_resp[i]; 253 p3_LG+p4_LG*exp(-p1_LG*(exp(-p1_LG*zed_slices)+ 254 p5_LG*zed_slices))+p6_LG*d); 251 255 response_sum_inner_LG += sing_resp_lowgain[i]; 256 252 257 } 253 258 … … 268 273 sigma = fwhm_resp_outer / 2.35 ; 269 274 x0 = 3*sigma ; 270 fadc_time_offset = trigger_delay-x0; // ns271 275 272 276 for (i = 0; i < fResponseSlicesFadc ; i++ ) … … 309 313 p6_LG = fPulseParametersLG[6]; 310 314 311 // Now define the time before trigger to read FADC signal when it312 // has to be written. Here FADC_SLICES_PER_NSEC (=0.3) is the value313 // for the 300 MHz MAGIC FADCs and must NOT be changed, even if you314 // use a faster sampling in the simulation (through the input card315 // command "fadc_GHz"), because this is just a conversion of units. The316 // parameters of the "pulpo" pulse shape were obtained with the 300 MHz317 // FADC and so we convert the time parameter to units of 3.3 ns slices318 // just to use the provided parametrization, and no matter what sampling319 // frequency we are simulating!320 321 fadc_time_offset = trigger_delay - p2 / FADC_SLICES_PER_NSEC; // ns322 323 315 for (i=0; i< fResponseSlicesFadc ; i++ ) 324 316 { … … 336 328 p4*exp(-p1*(exp(-p1*zed_slices)+ 337 329 p5*zed_slices))+p6*d); 330 response_sum_outer += sing_resp_outer[i]; 331 332 333 // Now the low gain: 338 334 339 335 zed_slices = x * FADC_SLICES_PER_NSEC - p2_LG; 340 336 d = (zed_slices>0)? 0.5 : -0.5; 337 341 338 sing_resp_outer_lowgain[i] = (Float_t) (p0_LG*exp(-p1_LG*(exp(-p1_LG*zed_slices)+zed_slices))+ 342 p3+p4*exp(-p1_LG*(exp(-p1_LG*zed_slices)+ 343 p5*zed_slices))+p6*d); 344 345 response_sum_outer += sing_resp_outer[i]; 339 p3_LG+p4_LG*exp(-p1_LG*(exp(-p1_LG*zed_slices)+ 340 p5_LG*zed_slices))+p6_LG*d); 346 341 response_sum_outer_LG += sing_resp_outer_lowgain[i]; 347 342 } … … 492 487 // We take the pulse height in the middle of FADC slices, we start in the 493 488 // first such point after the time "time" (=ichan in response bins). Each 494 // FADC slice corresponds to SUBBINS response bins (SUBBINS=5 by default).489 // FADC slice corresponds to SUBBINS response bins. 495 490 496 491 Int_t first_i = Int_t(SUBBINS/2) - ichan%(Int_t)SUBBINS;
Note:
See TracChangeset
for help on using the changeset viewer.