Changeset 17998 for branches/MarsISDCBranchBasedOn17887

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:

• branches/MarsISDCBranchBasedOn17887/ceres.rc (modified) (3 diffs)

branches/MarsISDCBranchBasedOn17887/ceres.rc

@@ -1 @@
-# ==========================================================================
-#############################################################################
-# ==========================================================================
-#                              General
-# ==========================================================================
-#############################################################################
-# ==========================================================================
+MLog.VerbosityLevel: 4
 
-# -------------------------------------------------------------------------
-# Use this if you want to setup the logging stream for the jobs
-# (overwrites command line options)
-# -------------------------------------------------------------------------
-MLog.VerbosityLevel: 4
-#MLog.DebugLevel:     1
-#MLog.NoColors:       yes
+BinningImpact.Raw: 40 0 100
 
-# ==========================================================================
-#############################################################################
-# ==========================================================================
-#                                   Ceres
-# ==========================================================================
-#############################################################################
-# ==========================================================================
+BinningTrigPos.Raw: 300 -25 275
 
-# -------------------------------------------------------------------------
-# Use this if you want to write the MJStar output somewhere
-# If you don't want it, it is written to the calibration output anyhow.
-# -------------------------------------------------------------------------
-#PathOut: .
-#PathIn: .
+RandomNumberGenerator: TRandom3
 
+Reflector.Constructor: MReflector
+Reflector.FileName: resmc/fact/fact-reflector-real_ceres.txt
+Reflector.SetSigmaPSF: 3.000000
 
-# -------------------------------------------------------------------------
-# Configure Eventloop
-# -------------------------------------------------------------------------
-#MaxEvents: 10000
-#Overwrite: yes,no
-
-
-# -------------------------------------------------------------------------
-# Use this to setup binnings. For more details see: MBinning::ReadEnv
-# -------------------------------------------------------------------------
-# BinningEnergy.Raw:      100      1    100000  log
-# BinningSize.Raw:        100      1  10000000  log
-BinningImpact.Raw:       40   0       1000
-# BinningHeight.Raw:       50      0        12
-# BinningAz.Raw:          360   -360       360
-# BinningZd.Raw:           70      0        70
-# BinningViewCone.Raw:    155      0        31
-BinningTrigPos.Raw:     300    -25       275
-# BinningEvtWidth.Raw:    150      0        25
-# BinningDist.Raw:        100      0       2.5
-# BinningDistC.Raw:       100      0       2.5
-# BinningThreshold.Raw:    35    0.9     90000  log
-# BinningEnergyEst.Raw:   100    0.9     90000  log
-
-
-# -------------------------------------------------------------------------
-# Initialize random number generator (see MJob::InitRandomNumberGenerator)
-# -------------------------------------------------------------------------
-RandomNumberGenerator: TRandom3
-#RandomNumberSeedValue: 0
-
-
-# -------------------------------------------------------------------------
-# Ceres general setup
-# -------------------------------------------------------------------------
-
-# Switch off the camera "electronics"
-#Camera: Off
-
-# Force the use of the "hardware" trigger for calibration data
-#ForceTrigger: Yes
-
-
-# -------------------------------------------------------------------------
-# Some setup for the atmosphere. The default should be well suited.
-# -------------------------------------------------------------------------
-#MSimAtmosphere.FileAerosols: resmc/fact/atmopshere-aerosols.txt
-#MSimAtmosphere.FileOzone:    resmc/fact/atmopshere-ozone.txt
-
-
-# -------------------------------------------------------------------------
-# Here you can control the poiting of the telescope. To switch on
-# off-target observations set a value for the distance !=0 [deg].
-# For details see MSimPointingPos
-# -------------------------------------------------------------------------
-#MSimPointingPos.OffTargetDistance: 0.4
-#MSimPointingPos.OffTargetPhi: -1
-
-
-# -------------------------------------------------------------------------
-# Setup the reflector and camera geometry
-# -------------------------------------------------------------------------
-Reflector.Constructor: MReflector
-
-# For the file definition see MReflector::ReadFile
-
-Reflector.FileName: resmc/fact/fact-reflector.txt
-# Units mm: ~0.02° 
-#Reflector.SetSigmaPSF: 1.8
-# Units mm: ~0.03°
-#Reflector.SetSigmaPSF: 2.6
-#Reflector.SetSigmaPSF: 1.66
-
-# new Value from 7.12.2012:
-# 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
-#Reflector.SetSigmaPSF: 1.3
-
-# new Value from 29.05.2013:
-# Set SigmaPSF to pixelsize
-Reflector.SetSigmaPSF: 10.0
-
-# --- DWARF ---
-# distance from mirror to camera [cm]
-#focal_distance     489.73
-# TRUE focal of the paraboloid (must correspond to that of the mirror positions defined below!)
-#paraboloid_focal   488.87
-
-#MGeomCam.Constructor:   MGeomCamDwarf(189.8,  9.5,  4.8887);
 MGeomCam.Constructor:   MGeomCamFACT();
 
-#MSimBundlePhotons.FileName: resmc/fact/dwarf-fact.txt
+MSimRandomPhotons.FrequencyFixed: 0.004000
+MSimAPD.Type:                   0
+MSimAPD.NumCells:               60
+MSimAPD.DeadTime:               3.000000
+MSimAPD.RecoveryTime:   8.750000
+MSimAPD.CrosstalkCoefficient:   0.130000
+MSimAPD.AfterpulseProb1:                0.140000
+MSimAPD.AfterpulseProb2:                0.110000
+MSimExcessNoise.ExcessNoise:    0.120000
 
-# Set the APD type (1: 30x30 <default>, 2: 60x60, 3:60x60(ct=15%))
-MSimAPD.Type: 4
-MSimExcessNoise.ExcessNoise: 0.12
-
-# -------------------------------------------------------------------------
-# Setup the absorption, conversion efficiency and angular acceptance
-# -------------------------------------------------------------------------
 MirrorReflectivity.FileName:           resmc/fact/MirrorReflectivity_Lustermann_FACT_bearbeitet.txt
 PhotonDetectionEfficiency.FileName:    resmc/fact/fact-pde-1.0V.txt
@@ -135 +28 @@
 ConesTransmission.FileName:            resmc/fact/Transmittance_1439Cones_FACT_bearbeitet.txt
 
-# -------------------------------------------------------------------------
-# Setup what in MMCS would be called "additional spot size"
-# -------------------------------------------------------------------------
-#MSimPSF.Sigma: -1
-#MSimReflector.DetectorMargin: 0
-
-# -------------------------------------------------------------------------
-# Setup the dark counts (FrequencyFixed) and the NSB noise per cm^2
-# -------------------------------------------------------------------------
-# Dark Counts per APD: ~4MHz
-#MSimRandomPhotons.FrequencyFixed: 0.004
-# NSB photon rate per cm^2 ~40MHz (folded with the cones' angular
-# acceptance and the wavelength acceptance of the camera (window, apd, etc)
-# 0.040 1/ns/cm^2 NSB-rate:
-#MSimRandomPhotons.FrequencyNSB:   0.025
-
-# FIXME: With a class describing the cones we could give NSB as
-# per sr and cm^2
-
-# -------------------------------------------------------------------------
-# Setup the trigger
-# -------------------------------------------------------------------------
-# This line could be omitted but then the discriminator would be
-# evaluated for all pixels not just for the pixels which are
-# later "connected" in the trigger (used in the coincidence map)
-# MSimTrigger.FileNameRouteAC:         resmc/fact/magic-trigger-1NN.txt
+MSimRandomPhotons.FileNameNSB: resmc/fact/night-sky-la-palma.txt
+MSimRandomPhotons.FrequencyNSB: 0.000000
 
 MSimTrigger.FileNameRouteAC: resmc/fact/fact-trigger-sum.txt
+MSimTrigger.DiscriminatorThreshold: -16180.125000
 
-# DiscriminatorThreshold is 9*offset(-1850)+21*gain(22.375)=16180.125
-MSimTrigger.DiscriminatorThreshold: -16180.125
-
-# Every Pixel(!) should see the same signal independant of its size
 MSimCalibrationSignal.NumPhotons: 24
 MSimCalibrationSignal.NumEvents:  1000
@@ -173 +39 @@
 IntendedPulsePos.Val: 25
 
-#PulseShape.Function.Name: exp(-(x/2)^2/2)
-#PulseShape.FileName: resmc/fact/dwarf-pulse.txt
+MRawRunHeader.SamplingFrequency: 2000
+MRawRunHeader.NumSamples:        300
+MRawRunHeader.NumBytesPerSample: 2
+MRawRunHeader.FadcResolution:    
 
-# -------------------------------------------------------------------------
-# Description how a pedestal/calibration signal is produced
-# -------------------------------------------------------------------------
+MSimCamera.DefaultOffset: -1850.000000
+MSimCamera.DefaultNoise:    2.812500
+MSimCamera.DefaultGain:    22.375000
 
-#MSimCalibrationSignal.NumEvents:  1000
-#MSimCalibrationSignal.NumPhotons;    5
-#MSimCalibrationSignal.TimeJitter;  1.0
+MSimCamera.fACFudgeFactor: 0.336700
+MSimCamera.fACTimeConstant: 20.000000
 
-
-# -------------------------------------------------------------------------
-# Setup the FADC
-# -------------------------------------------------------------------------
-
-MRawRunHeader.SamplingFrequency: 2000
-MRawRunHeader.NumSamples:         300
-MRawRunHeader.NumBytesPerSample:    2
-MRawRunHeader.FadcResolution:      12
-
-MSimCamera.DefaultOffset: -1850.0
-MSimCamera.DefaultNoise:    2.8125
-MSimCamera.DefaultGain:    22.375
-
-# Value for the fudgefactor in the calculation of the accoupling:
-MSimCamera.ACFudgeFactor: 0.3367
-MSimCamera.ACTimeConstant: 20
-
-#MSimReadout.fConversionFactor: 1
-
-# The number of sampling points is almost irrelevant because they
-# are equidistant, i.e. calculated and no search is necessary.
-# Nevertheless, you must make sure that there are enough points
-# to sample the function accuratly enough.
-# Attention: x in the function is given in slices, so if you change the sampling
-# frequency you have to change also this function
-PulseShape.Function.Name:  (1.626*(1-exp(-0.3803*0.5*x))*exp(-0.0649*0.5*x))
+PulseShape.Function.Name:  1.2384*(1-1/(1+exp(-(5.8-x)/2./1.07)))*exp((5.8-x)/2./19.4)
 PulseShape.Function.Npx:   300
-PulseShape.Function.Xmin:   0
+PulseShape.Function.Xmin:   -10
 PulseShape.Function.Xmax:   300
 
-# -------------------------------------------------------------------------
-# Setup an image cleaning on the pure gamma signal (without noise)
-# -------------------------------------------------------------------------
+Cut.Inverted: yes
+Cut.Condition: MHillas.fSize>10.000000
 
-#MImgCleanStd.CleanLevel1:   6.0
-#MImgCleanStd.CleanLevel2:   6.0
-
-# -------------------------------------------------------------------------
-# This is a cut executed after the calculation of the image parameters
-# -------------------------------------------------------------------------
-Cut.Inverted: Yes
-
-#Cut.Condition: MNewImagePar.fNumUsedPixels>3 && MHillas.fSize>6
-Cut.Condition: MHillas.fSize>10
-#&& MHillas.fSize>30
-#MSimTrigger.SimulateElectronics: Off
-
-# Does not trigger anyway
 ContEmpty3.Condition: MPhotonEvent.GetNumPhotons<10