| 1 | # ==========================================================================
|
|---|
| 2 | #############################################################################
|
|---|
| 3 | # ==========================================================================
|
|---|
| 4 | # General
|
|---|
| 5 | # ==========================================================================
|
|---|
| 6 | #############################################################################
|
|---|
| 7 | # ==========================================================================
|
|---|
| 8 |
|
|---|
| 9 | # -------------------------------------------------------------------------
|
|---|
| 10 | # Use this if you want to setup the logging stream for the jobs
|
|---|
| 11 | # (overwrites command line options)
|
|---|
| 12 | # -------------------------------------------------------------------------
|
|---|
| 13 | #MLog.VerbosityLevel: 2
|
|---|
| 14 | #MLog.DebugLevel: 1
|
|---|
| 15 | #MLog.NoColors: yes
|
|---|
| 16 |
|
|---|
| 17 | # ==========================================================================
|
|---|
| 18 | #############################################################################
|
|---|
| 19 | # ==========================================================================
|
|---|
| 20 | # Cuts
|
|---|
| 21 | # ==========================================================================
|
|---|
| 22 | #############################################################################
|
|---|
| 23 | # ==========================================================================
|
|---|
| 24 |
|
|---|
| 25 | # -------------------------------------------------------------------------
|
|---|
| 26 | # Use this if you want to write output to somewhere here you can give
|
|---|
| 27 | # the output path
|
|---|
| 28 | # -------------------------------------------------------------------------
|
|---|
| 29 | #MJCuts.PathOut: .
|
|---|
| 30 |
|
|---|
| 31 | # -------------------------------------------------------------------------
|
|---|
| 32 | # Use this to setup binnings. For more details see: MBinning::ReadEnv
|
|---|
| 33 | # -------------------------------------------------------------------------
|
|---|
| 34 | BinningAlpha.Raw: 36 0 90
|
|---|
| 35 | BinningFalseSource.Raw: 91 -1.5 1.5
|
|---|
| 36 | #BinningTheta.Raw: 50 0 60 cos
|
|---|
| 37 | BinningWidth.Raw: 50 0 1
|
|---|
| 38 | BinningLength.Raw: 25 0 1
|
|---|
| 39 | BinningDist.Raw: 25 0 2
|
|---|
| 40 | BinningMaxDist.Raw: 25 0 2
|
|---|
| 41 |
|
|---|
| 42 | # -------------------------------------------------------------------------
|
|---|
| 43 | # Using the starguider for pointing correction.
|
|---|
| 44 | # To switch it off use "MPointingDevCalc.MaxAbsDev: -1"
|
|---|
| 45 | # -------------------------------------------------------------------------
|
|---|
| 46 | #MPointingDevCalc.MaxAbsDev: 15
|
|---|
| 47 | #MPointingDevCalc.NumMinStars: 8
|
|---|
| 48 | #MPointingDevCalc.NsbLevel: 3.0
|
|---|
| 49 | #MPointingDevCalc.NsbMin: 30
|
|---|
| 50 | #MPointingDevCalc.NsbMax: 60
|
|---|
| 51 | #MPointingDevCalc.MaxAbsDev: 15
|
|---|
| 52 | #MPointingDevCalc.Dx -7
|
|---|
| 53 | #MPointingDevCalc.Dy 16
|
|---|
| 54 |
|
|---|
| 55 | # -------------------------------------------------------------------------
|
|---|
| 56 | # Setup fit mode and ranges
|
|---|
| 57 | # -------------------------------------------------------------------------
|
|---|
| 58 | #MAlphaFitter.SignalIntegralMax: 15
|
|---|
| 59 | #MAlphaFitter.SignalFitMax: 25
|
|---|
| 60 | #MAlphaFitter.BackgroundFitMin: 35
|
|---|
| 61 | #MAlphaFitter.BackgroundFitMax: 75
|
|---|
| 62 | #MAlphaFitter.ScaleMin: 35
|
|---|
| 63 | #MAlphaFitter.ScaleMax: 80
|
|---|
| 64 | #MAlphaFitter.PolynomOrder: 2
|
|---|
| 65 | #MAlphaFitter.ScaleMode: None
|
|---|
| 66 |
|
|---|
| 67 | # -------------------------------------------------------------------------
|
|---|
| 68 | # Define here which histogram class to use to determin the signal.
|
|---|
| 69 | # Currently availyble: MHAlpha, MHTheta <default>
|
|---|
| 70 | # -------------------------------------------------------------------------
|
|---|
| 71 | #MJCut.NameHist: MHAlpha
|
|---|
| 72 |
|
|---|
| 73 | # -------------------------------------------------------------------------
|
|---|
| 74 | # To get a meaningfull false source plot in wobble mode use magic cuts
|
|---|
| 75 | # and the following two lines.
|
|---|
| 76 | # -------------------------------------------------------------------------
|
|---|
| 77 | MJCut.NameHistFS: MHDisp
|
|---|
| 78 |
|
|---|
| 79 | MHDisp.Wobble: Yes
|
|---|
| 80 | MHDisp.Smearing: 0.06
|
|---|
| 81 | #MHDisp.ScaleMin: 0.325
|
|---|
| 82 | #MHDisp.ScaleMax: 0.475
|
|---|
| 83 |
|
|---|
| 84 | # -------------------------------------------------------------------------
|
|---|
| 85 | # Use this to set the number of bins corresponding to the signal region
|
|---|
| 86 | # and if you use a Anti-Theta cut the second one has to be enabled
|
|---|
| 87 | # -------------------------------------------------------------------------
|
|---|
| 88 | #MHPhi.NumBinsSignal: 2
|
|---|
| 89 | #MHPhi.UseAntiPhiCut: No
|
|---|
| 90 |
|
|---|
| 91 | # -------------------------------------------------------------------------
|
|---|
| 92 | # There are three cuts:
|
|---|
| 93 | # Cut0: After energy estimation, before writing summary file
|
|---|
| 94 | # Cut1: After writing summary file before filling false source plots
|
|---|
| 95 | # Cut2: After filling false source plots before filling alpha plots
|
|---|
| 96 | # Cut3: After filling alpha plots (eg. Alpha cuts) for image parameter
|
|---|
| 97 | # display
|
|---|
| 98 | # -------------------------------------------------------------------------
|
|---|
| 99 | Cut0.Inverted: Yes
|
|---|
| 100 | Cut1.Inverted: Yes
|
|---|
| 101 | Cut2.Inverted: Yes
|
|---|
| 102 | Cut3.Inverted: Yes
|
|---|
| 103 |
|
|---|
| 104 | # Get rid of triangular events...
|
|---|
| 105 | Cut0.Condition: ({0} || {1}) && {2} && {3} && {4} && {5} && {6}
|
|---|
| 106 | Cut0.0: MImagePar.fNumSatPixelsHG < 1
|
|---|
| 107 | Cut0.1: MHillas.GetArea*(MGeomCam.fConvMm2Deg^2) > (0.002*MImagePar.fNumSatPixelsHG) + 0.025
|
|---|
| 108 | Cut0.2: MNewImagePar.fNumUsedPixels>5
|
|---|
| 109 | Cut0.3: MNewImagePar.fLeakage1 < 0.3
|
|---|
| 110 | Cut0.4: MImagePar.fNumSatPixelsLG<2
|
|---|
| 111 | Cut0.5: MImagePar.fNumIslands<3
|
|---|
| 112 | Cut0.6: log10(MNewImagePar.fConc1) < (-0.371)*log10(MHillas.fSize) + 0.596
|
|---|
| 113 |
|
|---|
| 114 | # ---------- SETUP FOR WOBBLE-MODE -----------
|
|---|
| 115 | MAlphaFitter.ScaleMode: none
|
|---|
| 116 | MAlphaFitter.BackgroundFitMin: 0.09
|
|---|
| 117 | MAlphaFitter.BackgroundFitMax: 0.25
|
|---|
| 118 | MAlphaFitter.PolynomOrder: 1
|
|---|
| 119 | MAlphaFitter.SignalFunction: ThetaSq
|
|---|
| 120 |
|
|---|
| 121 | Cut1.Condition: MFMagicCuts
|
|---|
| 122 | Cut1.ThetaCut: Off
|
|---|
| 123 | Cut1.Param0: 1.221
|
|---|
| 124 | Cut1.Param1: 0.22
|
|---|
| 125 | Cut1.Param2: 0.222
|
|---|
| 126 | Cut1.Param3: 5.524
|
|---|
| 127 | Cut1.Param4: 0.0878
|
|---|
| 128 | Cut1.Param5: -0.07
|
|---|
| 129 | Cut1.Param6: 13.425
|
|---|
| 130 | Cut1.Param7: 1.04343
|
|---|
| 131 |
|
|---|
| 132 | # Be carefull: The spectrum might be done with a different binning
|
|---|
| 133 | # This might change the scale factor
|
|---|
| 134 | #MHThetaSq.NumBinsSignal: 3
|
|---|
| 135 | #MHThetaSq.NumBinsTotal: 75
|
|---|
| 136 |
|
|---|
| 137 | # -------------------------------------------------------------------------
|
|---|
| 138 | # Setup cuts in false source plot
|
|---|
| 139 | # -------------------------------------------------------------------------
|
|---|
| 140 | #MHFalseSource.DistMin: 0.55
|
|---|
| 141 | #MHFalseSource.DistMax: 0.55
|
|---|
| 142 | #MHFalseSource.DWMin: 0.55
|
|---|
| 143 | #MHFalseSource.DWMax: 0.55
|
|---|
| 144 |
|
|---|
| 145 | # -------------------------------------------------------------------------
|
|---|
| 146 | # Use this to executa a task (eg to calc hadronness) before Cut1
|
|---|
| 147 | # -------------------------------------------------------------------------
|
|---|
| 148 | # CalcHadronness: MRanForestCalc
|
|---|
| 149 | # CalcHadronness.FileName: rf-separation.root
|
|---|
| 150 | # CalcHadronness.NameOutput: Hadronness
|
|---|
| 151 | # CalcHadronness.Debug: No
|
|---|
| 152 |
|
|---|
| 153 | # -------------------------------------------------------------------------
|
|---|
| 154 | # Use this to executa a task (eg to calc disp) before Cut1
|
|---|
| 155 | # -------------------------------------------------------------------------
|
|---|
| 156 | # CalcDisp: MRanForestCalc
|
|---|
| 157 | # CalcDisp.FileName: rf-disp.root
|
|---|
| 158 | # CalcDisp.NameOutput: Disp
|
|---|
| 159 | # CalcDisp.Debug: No
|
|---|
| 160 | # Cut1.CalcDisp: No
|
|---|
