| 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: 1
|
|---|
| 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
|
|---|
