| 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: 2
|
|---|
| 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 |
|
|---|
| 111 |
|
|---|
| 112 | # --- DWARF ---
|
|---|
| 113 | # distance from mirror to camera [cm]
|
|---|
| 114 | #focal_distance 489.73
|
|---|
| 115 | # TRUE focal of the paraboloid (must correspond to that of the mirror positions defined below!)
|
|---|
| 116 | #paraboloid_focal 488.87
|
|---|
| 117 |
|
|---|
| 118 | #MGeomCam.Constructor: MGeomCamDwarf(189.8, 9.5, 4.8887);
|
|---|
| 119 | MGeomCam.Constructor: MGeomCamFACT();
|
|---|
| 120 |
|
|---|
| 121 | #MSimBundlePhotons.FileName: resmc/dwarf-fact.txt
|
|---|
| 122 |
|
|---|
| 123 | # Set the APD type (1: 30x30 <default>, 2: 60x60, 3:60x60(ct=15%))
|
|---|
| 124 | MSimAPD.Type: 4
|
|---|
| 125 | MSimExcessNoise.ExcessNoise: 0.12
|
|---|
| 126 |
|
|---|
| 127 | # -------------------------------------------------------------------------
|
|---|
| 128 | # Setup the absorption, conversion efficiency and angular acceptance
|
|---|
| 129 | # -------------------------------------------------------------------------
|
|---|
| 130 | MirrorReflectivity.FileName: /fhgfs/groups/app/mc_workshop/datacenter/setup/ceres/030/MirrorReflectivity_Lustermann_FACT_bearbeitet.txt
|
|---|
| 131 | PhotonDetectionEfficiency.FileName: /fhgfs/groups/app/mc_workshop/datacenter/setup/ceres/030/fact-pde-1.4V.txt
|
|---|
| 132 | ConesAngularAcceptance.FileName: /fhgfs/groups/app/mc_workshop/datacenter/setup/ceres/030/fact-cones-angular-acceptance.txt
|
|---|
| 133 | ConesTransmission.FileName: /fhgfs/groups/app/mc_workshop/datacenter/setup/ceres/030/Transmittance_1439Cones_FACT_bearbeitet.txt
|
|---|
| 134 | # -------------------------------------------------------------------------
|
|---|
| 135 | # Setup what in MMCS would be called "additional spot size"
|
|---|
| 136 | # -------------------------------------------------------------------------
|
|---|
| 137 | #MSimPSF.Sigma: -1
|
|---|
| 138 | #MSimReflector.DetectorMargin: 0
|
|---|
| 139 |
|
|---|
| 140 | # -------------------------------------------------------------------------
|
|---|
| 141 | # Setup the dark counts (FrequencyFixed) and the NSB noise per cm^2
|
|---|
| 142 | # -------------------------------------------------------------------------
|
|---|
| 143 | # Dark Counts per APD: ~4MHz
|
|---|
| 144 | #MSimRandomPhotons.FrequencyFixed: 0.004
|
|---|
| 145 | # NSB photon rate per cm^2 ~40MHz (folded with the cones' angular
|
|---|
| 146 | # acceptance and the wavelength acceptance of the camera (window, apd, etc)
|
|---|
| 147 | #MSimRandomPhotons.FrequencyNSB: 8.26e-4
|
|---|
| 148 |
|
|---|
| 149 | # FIXME: With a class describing the cones we could give NSB as
|
|---|
| 150 | # per sr and cm^2
|
|---|
| 151 |
|
|---|
| 152 | # -------------------------------------------------------------------------
|
|---|
| 153 | # Setup the trigger
|
|---|
| 154 | # -------------------------------------------------------------------------
|
|---|
| 155 | # This line could be omitted but then the discriminator would be
|
|---|
| 156 | # evaluated for all pixels not just for the pixels which are
|
|---|
| 157 | # later "connected" in the trigger (used in the coincidence map)
|
|---|
| 158 | # MSimTrigger.FileNameRouteAC: resmc/magic-trigger-1NN.txt
|
|---|
| 159 |
|
|---|
| 160 | # DiscriminatorThreshold is about 21 phe (->6888 DAC)
|
|---|
| 161 | MSimTrigger.FileNameRouteAC: /fhgfs/groups/app/mc_workshop/datacenter/setup/ceres/030/fact-trigger-sum.txt
|
|---|
| 162 | MSimTrigger.DiscriminatorThreshold: 6888
|
|---|
| 163 |
|
|---|
| 164 | # HEGRA like single pixel above threshold trigger
|
|---|
| 165 | #MSimTrigger.DiscriminatorThreshold: 10.00
|
|---|
| 166 | # Is this an option for a single pixel trigger?
|
|---|
| 167 | #MSimTrigger.DigitalSignalLength: -1
|
|---|
| 168 | #MSimTrigger.CoincidenceTime: 3
|
|---|
| 169 |
|
|---|
| 170 | # Magic like 2NN trigger (inner 396 pixel)
|
|---|
| 171 | #MSimTrigger.FileNameCoincidenceMap: resmc/magic-trigger-2NN.txt
|
|---|
| 172 | #MSimTrigger.DiscriminatorThreshold: 7.00
|
|---|
| 173 | #MSimTrigger.DigitalSignalLength: 8
|
|---|
| 174 | #MSimTrigger.CoincidenceTime: 0.5
|
|---|
| 175 |
|
|---|
| 176 | # Magic like 3NN trigger (inner 396 pixel)
|
|---|
| 177 | #MSimTrigger.FileNameCoincidenceMap: resmc/magic-trigger-3NN.txt
|
|---|
| 178 | #MSimTrigger.DiscriminatorThreshold: 6.0
|
|---|
| 179 | #MSimTrigger.DigitalSignalLength: 8
|
|---|
| 180 | #MSimTrigger.CoincidenceTime: 0.5
|
|---|
| 181 |
|
|---|
| 182 | # Magic like 4NN trigger (inner 396 pixel)
|
|---|
| 183 | #MSimTrigger.FileNameCoincidenceMap: resmc/magic-trigger-4NN.txt
|
|---|
| 184 | #MSimTrigger.DiscriminatorThreshold: 5.0
|
|---|
| 185 | #MSimTrigger.DigitalSignalLength: 8
|
|---|
| 186 | #MSimTrigger.CoincidenceTime: 0.5
|
|---|
| 187 |
|
|---|
| 188 | #MSimTrigger.FileNameRouteAC: resmc/fact-7map.txt
|
|---|
| 189 | #MSimTrigger.DiscriminatorThreshold: 14040
|
|---|
| 190 | #MSimTrigger.DigitalSignalLength: 8
|
|---|
| 191 | #MSimTrigger.CoincidenceTime: 0.5
|
|---|
| 192 |
|
|---|
| 193 | # Every Pixel(!) should see the same signal independant of its size
|
|---|
| 194 | MSimCalibrationSignal.NumPhotons: 24
|
|---|
| 195 | MSimCalibrationSignal.NumEvents: 1000
|
|---|
| 196 |
|
|---|
| 197 | IntendedPulsePos.Val: 30
|
|---|
| 198 |
|
|---|
| 199 | #PulseShape.Function.Name: exp(-(x/2)^2/2)
|
|---|
| 200 | #PulseShape.FileName: /fhgfs/groups/app/mc_workshop/datacenter/setup/ceres/030/template_rising_edge_maxprob_in_ns.txt
|
|---|
| 201 |
|
|---|
| 202 | # Pulse Shape:
|
|---|
| 203 | #
|
|---|
| 204 | # The number of sampling points is almost irrelevant because they
|
|---|
| 205 | # are equidistant, i.e. calculated and no search is necessary.
|
|---|
| 206 | # Nevertheless, you must make sure that there are enough points
|
|---|
| 207 | # to sample the function accuratly enough.
|
|---|
| 208 | #
|
|---|
| 209 | # The fitted function follows this expression:
|
|---|
| 210 | # A*(1-exp(-x/tau1))exp(-x/tau2)
|
|---|
| 211 | # A: Asymptotic maximum of the rising exp. function
|
|---|
| 212 | # tau1: time constant of the rising exp.function
|
|---|
| 213 | # tau2: time constant of the falling exp. function
|
|---|
| 214 | PulseShape.Function.Name: (1.626*(1-exp(-0.3803*x))*exp(-0.0649*x))
|
|---|
| 215 | PulseShape.Function.Npx: 250
|
|---|
| 216 | PulseShape.Function.Xmin: 0
|
|---|
| 217 | PulseShape.Function.Xmax: 100
|
|---|
| 218 |
|
|---|
| 219 | # -------------------------------------------------------------------------
|
|---|
| 220 | # Description how a pedestal/calibration signal is produced
|
|---|
| 221 | # -------------------------------------------------------------------------
|
|---|
| 222 |
|
|---|
| 223 | #MSimCalibrationSignal.NumEvents: 1000
|
|---|
| 224 | #MSimCalibrationSignal.NumPhotons; 5
|
|---|
| 225 | #MSimCalibrationSignal.TimeJitter; 1.0
|
|---|
| 226 |
|
|---|
| 227 |
|
|---|
| 228 | # -------------------------------------------------------------------------
|
|---|
| 229 | # Setup the FADC
|
|---|
| 230 | # -------------------------------------------------------------------------
|
|---|
| 231 |
|
|---|
| 232 | # Currently only 1GHz is supported
|
|---|
| 233 | MRawRunHeader.SamplingFrequency: 1000
|
|---|
| 234 | MRawRunHeader.NumSamples: 150
|
|---|
| 235 | MRawRunHeader.NumBytesPerSample: 2
|
|---|
| 236 | MRawRunHeader.FadcResolution: 12
|
|---|
| 237 |
|
|---|
| 238 | # Old values
|
|---|
| 239 | # Baseline at 15 like in MAGIC
|
|---|
| 240 | # 1.5 bit noise for a gain of 64
|
|---|
| 241 | # This allows a maximum of 200phe/pix
|
|---|
| 242 | #MSimCamera.DefaultOffset: 960
|
|---|
| 243 | #MSimCamera.DefaultNoise: 96
|
|---|
| 244 | #MSimCamera.DefaultGain: 256
|
|---|
| 245 |
|
|---|
| 246 | # New values
|
|---|
| 247 | #MSimCamera.DefaultOffset: 1000
|
|---|
| 248 | #MSimCamera.DefaultNoise: 30
|
|---|
| 249 | #MSimCamera.DefaultGain: 190
|
|---|
| 250 |
|
|---|
| 251 | # New values from 7.12.2012
|
|---|
| 252 | # gain of one phe is ~10mV, corresponds to 328 DAC (16 bit with an dynamic range of 2 V)
|
|---|
| 253 | # offset is 10 phe (->3280 DAC)
|
|---|
| 254 | # electronic noise is about 1.5 phe (->50 DAC)
|
|---|
| 255 | MSimCamera.DefaultOffset: 3280
|
|---|
| 256 | MSimCamera.DefaultNoise: 50
|
|---|
| 257 | MSimCamera.DefaultGain: 328
|
|---|
| 258 |
|
|---|
| 259 |
|
|---|
| 260 |
|
|---|
| 261 | #MSimReadout.fConversionFactor: 1
|
|---|
| 262 |
|
|---|
| 263 | # The number of sampling points is almost irrelevant because they
|
|---|
| 264 | # are equidistant, i.e. calculated and no search is necessary.
|
|---|
| 265 | # Nevertheless, you must make sure that there are enough points
|
|---|
| 266 | # to sample the function accuratly enough.
|
|---|
| 267 | #MPulseShape.Function.Name: 3.51308*exp(-1.4105*(exp(-0.5*x)+0.310942*x)) + 4.33521e-3*exp(-1.4105*(exp(-0.5*x)+6.09265e-3*x))
|
|---|
| 268 | #MPulseShape.Function.Npx: 250
|
|---|
| 269 | #MPulseShape.Function.Xmin: -5
|
|---|
| 270 | #MPulseShape.Function.Xmax: 45
|
|---|
| 271 |
|
|---|
| 272 | # -------------------------------------------------------------------------
|
|---|
| 273 | # Setup an image cleaning on the pure gamma signal (without noise)
|
|---|
| 274 | # -------------------------------------------------------------------------
|
|---|
| 275 |
|
|---|
| 276 | #MImgCleanStd.CleanLevel1: 6.0
|
|---|
| 277 | #MImgCleanStd.CleanLevel2: 6.0
|
|---|
| 278 |
|
|---|
| 279 | # -------------------------------------------------------------------------
|
|---|
| 280 | # This is a cut executed after the calculation of the image parameters
|
|---|
| 281 | # -------------------------------------------------------------------------
|
|---|
| 282 | Cut.Inverted: Yes
|
|---|
| 283 |
|
|---|
| 284 | #Cut.Condition: MNewImagePar.fNumUsedPixels>3 && MHillas.fSize>6
|
|---|
| 285 | Cut.Condition: MHillas.fSize>10
|
|---|
| 286 | #&& MHillas.fSize>30
|
|---|
| 287 | #MSimTrigger.SimulateElectronics: Off
|
|---|
| 288 |
|
|---|
| 289 | # Does not trigger anyway
|
|---|
| 290 | ContEmpty3.Condition: MPhotonEvent.GetNumPhotons<10
|
|---|
