Changeset 17998 for branches


Ignore:
Timestamp:
10/28/14 14:34:42 (10 years ago)
Author:
toscano
Message:
Ceres configuration file with the new pulse shape from the calibration paper
File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/MarsISDCBranchBasedOn17887/ceres.rc

    r17845 r17998  
    1 # ==========================================================================
    2 #############################################################################
    3 # ==========================================================================
    4 #                              General
    5 # ==========================================================================
    6 #############################################################################
    7 # ==========================================================================
     1MLog.VerbosityLevel: 4
    82
    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
     3BinningImpact.Raw: 40 0 100
    164
    17 # ==========================================================================
    18 #############################################################################
    19 # ==========================================================================
    20 #                                   Ceres
    21 # ==========================================================================
    22 #############################################################################
    23 # ==========================================================================
     5BinningTrigPos.Raw: 300 -25 275
    246
    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: .
     7RandomNumberGenerator: TRandom3
    318
     9Reflector.Constructor: MReflector
     10Reflector.FileName: resmc/fact/fact-reflector-real_ceres.txt
     11Reflector.SetSigmaPSF: 3.000000
    3212
    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/fact/atmopshere-aerosols.txt
    80 #MSimAtmosphere.FileOzone:    resmc/fact/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: resmc/fact/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);
    12113MGeomCam.Constructor:   MGeomCamFACT();
    12214
    123 #MSimBundlePhotons.FileName: resmc/fact/dwarf-fact.txt
     15MSimRandomPhotons.FrequencyFixed: 0.004000
     16MSimAPD.Type:                   0
     17MSimAPD.NumCells:               60
     18MSimAPD.DeadTime:               3.000000
     19MSimAPD.RecoveryTime:   8.750000
     20MSimAPD.CrosstalkCoefficient:   0.130000
     21MSimAPD.AfterpulseProb1:                0.140000
     22MSimAPD.AfterpulseProb2:                0.110000
     23MSimExcessNoise.ExcessNoise:    0.120000
    12424
    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 # -------------------------------------------------------------------------
    13225MirrorReflectivity.FileName:           resmc/fact/MirrorReflectivity_Lustermann_FACT_bearbeitet.txt
    13326PhotonDetectionEfficiency.FileName:    resmc/fact/fact-pde-1.0V.txt
     
    13528ConesTransmission.FileName:            resmc/fact/Transmittance_1439Cones_FACT_bearbeitet.txt
    13629
    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.FrequencyNSB:   0.025
    152 
    153 # FIXME: With a class describing the cones we could give NSB as
    154 # per sr and cm^2
    155 
    156 # -------------------------------------------------------------------------
    157 # Setup the trigger
    158 # -------------------------------------------------------------------------
    159 # This line could be omitted but then the discriminator would be
    160 # evaluated for all pixels not just for the pixels which are
    161 # later "connected" in the trigger (used in the coincidence map)
    162 # MSimTrigger.FileNameRouteAC:         resmc/fact/magic-trigger-1NN.txt
     30MSimRandomPhotons.FileNameNSB: resmc/fact/night-sky-la-palma.txt
     31MSimRandomPhotons.FrequencyNSB: 0.000000
    16332
    16433MSimTrigger.FileNameRouteAC: resmc/fact/fact-trigger-sum.txt
     34MSimTrigger.DiscriminatorThreshold: -16180.125000
    16535
    166 # DiscriminatorThreshold is 9*offset(-1850)+21*gain(22.375)=16180.125
    167 MSimTrigger.DiscriminatorThreshold: -16180.125
    168 
    169 # Every Pixel(!) should see the same signal independant of its size
    17036MSimCalibrationSignal.NumPhotons: 24
    17137MSimCalibrationSignal.NumEvents:  1000
     
    17339IntendedPulsePos.Val: 25
    17440
    175 #PulseShape.Function.Name: exp(-(x/2)^2/2)
    176 #PulseShape.FileName: resmc/fact/dwarf-pulse.txt
     41MRawRunHeader.SamplingFrequency: 2000
     42MRawRunHeader.NumSamples:        300
     43MRawRunHeader.NumBytesPerSample: 2
     44MRawRunHeader.FadcResolution:   
    17745
    178 # -------------------------------------------------------------------------
    179 # Description how a pedestal/calibration signal is produced
    180 # -------------------------------------------------------------------------
     46MSimCamera.DefaultOffset: -1850.000000
     47MSimCamera.DefaultNoise:    2.812500
     48MSimCamera.DefaultGain:    22.375000
    18149
    182 #MSimCalibrationSignal.NumEvents:  1000
    183 #MSimCalibrationSignal.NumPhotons;    5
    184 #MSimCalibrationSignal.TimeJitter;  1.0
     50MSimCamera.fACFudgeFactor: 0.336700
     51MSimCamera.fACTimeConstant: 20.000000
    18552
    186 
    187 # -------------------------------------------------------------------------
    188 # Setup the FADC
    189 # -------------------------------------------------------------------------
    190 
    191 MRawRunHeader.SamplingFrequency: 2000
    192 MRawRunHeader.NumSamples:         300
    193 MRawRunHeader.NumBytesPerSample:    2
    194 MRawRunHeader.FadcResolution:      12
    195 
    196 MSimCamera.DefaultOffset: -1850.0
    197 MSimCamera.DefaultNoise:    2.8125
    198 MSimCamera.DefaultGain:    22.375
    199 
    200 # Value for the fudgefactor in the calculation of the accoupling:
    201 MSimCamera.ACFudgeFactor: 0.3367
    202 MSimCamera.ACTimeConstant: 20
    203 
    204 #MSimReadout.fConversionFactor: 1
    205 
    206 # The number of sampling points is almost irrelevant because they
    207 # are equidistant, i.e. calculated and no search is necessary.
    208 # Nevertheless, you must make sure that there are enough points
    209 # to sample the function accuratly enough.
    210 # Attention: x in the function is given in slices, so if you change the sampling
    211 # frequency you have to change also this function
    212 PulseShape.Function.Name:  (1.626*(1-exp(-0.3803*0.5*x))*exp(-0.0649*0.5*x))
     53PulseShape.Function.Name:  1.2384*(1-1/(1+exp(-(5.8-x)/2./1.07)))*exp((5.8-x)/2./19.4)
    21354PulseShape.Function.Npx:   300
    214 PulseShape.Function.Xmin:   0
     55PulseShape.Function.Xmin:   -10
    21556PulseShape.Function.Xmax:   300
    21657
    217 # -------------------------------------------------------------------------
    218 # Setup an image cleaning on the pure gamma signal (without noise)
    219 # -------------------------------------------------------------------------
     58Cut.Inverted: yes
     59Cut.Condition: MHillas.fSize>10.000000
    22060
    221 #MImgCleanStd.CleanLevel1:   6.0
    222 #MImgCleanStd.CleanLevel2:   6.0
    223 
    224 # -------------------------------------------------------------------------
    225 # This is a cut executed after the calculation of the image parameters
    226 # -------------------------------------------------------------------------
    227 Cut.Inverted: Yes
    228 
    229 #Cut.Condition: MNewImagePar.fNumUsedPixels>3 && MHillas.fSize>6
    230 Cut.Condition: MHillas.fSize>10
    231 #&& MHillas.fSize>30
    232 #MSimTrigger.SimulateElectronics: Off
    233 
    234 # Does not trigger anyway
    23561ContEmpty3.Condition: MPhotonEvent.GetNumPhotons<10
Note: See TracChangeset for help on using the changeset viewer.