# ==========================================================================
#############################################################################
# ==========================================================================
#                              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.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:
#   Cut0: After energy estimation, 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
# -------------------------------------------------------------------------
Cut0.Inverted: Yes
Cut1.Inverted: Yes
Cut2.Inverted: Yes
Cut3.Inverted: Yes

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

# ---------- 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
# ThetaSq-Cut
Cut1.Param1:  0.22
# Area-Cut
Cut1.Param2:  0.222
Cut1.Param3:  5.524
Cut1.Param4:  0.0878

# To switch on random forest uncomment the following lines
#Cut1.HadronnessCut: Hadronness
#Cut1.Param8: 0.5
#Cut1.Param9: 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