Index: trunk/Mars/callisto.rc =================================================================== --- trunk/Mars/callisto.rc (revision 19871) +++ (revision ) @@ -1,458 +1,0 @@ -# ========================================================================== -############################################################################# -# ========================================================================== -# General -# ========================================================================== -############################################################################# -# ========================================================================== - -# ------------------------------------------------------------------------- -# Use this if you want to setup the logging stream for the jobs -# (overwrites command line options) -# ------------------------------------------------------------------------- -#MLog.VerbosityLevel: 2 -#MLog.DebugLevel: 1 -#MLog.NoColors: yes - -# ========================================================================== -############################################################################# -# ========================================================================== -# Calibration -# ========================================================================== -############################################################################# -# ========================================================================== - -# ------------------------------------------------------------------------- -# Use this if you want to write the MJPedestal output somewhere -# If you don't want it, it is written to the calibration output anyhow. -# ------------------------------------------------------------------------- -#MJPedestalC1.PathOut: . -#MJPedestalC2.PathOut: . -#MJPedestalC1.BadPixelsFile: mjobs/badpixels_0_559.rc -#MJPedestalC1.ReferenceFile: mjobs/pedestalref.rc -#MJPedestalC2.ReferenceFile: mjobs/pedestalref.rc -MJCalibration.ReferenceFile: mjobs/calibrationref_Nov04.rc -MJCalibration.MHCalibrationRelTimeCam.ReferenceFile: mjobs/calibrationref_Nov04.rc -MJCalibration.MHCalibrationChargeCam.ReferenceFile: mjobs/calibrationref_Nov04.rc - -# ------------------------------------------------------------------------- -# Use this to define where the calibration output is stored. The filename -# is created from the sequence number. If nothing is specified '.' is -# assumed. (overwrites command line options) -# ------------------------------------------------------------------------- -#MJCalibration.PathOut: calped - -# ------------------------------------------------------------------------- -# Define here which signal extractor you are going to use for the -# whole calibration chain -# ------------------------------------------------------------------------- -MJPedestalC1.ExtractSignal: MExtractTimeAndChargeDigitalFilter -# ------------------------------------------------------------------------- -# Define here an extractor which can be used for the December 04 data. -# ------------------------------------------------------------------------- -#MJPedestalC1.ExtractSignal: MExtractTimeAndChargeSpline, MExtractTimeAndChargeDigitalFilterPeakSearch, MExtractTimeAndChargeSlidingWindow -# ------------------------------------------------------------------------- -# Define here parameters valid for all extractors above -# The actual defaults used depend on the extractor in use. -# For details see the corresponding source code -# ------------------------------------------------------------------------- -#MJPedestalC1.ExtractSignal.HiGainFirst: 0 -#MJPedestalC1.ExtractSignal.HiGainLast: 15 -#MJPedestalC1.ExtractSignal.LoGainFirst: 3 -#MJPedestalC1.ExtractSignal.LoGainLast: 14 -#MJPedestalC1.ExtractSignal.SaturationLimit: 250 -#MJPedestalC1.ExtractSignal.OffsetLoGain: 0.95 -#MJPedestalC1.ExtractSignal.LoGainSwitch: 120 -#MJPedestalC1.ExtractSignal.LoGainStartShift: -2.5 -# ------------------------------------------------------------------------- -# Define here parameters valid for the spline: -# ------------------------------------------------------------------------- -#MJPedestalC1.ExtractSignal.Resolution: 0.05 -#MJPedestalC1.ExtractSignal.RiseTimeHiGain: 0.5 -#MJPedestalC1.ExtractSignal.FallTimeHiGain: 1.5 -#MJPedestalC1.ExtractSignal.LoGainStretch: 1.5 -#MJPedestalC1.ExtractSignal.ExtractionType: Integral -#MJPedestalC1.ExtractSignal.ExtractionType: Amplitude - -# ------------------------------------------------------------------------- -# In the case of the Digital filter you have to overwrite the default -# weights file depending on what you are extracting -# ------------------------------------------------------------------------- -MJPedestalC2.ExtractSignal.WeightsFile: msignal/calibration_weights_UV46.dat -MJPedestalY2.ExtractSignal.WeightsFile: msignal/cosmics_weights46.dat -MJPedestalY3.ExtractSignal.WeightsFile: msignal/cosmics_weights46.dat -MJCalibration.ExtractSignal.WeightsFile: msignal/calibration_weights_UV46.dat -MJCalibrateSignal.ExtractSignal.WeightsFile: msignal/cosmics_weights46.dat -MJCalibrateSignal.ExtractInterlaced.WeightsFile: msignal/calibration_weights_UV46.dat - -#MJPedestalC2.ExtractSignal.AutomaticWeightsFile: On -#MJPedestalY2.ExtractSignal.AutomaticWeightsFile: On -#MJPedestalY3.ExtractSignal.AutomaticWeightsFile: On -#MJCalibration.ExtractSignal.AutomaticWeightsFile: On -#MJCalibrateSignal.ExtractSignal.AutomaticWeightsFile: On -#MJCalibrateSignal.ExtractInterlaced.AutomaticWeightsFile: On - -# ------------------------------------------------------------------------- -# Configure MJCalibration -# ------------------------------------------------------------------------- - -# Switch on relative time calibration -#MJCalibration.RelTimeCalibration: Yes -# Type of used data format: raw, root, mc -#MJCalibration.DataType: Root -# Type of displayed plots: Full,DataCheck,Normal -#MJCalibration.Display: datacheck -# Write additinal debug output -#MJCalibration.Debug: No -# Use blind pixel -#MJCalibration.UseBlindPixel: No -# Use pin diode -#MJCalibration.UsePINDiode: No -# Use pixel check option -#MJCalibration.PixelCheck: Yes -#MJCalibTest.PixelCheck: Yes - -# ------------------------------------------------------------------------- -# Individual events check for calibration results -# ------------------------------------------------------------------------- - -#ContCosmics.MaxEmptyPixels: 0.2 -#ContCosmics.MaxExcludedFraction: 1 -#ContCosmics.MinAcceptedFraction: 0 - -# ------------------------------------------------------------------------- -# You can choose the pedestal extraction algorithm/task. To use -# MPedCalcPedRun use the lines below. Be carefull, a pedestal file -# could be (if no pedestal file available) a data-file, too. -# ------------------------------------------------------------------------- -#MJPedestalC1.ExtractPedestal: MPedCalcPedRun -#MJPedestalC1.ExtractPedestal.ExtractWinFirst: 0 -#MJPedestalC1.ExtractPedestal.ExtractWinSize: 6 -#MJPedestalC1.ExtractPedestal.PedestalUpdate: no -#MJPedestalC1.ExtractPedestal.CheckWinFirst: 0 -#MJPedestalC1.ExtractPedestal.CheckWinLast: 6 -#MJPedestalC1.ExtractPedestal.MaxSignalVar: 300 - -#MJPedestalC2.ExtractPedestal: MPedCalcPedRun -#MJPedestalC2.ExtractPedestal.ExtractWinFirst: 0 -#MJPedestalC2.ExtractPedestal.ExtractWinSize: 6 -#MJPedestalC2.ExtractPedestal.PedestalUpdate: no -#MJPedestalC2.ExtractPedestal.CheckWinFirst: 0 -#MJPedestalC2.ExtractPedestal.CheckWinSize: 29 - -# ------------------------------------------------------------------------- -# You can choose the pedestal extraction algorithm/task. To use -# MPedCalcFromLoGain use the lines below. Be carefull, a pedestal file -# could be (if no pedestal file available) a data-file, too. -# ------------------------------------------------------------------------- -#MJPedestalC1.ExtractPedestal: MPedCalcFromLoGain -#MJPedestalC1.ExtractPedestal.ExtractWinFirst: 0 -#MJPedestalC1.ExtractPedestal.ExtractWinSize: 6 -#MJPedestalC1.ExtractPedestal.PedestalUpdate: no -#MJPedestalC1.ExtractPedestal.CheckWinFirst: 0 -#MJPedestalC1.ExtractPedestal.CheckWinSize: 29 -#MJPedestalC1.ExtractPedestal.MaxSignalVar: 40 - -#MJPedestalC2.ExtractPedestal: MPedCalcFromLoGain -#MJPedestalC2.ExtractPedestal.ExtractWinFirst: 0 -#MJPedestalC2.ExtractPedestal.ExtractWinSize: 6 -#MJPedestalC2.ExtractPedestal.PedestalUpdate: no -#MJPedestalC2.ExtractPedestal.CheckWinFirst: 0 -#MJPedestalC2.ExtractPedestal.CheckWinSize: 29 -#MJPedestalC2.ExtractPedestal.MaxSignalVar: 40 - -# ------------------------------------------------------------------------- -# Configure MJPedestalC1 and MJPedestalC2 -# ------------------------------------------------------------------------- -# Possible data type: Raw,Root,MC -#MJPedestalC1.DataType: Root -#MJPedestalC2.DataType: Root -#MJPedestalC3.DataType: Root -# Type of displayed plots: Full,DataCheck,None -#MJPedestalC1.Display: datacheck -#MJPedestalC2.Display: datacheck -#MJPedestalC3.Display: datacheck -# Maximum number of event processed in the loop -#MJPedestalC1.MaxEvents: 1000 -#MJPedestalC2.MaxEvents: 1000 -# Allow to overwrite existing files with the output file -#MJPedestalC1.AllowOverwrite: No -#MJPedestalC2.AllowOverwrite: No -# Use data runs from the sequence instead of calibration runs -#MJPedestalC1.UseData: No -#MJPedestalC2.UseData: No -# Switch off DeadPixelCheck with "No" -#MJPedestalY2.DeadPixelCheck: yes - - -# ------------------------------------------------------------------------- -# Use this if you want to change the software low-gain switch in the -# calibration -# ------------------------------------------------------------------------- -MJCalibration.ExtractSignal.LoGainSwitch: 170 - -# ------------------------------------------------------------------------- -# Use this if you want to change the higain-vs. logain intercalibration -# ------------------------------------------------------------------------- -#MJCalibration.HiLoCalibFile: mjobs/hilocalib_df46.root - -# ------------------------------------------------------------------------- -# Use this to change the behaviour of the calibration -# ------------------------------------------------------------------------- -# Type if you set a colour explicitely from outside (TEST purposes!!) -#MJCalibration.MCalibColorSet.ExplicitColor: green,blue,uv,ct1 - -#MJCalibration.MCalibrationChargeCalc.ArrTimeRmsLimit: 5.0 -#MJCalibration.MCalibrationChargeCalc.ChargeLimit: 2.5 -#MJCalibration.MCalibrationChargeCalc.ChargeErrLimit: 0 -#MJCalibration.MCalibrationChargeCalc.ChargeRelErrLimit: 1 -#MJCalibration.MCalibrationChargeCalc.Debug: no - -#MJCalibration.MCalibrationChargeCalc.FFactorErrLimit: 4.5 -#MJCalibration.MCalibrationChargeCalc.LambdaErrLimit: 0.2 -#MJCalibration.MCalibrationChargeCalc.LambdaCheckLimit: 0.5 -#MJCalibration.MCalibrationChargeCalc.PheErrLowerLimit: 9.0 -#MJCalibration.MCalibrationChargeCalc.PheErrUpperLimit: 5.5 - -#MJCalibration.MHCalibrationChargeCam.Debug: no -#MJCalibration.MHCalibrationChargeCam.LoGain: yes -#MJCalibration.MHCalibrationChargeCam.Oscillations: yes -#MJCalibration.MHCalibrationChargeCam.SizeCheck: yes -#MJCalibration.MHCalibrationChargeCam.Averageing: yes -#MJCalibration.MHCalibrationChargeCam.HiGainNbins: 500 -#MJCalibration.MHCalibrationChargeCam.HiGainFirst: -98. -#MJCalibration.MHCalibrationChargeCam.HiGainLast: 1902. -#MJCalibration.MHCalibrationChargeCam.LoGainNbins: 500 -#MJCalibration.MHCalibrationChargeCam.LoGainFirst: -99. -#MJCalibration.MHCalibrationChargeCam.LoGainLast: 901. -#MJCalibration.MHCalibrationChargeCam.TimeLowerLimit: 1. -#MJCalibration.MHCalibrationChargeCam.TimeUpperLimit: 3. -#MJCalibration.MHCalibrationChargeCam.NumHiGainSaturationLimit: 0.02 -#MJCalibration.MHCalibrationChargeCam.NumLoGainSaturationLimit: 0.005 -#MJCalibration.MHCalibrationChargeCam.ProbLimit: 0.00000001 -#MJCalibration.MHCalibrationChargeCam.OverflowLimit: 0.005 -#MJCalibration.MHCalibrationChargeCam.PulserFrequency: 500 - -#MJCalibration.MHCalibrationRelTimeCam.Debug: no -#MJCalibration.MHCalibrationRelTimeCam.LoGain: no -#MJCalibration.MHCalibrationRelTimeCam.Oscillations: yes -#MJCalibration.MHCalibrationRelTimeCam.SizeCheck: yes -#MJCalibration.MHCalibrationRelTimeCam.Averageing: yes -#MJCalibration.MHCalibrationRelTimeCam.Nbins: 210 -#MJCalibration.MHCalibrationRelTimeCam.First: -5. -#MJCalibration.MHCalibrationRelTimeCam.Last: 10. -#MJCalibration.MHCalibrationRelTimeCam.NumHiGainSaturationLimit: 0.25 -#MJCalibration.MHCalibrationRelTimeCam.ProbLimit: 0.00001 -#MJCalibration.MHCalibrationRelTimeCam.OverflowLimit: 0.005 -#MJCalibration.MHCalibrationRelTimeCam.PulserFrequency: 500 - - -# ========================================================================== -############################################################################# -# ========================================================================== -# Signal extraction -# ========================================================================== -############################################################################# -# ========================================================================== - -# ------------------------------------------------------------------------- -# Use a OutputPath if you want to write the results to a file -# ------------------------------------------------------------------------- -#MJCalibrateSignal.PathOut: . -#MJPedestalY1.PathOut: . -#MJPedestalY2.PathOut: . -#MJPedestalY3.PathOut: . -#MJCalibrateSignal.PathIn: . - -# ------------------------------------------------------------------------- -# If you want to change the behaviour of the pedestal extraction... -# ------------------------------------------------------------------------- -# Possible data type: Raw,Root,MC -#MJPedestalY1.DataType: Root -#MJPedestalY2.DataType: Root -#MJPedestalY3.DataType: Root -#MJPedestalY1.UseData: Yes -#MJPedestalY2.UseData: Yes -#MJPedestalY3.UseData: Yes -#MJPedestalY1.MinEvents: 50 -#MJPedestalY2.MinEvents: 50 -#MJPedestalY3.MinEvents: 50 - -MJPedestalY1.MaxEvents: 1000 -MJPedestalY2.MaxEvents: 10000 -MJPedestalY3.MaxEvents: 1000 - -MJPedestalY1.MaxPedestals: 500 -MJPedestalY2.MaxPedestals: 500 -MJPedestalY3.MaxPedestals: 500 - -# ------------------------------------------------------------------------- -# When the starting pedestal is extracted from the first events of -# the data run also the average pulse position of these events -# is determined. If the pulse position check is switched on (default) -# the extraction window of the hi- and lo-gain extractor for the -# determined average pulse position is not allowed to exceed the -# extractors search region. In this scenario we assume that the -# start of the extraction window is at least similar to the -# rising edge of the pulse. -# To make sure that the search region and the extraction window -# does not touch increase the values. -# ------------------------------------------------------------------------- -#MJPedestalY2.PulsePosCheck: yes -#MJPedestalY2.ExtractWinLeft: 0.0 -#MJPedestalY2.ExtractWinRight: 0.0 - - -# ------------------------------------------------------------------------- -# Define the Pulse Position check parameters: -# ------------------------------------------------------------------------- -#MJPedestalY2.MHCalibrationPulseTimeCam.SaturationLimit: 245 -#MJPedestalY2.MHCalibrationPulseTimeCam.LowerSignalLimit: 85 -#MJPedestalY2.MHCalibrationPulseTimeCam.NumPixelsRequired: 3 -#MJPedestalY2.PixelCheck: no - -#MJPedestalY.ExtractPedestal: MPedCalcFromLoGain -#MJPedestalY.ExtractPedestal.PedestalUpdate: no - -# ------------------------------------------------------------------------- -# Setup pedestal extractor for running through data -# ------------------------------------------------------------------------- -#MJCalibrateSignal.MPedCalcFundamental.PedestalUpdate: yes -#MJCalibrateSignal.MPedCalcFundamental.NumDump: 500 -#MJCalibrateSignal.MPedCalcFundamental.NumEventsDump: 500 -#MJCalibrateSignal.MPedCalcFundamental.NumAreasDump: 500 -#MJCalibrateSignal.MPedCalcFundamental.NumSectorsDump: 500 -#MJCalibrateSignal.MPedCalcFundamental.ExtractWinFirst: 17 -#MJCalibrateSignal.MPedCalcFundamental.ExtractWinSize: 6 -#MJCalibrateSignal.MPedCalcFundamental.CheckWinFirst: 0 -#MJCalibrateSignal.MPedCalcFundamental.CheckWinSize: 29 -#MJCalibrateSignal.MPedCalcFundamental.MaxSignalVar: 40 - -#MJCalibrateSignal.MPedCalcFromExtractorRndm.PedestalUpdate: yes -#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumEventsDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumAreasDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumSectorsDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.ExtractWinFirst: 17 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.ExtractWinSize: 6 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.CheckWinFirst: 0 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.CheckWinSize: 29 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.MaxSignalVar: 40 - -#MJCalibrateSignal.MPedCalcFromExtractor.PedestalUpdate: yes -#MJCalibrateSignal.MPedCalcFromExtractor.NumDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractor.NumEventsDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractor.NumAreasDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractor.NumSectorsDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractor.ExtractWinFirst: 17 -#MJCalibrateSignal.MPedCalcFromExtractor.ExtractWinSize: 6 -#MJCalibrateSignal.MPedCalcFromExtractor.CheckWinFirst: 0 -#MJCalibrateSignal.MPedCalcFromExtractor.CheckWinSize: 29 -#MJCalibrateSignal.MPedCalcFromExtractor.MaxSignalVar: 40 - -# ------------------------------------------------------------------------- -# Configure the interlaced calibration updates -# ------------------------------------------------------------------------- -#MJCalibrateSignal.MCalibCalcFromPast.UpdateWithFFactorMethod: yes -#MJCalibrateSignal.MCalibCalcFromPast.NumEventsDump: 500 -#MJCalibrateSignal.MCalibCalcFromPast.UpdateNumPhes: yes -#MJCalibrateSignal.MCalibCalcFromPast.NumPhesDump: 10 -#MJCalibrateSignal.MCalibrationChargeCalc.Debug: no - -# ------------------------------------------------------------------------- -# Use this if you want to change the signal extractor -# ------------------------------------------------------------------------- -MJCalibrateSignal.ExtractSignal.LoGainSwitch: 150 -#MJCalibrateSignal.ExtractSignal.LoGainStretch: 1.5 -#MJCalibrateSignal.ExtractInterlaced.LoGainStretch: 1.5 - -# ------------------------------------------------------------------------- -# Setup calibration of data -# ------------------------------------------------------------------------- -# Type of used data format: raw,root,MC -#MJCalibrateSignal.DataType: Root -# Type if you set a colour explicitely from outside (only TESTS!!!) -#MJCalibrateSignal.MCalibColorSet.ExpicitColor: green,blue,uv,ct1 -#MJCalibrateSignal.MCalibrateData.PedestalFlag: Event -#MJCalibrateSignal.MCalibrateData.ScaleFactor: 1 -MJCalibrateSignal.MCalibrateData.FileNameScale: resources/calibration.rc -#MJCalibrateSignal.MCalibrateData.CalibrationMode: Default -#MJCalibrateSignal.MCalibrateData.CalibConvMinLimit: 0.01 -#MJCalibrateSignal.MCalibrateData.CalibConvMaxLimit: 5. -#MJCalibrateSignal.Interlaced: yes -#MJCalibrateSignal.RelTimesUpdate: yes - -# ------------------------------------------------------------------------- -# Setup level for determining of bad pixels -# ------------------------------------------------------------------------- -#MJCalibrateSignal.MBadPixelsCalc.PedestalLevel: 3.0 -#MJCalibrateSignal.MBadPixelsCalc.PedestalLevelVariance: 3.0 -#MJCalibrateSignal.MBadPixelsTreat.NumMinNeighbors: 3 -#MJCalibrateSignal.MBadPixelsTreat.MaxArrivalTimeDiff: 3.0 -#MJCalibrateSignal.MBadPixelsTreat.UseInterpolation: yes -#MJCalibrateSignal.MBadPixelsTreat.ProcessPedestalEvt: yes -#MJCalibrateSignal.MBadPixelsTreat.ProcessPedestalRun: no -#MJCalibrateSignal.MBadPixelsTreat.ProcessTimes: yes -#MJCalibrateSignal.MBadPixelsTreat.UseCentralPixel: no -#MJCalibrateSignal.MBadPixelsTreat.HardTreatment: no - -# ------------------------------------------------------------------------- -# Setup interlaced calibration events histogramming -# ------------------------------------------------------------------------- -#MJCalibrateSignal.MHCalibrationChargeCam.Debug: no -#MJCalibrateSignal.MHCalibrationChargeCam.LoGain: yes -MJCalibrateSignal.MHCalibrationChargeCam.Oscillations: no -#MJCalibrateSignal.MHCalibrationChargeCam.SizeCheck: yes -#MJCalibrateSignal.MHCalibrationChargeCam.Averageing: yes -#MJCalibrateSignal.MHCalibrationChargeCam.HiGainNbins: 500 -#MJCalibrateSignal.MHCalibrationChargeCam.HiGainFirst: -98. -#MJCalibrateSignal.MHCalibrationChargeCam.HiGainLast: 1902. -#MJCalibrateSignal.MHCalibrationChargeCam.HiGainLast: 1899.5 -MJCalibrateSignal.MHCalibrationChargeCam.LoGainNbins: 250 -MJCalibrateSignal.MHCalibrationChargeCam.LoGainFirst: -98 -MJCalibrateSignal.MHCalibrationChargeCam.LoGainLast: 902 -#MJCalibrateSignal.MHCalibrationChargeCam.TimeLowerLimit: 1. -#MJCalibrateSignal.MHCalibrationChargeCam.TimeUpperLimit: 3. -#MJCalibrateSignal.MHCalibrationChargeCam.NumHiGainSaturationLimit: 0.02 -#MJCalibrateSignal.MHCalibrationChargeCam.NumLoGainSaturationLimit: 0.005 -#MJCalibrateSignal.MHCalibrationChargeCam.ProbLimit: 0.0001 -#MJCalibrateSignal.MHCalibrationChargeCam.OverflowLimit: 0.005 -MJCalibrateSignal.MHCalibrationChargeCam.PulserFrequency: 50 - -# ------------------------------------------------------------------------- -# Setup relative time interlaced calibration -# ------------------------------------------------------------------------- -#MJCalibrateSignal.MHCalibrationRelTimeCam.Debug: no -#MJCalibrateSignal.MHCalibrationRelTimeCam.LoGain: no -MJCalibrateSignal.MHCalibrationRelTimeCam.Oscillations: no -#MJCalibrateSignal.MHCalibrationRelTimeCam.SizeCheck: yes -#MJCalibrateSignal.MHCalibrationRelTimeCam.Averageing: yes -#MJCalibrateSignal.MHCalibrationRelTimeCam.Nbins: 210 -#MJCalibrateSignal.MHCalibrationRelTimeCam.First: -5. -#MJCalibrateSignal.MHCalibrationRelTimeCam.Last: 10. -#MJCalibrateSignal.MHCalibrationRelTimeCam.NumHiGainSaturationLimit: 0.25 -#MJCalibrateSignal.MHCalibrationRelTimeCam.ProbLimit: 0.00001 -#MJCalibrateSignal.MHCalibrationRelTimeCam.OverflowLimit: 0.005 -MJCalibrateSignal.MHCalibrationRelTimeCam.PulserFrequency: 50 - -# ------------------------------------------------------------------------- -# Resource for movie encoding. For more details see IACT wiki -# ------------------------------------------------------------------------- -#MMoviePrepare.FirstSlice: 0 -#MMoviePrepare.LastSlice: 15 -#MovieFilter.Condition: MMovieData.fMax>5*MMovieData.fMedianPedestalRms -#MImgCleanStd.CleanLevel1: 8.0 -#MImgCleanStd.CleanLevel2: 4.5 -#MFEvtNumber.FileName: ganymed00001218.root -#MFEvtNumber.Selector: ThetaSquared.fVal<0.04 && DataType.fVal>0.5 -#MMovieWrite.TargetLength: 5 -#MMovieWrite.NumEvents: 10 -#MMovieWrite.Threshold: 3 -#MMovieWrite.Filename: movie.mpg - -# ------------------------------------------------------------------------- -# This is a special option for the datacheck. A task executed as last -# task in the event processing task list - in the standard software it -# is skipped -# ------------------------------------------------------------------------- -#MJCalibrateSignal.FinalFantasy: MHCamFlorian Index: trunk/Mars/callisto_Dec04Jan05.rc =================================================================== --- trunk/Mars/callisto_Dec04Jan05.rc (revision 19871) +++ (revision ) @@ -1,483 +1,0 @@ -# ========================================================================== -############################################################################# -# ========================================================================== -# General -# ========================================================================== -############################################################################# -# ========================================================================== - -# ------------------------------------------------------------------------- -# Use this if you want to setup the logging stream for the jobs -# (overwrites command line options) -# ------------------------------------------------------------------------- -#MLog.VerbosityLevel: 2 -#MLog.DebugLevel: 1 -#MLog.NoColors: yes - -# ========================================================================== -############################################################################# -# ========================================================================== -# Calibration -# ========================================================================== -############################################################################# -# ========================================================================== - -# ------------------------------------------------------------------------- -# Use this if you want to write the MJPedestal output somewhere -# If you don't want it, it is written to the calibration output anyhow. -# ------------------------------------------------------------------------- -#MJPedestalC1.PathOut: . -#MJPedestalC2.PathOut: . - -# ------------------------------------------------------------------------- -# Use this if you want to modify the exclusion of pixels from the beginning -# (default: Pixel 0 and 559) -# ------------------------------------------------------------------------- -MJPedestalC1.BadPixelsFile: mjobs/badpixels_0_559.rc - -# ------------------------------------------------------------------------- -# Use this if you want to modify files containing information about the -# reference lines in the status display -# ------------------------------------------------------------------------- -MJPedestalC1.ReferenceFile: mjobs/pedestalref_Dec04.rc -MJPedestalC2.ReferenceFile: mjobs/pedestalref_Dec04.rc -MJCalibration.ReferenceFile: mjobs/calibrationref_Dec04.rc -MJCalibration.MHCalibrationRelTimeCam.ReferenceFile: mjobs/calibrationref_Dec04.rc -MJCalibration.MHCalibrationChargeCam.ReferenceFile: mjobs/calibrationref_Dec04.rc - -# ------------------------------------------------------------------------- -# Use this to define where the calibration output is stored. The filename -# is created from the sequence number. If nothing is specified '.' is -# assumed. (overwrites command line options) -# ------------------------------------------------------------------------- -#MJCalibration.PathOut: calped - -# ------------------------------------------------------------------------- -# Use this to define where the program should search for the pedestal -# and calibration files defined in the sequence. To use the local -# directory use '.' If nothing is specified the default path in the -# datacenter is used. (overwrites command line options) -# ------------------------------------------------------------------------- -#MJPedestalC1.PathData: /data/MAGIC/Period016/ -#MJPedestalC2.PathData: /data/MAGIC/Period016/ -#MJCalibration.PathData: /data/MAGIC/Period016/ - -# ------------------------------------------------------------------------- -# Define here which signal extractor you are going to use for the -# whole calibration chain -# ------------------------------------------------------------------------- -MJPedestalC1.ExtractSignal: MExtractTimeAndChargeSpline -# ------------------------------------------------------------------------- -# Possibilities for a different extractor -# ------------------------------------------------------------------------- -#MJPedestalC1.ExtractSignal: MExtractTimeAndChargeDigitalFilter, MExtractTimeAndChargeDigitalFilterPeakSearch, MExtractTimeAndChargeSlidingWindow -# ------------------------------------------------------------------------- -# Define here parameters valid for all extractors above -# The actual defaults used depend on the extractor in use. -# For details see the corresponding source code -# ------------------------------------------------------------------------- -#MJPedestalC1.ExtractSignal.HiGainFirst: 0 -#MJPedestalC1.ExtractSignal.HiGainLast: 15 -#MJPedestalC1.ExtractSignal.LoGainFirst: 3 -#MJPedestalC1.ExtractSignal.LoGainLast: 14 -#MJPedestalC1.ExtractSignal.SaturationLimit: 250 -#MJPedestalC1.ExtractSignal.OffsetLoGain: 1.3 -#MJPedestalC1.ExtractSignal.LoGainSwitch: 120 -#MJPedestalC1.ExtractSignal.LoGainStartShift: -2.5 -# ------------------------------------------------------------------------- -# Define here parameters valid for the spline: -# ------------------------------------------------------------------------- -#MJPedestalC1.ExtractSignal.Resolution: 0.05 -MJPedestalC1.ExtractSignal.RiseTimeHiGain: 0.5 -MJPedestalC1.ExtractSignal.FallTimeHiGain: 0.5 -#MJPedestalC1.ExtractSignal.LoGainStretch: 1.5 -#MJPedestalC1.ExtractSignal.ExtractionType: Integral -#MJPedestalC1.ExtractSignal.ExtractionType: Amplitude - -# ------------------------------------------------------------------------- -# In the case of the Digital filter you have to overwrite the default -# weights file depending on what you are extracting -# ------------------------------------------------------------------------- -#MJPedestalC2.ExtractSignal.WeightsFile: msignal/calibration_weights_UV4.dat -#MJPedestalY2.ExtractSignal.WeightsFile: msignal/cosmics_weights4.dat -#MJPedestalY3.ExtractSignal.WeightsFile: msignal/cosmics_weights4.dat -#MJCalibration.ExtractSignal.WeightsFile: msignal/calibration_weights_UV4.dat -#MJCalibrateSignal.ExtractSignal.WeightsFile: msignal/cosmics_weights4.dat -#MJCalibrateSignal.ExtractInterlaced.WeightsFile: msignal/calibration_weights_UV4.dat -#MJCalibrateSignal.ExtractSignal.LoGainStretch: 1.5 -#MJCalibrateSignal.ExtractInterlaced.LoGainStretch: 1.5 - -# ------------------------------------------------------------------------- -# Configure MJCalibration -# ------------------------------------------------------------------------- - -# Switch on relative time calibration -#MJCalibration.RelTimeCalibration: Yes -# Set color to be used -#MJCalibration.Color: -# Type of used data format: raw, root, mc -#MJCalibration.DataType: Root -# Type of displayed plots: Full,DataCheck,Normal -#MJCalibration.Display: datacheck -# Write additinal debug output -#MJCalibration.Debug: No -# Use blind pixel -#MJCalibration.UseBlindPixel: No -# Use pin diode -#MJCalibration.UsePINDiode: No -# Use pixel check option -#MJCalibration.PixelCheck: No -#MJCalibTest.PixelCheck: No - -# ------------------------------------------------------------------------- -# Individual events check for calibration results -# ------------------------------------------------------------------------- - -#ContCosmics.MaxEmptyPixels: 0.2 -#ContCosmics.MaxExcludedFraction: 1 -#ContCosmics.MinAcceptedFraction: 0 - -# ------------------------------------------------------------------------- -# In case of MOON DATA, uncomment the following lines: -# (better: call callisto with option "-moon") -# ------------------------------------------------------------------------- -#MJPedestalC1.UseHists: yes -#MJPedestalC2.UseHists: yes -#MJPedestalY1.UseHists: yes -#MJPedestalY2.UseHists: yes -#MJPedestalY3.UseHists: yes - -# ------------------------------------------------------------------------- -# You can choose the pedestal extraction algorithm/task. To use -# MPedCalcPedRun use the lines below. Be carefull, a pedestal file -# could be (if no pedestal file available) a data-file, too. -# ------------------------------------------------------------------------- -#MJPedestalC1.ExtractPedestal: MPedCalcPedRun -#MJPedestalC1.ExtractPedestal.ExtractWinFirst: 0 -#MJPedestalC1.ExtractPedestal.ExtractWinSize: 6 -#MJPedestalC1.ExtractPedestal.PedestalUpdate: no -#MJPedestalC1.ExtractPedestal.CheckWinFirst: 0 -#MJPedestalC1.ExtractPedestal.CheckWinLast: 29 -#MJPedestalC1.ExtractPedestal.MaxSignalVar: 40 - -#MJPedestalC2.ExtractPedestal: MPedCalcPedRun -#MJPedestalC2.ExtractPedestal.ExtractWinFirst: 0 -#MJPedestalC2.ExtractPedestal.ExtractWinSize: 6 -#MJPedestalC2.ExtractPedestal.PedestalUpdate: no -#MJPedestalC2.ExtractPedestal.CheckWinFirst: 0 -#MJPedestalC2.ExtractPedestal.CheckWinLast: 29 -#MJPedestalC2.ExtractPedestal.MaxSignalVar: 40 - -# ------------------------------------------------------------------------- -# You can choose the pedestal extraction algorithm/task. To use -# MPedCalcFromLoGain use the lines below. Be carefull, a pedestal file -# could be (if no pedestal file available) a data-file, too. -# ------------------------------------------------------------------------- -#MJPedestalC1.ExtractPedestal: MPedCalcFromLoGain -#MJPedestalC1.ExtractPedestal.ExtractWinFirst: 0 -#MJPedestalC1.ExtractPedestal.ExtractWinSize: 6 -#MJPedestalC1.ExtractPedestal.PedestalUpdate: no -#MJPedestalC1.ExtractPedestal.CheckWinFirst: 0 -#MJPedestalC1.ExtractPedestal.CheckWinLast: 29 -#MJPedestalC1.ExtractPedestal.MaxSignalVar: 40 - -#MJPedestalC2.ExtractPedestal: MPedCalcFromLoGain -#MJPedestalC2.ExtractPedestal.ExtractWinFirst: 0 -#MJPedestalC2.ExtractPedestal.ExtractWinSize: 6 -#MJPedestalC2.ExtractPedestal.PedestalUpdate: no -#MJPedestalC2.ExtractPedestal.CheckWinFirst: 0 -#MJPedestalC2.ExtractPedestal.CheckWinLast: 29 -#MJPedestalC2.ExtractPedestal.MaxSignalVar: 40 - -# ------------------------------------------------------------------------- -# Configure MJPedestalC1 and MJPedestalC2 -# ------------------------------------------------------------------------- -# Possible data type: Raw,Root,MC -#MJPedestalC1.DataType: Root -#MJPedestalC2.DataType: Root -#MJPedestalC3.DataType: Root -# Type of displayed plots: Full,DataCheck,None -#MJPedestalC1.Display: datacheck -#MJPedestalC2.Display: datacheck -#MJPedestalC3.Display: datacheck -# Maximum number of event processed in the loop -#MJPedestalC1.MaxEvents: 1000 -#MJPedestalC2.MaxEvents: 1000 -# Allow to overwrite existing files with the output file -#MJPedestalC1.AllowOverwrite: No -#MJPedestalC2.AllowOverwrite: No -# Use data runs from the sequence instead of calibration runs -#MJPedestalC1.UseData: No -#MJPedestalC2.UseData: No - - -# ------------------------------------------------------------------------- -# Use this if you want to change the software low-gain switch -# ------------------------------------------------------------------------- -#MJCalibration.ExtractSignal.LoGainSwitch: 150 - -# ------------------------------------------------------------------------- -# Use this if you want to change the higain-vs. logain intercalibration -# ------------------------------------------------------------------------- -MJCalibration.HiLoCalibFile: resources/hilocalib_sp1.rc - -# ------------------------------------------------------------------------- -# Use this to change the behaviour of the calibration -# ------------------------------------------------------------------------- -# Type if you set a colour explicitely from outside (only for MC!!!) -#MJCalibration.MCalibColorSet.ExplicitColor: green,blue,uv,ct1 - -#MJCalibration.MCalibrationChargeCalc.ArrTimeRmsLimit: 5.0 -#MJCalibration.MCalibrationChargeCalc.ChargeLimit: 2.5 -#MJCalibration.MCalibrationChargeCalc.ChargeErrLimit: 0 -#MJCalibration.MCalibrationChargeCalc.ChargeRelErrLimit: 1 -#MJCalibration.MCalibrationChargeCalc.Debug: no - -#MJCalibration.MCalibrationChargeCalc.FFactorErrLimit: 4.5 -#MJCalibration.MCalibrationChargeCalc.LambdaErrLimit: 0.2 -#MJCalibration.MCalibrationChargeCalc.LambdaCheckLimit: 0.5 -#MJCalibration.MCalibrationChargeCalc.PheErrLowerLimit: 9.0 -#MJCalibration.MCalibrationChargeCalc.PheErrUpperLimit: 5.5 - -#MJCalibration.MHCalibrationChargeCam.Debug: no -#MJCalibration.MHCalibrationChargeCam.LoGain: yes -#MJCalibration.MHCalibrationChargeCam.Oscillations: yes -#MJCalibration.MHCalibrationChargeCam.SizeCheck: yes -#MJCalibration.MHCalibrationChargeCam.Averageing: yes -#MJCalibration.MHCalibrationChargeCam.HiGainNbins: 500 -#MJCalibration.MHCalibrationChargeCam.HiGainFirst: -100.125 -#MJCalibration.MHCalibrationChargeCam.HiGainLast: 1899.875 -#MJCalibration.MHCalibrationChargeCam.LoGainNbins: 500 -#MJCalibration.MHCalibrationChargeCam.LoGainFirst: -100.25 -#MJCalibration.MHCalibrationChargeCam.LoGainLast: 899.75 -#MJCalibration.MHCalibrationChargeCam.TimeLowerLimit: 1. -#MJCalibration.MHCalibrationChargeCam.TimeUpperLimit: 3. -#MJCalibration.MHCalibrationChargeCam.NumHiGainSaturationLimit: 0.02 -#MJCalibration.MHCalibrationChargeCam.NumLoGainSaturationLimit: 0.005 -MJCalibration.MHCalibrationChargeCam.ProbLimit: 1.0E-18 -#MJCalibration.MHCalibrationChargeCam.OverflowLimit: 0.005 -#MJCalibration.MHCalibrationChargeCam.PulserFrequency: 500 - -#MJCalibration.MHCalibrationRelTimeCam.Debug: no -#MJCalibration.MHCalibrationRelTimeCam.LoGain: no -#MJCalibration.MHCalibrationRelTimeCam.Oscillations: yes -#MJCalibration.MHCalibrationRelTimeCam.SizeCheck: yes -#MJCalibration.MHCalibrationRelTimeCam.Averageing: yes -#MJCalibration.MHCalibrationRelTimeCam.Nbins: 300 -#MJCalibration.MHCalibrationRelTimeCam.First: -4.975 -#MJCalibration.MHCalibrationRelTimeCam.Last: 10.025 -#MJCalibration.MHCalibrationRelTimeCam.NumHiGainSaturationLimit: 0.25 -#MJCalibration.MHCalibrationRelTimeCam.ProbLimit: 0.00001 -#MJCalibration.MHCalibrationRelTimeCam.OverflowLimit: 0.005 -#MJCalibration.MHCalibrationRelTimeCam.PulserFrequency: 500 - - -# ========================================================================== -############################################################################# -# ========================================================================== -# Signal extraction -# ========================================================================== -############################################################################# -# ========================================================================== - -# ------------------------------------------------------------------------- -# Use a OutputPath if you want to write the results to a file -# ------------------------------------------------------------------------- -#MJCalibrateSignal.PathOut: . -#MJPedestalY1.PathOut: . -#MJPedestalY2.PathOut: . -#MJPedestalY3.PathOut: . -#MJPedestalY1.PathData: /data/MAGIC/Period016/ -#MJPedestalY2.PathData: /data/MAGIC/Period016/ -#MJPedestalY3.PathData: /data/MAGIC/Period016/ -#MJCalibrateSignal.PathData: /data/MAGIC/Period016/ -#MJCalibrateSignal.PathIn: . - -# ------------------------------------------------------------------------- -# If you want to change the behaviour of the pedestal extraction... -# ------------------------------------------------------------------------- -# Possible data type: Raw,Root,MC -#MJPedestalY1.DataType: Root -#MJPedestalY2.DataType: Root -#MJPedestalY3.DataType: Root -#MJPedestalY1.UseData: Yes -#MJPedestalY2.UseData: Yes -#MJPedestalY3.UseData: Yes -#MJPedestalY1.MinEvents: 50 -#MJPedestalY2.MinEvents: 50 -#MJPedestalY3.MinEvents: 50 - -MJPedestalY1.MaxEvents: 1000 -MJPedestalY2.MaxEvents: 10000 -MJPedestalY3.MaxEvents: 1000 - -MJPedestalY1.MaxPedestals: 500 -MJPedestalY2.MaxPedestals: 500 -MJPedestalY3.MaxPedestals: 500 - -# ------------------------------------------------------------------------- -# When the starting pedestal is extracted from the first events of -# the data run also the average pulse position of these events -# is determined. If the pulse position check is switched on (default) -# the extraction window of the hi- and lo-gain extractor for the -# determined average pulse position is not allowed to exceed the -# extractors search region. In this scenario we assume that the -# start of the extraction window is at least similar to the -# rising edge of the pulse. -# To make sure that the search region and the extraction window -# does not touch increase the values. -# ------------------------------------------------------------------------- -#MJPedestalY2.PulsePosCheck: yes -#MJPedestalY2.ExtractWinLeft: 0.0 -#MJPedestalY2.ExtractWinRight: 0.0 - -# ------------------------------------------------------------------------- -# Define the Pulse Position check parameters: -# ------------------------------------------------------------------------- -#MJPedestalY2.MHCalibrationPulseTimeCam.SaturationLimit: 245 -#MJPedestalY2.MHCalibrationPulseTimeCam.LowerSignalLimit: 85 -#MJPedestalY2.MHCalibrationPulseTimeCam.NumPixelsRequired: 2 -#MJPedestalY2.PixelCheck: no - -#MJPedestalY.ExtractPedestal: MPedCalcFromLoGain -#MJPedestalY.ExtractPedestal.PedestalUpdate: no - -# ------------------------------------------------------------------------- -# Setup pedestal extractor for running through data -# ------------------------------------------------------------------------- -#MJCalibrateSignal.MPedCalcFundamental.PedestalUpdate: yes -#MJCalibrateSignal.MPedCalcFundamental.NumDump: 500 -#MJCalibrateSignal.MPedCalcFundamental.NumEventsDump: 500 -#MJCalibrateSignal.MPedCalcFundamental.NumAreasDump: 500 -#MJCalibrateSignal.MPedCalcFundamental.NumSectorsDump: 500 -#MJCalibrateSignal.MPedCalcFundamental.ExtractWinFirst: 17 -#MJCalibrateSignal.MPedCalcFundamental.ExtractWinSize: 6 -#MJCalibrateSignal.MPedCalcFundamental.CheckWinFirst: 0 -#MJCalibrateSignal.MPedCalcFundamental.CheckWinLast: 29 -#MJCalibrateSignal.MPedCalcFundamental.MaxSignalVar: 40 - -#MJCalibrateSignal.MPedCalcFromExtractorRndm.PedestalUpdate: yes -#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumEventsDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumAreasDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumSectorsDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.ExtractWinFirst: 17 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.ExtractWinSize: 6 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.CheckWinFirst: 0 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.CheckWinLast: 29 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.MaxSignalVar: 40 - -#MJCalibrateSignal.MPedCalcFromExtractor.PedestalUpdate: yes -#MJCalibrateSignal.MPedCalcFromExtractor.NumDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractor.NumEventsDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractor.NumAreasDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractor.NumSectorsDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractor.ExtractWinFirst: 17 -#MJCalibrateSignal.MPedCalcFromExtractor.ExtractWinSize: 6 -#MJCalibrateSignal.MPedCalcFromExtractor.CheckWinFirst: 0 -#MJCalibrateSignal.MPedCalcFromExtractor.CheckWinLast: 29 -#MJCalibrateSignal.MPedCalcFromExtractor.MaxSignalVar: 40 - -# ------------------------------------------------------------------------- -# Configure the interlaced calibration updates -# ------------------------------------------------------------------------- -#MJCalibrateSignal.MCalibCalcFromPast.UpdateWithFFactorMethod: yes -#MJCalibrateSignal.MCalibCalcFromPast.NumEventsDump: 500 -#MJCalibrateSignal.MCalibCalcFromPast.UpdateNumPhes: yes -#MJCalibrateSignal.MCalibCalcFromPast.NumPhesDump: 10 -#MJCalibrateSignal.MCalibrationChargeCalc.Debug: no - -# ------------------------------------------------------------------------- -# Use this if you want to change the signal extractor -# ------------------------------------------------------------------------- -MJCalibrateSignal.ExtractSignal.LoGainSwitch: 150 - -# ------------------------------------------------------------------------- -# Setup calibration of data -# ------------------------------------------------------------------------- -# Type of used data format: raw,root,MC -#MJCalibrateSignal.DataType: Root -# Type if you set a colour explicitely from outside (only for MC!!!) -#MJCalibrateSignal.MCalibColorSet.ExpicitColor: green,blue,uv,ct1 -#MJCalibrateSignal.MCalibrateData.PedestalFlag: Event -#MJCalibrateSignal.MCalibrateData.ScaleFactor: 0.85182358 -MJCalibrateSignal.MCalibrateData.FileNameScale: resources/calibration_spline.rc -#MJCalibrateSignal.MCalibrateData.CalibrationMode: Default -#MJCalibrateSignal.MCalibrateData.CalibConvMinLimit: 0.01 -#MJCalibrateSignal.MCalibrateData.CalibConvMaxLimit: 5. -#MJCalibrateSignal.Interlaced: yes -#MJCalibrateSignal.RelTimesUpdate: yes - -# ------------------------------------------------------------------------- -# Setup level for determining of bad pixels -# ------------------------------------------------------------------------- -#MJCalibrateSignal.MBadPixelsCalc.PedestalLevel: 3.0 -#MJCalibrateSignal.MBadPixelsCalc.PedestalLevelVariance: 3.0 -#MJCalibrateSignal.MBadPixelsTreat.NumMinNeighbors: 3 -#MJCalibrateSignal.MBadPixelsTreat.MaxArrivalTimeDiff: 3.0 -#MJCalibrateSignal.MBadPixelsTreat.UseInterpolation: yes -#MJCalibrateSignal.MBadPixelsTreat.ProcessPedestalEvt: yes -#MJCalibrateSignal.MBadPixelsTreat.ProcessPedestalRun: no -#MJCalibrateSignal.MBadPixelsTreat.ProcessTimes: yes -#MJCalibrateSignal.MBadPixelsTreat.UseCentralPixel: no -#MJCalibrateSignal.MBadPixelsTreat.HardTreatment: no - -# ------------------------------------------------------------------------- -# Setup interlaced calibration events histogramming -# ------------------------------------------------------------------------- -#MJCalibrateSignal.MHCalibrationChargeCam.Debug: no -#MJCalibrateSignal.MHCalibrationChargeCam.LoGain: yes -MJCalibrateSignal.MHCalibrationChargeCam.Oscillations: no -#MJCalibrateSignal.MHCalibrationChargeCam.SizeCheck: yes -#MJCalibrateSignal.MHCalibrationChargeCam.Averageing: yes -#MJCalibrateSignal.MHCalibrationChargeCam.HiGainNbins: 500 -#MJCalibrateSignal.MHCalibrationChargeCam.HiGainFirst: -100.5 -#MJCalibrateSignal.MHCalibrationChargeCam.HiGainLast: 1899.5 -MJCalibrateSignal.MHCalibrationChargeCam.LoGainNbins: 250 -MJCalibrateSignal.MHCalibrationChargeCam.LoGainFirst: -98 -MJCalibrateSignal.MHCalibrationChargeCam.LoGainLast: 902 -#MJCalibrateSignal.MHCalibrationChargeCam.TimeLowerLimit: 1. -#MJCalibrateSignal.MHCalibrationChargeCam.TimeUpperLimit: 3. -#MJCalibrateSignal.MHCalibrationChargeCam.NumHiGainSaturationLimit: 0.02 -#MJCalibrateSignal.MHCalibrationChargeCam.NumLoGainSaturationLimit: 0.005 -#MJCalibrateSignal.MHCalibrationChargeCam.ProbLimit: 0.0001 -#MJCalibrateSignal.MHCalibrationChargeCam.OverflowLimit: 0.005 -MJCalibrateSignal.MHCalibrationChargeCam.PulserFrequency: 50 - -# ------------------------------------------------------------------------- -# Setup relative time interlaced calibration -# ------------------------------------------------------------------------- -#MJCalibrateSignal.MHCalibrationRelTimeCam.Debug: no -#MJCalibrateSignal.MHCalibrationRelTimeCam.LoGain: no -MJCalibrateSignal.MHCalibrationRelTimeCam.Oscillations: no -#MJCalibrateSignal.MHCalibrationRelTimeCam.SizeCheck: yes -#MJCalibrateSignal.MHCalibrationRelTimeCam.Averageing: yes -#MJCalibrateSignal.MHCalibrationRelTimeCam.Nbins: 300 -#MJCalibrateSignal.MHCalibrationRelTimeCam.First: -4.975 -#MJCalibrateSignal.MHCalibrationRelTimeCam.Last: 10.025 -#MJCalibrateSignal.MHCalibrationRelTimeCam.NumHiGainSaturationLimit: 0.25 -#MJCalibrateSignal.MHCalibrationRelTimeCam.ProbLimit: 0.00001 -#MJCalibrateSignal.MHCalibrationRelTimeCam.OverflowLimit: 0.005 -MJCalibrateSignal.MHCalibrationRelTimeCam.PulserFrequency: 50 - -# ------------------------------------------------------------------------- -# Resource for movie encoding. For more details see IACT wiki -# ------------------------------------------------------------------------- -#MMoviePrepare.FirstSlice: 0 -#MMoviePrepare.LastSlice: 15 -#MovieFilter.Condition: MMovieData.fMax>5*MMovieData.fMedianPedestalRms -#MImgCleanStd.CleanLevel1: 8.0 -#MImgCleanStd.CleanLevel2: 4.5 -#MFEvtNumber.FileName: ganymed00001218.root -#MFEvtNumber.Selector: ThetaSquared.fVal<0.04 && DataType.fVal>0.5 -#MMovieWrite.TargetLength: 5 -#MMovieWrite.NumEvents: 500 -#MMovieWrite.Threshold: 3 -#MMovieWrite.FileName: filename.mpg - -# ------------------------------------------------------------------------- -# This is a special option for the datacheck. A task executed as last -# task in the event processing task list - in the standard software it -# is skipped -# ------------------------------------------------------------------------- -#MJCalibrateSignal.FinalFantasy: MHCamFlorian Index: trunk/Mars/callisto_MarApr05.rc =================================================================== --- trunk/Mars/callisto_MarApr05.rc (revision 19871) +++ (revision ) @@ -1,485 +1,0 @@ -# ========================================================================== -############################################################################# -# ========================================================================== -# General -# ========================================================================== -############################################################################# -# ========================================================================== - -# ------------------------------------------------------------------------- -# Use this if you want to setup the logging stream for the jobs -# (overwrites command line options) -# ------------------------------------------------------------------------- -#MLog.VerbosityLevel: 2 -#MLog.DebugLevel: 1 -#MLog.NoColors: yes - -# ========================================================================== -############################################################################# -# ========================================================================== -# Calibration -# ========================================================================== -############################################################################# -# ========================================================================== - -# ------------------------------------------------------------------------- -# Use this if you want to write the MJPedestal output somewhere -# If you don't want it, it is written to the calibration output anyhow. -# ------------------------------------------------------------------------- -#MJPedestalC1.PathOut: . -#MJPedestalC2.PathOut: . -#MJPedestalC1.BadPixelsFile: mjobs/badpixels_0_559.rc -#MJPedestalC1.ReferenceFile: mjobs/pedestalref.rc -#MJPedestalC2.ReferenceFile: mjobs/pedestalref.rc -MJCalibration.ReferenceFile: mjobs/calibrationref_Nov04.rc -MJCalibration.MHCalibrationRelTimeCam.ReferenceFile: mjobs/calibrationref_Nov04.rc -MJCalibration.MHCalibrationChargeCam.ReferenceFile: mjobs/calibrationref_Nov04.rc - -# ------------------------------------------------------------------------- -# Use this to define where the calibration output is stored. The filename -# is created from the sequence number. If nothing is specified '.' is -# assumed. (overwrites command line options) -# ------------------------------------------------------------------------- -#MJCalibration.PathOut: calped - -# ------------------------------------------------------------------------- -# Use this to define where the program should search for the pedestal -# and calibration files defined in the sequence. To use the local -# directory use '.' If nothing is specified the default path in the -# datacenter is used. (overwrites command line options) -# ------------------------------------------------------------------------- -#MJPedestalC1.PathData: /data/MAGIC/Period016/ -#MJPedestalC2.PathData: /data/MAGIC/Period016/ -#MJCalibration.PathData: /data/MAGIC/Period016/ - -# ------------------------------------------------------------------------- -# Define here which signal extractor you are going to use for the -# whole calibration chain -# ------------------------------------------------------------------------- -MJPedestalC1.ExtractSignal: MExtractTimeAndChargeDigitalFilter -# ------------------------------------------------------------------------- -# Define here an extractor which can be used for the December 04 data. -# ------------------------------------------------------------------------- -#MJPedestalC1.ExtractSignal: MExtractTimeAndChargeSpline, MExtractTimeAndChargeDigitalFilterPeakSearch, MExtractTimeAndChargeSlidingWindow -# ------------------------------------------------------------------------- -# Define here parameters valid for all extractors above -# ------------------------------------------------------------------------- -#MJPedestalC1.ExtractSignal.HiGainFirst: 0 -#MJPedestalC1.ExtractSignal.HiGainLast: 15 -#MJPedestalC1.ExtractSignal.LoGainFirst: 3 -#MJPedestalC1.ExtractSignal.LoGainLast: 14 -#MJPedestalC1.ExtractSignal.SaturationLimit: 250 -#MJPedestalC1.ExtractSignal.OffsetLoGain: 1.51 -#MJPedestalC1.ExtractSignal.LoGainSwitch: 150 -#MJPedestalC1.ExtractSignal.LoGainStartShift: -2.8 -# ------------------------------------------------------------------------- -# Define here parameters valid for the spline: -# ------------------------------------------------------------------------- -#MJPedestalC1.ExtractSignal.Resolution: 0.05 -#MJPedestalC1.ExtractSignal.RiseTimeHiGain: 0.5 -#MJPedestalC1.ExtractSignal.FallTimeHiGain: 1.5 -#MJPedestalC1.ExtractSignal.LoGainStretch: 1.5 -#MJPedestalC1.ExtractSignal.ExtractionType: Integral -#MJPedestalC1.ExtractSignal.ExtractionType: Amplitude -# ------------------------------------------------------------------------- -# Define here parameters valid for the sliding window: -# ------------------------------------------------------------------------- -#MJPedestalC1.ExtractSignal.HiGainWindowSize: 6 -#MJPedestalC1.ExtractSignal.LoGainWindowSize: 6 -# ------------------------------------------------------------------------- -# Define here parameters valid for the digital filter with peak search: -# ------------------------------------------------------------------------- -#MJPedestalC1.ExtractSignal.OffsetLeftFromPeak: 3 -#MJPedestalC1.ExtractSignal.OffsetRightFromPeak: 3 -#MJPedestalC1.ExtractSignal.PeakSearchWindowSize: 2 -#MJPedestalC1.ExtractSignal.HiGainFailureLimit: 10 -#MJPedestalC1.ExtractSignal.LoGainFailureLimit: 25 - -# ------------------------------------------------------------------------- -# In the case of the Digital filter you have to overwrite the default -# weights file depending on what you are extracting -# ------------------------------------------------------------------------- -MJPedestalC2.ExtractSignal.WeightsFile: msignal/calibration_weights_UV46.dat -MJPedestalY2.ExtractSignal.WeightsFile: msignal/cosmics_weights46.dat -MJPedestalY3.ExtractSignal.WeightsFile: msignal/cosmics_weights46.dat -MJCalibration.ExtractSignal.WeightsFile: msignal/calibration_weights_UV46.dat -MJCalibrateSignal.ExtractSignal.WeightsFile: msignal/cosmics_weights46.dat -MJCalibrateSignal.ExtractInterlaced.WeightsFile: msignal/calibration_weights_UV46.dat - -#MJPedestalC2.ExtractSignal.AutomaticWeightsFile: On -#MJPedestalY2.ExtractSignal.AutomaticWeightsFile: On -#MJPedestalY3.ExtractSignal.AutomaticWeightsFile: On -#MJCalibration.ExtractSignal.AutomaticWeightsFile: On -#MJCalibrateSignal.ExtractSignal.AutomaticWeightsFile: On -#MJCalibrateSignal.ExtractInterlaced.AutomaticWeightsFile: On - -# ------------------------------------------------------------------------- -# Configure MJCalibration -# ------------------------------------------------------------------------- - -# Switch on relative time calibration -#MJCalibration.RelTimeCalibration: Yes -# Type of used data format: raw, root, mc -#MJCalibration.DataType: Root -# Type of displayed plots: Full,DataCheck,Normal -#MJCalibration.Display: datacheck -# Write additinal debug output -#MJCalibration.Debug: No -# Use blind pixel -#MJCalibration.UseBlindPixel: No -# Use pin diode -#MJCalibration.UsePINDiode: No -# Use pixel check option -#MJCalibration.PixelCheck: Yes -#MJCalibTest.PixelCheck: Yes - -# ------------------------------------------------------------------------- -# Individual events check for calibration results -# ------------------------------------------------------------------------- - -#ContCosmics.MaxEmptyPixels: 0.2 -#ContCosmics.MaxExcludedFraction: 1 -#ContCosmics.MinAcceptedFraction: 0 - -# ------------------------------------------------------------------------- -# You can choose the pedestal extraction algorithm/task. To use -# MPedCalcPedRun use the lines below. Be carefull, a pedestal file -# could be (if no pedestal file available) a data-file, too. -# ------------------------------------------------------------------------- -#MJPedestalC1.ExtractPedestal: MPedCalcPedRun -#MJPedestalC1.ExtractPedestal.ExtractWinFirst: 0 -#MJPedestalC1.ExtractPedestal.ExtractWinSize: 6 -#MJPedestalC1.ExtractPedestal.PedestalUpdate: no -#MJPedestalC1.ExtractPedestal.CheckWinFirst: 0 -#MJPedestalC1.ExtractPedestal.CheckWinSize: 29 - -#MJPedestalC2.ExtractPedestal: MPedCalcPedRun -#MJPedestalC2.ExtractPedestal.ExtractWinFirst: 0 -#MJPedestalC2.ExtractPedestal.ExtractWinSize: 6 -#MJPedestalC2.ExtractPedestal.PedestalUpdate: no -#MJPedestalC2.ExtractPedestal.CheckWinFirst: 0 -#MJPedestalC2.ExtractPedestal.CheckWinSize: 29 - -# ------------------------------------------------------------------------- -# You can choose the pedestal extraction algorithm/task. To use -# MPedCalcFromLoGain use the lines below. Be carefull, a pedestal file -# could be (if no pedestal file available) a data-file, too. -# ------------------------------------------------------------------------- -#MJPedestalC1.ExtractPedestal: MPedCalcFromLoGain -#MJPedestalC1.ExtractPedestal.ExtractWinFirst: 0 -#MJPedestalC1.ExtractPedestal.ExtractWinSize: 6 -#MJPedestalC1.ExtractPedestal.PedestalUpdate: no -#MJPedestalC1.ExtractPedestal.CheckWinFirst: 0 -#MJPedestalC1.ExtractPedestal.CheckWinSize: 29 -#MJPedestalC1.ExtractPedestal.MaxSignalVar: 40 - -#MJPedestalC2.ExtractPedestal: MPedCalcFromLoGain -#MJPedestalC2.ExtractPedestal.ExtractWinFirst: 0 -#MJPedestalC2.ExtractPedestal.ExtractWinSize: 6 -#MJPedestalC2.ExtractPedestal.PedestalUpdate: no -#MJPedestalC2.ExtractPedestal.CheckWinFirst: 0 -#MJPedestalC2.ExtractPedestal.CheckWinSize: 29 -#MJPedestalC2.ExtractPedestal.MaxSignalVar: 40 - -# ------------------------------------------------------------------------- -# Configure MJPedestalC1 and MJPedestalC2 -# ------------------------------------------------------------------------- -# Possible data type: Raw,Root,MC -#MJPedestalC1.DataType: Root -#MJPedestalC2.DataType: Root -#MJPedestalC3.DataType: Root -# Type of displayed plots: Full,DataCheck,None -#MJPedestalC1.Display: datacheck -#MJPedestalC2.Display: datacheck -#MJPedestalC3.Display: datacheck -# Maximum number of event processed in the loop -#MJPedestalC1.MaxEvents: 1000 -#MJPedestalC2.MaxEvents: 1000 -# Allow to overwrite existing files with the output file -#MJPedestalC1.AllowOverwrite: No -#MJPedestalC2.AllowOverwrite: No -# Use data runs from the sequence instead of calibration runs -#MJPedestalC1.UseData: No -#MJPedestalC2.UseData: No -# Switch off DeadPixelCheck with "No" -#MJPedestalY2.DeadPixelCheck: yes - - -# ------------------------------------------------------------------------- -# Use this if you want to change the software low-gain switch in the -# calibration -# ------------------------------------------------------------------------- -MJCalibration.ExtractSignal.LoGainSwitch: 170 - -# ------------------------------------------------------------------------- -# Use this if you want to change the higain-vs. logain intercalibration -# ------------------------------------------------------------------------- -#MJCalibration.HiLoCalibFile: mjobs/hilocalib_df46.root -#MJCalibration.HiLoCalibration: yes - -# ------------------------------------------------------------------------- -# Use this if you want to change the time extractor for the calibration -# and automatically the data extraction -# ------------------------------------------------------------------------- -#MJCalibration.ExtractTime: -#MJCalibration.ExtractTime: MExtractTimeFastSpline -#MJCalibration.ExtractTime.HiGainFirst: 0 -#MJCalibration.ExtractTime.HiGainLast: 7 -#MJCalibration.ExtractTime.LoGainFirst: 3 -#MJCalibration.ExtractTime.LoGainLast: 8 -#MJCalibration.ExtractTime.WindowSizeHiGain: 6 -#MJCalibration.ExtractTime.WindowSizeLoGain: 6 -#MJCalibration.ExtractTime.OffsetLoGain: 1.51 - -# ------------------------------------------------------------------------- -# Use this to change the behaviour of the calibration -# ------------------------------------------------------------------------- -# Type if you set a colour explicitely from outside (TEST purposes!!) -#MJCalibration.MCalibColorSet.ExplicitColor: green,blue,uv,ct1 - -#MJCalibration.MCalibrationChargeCalc.ArrTimeRmsLimit: 5.0 -#MJCalibration.MCalibrationChargeCalc.ChargeLimit: 2.5 -#MJCalibration.MCalibrationChargeCalc.ChargeErrLimit: 0 -#MJCalibration.MCalibrationChargeCalc.ChargeRelErrLimit: 1 -#MJCalibration.MCalibrationChargeCalc.Debug: no - -#MJCalibration.MCalibrationChargeCalc.FFactorErrLimit: 4.5 -#MJCalibration.MCalibrationChargeCalc.LambdaErrLimit: 0.2 -#MJCalibration.MCalibrationChargeCalc.LambdaCheckLimit: 0.5 -#MJCalibration.MCalibrationChargeCalc.PheErrLowerLimit: 9.0 -#MJCalibration.MCalibrationChargeCalc.PheErrUpperLimit: 5.5 - -#MJCalibration.MHCalibrationChargeCam.Debug: no -#MJCalibration.MHCalibrationChargeCam.LoGain: yes -#MJCalibration.MHCalibrationChargeCam.Oscillations: yes -#MJCalibration.MHCalibrationChargeCam.SizeCheck: yes -#MJCalibration.MHCalibrationChargeCam.Averageing: yes -#MJCalibration.MHCalibrationChargeCam.HiGainNbins: 500 -#MJCalibration.MHCalibrationChargeCam.HiGainFirst: -98. -#MJCalibration.MHCalibrationChargeCam.HiGainLast: 1902. -#MJCalibration.MHCalibrationChargeCam.LoGainNbins: 500 -#MJCalibration.MHCalibrationChargeCam.LoGainFirst: -99. -#MJCalibration.MHCalibrationChargeCam.LoGainLast: 901. -#MJCalibration.MHCalibrationChargeCam.TimeLowerLimit: 1. -#MJCalibration.MHCalibrationChargeCam.TimeUpperLimit: 3. -#MJCalibration.MHCalibrationChargeCam.NumHiGainSaturationLimit: 0.02 -#MJCalibration.MHCalibrationChargeCam.NumLoGainSaturationLimit: 0.005 -#MJCalibration.MHCalibrationChargeCam.ProbLimit: 0.00000001 -#MJCalibration.MHCalibrationChargeCam.OverflowLimit: 0.005 -#MJCalibration.MHCalibrationChargeCam.PulserFrequency: 500 - - -#MJCalibration.MHCalibrationRelTimeCam.Debug: no -#MJCalibration.MHCalibrationRelTimeCam.LoGain: no -#MJCalibration.MHCalibrationRelTimeCam.Oscillations: yes -#MJCalibration.MHCalibrationRelTimeCam.SizeCheck: yes -#MJCalibration.MHCalibrationRelTimeCam.Averageing: yes -#MJCalibration.MHCalibrationRelTimeCam.HiGainNbins: 210 -#MJCalibration.MHCalibrationRelTimeCam.HiGainFirst: -5. -#MJCalibration.MHCalibrationRelTimeCam.HiGainLast: 10. -#MJCalibration.MHCalibrationRelTimeCam.NumHiGainSaturationLimit: 0.25 -#MJCalibration.MHCalibrationRelTimeCam.ProbLimit: 0.00001 -#MJCalibration.MHCalibrationRelTimeCam.OverflowLimit: 0.005 -#MJCalibration.MHCalibrationRelTimeCam.PulserFrequency: 500 - - -# ========================================================================== -############################################################################# -# ========================================================================== -# Signal extraction -# ========================================================================== -############################################################################# -# ========================================================================== - -# ------------------------------------------------------------------------- -# Use a OutputPath if you want to write the results to a file -# ------------------------------------------------------------------------- -#MJCalibrateSignal.PathOut: . -#MJPedestalY1.PathOut: . -#MJPedestalY2.PathOut: . -#MJPedestalY3.PathOut: . -#MJPedestalY1.PathData: /data/MAGIC/Period016/ -#MJPedestalY2.PathData: /data/MAGIC/Period016/ -#MJPedestalY3.PathData: /data/MAGIC/Period016/ -#MJCalibrateSignal.PathData: /data/MAGIC/Period016/ -#MJCalibrateSignal.PathIn: . - -# ------------------------------------------------------------------------- -# If you want to change the behaviour of the pedestal extraction... -# ------------------------------------------------------------------------- -# Possible data type: Raw,Root,MC -#MJPedestalY1.DataType: Root -#MJPedestalY2.DataType: Root -#MJPedestalY3.DataType: Root -#MJPedestalY1.UseData: Yes -#MJPedestalY2.UseData: Yes -#MJPedestalY3.UseData: Yes - -MJPedestalY1.MaxEvents: 1000 -MJPedestalY2.MaxEvents: 10000 -MJPedestalY3.MaxEvents: 1000 - -MJPedestalY1.MaxPedestals: 500 -MJPedestalY2.MaxPedestals: 500 -MJPedestalY3.MaxPedestals: 500 - -# ------------------------------------------------------------------------- -# Use Pulse Position check to define the extraction ranges for the data? -# ------------------------------------------------------------------------- -#MJPedestalY2.PulsePosCheck: yes -#MJPedestalY2.ExtractWinLeft: 0.0 -#MJPedestalY2.ExtractWinRight: 0.5 - -# ------------------------------------------------------------------------- -# Define the Pulse Position check parameters: -# ------------------------------------------------------------------------- -#MJPedestalY2.MHCalibrationPulseTimeCam.SaturationLimit: 255 -#MJPedestalY2.MHCalibrationPulseTimeCam.LowerSignalLimit: 100 -#MJPedestalY2.MHCalibrationPulseTimeCam.NumPixelsRequired: 2 -#MJPedestalY2.PixelCheck: no - -#MJPedestalY.ExtractPedestal: MPedCalcFromLoGain -#MJPedestalY.ExtractPedestal.PedestalUpdate: no - -# ------------------------------------------------------------------------- -# Setup pedestal extractor for running through data -# ------------------------------------------------------------------------- -#MJCalibrateSignal.MPedCalcFundamental.PedestalUpdate: yes -#MJCalibrateSignal.MPedCalcFundamental.NumDump: 500 -#MJCalibrateSignal.MPedCalcFundamental.NumEventsDump: 500 -#MJCalibrateSignal.MPedCalcFundamental.NumAreasDump: 500 -#MJCalibrateSignal.MPedCalcFundamental.NumSectorsDump: 500 -#MJCalibrateSignal.MPedCalcFundamental.ExtractWinFirst: 17 -#MJCalibrateSignal.MPedCalcFundamental.ExtractWinSize: 6 -#MJCalibrateSignal.MPedCalcFundamental.CheckWinFirst: 0 -#MJCalibrateSignal.MPedCalcFundamental.CheckWinSize: 29 -#MJCalibrateSignal.MPedCalcFundamental.MaxSignalVar: 40 - -#MJCalibrateSignal.MPedCalcFromExtractorRndm.PedestalUpdate: yes -#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumEventsDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumAreasDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumSectorsDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.ExtractWinFirst: 17 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.ExtractWinSize: 6 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.CheckWinFirst: 0 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.CheckWinSize: 29 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.MaxSignalVar: 40 - -#MJCalibrateSignal.MPedCalcFromExtractor.PedestalUpdate: yes -#MJCalibrateSignal.MPedCalcFromExtractor.NumDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractor.NumEventsDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractor.NumAreasDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractor.NumSectorsDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractor.ExtractWinFirst: 17 -#MJCalibrateSignal.MPedCalcFromExtractor.ExtractWinSize: 6 -#MJCalibrateSignal.MPedCalcFromExtractor.CheckWinFirst: 0 -#MJCalibrateSignal.MPedCalcFromExtractor.CheckWinSize: 29 -#MJCalibrateSignal.MPedCalcFromExtractor.MaxSignalVar: 40 - -# ------------------------------------------------------------------------- -# Configure the interlaced calibration updates -# ------------------------------------------------------------------------- -#MJCalibrateSignal.MCalibCalcFromPast.UpdateWithFFactorMethod: yes -#MJCalibrateSignal.MCalibCalcFromPast.NumEventsDump: 500 -#MJCalibrateSignal.MCalibCalcFromPast.UpdateNumPhes: yes -#MJCalibrateSignal.MCalibCalcFromPast.NumPhesDump: 10 -#MJCalibrateSignal.MCalibrationChargeCalc.Debug: no - -# ------------------------------------------------------------------------- -# Use this if you want to change the signal extractor -# ------------------------------------------------------------------------- -MJCalibrateSignal.ExtractSignal.LoGainSwitch: 150 -#MJCalibrateSignal.ExtractSignal.LoGainStretch: 1.5 -#MJCalibrateSignal.ExtractInterlaced.LoGainStretch: 1.5 - -# ------------------------------------------------------------------------- -# Setup calibration of data -# ------------------------------------------------------------------------- -# Type of used data format: raw,root,MC -#MJCalibrateSignal.DataType: Root -# Type if you set a colour explicitely from outside (only TESTS!!!) -#MJCalibrateSignal.MCalibColorSet.ExpicitColor: green,blue,uv,ct1 -#MJCalibrateSignal.MCalibrateData.PedestalFlag: Event -#MJCalibrateSignal.MCalibrateData.ScaleFactor: 1 -MJCalibrateSignal.MCalibrateData.FileNameScale: resources/calibration.rc -#MJCalibrateSignal.MCalibrateData.CalibrationMode: Default -#MJCalibrateSignal.MCalibrateData.CalibConvMinLimit: 0.01 -#MJCalibrateSignal.MCalibrateData.CalibConvMaxLimit: 5. -#MJCalibrateSignal.Interlaced: yes -#MJCalibrateSignal.RelTimesUpdate: yes - -# ------------------------------------------------------------------------- -# Setup level for determining of bad pixels -# ------------------------------------------------------------------------- -#MJCalibrateSignal.MBadPixelsCalc.PedestalLevel: 3.0 -#MJCalibrateSignal.MBadPixelsCalc.PedestalLevelVariance: 3.0 -#MJCalibrateSignal.MBadPixelsTreat.NumMinNeighbors: 3 -#MJCalibrateSignal.MBadPixelsTreat.MaxArrivalTimeDiff: 1.0 -#MJCalibrateSignal.MBadPixelsTreat.UseInterpolation: yes -#MJCalibrateSignal.MBadPixelsTreat.ProcessPedestalEvt: yes -#MJCalibrateSignal.MBadPixelsTreat.ProcessPedestalRun: no -#MJCalibrateSignal.MBadPixelsTreat.ProcessTimes: yes -#MJCalibrateSignal.MBadPixelsTreat.UseCentralPixel: no -#MJCalibrateSignal.MBadPixelsTreat.HardTreatment: no - -# ------------------------------------------------------------------------- -# Setup interlaced calibration events histogramming -# ------------------------------------------------------------------------- -#MJCalibrateSignal.MHCalibrationChargeCam.Debug: no -#MJCalibrateSignal.MHCalibrationChargeCam.LoGain: yes -MJCalibrateSignal.MHCalibrationChargeCam.Oscillations: no -#MJCalibrateSignal.MHCalibrationChargeCam.SizeCheck: yes -#MJCalibrateSignal.MHCalibrationChargeCam.Averageing: yes -#MJCalibrateSignal.MHCalibrationChargeCam.HiGainNbins: 500 -#MJCalibrateSignal.MHCalibrationChargeCam.HiGainFirst: -98. -#MJCalibrateSignal.MHCalibrationChargeCam.HiGainLast: 1902. -#MJCalibrateSignal.MHCalibrationChargeCam.HiGainLast: 1899.5 -MJCalibrateSignal.MHCalibrationChargeCam.LoGainNbins: 250 -MJCalibrateSignal.MHCalibrationChargeCam.LoGainFirst: -98 -MJCalibrateSignal.MHCalibrationChargeCam.LoGainLast: 902 -#MJCalibrateSignal.MHCalibrationChargeCam.TimeLowerLimit: 1. -#MJCalibrateSignal.MHCalibrationChargeCam.TimeUpperLimit: 3. -#MJCalibrateSignal.MHCalibrationChargeCam.NumHiGainSaturationLimit: 0.02 -#MJCalibrateSignal.MHCalibrationChargeCam.NumLoGainSaturationLimit: 0.005 -#MJCalibrateSignal.MHCalibrationChargeCam.ProbLimit: 0.0001 -#MJCalibrateSignal.MHCalibrationChargeCam.OverflowLimit: 0.005 -MJCalibrateSignal.MHCalibrationChargeCam.PulserFrequency: 50 - -# ------------------------------------------------------------------------- -# Setup relative time interlaced calibration -# ------------------------------------------------------------------------- -#MJCalibrateSignal.MHCalibrationRelTimeCam.Debug: no -#MJCalibrateSignal.MHCalibrationRelTimeCam.LoGain: no -MJCalibrateSignal.MHCalibrationRelTimeCam.Oscillations: no -#MJCalibrateSignal.MHCalibrationRelTimeCam.SizeCheck: yes -#MJCalibrateSignal.MHCalibrationRelTimeCam.Averageing: yes -#MJCalibrateSignal.MHCalibrationRelTimeCam.HiGainNbins: 210 -#MJCalibrateSignal.MHCalibrationRelTimeCam.HiGainFirst: -5. -#MJCalibrateSignal.MHCalibrationRelTimeCam.HiGainLast: 10. -#MJCalibrateSignal.MHCalibrationRelTimeCam.NumHiGainSaturationLimit: 0.25 -#MJCalibrateSignal.MHCalibrationRelTimeCam.ProbLimit: 0.00001 -#MJCalibrateSignal.MHCalibrationRelTimeCam.OverflowLimit: 0.005 -MJCalibrateSignal.MHCalibrationRelTimeCam.PulserFrequency: 50 - -# ------------------------------------------------------------------------- -# Resource for movie encoding. For more details see IACT wiki -# ------------------------------------------------------------------------- -#MMoviePrepare.FirstSlice: 0 -#MMoviePrepare.LastSlice: 15 -#MovieFilter.Condition: MMovieData.fMax>5*MMovieData.fMedianPedestalRms -#MImgCleanStd.CleanLevel1: 8.0 -#MImgCleanStd.CleanLevel2: 4.5 -#MFEvtNumber.FileName: ganymed00001218.root -#MFEvtNumber.Selector: ThetaSquared.fVal<0.04 && DataType.fVal>0.5 -#MMovieWrite.TargetLength: 5 -#MMovieWrite.NumEvents: 500 -#MMovieWrite.Threshold: 3 -#MMovieWrite.FileName: filename.mpg - -# ------------------------------------------------------------------------- -# This is a special option for the datacheck. A task executed as last -# task in the event processing task list - in the standard software it -# is skipped -# ------------------------------------------------------------------------- -#MJCalibrateSignal.FinalFantasy: MHCamFlorian Index: trunk/Mars/callisto_mux.rc =================================================================== --- trunk/Mars/callisto_mux.rc (revision 19871) +++ (revision ) @@ -1,479 +1,0 @@ -# ========================================================================== -############################################################################# -# ========================================================================== -# General -# ========================================================================== -############################################################################# -# ========================================================================== - -# ------------------------------------------------------------------------- -# Use this if you want to setup the logging stream for the jobs -# (overwrites command line options) -# ------------------------------------------------------------------------- -#MLog.VerbosityLevel: 2 -#MLog.DebugLevel: 1 -#MLog.NoColors: yes - -# ========================================================================== -############################################################################# -# ========================================================================== -# Calibration -# ========================================================================== -############################################################################# -# ========================================================================== - -# ------------------------------------------------------------------------- -# Use this if you want to write the MJPedestal output somewhere -# If you don't want it, it is written to the calibration output anyhow. -# ------------------------------------------------------------------------- -#MJPedestalC1.PathOut: . -#MJPedestalC2.PathOut: . - -# ------------------------------------------------------------------------- -# Use this if you want to modify the exclusion of pixels from the beginning -# (default: Pixel 0 and 559) -# ------------------------------------------------------------------------- -MJPedestalC1.BadPixelsFile: mjobs/badpixels_0_559.rc - -# ------------------------------------------------------------------------- -# Use this if you want to modify files containing information about the -# reference lines in the status display -# ------------------------------------------------------------------------- -MJPedestalC1.ReferenceFile: resources/pedestalref_mux.rc -MJPedestalC2.ReferenceFile: resources/pedestalref_mux.rc -MJCalibration.ReferenceFile: resources/calibrationref_mux.rc -MJCalibration.MHCalibrationRelTimeCam.ReferenceFile: resources/calibrationref_mux.rc -MJCalibration.MHCalibrationChargeCam.ReferenceFile: resources/calibrationref_mux.rc - -# ------------------------------------------------------------------------- -# Use this to define where the calibration output is stored. The filename -# is created from the sequence number. If nothing is specified '.' is -# assumed. (overwrites command line options) -# ------------------------------------------------------------------------- -#MJCalibration.PathOut: calped - -# ------------------------------------------------------------------------- -# Use this to define where the program should search for the pedestal -# and calibration files defined in the sequence. To use the local -# directory use '.' If nothing is specified the default path in the -# datacenter is used. (overwrites command line options) -# ------------------------------------------------------------------------- -#MJPedestalC1.PathData: /data/MAGIC/Period016/ -#MJPedestalC2.PathData: /data/MAGIC/Period016/ -#MJCalibration.PathData: /data/MAGIC/Period016/ - -# ------------------------------------------------------------------------- -# Define here which signal extractor you are going to use for the -# whole calibration chain -# ------------------------------------------------------------------------- -MJPedestalC1.ExtractSignal: MExtractTimeAndChargeSpline -# ------------------------------------------------------------------------- -# Possibilities for a different extractor -# ------------------------------------------------------------------------- -#MJPedestalC1.ExtractSignal: MExtractTimeAndChargeDigitalFilter, MExtractTimeAndChargeDigitalFilterPeakSearch, MExtractTimeAndChargeSlidingWindow -# ------------------------------------------------------------------------- -# Define here parameters valid for all extractors above -# The actual defaults used depend on the extractor in use. -# For details see the corresponding source code -# ------------------------------------------------------------------------- -MJPedestalC1.ExtractSignal.HiGainFirst: 0 -MJPedestalC1.ExtractSignal.HiGainLast: 49 -#MJPedestalC1.ExtractSignal.LoGainFirst: 0 -#MJPedestalC1.ExtractSignal.LoGainLast: 0 -MJPedestalC1.ExtractSignal.SaturationLimit: 252 -#MJPedestalC1.ExtractSignal.OffsetLoGain: 1.51 -#MJPedestalC1.ExtractSignal.LoGainSwitch: 150 -#MJPedestalC1.ExtractSignal.LoGainStartShift: -2.8 -# ------------------------------------------------------------------------- -# Define here parameters valid for the spline: -# ------------------------------------------------------------------------- -#MJPedestalC1.ExtractSignal.Resolution: 0.05 -MJPedestalC1.ExtractSignal.RiseTimeHiGain: 4.0 -MJPedestalC1.ExtractSignal.FallTimeHiGain: 5.0 -#MJPedestalC1.ExtractSignal.LoGainStretch: 1.5 -#MJPedestalC1.ExtractSignal.ExtractionType: Integral -#MJPedestalC1.ExtractSignal.ExtractionType: Amplitude - -# ------------------------------------------------------------------------- -# In the case of the Digital filter you have to overwrite the default -# weights file depending on what you are extracting -# ------------------------------------------------------------------------- -#MJPedestalC2.ExtractSignal.WeightsFile: msignal/calibration_weights_UV4.dat -#MJPedestalY2.ExtractSignal.WeightsFile: msignal/cosmics_weights4.dat -#MJPedestalY3.ExtractSignal.WeightsFile: msignal/cosmics_weights4.dat -#MJCalibration.ExtractSignal.WeightsFile: msignal/calibration_weights_UV4.dat -#MJCalibrateSignal.ExtractSignal.WeightsFile: msignal/cosmics_weights4.dat -#MJCalibrateSignal.ExtractInterlaced.WeightsFile: msignal/calibration_weights_UV4.dat -#MJCalibrateSignal.ExtractSignal.LoGainStretch: 1.5 -#MJCalibrateSignal.ExtractInterlaced.LoGainStretch: 1.5 - -# ------------------------------------------------------------------------- -# Configure MJCalibration -# ------------------------------------------------------------------------- - -# Switch on relative time calibration -#MJCalibration.RelTimeCalibration: Yes -# Set color to be used -#MJCalibration.Color: -# Type of used data format: raw, root, mc -#MJCalibration.DataType: Root -# Type of displayed plots: Full,DataCheck,Normal -#MJCalibration.Display: datacheck -# Write additinal debug output -#MJCalibration.Debug: No -# Use blind pixel -#MJCalibration.UseBlindPixel: No -# Use pin diode -#MJCalibration.UsePINDiode: No -# Use pixel check option -#MJCalibration.PixelCheck: No -#MJCalibTest.PixelCheck: No - -# ------------------------------------------------------------------------- -# Individual events check for calibration results -# ------------------------------------------------------------------------- - -#ContCosmics.MaxEmptyPixels: 0.2 -#ContCosmics.MaxExcludedFraction: 1 -#ContCosmics.MinAcceptedFraction: 0 - -# ------------------------------------------------------------------------- -# In case of MOON DATA, uncomment the following lines: -# (better: call callisto with option "-moon") -# ------------------------------------------------------------------------- -#MJPedestalC1.UseHists: yes -#MJPedestalC2.UseHists: yes -#MJPedestalY1.UseHists: yes -#MJPedestalY2.UseHists: yes -#MJPedestalY3.UseHists: yes - -# ------------------------------------------------------------------------- -# You can choose the pedestal extraction algorithm/task. To use -# MPedCalcPedRun use the lines below. Be carefull, a pedestal file -# could be (if no pedestal file available) a data-file, too. -# ------------------------------------------------------------------------- -#MJPedestalC1.ExtractPedestal: MPedCalcPedRun -#MJPedestalC1.ExtractPedestal.ExtractWinFirst: 0 -#MJPedestalC1.ExtractPedestal.ExtractWinSize: 6 -#MJPedestalC1.ExtractPedestal.PedestalUpdate: no -#MJPedestalC1.ExtractPedestal.CheckWinFirst: 0 -#MJPedestalC1.ExtractPedestal.CheckWinSize: 29 - -#MJPedestalC2.ExtractPedestal: MPedCalcPedRun -#MJPedestalC2.ExtractPedestal.ExtractWinFirst: 0 -#MJPedestalC2.ExtractPedestal.ExtractWinSize: 6 -#MJPedestalC2.ExtractPedestal.PedestalUpdate: no -#MJPedestalC2.ExtractPedestal.CheckWinFirst: 0 -#MJPedestalC2.ExtractPedestal.CheckWinSize: 29 - -# ------------------------------------------------------------------------- -# You can choose the pedestal extraction algorithm/task. To use -# MPedCalcFromLoGain use the lines below. Be carefull, a pedestal file -# could be (if no pedestal file available) a data-file, too. -# ------------------------------------------------------------------------- -#MJPedestalC1.ExtractPedestal: MPedCalcFromLoGain -#MJPedestalC1.ExtractPedestal.ExtractWinFirst: 0 -#MJPedestalC1.ExtractPedestal.ExtractWinSize: 6 -#MJPedestalC1.ExtractPedestal.PedestalUpdate: no -#MJPedestalC1.ExtractPedestal.CheckWinFirst: 0 -#MJPedestalC1.ExtractPedestal.CheckWinSize: 29 -#MJPedestalC1.ExtractPedestal.MaxSignalVar: 40 - -#MJPedestalC2.ExtractPedestal: MPedCalcFromLoGain -#MJPedestalC2.ExtractPedestal.ExtractWinFirst: 0 -#MJPedestalC2.ExtractPedestal.ExtractWinSize: 6 -#MJPedestalC2.ExtractPedestal.PedestalUpdate: no -#MJPedestalC2.ExtractPedestal.CheckWinFirst: 0 -#MJPedestalC2.ExtractPedestal.CheckWinSize: 29 -#MJPedestalC2.ExtractPedestal.MaxSignalVar: 40 - -# ------------------------------------------------------------------------- -# Configure MJPedestalC1 and MJPedestalC2 -# ------------------------------------------------------------------------- -# Possible data type: Raw,Root,MC -#MJPedestalC1.DataType: Root -#MJPedestalC2.DataType: Root -#MJPedestalC3.DataType: Root -# Type of displayed plots: Full,DataCheck,None -#MJPedestalC1.Display: datacheck -#MJPedestalC2.Display: datacheck -#MJPedestalC3.Display: datacheck -# Maximum number of event processed in the loop -#MJPedestalC1.MaxEvents: 1000 -#MJPedestalC2.MaxEvents: 1000 -# Allow to overwrite existing files with the output file -#MJPedestalC1.AllowOverwrite: No -#MJPedestalC2.AllowOverwrite: No -# Use data runs from the sequence instead of calibration runs -#MJPedestalC1.UseData: No -#MJPedestalC2.UseData: No - -# ------------------------------------------------------------------------- -# Use this if you want to change the software low-gain switch -# ------------------------------------------------------------------------- -#MJCalibration.ExtractSignal.LoGainSwitch: 150 - -# ------------------------------------------------------------------------- -# Use this if you want to change the higain-vs. logain intercalibration -# ------------------------------------------------------------------------- -#MJCalibration.HiLoCalibFile: resources/hilocalib_sp1.rc - -# ------------------------------------------------------------------------- -# Use this to change the behaviour of the calibration -# ------------------------------------------------------------------------- -# Type if you set a colour explicitely from outside (only for MC!!!) -#MJCalibration.MCalibColorSet.ExplicitColor: green,blue,uv,ct1 - -#MJCalibration.MCalibrationChargeCalc.ArrTimeRmsLimit: 5.0 -#MJCalibration.MCalibrationChargeCalc.ChargeLimit: 2.5 -#MJCalibration.MCalibrationChargeCalc.ChargeErrLimit: 0 -#MJCalibration.MCalibrationChargeCalc.ChargeRelErrLimit: 1 -#MJCalibration.MCalibrationChargeCalc.Debug: no - -#MJCalibration.MCalibrationChargeCalc.FFactorErrLimit: 4.5 -#MJCalibration.MCalibrationChargeCalc.LambdaErrLimit: 0.2 -#MJCalibration.MCalibrationChargeCalc.LambdaCheckLimit: 0.5 -#MJCalibration.MCalibrationChargeCalc.PheErrLowerLimit: 9.0 -#MJCalibration.MCalibrationChargeCalc.PheErrUpperLimit: 5.5 - -#MJCalibration.MHCalibrationChargeCam.Debug: no -#MJCalibration.MHCalibrationChargeCam.LoGain: yes -#MJCalibration.MHCalibrationChargeCam.Oscillations: yes -#MJCalibration.MHCalibrationChargeCam.SizeCheck: yes -#MJCalibration.MHCalibrationChargeCam.Averageing: yes -MJCalibration.MHCalibrationChargeCam.HiGainNbins: 800 -MJCalibration.MHCalibrationChargeCam.HiGainFirst: -0.5 -MJCalibration.MHCalibrationChargeCam.HiGainLast: 3999.5 -#MJCalibration.MHCalibrationChargeCam.LoGainNbins: 500 -#MJCalibration.MHCalibrationChargeCam.LoGainFirst: -100.25 -#MJCalibration.MHCalibrationChargeCam.LoGainLast: 899.75 -#MJCalibration.MHCalibrationChargeCam.TimeLowerLimit: 1. -#MJCalibration.MHCalibrationChargeCam.TimeUpperLimit: 3. -#MJCalibration.MHCalibrationChargeCam.NumHiGainSaturationLimit: 0.02 -#MJCalibration.MHCalibrationChargeCam.NumLoGainSaturationLimit: 0.005 -MJCalibration.MHCalibrationChargeCam.ProbLimit: 1.0E-18 -#MJCalibration.MHCalibrationChargeCam.OverflowLimit: 0.005 -#MJCalibration.MHCalibrationChargeCam.PulserFrequency: 500 - -#MJCalibration.MHCalibrationRelTimeCam.Debug: no -#MJCalibration.MHCalibrationRelTimeCam.LoGain: no -#MJCalibration.MHCalibrationRelTimeCam.Oscillations: yes -#MJCalibration.MHCalibrationRelTimeCam.SizeCheck: yes -#MJCalibration.MHCalibrationRelTimeCam.Averageing: yes -MJCalibration.MHCalibrationRelTimeCam.Nbins: 751 -MJCalibration.MHCalibrationRelTimeCam.First: -60 -MJCalibration.MHCalibrationRelTimeCam.Last: 60 -#MJCalibration.MHCalibrationRelTimeCam.NumHiGainSaturationLimit: 0.25 -#MJCalibration.MHCalibrationRelTimeCam.ProbLimit: 0.00001 -#MJCalibration.MHCalibrationRelTimeCam.OverflowLimit: 0.005 -#MJCalibration.MHCalibrationRelTimeCam.PulserFrequency: 500 - -# ========================================================================== -############################################################################# -# ========================================================================== -# Signal extraction -# ========================================================================== -############################################################################# -# ========================================================================== - -# ------------------------------------------------------------------------- -# Use a OutputPath if you want to write the results to a file -# ------------------------------------------------------------------------- -#MJCalibrateSignal.PathOut: . -#MJPedestalY1.PathOut: . -#MJPedestalY2.PathOut: . -#MJPedestalY3.PathOut: . -#MJPedestalY1.PathData: /data/MAGIC/Period016/ -#MJPedestalY2.PathData: /data/MAGIC/Period016/ -#MJPedestalY3.PathData: /data/MAGIC/Period016/ -#MJCalibrateSignal.PathData: /data/MAGIC/Period016/ -#MJCalibrateSignal.PathIn: . - -# ------------------------------------------------------------------------- -# If you want to change the behaviour of the pedestal extraction... -# ------------------------------------------------------------------------- -# Possible data type: Raw,Root,MC -#MJPedestalY1.DataType: Root -#MJPedestalY2.DataType: Root -#MJPedestalY3.DataType: Root -#MJPedestalY1.UseData: Yes -#MJPedestalY2.UseData: Yes -#MJPedestalY3.UseData: Yes -#MJPedestalY1.MinEvents: 50 -#MJPedestalY2.MinEvents: 50 -#MJPedestalY3.MinEvents: 50 - -MJPedestalY1.MaxEvents: 50000 -MJPedestalY2.MaxEvents: 100000 -MJPedestalY3.MaxEvents: 50000 - -MJPedestalY1.MaxPedestals: 500 -MJPedestalY2.MaxPedestals: 500 -MJPedestalY3.MaxPedestals: 500 - -# ------------------------------------------------------------------------- -# When the starting pedestal is extracted from the first events of -# the data run also the average pulse position of these events -# is determined. If the pulse position check is switched on (default) -# the extraction window of the hi- and lo-gain extractor for the -# determined average pulse position is not allowed to exceed the -# extractors search region. In this scenario we assume that the -# start of the extraction window is at least similar to the -# rising edge of the pulse. -# To make sure that the search region and the extraction window -# does not touch increase the values. -# ------------------------------------------------------------------------- -#MJPedestalY2.PulsePosCheck: yes -#MJPedestalY2.ExtractWinLeft: 0.0 -#MJPedestalY2.ExtractWinRight: 0.0 - -# ------------------------------------------------------------------------- -# Define the Pulse Position check parameters: -# ------------------------------------------------------------------------- -MJPedestalY2.MHCalibrationPulseTimeCam.SaturationLimit: 200 -MJPedestalY2.MHCalibrationPulseTimeCam.LowerSignalLimit: 35 -#MJPedestalY2.MHCalibrationPulseTimeCam.NumPixelsRequired: 3 -#MJPedestalY2.PixelCheck: no - -#MJPedestalY.ExtractPedestal: MPedCalcFromLoGain -#MJPedestalY.ExtractPedestal.PedestalUpdate: no - -# ------------------------------------------------------------------------- -# Setup pedestal extractor for running through data -# ------------------------------------------------------------------------- -#MJCalibrateSignal.MPedCalcFundamental.PedestalUpdate: yes -#MJCalibrateSignal.MPedCalcFundamental.NumDump: 500 -#MJCalibrateSignal.MPedCalcFundamental.NumEventsDump: 500 -#MJCalibrateSignal.MPedCalcFundamental.NumAreasDump: 500 -#MJCalibrateSignal.MPedCalcFundamental.NumSectorsDump: 500 -#MJCalibrateSignal.MPedCalcFundamental.ExtractWinFirst: 17 -#MJCalibrateSignal.MPedCalcFundamental.ExtractWinSize: 6 -#MJCalibrateSignal.MPedCalcFundamental.CheckWinFirst: 0 -#MJCalibrateSignal.MPedCalcFundamental.CheckWinSize: 29 -#MJCalibrateSignal.MPedCalcFundamental.MaxSignalVar: 40 - -#MJCalibrateSignal.MPedCalcFromExtractorRndm.PedestalUpdate: yes -#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumEventsDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumAreasDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumSectorsDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.ExtractWinFirst: 17 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.ExtractWinSize: 6 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.CheckWinFirst: 0 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.CheckWinSize: 29 -#MJCalibrateSignal.MPedCalcFromExtractorRndm.MaxSignalVar: 40 - -#MJCalibrateSignal.MPedCalcFromExtractor.PedestalUpdate: yes -#MJCalibrateSignal.MPedCalcFromExtractor.NumDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractor.NumEventsDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractor.NumAreasDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractor.NumSectorsDump: 500 -#MJCalibrateSignal.MPedCalcFromExtractor.ExtractWinFirst: 17 -#MJCalibrateSignal.MPedCalcFromExtractor.ExtractWinSize: 6 -#MJCalibrateSignal.MPedCalcFromExtractor.CheckWinFirst: 0 -#MJCalibrateSignal.MPedCalcFromExtractor.CheckWinSize: 29 -#MJCalibrateSignal.MPedCalcFromExtractor.MaxSignalVar: 40 - -# ------------------------------------------------------------------------- -# Configure the interlaced calibration updates -# ------------------------------------------------------------------------- -#MJCalibrateSignal.MCalibCalcFromPast.UpdateWithFFactorMethod: yes -#MJCalibrateSignal.MCalibCalcFromPast.NumEventsDump: 500 -#MJCalibrateSignal.MCalibCalcFromPast.UpdateNumPhes: yes -#MJCalibrateSignal.MCalibCalcFromPast.NumPhesDump: 10 -#MJCalibrateSignal.MCalibrationChargeCalc.Debug: no - -# ------------------------------------------------------------------------- -# Use this if you want to change the signal extractor -# ------------------------------------------------------------------------- -#MJCalibrateSignal.ExtractSignal.LoGainSwitch: 150 - -# ------------------------------------------------------------------------- -# Setup calibration of data -# ------------------------------------------------------------------------- -# Type of used data format: raw,root,MC -#MJCalibrateSignal.DataType: Root -# Type if you set a colour explicitely from outside (only for MC!!!) -#MJCalibrateSignal.MCalibColorSet.ExpicitColor: green,blue,uv,ct1 -#MJCalibrateSignal.MCalibrateData.PedestalFlag: Event -#MJCalibrateSignal.MCalibrateData.ScaleFactor: 0.85182358 -MJCalibrateSignal.MCalibrateData.FileNameScale: resources/calibration_spline.rc -#MJCalibrateSignal.MCalibrateData.CalibrationMode: Default -#MJCalibrateSignal.MCalibrateData.CalibConvMinLimit: 0.01 -#MJCalibrateSignal.MCalibrateData.CalibConvMaxLimit: 5. -#MJCalibrateSignal.Interlaced: yes -#MJCalibrateSignal.RelTimesUpdate: yes - -# ------------------------------------------------------------------------- -# Setup level for determining of bad pixels -# ------------------------------------------------------------------------- -#MJCalibrateSignal.MBadPixelsCalc.PedestalLevel: 3.0 -#MJCalibrateSignal.MBadPixelsCalc.PedestalLevelVariance: 3.0 -#MJCalibrateSignal.MBadPixelsTreat.NumMinNeighbors: 3 -#MJCalibrateSignal.MBadPixelsTreat.MaxArrivalTimeDiff: 3.0 -#MJCalibrateSignal.MBadPixelsTreat.UseInterpolation: yes -#MJCalibrateSignal.MBadPixelsTreat.ProcessPedestalEvt: yes -#MJCalibrateSignal.MBadPixelsTreat.ProcessPedestalRun: no -#MJCalibrateSignal.MBadPixelsTreat.ProcessTimes: yes -#MJCalibrateSignal.MBadPixelsTreat.UseCentralPixel: no -#MJCalibrateSignal.MBadPixelsTreat.HardTreatment: no - -# ------------------------------------------------------------------------- -# Setup interlaced calibration events histogramming -# ------------------------------------------------------------------------- -#MJCalibrateSignal.MHCalibrationChargeCam.Debug: no -#MJCalibrateSignal.MHCalibrationChargeCam.LoGain: yes -MJCalibrateSignal.MHCalibrationChargeCam.Oscillations: no -#MJCalibrateSignal.MHCalibrationChargeCam.SizeCheck: yes -#MJCalibrateSignal.MHCalibrationChargeCam.Averageing: yes -#MJCalibrateSignal.MHCalibrationChargeCam.HiGainNbins: 500 -#MJCalibrateSignal.MHCalibrationChargeCam.HiGainFirst: -100.5 -#MJCalibrateSignal.MHCalibrationChargeCam.HiGainLast: 1899.5 -MJCalibrateSignal.MHCalibrationChargeCam.LoGainNbins: 800 -MJCalibrateSignal.MHCalibrationChargeCam.LoGainFirst: -0.5 -MJCalibrateSignal.MHCalibrationChargeCam.LoGainLast: 3999.5 -#MJCalibrateSignal.MHCalibrationChargeCam.TimeLowerLimit: 1. -#MJCalibrateSignal.MHCalibrationChargeCam.TimeUpperLimit: 3. -#MJCalibrateSignal.MHCalibrationChargeCam.NumHiGainSaturationLimit: 0.02 -#MJCalibrateSignal.MHCalibrationChargeCam.NumLoGainSaturationLimit: 0.005 -#MJCalibrateSignal.MHCalibrationChargeCam.ProbLimit: 0.0001 -#MJCalibrateSignal.MHCalibrationChargeCam.OverflowLimit: 0.005 -MJCalibrateSignal.MHCalibrationChargeCam.PulserFrequency: 25 - -# ------------------------------------------------------------------------- -# Setup relative time interlaced calibration -# ------------------------------------------------------------------------- -#MJCalibrateSignal.MHCalibrationRelTimeCam.Debug: no -#MJCalibrateSignal.MHCalibrationRelTimeCam.LoGain: no -MJCalibrateSignal.MHCalibrationRelTimeCam.Oscillations: no -#MJCalibrateSignal.MHCalibrationRelTimeCam.SizeCheck: yes -#MJCalibrateSignal.MHCalibrationRelTimeCam.Averageing: yes -MJCalibrateSignal.MHCalibrationRelTimeCam.Nbins: 751 -MJCalibrateSignal.MHCalibrationRelTimeCam.First: -60 -MJCalibrateSignal.MHCalibrationRelTimeCam.Last: 60 -#MJCalibrateSignal.MHCalibrationRelTimeCam.NumHiGainSaturationLimit: 0.25 -#MJCalibrateSignal.MHCalibrationRelTimeCam.ProbLimit: 0.00001 -#MJCalibrateSignal.MHCalibrationRelTimeCam.OverflowLimit: 0.005 -MJCalibrateSignal.MHCalibrationRelTimeCam.PulserFrequency: 25 - -# ------------------------------------------------------------------------- -# Resource for movie encoding. For more details see IACT wiki -# ------------------------------------------------------------------------- -#MMoviePrepare.FirstSlice: 10 -#MMoviePrepare.LastSlice: 50 -#MovieFilter.Condition: MMovieData.fMax>5*MMovieData.fMedianPedestalRms -#MImgCleanStd.CleanLevel1: 8.0 -#MImgCleanStd.CleanLevel2: 4.5 -#MFEvtNumber.FileName: ganymed00001218.root -#MFEvtNumber.Selector: ThetaSquared.fVal<0.04 && DataType.fVal>0.5 -#MMovieWrite.TargetLength: 5 -#MMovieWrite.NumEvents: 500 -#MMovieWrite.Threshold: 3 -#MMovieWrite.FileName: filename.mpg - -# ------------------------------------------------------------------------- -# This is a special option for the datacheck. A task executed as last -# task in the event processing task list - in the standard software it -# is skipped -# ------------------------------------------------------------------------- -#MJCalibrateSignal.FinalFantasy: MHCamFlorian Index: trunk/Mars/ganymed.rc =================================================================== --- trunk/Mars/ganymed.rc (revision 19871) +++ (revision ) @@ -1,221 +1,0 @@ -# ========================================================================== -############################################################################# -# ========================================================================== -# General -# ========================================================================== -############################################################################# -# ========================================================================== - -# ------------------------------------------------------------------------- -# Use this if you want to setup the logging stream for the jobs -# (overwrites command line options) -# ------------------------------------------------------------------------- -#MLog.VerbosityLevel: 2 -#MLog.DebugLevel: 1 -#MLog.NoColors: yes - -# ========================================================================== -############################################################################# -# ========================================================================== -# Cuts -# ========================================================================== -############################################################################# -# ========================================================================== - -# ------------------------------------------------------------------------- -# Use this if you want to write output to somewhere here you can give -# the output path -# ------------------------------------------------------------------------- -#MJCuts.PathOut: . - -# ------------------------------------------------------------------------- -# Use this to setup binnings. For more details see: MBinning::ReadEnv -# ------------------------------------------------------------------------- -BinningAlpha.Raw: 36 0 90 -BinningFalseSource.Raw: 91 -1.5 1.5 -BinningWidth.Raw: 50 0 1 -BinningLength.Raw: 25 0 1 -BinningDist.Raw: 25 0 2 -BinningAsym.Raw: 51 -0.7 0.7 -BinningM3Long.Raw: 51 -0.7 0.7 -BinningM3Trans.Raw: 51 -0.7 0.7 -BinningSlope.Raw: 51 -7.5 7.5 - -# ------------------------------------------------------------------------- -# Using the starguider for pointing correction. -# To switch it off use "MPointingDevCalc.MaxAbsDev: -1" -# For details see the class reference of MPointingDevCalc -# ------------------------------------------------------------------------- -#MPointingDevCalc.MaxAbsDev: 15 -#MPointingDevCalc.NumMinStars: 8 -#MPointingDevCalc.NsbLevel: 3.0 -#MPointingDevCalc.NsbMin: 30 -#MPointingDevCalc.NsbMax: 60 -#MPointingDevCalc.Dx -0.001 -#MPointingDevCalc.Dy -0.004 -#MPointingDevCalc.MaxAge 1.0 - -MPointingDevCalc.PointingModels: 1 81855 83722 89180 267253 291104 328198 1000534 - -# ------------------------------------------------------------------------- -# Setup misfocussing correction (see class reference for MSrcPosCorrect) -# ------------------------------------------------------------------------- -#MSrcPosCorrect.Dx: 0 -#MSrcPosCorrect.Dy: 0 - -# ------------------------------------------------------------------------- -# Setup fit mode and ranges -# ------------------------------------------------------------------------- -#MAlphaFitter.SignalIntegralMax: 15 -#MAlphaFitter.SignalFitMax: 25 -#MAlphaFitter.BackgroundFitMin: 35 -#MAlphaFitter.BackgroundFitMax: 75 -#MAlphaFitter.ScaleMin: 35 -#MAlphaFitter.ScaleMax: 80 -#MAlphaFitter.PolynomOrder: 2 -#MAlphaFitter.ScaleMode: None - -# ------------------------------------------------------------------------- -# Define here which histogram class to use to determin the signal. -# Currently availyble: MHAlpha, MHTheta -# ------------------------------------------------------------------------- -#MJCut.NameHist: MHAlpha - -# Be carefull: The spectrum might be done with a different binning -# This might change the scale factor -#MHThetaSq.NumBinsSignal: 3 -#MHThetaSq.NumBinsTotal: 75 - -# ------------------------------------------------------------------------- -# To get a meaningfull false source plot in wobble mode use magic cuts -# and the following two lines. -# ------------------------------------------------------------------------- -MJCut.NameHistFS: MHDisp - -MHDisp.Smearing: 0.06 -#MHDisp.Wobble: Yes -#MHDisp.ScaleMin: 0.325 -#MHDisp.ScaleMax: 0.475 - -# ------------------------------------------------------------------------- -# To get/use a random source position in the off data with the distribution -# in the camera of the on-data -# ------------------------------------------------------------------------- -#MJCut.RandomSourcePosition: Yes - -# ------------------------------------------------------------------------- -# There are three cuts: -# CutQ: Cut after reading file -# Cut0: After calculation of src dependent parameters, energy estimation -# and disp calculation, before writing summary file -# Cut1: After writing summary file before filling false source plots -# Cut2: After filling false source plots before filling alpha plots -# Cut3: After filling alpha plots (eg. Alpha cuts) for image parameter -# display, eg. write only events after alpha cut to ganymed file -# ------------------------------------------------------------------------- -CutQ.Inverted: Yes -Cut0.Inverted: Yes -Cut1.Inverted: Yes -Cut2.Inverted: Yes -Cut3.Inverted: Yes - -# Quality-cuts (remove unphysical and obvious non-gamma events) -CutQ.Condition: {0} && {1} && {2} && {3} && {4} && {5} -CutQ.0: MImagePar.fNumIslands<3 -CutQ.1: MNewImagePar.fNumUsedPixels>5 -CutQ.2: MNewImagePar.fLeakage1 < 0.3 -CutQ.3: log10(MNewImagePar.fConcCOG)<-0.45+0.08*(log10(MHillas.fSize)<3.9)*(log10(MHillas.fSize)-3.9)^2 -CutQ.4: log10(MNewImagePar.fConc1) <-0.75+0.10*(log10(MHillas.fSize)<3.8)*(log10(MHillas.fSize)-3.8)^2 -CutQ.5: MHillas.fLength>-3.6*(log10(MHillas.fSize)-6.0)^2+70 - -# ------------------------------------------------------------------------- -# Such a cut could be used to exclude strange pointing positions. like -# the camera center in wobble mode. Note, that a cut in the poiting -# position directly influences the effective observation time and should -# only be applied if the events are lost anyhow. In both cases you must -# correct for the loss of observation time calculating your flux! -# ------------------------------------------------------------------------- -#Cut0.Condition: MSrcPosCam.GetDist*MGeomCam.fConvMm2Deg>0.3 - -# ---------- SETUP FOR ONOFF-MODE ----------- -MAlphaFitter.ScaleMin: 0.137 -MAlphaFitter.ScaleMax: 0.640 -MAlphaFitter.BackgroundFitMin: 0.137 -MAlphaFitter.BackgroundFitMax: 1.000 -MAlphaFitter.PolynomOrder: 1 -MAlphaFitter.ScaleMode: OffRegion -MAlphaFitter.SignalFunction: ThetaSq - -# For more details on the magic cuts see its class reference -Cut1.Condition: MFMagicCuts -Cut1.ThetaCut: None -# Parametrization of Disp -Cut1.Param0: 1.15136 -Cut1.Param8: 0.0681437 -Cut1.Param9: 2.62932 -Cut1.Param10: 1.51279 -Cut1.Param11: 0.0507821 -# Parametrization for sign of disp (m3long, slope) -Cut1.Param5: -0.07 -Cut1.Param6: 7.2 -Cut1.Param7: 0.5 -# ThetaSq-Cut -Cut1.Param1: 0.215 -# Area-Cut -Cut1.Param2: 0.215468 -Cut1.Param3: 5.63973 -Cut1.Param4: 0.0836169 - -# ---------- SETUP FOR WOBBLE-MODE ----------- -#MAlphaFitter.ScaleMode: none -#MAlphaFitter.BackgroundFitMin: 0.09 -#MAlphaFitter.BackgroundFitMax: 0.25 -#MAlphaFitter.PolynomOrder: 1 -#MAlphaFitter.SignalFunction: ThetaSq - -# ---------- SETUP TO USE RANDOM FOREST ----------- - -# To switch on random forest uncomment the following lines -#Cut1.HadronnessCut: Hadronness -#Cut1.Param13: 0.5 -#Cut1.Param14: 100 - -# To use both, hadronness and area cut do -#Cut1.HadronnessCut: All - -# ------------------------------------------------------------------------- -# This cut is applied exclusively to the plots in the tab CutT. It is -# resonable to show the events within your signal region here to -# see if your signal region is clean from strange events (e.g. sparks) -# ------------------------------------------------------------------------- -CutT.Condition: ThetaSquared.fVal < ThetaSquaredCut.fVal - -# ------------------------------------------------------------------------- -# Use this to executa a task (eg to calc hadronness) before Cut1 -# ------------------------------------------------------------------------- -# CalcHadronness: MRanForestCalc -# CalcHadronness.FileName: rf-separation.root -# CalcHadronness.NameOutput: Hadronness -# CalcHadronness.Debug: No - -# ------------------------------------------------------------------------- -# Use this to executa a task (eg to calc disp) before Cut1 -# ------------------------------------------------------------------------- -# CalcDisp: MRanForestCalc -# CalcDisp.FileName: rf-disp.root -# CalcDisp.NameOutput: Disp -# CalcDisp.Debug: No -# Cut1.CalcDisp: No - -# CalcDisp: MParameterCalc -# CalcDisp.Rule: 1.15*(1-MHillas.fWidth/MHillas.fLength) -# CalcDisp.NameOutput: Disp -# Cut1.CalcDisp: No - -# ------------------------------------------------------------------------- -# Use this to executa a task (eg to estimate energy) after all Cuts -# ------------------------------------------------------------------------- -# EstimateEnergy: MRanForestCalc -# EstimateEnergy.NameOutput: MEnergyEst -# EstimateEnergy.FileName: ranforest.root -# EstimateEnergy.Debug: No Index: trunk/Mars/ganymed_onoff.rc =================================================================== --- trunk/Mars/ganymed_onoff.rc (revision 19871) +++ (revision ) @@ -1,211 +1,0 @@ -# ========================================================================== -############################################################################# -# ========================================================================== -# General -# ========================================================================== -############################################################################# -# ========================================================================== - -# ------------------------------------------------------------------------- -# Use this if you want to setup the logging stream for the jobs -# (overwrites command line options) -# ------------------------------------------------------------------------- -#MLog.VerbosityLevel: 2 -#MLog.DebugLevel: 1 -#MLog.NoColors: yes - -# ========================================================================== -############################################################################# -# ========================================================================== -# Cuts -# ========================================================================== -############################################################################# -# ========================================================================== - -# ------------------------------------------------------------------------- -# Use this if you want to write output to somewhere here you can give -# the output path -# ------------------------------------------------------------------------- -#MJCuts.PathOut: . - -# ------------------------------------------------------------------------- -# Use this to setup binnings. For more details see: MBinning::ReadEnv -# ------------------------------------------------------------------------- -BinningAlpha.Raw: 36 0 90 -BinningFalseSource.Raw: 91 -1.5 1.5 -BinningWidth.Raw: 50 0 1 -BinningLength.Raw: 25 0 1 -BinningDist.Raw: 25 0 2 -BinningAsym.Raw: 51 -0.7 0.7 -BinningM3Long.Raw: 51 -0.7 0.7 -BinningM3Trans.Raw: 51 -0.7 0.7 -BinningSlope.Raw: 51 -7.5 7.5 - -# ------------------------------------------------------------------------- -# Using the starguider for pointing correction. -# To switch it off use "MPointingDevCalc.MaxAbsDev: -1" -# For details see the class reference of MPointingDevCalc -# ------------------------------------------------------------------------- -#MPointingDevCalc.MaxAbsDev: 15 -#MPointingDevCalc.NumMinStars: 8 -#MPointingDevCalc.NsbLevel: 3.0 -#MPointingDevCalc.NsbMin: 30 -#MPointingDevCalc.NsbMax: 60 -#MPointingDevCalc.Dx -0.001 -#MPointingDevCalc.Dy -0.004 -#MPointingDevCalc.MaxAge 1.0 - -MPointingDevCalc.PointingModels: 1 81855 83722 89180 267253 291104 328198 1000534 - -# ------------------------------------------------------------------------- -# Setup misfocussing correction (see class reference for MSrcPosCorrect) -# ------------------------------------------------------------------------- -#MSrcPosCorrect.Dx: 0 -#MSrcPosCorrect.Dy: 0 - -# ------------------------------------------------------------------------- -# Setup fit mode and ranges -# ------------------------------------------------------------------------- -#MAlphaFitter.SignalIntegralMax: 15 -#MAlphaFitter.SignalFitMax: 25 -#MAlphaFitter.BackgroundFitMin: 35 -#MAlphaFitter.BackgroundFitMax: 75 -#MAlphaFitter.ScaleMin: 35 -#MAlphaFitter.ScaleMax: 80 -#MAlphaFitter.PolynomOrder: 2 -#MAlphaFitter.ScaleMode: None - -# ------------------------------------------------------------------------- -# Define here which histogram class to use to determin the signal. -# Currently availyble: MHAlpha, MHTheta -# ------------------------------------------------------------------------- -#MJCut.NameHist: MHAlpha - -# Be carefull: The spectrum might be done with a different binning -# This might change the scale factor -#MHThetaSq.NumBinsSignal: 3 -#MHThetaSq.NumBinsTotal: 75 - -# ------------------------------------------------------------------------- -# To get a meaningfull false source plot in wobble mode use magic cuts -# and the following two lines. -# ------------------------------------------------------------------------- -MJCut.NameHistFS: MHDisp - -MHDisp.Smearing: 0.06 -#MHDisp.ScaleMin: 0.325 -#MHDisp.ScaleMax: 0.475 - -# ------------------------------------------------------------------------- -# To get use a random source position in the off data with the distribution -# in the camera of the on-data -# ------------------------------------------------------------------------- -MJCut.RandomSourcePosition: Yes - -# ------------------------------------------------------------------------- -# There are three cuts: -# CutQ: Cut after reading file -# Cut0: After calculation of src dependent parameters, energy estimation -# and disp calculation, before writing summary file -# Cut1: After writing summary file before filling false source plots -# Cut2: After filling false source plots before filling alpha plots -# Cut3: After filling alpha plots (eg. Alpha cuts) for image parameter -# display, eg. write only events after alpha cut to ganymed file -# ------------------------------------------------------------------------- -CutQ.Inverted: Yes -Cut0.Inverted: Yes -Cut1.Inverted: Yes -Cut2.Inverted: Yes -Cut3.Inverted: Yes - -# Quality-cuts (remove unphysical and obvious non-gamma events) -CutQ.Condition: {0} && {1} && {2} && {3} && {4} && {5} -CutQ.0: MImagePar.fNumIslands<3 -CutQ.1: MNewImagePar.fNumUsedPixels>5 -CutQ.2: MNewImagePar.fLeakage1 < 0.3 -CutQ.3: log10(MNewImagePar.fConcCOG)<-0.45+0.08*(log10(MHillas.fSize)<3.9)*(log10(MHillas.fSize)-3.9)^2 -CutQ.4: log10(MNewImagePar.fConc1) <-0.75+0.10*(log10(MHillas.fSize)<3.8)*(log10(MHillas.fSize)-3.8)^2 -CutQ.5: MHillas.fLength>-3.6*(log10(MHillas.fSize)-6.0)^2+70 - -# ------------------------------------------------------------------------- -# Such a cut could be used to exclude strange pointing positions. like -# the camera center in wobble mode. Note, that a cut in the poiting -# position directly influences the effective observation time and should -# only be applied if the events are lost anyhow. In both cases you must -# correct for the loss of observation time calculating your flux! -# ------------------------------------------------------------------------- -#Cut0.Condition: MSrcPosCam.GetDist*MGeomCam.fConvMm2Deg>0.3 - -# ---------- SETUP FOR ONOFF-MODE ----------- -MAlphaFitter.ScaleMin: 0.137 -MAlphaFitter.ScaleMax: 0.640 -MAlphaFitter.BackgroundFitMin: 0.137 -MAlphaFitter.BackgroundFitMax: 1.000 -MAlphaFitter.PolynomOrder: 1 -MAlphaFitter.ScaleMode: OffRegion -MAlphaFitter.SignalFunction: ThetaSq - -# For more details on the magic cuts see its class reference -Cut1.Condition: MFMagicCuts -Cut1.ThetaCut: None -# Parametrization of Disp -Cut1.Param0: 1.15136 -Cut1.Param8: 0.0681437 -Cut1.Param9: 2.62932 -Cut1.Param10: 1.51279 -Cut1.Param11: 0.0507821 -# Parametrization for sign of disp (m3long, slope) -Cut1.Param5: -0.07 -Cut1.Param6: 7.2 -Cut1.Param7: 0.5 -# ThetaSq-Cut -Cut1.Param1: 0.215 -# Area-Cut -Cut1.Param2: 0.215468 -Cut1.Param3: 5.63973 -Cut1.Param4: 0.0836169 - -# To switch on random forest uncomment the following lines -#Cut1.HadronnessCut: Hadronness -#Cut1.Param13: 0.5 -#Cut1.Param14: 100 - -# To use both, hadronness and area cut do -#Cut1.HadronnessCut: All - -# ------------------------------------------------------------------------- -# This cut is applied exclusively to the plots in the tab CutT. It is -# resonable to show the events within your signal region here to -# see if your signal region is clean from strange events (e.g. sparks) -# ------------------------------------------------------------------------- -CutT.Condition: ThetaSquared.fVal < ThetaSquaredCut.fVal - -# ------------------------------------------------------------------------- -# Use this to executa a task (eg to calc hadronness) before Cut1 -# ------------------------------------------------------------------------- -# CalcHadronness: MRanForestCalc -# CalcHadronness.FileName: rf-separation.root -# CalcHadronness.NameOutput: Hadronness -# CalcHadronness.Debug: No - -# ------------------------------------------------------------------------- -# Use this to executa a task (eg to calc disp) before Cut1 -# ------------------------------------------------------------------------- -# CalcDisp: MRanForestCalc -# CalcDisp.FileName: rf-disp.root -# CalcDisp.NameOutput: Disp -# CalcDisp.Debug: No -# Cut1.CalcDisp: No - -# CalcDisp: MParameterCalc -# CalcDisp.Rule: 1.15*(1-MHillas.fWidth/MHillas.fLength) -# CalcDisp.NameOutput: Disp -# Cut1.CalcDisp: No - -# ------------------------------------------------------------------------- -# Use this to executa a task (eg to estimate energy) after all Cuts -# ------------------------------------------------------------------------- -# EstimateEnergy: MRanForestCalc -# EstimateEnergy.NameOutput: MEnergyEst -# EstimateEnergy.FileName: ranforest.root -# EstimateEnergy.Debug: No Index: trunk/Mars/ganymed_wobble.rc =================================================================== --- trunk/Mars/ganymed_wobble.rc (revision 19871) +++ (revision ) @@ -1,211 +1,0 @@ -# ========================================================================== -############################################################################# -# ========================================================================== -# General -# ========================================================================== -############################################################################# -# ========================================================================== - -# ------------------------------------------------------------------------- -# Use this if you want to setup the logging stream for the jobs -# (overwrites command line options) -# ------------------------------------------------------------------------- -#MLog.VerbosityLevel: 2 -#MLog.DebugLevel: 1 -#MLog.NoColors: yes - -# ========================================================================== -############################################################################# -# ========================================================================== -# Cuts -# ========================================================================== -############################################################################# -# ========================================================================== - -# ------------------------------------------------------------------------- -# Use this if you want to write output to somewhere here you can give -# the output path -# ------------------------------------------------------------------------- -#MJCuts.PathOut: . - -# ------------------------------------------------------------------------- -# Use this to setup binnings. For more details see: MBinning::ReadEnv -# ------------------------------------------------------------------------- -BinningAlpha.Raw: 36 0 90 -BinningFalseSource.Raw: 91 -1.5 1.5 -BinningWidth.Raw: 50 0 1 -BinningLength.Raw: 25 0 1 -BinningDist.Raw: 25 0 2 -BinningAsym.Raw: 51 -0.7 0.7 -BinningM3Long.Raw: 51 -0.7 0.7 -BinningM3Trans.Raw: 51 -0.7 0.7 -BinningSlope.Raw: 51 -7.5 7.5 - -# ------------------------------------------------------------------------- -# Using the starguider for pointing correction. -# To switch it off use "MPointingDevCalc.MaxAbsDev: -1" -# For details see the class reference of MPointingDevCalc -# ------------------------------------------------------------------------- -#MPointingDevCalc.MaxAbsDev: 15 -#MPointingDevCalc.NumMinStars: 8 -#MPointingDevCalc.NsbLevel: 3.0 -#MPointingDevCalc.NsbMin: 30 -#MPointingDevCalc.NsbMax: 60 -#MPointingDevCalc.Dx -0.001 -#MPointingDevCalc.Dy -0.004 -#MPointingDevCalc.MaxAge 1.0 - -MPointingDevCalc.PointingModels: 1 81855 83722 89180 267253 291104 328198 1000534 - -# ------------------------------------------------------------------------- -# Setup misfocussing correction (see class reference for MSrcPosCorrect) -# ------------------------------------------------------------------------- -#MSrcPosCorrect.Dx: 0 -#MSrcPosCorrect.Dy: 0 - -# ------------------------------------------------------------------------- -# Setup fit mode and ranges -# ------------------------------------------------------------------------- -#MAlphaFitter.SignalIntegralMax: 15 -#MAlphaFitter.SignalFitMax: 25 -#MAlphaFitter.BackgroundFitMin: 35 -#MAlphaFitter.BackgroundFitMax: 75 -#MAlphaFitter.ScaleMin: 35 -#MAlphaFitter.ScaleMax: 80 -#MAlphaFitter.PolynomOrder: 2 -#MAlphaFitter.ScaleMode: None - -# ------------------------------------------------------------------------- -# Define here which histogram class to use to determin the signal. -# Currently availyble: MHAlpha, MHTheta -# ------------------------------------------------------------------------- -#MJCut.NameHist: MHAlpha - -# ------------------------------------------------------------------------- -# Define here the number of off-source regions used in wobble mode -# ------------------------------------------------------------------------- -#MJCut.NumOffSourcePositions: 3 - -# ------------------------------------------------------------------------- -# To get a meaningfull false source plot in wobble mode use magic cuts -# and the following two lines. -# ------------------------------------------------------------------------- -MJCut.NameHistFS: MHDisp - -MHDisp.Wobble: Yes -MHDisp.Smearing: 0.06 -#MHDisp.ScaleMin: 0.325 -#MHDisp.ScaleMax: 0.475 - -# ------------------------------------------------------------------------- -# Use this to set the number of bins corresponding to the signal region -# and if you use a Anti-Theta cut the second one has to be enabled -# ------------------------------------------------------------------------- -#MHPhi.NumBinsSignal: 2 -#MHPhi.UseAntiPhiCut: No - -# ------------------------------------------------------------------------- -# There are three cuts: -# CutQ: Cut after reading file -# Cut0: After calculation of src dependent parameters, energy estimation -# and disp calculation, before writing summary file -# Cut1: After writing summary file before filling false source plots -# Cut2: After filling false source plots before filling alpha plots -# Cut3: After filling alpha plots (eg. Alpha cuts) for image parameter -# display, eg. write only events after alpha cut to ganymed file -# ------------------------------------------------------------------------- -CutQ.Inverted: Yes -Cut0.Inverted: Yes -Cut1.Inverted: Yes -Cut2.Inverted: Yes -Cut3.Inverted: Yes - -# Quality-cuts (remove unphysical and obvious non-gamma events) -CutQ.Condition: {0} && {1} && {2} && {3} && {4} && {5} -CutQ.0: MImagePar.fNumIslands<3 -CutQ.1: MNewImagePar.fNumUsedPixels>5 -CutQ.2: MNewImagePar.fLeakage1 < 0.3 -CutQ.3: log10(MNewImagePar.fConcCOG)<-0.45+0.08*(log10(MHillas.fSize)<3.9)*(log10(MHillas.fSize)-3.9)^2 -CutQ.4: log10(MNewImagePar.fConc1) <-0.75+0.10*(log10(MHillas.fSize)<3.8)*(log10(MHillas.fSize)-3.8)^2 -CutQ.5: MHillas.fLength>-3.6*(log10(MHillas.fSize)-6.0)^2+70 - -# ------------------------------------------------------------------------- -# Such a cut could be used to exclude strange pointing positions. like -# the camera center in wobble mode. Note, that a cut in the poiting -# position directly influences the effective observation time and should -# only be applied if the events are lost anyhow. In both cases you must -# correct for the loss of observation time calculating your flux! -# ------------------------------------------------------------------------- -#Cut0.Condition: MSrcPosCam.GetDist*MGeomCam.fConvMm2Deg>0.3 - -# ---------- SETUP FOR WOBBLE-MODE ----------- -#MAlphaFitter.ScaleMode: none -MAlphaFitter.BackgroundFitMin: 0.09 -MAlphaFitter.BackgroundFitMax: 0.25 -MAlphaFitter.PolynomOrder: 1 -MAlphaFitter.SignalFunction: ThetaSq - -# For more details on the magic cuts see its class reference -Cut1.Condition: MFMagicCuts -Cut1.ThetaCut: None -# Parametrization of Disp -Cut1.Param0: 1.15136 -Cut1.Param8: 0.0681437 -Cut1.Param9: 2.62932 -Cut1.Param10: 1.51279 -Cut1.Param11: 0.0507821 -# Parametrization for sign of disp (m3long, slope) -Cut1.Param5: -0.07 -Cut1.Param6: 7.2 -Cut1.Param7: 0.5 -# ThetaSq-Cut -Cut1.Param1: 0.215 -# Area-Cut -Cut1.Param2: 0.215468 -Cut1.Param3: 5.63973 -Cut1.Param4: 0.0836169 - -# To switch on random forest uncomment the following lines -#Cut1.HadronnessCut: Hadronness -#Cut1.Param13: 0.5 -#Cut1.Param14: 100 - -# To use both, hadronness and area cut do -#Cut1.HadronnessCut: All - -# ------------------------------------------------------------------------- -# This cut is applied exclusively to the plots in the tab CutT. It is -# resonable to show the events within your signal region here to -# see if your signal region is clean from strange events (e.g. sparks) -# ------------------------------------------------------------------------- -CutT.Condition: ThetaSquared.fVal < ThetaSquaredCut.fVal - -# ------------------------------------------------------------------------- -# Use this to executa a task (eg to calc hadronness) before Cut1 -# ------------------------------------------------------------------------- -# CalcHadronness: MRanForestCalc -# CalcHadronness.FileName: rf-separation.root -# CalcHadronness.NameOutput: Hadronness -# CalcHadronness.Debug: No - -# ------------------------------------------------------------------------- -# Use this to executa a task (eg to calc disp) before Cut1 -# ------------------------------------------------------------------------- -# CalcDisp: MRanForestCalc -# CalcDisp.FileName: rf-disp.root -# CalcDisp.NameOutput: Disp -# CalcDisp.Debug: No -# Cut1.CalcDisp: No - -# CalcDisp: MParameterCalc -# CalcDisp.Rule: 1.15*(1-MHillas.fWidth/MHillas.fLength) -# CalcDisp.NameOutput: Disp -# Cut1.CalcDisp: No - -# ------------------------------------------------------------------------- -# Use this to executa a task (eg to estimate energy) after all Cuts -# ------------------------------------------------------------------------- -# EstimateEnergy: MRanForestCalc -# EstimateEnergy.NameOutput: MEnergyEst -# EstimateEnergy.FileName: ranforest.root -# EstimateEnergy.Debug: No Index: trunk/Mars/mars.rc =================================================================== --- trunk/Mars/mars.rc (revision 19871) +++ (revision ) @@ -1,106 +1,0 @@ -# ========================================================================== -############################################################################# -# ========================================================================== -# -# This is a resource file for the Mars executable. You can here setup -# setup the tasks of the eventloops executed in the Mars executable - -# if supported by these tasks. -# -# ========================================================================== -############################################################################# -# ========================================================================== - -# ------------------------------------------------------------------------- -# Use this if you want to setup the logging stream for the jobs -# (overwrites command line options) -# ------------------------------------------------------------------------- -#MLog.VerbosityLevel: 2 -#MLog.DebugLevel: 1 -#MLog.NoColors: yes - -# ------------------------------------------------------------------------- -# Define here which signal extractor you are going to use for the -# whole calibration chain -# ------------------------------------------------------------------------- -#ExtractSignal: MExtractTimeAndChargeDigitalFilter -#ExtractSignal.WeightsFile: msignal/MC_weights46.dat -# ------------------------------------------------------------------------- -# Define here parameters valid for all extractors above -# ------------------------------------------------------------------------- -#ExtractSignal.LoGainSwitch: 120 -#ExtractSignal.LoGainStartShift: -2.5 -# ------------------------------------------------------------------------- -# Define here parameters valid for the spline: -# ------------------------------------------------------------------------- -#ExtractSignal.HiGainFirst: 0 -#ExtractSignal.HiGainLast: 14 -#ExtractSignal.LoGainFirst: 1 -#ExtractSignal.LoGainLast: 14 -#ExtractSignal.Resolution: 0.05 -#ExtractSignal.RiseTimeHiGain: 0.64 -#ExtractSignal.FallTimeHiGain: 0.76 -#ExtractSignal.LoGainStretch: 1.5 -#ExtractSignal.OffsetLoGain: 1.3 -#ExtractSignal.ExtractionType: Integral -#ExtractSignal.ExtractionType: Amplitude -# ------------------------------------------------------------------------- -# Define here parameters valid for the Digital Filter: -# ------------------------------------------------------------------------- -#ExtractSignal.HiGainFirst: 0 -#ExtractSignal.HiGainLast: 16 -#ExtractSignal.LoGainFirst: 1 -#ExtractSignal.LoGainLast: 14 -#ExtractSignal.OffsetLoGain: 0.95 - -# ========================================================================= -# To look into MUX data use the following setup -# ========================================================================= -#ExtractSignal: MExtractTimeAndChargeSpline -#ExtractSignal.RiseTimeHiGain: 4.3 -#ExtractSignal.FallTimeHiGain: 5.1 -#ExtractSignal.HiGainFirst: 0 -#ExtractSignal.HiGainLast: 49 - -# ------------------------------------------------------------------------- -# In the case of the Digital filter you have to overwrite the default -# weights file depending on what you are extracting -# ------------------------------------------------------------------------- -#ExtractSignal.WeightsFile: msignal/cosmics_weights.dat - -# ------------------------------------------------------------------------- -# Setup level for determining of bad pixels -# ------------------------------------------------------------------------- -#MBadPixelsCalc.PedestalLevel: 3.0 -#MBadPixelsCalc.PedestalLevelVariance: 3.0 -#MBadPixelsTreat.NumMinNeighbors: 3 -#MBadPixelsTreat.UseInterpolation: yes -#MBadPixelsTreat.ProcessPedestalEvt: yes -#MBadPixelsTreat.ProcessPedestalRun: no -#MBadPixelsTreat.ProcessTimes: yes -#MBadPixelsTreat.UseCentralPixel: no -#MBadPixelsTreat.HardTreatment: no - -# ------------------------------------------------------------------------- -# This is how you can setup the Image cleaning used in the Camera Display -# e.g. Standard, Absolute, Scaled, Time -# ------------------------------------------------------------------------- -MImgCleanStd.CleanMethod: absolute -#MImgCleanStd.CleanLevel0: 6.0 -MImgCleanStd.CleanLevel1: 6.0 -MImgCleanStd.CleanLevel2: 3.0 -MImgCleanStd.TimeLevel1: 17.5 -MImgCleanStd.TimeLevel2: 17.5 -MImgCleanStd.PostCleanType: 3 - -#MImgCleanStd.CleanRings: 1 -#MImgCleanStd.KeepSinglePixels: No - -# ------------------------------------------------------------------------- -# Setup a cut to choose the events which should be displayed -# (Use this cuts to display only good Muon candidates) -# ------------------------------------------------------------------------- -#Cut.Inverted: Yes -#Cut.Condition: {0} && {1} && {2} -#Cut.0: MHillas.fSize>150 -#Cut.1: (MMuonSearchPar.fRadius>180) && (MMuonSearchPar.fRadius<400) && (MMuonSearchPar.fDeviation<45) -#Cut.2: (MMuonCalibPar.fArcPhi>190) && (MMuonSearchPar.fDeviation<35) && (MMuonCalibPar.fArcWidth<0.20) && (MMuonCalibPar.fArcWidth>0.04) Index: trunk/Mars/resources/callisto.rc =================================================================== --- trunk/Mars/resources/callisto.rc (revision 19872) +++ trunk/Mars/resources/callisto.rc (revision 19872) @@ -0,0 +1,458 @@ +# ========================================================================== +############################################################################# +# ========================================================================== +# General +# ========================================================================== +############################################################################# +# ========================================================================== + +# ------------------------------------------------------------------------- +# Use this if you want to setup the logging stream for the jobs +# (overwrites command line options) +# ------------------------------------------------------------------------- +#MLog.VerbosityLevel: 2 +#MLog.DebugLevel: 1 +#MLog.NoColors: yes + +# ========================================================================== +############################################################################# +# ========================================================================== +# Calibration +# ========================================================================== +############################################################################# +# ========================================================================== + +# ------------------------------------------------------------------------- +# Use this if you want to write the MJPedestal output somewhere +# If you don't want it, it is written to the calibration output anyhow. +# ------------------------------------------------------------------------- +#MJPedestalC1.PathOut: . +#MJPedestalC2.PathOut: . +#MJPedestalC1.BadPixelsFile: mjobs/badpixels_0_559.rc +#MJPedestalC1.ReferenceFile: mjobs/pedestalref.rc +#MJPedestalC2.ReferenceFile: mjobs/pedestalref.rc +MJCalibration.ReferenceFile: mjobs/calibrationref_Nov04.rc +MJCalibration.MHCalibrationRelTimeCam.ReferenceFile: mjobs/calibrationref_Nov04.rc +MJCalibration.MHCalibrationChargeCam.ReferenceFile: mjobs/calibrationref_Nov04.rc + +# ------------------------------------------------------------------------- +# Use this to define where the calibration output is stored. The filename +# is created from the sequence number. If nothing is specified '.' is +# assumed. (overwrites command line options) +# ------------------------------------------------------------------------- +#MJCalibration.PathOut: calped + +# ------------------------------------------------------------------------- +# Define here which signal extractor you are going to use for the +# whole calibration chain +# ------------------------------------------------------------------------- +MJPedestalC1.ExtractSignal: MExtractTimeAndChargeDigitalFilter +# ------------------------------------------------------------------------- +# Define here an extractor which can be used for the December 04 data. +# ------------------------------------------------------------------------- +#MJPedestalC1.ExtractSignal: MExtractTimeAndChargeSpline, MExtractTimeAndChargeDigitalFilterPeakSearch, MExtractTimeAndChargeSlidingWindow +# ------------------------------------------------------------------------- +# Define here parameters valid for all extractors above +# The actual defaults used depend on the extractor in use. +# For details see the corresponding source code +# ------------------------------------------------------------------------- +#MJPedestalC1.ExtractSignal.HiGainFirst: 0 +#MJPedestalC1.ExtractSignal.HiGainLast: 15 +#MJPedestalC1.ExtractSignal.LoGainFirst: 3 +#MJPedestalC1.ExtractSignal.LoGainLast: 14 +#MJPedestalC1.ExtractSignal.SaturationLimit: 250 +#MJPedestalC1.ExtractSignal.OffsetLoGain: 0.95 +#MJPedestalC1.ExtractSignal.LoGainSwitch: 120 +#MJPedestalC1.ExtractSignal.LoGainStartShift: -2.5 +# ------------------------------------------------------------------------- +# Define here parameters valid for the spline: +# ------------------------------------------------------------------------- +#MJPedestalC1.ExtractSignal.Resolution: 0.05 +#MJPedestalC1.ExtractSignal.RiseTimeHiGain: 0.5 +#MJPedestalC1.ExtractSignal.FallTimeHiGain: 1.5 +#MJPedestalC1.ExtractSignal.LoGainStretch: 1.5 +#MJPedestalC1.ExtractSignal.ExtractionType: Integral +#MJPedestalC1.ExtractSignal.ExtractionType: Amplitude + +# ------------------------------------------------------------------------- +# In the case of the Digital filter you have to overwrite the default +# weights file depending on what you are extracting +# ------------------------------------------------------------------------- +MJPedestalC2.ExtractSignal.WeightsFile: msignal/calibration_weights_UV46.dat +MJPedestalY2.ExtractSignal.WeightsFile: msignal/cosmics_weights46.dat +MJPedestalY3.ExtractSignal.WeightsFile: msignal/cosmics_weights46.dat +MJCalibration.ExtractSignal.WeightsFile: msignal/calibration_weights_UV46.dat +MJCalibrateSignal.ExtractSignal.WeightsFile: msignal/cosmics_weights46.dat +MJCalibrateSignal.ExtractInterlaced.WeightsFile: msignal/calibration_weights_UV46.dat + +#MJPedestalC2.ExtractSignal.AutomaticWeightsFile: On +#MJPedestalY2.ExtractSignal.AutomaticWeightsFile: On +#MJPedestalY3.ExtractSignal.AutomaticWeightsFile: On +#MJCalibration.ExtractSignal.AutomaticWeightsFile: On +#MJCalibrateSignal.ExtractSignal.AutomaticWeightsFile: On +#MJCalibrateSignal.ExtractInterlaced.AutomaticWeightsFile: On + +# ------------------------------------------------------------------------- +# Configure MJCalibration +# ------------------------------------------------------------------------- + +# Switch on relative time calibration +#MJCalibration.RelTimeCalibration: Yes +# Type of used data format: raw, root, mc +#MJCalibration.DataType: Root +# Type of displayed plots: Full,DataCheck,Normal +#MJCalibration.Display: datacheck +# Write additinal debug output +#MJCalibration.Debug: No +# Use blind pixel +#MJCalibration.UseBlindPixel: No +# Use pin diode +#MJCalibration.UsePINDiode: No +# Use pixel check option +#MJCalibration.PixelCheck: Yes +#MJCalibTest.PixelCheck: Yes + +# ------------------------------------------------------------------------- +# Individual events check for calibration results +# ------------------------------------------------------------------------- + +#ContCosmics.MaxEmptyPixels: 0.2 +#ContCosmics.MaxExcludedFraction: 1 +#ContCosmics.MinAcceptedFraction: 0 + +# ------------------------------------------------------------------------- +# You can choose the pedestal extraction algorithm/task. To use +# MPedCalcPedRun use the lines below. Be carefull, a pedestal file +# could be (if no pedestal file available) a data-file, too. +# ------------------------------------------------------------------------- +#MJPedestalC1.ExtractPedestal: MPedCalcPedRun +#MJPedestalC1.ExtractPedestal.ExtractWinFirst: 0 +#MJPedestalC1.ExtractPedestal.ExtractWinSize: 6 +#MJPedestalC1.ExtractPedestal.PedestalUpdate: no +#MJPedestalC1.ExtractPedestal.CheckWinFirst: 0 +#MJPedestalC1.ExtractPedestal.CheckWinLast: 6 +#MJPedestalC1.ExtractPedestal.MaxSignalVar: 300 + +#MJPedestalC2.ExtractPedestal: MPedCalcPedRun +#MJPedestalC2.ExtractPedestal.ExtractWinFirst: 0 +#MJPedestalC2.ExtractPedestal.ExtractWinSize: 6 +#MJPedestalC2.ExtractPedestal.PedestalUpdate: no +#MJPedestalC2.ExtractPedestal.CheckWinFirst: 0 +#MJPedestalC2.ExtractPedestal.CheckWinSize: 29 + +# ------------------------------------------------------------------------- +# You can choose the pedestal extraction algorithm/task. To use +# MPedCalcFromLoGain use the lines below. Be carefull, a pedestal file +# could be (if no pedestal file available) a data-file, too. +# ------------------------------------------------------------------------- +#MJPedestalC1.ExtractPedestal: MPedCalcFromLoGain +#MJPedestalC1.ExtractPedestal.ExtractWinFirst: 0 +#MJPedestalC1.ExtractPedestal.ExtractWinSize: 6 +#MJPedestalC1.ExtractPedestal.PedestalUpdate: no +#MJPedestalC1.ExtractPedestal.CheckWinFirst: 0 +#MJPedestalC1.ExtractPedestal.CheckWinSize: 29 +#MJPedestalC1.ExtractPedestal.MaxSignalVar: 40 + +#MJPedestalC2.ExtractPedestal: MPedCalcFromLoGain +#MJPedestalC2.ExtractPedestal.ExtractWinFirst: 0 +#MJPedestalC2.ExtractPedestal.ExtractWinSize: 6 +#MJPedestalC2.ExtractPedestal.PedestalUpdate: no +#MJPedestalC2.ExtractPedestal.CheckWinFirst: 0 +#MJPedestalC2.ExtractPedestal.CheckWinSize: 29 +#MJPedestalC2.ExtractPedestal.MaxSignalVar: 40 + +# ------------------------------------------------------------------------- +# Configure MJPedestalC1 and MJPedestalC2 +# ------------------------------------------------------------------------- +# Possible data type: Raw,Root,MC +#MJPedestalC1.DataType: Root +#MJPedestalC2.DataType: Root +#MJPedestalC3.DataType: Root +# Type of displayed plots: Full,DataCheck,None +#MJPedestalC1.Display: datacheck +#MJPedestalC2.Display: datacheck +#MJPedestalC3.Display: datacheck +# Maximum number of event processed in the loop +#MJPedestalC1.MaxEvents: 1000 +#MJPedestalC2.MaxEvents: 1000 +# Allow to overwrite existing files with the output file +#MJPedestalC1.AllowOverwrite: No +#MJPedestalC2.AllowOverwrite: No +# Use data runs from the sequence instead of calibration runs +#MJPedestalC1.UseData: No +#MJPedestalC2.UseData: No +# Switch off DeadPixelCheck with "No" +#MJPedestalY2.DeadPixelCheck: yes + + +# ------------------------------------------------------------------------- +# Use this if you want to change the software low-gain switch in the +# calibration +# ------------------------------------------------------------------------- +MJCalibration.ExtractSignal.LoGainSwitch: 170 + +# ------------------------------------------------------------------------- +# Use this if you want to change the higain-vs. logain intercalibration +# ------------------------------------------------------------------------- +#MJCalibration.HiLoCalibFile: mjobs/hilocalib_df46.root + +# ------------------------------------------------------------------------- +# Use this to change the behaviour of the calibration +# ------------------------------------------------------------------------- +# Type if you set a colour explicitely from outside (TEST purposes!!) +#MJCalibration.MCalibColorSet.ExplicitColor: green,blue,uv,ct1 + +#MJCalibration.MCalibrationChargeCalc.ArrTimeRmsLimit: 5.0 +#MJCalibration.MCalibrationChargeCalc.ChargeLimit: 2.5 +#MJCalibration.MCalibrationChargeCalc.ChargeErrLimit: 0 +#MJCalibration.MCalibrationChargeCalc.ChargeRelErrLimit: 1 +#MJCalibration.MCalibrationChargeCalc.Debug: no + +#MJCalibration.MCalibrationChargeCalc.FFactorErrLimit: 4.5 +#MJCalibration.MCalibrationChargeCalc.LambdaErrLimit: 0.2 +#MJCalibration.MCalibrationChargeCalc.LambdaCheckLimit: 0.5 +#MJCalibration.MCalibrationChargeCalc.PheErrLowerLimit: 9.0 +#MJCalibration.MCalibrationChargeCalc.PheErrUpperLimit: 5.5 + +#MJCalibration.MHCalibrationChargeCam.Debug: no +#MJCalibration.MHCalibrationChargeCam.LoGain: yes +#MJCalibration.MHCalibrationChargeCam.Oscillations: yes +#MJCalibration.MHCalibrationChargeCam.SizeCheck: yes +#MJCalibration.MHCalibrationChargeCam.Averageing: yes +#MJCalibration.MHCalibrationChargeCam.HiGainNbins: 500 +#MJCalibration.MHCalibrationChargeCam.HiGainFirst: -98. +#MJCalibration.MHCalibrationChargeCam.HiGainLast: 1902. +#MJCalibration.MHCalibrationChargeCam.LoGainNbins: 500 +#MJCalibration.MHCalibrationChargeCam.LoGainFirst: -99. +#MJCalibration.MHCalibrationChargeCam.LoGainLast: 901. +#MJCalibration.MHCalibrationChargeCam.TimeLowerLimit: 1. +#MJCalibration.MHCalibrationChargeCam.TimeUpperLimit: 3. +#MJCalibration.MHCalibrationChargeCam.NumHiGainSaturationLimit: 0.02 +#MJCalibration.MHCalibrationChargeCam.NumLoGainSaturationLimit: 0.005 +#MJCalibration.MHCalibrationChargeCam.ProbLimit: 0.00000001 +#MJCalibration.MHCalibrationChargeCam.OverflowLimit: 0.005 +#MJCalibration.MHCalibrationChargeCam.PulserFrequency: 500 + +#MJCalibration.MHCalibrationRelTimeCam.Debug: no +#MJCalibration.MHCalibrationRelTimeCam.LoGain: no +#MJCalibration.MHCalibrationRelTimeCam.Oscillations: yes +#MJCalibration.MHCalibrationRelTimeCam.SizeCheck: yes +#MJCalibration.MHCalibrationRelTimeCam.Averageing: yes +#MJCalibration.MHCalibrationRelTimeCam.Nbins: 210 +#MJCalibration.MHCalibrationRelTimeCam.First: -5. +#MJCalibration.MHCalibrationRelTimeCam.Last: 10. +#MJCalibration.MHCalibrationRelTimeCam.NumHiGainSaturationLimit: 0.25 +#MJCalibration.MHCalibrationRelTimeCam.ProbLimit: 0.00001 +#MJCalibration.MHCalibrationRelTimeCam.OverflowLimit: 0.005 +#MJCalibration.MHCalibrationRelTimeCam.PulserFrequency: 500 + + +# ========================================================================== +############################################################################# +# ========================================================================== +# Signal extraction +# ========================================================================== +############################################################################# +# ========================================================================== + +# ------------------------------------------------------------------------- +# Use a OutputPath if you want to write the results to a file +# ------------------------------------------------------------------------- +#MJCalibrateSignal.PathOut: . +#MJPedestalY1.PathOut: . +#MJPedestalY2.PathOut: . +#MJPedestalY3.PathOut: . +#MJCalibrateSignal.PathIn: . + +# ------------------------------------------------------------------------- +# If you want to change the behaviour of the pedestal extraction... +# ------------------------------------------------------------------------- +# Possible data type: Raw,Root,MC +#MJPedestalY1.DataType: Root +#MJPedestalY2.DataType: Root +#MJPedestalY3.DataType: Root +#MJPedestalY1.UseData: Yes +#MJPedestalY2.UseData: Yes +#MJPedestalY3.UseData: Yes +#MJPedestalY1.MinEvents: 50 +#MJPedestalY2.MinEvents: 50 +#MJPedestalY3.MinEvents: 50 + +MJPedestalY1.MaxEvents: 1000 +MJPedestalY2.MaxEvents: 10000 +MJPedestalY3.MaxEvents: 1000 + +MJPedestalY1.MaxPedestals: 500 +MJPedestalY2.MaxPedestals: 500 +MJPedestalY3.MaxPedestals: 500 + +# ------------------------------------------------------------------------- +# When the starting pedestal is extracted from the first events of +# the data run also the average pulse position of these events +# is determined. If the pulse position check is switched on (default) +# the extraction window of the hi- and lo-gain extractor for the +# determined average pulse position is not allowed to exceed the +# extractors search region. In this scenario we assume that the +# start of the extraction window is at least similar to the +# rising edge of the pulse. +# To make sure that the search region and the extraction window +# does not touch increase the values. +# ------------------------------------------------------------------------- +#MJPedestalY2.PulsePosCheck: yes +#MJPedestalY2.ExtractWinLeft: 0.0 +#MJPedestalY2.ExtractWinRight: 0.0 + + +# ------------------------------------------------------------------------- +# Define the Pulse Position check parameters: +# ------------------------------------------------------------------------- +#MJPedestalY2.MHCalibrationPulseTimeCam.SaturationLimit: 245 +#MJPedestalY2.MHCalibrationPulseTimeCam.LowerSignalLimit: 85 +#MJPedestalY2.MHCalibrationPulseTimeCam.NumPixelsRequired: 3 +#MJPedestalY2.PixelCheck: no + +#MJPedestalY.ExtractPedestal: MPedCalcFromLoGain +#MJPedestalY.ExtractPedestal.PedestalUpdate: no + +# ------------------------------------------------------------------------- +# Setup pedestal extractor for running through data +# ------------------------------------------------------------------------- +#MJCalibrateSignal.MPedCalcFundamental.PedestalUpdate: yes +#MJCalibrateSignal.MPedCalcFundamental.NumDump: 500 +#MJCalibrateSignal.MPedCalcFundamental.NumEventsDump: 500 +#MJCalibrateSignal.MPedCalcFundamental.NumAreasDump: 500 +#MJCalibrateSignal.MPedCalcFundamental.NumSectorsDump: 500 +#MJCalibrateSignal.MPedCalcFundamental.ExtractWinFirst: 17 +#MJCalibrateSignal.MPedCalcFundamental.ExtractWinSize: 6 +#MJCalibrateSignal.MPedCalcFundamental.CheckWinFirst: 0 +#MJCalibrateSignal.MPedCalcFundamental.CheckWinSize: 29 +#MJCalibrateSignal.MPedCalcFundamental.MaxSignalVar: 40 + +#MJCalibrateSignal.MPedCalcFromExtractorRndm.PedestalUpdate: yes +#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumEventsDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumAreasDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumSectorsDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.ExtractWinFirst: 17 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.ExtractWinSize: 6 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.CheckWinFirst: 0 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.CheckWinSize: 29 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.MaxSignalVar: 40 + +#MJCalibrateSignal.MPedCalcFromExtractor.PedestalUpdate: yes +#MJCalibrateSignal.MPedCalcFromExtractor.NumDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractor.NumEventsDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractor.NumAreasDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractor.NumSectorsDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractor.ExtractWinFirst: 17 +#MJCalibrateSignal.MPedCalcFromExtractor.ExtractWinSize: 6 +#MJCalibrateSignal.MPedCalcFromExtractor.CheckWinFirst: 0 +#MJCalibrateSignal.MPedCalcFromExtractor.CheckWinSize: 29 +#MJCalibrateSignal.MPedCalcFromExtractor.MaxSignalVar: 40 + +# ------------------------------------------------------------------------- +# Configure the interlaced calibration updates +# ------------------------------------------------------------------------- +#MJCalibrateSignal.MCalibCalcFromPast.UpdateWithFFactorMethod: yes +#MJCalibrateSignal.MCalibCalcFromPast.NumEventsDump: 500 +#MJCalibrateSignal.MCalibCalcFromPast.UpdateNumPhes: yes +#MJCalibrateSignal.MCalibCalcFromPast.NumPhesDump: 10 +#MJCalibrateSignal.MCalibrationChargeCalc.Debug: no + +# ------------------------------------------------------------------------- +# Use this if you want to change the signal extractor +# ------------------------------------------------------------------------- +MJCalibrateSignal.ExtractSignal.LoGainSwitch: 150 +#MJCalibrateSignal.ExtractSignal.LoGainStretch: 1.5 +#MJCalibrateSignal.ExtractInterlaced.LoGainStretch: 1.5 + +# ------------------------------------------------------------------------- +# Setup calibration of data +# ------------------------------------------------------------------------- +# Type of used data format: raw,root,MC +#MJCalibrateSignal.DataType: Root +# Type if you set a colour explicitely from outside (only TESTS!!!) +#MJCalibrateSignal.MCalibColorSet.ExpicitColor: green,blue,uv,ct1 +#MJCalibrateSignal.MCalibrateData.PedestalFlag: Event +#MJCalibrateSignal.MCalibrateData.ScaleFactor: 1 +MJCalibrateSignal.MCalibrateData.FileNameScale: resources/calibration.rc +#MJCalibrateSignal.MCalibrateData.CalibrationMode: Default +#MJCalibrateSignal.MCalibrateData.CalibConvMinLimit: 0.01 +#MJCalibrateSignal.MCalibrateData.CalibConvMaxLimit: 5. +#MJCalibrateSignal.Interlaced: yes +#MJCalibrateSignal.RelTimesUpdate: yes + +# ------------------------------------------------------------------------- +# Setup level for determining of bad pixels +# ------------------------------------------------------------------------- +#MJCalibrateSignal.MBadPixelsCalc.PedestalLevel: 3.0 +#MJCalibrateSignal.MBadPixelsCalc.PedestalLevelVariance: 3.0 +#MJCalibrateSignal.MBadPixelsTreat.NumMinNeighbors: 3 +#MJCalibrateSignal.MBadPixelsTreat.MaxArrivalTimeDiff: 3.0 +#MJCalibrateSignal.MBadPixelsTreat.UseInterpolation: yes +#MJCalibrateSignal.MBadPixelsTreat.ProcessPedestalEvt: yes +#MJCalibrateSignal.MBadPixelsTreat.ProcessPedestalRun: no +#MJCalibrateSignal.MBadPixelsTreat.ProcessTimes: yes +#MJCalibrateSignal.MBadPixelsTreat.UseCentralPixel: no +#MJCalibrateSignal.MBadPixelsTreat.HardTreatment: no + +# ------------------------------------------------------------------------- +# Setup interlaced calibration events histogramming +# ------------------------------------------------------------------------- +#MJCalibrateSignal.MHCalibrationChargeCam.Debug: no +#MJCalibrateSignal.MHCalibrationChargeCam.LoGain: yes +MJCalibrateSignal.MHCalibrationChargeCam.Oscillations: no +#MJCalibrateSignal.MHCalibrationChargeCam.SizeCheck: yes +#MJCalibrateSignal.MHCalibrationChargeCam.Averageing: yes +#MJCalibrateSignal.MHCalibrationChargeCam.HiGainNbins: 500 +#MJCalibrateSignal.MHCalibrationChargeCam.HiGainFirst: -98. +#MJCalibrateSignal.MHCalibrationChargeCam.HiGainLast: 1902. +#MJCalibrateSignal.MHCalibrationChargeCam.HiGainLast: 1899.5 +MJCalibrateSignal.MHCalibrationChargeCam.LoGainNbins: 250 +MJCalibrateSignal.MHCalibrationChargeCam.LoGainFirst: -98 +MJCalibrateSignal.MHCalibrationChargeCam.LoGainLast: 902 +#MJCalibrateSignal.MHCalibrationChargeCam.TimeLowerLimit: 1. +#MJCalibrateSignal.MHCalibrationChargeCam.TimeUpperLimit: 3. +#MJCalibrateSignal.MHCalibrationChargeCam.NumHiGainSaturationLimit: 0.02 +#MJCalibrateSignal.MHCalibrationChargeCam.NumLoGainSaturationLimit: 0.005 +#MJCalibrateSignal.MHCalibrationChargeCam.ProbLimit: 0.0001 +#MJCalibrateSignal.MHCalibrationChargeCam.OverflowLimit: 0.005 +MJCalibrateSignal.MHCalibrationChargeCam.PulserFrequency: 50 + +# ------------------------------------------------------------------------- +# Setup relative time interlaced calibration +# ------------------------------------------------------------------------- +#MJCalibrateSignal.MHCalibrationRelTimeCam.Debug: no +#MJCalibrateSignal.MHCalibrationRelTimeCam.LoGain: no +MJCalibrateSignal.MHCalibrationRelTimeCam.Oscillations: no +#MJCalibrateSignal.MHCalibrationRelTimeCam.SizeCheck: yes +#MJCalibrateSignal.MHCalibrationRelTimeCam.Averageing: yes +#MJCalibrateSignal.MHCalibrationRelTimeCam.Nbins: 210 +#MJCalibrateSignal.MHCalibrationRelTimeCam.First: -5. +#MJCalibrateSignal.MHCalibrationRelTimeCam.Last: 10. +#MJCalibrateSignal.MHCalibrationRelTimeCam.NumHiGainSaturationLimit: 0.25 +#MJCalibrateSignal.MHCalibrationRelTimeCam.ProbLimit: 0.00001 +#MJCalibrateSignal.MHCalibrationRelTimeCam.OverflowLimit: 0.005 +MJCalibrateSignal.MHCalibrationRelTimeCam.PulserFrequency: 50 + +# ------------------------------------------------------------------------- +# Resource for movie encoding. For more details see IACT wiki +# ------------------------------------------------------------------------- +#MMoviePrepare.FirstSlice: 0 +#MMoviePrepare.LastSlice: 15 +#MovieFilter.Condition: MMovieData.fMax>5*MMovieData.fMedianPedestalRms +#MImgCleanStd.CleanLevel1: 8.0 +#MImgCleanStd.CleanLevel2: 4.5 +#MFEvtNumber.FileName: ganymed00001218.root +#MFEvtNumber.Selector: ThetaSquared.fVal<0.04 && DataType.fVal>0.5 +#MMovieWrite.TargetLength: 5 +#MMovieWrite.NumEvents: 10 +#MMovieWrite.Threshold: 3 +#MMovieWrite.Filename: movie.mpg + +# ------------------------------------------------------------------------- +# This is a special option for the datacheck. A task executed as last +# task in the event processing task list - in the standard software it +# is skipped +# ------------------------------------------------------------------------- +#MJCalibrateSignal.FinalFantasy: MHCamFlorian Index: trunk/Mars/resources/callisto_Dec04Jan05.rc =================================================================== --- trunk/Mars/resources/callisto_Dec04Jan05.rc (revision 19872) +++ trunk/Mars/resources/callisto_Dec04Jan05.rc (revision 19872) @@ -0,0 +1,483 @@ +# ========================================================================== +############################################################################# +# ========================================================================== +# General +# ========================================================================== +############################################################################# +# ========================================================================== + +# ------------------------------------------------------------------------- +# Use this if you want to setup the logging stream for the jobs +# (overwrites command line options) +# ------------------------------------------------------------------------- +#MLog.VerbosityLevel: 2 +#MLog.DebugLevel: 1 +#MLog.NoColors: yes + +# ========================================================================== +############################################################################# +# ========================================================================== +# Calibration +# ========================================================================== +############################################################################# +# ========================================================================== + +# ------------------------------------------------------------------------- +# Use this if you want to write the MJPedestal output somewhere +# If you don't want it, it is written to the calibration output anyhow. +# ------------------------------------------------------------------------- +#MJPedestalC1.PathOut: . +#MJPedestalC2.PathOut: . + +# ------------------------------------------------------------------------- +# Use this if you want to modify the exclusion of pixels from the beginning +# (default: Pixel 0 and 559) +# ------------------------------------------------------------------------- +MJPedestalC1.BadPixelsFile: mjobs/badpixels_0_559.rc + +# ------------------------------------------------------------------------- +# Use this if you want to modify files containing information about the +# reference lines in the status display +# ------------------------------------------------------------------------- +MJPedestalC1.ReferenceFile: mjobs/pedestalref_Dec04.rc +MJPedestalC2.ReferenceFile: mjobs/pedestalref_Dec04.rc +MJCalibration.ReferenceFile: mjobs/calibrationref_Dec04.rc +MJCalibration.MHCalibrationRelTimeCam.ReferenceFile: mjobs/calibrationref_Dec04.rc +MJCalibration.MHCalibrationChargeCam.ReferenceFile: mjobs/calibrationref_Dec04.rc + +# ------------------------------------------------------------------------- +# Use this to define where the calibration output is stored. The filename +# is created from the sequence number. If nothing is specified '.' is +# assumed. (overwrites command line options) +# ------------------------------------------------------------------------- +#MJCalibration.PathOut: calped + +# ------------------------------------------------------------------------- +# Use this to define where the program should search for the pedestal +# and calibration files defined in the sequence. To use the local +# directory use '.' If nothing is specified the default path in the +# datacenter is used. (overwrites command line options) +# ------------------------------------------------------------------------- +#MJPedestalC1.PathData: /data/MAGIC/Period016/ +#MJPedestalC2.PathData: /data/MAGIC/Period016/ +#MJCalibration.PathData: /data/MAGIC/Period016/ + +# ------------------------------------------------------------------------- +# Define here which signal extractor you are going to use for the +# whole calibration chain +# ------------------------------------------------------------------------- +MJPedestalC1.ExtractSignal: MExtractTimeAndChargeSpline +# ------------------------------------------------------------------------- +# Possibilities for a different extractor +# ------------------------------------------------------------------------- +#MJPedestalC1.ExtractSignal: MExtractTimeAndChargeDigitalFilter, MExtractTimeAndChargeDigitalFilterPeakSearch, MExtractTimeAndChargeSlidingWindow +# ------------------------------------------------------------------------- +# Define here parameters valid for all extractors above +# The actual defaults used depend on the extractor in use. +# For details see the corresponding source code +# ------------------------------------------------------------------------- +#MJPedestalC1.ExtractSignal.HiGainFirst: 0 +#MJPedestalC1.ExtractSignal.HiGainLast: 15 +#MJPedestalC1.ExtractSignal.LoGainFirst: 3 +#MJPedestalC1.ExtractSignal.LoGainLast: 14 +#MJPedestalC1.ExtractSignal.SaturationLimit: 250 +#MJPedestalC1.ExtractSignal.OffsetLoGain: 1.3 +#MJPedestalC1.ExtractSignal.LoGainSwitch: 120 +#MJPedestalC1.ExtractSignal.LoGainStartShift: -2.5 +# ------------------------------------------------------------------------- +# Define here parameters valid for the spline: +# ------------------------------------------------------------------------- +#MJPedestalC1.ExtractSignal.Resolution: 0.05 +MJPedestalC1.ExtractSignal.RiseTimeHiGain: 0.5 +MJPedestalC1.ExtractSignal.FallTimeHiGain: 0.5 +#MJPedestalC1.ExtractSignal.LoGainStretch: 1.5 +#MJPedestalC1.ExtractSignal.ExtractionType: Integral +#MJPedestalC1.ExtractSignal.ExtractionType: Amplitude + +# ------------------------------------------------------------------------- +# In the case of the Digital filter you have to overwrite the default +# weights file depending on what you are extracting +# ------------------------------------------------------------------------- +#MJPedestalC2.ExtractSignal.WeightsFile: msignal/calibration_weights_UV4.dat +#MJPedestalY2.ExtractSignal.WeightsFile: msignal/cosmics_weights4.dat +#MJPedestalY3.ExtractSignal.WeightsFile: msignal/cosmics_weights4.dat +#MJCalibration.ExtractSignal.WeightsFile: msignal/calibration_weights_UV4.dat +#MJCalibrateSignal.ExtractSignal.WeightsFile: msignal/cosmics_weights4.dat +#MJCalibrateSignal.ExtractInterlaced.WeightsFile: msignal/calibration_weights_UV4.dat +#MJCalibrateSignal.ExtractSignal.LoGainStretch: 1.5 +#MJCalibrateSignal.ExtractInterlaced.LoGainStretch: 1.5 + +# ------------------------------------------------------------------------- +# Configure MJCalibration +# ------------------------------------------------------------------------- + +# Switch on relative time calibration +#MJCalibration.RelTimeCalibration: Yes +# Set color to be used +#MJCalibration.Color: +# Type of used data format: raw, root, mc +#MJCalibration.DataType: Root +# Type of displayed plots: Full,DataCheck,Normal +#MJCalibration.Display: datacheck +# Write additinal debug output +#MJCalibration.Debug: No +# Use blind pixel +#MJCalibration.UseBlindPixel: No +# Use pin diode +#MJCalibration.UsePINDiode: No +# Use pixel check option +#MJCalibration.PixelCheck: No +#MJCalibTest.PixelCheck: No + +# ------------------------------------------------------------------------- +# Individual events check for calibration results +# ------------------------------------------------------------------------- + +#ContCosmics.MaxEmptyPixels: 0.2 +#ContCosmics.MaxExcludedFraction: 1 +#ContCosmics.MinAcceptedFraction: 0 + +# ------------------------------------------------------------------------- +# In case of MOON DATA, uncomment the following lines: +# (better: call callisto with option "-moon") +# ------------------------------------------------------------------------- +#MJPedestalC1.UseHists: yes +#MJPedestalC2.UseHists: yes +#MJPedestalY1.UseHists: yes +#MJPedestalY2.UseHists: yes +#MJPedestalY3.UseHists: yes + +# ------------------------------------------------------------------------- +# You can choose the pedestal extraction algorithm/task. To use +# MPedCalcPedRun use the lines below. Be carefull, a pedestal file +# could be (if no pedestal file available) a data-file, too. +# ------------------------------------------------------------------------- +#MJPedestalC1.ExtractPedestal: MPedCalcPedRun +#MJPedestalC1.ExtractPedestal.ExtractWinFirst: 0 +#MJPedestalC1.ExtractPedestal.ExtractWinSize: 6 +#MJPedestalC1.ExtractPedestal.PedestalUpdate: no +#MJPedestalC1.ExtractPedestal.CheckWinFirst: 0 +#MJPedestalC1.ExtractPedestal.CheckWinLast: 29 +#MJPedestalC1.ExtractPedestal.MaxSignalVar: 40 + +#MJPedestalC2.ExtractPedestal: MPedCalcPedRun +#MJPedestalC2.ExtractPedestal.ExtractWinFirst: 0 +#MJPedestalC2.ExtractPedestal.ExtractWinSize: 6 +#MJPedestalC2.ExtractPedestal.PedestalUpdate: no +#MJPedestalC2.ExtractPedestal.CheckWinFirst: 0 +#MJPedestalC2.ExtractPedestal.CheckWinLast: 29 +#MJPedestalC2.ExtractPedestal.MaxSignalVar: 40 + +# ------------------------------------------------------------------------- +# You can choose the pedestal extraction algorithm/task. To use +# MPedCalcFromLoGain use the lines below. Be carefull, a pedestal file +# could be (if no pedestal file available) a data-file, too. +# ------------------------------------------------------------------------- +#MJPedestalC1.ExtractPedestal: MPedCalcFromLoGain +#MJPedestalC1.ExtractPedestal.ExtractWinFirst: 0 +#MJPedestalC1.ExtractPedestal.ExtractWinSize: 6 +#MJPedestalC1.ExtractPedestal.PedestalUpdate: no +#MJPedestalC1.ExtractPedestal.CheckWinFirst: 0 +#MJPedestalC1.ExtractPedestal.CheckWinLast: 29 +#MJPedestalC1.ExtractPedestal.MaxSignalVar: 40 + +#MJPedestalC2.ExtractPedestal: MPedCalcFromLoGain +#MJPedestalC2.ExtractPedestal.ExtractWinFirst: 0 +#MJPedestalC2.ExtractPedestal.ExtractWinSize: 6 +#MJPedestalC2.ExtractPedestal.PedestalUpdate: no +#MJPedestalC2.ExtractPedestal.CheckWinFirst: 0 +#MJPedestalC2.ExtractPedestal.CheckWinLast: 29 +#MJPedestalC2.ExtractPedestal.MaxSignalVar: 40 + +# ------------------------------------------------------------------------- +# Configure MJPedestalC1 and MJPedestalC2 +# ------------------------------------------------------------------------- +# Possible data type: Raw,Root,MC +#MJPedestalC1.DataType: Root +#MJPedestalC2.DataType: Root +#MJPedestalC3.DataType: Root +# Type of displayed plots: Full,DataCheck,None +#MJPedestalC1.Display: datacheck +#MJPedestalC2.Display: datacheck +#MJPedestalC3.Display: datacheck +# Maximum number of event processed in the loop +#MJPedestalC1.MaxEvents: 1000 +#MJPedestalC2.MaxEvents: 1000 +# Allow to overwrite existing files with the output file +#MJPedestalC1.AllowOverwrite: No +#MJPedestalC2.AllowOverwrite: No +# Use data runs from the sequence instead of calibration runs +#MJPedestalC1.UseData: No +#MJPedestalC2.UseData: No + + +# ------------------------------------------------------------------------- +# Use this if you want to change the software low-gain switch +# ------------------------------------------------------------------------- +#MJCalibration.ExtractSignal.LoGainSwitch: 150 + +# ------------------------------------------------------------------------- +# Use this if you want to change the higain-vs. logain intercalibration +# ------------------------------------------------------------------------- +MJCalibration.HiLoCalibFile: resources/hilocalib_sp1.rc + +# ------------------------------------------------------------------------- +# Use this to change the behaviour of the calibration +# ------------------------------------------------------------------------- +# Type if you set a colour explicitely from outside (only for MC!!!) +#MJCalibration.MCalibColorSet.ExplicitColor: green,blue,uv,ct1 + +#MJCalibration.MCalibrationChargeCalc.ArrTimeRmsLimit: 5.0 +#MJCalibration.MCalibrationChargeCalc.ChargeLimit: 2.5 +#MJCalibration.MCalibrationChargeCalc.ChargeErrLimit: 0 +#MJCalibration.MCalibrationChargeCalc.ChargeRelErrLimit: 1 +#MJCalibration.MCalibrationChargeCalc.Debug: no + +#MJCalibration.MCalibrationChargeCalc.FFactorErrLimit: 4.5 +#MJCalibration.MCalibrationChargeCalc.LambdaErrLimit: 0.2 +#MJCalibration.MCalibrationChargeCalc.LambdaCheckLimit: 0.5 +#MJCalibration.MCalibrationChargeCalc.PheErrLowerLimit: 9.0 +#MJCalibration.MCalibrationChargeCalc.PheErrUpperLimit: 5.5 + +#MJCalibration.MHCalibrationChargeCam.Debug: no +#MJCalibration.MHCalibrationChargeCam.LoGain: yes +#MJCalibration.MHCalibrationChargeCam.Oscillations: yes +#MJCalibration.MHCalibrationChargeCam.SizeCheck: yes +#MJCalibration.MHCalibrationChargeCam.Averageing: yes +#MJCalibration.MHCalibrationChargeCam.HiGainNbins: 500 +#MJCalibration.MHCalibrationChargeCam.HiGainFirst: -100.125 +#MJCalibration.MHCalibrationChargeCam.HiGainLast: 1899.875 +#MJCalibration.MHCalibrationChargeCam.LoGainNbins: 500 +#MJCalibration.MHCalibrationChargeCam.LoGainFirst: -100.25 +#MJCalibration.MHCalibrationChargeCam.LoGainLast: 899.75 +#MJCalibration.MHCalibrationChargeCam.TimeLowerLimit: 1. +#MJCalibration.MHCalibrationChargeCam.TimeUpperLimit: 3. +#MJCalibration.MHCalibrationChargeCam.NumHiGainSaturationLimit: 0.02 +#MJCalibration.MHCalibrationChargeCam.NumLoGainSaturationLimit: 0.005 +MJCalibration.MHCalibrationChargeCam.ProbLimit: 1.0E-18 +#MJCalibration.MHCalibrationChargeCam.OverflowLimit: 0.005 +#MJCalibration.MHCalibrationChargeCam.PulserFrequency: 500 + +#MJCalibration.MHCalibrationRelTimeCam.Debug: no +#MJCalibration.MHCalibrationRelTimeCam.LoGain: no +#MJCalibration.MHCalibrationRelTimeCam.Oscillations: yes +#MJCalibration.MHCalibrationRelTimeCam.SizeCheck: yes +#MJCalibration.MHCalibrationRelTimeCam.Averageing: yes +#MJCalibration.MHCalibrationRelTimeCam.Nbins: 300 +#MJCalibration.MHCalibrationRelTimeCam.First: -4.975 +#MJCalibration.MHCalibrationRelTimeCam.Last: 10.025 +#MJCalibration.MHCalibrationRelTimeCam.NumHiGainSaturationLimit: 0.25 +#MJCalibration.MHCalibrationRelTimeCam.ProbLimit: 0.00001 +#MJCalibration.MHCalibrationRelTimeCam.OverflowLimit: 0.005 +#MJCalibration.MHCalibrationRelTimeCam.PulserFrequency: 500 + + +# ========================================================================== +############################################################################# +# ========================================================================== +# Signal extraction +# ========================================================================== +############################################################################# +# ========================================================================== + +# ------------------------------------------------------------------------- +# Use a OutputPath if you want to write the results to a file +# ------------------------------------------------------------------------- +#MJCalibrateSignal.PathOut: . +#MJPedestalY1.PathOut: . +#MJPedestalY2.PathOut: . +#MJPedestalY3.PathOut: . +#MJPedestalY1.PathData: /data/MAGIC/Period016/ +#MJPedestalY2.PathData: /data/MAGIC/Period016/ +#MJPedestalY3.PathData: /data/MAGIC/Period016/ +#MJCalibrateSignal.PathData: /data/MAGIC/Period016/ +#MJCalibrateSignal.PathIn: . + +# ------------------------------------------------------------------------- +# If you want to change the behaviour of the pedestal extraction... +# ------------------------------------------------------------------------- +# Possible data type: Raw,Root,MC +#MJPedestalY1.DataType: Root +#MJPedestalY2.DataType: Root +#MJPedestalY3.DataType: Root +#MJPedestalY1.UseData: Yes +#MJPedestalY2.UseData: Yes +#MJPedestalY3.UseData: Yes +#MJPedestalY1.MinEvents: 50 +#MJPedestalY2.MinEvents: 50 +#MJPedestalY3.MinEvents: 50 + +MJPedestalY1.MaxEvents: 1000 +MJPedestalY2.MaxEvents: 10000 +MJPedestalY3.MaxEvents: 1000 + +MJPedestalY1.MaxPedestals: 500 +MJPedestalY2.MaxPedestals: 500 +MJPedestalY3.MaxPedestals: 500 + +# ------------------------------------------------------------------------- +# When the starting pedestal is extracted from the first events of +# the data run also the average pulse position of these events +# is determined. If the pulse position check is switched on (default) +# the extraction window of the hi- and lo-gain extractor for the +# determined average pulse position is not allowed to exceed the +# extractors search region. In this scenario we assume that the +# start of the extraction window is at least similar to the +# rising edge of the pulse. +# To make sure that the search region and the extraction window +# does not touch increase the values. +# ------------------------------------------------------------------------- +#MJPedestalY2.PulsePosCheck: yes +#MJPedestalY2.ExtractWinLeft: 0.0 +#MJPedestalY2.ExtractWinRight: 0.0 + +# ------------------------------------------------------------------------- +# Define the Pulse Position check parameters: +# ------------------------------------------------------------------------- +#MJPedestalY2.MHCalibrationPulseTimeCam.SaturationLimit: 245 +#MJPedestalY2.MHCalibrationPulseTimeCam.LowerSignalLimit: 85 +#MJPedestalY2.MHCalibrationPulseTimeCam.NumPixelsRequired: 2 +#MJPedestalY2.PixelCheck: no + +#MJPedestalY.ExtractPedestal: MPedCalcFromLoGain +#MJPedestalY.ExtractPedestal.PedestalUpdate: no + +# ------------------------------------------------------------------------- +# Setup pedestal extractor for running through data +# ------------------------------------------------------------------------- +#MJCalibrateSignal.MPedCalcFundamental.PedestalUpdate: yes +#MJCalibrateSignal.MPedCalcFundamental.NumDump: 500 +#MJCalibrateSignal.MPedCalcFundamental.NumEventsDump: 500 +#MJCalibrateSignal.MPedCalcFundamental.NumAreasDump: 500 +#MJCalibrateSignal.MPedCalcFundamental.NumSectorsDump: 500 +#MJCalibrateSignal.MPedCalcFundamental.ExtractWinFirst: 17 +#MJCalibrateSignal.MPedCalcFundamental.ExtractWinSize: 6 +#MJCalibrateSignal.MPedCalcFundamental.CheckWinFirst: 0 +#MJCalibrateSignal.MPedCalcFundamental.CheckWinLast: 29 +#MJCalibrateSignal.MPedCalcFundamental.MaxSignalVar: 40 + +#MJCalibrateSignal.MPedCalcFromExtractorRndm.PedestalUpdate: yes +#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumEventsDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumAreasDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumSectorsDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.ExtractWinFirst: 17 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.ExtractWinSize: 6 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.CheckWinFirst: 0 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.CheckWinLast: 29 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.MaxSignalVar: 40 + +#MJCalibrateSignal.MPedCalcFromExtractor.PedestalUpdate: yes +#MJCalibrateSignal.MPedCalcFromExtractor.NumDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractor.NumEventsDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractor.NumAreasDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractor.NumSectorsDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractor.ExtractWinFirst: 17 +#MJCalibrateSignal.MPedCalcFromExtractor.ExtractWinSize: 6 +#MJCalibrateSignal.MPedCalcFromExtractor.CheckWinFirst: 0 +#MJCalibrateSignal.MPedCalcFromExtractor.CheckWinLast: 29 +#MJCalibrateSignal.MPedCalcFromExtractor.MaxSignalVar: 40 + +# ------------------------------------------------------------------------- +# Configure the interlaced calibration updates +# ------------------------------------------------------------------------- +#MJCalibrateSignal.MCalibCalcFromPast.UpdateWithFFactorMethod: yes +#MJCalibrateSignal.MCalibCalcFromPast.NumEventsDump: 500 +#MJCalibrateSignal.MCalibCalcFromPast.UpdateNumPhes: yes +#MJCalibrateSignal.MCalibCalcFromPast.NumPhesDump: 10 +#MJCalibrateSignal.MCalibrationChargeCalc.Debug: no + +# ------------------------------------------------------------------------- +# Use this if you want to change the signal extractor +# ------------------------------------------------------------------------- +MJCalibrateSignal.ExtractSignal.LoGainSwitch: 150 + +# ------------------------------------------------------------------------- +# Setup calibration of data +# ------------------------------------------------------------------------- +# Type of used data format: raw,root,MC +#MJCalibrateSignal.DataType: Root +# Type if you set a colour explicitely from outside (only for MC!!!) +#MJCalibrateSignal.MCalibColorSet.ExpicitColor: green,blue,uv,ct1 +#MJCalibrateSignal.MCalibrateData.PedestalFlag: Event +#MJCalibrateSignal.MCalibrateData.ScaleFactor: 0.85182358 +MJCalibrateSignal.MCalibrateData.FileNameScale: resources/calibration_spline.rc +#MJCalibrateSignal.MCalibrateData.CalibrationMode: Default +#MJCalibrateSignal.MCalibrateData.CalibConvMinLimit: 0.01 +#MJCalibrateSignal.MCalibrateData.CalibConvMaxLimit: 5. +#MJCalibrateSignal.Interlaced: yes +#MJCalibrateSignal.RelTimesUpdate: yes + +# ------------------------------------------------------------------------- +# Setup level for determining of bad pixels +# ------------------------------------------------------------------------- +#MJCalibrateSignal.MBadPixelsCalc.PedestalLevel: 3.0 +#MJCalibrateSignal.MBadPixelsCalc.PedestalLevelVariance: 3.0 +#MJCalibrateSignal.MBadPixelsTreat.NumMinNeighbors: 3 +#MJCalibrateSignal.MBadPixelsTreat.MaxArrivalTimeDiff: 3.0 +#MJCalibrateSignal.MBadPixelsTreat.UseInterpolation: yes +#MJCalibrateSignal.MBadPixelsTreat.ProcessPedestalEvt: yes +#MJCalibrateSignal.MBadPixelsTreat.ProcessPedestalRun: no +#MJCalibrateSignal.MBadPixelsTreat.ProcessTimes: yes +#MJCalibrateSignal.MBadPixelsTreat.UseCentralPixel: no +#MJCalibrateSignal.MBadPixelsTreat.HardTreatment: no + +# ------------------------------------------------------------------------- +# Setup interlaced calibration events histogramming +# ------------------------------------------------------------------------- +#MJCalibrateSignal.MHCalibrationChargeCam.Debug: no +#MJCalibrateSignal.MHCalibrationChargeCam.LoGain: yes +MJCalibrateSignal.MHCalibrationChargeCam.Oscillations: no +#MJCalibrateSignal.MHCalibrationChargeCam.SizeCheck: yes +#MJCalibrateSignal.MHCalibrationChargeCam.Averageing: yes +#MJCalibrateSignal.MHCalibrationChargeCam.HiGainNbins: 500 +#MJCalibrateSignal.MHCalibrationChargeCam.HiGainFirst: -100.5 +#MJCalibrateSignal.MHCalibrationChargeCam.HiGainLast: 1899.5 +MJCalibrateSignal.MHCalibrationChargeCam.LoGainNbins: 250 +MJCalibrateSignal.MHCalibrationChargeCam.LoGainFirst: -98 +MJCalibrateSignal.MHCalibrationChargeCam.LoGainLast: 902 +#MJCalibrateSignal.MHCalibrationChargeCam.TimeLowerLimit: 1. +#MJCalibrateSignal.MHCalibrationChargeCam.TimeUpperLimit: 3. +#MJCalibrateSignal.MHCalibrationChargeCam.NumHiGainSaturationLimit: 0.02 +#MJCalibrateSignal.MHCalibrationChargeCam.NumLoGainSaturationLimit: 0.005 +#MJCalibrateSignal.MHCalibrationChargeCam.ProbLimit: 0.0001 +#MJCalibrateSignal.MHCalibrationChargeCam.OverflowLimit: 0.005 +MJCalibrateSignal.MHCalibrationChargeCam.PulserFrequency: 50 + +# ------------------------------------------------------------------------- +# Setup relative time interlaced calibration +# ------------------------------------------------------------------------- +#MJCalibrateSignal.MHCalibrationRelTimeCam.Debug: no +#MJCalibrateSignal.MHCalibrationRelTimeCam.LoGain: no +MJCalibrateSignal.MHCalibrationRelTimeCam.Oscillations: no +#MJCalibrateSignal.MHCalibrationRelTimeCam.SizeCheck: yes +#MJCalibrateSignal.MHCalibrationRelTimeCam.Averageing: yes +#MJCalibrateSignal.MHCalibrationRelTimeCam.Nbins: 300 +#MJCalibrateSignal.MHCalibrationRelTimeCam.First: -4.975 +#MJCalibrateSignal.MHCalibrationRelTimeCam.Last: 10.025 +#MJCalibrateSignal.MHCalibrationRelTimeCam.NumHiGainSaturationLimit: 0.25 +#MJCalibrateSignal.MHCalibrationRelTimeCam.ProbLimit: 0.00001 +#MJCalibrateSignal.MHCalibrationRelTimeCam.OverflowLimit: 0.005 +MJCalibrateSignal.MHCalibrationRelTimeCam.PulserFrequency: 50 + +# ------------------------------------------------------------------------- +# Resource for movie encoding. For more details see IACT wiki +# ------------------------------------------------------------------------- +#MMoviePrepare.FirstSlice: 0 +#MMoviePrepare.LastSlice: 15 +#MovieFilter.Condition: MMovieData.fMax>5*MMovieData.fMedianPedestalRms +#MImgCleanStd.CleanLevel1: 8.0 +#MImgCleanStd.CleanLevel2: 4.5 +#MFEvtNumber.FileName: ganymed00001218.root +#MFEvtNumber.Selector: ThetaSquared.fVal<0.04 && DataType.fVal>0.5 +#MMovieWrite.TargetLength: 5 +#MMovieWrite.NumEvents: 500 +#MMovieWrite.Threshold: 3 +#MMovieWrite.FileName: filename.mpg + +# ------------------------------------------------------------------------- +# This is a special option for the datacheck. A task executed as last +# task in the event processing task list - in the standard software it +# is skipped +# ------------------------------------------------------------------------- +#MJCalibrateSignal.FinalFantasy: MHCamFlorian Index: trunk/Mars/resources/callisto_MarApr05.rc =================================================================== --- trunk/Mars/resources/callisto_MarApr05.rc (revision 19872) +++ trunk/Mars/resources/callisto_MarApr05.rc (revision 19872) @@ -0,0 +1,485 @@ +# ========================================================================== +############################################################################# +# ========================================================================== +# General +# ========================================================================== +############################################################################# +# ========================================================================== + +# ------------------------------------------------------------------------- +# Use this if you want to setup the logging stream for the jobs +# (overwrites command line options) +# ------------------------------------------------------------------------- +#MLog.VerbosityLevel: 2 +#MLog.DebugLevel: 1 +#MLog.NoColors: yes + +# ========================================================================== +############################################################################# +# ========================================================================== +# Calibration +# ========================================================================== +############################################################################# +# ========================================================================== + +# ------------------------------------------------------------------------- +# Use this if you want to write the MJPedestal output somewhere +# If you don't want it, it is written to the calibration output anyhow. +# ------------------------------------------------------------------------- +#MJPedestalC1.PathOut: . +#MJPedestalC2.PathOut: . +#MJPedestalC1.BadPixelsFile: mjobs/badpixels_0_559.rc +#MJPedestalC1.ReferenceFile: mjobs/pedestalref.rc +#MJPedestalC2.ReferenceFile: mjobs/pedestalref.rc +MJCalibration.ReferenceFile: mjobs/calibrationref_Nov04.rc +MJCalibration.MHCalibrationRelTimeCam.ReferenceFile: mjobs/calibrationref_Nov04.rc +MJCalibration.MHCalibrationChargeCam.ReferenceFile: mjobs/calibrationref_Nov04.rc + +# ------------------------------------------------------------------------- +# Use this to define where the calibration output is stored. The filename +# is created from the sequence number. If nothing is specified '.' is +# assumed. (overwrites command line options) +# ------------------------------------------------------------------------- +#MJCalibration.PathOut: calped + +# ------------------------------------------------------------------------- +# Use this to define where the program should search for the pedestal +# and calibration files defined in the sequence. To use the local +# directory use '.' If nothing is specified the default path in the +# datacenter is used. (overwrites command line options) +# ------------------------------------------------------------------------- +#MJPedestalC1.PathData: /data/MAGIC/Period016/ +#MJPedestalC2.PathData: /data/MAGIC/Period016/ +#MJCalibration.PathData: /data/MAGIC/Period016/ + +# ------------------------------------------------------------------------- +# Define here which signal extractor you are going to use for the +# whole calibration chain +# ------------------------------------------------------------------------- +MJPedestalC1.ExtractSignal: MExtractTimeAndChargeDigitalFilter +# ------------------------------------------------------------------------- +# Define here an extractor which can be used for the December 04 data. +# ------------------------------------------------------------------------- +#MJPedestalC1.ExtractSignal: MExtractTimeAndChargeSpline, MExtractTimeAndChargeDigitalFilterPeakSearch, MExtractTimeAndChargeSlidingWindow +# ------------------------------------------------------------------------- +# Define here parameters valid for all extractors above +# ------------------------------------------------------------------------- +#MJPedestalC1.ExtractSignal.HiGainFirst: 0 +#MJPedestalC1.ExtractSignal.HiGainLast: 15 +#MJPedestalC1.ExtractSignal.LoGainFirst: 3 +#MJPedestalC1.ExtractSignal.LoGainLast: 14 +#MJPedestalC1.ExtractSignal.SaturationLimit: 250 +#MJPedestalC1.ExtractSignal.OffsetLoGain: 1.51 +#MJPedestalC1.ExtractSignal.LoGainSwitch: 150 +#MJPedestalC1.ExtractSignal.LoGainStartShift: -2.8 +# ------------------------------------------------------------------------- +# Define here parameters valid for the spline: +# ------------------------------------------------------------------------- +#MJPedestalC1.ExtractSignal.Resolution: 0.05 +#MJPedestalC1.ExtractSignal.RiseTimeHiGain: 0.5 +#MJPedestalC1.ExtractSignal.FallTimeHiGain: 1.5 +#MJPedestalC1.ExtractSignal.LoGainStretch: 1.5 +#MJPedestalC1.ExtractSignal.ExtractionType: Integral +#MJPedestalC1.ExtractSignal.ExtractionType: Amplitude +# ------------------------------------------------------------------------- +# Define here parameters valid for the sliding window: +# ------------------------------------------------------------------------- +#MJPedestalC1.ExtractSignal.HiGainWindowSize: 6 +#MJPedestalC1.ExtractSignal.LoGainWindowSize: 6 +# ------------------------------------------------------------------------- +# Define here parameters valid for the digital filter with peak search: +# ------------------------------------------------------------------------- +#MJPedestalC1.ExtractSignal.OffsetLeftFromPeak: 3 +#MJPedestalC1.ExtractSignal.OffsetRightFromPeak: 3 +#MJPedestalC1.ExtractSignal.PeakSearchWindowSize: 2 +#MJPedestalC1.ExtractSignal.HiGainFailureLimit: 10 +#MJPedestalC1.ExtractSignal.LoGainFailureLimit: 25 + +# ------------------------------------------------------------------------- +# In the case of the Digital filter you have to overwrite the default +# weights file depending on what you are extracting +# ------------------------------------------------------------------------- +MJPedestalC2.ExtractSignal.WeightsFile: msignal/calibration_weights_UV46.dat +MJPedestalY2.ExtractSignal.WeightsFile: msignal/cosmics_weights46.dat +MJPedestalY3.ExtractSignal.WeightsFile: msignal/cosmics_weights46.dat +MJCalibration.ExtractSignal.WeightsFile: msignal/calibration_weights_UV46.dat +MJCalibrateSignal.ExtractSignal.WeightsFile: msignal/cosmics_weights46.dat +MJCalibrateSignal.ExtractInterlaced.WeightsFile: msignal/calibration_weights_UV46.dat + +#MJPedestalC2.ExtractSignal.AutomaticWeightsFile: On +#MJPedestalY2.ExtractSignal.AutomaticWeightsFile: On +#MJPedestalY3.ExtractSignal.AutomaticWeightsFile: On +#MJCalibration.ExtractSignal.AutomaticWeightsFile: On +#MJCalibrateSignal.ExtractSignal.AutomaticWeightsFile: On +#MJCalibrateSignal.ExtractInterlaced.AutomaticWeightsFile: On + +# ------------------------------------------------------------------------- +# Configure MJCalibration +# ------------------------------------------------------------------------- + +# Switch on relative time calibration +#MJCalibration.RelTimeCalibration: Yes +# Type of used data format: raw, root, mc +#MJCalibration.DataType: Root +# Type of displayed plots: Full,DataCheck,Normal +#MJCalibration.Display: datacheck +# Write additinal debug output +#MJCalibration.Debug: No +# Use blind pixel +#MJCalibration.UseBlindPixel: No +# Use pin diode +#MJCalibration.UsePINDiode: No +# Use pixel check option +#MJCalibration.PixelCheck: Yes +#MJCalibTest.PixelCheck: Yes + +# ------------------------------------------------------------------------- +# Individual events check for calibration results +# ------------------------------------------------------------------------- + +#ContCosmics.MaxEmptyPixels: 0.2 +#ContCosmics.MaxExcludedFraction: 1 +#ContCosmics.MinAcceptedFraction: 0 + +# ------------------------------------------------------------------------- +# You can choose the pedestal extraction algorithm/task. To use +# MPedCalcPedRun use the lines below. Be carefull, a pedestal file +# could be (if no pedestal file available) a data-file, too. +# ------------------------------------------------------------------------- +#MJPedestalC1.ExtractPedestal: MPedCalcPedRun +#MJPedestalC1.ExtractPedestal.ExtractWinFirst: 0 +#MJPedestalC1.ExtractPedestal.ExtractWinSize: 6 +#MJPedestalC1.ExtractPedestal.PedestalUpdate: no +#MJPedestalC1.ExtractPedestal.CheckWinFirst: 0 +#MJPedestalC1.ExtractPedestal.CheckWinSize: 29 + +#MJPedestalC2.ExtractPedestal: MPedCalcPedRun +#MJPedestalC2.ExtractPedestal.ExtractWinFirst: 0 +#MJPedestalC2.ExtractPedestal.ExtractWinSize: 6 +#MJPedestalC2.ExtractPedestal.PedestalUpdate: no +#MJPedestalC2.ExtractPedestal.CheckWinFirst: 0 +#MJPedestalC2.ExtractPedestal.CheckWinSize: 29 + +# ------------------------------------------------------------------------- +# You can choose the pedestal extraction algorithm/task. To use +# MPedCalcFromLoGain use the lines below. Be carefull, a pedestal file +# could be (if no pedestal file available) a data-file, too. +# ------------------------------------------------------------------------- +#MJPedestalC1.ExtractPedestal: MPedCalcFromLoGain +#MJPedestalC1.ExtractPedestal.ExtractWinFirst: 0 +#MJPedestalC1.ExtractPedestal.ExtractWinSize: 6 +#MJPedestalC1.ExtractPedestal.PedestalUpdate: no +#MJPedestalC1.ExtractPedestal.CheckWinFirst: 0 +#MJPedestalC1.ExtractPedestal.CheckWinSize: 29 +#MJPedestalC1.ExtractPedestal.MaxSignalVar: 40 + +#MJPedestalC2.ExtractPedestal: MPedCalcFromLoGain +#MJPedestalC2.ExtractPedestal.ExtractWinFirst: 0 +#MJPedestalC2.ExtractPedestal.ExtractWinSize: 6 +#MJPedestalC2.ExtractPedestal.PedestalUpdate: no +#MJPedestalC2.ExtractPedestal.CheckWinFirst: 0 +#MJPedestalC2.ExtractPedestal.CheckWinSize: 29 +#MJPedestalC2.ExtractPedestal.MaxSignalVar: 40 + +# ------------------------------------------------------------------------- +# Configure MJPedestalC1 and MJPedestalC2 +# ------------------------------------------------------------------------- +# Possible data type: Raw,Root,MC +#MJPedestalC1.DataType: Root +#MJPedestalC2.DataType: Root +#MJPedestalC3.DataType: Root +# Type of displayed plots: Full,DataCheck,None +#MJPedestalC1.Display: datacheck +#MJPedestalC2.Display: datacheck +#MJPedestalC3.Display: datacheck +# Maximum number of event processed in the loop +#MJPedestalC1.MaxEvents: 1000 +#MJPedestalC2.MaxEvents: 1000 +# Allow to overwrite existing files with the output file +#MJPedestalC1.AllowOverwrite: No +#MJPedestalC2.AllowOverwrite: No +# Use data runs from the sequence instead of calibration runs +#MJPedestalC1.UseData: No +#MJPedestalC2.UseData: No +# Switch off DeadPixelCheck with "No" +#MJPedestalY2.DeadPixelCheck: yes + + +# ------------------------------------------------------------------------- +# Use this if you want to change the software low-gain switch in the +# calibration +# ------------------------------------------------------------------------- +MJCalibration.ExtractSignal.LoGainSwitch: 170 + +# ------------------------------------------------------------------------- +# Use this if you want to change the higain-vs. logain intercalibration +# ------------------------------------------------------------------------- +#MJCalibration.HiLoCalibFile: mjobs/hilocalib_df46.root +#MJCalibration.HiLoCalibration: yes + +# ------------------------------------------------------------------------- +# Use this if you want to change the time extractor for the calibration +# and automatically the data extraction +# ------------------------------------------------------------------------- +#MJCalibration.ExtractTime: +#MJCalibration.ExtractTime: MExtractTimeFastSpline +#MJCalibration.ExtractTime.HiGainFirst: 0 +#MJCalibration.ExtractTime.HiGainLast: 7 +#MJCalibration.ExtractTime.LoGainFirst: 3 +#MJCalibration.ExtractTime.LoGainLast: 8 +#MJCalibration.ExtractTime.WindowSizeHiGain: 6 +#MJCalibration.ExtractTime.WindowSizeLoGain: 6 +#MJCalibration.ExtractTime.OffsetLoGain: 1.51 + +# ------------------------------------------------------------------------- +# Use this to change the behaviour of the calibration +# ------------------------------------------------------------------------- +# Type if you set a colour explicitely from outside (TEST purposes!!) +#MJCalibration.MCalibColorSet.ExplicitColor: green,blue,uv,ct1 + +#MJCalibration.MCalibrationChargeCalc.ArrTimeRmsLimit: 5.0 +#MJCalibration.MCalibrationChargeCalc.ChargeLimit: 2.5 +#MJCalibration.MCalibrationChargeCalc.ChargeErrLimit: 0 +#MJCalibration.MCalibrationChargeCalc.ChargeRelErrLimit: 1 +#MJCalibration.MCalibrationChargeCalc.Debug: no + +#MJCalibration.MCalibrationChargeCalc.FFactorErrLimit: 4.5 +#MJCalibration.MCalibrationChargeCalc.LambdaErrLimit: 0.2 +#MJCalibration.MCalibrationChargeCalc.LambdaCheckLimit: 0.5 +#MJCalibration.MCalibrationChargeCalc.PheErrLowerLimit: 9.0 +#MJCalibration.MCalibrationChargeCalc.PheErrUpperLimit: 5.5 + +#MJCalibration.MHCalibrationChargeCam.Debug: no +#MJCalibration.MHCalibrationChargeCam.LoGain: yes +#MJCalibration.MHCalibrationChargeCam.Oscillations: yes +#MJCalibration.MHCalibrationChargeCam.SizeCheck: yes +#MJCalibration.MHCalibrationChargeCam.Averageing: yes +#MJCalibration.MHCalibrationChargeCam.HiGainNbins: 500 +#MJCalibration.MHCalibrationChargeCam.HiGainFirst: -98. +#MJCalibration.MHCalibrationChargeCam.HiGainLast: 1902. +#MJCalibration.MHCalibrationChargeCam.LoGainNbins: 500 +#MJCalibration.MHCalibrationChargeCam.LoGainFirst: -99. +#MJCalibration.MHCalibrationChargeCam.LoGainLast: 901. +#MJCalibration.MHCalibrationChargeCam.TimeLowerLimit: 1. +#MJCalibration.MHCalibrationChargeCam.TimeUpperLimit: 3. +#MJCalibration.MHCalibrationChargeCam.NumHiGainSaturationLimit: 0.02 +#MJCalibration.MHCalibrationChargeCam.NumLoGainSaturationLimit: 0.005 +#MJCalibration.MHCalibrationChargeCam.ProbLimit: 0.00000001 +#MJCalibration.MHCalibrationChargeCam.OverflowLimit: 0.005 +#MJCalibration.MHCalibrationChargeCam.PulserFrequency: 500 + + +#MJCalibration.MHCalibrationRelTimeCam.Debug: no +#MJCalibration.MHCalibrationRelTimeCam.LoGain: no +#MJCalibration.MHCalibrationRelTimeCam.Oscillations: yes +#MJCalibration.MHCalibrationRelTimeCam.SizeCheck: yes +#MJCalibration.MHCalibrationRelTimeCam.Averageing: yes +#MJCalibration.MHCalibrationRelTimeCam.HiGainNbins: 210 +#MJCalibration.MHCalibrationRelTimeCam.HiGainFirst: -5. +#MJCalibration.MHCalibrationRelTimeCam.HiGainLast: 10. +#MJCalibration.MHCalibrationRelTimeCam.NumHiGainSaturationLimit: 0.25 +#MJCalibration.MHCalibrationRelTimeCam.ProbLimit: 0.00001 +#MJCalibration.MHCalibrationRelTimeCam.OverflowLimit: 0.005 +#MJCalibration.MHCalibrationRelTimeCam.PulserFrequency: 500 + + +# ========================================================================== +############################################################################# +# ========================================================================== +# Signal extraction +# ========================================================================== +############################################################################# +# ========================================================================== + +# ------------------------------------------------------------------------- +# Use a OutputPath if you want to write the results to a file +# ------------------------------------------------------------------------- +#MJCalibrateSignal.PathOut: . +#MJPedestalY1.PathOut: . +#MJPedestalY2.PathOut: . +#MJPedestalY3.PathOut: . +#MJPedestalY1.PathData: /data/MAGIC/Period016/ +#MJPedestalY2.PathData: /data/MAGIC/Period016/ +#MJPedestalY3.PathData: /data/MAGIC/Period016/ +#MJCalibrateSignal.PathData: /data/MAGIC/Period016/ +#MJCalibrateSignal.PathIn: . + +# ------------------------------------------------------------------------- +# If you want to change the behaviour of the pedestal extraction... +# ------------------------------------------------------------------------- +# Possible data type: Raw,Root,MC +#MJPedestalY1.DataType: Root +#MJPedestalY2.DataType: Root +#MJPedestalY3.DataType: Root +#MJPedestalY1.UseData: Yes +#MJPedestalY2.UseData: Yes +#MJPedestalY3.UseData: Yes + +MJPedestalY1.MaxEvents: 1000 +MJPedestalY2.MaxEvents: 10000 +MJPedestalY3.MaxEvents: 1000 + +MJPedestalY1.MaxPedestals: 500 +MJPedestalY2.MaxPedestals: 500 +MJPedestalY3.MaxPedestals: 500 + +# ------------------------------------------------------------------------- +# Use Pulse Position check to define the extraction ranges for the data? +# ------------------------------------------------------------------------- +#MJPedestalY2.PulsePosCheck: yes +#MJPedestalY2.ExtractWinLeft: 0.0 +#MJPedestalY2.ExtractWinRight: 0.5 + +# ------------------------------------------------------------------------- +# Define the Pulse Position check parameters: +# ------------------------------------------------------------------------- +#MJPedestalY2.MHCalibrationPulseTimeCam.SaturationLimit: 255 +#MJPedestalY2.MHCalibrationPulseTimeCam.LowerSignalLimit: 100 +#MJPedestalY2.MHCalibrationPulseTimeCam.NumPixelsRequired: 2 +#MJPedestalY2.PixelCheck: no + +#MJPedestalY.ExtractPedestal: MPedCalcFromLoGain +#MJPedestalY.ExtractPedestal.PedestalUpdate: no + +# ------------------------------------------------------------------------- +# Setup pedestal extractor for running through data +# ------------------------------------------------------------------------- +#MJCalibrateSignal.MPedCalcFundamental.PedestalUpdate: yes +#MJCalibrateSignal.MPedCalcFundamental.NumDump: 500 +#MJCalibrateSignal.MPedCalcFundamental.NumEventsDump: 500 +#MJCalibrateSignal.MPedCalcFundamental.NumAreasDump: 500 +#MJCalibrateSignal.MPedCalcFundamental.NumSectorsDump: 500 +#MJCalibrateSignal.MPedCalcFundamental.ExtractWinFirst: 17 +#MJCalibrateSignal.MPedCalcFundamental.ExtractWinSize: 6 +#MJCalibrateSignal.MPedCalcFundamental.CheckWinFirst: 0 +#MJCalibrateSignal.MPedCalcFundamental.CheckWinSize: 29 +#MJCalibrateSignal.MPedCalcFundamental.MaxSignalVar: 40 + +#MJCalibrateSignal.MPedCalcFromExtractorRndm.PedestalUpdate: yes +#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumEventsDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumAreasDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumSectorsDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.ExtractWinFirst: 17 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.ExtractWinSize: 6 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.CheckWinFirst: 0 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.CheckWinSize: 29 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.MaxSignalVar: 40 + +#MJCalibrateSignal.MPedCalcFromExtractor.PedestalUpdate: yes +#MJCalibrateSignal.MPedCalcFromExtractor.NumDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractor.NumEventsDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractor.NumAreasDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractor.NumSectorsDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractor.ExtractWinFirst: 17 +#MJCalibrateSignal.MPedCalcFromExtractor.ExtractWinSize: 6 +#MJCalibrateSignal.MPedCalcFromExtractor.CheckWinFirst: 0 +#MJCalibrateSignal.MPedCalcFromExtractor.CheckWinSize: 29 +#MJCalibrateSignal.MPedCalcFromExtractor.MaxSignalVar: 40 + +# ------------------------------------------------------------------------- +# Configure the interlaced calibration updates +# ------------------------------------------------------------------------- +#MJCalibrateSignal.MCalibCalcFromPast.UpdateWithFFactorMethod: yes +#MJCalibrateSignal.MCalibCalcFromPast.NumEventsDump: 500 +#MJCalibrateSignal.MCalibCalcFromPast.UpdateNumPhes: yes +#MJCalibrateSignal.MCalibCalcFromPast.NumPhesDump: 10 +#MJCalibrateSignal.MCalibrationChargeCalc.Debug: no + +# ------------------------------------------------------------------------- +# Use this if you want to change the signal extractor +# ------------------------------------------------------------------------- +MJCalibrateSignal.ExtractSignal.LoGainSwitch: 150 +#MJCalibrateSignal.ExtractSignal.LoGainStretch: 1.5 +#MJCalibrateSignal.ExtractInterlaced.LoGainStretch: 1.5 + +# ------------------------------------------------------------------------- +# Setup calibration of data +# ------------------------------------------------------------------------- +# Type of used data format: raw,root,MC +#MJCalibrateSignal.DataType: Root +# Type if you set a colour explicitely from outside (only TESTS!!!) +#MJCalibrateSignal.MCalibColorSet.ExpicitColor: green,blue,uv,ct1 +#MJCalibrateSignal.MCalibrateData.PedestalFlag: Event +#MJCalibrateSignal.MCalibrateData.ScaleFactor: 1 +MJCalibrateSignal.MCalibrateData.FileNameScale: resources/calibration.rc +#MJCalibrateSignal.MCalibrateData.CalibrationMode: Default +#MJCalibrateSignal.MCalibrateData.CalibConvMinLimit: 0.01 +#MJCalibrateSignal.MCalibrateData.CalibConvMaxLimit: 5. +#MJCalibrateSignal.Interlaced: yes +#MJCalibrateSignal.RelTimesUpdate: yes + +# ------------------------------------------------------------------------- +# Setup level for determining of bad pixels +# ------------------------------------------------------------------------- +#MJCalibrateSignal.MBadPixelsCalc.PedestalLevel: 3.0 +#MJCalibrateSignal.MBadPixelsCalc.PedestalLevelVariance: 3.0 +#MJCalibrateSignal.MBadPixelsTreat.NumMinNeighbors: 3 +#MJCalibrateSignal.MBadPixelsTreat.MaxArrivalTimeDiff: 1.0 +#MJCalibrateSignal.MBadPixelsTreat.UseInterpolation: yes +#MJCalibrateSignal.MBadPixelsTreat.ProcessPedestalEvt: yes +#MJCalibrateSignal.MBadPixelsTreat.ProcessPedestalRun: no +#MJCalibrateSignal.MBadPixelsTreat.ProcessTimes: yes +#MJCalibrateSignal.MBadPixelsTreat.UseCentralPixel: no +#MJCalibrateSignal.MBadPixelsTreat.HardTreatment: no + +# ------------------------------------------------------------------------- +# Setup interlaced calibration events histogramming +# ------------------------------------------------------------------------- +#MJCalibrateSignal.MHCalibrationChargeCam.Debug: no +#MJCalibrateSignal.MHCalibrationChargeCam.LoGain: yes +MJCalibrateSignal.MHCalibrationChargeCam.Oscillations: no +#MJCalibrateSignal.MHCalibrationChargeCam.SizeCheck: yes +#MJCalibrateSignal.MHCalibrationChargeCam.Averageing: yes +#MJCalibrateSignal.MHCalibrationChargeCam.HiGainNbins: 500 +#MJCalibrateSignal.MHCalibrationChargeCam.HiGainFirst: -98. +#MJCalibrateSignal.MHCalibrationChargeCam.HiGainLast: 1902. +#MJCalibrateSignal.MHCalibrationChargeCam.HiGainLast: 1899.5 +MJCalibrateSignal.MHCalibrationChargeCam.LoGainNbins: 250 +MJCalibrateSignal.MHCalibrationChargeCam.LoGainFirst: -98 +MJCalibrateSignal.MHCalibrationChargeCam.LoGainLast: 902 +#MJCalibrateSignal.MHCalibrationChargeCam.TimeLowerLimit: 1. +#MJCalibrateSignal.MHCalibrationChargeCam.TimeUpperLimit: 3. +#MJCalibrateSignal.MHCalibrationChargeCam.NumHiGainSaturationLimit: 0.02 +#MJCalibrateSignal.MHCalibrationChargeCam.NumLoGainSaturationLimit: 0.005 +#MJCalibrateSignal.MHCalibrationChargeCam.ProbLimit: 0.0001 +#MJCalibrateSignal.MHCalibrationChargeCam.OverflowLimit: 0.005 +MJCalibrateSignal.MHCalibrationChargeCam.PulserFrequency: 50 + +# ------------------------------------------------------------------------- +# Setup relative time interlaced calibration +# ------------------------------------------------------------------------- +#MJCalibrateSignal.MHCalibrationRelTimeCam.Debug: no +#MJCalibrateSignal.MHCalibrationRelTimeCam.LoGain: no +MJCalibrateSignal.MHCalibrationRelTimeCam.Oscillations: no +#MJCalibrateSignal.MHCalibrationRelTimeCam.SizeCheck: yes +#MJCalibrateSignal.MHCalibrationRelTimeCam.Averageing: yes +#MJCalibrateSignal.MHCalibrationRelTimeCam.HiGainNbins: 210 +#MJCalibrateSignal.MHCalibrationRelTimeCam.HiGainFirst: -5. +#MJCalibrateSignal.MHCalibrationRelTimeCam.HiGainLast: 10. +#MJCalibrateSignal.MHCalibrationRelTimeCam.NumHiGainSaturationLimit: 0.25 +#MJCalibrateSignal.MHCalibrationRelTimeCam.ProbLimit: 0.00001 +#MJCalibrateSignal.MHCalibrationRelTimeCam.OverflowLimit: 0.005 +MJCalibrateSignal.MHCalibrationRelTimeCam.PulserFrequency: 50 + +# ------------------------------------------------------------------------- +# Resource for movie encoding. For more details see IACT wiki +# ------------------------------------------------------------------------- +#MMoviePrepare.FirstSlice: 0 +#MMoviePrepare.LastSlice: 15 +#MovieFilter.Condition: MMovieData.fMax>5*MMovieData.fMedianPedestalRms +#MImgCleanStd.CleanLevel1: 8.0 +#MImgCleanStd.CleanLevel2: 4.5 +#MFEvtNumber.FileName: ganymed00001218.root +#MFEvtNumber.Selector: ThetaSquared.fVal<0.04 && DataType.fVal>0.5 +#MMovieWrite.TargetLength: 5 +#MMovieWrite.NumEvents: 500 +#MMovieWrite.Threshold: 3 +#MMovieWrite.FileName: filename.mpg + +# ------------------------------------------------------------------------- +# This is a special option for the datacheck. A task executed as last +# task in the event processing task list - in the standard software it +# is skipped +# ------------------------------------------------------------------------- +#MJCalibrateSignal.FinalFantasy: MHCamFlorian Index: trunk/Mars/resources/callisto_mux.rc =================================================================== --- trunk/Mars/resources/callisto_mux.rc (revision 19872) +++ trunk/Mars/resources/callisto_mux.rc (revision 19872) @@ -0,0 +1,479 @@ +# ========================================================================== +############################################################################# +# ========================================================================== +# General +# ========================================================================== +############################################################################# +# ========================================================================== + +# ------------------------------------------------------------------------- +# Use this if you want to setup the logging stream for the jobs +# (overwrites command line options) +# ------------------------------------------------------------------------- +#MLog.VerbosityLevel: 2 +#MLog.DebugLevel: 1 +#MLog.NoColors: yes + +# ========================================================================== +############################################################################# +# ========================================================================== +# Calibration +# ========================================================================== +############################################################################# +# ========================================================================== + +# ------------------------------------------------------------------------- +# Use this if you want to write the MJPedestal output somewhere +# If you don't want it, it is written to the calibration output anyhow. +# ------------------------------------------------------------------------- +#MJPedestalC1.PathOut: . +#MJPedestalC2.PathOut: . + +# ------------------------------------------------------------------------- +# Use this if you want to modify the exclusion of pixels from the beginning +# (default: Pixel 0 and 559) +# ------------------------------------------------------------------------- +MJPedestalC1.BadPixelsFile: mjobs/badpixels_0_559.rc + +# ------------------------------------------------------------------------- +# Use this if you want to modify files containing information about the +# reference lines in the status display +# ------------------------------------------------------------------------- +MJPedestalC1.ReferenceFile: resources/pedestalref_mux.rc +MJPedestalC2.ReferenceFile: resources/pedestalref_mux.rc +MJCalibration.ReferenceFile: resources/calibrationref_mux.rc +MJCalibration.MHCalibrationRelTimeCam.ReferenceFile: resources/calibrationref_mux.rc +MJCalibration.MHCalibrationChargeCam.ReferenceFile: resources/calibrationref_mux.rc + +# ------------------------------------------------------------------------- +# Use this to define where the calibration output is stored. The filename +# is created from the sequence number. If nothing is specified '.' is +# assumed. (overwrites command line options) +# ------------------------------------------------------------------------- +#MJCalibration.PathOut: calped + +# ------------------------------------------------------------------------- +# Use this to define where the program should search for the pedestal +# and calibration files defined in the sequence. To use the local +# directory use '.' If nothing is specified the default path in the +# datacenter is used. (overwrites command line options) +# ------------------------------------------------------------------------- +#MJPedestalC1.PathData: /data/MAGIC/Period016/ +#MJPedestalC2.PathData: /data/MAGIC/Period016/ +#MJCalibration.PathData: /data/MAGIC/Period016/ + +# ------------------------------------------------------------------------- +# Define here which signal extractor you are going to use for the +# whole calibration chain +# ------------------------------------------------------------------------- +MJPedestalC1.ExtractSignal: MExtractTimeAndChargeSpline +# ------------------------------------------------------------------------- +# Possibilities for a different extractor +# ------------------------------------------------------------------------- +#MJPedestalC1.ExtractSignal: MExtractTimeAndChargeDigitalFilter, MExtractTimeAndChargeDigitalFilterPeakSearch, MExtractTimeAndChargeSlidingWindow +# ------------------------------------------------------------------------- +# Define here parameters valid for all extractors above +# The actual defaults used depend on the extractor in use. +# For details see the corresponding source code +# ------------------------------------------------------------------------- +MJPedestalC1.ExtractSignal.HiGainFirst: 0 +MJPedestalC1.ExtractSignal.HiGainLast: 49 +#MJPedestalC1.ExtractSignal.LoGainFirst: 0 +#MJPedestalC1.ExtractSignal.LoGainLast: 0 +MJPedestalC1.ExtractSignal.SaturationLimit: 252 +#MJPedestalC1.ExtractSignal.OffsetLoGain: 1.51 +#MJPedestalC1.ExtractSignal.LoGainSwitch: 150 +#MJPedestalC1.ExtractSignal.LoGainStartShift: -2.8 +# ------------------------------------------------------------------------- +# Define here parameters valid for the spline: +# ------------------------------------------------------------------------- +#MJPedestalC1.ExtractSignal.Resolution: 0.05 +MJPedestalC1.ExtractSignal.RiseTimeHiGain: 4.0 +MJPedestalC1.ExtractSignal.FallTimeHiGain: 5.0 +#MJPedestalC1.ExtractSignal.LoGainStretch: 1.5 +#MJPedestalC1.ExtractSignal.ExtractionType: Integral +#MJPedestalC1.ExtractSignal.ExtractionType: Amplitude + +# ------------------------------------------------------------------------- +# In the case of the Digital filter you have to overwrite the default +# weights file depending on what you are extracting +# ------------------------------------------------------------------------- +#MJPedestalC2.ExtractSignal.WeightsFile: msignal/calibration_weights_UV4.dat +#MJPedestalY2.ExtractSignal.WeightsFile: msignal/cosmics_weights4.dat +#MJPedestalY3.ExtractSignal.WeightsFile: msignal/cosmics_weights4.dat +#MJCalibration.ExtractSignal.WeightsFile: msignal/calibration_weights_UV4.dat +#MJCalibrateSignal.ExtractSignal.WeightsFile: msignal/cosmics_weights4.dat +#MJCalibrateSignal.ExtractInterlaced.WeightsFile: msignal/calibration_weights_UV4.dat +#MJCalibrateSignal.ExtractSignal.LoGainStretch: 1.5 +#MJCalibrateSignal.ExtractInterlaced.LoGainStretch: 1.5 + +# ------------------------------------------------------------------------- +# Configure MJCalibration +# ------------------------------------------------------------------------- + +# Switch on relative time calibration +#MJCalibration.RelTimeCalibration: Yes +# Set color to be used +#MJCalibration.Color: +# Type of used data format: raw, root, mc +#MJCalibration.DataType: Root +# Type of displayed plots: Full,DataCheck,Normal +#MJCalibration.Display: datacheck +# Write additinal debug output +#MJCalibration.Debug: No +# Use blind pixel +#MJCalibration.UseBlindPixel: No +# Use pin diode +#MJCalibration.UsePINDiode: No +# Use pixel check option +#MJCalibration.PixelCheck: No +#MJCalibTest.PixelCheck: No + +# ------------------------------------------------------------------------- +# Individual events check for calibration results +# ------------------------------------------------------------------------- + +#ContCosmics.MaxEmptyPixels: 0.2 +#ContCosmics.MaxExcludedFraction: 1 +#ContCosmics.MinAcceptedFraction: 0 + +# ------------------------------------------------------------------------- +# In case of MOON DATA, uncomment the following lines: +# (better: call callisto with option "-moon") +# ------------------------------------------------------------------------- +#MJPedestalC1.UseHists: yes +#MJPedestalC2.UseHists: yes +#MJPedestalY1.UseHists: yes +#MJPedestalY2.UseHists: yes +#MJPedestalY3.UseHists: yes + +# ------------------------------------------------------------------------- +# You can choose the pedestal extraction algorithm/task. To use +# MPedCalcPedRun use the lines below. Be carefull, a pedestal file +# could be (if no pedestal file available) a data-file, too. +# ------------------------------------------------------------------------- +#MJPedestalC1.ExtractPedestal: MPedCalcPedRun +#MJPedestalC1.ExtractPedestal.ExtractWinFirst: 0 +#MJPedestalC1.ExtractPedestal.ExtractWinSize: 6 +#MJPedestalC1.ExtractPedestal.PedestalUpdate: no +#MJPedestalC1.ExtractPedestal.CheckWinFirst: 0 +#MJPedestalC1.ExtractPedestal.CheckWinSize: 29 + +#MJPedestalC2.ExtractPedestal: MPedCalcPedRun +#MJPedestalC2.ExtractPedestal.ExtractWinFirst: 0 +#MJPedestalC2.ExtractPedestal.ExtractWinSize: 6 +#MJPedestalC2.ExtractPedestal.PedestalUpdate: no +#MJPedestalC2.ExtractPedestal.CheckWinFirst: 0 +#MJPedestalC2.ExtractPedestal.CheckWinSize: 29 + +# ------------------------------------------------------------------------- +# You can choose the pedestal extraction algorithm/task. To use +# MPedCalcFromLoGain use the lines below. Be carefull, a pedestal file +# could be (if no pedestal file available) a data-file, too. +# ------------------------------------------------------------------------- +#MJPedestalC1.ExtractPedestal: MPedCalcFromLoGain +#MJPedestalC1.ExtractPedestal.ExtractWinFirst: 0 +#MJPedestalC1.ExtractPedestal.ExtractWinSize: 6 +#MJPedestalC1.ExtractPedestal.PedestalUpdate: no +#MJPedestalC1.ExtractPedestal.CheckWinFirst: 0 +#MJPedestalC1.ExtractPedestal.CheckWinSize: 29 +#MJPedestalC1.ExtractPedestal.MaxSignalVar: 40 + +#MJPedestalC2.ExtractPedestal: MPedCalcFromLoGain +#MJPedestalC2.ExtractPedestal.ExtractWinFirst: 0 +#MJPedestalC2.ExtractPedestal.ExtractWinSize: 6 +#MJPedestalC2.ExtractPedestal.PedestalUpdate: no +#MJPedestalC2.ExtractPedestal.CheckWinFirst: 0 +#MJPedestalC2.ExtractPedestal.CheckWinSize: 29 +#MJPedestalC2.ExtractPedestal.MaxSignalVar: 40 + +# ------------------------------------------------------------------------- +# Configure MJPedestalC1 and MJPedestalC2 +# ------------------------------------------------------------------------- +# Possible data type: Raw,Root,MC +#MJPedestalC1.DataType: Root +#MJPedestalC2.DataType: Root +#MJPedestalC3.DataType: Root +# Type of displayed plots: Full,DataCheck,None +#MJPedestalC1.Display: datacheck +#MJPedestalC2.Display: datacheck +#MJPedestalC3.Display: datacheck +# Maximum number of event processed in the loop +#MJPedestalC1.MaxEvents: 1000 +#MJPedestalC2.MaxEvents: 1000 +# Allow to overwrite existing files with the output file +#MJPedestalC1.AllowOverwrite: No +#MJPedestalC2.AllowOverwrite: No +# Use data runs from the sequence instead of calibration runs +#MJPedestalC1.UseData: No +#MJPedestalC2.UseData: No + +# ------------------------------------------------------------------------- +# Use this if you want to change the software low-gain switch +# ------------------------------------------------------------------------- +#MJCalibration.ExtractSignal.LoGainSwitch: 150 + +# ------------------------------------------------------------------------- +# Use this if you want to change the higain-vs. logain intercalibration +# ------------------------------------------------------------------------- +#MJCalibration.HiLoCalibFile: resources/hilocalib_sp1.rc + +# ------------------------------------------------------------------------- +# Use this to change the behaviour of the calibration +# ------------------------------------------------------------------------- +# Type if you set a colour explicitely from outside (only for MC!!!) +#MJCalibration.MCalibColorSet.ExplicitColor: green,blue,uv,ct1 + +#MJCalibration.MCalibrationChargeCalc.ArrTimeRmsLimit: 5.0 +#MJCalibration.MCalibrationChargeCalc.ChargeLimit: 2.5 +#MJCalibration.MCalibrationChargeCalc.ChargeErrLimit: 0 +#MJCalibration.MCalibrationChargeCalc.ChargeRelErrLimit: 1 +#MJCalibration.MCalibrationChargeCalc.Debug: no + +#MJCalibration.MCalibrationChargeCalc.FFactorErrLimit: 4.5 +#MJCalibration.MCalibrationChargeCalc.LambdaErrLimit: 0.2 +#MJCalibration.MCalibrationChargeCalc.LambdaCheckLimit: 0.5 +#MJCalibration.MCalibrationChargeCalc.PheErrLowerLimit: 9.0 +#MJCalibration.MCalibrationChargeCalc.PheErrUpperLimit: 5.5 + +#MJCalibration.MHCalibrationChargeCam.Debug: no +#MJCalibration.MHCalibrationChargeCam.LoGain: yes +#MJCalibration.MHCalibrationChargeCam.Oscillations: yes +#MJCalibration.MHCalibrationChargeCam.SizeCheck: yes +#MJCalibration.MHCalibrationChargeCam.Averageing: yes +MJCalibration.MHCalibrationChargeCam.HiGainNbins: 800 +MJCalibration.MHCalibrationChargeCam.HiGainFirst: -0.5 +MJCalibration.MHCalibrationChargeCam.HiGainLast: 3999.5 +#MJCalibration.MHCalibrationChargeCam.LoGainNbins: 500 +#MJCalibration.MHCalibrationChargeCam.LoGainFirst: -100.25 +#MJCalibration.MHCalibrationChargeCam.LoGainLast: 899.75 +#MJCalibration.MHCalibrationChargeCam.TimeLowerLimit: 1. +#MJCalibration.MHCalibrationChargeCam.TimeUpperLimit: 3. +#MJCalibration.MHCalibrationChargeCam.NumHiGainSaturationLimit: 0.02 +#MJCalibration.MHCalibrationChargeCam.NumLoGainSaturationLimit: 0.005 +MJCalibration.MHCalibrationChargeCam.ProbLimit: 1.0E-18 +#MJCalibration.MHCalibrationChargeCam.OverflowLimit: 0.005 +#MJCalibration.MHCalibrationChargeCam.PulserFrequency: 500 + +#MJCalibration.MHCalibrationRelTimeCam.Debug: no +#MJCalibration.MHCalibrationRelTimeCam.LoGain: no +#MJCalibration.MHCalibrationRelTimeCam.Oscillations: yes +#MJCalibration.MHCalibrationRelTimeCam.SizeCheck: yes +#MJCalibration.MHCalibrationRelTimeCam.Averageing: yes +MJCalibration.MHCalibrationRelTimeCam.Nbins: 751 +MJCalibration.MHCalibrationRelTimeCam.First: -60 +MJCalibration.MHCalibrationRelTimeCam.Last: 60 +#MJCalibration.MHCalibrationRelTimeCam.NumHiGainSaturationLimit: 0.25 +#MJCalibration.MHCalibrationRelTimeCam.ProbLimit: 0.00001 +#MJCalibration.MHCalibrationRelTimeCam.OverflowLimit: 0.005 +#MJCalibration.MHCalibrationRelTimeCam.PulserFrequency: 500 + +# ========================================================================== +############################################################################# +# ========================================================================== +# Signal extraction +# ========================================================================== +############################################################################# +# ========================================================================== + +# ------------------------------------------------------------------------- +# Use a OutputPath if you want to write the results to a file +# ------------------------------------------------------------------------- +#MJCalibrateSignal.PathOut: . +#MJPedestalY1.PathOut: . +#MJPedestalY2.PathOut: . +#MJPedestalY3.PathOut: . +#MJPedestalY1.PathData: /data/MAGIC/Period016/ +#MJPedestalY2.PathData: /data/MAGIC/Period016/ +#MJPedestalY3.PathData: /data/MAGIC/Period016/ +#MJCalibrateSignal.PathData: /data/MAGIC/Period016/ +#MJCalibrateSignal.PathIn: . + +# ------------------------------------------------------------------------- +# If you want to change the behaviour of the pedestal extraction... +# ------------------------------------------------------------------------- +# Possible data type: Raw,Root,MC +#MJPedestalY1.DataType: Root +#MJPedestalY2.DataType: Root +#MJPedestalY3.DataType: Root +#MJPedestalY1.UseData: Yes +#MJPedestalY2.UseData: Yes +#MJPedestalY3.UseData: Yes +#MJPedestalY1.MinEvents: 50 +#MJPedestalY2.MinEvents: 50 +#MJPedestalY3.MinEvents: 50 + +MJPedestalY1.MaxEvents: 50000 +MJPedestalY2.MaxEvents: 100000 +MJPedestalY3.MaxEvents: 50000 + +MJPedestalY1.MaxPedestals: 500 +MJPedestalY2.MaxPedestals: 500 +MJPedestalY3.MaxPedestals: 500 + +# ------------------------------------------------------------------------- +# When the starting pedestal is extracted from the first events of +# the data run also the average pulse position of these events +# is determined. If the pulse position check is switched on (default) +# the extraction window of the hi- and lo-gain extractor for the +# determined average pulse position is not allowed to exceed the +# extractors search region. In this scenario we assume that the +# start of the extraction window is at least similar to the +# rising edge of the pulse. +# To make sure that the search region and the extraction window +# does not touch increase the values. +# ------------------------------------------------------------------------- +#MJPedestalY2.PulsePosCheck: yes +#MJPedestalY2.ExtractWinLeft: 0.0 +#MJPedestalY2.ExtractWinRight: 0.0 + +# ------------------------------------------------------------------------- +# Define the Pulse Position check parameters: +# ------------------------------------------------------------------------- +MJPedestalY2.MHCalibrationPulseTimeCam.SaturationLimit: 200 +MJPedestalY2.MHCalibrationPulseTimeCam.LowerSignalLimit: 35 +#MJPedestalY2.MHCalibrationPulseTimeCam.NumPixelsRequired: 3 +#MJPedestalY2.PixelCheck: no + +#MJPedestalY.ExtractPedestal: MPedCalcFromLoGain +#MJPedestalY.ExtractPedestal.PedestalUpdate: no + +# ------------------------------------------------------------------------- +# Setup pedestal extractor for running through data +# ------------------------------------------------------------------------- +#MJCalibrateSignal.MPedCalcFundamental.PedestalUpdate: yes +#MJCalibrateSignal.MPedCalcFundamental.NumDump: 500 +#MJCalibrateSignal.MPedCalcFundamental.NumEventsDump: 500 +#MJCalibrateSignal.MPedCalcFundamental.NumAreasDump: 500 +#MJCalibrateSignal.MPedCalcFundamental.NumSectorsDump: 500 +#MJCalibrateSignal.MPedCalcFundamental.ExtractWinFirst: 17 +#MJCalibrateSignal.MPedCalcFundamental.ExtractWinSize: 6 +#MJCalibrateSignal.MPedCalcFundamental.CheckWinFirst: 0 +#MJCalibrateSignal.MPedCalcFundamental.CheckWinSize: 29 +#MJCalibrateSignal.MPedCalcFundamental.MaxSignalVar: 40 + +#MJCalibrateSignal.MPedCalcFromExtractorRndm.PedestalUpdate: yes +#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumEventsDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumAreasDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.NumSectorsDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.ExtractWinFirst: 17 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.ExtractWinSize: 6 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.CheckWinFirst: 0 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.CheckWinSize: 29 +#MJCalibrateSignal.MPedCalcFromExtractorRndm.MaxSignalVar: 40 + +#MJCalibrateSignal.MPedCalcFromExtractor.PedestalUpdate: yes +#MJCalibrateSignal.MPedCalcFromExtractor.NumDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractor.NumEventsDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractor.NumAreasDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractor.NumSectorsDump: 500 +#MJCalibrateSignal.MPedCalcFromExtractor.ExtractWinFirst: 17 +#MJCalibrateSignal.MPedCalcFromExtractor.ExtractWinSize: 6 +#MJCalibrateSignal.MPedCalcFromExtractor.CheckWinFirst: 0 +#MJCalibrateSignal.MPedCalcFromExtractor.CheckWinSize: 29 +#MJCalibrateSignal.MPedCalcFromExtractor.MaxSignalVar: 40 + +# ------------------------------------------------------------------------- +# Configure the interlaced calibration updates +# ------------------------------------------------------------------------- +#MJCalibrateSignal.MCalibCalcFromPast.UpdateWithFFactorMethod: yes +#MJCalibrateSignal.MCalibCalcFromPast.NumEventsDump: 500 +#MJCalibrateSignal.MCalibCalcFromPast.UpdateNumPhes: yes +#MJCalibrateSignal.MCalibCalcFromPast.NumPhesDump: 10 +#MJCalibrateSignal.MCalibrationChargeCalc.Debug: no + +# ------------------------------------------------------------------------- +# Use this if you want to change the signal extractor +# ------------------------------------------------------------------------- +#MJCalibrateSignal.ExtractSignal.LoGainSwitch: 150 + +# ------------------------------------------------------------------------- +# Setup calibration of data +# ------------------------------------------------------------------------- +# Type of used data format: raw,root,MC +#MJCalibrateSignal.DataType: Root +# Type if you set a colour explicitely from outside (only for MC!!!) +#MJCalibrateSignal.MCalibColorSet.ExpicitColor: green,blue,uv,ct1 +#MJCalibrateSignal.MCalibrateData.PedestalFlag: Event +#MJCalibrateSignal.MCalibrateData.ScaleFactor: 0.85182358 +MJCalibrateSignal.MCalibrateData.FileNameScale: resources/calibration_spline.rc +#MJCalibrateSignal.MCalibrateData.CalibrationMode: Default +#MJCalibrateSignal.MCalibrateData.CalibConvMinLimit: 0.01 +#MJCalibrateSignal.MCalibrateData.CalibConvMaxLimit: 5. +#MJCalibrateSignal.Interlaced: yes +#MJCalibrateSignal.RelTimesUpdate: yes + +# ------------------------------------------------------------------------- +# Setup level for determining of bad pixels +# ------------------------------------------------------------------------- +#MJCalibrateSignal.MBadPixelsCalc.PedestalLevel: 3.0 +#MJCalibrateSignal.MBadPixelsCalc.PedestalLevelVariance: 3.0 +#MJCalibrateSignal.MBadPixelsTreat.NumMinNeighbors: 3 +#MJCalibrateSignal.MBadPixelsTreat.MaxArrivalTimeDiff: 3.0 +#MJCalibrateSignal.MBadPixelsTreat.UseInterpolation: yes +#MJCalibrateSignal.MBadPixelsTreat.ProcessPedestalEvt: yes +#MJCalibrateSignal.MBadPixelsTreat.ProcessPedestalRun: no +#MJCalibrateSignal.MBadPixelsTreat.ProcessTimes: yes +#MJCalibrateSignal.MBadPixelsTreat.UseCentralPixel: no +#MJCalibrateSignal.MBadPixelsTreat.HardTreatment: no + +# ------------------------------------------------------------------------- +# Setup interlaced calibration events histogramming +# ------------------------------------------------------------------------- +#MJCalibrateSignal.MHCalibrationChargeCam.Debug: no +#MJCalibrateSignal.MHCalibrationChargeCam.LoGain: yes +MJCalibrateSignal.MHCalibrationChargeCam.Oscillations: no +#MJCalibrateSignal.MHCalibrationChargeCam.SizeCheck: yes +#MJCalibrateSignal.MHCalibrationChargeCam.Averageing: yes +#MJCalibrateSignal.MHCalibrationChargeCam.HiGainNbins: 500 +#MJCalibrateSignal.MHCalibrationChargeCam.HiGainFirst: -100.5 +#MJCalibrateSignal.MHCalibrationChargeCam.HiGainLast: 1899.5 +MJCalibrateSignal.MHCalibrationChargeCam.LoGainNbins: 800 +MJCalibrateSignal.MHCalibrationChargeCam.LoGainFirst: -0.5 +MJCalibrateSignal.MHCalibrationChargeCam.LoGainLast: 3999.5 +#MJCalibrateSignal.MHCalibrationChargeCam.TimeLowerLimit: 1. +#MJCalibrateSignal.MHCalibrationChargeCam.TimeUpperLimit: 3. +#MJCalibrateSignal.MHCalibrationChargeCam.NumHiGainSaturationLimit: 0.02 +#MJCalibrateSignal.MHCalibrationChargeCam.NumLoGainSaturationLimit: 0.005 +#MJCalibrateSignal.MHCalibrationChargeCam.ProbLimit: 0.0001 +#MJCalibrateSignal.MHCalibrationChargeCam.OverflowLimit: 0.005 +MJCalibrateSignal.MHCalibrationChargeCam.PulserFrequency: 25 + +# ------------------------------------------------------------------------- +# Setup relative time interlaced calibration +# ------------------------------------------------------------------------- +#MJCalibrateSignal.MHCalibrationRelTimeCam.Debug: no +#MJCalibrateSignal.MHCalibrationRelTimeCam.LoGain: no +MJCalibrateSignal.MHCalibrationRelTimeCam.Oscillations: no +#MJCalibrateSignal.MHCalibrationRelTimeCam.SizeCheck: yes +#MJCalibrateSignal.MHCalibrationRelTimeCam.Averageing: yes +MJCalibrateSignal.MHCalibrationRelTimeCam.Nbins: 751 +MJCalibrateSignal.MHCalibrationRelTimeCam.First: -60 +MJCalibrateSignal.MHCalibrationRelTimeCam.Last: 60 +#MJCalibrateSignal.MHCalibrationRelTimeCam.NumHiGainSaturationLimit: 0.25 +#MJCalibrateSignal.MHCalibrationRelTimeCam.ProbLimit: 0.00001 +#MJCalibrateSignal.MHCalibrationRelTimeCam.OverflowLimit: 0.005 +MJCalibrateSignal.MHCalibrationRelTimeCam.PulserFrequency: 25 + +# ------------------------------------------------------------------------- +# Resource for movie encoding. For more details see IACT wiki +# ------------------------------------------------------------------------- +#MMoviePrepare.FirstSlice: 10 +#MMoviePrepare.LastSlice: 50 +#MovieFilter.Condition: MMovieData.fMax>5*MMovieData.fMedianPedestalRms +#MImgCleanStd.CleanLevel1: 8.0 +#MImgCleanStd.CleanLevel2: 4.5 +#MFEvtNumber.FileName: ganymed00001218.root +#MFEvtNumber.Selector: ThetaSquared.fVal<0.04 && DataType.fVal>0.5 +#MMovieWrite.TargetLength: 5 +#MMovieWrite.NumEvents: 500 +#MMovieWrite.Threshold: 3 +#MMovieWrite.FileName: filename.mpg + +# ------------------------------------------------------------------------- +# This is a special option for the datacheck. A task executed as last +# task in the event processing task list - in the standard software it +# is skipped +# ------------------------------------------------------------------------- +#MJCalibrateSignal.FinalFantasy: MHCamFlorian Index: trunk/Mars/resources/ganymed.rc =================================================================== --- trunk/Mars/resources/ganymed.rc (revision 19872) +++ trunk/Mars/resources/ganymed.rc (revision 19872) @@ -0,0 +1,221 @@ +# ========================================================================== +############################################################################# +# ========================================================================== +# General +# ========================================================================== +############################################################################# +# ========================================================================== + +# ------------------------------------------------------------------------- +# Use this if you want to setup the logging stream for the jobs +# (overwrites command line options) +# ------------------------------------------------------------------------- +#MLog.VerbosityLevel: 2 +#MLog.DebugLevel: 1 +#MLog.NoColors: yes + +# ========================================================================== +############################################################################# +# ========================================================================== +# Cuts +# ========================================================================== +############################################################################# +# ========================================================================== + +# ------------------------------------------------------------------------- +# Use this if you want to write output to somewhere here you can give +# the output path +# ------------------------------------------------------------------------- +#MJCuts.PathOut: . + +# ------------------------------------------------------------------------- +# Use this to setup binnings. For more details see: MBinning::ReadEnv +# ------------------------------------------------------------------------- +BinningAlpha.Raw: 36 0 90 +BinningFalseSource.Raw: 91 -1.5 1.5 +BinningWidth.Raw: 50 0 1 +BinningLength.Raw: 25 0 1 +BinningDist.Raw: 25 0 2 +BinningAsym.Raw: 51 -0.7 0.7 +BinningM3Long.Raw: 51 -0.7 0.7 +BinningM3Trans.Raw: 51 -0.7 0.7 +BinningSlope.Raw: 51 -7.5 7.5 + +# ------------------------------------------------------------------------- +# Using the starguider for pointing correction. +# To switch it off use "MPointingDevCalc.MaxAbsDev: -1" +# For details see the class reference of MPointingDevCalc +# ------------------------------------------------------------------------- +#MPointingDevCalc.MaxAbsDev: 15 +#MPointingDevCalc.NumMinStars: 8 +#MPointingDevCalc.NsbLevel: 3.0 +#MPointingDevCalc.NsbMin: 30 +#MPointingDevCalc.NsbMax: 60 +#MPointingDevCalc.Dx -0.001 +#MPointingDevCalc.Dy -0.004 +#MPointingDevCalc.MaxAge 1.0 + +MPointingDevCalc.PointingModels: 1 81855 83722 89180 267253 291104 328198 1000534 + +# ------------------------------------------------------------------------- +# Setup misfocussing correction (see class reference for MSrcPosCorrect) +# ------------------------------------------------------------------------- +#MSrcPosCorrect.Dx: 0 +#MSrcPosCorrect.Dy: 0 + +# ------------------------------------------------------------------------- +# Setup fit mode and ranges +# ------------------------------------------------------------------------- +#MAlphaFitter.SignalIntegralMax: 15 +#MAlphaFitter.SignalFitMax: 25 +#MAlphaFitter.BackgroundFitMin: 35 +#MAlphaFitter.BackgroundFitMax: 75 +#MAlphaFitter.ScaleMin: 35 +#MAlphaFitter.ScaleMax: 80 +#MAlphaFitter.PolynomOrder: 2 +#MAlphaFitter.ScaleMode: None + +# ------------------------------------------------------------------------- +# Define here which histogram class to use to determin the signal. +# Currently availyble: MHAlpha, MHTheta +# ------------------------------------------------------------------------- +#MJCut.NameHist: MHAlpha + +# Be carefull: The spectrum might be done with a different binning +# This might change the scale factor +#MHThetaSq.NumBinsSignal: 3 +#MHThetaSq.NumBinsTotal: 75 + +# ------------------------------------------------------------------------- +# To get a meaningfull false source plot in wobble mode use magic cuts +# and the following two lines. +# ------------------------------------------------------------------------- +MJCut.NameHistFS: MHDisp + +MHDisp.Smearing: 0.06 +#MHDisp.Wobble: Yes +#MHDisp.ScaleMin: 0.325 +#MHDisp.ScaleMax: 0.475 + +# ------------------------------------------------------------------------- +# To get/use a random source position in the off data with the distribution +# in the camera of the on-data +# ------------------------------------------------------------------------- +#MJCut.RandomSourcePosition: Yes + +# ------------------------------------------------------------------------- +# There are three cuts: +# CutQ: Cut after reading file +# Cut0: After calculation of src dependent parameters, energy estimation +# and disp calculation, before writing summary file +# Cut1: After writing summary file before filling false source plots +# Cut2: After filling false source plots before filling alpha plots +# Cut3: After filling alpha plots (eg. Alpha cuts) for image parameter +# display, eg. write only events after alpha cut to ganymed file +# ------------------------------------------------------------------------- +CutQ.Inverted: Yes +Cut0.Inverted: Yes +Cut1.Inverted: Yes +Cut2.Inverted: Yes +Cut3.Inverted: Yes + +# Quality-cuts (remove unphysical and obvious non-gamma events) +CutQ.Condition: {0} && {1} && {2} && {3} && {4} && {5} +CutQ.0: MImagePar.fNumIslands<3 +CutQ.1: MNewImagePar.fNumUsedPixels>5 +CutQ.2: MNewImagePar.fLeakage1 < 0.3 +CutQ.3: log10(MNewImagePar.fConcCOG)<-0.45+0.08*(log10(MHillas.fSize)<3.9)*(log10(MHillas.fSize)-3.9)^2 +CutQ.4: log10(MNewImagePar.fConc1) <-0.75+0.10*(log10(MHillas.fSize)<3.8)*(log10(MHillas.fSize)-3.8)^2 +CutQ.5: MHillas.fLength>-3.6*(log10(MHillas.fSize)-6.0)^2+70 + +# ------------------------------------------------------------------------- +# Such a cut could be used to exclude strange pointing positions. like +# the camera center in wobble mode. Note, that a cut in the poiting +# position directly influences the effective observation time and should +# only be applied if the events are lost anyhow. In both cases you must +# correct for the loss of observation time calculating your flux! +# ------------------------------------------------------------------------- +#Cut0.Condition: MSrcPosCam.GetDist*MGeomCam.fConvMm2Deg>0.3 + +# ---------- SETUP FOR ONOFF-MODE ----------- +MAlphaFitter.ScaleMin: 0.137 +MAlphaFitter.ScaleMax: 0.640 +MAlphaFitter.BackgroundFitMin: 0.137 +MAlphaFitter.BackgroundFitMax: 1.000 +MAlphaFitter.PolynomOrder: 1 +MAlphaFitter.ScaleMode: OffRegion +MAlphaFitter.SignalFunction: ThetaSq + +# For more details on the magic cuts see its class reference +Cut1.Condition: MFMagicCuts +Cut1.ThetaCut: None +# Parametrization of Disp +Cut1.Param0: 1.15136 +Cut1.Param8: 0.0681437 +Cut1.Param9: 2.62932 +Cut1.Param10: 1.51279 +Cut1.Param11: 0.0507821 +# Parametrization for sign of disp (m3long, slope) +Cut1.Param5: -0.07 +Cut1.Param6: 7.2 +Cut1.Param7: 0.5 +# ThetaSq-Cut +Cut1.Param1: 0.215 +# Area-Cut +Cut1.Param2: 0.215468 +Cut1.Param3: 5.63973 +Cut1.Param4: 0.0836169 + +# ---------- SETUP FOR WOBBLE-MODE ----------- +#MAlphaFitter.ScaleMode: none +#MAlphaFitter.BackgroundFitMin: 0.09 +#MAlphaFitter.BackgroundFitMax: 0.25 +#MAlphaFitter.PolynomOrder: 1 +#MAlphaFitter.SignalFunction: ThetaSq + +# ---------- SETUP TO USE RANDOM FOREST ----------- + +# To switch on random forest uncomment the following lines +#Cut1.HadronnessCut: Hadronness +#Cut1.Param13: 0.5 +#Cut1.Param14: 100 + +# To use both, hadronness and area cut do +#Cut1.HadronnessCut: All + +# ------------------------------------------------------------------------- +# This cut is applied exclusively to the plots in the tab CutT. It is +# resonable to show the events within your signal region here to +# see if your signal region is clean from strange events (e.g. sparks) +# ------------------------------------------------------------------------- +CutT.Condition: ThetaSquared.fVal < ThetaSquaredCut.fVal + +# ------------------------------------------------------------------------- +# Use this to executa a task (eg to calc hadronness) before Cut1 +# ------------------------------------------------------------------------- +# CalcHadronness: MRanForestCalc +# CalcHadronness.FileName: rf-separation.root +# CalcHadronness.NameOutput: Hadronness +# CalcHadronness.Debug: No + +# ------------------------------------------------------------------------- +# Use this to executa a task (eg to calc disp) before Cut1 +# ------------------------------------------------------------------------- +# CalcDisp: MRanForestCalc +# CalcDisp.FileName: rf-disp.root +# CalcDisp.NameOutput: Disp +# CalcDisp.Debug: No +# Cut1.CalcDisp: No + +# CalcDisp: MParameterCalc +# CalcDisp.Rule: 1.15*(1-MHillas.fWidth/MHillas.fLength) +# CalcDisp.NameOutput: Disp +# Cut1.CalcDisp: No + +# ------------------------------------------------------------------------- +# Use this to executa a task (eg to estimate energy) after all Cuts +# ------------------------------------------------------------------------- +# EstimateEnergy: MRanForestCalc +# EstimateEnergy.NameOutput: MEnergyEst +# EstimateEnergy.FileName: ranforest.root +# EstimateEnergy.Debug: No Index: trunk/Mars/resources/ganymed_onoff.rc =================================================================== --- trunk/Mars/resources/ganymed_onoff.rc (revision 19872) +++ trunk/Mars/resources/ganymed_onoff.rc (revision 19872) @@ -0,0 +1,211 @@ +# ========================================================================== +############################################################################# +# ========================================================================== +# General +# ========================================================================== +############################################################################# +# ========================================================================== + +# ------------------------------------------------------------------------- +# Use this if you want to setup the logging stream for the jobs +# (overwrites command line options) +# ------------------------------------------------------------------------- +#MLog.VerbosityLevel: 2 +#MLog.DebugLevel: 1 +#MLog.NoColors: yes + +# ========================================================================== +############################################################################# +# ========================================================================== +# Cuts +# ========================================================================== +############################################################################# +# ========================================================================== + +# ------------------------------------------------------------------------- +# Use this if you want to write output to somewhere here you can give +# the output path +# ------------------------------------------------------------------------- +#MJCuts.PathOut: . + +# ------------------------------------------------------------------------- +# Use this to setup binnings. For more details see: MBinning::ReadEnv +# ------------------------------------------------------------------------- +BinningAlpha.Raw: 36 0 90 +BinningFalseSource.Raw: 91 -1.5 1.5 +BinningWidth.Raw: 50 0 1 +BinningLength.Raw: 25 0 1 +BinningDist.Raw: 25 0 2 +BinningAsym.Raw: 51 -0.7 0.7 +BinningM3Long.Raw: 51 -0.7 0.7 +BinningM3Trans.Raw: 51 -0.7 0.7 +BinningSlope.Raw: 51 -7.5 7.5 + +# ------------------------------------------------------------------------- +# Using the starguider for pointing correction. +# To switch it off use "MPointingDevCalc.MaxAbsDev: -1" +# For details see the class reference of MPointingDevCalc +# ------------------------------------------------------------------------- +#MPointingDevCalc.MaxAbsDev: 15 +#MPointingDevCalc.NumMinStars: 8 +#MPointingDevCalc.NsbLevel: 3.0 +#MPointingDevCalc.NsbMin: 30 +#MPointingDevCalc.NsbMax: 60 +#MPointingDevCalc.Dx -0.001 +#MPointingDevCalc.Dy -0.004 +#MPointingDevCalc.MaxAge 1.0 + +MPointingDevCalc.PointingModels: 1 81855 83722 89180 267253 291104 328198 1000534 + +# ------------------------------------------------------------------------- +# Setup misfocussing correction (see class reference for MSrcPosCorrect) +# ------------------------------------------------------------------------- +#MSrcPosCorrect.Dx: 0 +#MSrcPosCorrect.Dy: 0 + +# ------------------------------------------------------------------------- +# Setup fit mode and ranges +# ------------------------------------------------------------------------- +#MAlphaFitter.SignalIntegralMax: 15 +#MAlphaFitter.SignalFitMax: 25 +#MAlphaFitter.BackgroundFitMin: 35 +#MAlphaFitter.BackgroundFitMax: 75 +#MAlphaFitter.ScaleMin: 35 +#MAlphaFitter.ScaleMax: 80 +#MAlphaFitter.PolynomOrder: 2 +#MAlphaFitter.ScaleMode: None + +# ------------------------------------------------------------------------- +# Define here which histogram class to use to determin the signal. +# Currently availyble: MHAlpha, MHTheta +# ------------------------------------------------------------------------- +#MJCut.NameHist: MHAlpha + +# Be carefull: The spectrum might be done with a different binning +# This might change the scale factor +#MHThetaSq.NumBinsSignal: 3 +#MHThetaSq.NumBinsTotal: 75 + +# ------------------------------------------------------------------------- +# To get a meaningfull false source plot in wobble mode use magic cuts +# and the following two lines. +# ------------------------------------------------------------------------- +MJCut.NameHistFS: MHDisp + +MHDisp.Smearing: 0.06 +#MHDisp.ScaleMin: 0.325 +#MHDisp.ScaleMax: 0.475 + +# ------------------------------------------------------------------------- +# To get use a random source position in the off data with the distribution +# in the camera of the on-data +# ------------------------------------------------------------------------- +MJCut.RandomSourcePosition: Yes + +# ------------------------------------------------------------------------- +# There are three cuts: +# CutQ: Cut after reading file +# Cut0: After calculation of src dependent parameters, energy estimation +# and disp calculation, before writing summary file +# Cut1: After writing summary file before filling false source plots +# Cut2: After filling false source plots before filling alpha plots +# Cut3: After filling alpha plots (eg. Alpha cuts) for image parameter +# display, eg. write only events after alpha cut to ganymed file +# ------------------------------------------------------------------------- +CutQ.Inverted: Yes +Cut0.Inverted: Yes +Cut1.Inverted: Yes +Cut2.Inverted: Yes +Cut3.Inverted: Yes + +# Quality-cuts (remove unphysical and obvious non-gamma events) +CutQ.Condition: {0} && {1} && {2} && {3} && {4} && {5} +CutQ.0: MImagePar.fNumIslands<3 +CutQ.1: MNewImagePar.fNumUsedPixels>5 +CutQ.2: MNewImagePar.fLeakage1 < 0.3 +CutQ.3: log10(MNewImagePar.fConcCOG)<-0.45+0.08*(log10(MHillas.fSize)<3.9)*(log10(MHillas.fSize)-3.9)^2 +CutQ.4: log10(MNewImagePar.fConc1) <-0.75+0.10*(log10(MHillas.fSize)<3.8)*(log10(MHillas.fSize)-3.8)^2 +CutQ.5: MHillas.fLength>-3.6*(log10(MHillas.fSize)-6.0)^2+70 + +# ------------------------------------------------------------------------- +# Such a cut could be used to exclude strange pointing positions. like +# the camera center in wobble mode. Note, that a cut in the poiting +# position directly influences the effective observation time and should +# only be applied if the events are lost anyhow. In both cases you must +# correct for the loss of observation time calculating your flux! +# ------------------------------------------------------------------------- +#Cut0.Condition: MSrcPosCam.GetDist*MGeomCam.fConvMm2Deg>0.3 + +# ---------- SETUP FOR ONOFF-MODE ----------- +MAlphaFitter.ScaleMin: 0.137 +MAlphaFitter.ScaleMax: 0.640 +MAlphaFitter.BackgroundFitMin: 0.137 +MAlphaFitter.BackgroundFitMax: 1.000 +MAlphaFitter.PolynomOrder: 1 +MAlphaFitter.ScaleMode: OffRegion +MAlphaFitter.SignalFunction: ThetaSq + +# For more details on the magic cuts see its class reference +Cut1.Condition: MFMagicCuts +Cut1.ThetaCut: None +# Parametrization of Disp +Cut1.Param0: 1.15136 +Cut1.Param8: 0.0681437 +Cut1.Param9: 2.62932 +Cut1.Param10: 1.51279 +Cut1.Param11: 0.0507821 +# Parametrization for sign of disp (m3long, slope) +Cut1.Param5: -0.07 +Cut1.Param6: 7.2 +Cut1.Param7: 0.5 +# ThetaSq-Cut +Cut1.Param1: 0.215 +# Area-Cut +Cut1.Param2: 0.215468 +Cut1.Param3: 5.63973 +Cut1.Param4: 0.0836169 + +# To switch on random forest uncomment the following lines +#Cut1.HadronnessCut: Hadronness +#Cut1.Param13: 0.5 +#Cut1.Param14: 100 + +# To use both, hadronness and area cut do +#Cut1.HadronnessCut: All + +# ------------------------------------------------------------------------- +# This cut is applied exclusively to the plots in the tab CutT. It is +# resonable to show the events within your signal region here to +# see if your signal region is clean from strange events (e.g. sparks) +# ------------------------------------------------------------------------- +CutT.Condition: ThetaSquared.fVal < ThetaSquaredCut.fVal + +# ------------------------------------------------------------------------- +# Use this to executa a task (eg to calc hadronness) before Cut1 +# ------------------------------------------------------------------------- +# CalcHadronness: MRanForestCalc +# CalcHadronness.FileName: rf-separation.root +# CalcHadronness.NameOutput: Hadronness +# CalcHadronness.Debug: No + +# ------------------------------------------------------------------------- +# Use this to executa a task (eg to calc disp) before Cut1 +# ------------------------------------------------------------------------- +# CalcDisp: MRanForestCalc +# CalcDisp.FileName: rf-disp.root +# CalcDisp.NameOutput: Disp +# CalcDisp.Debug: No +# Cut1.CalcDisp: No + +# CalcDisp: MParameterCalc +# CalcDisp.Rule: 1.15*(1-MHillas.fWidth/MHillas.fLength) +# CalcDisp.NameOutput: Disp +# Cut1.CalcDisp: No + +# ------------------------------------------------------------------------- +# Use this to executa a task (eg to estimate energy) after all Cuts +# ------------------------------------------------------------------------- +# EstimateEnergy: MRanForestCalc +# EstimateEnergy.NameOutput: MEnergyEst +# EstimateEnergy.FileName: ranforest.root +# EstimateEnergy.Debug: No Index: trunk/Mars/resources/ganymed_wobble.rc =================================================================== --- trunk/Mars/resources/ganymed_wobble.rc (revision 19872) +++ trunk/Mars/resources/ganymed_wobble.rc (revision 19872) @@ -0,0 +1,211 @@ +# ========================================================================== +############################################################################# +# ========================================================================== +# General +# ========================================================================== +############################################################################# +# ========================================================================== + +# ------------------------------------------------------------------------- +# Use this if you want to setup the logging stream for the jobs +# (overwrites command line options) +# ------------------------------------------------------------------------- +#MLog.VerbosityLevel: 2 +#MLog.DebugLevel: 1 +#MLog.NoColors: yes + +# ========================================================================== +############################################################################# +# ========================================================================== +# Cuts +# ========================================================================== +############################################################################# +# ========================================================================== + +# ------------------------------------------------------------------------- +# Use this if you want to write output to somewhere here you can give +# the output path +# ------------------------------------------------------------------------- +#MJCuts.PathOut: . + +# ------------------------------------------------------------------------- +# Use this to setup binnings. For more details see: MBinning::ReadEnv +# ------------------------------------------------------------------------- +BinningAlpha.Raw: 36 0 90 +BinningFalseSource.Raw: 91 -1.5 1.5 +BinningWidth.Raw: 50 0 1 +BinningLength.Raw: 25 0 1 +BinningDist.Raw: 25 0 2 +BinningAsym.Raw: 51 -0.7 0.7 +BinningM3Long.Raw: 51 -0.7 0.7 +BinningM3Trans.Raw: 51 -0.7 0.7 +BinningSlope.Raw: 51 -7.5 7.5 + +# ------------------------------------------------------------------------- +# Using the starguider for pointing correction. +# To switch it off use "MPointingDevCalc.MaxAbsDev: -1" +# For details see the class reference of MPointingDevCalc +# ------------------------------------------------------------------------- +#MPointingDevCalc.MaxAbsDev: 15 +#MPointingDevCalc.NumMinStars: 8 +#MPointingDevCalc.NsbLevel: 3.0 +#MPointingDevCalc.NsbMin: 30 +#MPointingDevCalc.NsbMax: 60 +#MPointingDevCalc.Dx -0.001 +#MPointingDevCalc.Dy -0.004 +#MPointingDevCalc.MaxAge 1.0 + +MPointingDevCalc.PointingModels: 1 81855 83722 89180 267253 291104 328198 1000534 + +# ------------------------------------------------------------------------- +# Setup misfocussing correction (see class reference for MSrcPosCorrect) +# ------------------------------------------------------------------------- +#MSrcPosCorrect.Dx: 0 +#MSrcPosCorrect.Dy: 0 + +# ------------------------------------------------------------------------- +# Setup fit mode and ranges +# ------------------------------------------------------------------------- +#MAlphaFitter.SignalIntegralMax: 15 +#MAlphaFitter.SignalFitMax: 25 +#MAlphaFitter.BackgroundFitMin: 35 +#MAlphaFitter.BackgroundFitMax: 75 +#MAlphaFitter.ScaleMin: 35 +#MAlphaFitter.ScaleMax: 80 +#MAlphaFitter.PolynomOrder: 2 +#MAlphaFitter.ScaleMode: None + +# ------------------------------------------------------------------------- +# Define here which histogram class to use to determin the signal. +# Currently availyble: MHAlpha, MHTheta +# ------------------------------------------------------------------------- +#MJCut.NameHist: MHAlpha + +# ------------------------------------------------------------------------- +# Define here the number of off-source regions used in wobble mode +# ------------------------------------------------------------------------- +#MJCut.NumOffSourcePositions: 3 + +# ------------------------------------------------------------------------- +# To get a meaningfull false source plot in wobble mode use magic cuts +# and the following two lines. +# ------------------------------------------------------------------------- +MJCut.NameHistFS: MHDisp + +MHDisp.Wobble: Yes +MHDisp.Smearing: 0.06 +#MHDisp.ScaleMin: 0.325 +#MHDisp.ScaleMax: 0.475 + +# ------------------------------------------------------------------------- +# Use this to set the number of bins corresponding to the signal region +# and if you use a Anti-Theta cut the second one has to be enabled +# ------------------------------------------------------------------------- +#MHPhi.NumBinsSignal: 2 +#MHPhi.UseAntiPhiCut: No + +# ------------------------------------------------------------------------- +# There are three cuts: +# CutQ: Cut after reading file +# Cut0: After calculation of src dependent parameters, energy estimation +# and disp calculation, before writing summary file +# Cut1: After writing summary file before filling false source plots +# Cut2: After filling false source plots before filling alpha plots +# Cut3: After filling alpha plots (eg. Alpha cuts) for image parameter +# display, eg. write only events after alpha cut to ganymed file +# ------------------------------------------------------------------------- +CutQ.Inverted: Yes +Cut0.Inverted: Yes +Cut1.Inverted: Yes +Cut2.Inverted: Yes +Cut3.Inverted: Yes + +# Quality-cuts (remove unphysical and obvious non-gamma events) +CutQ.Condition: {0} && {1} && {2} && {3} && {4} && {5} +CutQ.0: MImagePar.fNumIslands<3 +CutQ.1: MNewImagePar.fNumUsedPixels>5 +CutQ.2: MNewImagePar.fLeakage1 < 0.3 +CutQ.3: log10(MNewImagePar.fConcCOG)<-0.45+0.08*(log10(MHillas.fSize)<3.9)*(log10(MHillas.fSize)-3.9)^2 +CutQ.4: log10(MNewImagePar.fConc1) <-0.75+0.10*(log10(MHillas.fSize)<3.8)*(log10(MHillas.fSize)-3.8)^2 +CutQ.5: MHillas.fLength>-3.6*(log10(MHillas.fSize)-6.0)^2+70 + +# ------------------------------------------------------------------------- +# Such a cut could be used to exclude strange pointing positions. like +# the camera center in wobble mode. Note, that a cut in the poiting +# position directly influences the effective observation time and should +# only be applied if the events are lost anyhow. In both cases you must +# correct for the loss of observation time calculating your flux! +# ------------------------------------------------------------------------- +#Cut0.Condition: MSrcPosCam.GetDist*MGeomCam.fConvMm2Deg>0.3 + +# ---------- SETUP FOR WOBBLE-MODE ----------- +#MAlphaFitter.ScaleMode: none +MAlphaFitter.BackgroundFitMin: 0.09 +MAlphaFitter.BackgroundFitMax: 0.25 +MAlphaFitter.PolynomOrder: 1 +MAlphaFitter.SignalFunction: ThetaSq + +# For more details on the magic cuts see its class reference +Cut1.Condition: MFMagicCuts +Cut1.ThetaCut: None +# Parametrization of Disp +Cut1.Param0: 1.15136 +Cut1.Param8: 0.0681437 +Cut1.Param9: 2.62932 +Cut1.Param10: 1.51279 +Cut1.Param11: 0.0507821 +# Parametrization for sign of disp (m3long, slope) +Cut1.Param5: -0.07 +Cut1.Param6: 7.2 +Cut1.Param7: 0.5 +# ThetaSq-Cut +Cut1.Param1: 0.215 +# Area-Cut +Cut1.Param2: 0.215468 +Cut1.Param3: 5.63973 +Cut1.Param4: 0.0836169 + +# To switch on random forest uncomment the following lines +#Cut1.HadronnessCut: Hadronness +#Cut1.Param13: 0.5 +#Cut1.Param14: 100 + +# To use both, hadronness and area cut do +#Cut1.HadronnessCut: All + +# ------------------------------------------------------------------------- +# This cut is applied exclusively to the plots in the tab CutT. It is +# resonable to show the events within your signal region here to +# see if your signal region is clean from strange events (e.g. sparks) +# ------------------------------------------------------------------------- +CutT.Condition: ThetaSquared.fVal < ThetaSquaredCut.fVal + +# ------------------------------------------------------------------------- +# Use this to executa a task (eg to calc hadronness) before Cut1 +# ------------------------------------------------------------------------- +# CalcHadronness: MRanForestCalc +# CalcHadronness.FileName: rf-separation.root +# CalcHadronness.NameOutput: Hadronness +# CalcHadronness.Debug: No + +# ------------------------------------------------------------------------- +# Use this to executa a task (eg to calc disp) before Cut1 +# ------------------------------------------------------------------------- +# CalcDisp: MRanForestCalc +# CalcDisp.FileName: rf-disp.root +# CalcDisp.NameOutput: Disp +# CalcDisp.Debug: No +# Cut1.CalcDisp: No + +# CalcDisp: MParameterCalc +# CalcDisp.Rule: 1.15*(1-MHillas.fWidth/MHillas.fLength) +# CalcDisp.NameOutput: Disp +# Cut1.CalcDisp: No + +# ------------------------------------------------------------------------- +# Use this to executa a task (eg to estimate energy) after all Cuts +# ------------------------------------------------------------------------- +# EstimateEnergy: MRanForestCalc +# EstimateEnergy.NameOutput: MEnergyEst +# EstimateEnergy.FileName: ranforest.root +# EstimateEnergy.Debug: No Index: trunk/Mars/resources/mars.rc =================================================================== --- trunk/Mars/resources/mars.rc (revision 19872) +++ trunk/Mars/resources/mars.rc (revision 19872) @@ -0,0 +1,106 @@ +# ========================================================================== +############################################################################# +# ========================================================================== +# +# This is a resource file for the Mars executable. You can here setup +# setup the tasks of the eventloops executed in the Mars executable - +# if supported by these tasks. +# +# ========================================================================== +############################################################################# +# ========================================================================== + +# ------------------------------------------------------------------------- +# Use this if you want to setup the logging stream for the jobs +# (overwrites command line options) +# ------------------------------------------------------------------------- +#MLog.VerbosityLevel: 2 +#MLog.DebugLevel: 1 +#MLog.NoColors: yes + +# ------------------------------------------------------------------------- +# Define here which signal extractor you are going to use for the +# whole calibration chain +# ------------------------------------------------------------------------- +#ExtractSignal: MExtractTimeAndChargeDigitalFilter +#ExtractSignal.WeightsFile: msignal/MC_weights46.dat +# ------------------------------------------------------------------------- +# Define here parameters valid for all extractors above +# ------------------------------------------------------------------------- +#ExtractSignal.LoGainSwitch: 120 +#ExtractSignal.LoGainStartShift: -2.5 +# ------------------------------------------------------------------------- +# Define here parameters valid for the spline: +# ------------------------------------------------------------------------- +#ExtractSignal.HiGainFirst: 0 +#ExtractSignal.HiGainLast: 14 +#ExtractSignal.LoGainFirst: 1 +#ExtractSignal.LoGainLast: 14 +#ExtractSignal.Resolution: 0.05 +#ExtractSignal.RiseTimeHiGain: 0.64 +#ExtractSignal.FallTimeHiGain: 0.76 +#ExtractSignal.LoGainStretch: 1.5 +#ExtractSignal.OffsetLoGain: 1.3 +#ExtractSignal.ExtractionType: Integral +#ExtractSignal.ExtractionType: Amplitude +# ------------------------------------------------------------------------- +# Define here parameters valid for the Digital Filter: +# ------------------------------------------------------------------------- +#ExtractSignal.HiGainFirst: 0 +#ExtractSignal.HiGainLast: 16 +#ExtractSignal.LoGainFirst: 1 +#ExtractSignal.LoGainLast: 14 +#ExtractSignal.OffsetLoGain: 0.95 + +# ========================================================================= +# To look into MUX data use the following setup +# ========================================================================= +#ExtractSignal: MExtractTimeAndChargeSpline +#ExtractSignal.RiseTimeHiGain: 4.3 +#ExtractSignal.FallTimeHiGain: 5.1 +#ExtractSignal.HiGainFirst: 0 +#ExtractSignal.HiGainLast: 49 + +# ------------------------------------------------------------------------- +# In the case of the Digital filter you have to overwrite the default +# weights file depending on what you are extracting +# ------------------------------------------------------------------------- +#ExtractSignal.WeightsFile: msignal/cosmics_weights.dat + +# ------------------------------------------------------------------------- +# Setup level for determining of bad pixels +# ------------------------------------------------------------------------- +#MBadPixelsCalc.PedestalLevel: 3.0 +#MBadPixelsCalc.PedestalLevelVariance: 3.0 +#MBadPixelsTreat.NumMinNeighbors: 3 +#MBadPixelsTreat.UseInterpolation: yes +#MBadPixelsTreat.ProcessPedestalEvt: yes +#MBadPixelsTreat.ProcessPedestalRun: no +#MBadPixelsTreat.ProcessTimes: yes +#MBadPixelsTreat.UseCentralPixel: no +#MBadPixelsTreat.HardTreatment: no + +# ------------------------------------------------------------------------- +# This is how you can setup the Image cleaning used in the Camera Display +# e.g. Standard, Absolute, Scaled, Time +# ------------------------------------------------------------------------- +MImgCleanStd.CleanMethod: absolute +#MImgCleanStd.CleanLevel0: 6.0 +MImgCleanStd.CleanLevel1: 6.0 +MImgCleanStd.CleanLevel2: 3.0 +MImgCleanStd.TimeLevel1: 17.5 +MImgCleanStd.TimeLevel2: 17.5 +MImgCleanStd.PostCleanType: 3 + +#MImgCleanStd.CleanRings: 1 +#MImgCleanStd.KeepSinglePixels: No + +# ------------------------------------------------------------------------- +# Setup a cut to choose the events which should be displayed +# (Use this cuts to display only good Muon candidates) +# ------------------------------------------------------------------------- +#Cut.Inverted: Yes +#Cut.Condition: {0} && {1} && {2} +#Cut.0: MHillas.fSize>150 +#Cut.1: (MMuonSearchPar.fRadius>180) && (MMuonSearchPar.fRadius<400) && (MMuonSearchPar.fDeviation<45) +#Cut.2: (MMuonCalibPar.fArcPhi>190) && (MMuonSearchPar.fDeviation<35) && (MMuonCalibPar.fArcWidth<0.20) && (MMuonCalibPar.fArcWidth>0.04) Index: trunk/Mars/resources/sponde.rc =================================================================== --- trunk/Mars/resources/sponde.rc (revision 19872) +++ trunk/Mars/resources/sponde.rc (revision 19872) @@ -0,0 +1,89 @@ +# -------------------- Define your energy estimator ------------------- + +# Manually optimized (for example using trainenergy.C) +EstimateEnergy.Rule: MHillas.fSize^0.906*pow(10, -0.831*(MHillasSrc.fDist*3.37e-3-0.161)^2*(MHillasSrc.fDist*3.37e-3>0.161) + 2.21*(MHillasSrc.fDist*3.37e-3-0.493)^2*(MHillasSrc.fDist*3.37e-3>0.493) + 8.18e-6*MPointingPos.fZd^2.82 - 0.1*MNewImagePar.fLeakage1^0.2) + +# You can use instead the Random Forest energy estimator +# (for example trained by trainenergy.C) +#EstimateEnergy: MRanForestCalc +#EstimateEnergy.NameOutput: MEnergyEst +#EstimateEnergy.FileName: rf-energy.root +#EstimateEnergy.Debug: No + +# ------------------------ Define your binnings ----------------------- + +# The energy binning is mandatory +BinningEnergyEst.Raw: 18 53 35800 log + +# Binnings used for the Energy- and Disp-tab +#BinningSize.Raw: 50 10 100000 log +#BinningLeakage.Raw: 50 0 0.3 lin +#BinningEnergyResidual.Raw: 51 -1 1 lin +#BinningResidualDist.Raw: 51 -1 1 lin + +# A change of the theta binning must be done in ganymed (see below) :( +#BinningTheta.Raw: 101 -0.005 0.505 asin + +# --------------- Setup the weighting for your spectrum --------------- +# For more details see McSpectrumWeight::ReadEnv + +# A simple power law with a differential spectral index of -4 +#MMcSpectrumWeight.NewSlope: -4 + +# Give a formula for the spectrum. Capital X is a place-holder +# for the energy (eg. the crab-spectrum from our publication) +#MMcSpectrumWeight.Formula: pow(X/300, -2.31-0.26*log10(X/300)) + +# ----------------------- Additional setup ---------------------------- + +# In the loop filling the ThetaSq histogram the theta cut is switched off +# and the theta-cut is implicitly done in the histogram. Reading the +# Monte Carlo for determining the collection area there is no theta- +# histogram involved, thus you need the theta cut switched on +Spectrum.Cut1.ThetaCut: On + +# You may want to change the scale region to a different region than +# previously used in your ganymed-file. This might be usefull for +# on/off-data. Note, that the result can also depend on what else +# was setup for your MAlphaFitter in ganymed. +#MAlphaFitter.ScaleMin: 0.137 +#MAlphaFitter.ScaleMax: 0.500 + +# You may use this to change the scale and/or fit region for +# determination of the excess vs energy (can be usefull in On/Off mode) +#MAlphaFitter.ScaleMode: Fixed + +# You may want to allow scaling of each energy bin individually. (This +# can be useful if on- and off- data in wobble mode do not agree well) +#MAlphaFitter.ScaleMode: background +#MAlphaFitter.BackgroundFitMin: 0.09 +#MAlphaFitter.BackgroundFitMax: 0.25 + +# ------------------------- Additional cut ---------------------------- + +# You can apply a cut in addition to all other cuts to your data and +# Monte Carlos. This can be usefull for example to setup a second +# MFMagicCuts with a tighter behaviour assuming that your ganymed +# was written with a really loose cut. +#CutS.Inverted: Yes +#CutS.Condition: MFMagicCuts +#CutS.ThetaCut: None +#CutS.HadronnesCut: Area +# Parametrization of Area +#CutS.Param2: 0.215468 +#CutS.Param3: 5.63973 +#CutS.Param4: 0.0836169 + +# --------------------------- Other binnings -------------------------- + +# Binnings that are used by sponde but can not/shoud not be changed +# because they are taken from the ganymed file +#BinningTheta +#BinningWidth +#BinningLength +#BinningDist +#BinningM3Long +#BinningM3Trans +#BinningSlope +#BinningAsym +#BinningConc1 Index: trunk/Mars/resources/sponde_onoff.rc =================================================================== --- trunk/Mars/resources/sponde_onoff.rc (revision 19872) +++ trunk/Mars/resources/sponde_onoff.rc (revision 19872) @@ -0,0 +1,89 @@ +# -------------------- Define your energy estimator ------------------- + +# Manually optimized (for example using trainenergy.C) +EstimateEnergy.Rule: MHillas.fSize^0.906*pow(10, -0.831*(MHillasSrc.fDist*3.37e-3-0.161)^2*(MHillasSrc.fDist*3.37e-3>0.161) + 2.21*(MHillasSrc.fDist*3.37e-3-0.493)^2*(MHillasSrc.fDist*3.37e-3>0.493) + 8.18e-6*MPointingPos.fZd^2.82 - 0.1*MNewImagePar.fLeakage1^0.2) + +# You can use instead the Random Forest energy estimator +# (for example trained by trainenergy.C) +#EstimateEnergy: MRanForestCalc +#EstimateEnergy.NameOutput: MEnergyEst +#EstimateEnergy.FileName: rf-energy.root +#EstimateEnergy.Debug: No + +# ------------------------ Define your binnings ----------------------- + +# The energy binning is mandatory +BinningEnergyEst.Raw: 18 53 35800 log + +# Binnings used for the Energy- and Disp-tab +#BinningSize.Raw: 50 10 100000 log +#BinningLeakage.Raw: 50 0 0.3 lin +#BinningEnergyResidual.Raw: 51 -1 1 lin +#BinningResidualDist.Raw: 51 -1 1 lin + +# A change of the theta binning must be done in ganymed (see below) :( +#BinningTheta.Raw: 101 -0.005 0.505 asin + +# --------------- Setup the weighting for your spectrum --------------- +# For more details see MMcSpectrumWeight::ReadEnv + +# A simple power law with a differential spectral index of -4 +#MMcSpectrumWeight.NewSlope: -4 + +# Give a formula for the spectrum. Capital X is a place-holder +# for the energy (eg. the crab-spectrum from our publication) +#MMcSpectrumWeight.Formula: pow(X/300, -2.31-0.26*log10(X/300)) + +# ----------------------- Additional setup ---------------------------- + +# In the loop filling the ThetaSq histogram the theta cut is switched off +# and the theta-cut is implicitly done in the histogram. Reading the +# Monte Carlo for determining the collection area there is no theta- +# histogram involved, thus you need the theta cut switched on +Spectrum.Cut1.ThetaCut: On + +# You may want to change the scale region to a different region than +# previously used in your ganymed-file. This might be usefull for +# on/off-data. Note, that the result can also depend on what else +# was setup for your MAlphaFitter in ganymed. +#MAlphaFitter.ScaleMin: 0.137 +#MAlphaFitter.ScaleMax: 0.500 + +# You may use this to change the scale and/or fit region for +# determination of the excess vs energy (can be usefull in On/Off mode) +#MAlphaFitter.ScaleMode: Fixed + +# You may want to allow scaling of each energy bin individually. (This +# can be useful if on- and off- data in wobble mode do not agree well) +#MAlphaFitter.ScaleMode: background +#MAlphaFitter.BackgroundFitMin: 0.09 +#MAlphaFitter.BackgroundFitMax: 0.25 + +# ------------------------- Additional cut ---------------------------- + +# You can apply a cut in addition to all other cuts to your data and +# MOnte Carlos. This can be usefull for example to setup a second +# MFMagicCuts with a tighter behaviour assuming that your ganymed +# was written with a really loose cut. +#CutS.Inverted: Yes +#CutS.Condition: MFMagicCuts +#CutS.ThetaCut: None +#CutS.HadronnesCut: Area +# Parametrization of Area +#CutS.Param2: 0.215468 +#CutS.Param3: 5.63973 +#CutS.Param4: 0.0836169 + +# --------------------------- Other binnings -------------------------- + +# Binnings that are used by sponde but can not/shoud not be changed +# because they are taken from the ganymed file +#BinningTheta +#BinningWidth +#BinningLength +#BinningDist +#BinningM3Long +#BinningM3Trans +#BinningSlope +#BinningAsym +#BinningConc1 Index: trunk/Mars/resources/star.rc =================================================================== --- trunk/Mars/resources/star.rc (revision 19872) +++ trunk/Mars/resources/star.rc (revision 19872) @@ -0,0 +1,96 @@ +# ========================================================================== +############################################################################# +# ========================================================================== +# General +# ========================================================================== +############################################################################# +# ========================================================================== + +# ------------------------------------------------------------------------- +# Use this if you want to setup the logging stream for the jobs +# (overwrites command line options) +# ------------------------------------------------------------------------- +#MLog.VerbosityLevel: 2 +#MLog.DebugLevel: 1 +#MLog.NoColors: yes + +# ========================================================================== +############################################################################# +# ========================================================================== +# Star +# ========================================================================== +############################################################################# +# ========================================================================== + +# ------------------------------------------------------------------------- +# 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. +# ------------------------------------------------------------------------- +#MJStar.PathOut: . +#MJStar.PathIn: . + +# ------------------------------------------------------------------------- +# Configure MJStar +# ------------------------------------------------------------------------- +#MJStar.MaxEvents: 10000 +#MJStar.Overwrite: yes,no + +# ------------------------------------------------------------------------- +# Use this if you want to change the parameters of the software trigger +# ------------------------------------------------------------------------- +#FilterSwTrigger.Threshold: 5 +#FilterSwTrigger.NumNeighbors: 4 +#FilterSwTrigger.TimeWindow: 0.5 + +# ------------------------------------------------------------------------- +# Setup the image cleaning here +# ------------------------------------------------------------------------- +MImgCleanStd.CleanMethod: absolute +#MImgCleanStd.CleanLevel0: 6.0 +MImgCleanStd.CleanLevel1: 6.0 +MImgCleanStd.CleanLevel2: 3.0 +MImgCleanStd.TimeLevel1: 17.5 +MImgCleanStd.TimeLevel2: 17.5 +MImgCleanStd.PostCleanType: 3 + +#MImgCleanStd.CleanRings: 3 +#MImgCleanStd.KeepIsolatedPixels: Yes +#MImgCleanStd.RecoverIsolatedPixles: Yes + +# ------------------------------------------------------------------------- +# Setup the cut applied to the Sparkless tab +# ------------------------------------------------------------------------- + +# This is the simple default +#SparkCut.Condition: log10(MNewImagePar.fConc1) < (-0.371)*log10(MHillas.fSize) + 0.596 + +# This is the more complicated cut from ganymed.rc +#SparkCut.Condition: ({0} || {1}) && {2} +#SparkCut.0: MImagePar.fNumSatPixelsHG < 1 +#SparkCut.1: MHillas.GetArea*(MGeomCam.fConvMm2Deg^2) > (0.003*MImagePar.fNumSatPixelsHG) + 0.0325 +#SparkCut.2: log10(MNewImagePar.fConc1) < (-0.371)*log10(MHillas.fSize) + 0.596 + +# ------------------------------------------------------------------------- +# setup the image parameter calculation here +# ------------------------------------------------------------------------- +#MHillasCalc.IdxIsland: 0 + +# ------------------------------------------------------------------------- +# Setup or switch off the muon analysis here +# ------------------------------------------------------------------------- +#MJStar.MuonAnalysis: On + +#MMuonSetup.Margin: 0.2 +#MMuonSetup.ThresholdArcPhi: 30 +#MMuonSetup.ThresholdArcWidth: 2 + +#BinningRadius.Raw: 20 0.5 1.5 +#BinningArcWidth.Raw: 60 0.0 0.3 +#BinningRingBroadening.Raw: 20 0.5 1.5 +#BinningSizeVsArcRadius.Raw: 20 0.5 1.5 +#BinningsArcPhi.Raw: 20 -180 180 +#BinningsMuonWidth.Raw: 28 0.3 1.7 + +#MuonPreCut.Condition: MHillas.fSize>150 +#MuonSearcCut.Condition: (MMuonSearchPar.fRadius>180) && (MMuonSearchPar.fRadius<400) && (MMuonSearchPar.fDeviation<45) +#MuonFinalCut.Condition: (MMuonCalibPar.fArcPhi>190) && (MMuonSearchPar.fDeviation<35) && (MMuonCalibPar.fArcWidth<0.20) && (MMuonCalibPar.fArcWidth>0.04) Index: trunk/Mars/resources/star_mux.rc =================================================================== --- trunk/Mars/resources/star_mux.rc (revision 19872) +++ trunk/Mars/resources/star_mux.rc (revision 19872) @@ -0,0 +1,101 @@ +# ========================================================================== +############################################################################# +# ========================================================================== +# General +# ========================================================================== +############################################################################# +# ========================================================================== + +# ------------------------------------------------------------------------- +# Use this if you want to setup the logging stream for the jobs +# (overwrites command line options) +# ------------------------------------------------------------------------- +#MLog.VerbosityLevel: 2 +#MLog.DebugLevel: 1 +#MLog.NoColors: yes + +# ========================================================================== +############################################################################# +# ========================================================================== +# Star +# ========================================================================== +############################################################################# +# ========================================================================== + +# ------------------------------------------------------------------------- +# 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. +# ------------------------------------------------------------------------- +#MJStar.PathOut: . +#MJStar.PathIn: . + +# ------------------------------------------------------------------------- +# Configure MJStar +# ------------------------------------------------------------------------- +#MJStar.MaxEvents: 10000 +#MJStar.Overwrite: yes,no + +# ------------------------------------------------------------------------- +# Use this if you want to change the parameters of the software trigger +# ------------------------------------------------------------------------- +#FilterSwTrigger.Threshold: 5 +#FilterSwTrigger.NumNeighbors: 4 +#FilterSwTrigger.TimeWindow: 0.5 + +# ------------------------------------------------------------------------- +# It seems that the MUX FADCs have s lightly higher dead time... +# ------------------------------------------------------------------------- +MHEffectiveOnTime.FirstBin: 4 + +# ------------------------------------------------------------------------- +# Setup the image cleaning here +# ------------------------------------------------------------------------- +MImgCleanStd.CleanMethod: absolute +#MImgCleanStd.CleanLevel0: 6.0 +MImgCleanStd.CleanLevel1: 6.0 +MImgCleanStd.CleanLevel2: 3.0 +MImgCleanStd.TimeLevel1: 1.75 +MImgCleanStd.TimeLevel2: 1.75 +MImgCleanStd.PostCleanType: 3 + +#MImgCleanStd.CleanRings: 3 +#MImgCleanStd.KeepIsolatedPixels: Yes +#MImgCleanStd.RecoverIsolatedPixles: Yes + +# ------------------------------------------------------------------------- +# Setup the cut applied to the Sparkless tab +# ------------------------------------------------------------------------- + +# This is the simple default +#SparkCut.Condition: log10(MNewImagePar.fConc1) < (-0.371)*log10(MHillas.fSize) + 0.596 + +# This is the more complicated cut from ganymed.rc +#SparkCut.Condition: ({0} || {1}) && {2} +#SparkCut.0: MImagePar.fNumSatPixelsHG < 1 +#SparkCut.1: MHillas.GetArea*(MGeomCam.fConvMm2Deg^2) > (0.003*MImagePar.fNumSatPixelsHG) + 0.0325 +#SparkCut.2: log10(MNewImagePar.fConc1) < (-0.371)*log10(MHillas.fSize) + 0.596 + +# ------------------------------------------------------------------------- +# setup the image parameter calculation here +# ------------------------------------------------------------------------- +#MHillasCalc.IdxIsland: 0 + +# ------------------------------------------------------------------------- +# Setup or switch off the muon analysis here +# ------------------------------------------------------------------------- +#MJStar.MuonAnalysis: On + +#MMuonSetup.Margin: 0.2 +#MMuonSetup.ThresholdArcPhi: 30 +#MMuonSetup.ThresholdArcWidth: 2 + +#BinningRadius.Raw: 20 0.5 1.5 +#BinningArcWidth.Raw: 60 0.0 0.3 +#BinningRingBroadening.Raw: 20 0.5 1.5 +#BinningSizeVsArcRadius.Raw: 20 0.5 1.5 +#BinningsArcPhi.Raw: 20 -180 180 +#BinningsMuonWidth.Raw: 28 0.3 1.7 + +#MuonPreCut.Condition: MHillas.fSize>150 +#MuonSearcCut.Condition: (MMuonSearchPar.fRadius>180) && (MMuonSearchPar.fRadius<400) && (MMuonSearchPar.fDeviation<45) +#MuonFinalCut.Condition: (MMuonCalibPar.fArcPhi>190) && (MMuonSearchPar.fDeviation<35) && (MMuonCalibPar.fArcWidth<0.20) && (MMuonCalibPar.fArcWidth>0.04) Index: trunk/Mars/sponde.rc =================================================================== --- trunk/Mars/sponde.rc (revision 19871) +++ (revision ) @@ -1,89 +1,0 @@ -# -------------------- Define your energy estimator ------------------- - -# Manually optimized (for example using trainenergy.C) -EstimateEnergy.Rule: MHillas.fSize^0.906*pow(10, -0.831*(MHillasSrc.fDist*3.37e-3-0.161)^2*(MHillasSrc.fDist*3.37e-3>0.161) + 2.21*(MHillasSrc.fDist*3.37e-3-0.493)^2*(MHillasSrc.fDist*3.37e-3>0.493) + 8.18e-6*MPointingPos.fZd^2.82 - 0.1*MNewImagePar.fLeakage1^0.2) - -# You can use instead the Random Forest energy estimator -# (for example trained by trainenergy.C) -#EstimateEnergy: MRanForestCalc -#EstimateEnergy.NameOutput: MEnergyEst -#EstimateEnergy.FileName: rf-energy.root -#EstimateEnergy.Debug: No - -# ------------------------ Define your binnings ----------------------- - -# The energy binning is mandatory -BinningEnergyEst.Raw: 18 53 35800 log - -# Binnings used for the Energy- and Disp-tab -#BinningSize.Raw: 50 10 100000 log -#BinningLeakage.Raw: 50 0 0.3 lin -#BinningEnergyResidual.Raw: 51 -1 1 lin -#BinningResidualDist.Raw: 51 -1 1 lin - -# A change of the theta binning must be done in ganymed (see below) :( -#BinningTheta.Raw: 101 -0.005 0.505 asin - -# --------------- Setup the weighting for your spectrum --------------- -# For more details see McSpectrumWeight::ReadEnv - -# A simple power law with a differential spectral index of -4 -#MMcSpectrumWeight.NewSlope: -4 - -# Give a formula for the spectrum. Capital X is a place-holder -# for the energy (eg. the crab-spectrum from our publication) -#MMcSpectrumWeight.Formula: pow(X/300, -2.31-0.26*log10(X/300)) - -# ----------------------- Additional setup ---------------------------- - -# In the loop filling the ThetaSq histogram the theta cut is switched off -# and the theta-cut is implicitly done in the histogram. Reading the -# Monte Carlo for determining the collection area there is no theta- -# histogram involved, thus you need the theta cut switched on -Spectrum.Cut1.ThetaCut: On - -# You may want to change the scale region to a different region than -# previously used in your ganymed-file. This might be usefull for -# on/off-data. Note, that the result can also depend on what else -# was setup for your MAlphaFitter in ganymed. -#MAlphaFitter.ScaleMin: 0.137 -#MAlphaFitter.ScaleMax: 0.500 - -# You may use this to change the scale and/or fit region for -# determination of the excess vs energy (can be usefull in On/Off mode) -#MAlphaFitter.ScaleMode: Fixed - -# You may want to allow scaling of each energy bin individually. (This -# can be useful if on- and off- data in wobble mode do not agree well) -#MAlphaFitter.ScaleMode: background -#MAlphaFitter.BackgroundFitMin: 0.09 -#MAlphaFitter.BackgroundFitMax: 0.25 - -# ------------------------- Additional cut ---------------------------- - -# You can apply a cut in addition to all other cuts to your data and -# Monte Carlos. This can be usefull for example to setup a second -# MFMagicCuts with a tighter behaviour assuming that your ganymed -# was written with a really loose cut. -#CutS.Inverted: Yes -#CutS.Condition: MFMagicCuts -#CutS.ThetaCut: None -#CutS.HadronnesCut: Area -# Parametrization of Area -#CutS.Param2: 0.215468 -#CutS.Param3: 5.63973 -#CutS.Param4: 0.0836169 - -# --------------------------- Other binnings -------------------------- - -# Binnings that are used by sponde but can not/shoud not be changed -# because they are taken from the ganymed file -#BinningTheta -#BinningWidth -#BinningLength -#BinningDist -#BinningM3Long -#BinningM3Trans -#BinningSlope -#BinningAsym -#BinningConc1 Index: trunk/Mars/sponde_onoff.rc =================================================================== --- trunk/Mars/sponde_onoff.rc (revision 19871) +++ (revision ) @@ -1,89 +1,0 @@ -# -------------------- Define your energy estimator ------------------- - -# Manually optimized (for example using trainenergy.C) -EstimateEnergy.Rule: MHillas.fSize^0.906*pow(10, -0.831*(MHillasSrc.fDist*3.37e-3-0.161)^2*(MHillasSrc.fDist*3.37e-3>0.161) + 2.21*(MHillasSrc.fDist*3.37e-3-0.493)^2*(MHillasSrc.fDist*3.37e-3>0.493) + 8.18e-6*MPointingPos.fZd^2.82 - 0.1*MNewImagePar.fLeakage1^0.2) - -# You can use instead the Random Forest energy estimator -# (for example trained by trainenergy.C) -#EstimateEnergy: MRanForestCalc -#EstimateEnergy.NameOutput: MEnergyEst -#EstimateEnergy.FileName: rf-energy.root -#EstimateEnergy.Debug: No - -# ------------------------ Define your binnings ----------------------- - -# The energy binning is mandatory -BinningEnergyEst.Raw: 18 53 35800 log - -# Binnings used for the Energy- and Disp-tab -#BinningSize.Raw: 50 10 100000 log -#BinningLeakage.Raw: 50 0 0.3 lin -#BinningEnergyResidual.Raw: 51 -1 1 lin -#BinningResidualDist.Raw: 51 -1 1 lin - -# A change of the theta binning must be done in ganymed (see below) :( -#BinningTheta.Raw: 101 -0.005 0.505 asin - -# --------------- Setup the weighting for your spectrum --------------- -# For more details see MMcSpectrumWeight::ReadEnv - -# A simple power law with a differential spectral index of -4 -#MMcSpectrumWeight.NewSlope: -4 - -# Give a formula for the spectrum. Capital X is a place-holder -# for the energy (eg. the crab-spectrum from our publication) -#MMcSpectrumWeight.Formula: pow(X/300, -2.31-0.26*log10(X/300)) - -# ----------------------- Additional setup ---------------------------- - -# In the loop filling the ThetaSq histogram the theta cut is switched off -# and the theta-cut is implicitly done in the histogram. Reading the -# Monte Carlo for determining the collection area there is no theta- -# histogram involved, thus you need the theta cut switched on -Spectrum.Cut1.ThetaCut: On - -# You may want to change the scale region to a different region than -# previously used in your ganymed-file. This might be usefull for -# on/off-data. Note, that the result can also depend on what else -# was setup for your MAlphaFitter in ganymed. -#MAlphaFitter.ScaleMin: 0.137 -#MAlphaFitter.ScaleMax: 0.500 - -# You may use this to change the scale and/or fit region for -# determination of the excess vs energy (can be usefull in On/Off mode) -#MAlphaFitter.ScaleMode: Fixed - -# You may want to allow scaling of each energy bin individually. (This -# can be useful if on- and off- data in wobble mode do not agree well) -#MAlphaFitter.ScaleMode: background -#MAlphaFitter.BackgroundFitMin: 0.09 -#MAlphaFitter.BackgroundFitMax: 0.25 - -# ------------------------- Additional cut ---------------------------- - -# You can apply a cut in addition to all other cuts to your data and -# MOnte Carlos. This can be usefull for example to setup a second -# MFMagicCuts with a tighter behaviour assuming that your ganymed -# was written with a really loose cut. -#CutS.Inverted: Yes -#CutS.Condition: MFMagicCuts -#CutS.ThetaCut: None -#CutS.HadronnesCut: Area -# Parametrization of Area -#CutS.Param2: 0.215468 -#CutS.Param3: 5.63973 -#CutS.Param4: 0.0836169 - -# --------------------------- Other binnings -------------------------- - -# Binnings that are used by sponde but can not/shoud not be changed -# because they are taken from the ganymed file -#BinningTheta -#BinningWidth -#BinningLength -#BinningDist -#BinningM3Long -#BinningM3Trans -#BinningSlope -#BinningAsym -#BinningConc1 Index: trunk/Mars/star.rc =================================================================== --- trunk/Mars/star.rc (revision 19871) +++ (revision ) @@ -1,96 +1,0 @@ -# ========================================================================== -############################################################################# -# ========================================================================== -# General -# ========================================================================== -############################################################################# -# ========================================================================== - -# ------------------------------------------------------------------------- -# Use this if you want to setup the logging stream for the jobs -# (overwrites command line options) -# ------------------------------------------------------------------------- -#MLog.VerbosityLevel: 2 -#MLog.DebugLevel: 1 -#MLog.NoColors: yes - -# ========================================================================== -############################################################################# -# ========================================================================== -# Star -# ========================================================================== -############################################################################# -# ========================================================================== - -# ------------------------------------------------------------------------- -# 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. -# ------------------------------------------------------------------------- -#MJStar.PathOut: . -#MJStar.PathIn: . - -# ------------------------------------------------------------------------- -# Configure MJStar -# ------------------------------------------------------------------------- -#MJStar.MaxEvents: 10000 -#MJStar.Overwrite: yes,no - -# ------------------------------------------------------------------------- -# Use this if you want to change the parameters of the software trigger -# ------------------------------------------------------------------------- -#FilterSwTrigger.Threshold: 5 -#FilterSwTrigger.NumNeighbors: 4 -#FilterSwTrigger.TimeWindow: 0.5 - -# ------------------------------------------------------------------------- -# Setup the image cleaning here -# ------------------------------------------------------------------------- -MImgCleanStd.CleanMethod: absolute -#MImgCleanStd.CleanLevel0: 6.0 -MImgCleanStd.CleanLevel1: 6.0 -MImgCleanStd.CleanLevel2: 3.0 -MImgCleanStd.TimeLevel1: 17.5 -MImgCleanStd.TimeLevel2: 17.5 -MImgCleanStd.PostCleanType: 3 - -#MImgCleanStd.CleanRings: 3 -#MImgCleanStd.KeepIsolatedPixels: Yes -#MImgCleanStd.RecoverIsolatedPixles: Yes - -# ------------------------------------------------------------------------- -# Setup the cut applied to the Sparkless tab -# ------------------------------------------------------------------------- - -# This is the simple default -#SparkCut.Condition: log10(MNewImagePar.fConc1) < (-0.371)*log10(MHillas.fSize) + 0.596 - -# This is the more complicated cut from ganymed.rc -#SparkCut.Condition: ({0} || {1}) && {2} -#SparkCut.0: MImagePar.fNumSatPixelsHG < 1 -#SparkCut.1: MHillas.GetArea*(MGeomCam.fConvMm2Deg^2) > (0.003*MImagePar.fNumSatPixelsHG) + 0.0325 -#SparkCut.2: log10(MNewImagePar.fConc1) < (-0.371)*log10(MHillas.fSize) + 0.596 - -# ------------------------------------------------------------------------- -# setup the image parameter calculation here -# ------------------------------------------------------------------------- -#MHillasCalc.IdxIsland: 0 - -# ------------------------------------------------------------------------- -# Setup or switch off the muon analysis here -# ------------------------------------------------------------------------- -#MJStar.MuonAnalysis: On - -#MMuonSetup.Margin: 0.2 -#MMuonSetup.ThresholdArcPhi: 30 -#MMuonSetup.ThresholdArcWidth: 2 - -#BinningRadius.Raw: 20 0.5 1.5 -#BinningArcWidth.Raw: 60 0.0 0.3 -#BinningRingBroadening.Raw: 20 0.5 1.5 -#BinningSizeVsArcRadius.Raw: 20 0.5 1.5 -#BinningsArcPhi.Raw: 20 -180 180 -#BinningsMuonWidth.Raw: 28 0.3 1.7 - -#MuonPreCut.Condition: MHillas.fSize>150 -#MuonSearcCut.Condition: (MMuonSearchPar.fRadius>180) && (MMuonSearchPar.fRadius<400) && (MMuonSearchPar.fDeviation<45) -#MuonFinalCut.Condition: (MMuonCalibPar.fArcPhi>190) && (MMuonSearchPar.fDeviation<35) && (MMuonCalibPar.fArcWidth<0.20) && (MMuonCalibPar.fArcWidth>0.04) Index: trunk/Mars/star_mux.rc =================================================================== --- trunk/Mars/star_mux.rc (revision 19871) +++ (revision ) @@ -1,101 +1,0 @@ -# ========================================================================== -############################################################################# -# ========================================================================== -# General -# ========================================================================== -############################################################################# -# ========================================================================== - -# ------------------------------------------------------------------------- -# Use this if you want to setup the logging stream for the jobs -# (overwrites command line options) -# ------------------------------------------------------------------------- -#MLog.VerbosityLevel: 2 -#MLog.DebugLevel: 1 -#MLog.NoColors: yes - -# ========================================================================== -############################################################################# -# ========================================================================== -# Star -# ========================================================================== -############################################################################# -# ========================================================================== - -# ------------------------------------------------------------------------- -# 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. -# ------------------------------------------------------------------------- -#MJStar.PathOut: . -#MJStar.PathIn: . - -# ------------------------------------------------------------------------- -# Configure MJStar -# ------------------------------------------------------------------------- -#MJStar.MaxEvents: 10000 -#MJStar.Overwrite: yes,no - -# ------------------------------------------------------------------------- -# Use this if you want to change the parameters of the software trigger -# ------------------------------------------------------------------------- -#FilterSwTrigger.Threshold: 5 -#FilterSwTrigger.NumNeighbors: 4 -#FilterSwTrigger.TimeWindow: 0.5 - -# ------------------------------------------------------------------------- -# It seems that the MUX FADCs have s lightly higher dead time... -# ------------------------------------------------------------------------- -MHEffectiveOnTime.FirstBin: 4 - -# ------------------------------------------------------------------------- -# Setup the image cleaning here -# ------------------------------------------------------------------------- -MImgCleanStd.CleanMethod: absolute -#MImgCleanStd.CleanLevel0: 6.0 -MImgCleanStd.CleanLevel1: 6.0 -MImgCleanStd.CleanLevel2: 3.0 -MImgCleanStd.TimeLevel1: 1.75 -MImgCleanStd.TimeLevel2: 1.75 -MImgCleanStd.PostCleanType: 3 - -#MImgCleanStd.CleanRings: 3 -#MImgCleanStd.KeepIsolatedPixels: Yes -#MImgCleanStd.RecoverIsolatedPixles: Yes - -# ------------------------------------------------------------------------- -# Setup the cut applied to the Sparkless tab -# ------------------------------------------------------------------------- - -# This is the simple default -#SparkCut.Condition: log10(MNewImagePar.fConc1) < (-0.371)*log10(MHillas.fSize) + 0.596 - -# This is the more complicated cut from ganymed.rc -#SparkCut.Condition: ({0} || {1}) && {2} -#SparkCut.0: MImagePar.fNumSatPixelsHG < 1 -#SparkCut.1: MHillas.GetArea*(MGeomCam.fConvMm2Deg^2) > (0.003*MImagePar.fNumSatPixelsHG) + 0.0325 -#SparkCut.2: log10(MNewImagePar.fConc1) < (-0.371)*log10(MHillas.fSize) + 0.596 - -# ------------------------------------------------------------------------- -# setup the image parameter calculation here -# ------------------------------------------------------------------------- -#MHillasCalc.IdxIsland: 0 - -# ------------------------------------------------------------------------- -# Setup or switch off the muon analysis here -# ------------------------------------------------------------------------- -#MJStar.MuonAnalysis: On - -#MMuonSetup.Margin: 0.2 -#MMuonSetup.ThresholdArcPhi: 30 -#MMuonSetup.ThresholdArcWidth: 2 - -#BinningRadius.Raw: 20 0.5 1.5 -#BinningArcWidth.Raw: 60 0.0 0.3 -#BinningRingBroadening.Raw: 20 0.5 1.5 -#BinningSizeVsArcRadius.Raw: 20 0.5 1.5 -#BinningsArcPhi.Raw: 20 -180 180 -#BinningsMuonWidth.Raw: 28 0.3 1.7 - -#MuonPreCut.Condition: MHillas.fSize>150 -#MuonSearcCut.Condition: (MMuonSearchPar.fRadius>180) && (MMuonSearchPar.fRadius<400) && (MMuonSearchPar.fDeviation<45) -#MuonFinalCut.Condition: (MMuonCalibPar.fArcPhi>190) && (MMuonSearchPar.fDeviation<35) && (MMuonCalibPar.fArcWidth<0.20) && (MMuonCalibPar.fArcWidth>0.04)