1 | # ==========================================================================
|
---|
2 | #############################################################################
|
---|
3 | # ==========================================================================
|
---|
4 | # General
|
---|
5 | # ==========================================================================
|
---|
6 | #############################################################################
|
---|
7 | # ==========================================================================
|
---|
8 |
|
---|
9 | # -------------------------------------------------------------------------
|
---|
10 | # Use this if you want to setup the logging stream for the jobs
|
---|
11 | # (overwrites command line options)
|
---|
12 | # -------------------------------------------------------------------------
|
---|
13 | MLog.VerbosityLevel: 4
|
---|
14 | #MLog.DebugLevel: 1
|
---|
15 | #MLog.NoColors: yes
|
---|
16 |
|
---|
17 | # ==========================================================================
|
---|
18 | #############################################################################
|
---|
19 | # ==========================================================================
|
---|
20 | # Ceres
|
---|
21 | # ==========================================================================
|
---|
22 | #############################################################################
|
---|
23 | # ==========================================================================
|
---|
24 |
|
---|
25 | # -------------------------------------------------------------------------
|
---|
26 | # Use this if you want to write the MJStar output somewhere
|
---|
27 | # If you don't want it, it is written to the calibration output anyhow.
|
---|
28 | # -------------------------------------------------------------------------
|
---|
29 | #PathOut: .
|
---|
30 | #PathIn: .
|
---|
31 |
|
---|
32 |
|
---|
33 | # -------------------------------------------------------------------------
|
---|
34 | # Configure Eventloop
|
---|
35 | # -------------------------------------------------------------------------
|
---|
36 | #MaxEvents: 10000
|
---|
37 | #Overwrite: yes,no
|
---|
38 |
|
---|
39 |
|
---|
40 | # -------------------------------------------------------------------------
|
---|
41 | # Use this to setup binnings. For more details see: MBinning::ReadEnv
|
---|
42 | # -------------------------------------------------------------------------
|
---|
43 | # BinningEnergy.Raw: 100 1 100000 log
|
---|
44 | # BinningSize.Raw: 100 1 10000000 log
|
---|
45 | BinningImpact.Raw: 40 0 1000
|
---|
46 | # BinningHeight.Raw: 50 0 12
|
---|
47 | # BinningAz.Raw: 360 -360 360
|
---|
48 | # BinningZd.Raw: 70 0 70
|
---|
49 | # BinningViewCone.Raw: 155 0 31
|
---|
50 | BinningTrigPos.Raw: 300 -25 275
|
---|
51 | # BinningEvtWidth.Raw: 150 0 25
|
---|
52 | # BinningDist.Raw: 100 0 2.5
|
---|
53 | # BinningDistC.Raw: 100 0 2.5
|
---|
54 | # BinningThreshold.Raw: 35 0.9 90000 log
|
---|
55 | # BinningEnergyEst.Raw: 100 0.9 90000 log
|
---|
56 |
|
---|
57 |
|
---|
58 | # -------------------------------------------------------------------------
|
---|
59 | # Initialize random number generator (see MJob::InitRandomNumberGenerator)
|
---|
60 | # -------------------------------------------------------------------------
|
---|
61 | RandomNumberGenerator: TRandom3
|
---|
62 | #RandomNumberSeedValue: 0
|
---|
63 |
|
---|
64 |
|
---|
65 | # -------------------------------------------------------------------------
|
---|
66 | # Ceres general setup
|
---|
67 | # -------------------------------------------------------------------------
|
---|
68 |
|
---|
69 | # Switch off the camera "electronics"
|
---|
70 | #Camera: Off
|
---|
71 |
|
---|
72 | # Force the use of the "hardware" trigger for calibration data
|
---|
73 | #ForceTrigger: Yes
|
---|
74 |
|
---|
75 |
|
---|
76 | # -------------------------------------------------------------------------
|
---|
77 | # Some setup for the atmosphere. The default should be well suited.
|
---|
78 | # -------------------------------------------------------------------------
|
---|
79 | #MSimAtmosphere.FileAerosols: resmc/atmopshere-aerosols.txt
|
---|
80 | #MSimAtmosphere.FileOzone: resmc/atmopshere-ozone.txt
|
---|
81 |
|
---|
82 |
|
---|
83 | # -------------------------------------------------------------------------
|
---|
84 | # Here you can control the poiting of the telescope. To switch on
|
---|
85 | # off-target observations set a value for the distance !=0 [deg].
|
---|
86 | # For details see MSimPointingPos
|
---|
87 | # -------------------------------------------------------------------------
|
---|
88 | #MSimPointingPos.OffTargetDistance: 0.4
|
---|
89 | #MSimPointingPos.OffTargetPhi: -1
|
---|
90 |
|
---|
91 |
|
---|
92 | # -------------------------------------------------------------------------
|
---|
93 | # Setup the reflector and camera geometry
|
---|
94 | # -------------------------------------------------------------------------
|
---|
95 | Reflector.Constructor: MReflector
|
---|
96 |
|
---|
97 | # For the file definition see MReflector::ReadFile
|
---|
98 |
|
---|
99 | Reflector.FileName: /fhgfs/groups/app/mc_workshop/datacenter/setup/ceres/030/fact-reflector.txt
|
---|
100 | # Units mm: ~0.02°
|
---|
101 | #Reflector.SetSigmaPSF: 1.8
|
---|
102 | # Units mm: ~0.03°
|
---|
103 | #Reflector.SetSigmaPSF: 2.6
|
---|
104 | #Reflector.SetSigmaPSF: 1.66
|
---|
105 |
|
---|
106 | # new Value from 7.12.2012:
|
---|
107 | # 0.95 of the reflected light was inside a circle with a radius of 2.25mm. This corresponds to a SigmaPSF of 1.3 mm
|
---|
108 | #Reflector.SetSigmaPSF: 1.3
|
---|
109 |
|
---|
110 | # new Value from 29.05.2013:
|
---|
111 | # Set SigmaPSF to pixelsize
|
---|
112 | Reflector.SetSigmaPSF: 10.0
|
---|
113 |
|
---|
114 | # --- DWARF ---
|
---|
115 | # distance from mirror to camera [cm]
|
---|
116 | #focal_distance 489.73
|
---|
117 | # TRUE focal of the paraboloid (must correspond to that of the mirror positions defined below!)
|
---|
118 | #paraboloid_focal 488.87
|
---|
119 |
|
---|
120 | #MGeomCam.Constructor: MGeomCamDwarf(189.8, 9.5, 4.8887);
|
---|
121 | MGeomCam.Constructor: MGeomCamFACT();
|
---|
122 |
|
---|
123 | #MSimBundlePhotons.FileName: resmc/dwarf-fact.txt
|
---|
124 |
|
---|
125 | # Set the APD type (1: 30x30 <default>, 2: 60x60, 3:60x60(ct=15%))
|
---|
126 | MSimAPD.Type: 4
|
---|
127 | MSimExcessNoise.ExcessNoise: 0.12
|
---|
128 |
|
---|
129 | # -------------------------------------------------------------------------
|
---|
130 | # Setup the absorption, conversion efficiency and angular acceptance
|
---|
131 | # -------------------------------------------------------------------------
|
---|
132 | MirrorReflectivity.FileName: /fhgfs/groups/app/mc_workshop/datacenter/setup/ceres/030/MirrorReflectivity_Lustermann_FACT_bearbeitet.txt
|
---|
133 | PhotonDetectionEfficiency.FileName: /fhgfs/groups/app/mc_workshop/datacenter/setup/ceres/030/fact-pde-1.0V.txt
|
---|
134 | ConesAngularAcceptance.FileName: /fhgfs/groups/app/mc_workshop/datacenter/setup/ceres/030/fact-cones-angular-acceptance.txt
|
---|
135 | ConesTransmission.FileName: /fhgfs/groups/app/mc_workshop/datacenter/setup/ceres/030/Transmittance_1439Cones_FACT_bearbeitet.txt
|
---|
136 |
|
---|
137 | # -------------------------------------------------------------------------
|
---|
138 | # Setup what in MMCS would be called "additional spot size"
|
---|
139 | # -------------------------------------------------------------------------
|
---|
140 | #MSimPSF.Sigma: -1
|
---|
141 | #MSimReflector.DetectorMargin: 0
|
---|
142 |
|
---|
143 | # -------------------------------------------------------------------------
|
---|
144 | # Setup the dark counts (FrequencyFixed) and the NSB noise per cm^2
|
---|
145 | # -------------------------------------------------------------------------
|
---|
146 | # Dark Counts per APD: ~4MHz
|
---|
147 | #MSimRandomPhotons.FrequencyFixed: 0.004
|
---|
148 | # NSB photon rate per cm^2 ~40MHz (folded with the cones' angular
|
---|
149 | # acceptance and the wavelength acceptance of the camera (window, apd, etc)
|
---|
150 | # 0.040 1/ns/cm^2 NSB-rate:
|
---|
151 | MSimRandomPhotons.FileNameNSB:
|
---|
152 | MSimRandomPhotons.FrequencyNSB: 0.0
|
---|
153 |
|
---|
154 | # FIXME: With a class describing the cones we could give NSB as
|
---|
155 | # per sr and cm^2
|
---|
156 |
|
---|
157 | # -------------------------------------------------------------------------
|
---|
158 | # Setup the trigger
|
---|
159 | # -------------------------------------------------------------------------
|
---|
160 | # This line could be omitted but then the discriminator would be
|
---|
161 | # evaluated for all pixels not just for the pixels which are
|
---|
162 | # later "connected" in the trigger (used in the coincidence map)
|
---|
163 | # MSimTrigger.FileNameRouteAC: resmc/magic-trigger-1NN.txt
|
---|
164 |
|
---|
165 | MSimTrigger.FileNameRouteAC: /fhgfs/groups/app/mc_workshop/datacenter/setup/ceres/030/fact-trigger-sum.txt
|
---|
166 | #MSimTrigger.DiscriminatorThreshold: 13416
|
---|
167 |
|
---|
168 | # new Value from 29.01.2013: DiscriminatorThreshold is 9*offset(3280)+21*gain(328)=36408
|
---|
169 | MSimTrigger.DiscriminatorThreshold: 36408
|
---|
170 |
|
---|
171 | # new Value from 28.08.2013: DiscriminatorThreshold is 9*offset(3720)+21*gain(372)=41292
|
---|
172 | MSimTrigger.DiscriminatorThreshold: 41292
|
---|
173 |
|
---|
174 | # HEGRA like single pixel above threshold trigger
|
---|
175 | #MSimTrigger.DiscriminatorThreshold: 10.00
|
---|
176 | # Is this an option for a single pixel trigger?
|
---|
177 | #MSimTrigger.DigitalSignalLength: -1
|
---|
178 | #MSimTrigger.CoincidenceTime: 3
|
---|
179 |
|
---|
180 | # Magic like 2NN trigger (inner 396 pixel)
|
---|
181 | #MSimTrigger.FileNameCoincidenceMap: resmc/magic-trigger-2NN.txt
|
---|
182 | #MSimTrigger.DiscriminatorThreshold: 7.00
|
---|
183 | #MSimTrigger.DigitalSignalLength: 8
|
---|
184 | #MSimTrigger.CoincidenceTime: 0.5
|
---|
185 |
|
---|
186 | # Magic like 3NN trigger (inner 396 pixel)
|
---|
187 | #MSimTrigger.FileNameCoincidenceMap: resmc/magic-trigger-3NN.txt
|
---|
188 | #MSimTrigger.DiscriminatorThreshold: 6.0
|
---|
189 | #MSimTrigger.DigitalSignalLength: 8
|
---|
190 | #MSimTrigger.CoincidenceTime: 0.5
|
---|
191 |
|
---|
192 | # Magic like 4NN trigger (inner 396 pixel)
|
---|
193 | #MSimTrigger.FileNameCoincidenceMap: resmc/magic-trigger-4NN.txt
|
---|
194 | #MSimTrigger.DiscriminatorThreshold: 5.0
|
---|
195 | #MSimTrigger.DigitalSignalLength: 8
|
---|
196 | #MSimTrigger.CoincidenceTime: 0.5
|
---|
197 |
|
---|
198 | #MSimTrigger.FileNameRouteAC: resmc/fact-7map.txt
|
---|
199 | #MSimTrigger.DiscriminatorThreshold: 14040
|
---|
200 | #MSimTrigger.DigitalSignalLength: 8
|
---|
201 | #MSimTrigger.CoincidenceTime: 0.5
|
---|
202 |
|
---|
203 | # Every Pixel(!) should see the same signal independant of its size
|
---|
204 | MSimCalibrationSignal.NumPhotons: 24
|
---|
205 | MSimCalibrationSignal.NumEvents: 1000
|
---|
206 |
|
---|
207 | IntendedPulsePos.Val: 30
|
---|
208 |
|
---|
209 | #PulseShape.Function.Name: exp(-(x/2)^2/2)
|
---|
210 | #PulseShape.FileName: /fhgfs/groups/app/mc_workshop/datacenter/setup/ceres/030/dwarf-pulse.txt
|
---|
211 |
|
---|
212 | # -------------------------------------------------------------------------
|
---|
213 | # Description how a pedestal/calibration signal is produced
|
---|
214 | # -------------------------------------------------------------------------
|
---|
215 |
|
---|
216 | #MSimCalibrationSignal.NumEvents: 1000
|
---|
217 | #MSimCalibrationSignal.NumPhotons; 5
|
---|
218 | #MSimCalibrationSignal.TimeJitter; 1.0
|
---|
219 |
|
---|
220 |
|
---|
221 | # -------------------------------------------------------------------------
|
---|
222 | # Setup the FADC
|
---|
223 | # -------------------------------------------------------------------------
|
---|
224 |
|
---|
225 | MRawRunHeader.SamplingFrequency: 2000
|
---|
226 | MRawRunHeader.NumSamples: 1024
|
---|
227 | MRawRunHeader.NumBytesPerSample: 2
|
---|
228 | MRawRunHeader.FadcResolution: 12
|
---|
229 |
|
---|
230 | # Old values
|
---|
231 | # Baseline at 15 like in MAGIC
|
---|
232 | # 1.5 bit noise for a gain of 64
|
---|
233 | # This allows a maximum of 200phe/pix
|
---|
234 | #MSimCamera.DefaultOffset: 960
|
---|
235 | #MSimCamera.DefaultNoise: 96
|
---|
236 | #MSimCamera.DefaultGain: 256
|
---|
237 |
|
---|
238 | # New values
|
---|
239 | #MSimCamera.DefaultOffset: 1000
|
---|
240 | #MSimCamera.DefaultNoise: 30
|
---|
241 | #MSimCamera.DefaultGain: 190
|
---|
242 |
|
---|
243 | # New values from 7.12.2012
|
---|
244 | # gain of one pe is ~10mV, corresponds to 328 DAC (16 bit with a dynamic range of 2 V)
|
---|
245 | # offset is 10 pe (->3280 DAC)
|
---|
246 | # electronic noise is about 0.15 pe (->50 DAC)
|
---|
247 | #MSimCamera.DefaultOffset: 3280
|
---|
248 | #MSimCamera.DefaultNoise: 50
|
---|
249 | #MSimCamera.DefaultGain: 328
|
---|
250 |
|
---|
251 | # New values from 28.08.2013
|
---|
252 | # Gain (integral over 30 slices) of one pe is 123.5 mV ns in real data.
|
---|
253 | # Gain(integral over 30 slices over the PulseShape is 10.8757 1ns.
|
---|
254 | # The scaling factor (DefaultGain) should be 123.5 / 10.8757 mV = 11,3556 mV
|
---|
255 | # this corresponds to 372 ADC (16 bit with a dynamic range of 2 V)
|
---|
256 | # offset is 10 pe (->3720 ADC)
|
---|
257 | # electronic noise is about 1.4 mV (->45 ADC)
|
---|
258 | MSimCamera.DefaultOffset: 3720
|
---|
259 | MSimCamera.DefaultNoise: 45
|
---|
260 | MSimCamera.DefaultGain: 358
|
---|
261 |
|
---|
262 | # Value for the fudgefactor in the calculation of the accoupling:
|
---|
263 | MSimCamera.fACFudgeFactor: 0.3367
|
---|
264 |
|
---|
265 | #MSimReadout.fConversionFactor: 1
|
---|
266 |
|
---|
267 | # The number of sampling points is almost irrelevant because they
|
---|
268 | # are equidistant, i.e. calculated and no search is necessary.
|
---|
269 | # Nevertheless, you must make sure that there are enough points
|
---|
270 | # to sample the function accuratly enough.
|
---|
271 | # Attention: x in the function is given in slices, so if you change the sampling
|
---|
272 | # frequency you have to change also this function
|
---|
273 | PulseShape.Function.Name: (1.626*(1-exp(-0.3803*0.5*x))*exp(-0.0649*0.5*x))
|
---|
274 | PulseShape.Function.Npx: 300
|
---|
275 | PulseShape.Function.Xmin: 0
|
---|
276 | PulseShape.Function.Xmax: 300
|
---|
277 |
|
---|
278 | # -------------------------------------------------------------------------
|
---|
279 | # Setup an image cleaning on the pure gamma signal (without noise)
|
---|
280 | # -------------------------------------------------------------------------
|
---|
281 |
|
---|
282 | #MImgCleanStd.CleanLevel1: 6.0
|
---|
283 | #MImgCleanStd.CleanLevel2: 6.0
|
---|
284 |
|
---|
285 | # -------------------------------------------------------------------------
|
---|
286 | # This is a cut executed after the calculation of the image parameters
|
---|
287 | # -------------------------------------------------------------------------
|
---|
288 | Cut.Inverted: Yes
|
---|
289 |
|
---|
290 | #Cut.Condition: MNewImagePar.fNumUsedPixels>3 && MHillas.fSize>6
|
---|
291 | Cut.Condition: MHillas.fSize>10
|
---|
292 | #&& MHillas.fSize>30
|
---|
293 | #MSimTrigger.SimulateElectronics: Off
|
---|
294 |
|
---|
295 | # Does not trigger anyway
|
---|
296 | ContEmpty3.Condition: MPhotonEvent.GetNumPhotons<10
|
---|