# ==========================================================================
#############################################################################
# ==========================================================================
#                              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
#BinningTheta.Raw:       50    0    60      cos
BinningWidth.Raw:        50    0    1
BinningLength.Raw:       25    0    1
BinningDist.Raw:         25    0    2
BinningMaxDist.Raw:      25    0    2

# -------------------------------------------------------------------------
# Using the starguider for pointing correction. 
# To switch it off use "MPointingDevCalc.MaxAbsDev: -1"
# -------------------------------------------------------------------------
#MPointingDevCalc.MaxAbsDev:   15
#MPointingDevCalc.NumMinStars: 8
#MPointingDevCalc.NsbLevel:    3.0
#MPointingDevCalc.NsbMin:      30
#MPointingDevCalc.NsbMax:      60
#MPointingDevCalc.Dx           -7
#MPointingDevCalc.Dy           16
#MPointingDevCalc.MaxAge       1.0

# -------------------------------------------------------------------------
# Setup the misfocussing correction
# -------------------------------------------------------------------------
#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 <default>
# -------------------------------------------------------------------------
#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 dependant 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

# Get rid of triangular events...
CutQ.Condition: ({0} || {1}) && {2} && {3} && {4} && {5} && {6}
CutQ.0: MImagePar.fNumSatPixelsHG < 1
CutQ.1: MHillas.GetArea*(MGeomCam.fConvMm2Deg^2) > (0.003*MImagePar.fNumSatPixelsHG) + 0.0325
CutQ.2: MNewImagePar.fNumUsedPixels>5
CutQ.3: MNewImagePar.fLeakage1 < 0.3
CutQ.4: MImagePar.fNumSatPixelsLG<2
CutQ.5: MImagePar.fNumIslands<3
CutQ.6: log10(MNewImagePar.fConc1) < (-0.371)*log10(MHillas.fSize) + 0.596

# ---------- 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.3245
Cut1.Param6:  8.2957
Cut1.Param7:  0.8677
# Parametrization for sign of disp (m3long)
Cut1.Param5: -0.07
# Parametrization for sign of disp (slope)
Cut1.Param8:  7.2
Cut1.Param9:  0.5
# ThetaSq-Cut
Cut1.Param1:  0.22
# Area-Cut
Cut1.Param2:  0.215
Cut1.Param3:  5.486
Cut1.Param4:  0.0897

# ---------- 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: Off
# Parametrization of Disp
#Cut1.Param0:  1.221
#Cut1.Param6:  13.425
#Cut1.Param7:  1.04343
# Parametrization for sign of disp (m3long)
#Cut1.Param5: -0.07
# Parametrization for sign of disp (slope)
#Cut1.Param8:  7.2
#Cut1.Param9:  0.5
# ThetaSq-Cut
#Cut1.Param1:  0.22
# Area-Cut
#Cut1.Param2:  0.222
#Cut1.Param3:  5.524
#Cut1.Param4:  0.0878

# ---------- SETUP TO USE RANDOM FOREST -----------

# To switch on random forest uncomment the following lines
#Cut1.HadronnessCut: Hadronness
#Cut1.Param10: 0.5
#Cut1.Param11: 100

# To use both, hadronness and area cut do
#Cut1.HadronnessCut: All


# -------------------------------------------------------------------------
# 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

# -------------------------------------------------------------------------
# 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