| 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 | BinningWidth.Raw: 50 0 1
|
|---|
| 37 | BinningLength.Raw: 25 0 1
|
|---|
| 38 | BinningDist.Raw: 25 0 2
|
|---|
| 39 | BinningAsym.Raw: 51 -0.7 0.7
|
|---|
| 40 | BinningM3Long.Raw: 51 -0.7 0.7
|
|---|
| 41 | BinningM3Trans.Raw: 51 -0.7 0.7
|
|---|
| 42 | BinningSlope.Raw: 51 -7.5 7.5
|
|---|
| 43 |
|
|---|
| 44 | # -------------------------------------------------------------------------
|
|---|
| 45 | # Using the starguider for pointing correction.
|
|---|
| 46 | # To switch it off use "MPointingDevCalc.MaxAbsDev: -1"
|
|---|
| 47 | # For details see the class reference of MPointingDevCalc
|
|---|
| 48 | # -------------------------------------------------------------------------
|
|---|
| 49 | #MPointingDevCalc.MaxAbsDev: 15
|
|---|
| 50 | #MPointingDevCalc.NumMinStars: 8
|
|---|
| 51 | #MPointingDevCalc.NsbLevel: 3.0
|
|---|
| 52 | #MPointingDevCalc.NsbMin: 30
|
|---|
| 53 | #MPointingDevCalc.NsbMax: 60
|
|---|
| 54 | #MPointingDevCalc.Dx -0.001
|
|---|
| 55 | #MPointingDevCalc.Dy -0.004
|
|---|
| 56 | #MPointingDevCalc.MaxAge 1.0
|
|---|
| 57 |
|
|---|
| 58 | MPointingDevCalc.PointingModels: 1 81855 83722 89180 267253 291104
|
|---|
| 59 |
|
|---|
| 60 | # -------------------------------------------------------------------------
|
|---|
| 61 | # Setup misfocussing correction (see class reference for MSrcPosCorrect)
|
|---|
| 62 | # -------------------------------------------------------------------------
|
|---|
| 63 | #MSrcPosCorrect.Dx: 0
|
|---|
| 64 | #MSrcPosCorrect.Dy: 0
|
|---|
| 65 |
|
|---|
| 66 | # -------------------------------------------------------------------------
|
|---|
| 67 | # Setup fit mode and ranges
|
|---|
| 68 | # -------------------------------------------------------------------------
|
|---|
| 69 | #MAlphaFitter.SignalIntegralMax: 15
|
|---|
| 70 | #MAlphaFitter.SignalFitMax: 25
|
|---|
| 71 | #MAlphaFitter.BackgroundFitMin: 35
|
|---|
| 72 | #MAlphaFitter.BackgroundFitMax: 75
|
|---|
| 73 | #MAlphaFitter.ScaleMin: 35
|
|---|
| 74 | #MAlphaFitter.ScaleMax: 80
|
|---|
| 75 | #MAlphaFitter.PolynomOrder: 2
|
|---|
| 76 | #MAlphaFitter.ScaleMode: None
|
|---|
| 77 |
|
|---|
| 78 | # -------------------------------------------------------------------------
|
|---|
| 79 | # Define here which histogram class to use to determin the signal.
|
|---|
| 80 | # Currently availyble: MHAlpha, MHTheta <default>
|
|---|
| 81 | # -------------------------------------------------------------------------
|
|---|
| 82 | #MJCut.NameHist: MHAlpha
|
|---|
| 83 |
|
|---|
| 84 | # Be carefull: The spectrum might be done with a different binning
|
|---|
| 85 | # This might change the scale factor
|
|---|
| 86 | #MHThetaSq.NumBinsSignal: 3
|
|---|
| 87 | #MHThetaSq.NumBinsTotal: 75
|
|---|
| 88 |
|
|---|
| 89 | # -------------------------------------------------------------------------
|
|---|
| 90 | # To get a meaningfull false source plot in wobble mode use magic cuts
|
|---|
| 91 | # and the following two lines.
|
|---|
| 92 | # -------------------------------------------------------------------------
|
|---|
| 93 | MJCut.NameHistFS: MHDisp
|
|---|
| 94 |
|
|---|
| 95 | MHDisp.Smearing: 0.06
|
|---|
| 96 | #MHDisp.Wobble: Yes
|
|---|
| 97 | #MHDisp.ScaleMin: 0.325
|
|---|
| 98 | #MHDisp.ScaleMax: 0.475
|
|---|
| 99 |
|
|---|
| 100 | # -------------------------------------------------------------------------
|
|---|
| 101 | # To get/use a random source position in the off data with the distribution
|
|---|
| 102 | # in the camera of the on-data
|
|---|
| 103 | # -------------------------------------------------------------------------
|
|---|
| 104 | #MJCut.RandomSourcePosition: Yes
|
|---|
| 105 |
|
|---|
| 106 | # -------------------------------------------------------------------------
|
|---|
| 107 | # There are three cuts:
|
|---|
| 108 | # CutQ: Cut after reading file
|
|---|
| 109 | # Cut0: After calculation of src dependent parameters, energy estimation
|
|---|
| 110 | # and disp calculation, before writing summary file
|
|---|
| 111 | # Cut1: After writing summary file before filling false source plots
|
|---|
| 112 | # Cut2: After filling false source plots before filling alpha plots
|
|---|
| 113 | # Cut3: After filling alpha plots (eg. Alpha cuts) for image parameter
|
|---|
| 114 | # display, eg. write only events after alpha cut to ganymed file
|
|---|
| 115 | # -------------------------------------------------------------------------
|
|---|
| 116 | CutQ.Inverted: Yes
|
|---|
| 117 | Cut0.Inverted: Yes
|
|---|
| 118 | Cut1.Inverted: Yes
|
|---|
| 119 | Cut2.Inverted: Yes
|
|---|
| 120 | Cut3.Inverted: Yes
|
|---|
| 121 |
|
|---|
| 122 | # Quality-cuts (remove unphysical and obvious non-gamma events)
|
|---|
| 123 | CutQ.Condition: {0} && {1} && {2} && {3} && {4} && {5}
|
|---|
| 124 | CutQ.0: MImagePar.fNumIslands<3
|
|---|
| 125 | CutQ.1: MNewImagePar.fNumUsedPixels>5
|
|---|
| 126 | CutQ.2: MNewImagePar.fLeakage1 < 0.3
|
|---|
| 127 | CutQ.3: log10(MNewImagePar.fConcCOG)<-0.45+0.08*(log10(MHillas.fSize)<3.9)*(log10(MHillas.fSize)-3.9)^2
|
|---|
| 128 | CutQ.4: log10(MNewImagePar.fConc1) <-0.75+0.10*(log10(MHillas.fSize)<3.8)*(log10(MHillas.fSize)-3.8)^2
|
|---|
| 129 | CutQ.5: MHillas.fLength>-5*(log10(MHillas.fSize)-4.7)^2+50
|
|---|
| 130 |
|
|---|
| 131 | # -------------------------------------------------------------------------
|
|---|
| 132 | # Such a cut could be used to exclude strange pointing positions. like
|
|---|
| 133 | # the camera center in wobble mode. Note, that a cut in the poiting
|
|---|
| 134 | # position directly influences the effective observation time and should
|
|---|
| 135 | # only be applied if the events are lost anyhow. In both cases you must
|
|---|
| 136 | # correct for the loss of observation time calculating your flux!
|
|---|
| 137 | # -------------------------------------------------------------------------
|
|---|
| 138 | #Cut0.Condition: MSrcPosCam.GetDist*MGeomCam.fConvMm2Deg>0.3
|
|---|
| 139 |
|
|---|
| 140 | # ---------- SETUP FOR ONOFF-MODE -----------
|
|---|
| 141 | MAlphaFitter.ScaleMin: 0.137
|
|---|
| 142 | MAlphaFitter.ScaleMax: 0.640
|
|---|
| 143 | MAlphaFitter.BackgroundFitMin: 0.137
|
|---|
| 144 | MAlphaFitter.BackgroundFitMax: 1.000
|
|---|
| 145 | MAlphaFitter.PolynomOrder: 1
|
|---|
| 146 | MAlphaFitter.ScaleMode: OffRegion
|
|---|
| 147 | MAlphaFitter.SignalFunction: ThetaSq
|
|---|
| 148 |
|
|---|
| 149 | # For more details on the magic cuts see its class reference
|
|---|
| 150 | Cut1.Condition: MFMagicCuts
|
|---|
| 151 | Cut1.ThetaCut: None
|
|---|
| 152 | # Parametrization of Disp
|
|---|
| 153 | Cut1.Param0: 1.15136
|
|---|
| 154 | Cut1.Param8: 0.0681437
|
|---|
| 155 | Cut1.Param9: 2.62932
|
|---|
| 156 | Cut1.Param10: 1.51279
|
|---|
| 157 | Cut1.Param11: 0.0507821
|
|---|
| 158 | # Parametrization for sign of disp (m3long, slope)
|
|---|
| 159 | Cut1.Param5: -0.07
|
|---|
| 160 | Cut1.Param6: 7.2
|
|---|
| 161 | Cut1.Param7: 0.5
|
|---|
| 162 | # ThetaSq-Cut
|
|---|
| 163 | Cut1.Param1: 0.215
|
|---|
| 164 | # Area-Cut
|
|---|
| 165 | Cut1.Param2: 0.215468
|
|---|
| 166 | Cut1.Param3: 5.63973
|
|---|
| 167 | Cut1.Param4: 0.0836169
|
|---|
| 168 |
|
|---|
| 169 | # ---------- SETUP FOR WOBBLE-MODE -----------
|
|---|
| 170 | #MAlphaFitter.ScaleMode: none
|
|---|
| 171 | #MAlphaFitter.BackgroundFitMin: 0.09
|
|---|
| 172 | #MAlphaFitter.BackgroundFitMax: 0.25
|
|---|
| 173 | #MAlphaFitter.PolynomOrder: 1
|
|---|
| 174 | #MAlphaFitter.SignalFunction: ThetaSq
|
|---|
| 175 |
|
|---|
| 176 | # ---------- SETUP TO USE RANDOM FOREST -----------
|
|---|
| 177 |
|
|---|
| 178 | # To switch on random forest uncomment the following lines
|
|---|
| 179 | #Cut1.HadronnessCut: Hadronness
|
|---|
| 180 | #Cut1.Param14: 0.5
|
|---|
| 181 | #Cut1.Param15: 100
|
|---|
| 182 |
|
|---|
| 183 | # To use both, hadronness and area cut do
|
|---|
| 184 | #Cut1.HadronnessCut: All
|
|---|
| 185 |
|
|---|
| 186 | # -------------------------------------------------------------------------
|
|---|
| 187 | # This cut is applied exclusively to the plots in the tab CutT. It is
|
|---|
| 188 | # resonable to show the events within your signal region here to
|
|---|
| 189 | # see if your signal region is clean from strange events (e.g. sparks)
|
|---|
| 190 | # -------------------------------------------------------------------------
|
|---|
| 191 | CutT.Condition: ThetaSquared.fVal < ThetaSquaredCut.fVal
|
|---|
| 192 |
|
|---|
| 193 | # -------------------------------------------------------------------------
|
|---|
| 194 | # Use this to executa a task (eg to calc hadronness) before Cut1
|
|---|
| 195 | # -------------------------------------------------------------------------
|
|---|
| 196 | # CalcHadronness: MRanForestCalc
|
|---|
| 197 | # CalcHadronness.FileName: rf-separation.root
|
|---|
| 198 | # CalcHadronness.NameOutput: Hadronness
|
|---|
| 199 | # CalcHadronness.Debug: No
|
|---|
| 200 |
|
|---|
| 201 | # -------------------------------------------------------------------------
|
|---|
| 202 | # Use this to executa a task (eg to calc disp) before Cut1
|
|---|
| 203 | # -------------------------------------------------------------------------
|
|---|
| 204 | # CalcDisp: MRanForestCalc
|
|---|
| 205 | # CalcDisp.FileName: rf-disp.root
|
|---|
| 206 | # CalcDisp.NameOutput: Disp
|
|---|
| 207 | # CalcDisp.Debug: No
|
|---|
| 208 | # Cut1.CalcDisp: No
|
|---|
| 209 |
|
|---|
| 210 | # CalcDisp: MParameterCalc
|
|---|
| 211 | # CalcDisp.Rule: 1.15*(1-MHillas.fWidth/MHillas.fLength)
|
|---|
| 212 | # CalcDisp.NameOutput: Disp
|
|---|
| 213 | # Cut1.CalcDisp: No
|
|---|
| 214 |
|
|---|
| 215 | # -------------------------------------------------------------------------
|
|---|
| 216 | # Use this to executa a task (eg to estimate energy) after all Cuts
|
|---|
| 217 | # -------------------------------------------------------------------------
|
|---|
| 218 | # EstimateEnergy: MRanForestCalc
|
|---|
| 219 | # EstimateEnergy.NameOutput: MEnergyEst
|
|---|
| 220 | # EstimateEnergy.FileName: ranforest.root
|
|---|
| 221 | # EstimateEnergy.Debug: No
|
|---|
