| 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.Dx -7
|
|---|
| 52 | #MPointingDevCalc.Dy 16
|
|---|
| 53 | #MPointingDevCalc.MaxAge 1.0
|
|---|
| 54 |
|
|---|
| 55 | # -------------------------------------------------------------------------
|
|---|
| 56 | # Setup the misfocussing correction
|
|---|
| 57 | # -------------------------------------------------------------------------
|
|---|
| 58 | #MSrcPosCorrect.Dx: 0
|
|---|
| 59 | #MSrcPosCorrect.Dy: 0
|
|---|
| 60 |
|
|---|
| 61 | # -------------------------------------------------------------------------
|
|---|
| 62 | # Setup fit mode and ranges
|
|---|
| 63 | # -------------------------------------------------------------------------
|
|---|
| 64 | #MAlphaFitter.SignalIntegralMax: 15
|
|---|
| 65 | #MAlphaFitter.SignalFitMax: 25
|
|---|
| 66 | #MAlphaFitter.BackgroundFitMin: 35
|
|---|
| 67 | #MAlphaFitter.BackgroundFitMax: 75
|
|---|
| 68 | #MAlphaFitter.ScaleMin: 35
|
|---|
| 69 | #MAlphaFitter.ScaleMax: 80
|
|---|
| 70 | #MAlphaFitter.PolynomOrder: 2
|
|---|
| 71 | #MAlphaFitter.ScaleMode: None
|
|---|
| 72 |
|
|---|
| 73 | # -------------------------------------------------------------------------
|
|---|
| 74 | # Define here which histogram class to use to determin the signal.
|
|---|
| 75 | # Currently availyble: MHAlpha, MHTheta <default>
|
|---|
| 76 | # -------------------------------------------------------------------------
|
|---|
| 77 | #MJCut.NameHist: MHAlpha
|
|---|
| 78 |
|
|---|
| 79 | # Be carefull: The spectrum might be done with a different binning
|
|---|
| 80 | # This might change the scale factor
|
|---|
| 81 | #MHThetaSq.NumBinsSignal: 3
|
|---|
| 82 | #MHThetaSq.NumBinsTotal: 75
|
|---|
| 83 |
|
|---|
| 84 | # -------------------------------------------------------------------------
|
|---|
| 85 | # To get a meaningfull false source plot in wobble mode use magic cuts
|
|---|
| 86 | # and the following two lines.
|
|---|
| 87 | # -------------------------------------------------------------------------
|
|---|
| 88 | MJCut.NameHistFS: MHDisp
|
|---|
| 89 |
|
|---|
| 90 | MHDisp.Smearing: 0.06
|
|---|
| 91 | #MHDisp.Wobble: Yes
|
|---|
| 92 | #MHDisp.ScaleMin: 0.325
|
|---|
| 93 | #MHDisp.ScaleMax: 0.475
|
|---|
| 94 |
|
|---|
| 95 | # -------------------------------------------------------------------------
|
|---|
| 96 | # To get use a random source position in the off data with the distribution
|
|---|
| 97 | # in the camera of the on-data
|
|---|
| 98 | # -------------------------------------------------------------------------
|
|---|
| 99 | #MJCut.RandomSourcePosition: Yes
|
|---|
| 100 |
|
|---|
| 101 | # -------------------------------------------------------------------------
|
|---|
| 102 | # There are three cuts:
|
|---|
| 103 | # CutQ: Cut after reading file
|
|---|
| 104 | # Cut0: After calculation of src dependant parameters, energy estimation
|
|---|
| 105 | # and disp calculation, before writing summary file
|
|---|
| 106 | # Cut1: After writing summary file before filling false source plots
|
|---|
| 107 | # Cut2: After filling false source plots before filling alpha plots
|
|---|
| 108 | # Cut3: After filling alpha plots (eg. Alpha cuts) for image parameter
|
|---|
| 109 | # display, eg. write only events after alpha cut to ganymed file
|
|---|
| 110 | # -------------------------------------------------------------------------
|
|---|
| 111 | CutQ.Inverted: Yes
|
|---|
| 112 | Cut0.Inverted: Yes
|
|---|
| 113 | Cut1.Inverted: Yes
|
|---|
| 114 | Cut2.Inverted: Yes
|
|---|
| 115 | Cut3.Inverted: Yes
|
|---|
| 116 |
|
|---|
| 117 | # Get rid of triangular events...
|
|---|
| 118 | CutQ.Condition: ({0} || {1}) && {2} && {3} && {4} && {5} && {6}
|
|---|
| 119 | CutQ.0: MImagePar.fNumSatPixelsHG < 1
|
|---|
| 120 | CutQ.1: MHillas.GetArea*(MGeomCam.fConvMm2Deg^2) > (0.003*MImagePar.fNumSatPixelsHG) + 0.0325
|
|---|
| 121 | CutQ.2: MNewImagePar.fNumUsedPixels>5
|
|---|
| 122 | CutQ.3: MNewImagePar.fLeakage1 < 0.3
|
|---|
| 123 | CutQ.4: MImagePar.fNumSatPixelsLG<2
|
|---|
| 124 | CutQ.5: MImagePar.fNumIslands<3
|
|---|
| 125 | CutQ.6: log10(MNewImagePar.fConc1) < (-0.371)*log10(MHillas.fSize) + 0.596
|
|---|
| 126 |
|
|---|
| 127 | # ---------- SETUP FOR ONOFF-MODE -----------
|
|---|
| 128 | MAlphaFitter.ScaleMin: 0.137
|
|---|
| 129 | MAlphaFitter.ScaleMax: 0.640
|
|---|
| 130 | MAlphaFitter.BackgroundFitMin: 0.137
|
|---|
| 131 | MAlphaFitter.BackgroundFitMax: 1.000
|
|---|
| 132 | MAlphaFitter.PolynomOrder: 1
|
|---|
| 133 | MAlphaFitter.ScaleMode: OffRegion
|
|---|
| 134 | MAlphaFitter.SignalFunction: ThetaSq
|
|---|
| 135 |
|
|---|
| 136 | # For more details on the magic cuts see its class reference
|
|---|
| 137 | Cut1.Condition: MFMagicCuts
|
|---|
| 138 | Cut1.ThetaCut: None
|
|---|
| 139 | # Parametrization of Disp
|
|---|
| 140 | Cut1.Param0: 1.3245
|
|---|
| 141 | Cut1.Param6: 8.2957
|
|---|
| 142 | Cut1.Param7: 0.8677
|
|---|
| 143 | # Parametrization for sign of disp (m3long)
|
|---|
| 144 | Cut1.Param5: -0.07
|
|---|
| 145 | # Parametrization for sign of disp (slope)
|
|---|
| 146 | Cut1.Param8: 7.2
|
|---|
| 147 | Cut1.Param9: 0.5
|
|---|
| 148 | # ThetaSq-Cut
|
|---|
| 149 | Cut1.Param1: 0.22
|
|---|
| 150 | # Area-Cut
|
|---|
| 151 | Cut1.Param2: 0.215
|
|---|
| 152 | Cut1.Param3: 5.486
|
|---|
| 153 | Cut1.Param4: 0.0897
|
|---|
| 154 |
|
|---|
| 155 | # ---------- SETUP FOR WOBBLE-MODE -----------
|
|---|
| 156 | #MAlphaFitter.ScaleMode: none
|
|---|
| 157 | #MAlphaFitter.BackgroundFitMin: 0.09
|
|---|
| 158 | #MAlphaFitter.BackgroundFitMax: 0.25
|
|---|
| 159 | #MAlphaFitter.PolynomOrder: 1
|
|---|
| 160 | #MAlphaFitter.SignalFunction: ThetaSq
|
|---|
| 161 |
|
|---|
| 162 | # For more details on the magic cuts see its class reference
|
|---|
| 163 | #Cut1.Condition: MFMagicCuts
|
|---|
| 164 | #Cut1.ThetaCut: Off
|
|---|
| 165 | # Parametrization of Disp
|
|---|
| 166 | #Cut1.Param0: 1.221
|
|---|
| 167 | #Cut1.Param6: 13.425
|
|---|
| 168 | #Cut1.Param7: 1.04343
|
|---|
| 169 | # Parametrization for sign of disp (m3long)
|
|---|
| 170 | #Cut1.Param5: -0.07
|
|---|
| 171 | # Parametrization for sign of disp (slope)
|
|---|
| 172 | #Cut1.Param8: 7.2
|
|---|
| 173 | #Cut1.Param9: 0.5
|
|---|
| 174 | # ThetaSq-Cut
|
|---|
| 175 | #Cut1.Param1: 0.22
|
|---|
| 176 | # Area-Cut
|
|---|
| 177 | #Cut1.Param2: 0.222
|
|---|
| 178 | #Cut1.Param3: 5.524
|
|---|
| 179 | #Cut1.Param4: 0.0878
|
|---|
| 180 |
|
|---|
| 181 | # ---------- SETUP TO USE RANDOM FOREST -----------
|
|---|
| 182 |
|
|---|
| 183 | # To switch on random forest uncomment the following lines
|
|---|
| 184 | #Cut1.HadronnessCut: Hadronness
|
|---|
| 185 | #Cut1.Param10: 0.5
|
|---|
| 186 | #Cut1.Param11: 100
|
|---|
| 187 |
|
|---|
| 188 | # To use both, hadronness and area cut do
|
|---|
| 189 | #Cut1.HadronnessCut: All
|
|---|
| 190 |
|
|---|
| 191 |
|
|---|
| 192 | # -------------------------------------------------------------------------
|
|---|
| 193 | # Use this to executa a task (eg to calc hadronness) before Cut1
|
|---|
| 194 | # -------------------------------------------------------------------------
|
|---|
| 195 | # CalcHadronness: MRanForestCalc
|
|---|
| 196 | # CalcHadronness.FileName: rf-separation.root
|
|---|
| 197 | # CalcHadronness.NameOutput: Hadronness
|
|---|
| 198 | # CalcHadronness.Debug: No
|
|---|
| 199 |
|
|---|
| 200 | # -------------------------------------------------------------------------
|
|---|
| 201 | # Use this to executa a task (eg to calc disp) before Cut1
|
|---|
| 202 | # -------------------------------------------------------------------------
|
|---|
| 203 | # CalcDisp: MRanForestCalc
|
|---|
| 204 | # CalcDisp.FileName: rf-disp.root
|
|---|
| 205 | # CalcDisp.NameOutput: Disp
|
|---|
| 206 | # CalcDisp.Debug: No
|
|---|
| 207 | # Cut1.CalcDisp: No
|
|---|
| 208 |
|
|---|
| 209 | # -------------------------------------------------------------------------
|
|---|
| 210 | # Use this to executa a task (eg to estimate energy) after all Cuts
|
|---|
| 211 | # -------------------------------------------------------------------------
|
|---|
| 212 | # EstimateEnergy: MRanForestCalc
|
|---|
| 213 | # EstimateEnergy.NameOutput: MEnergyEst
|
|---|
| 214 | # EstimateEnergy.FileName: ranforest.root
|
|---|
| 215 | # EstimateEnergy.Debug: No
|
|---|
