| 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 | # -------------------------------------------------------------------------
|
|---|
| 48 | #MPointingDevCalc.MaxAbsDev: 15
|
|---|
| 49 | #MPointingDevCalc.NumMinStars: 8
|
|---|
| 50 | #MPointingDevCalc.NsbLevel: 3.0
|
|---|
| 51 | #MPointingDevCalc.NsbMin: 30
|
|---|
| 52 | #MPointingDevCalc.NsbMax: 60
|
|---|
| 53 | #MPointingDevCalc.Dx -7
|
|---|
| 54 | #MPointingDevCalc.Dy 16
|
|---|
| 55 | #MPointingDevCalc.MaxAge 1.0
|
|---|
| 56 |
|
|---|
| 57 | # -------------------------------------------------------------------------
|
|---|
| 58 | # Setup the misfocussing correction
|
|---|
| 59 | # -------------------------------------------------------------------------
|
|---|
| 60 | #MSrcPosCorrect.Dx: 0
|
|---|
| 61 | #MSrcPosCorrect.Dy: 0
|
|---|
| 62 |
|
|---|
| 63 | # -------------------------------------------------------------------------
|
|---|
| 64 | # Setup fit mode and ranges
|
|---|
| 65 | # -------------------------------------------------------------------------
|
|---|
| 66 | #MAlphaFitter.SignalIntegralMax: 15
|
|---|
| 67 | #MAlphaFitter.SignalFitMax: 25
|
|---|
| 68 | #MAlphaFitter.BackgroundFitMin: 35
|
|---|
| 69 | #MAlphaFitter.BackgroundFitMax: 75
|
|---|
| 70 | #MAlphaFitter.ScaleMin: 35
|
|---|
| 71 | #MAlphaFitter.ScaleMax: 80
|
|---|
| 72 | #MAlphaFitter.PolynomOrder: 2
|
|---|
| 73 | #MAlphaFitter.ScaleMode: None
|
|---|
| 74 |
|
|---|
| 75 | # -------------------------------------------------------------------------
|
|---|
| 76 | # Define here which histogram class to use to determin the signal.
|
|---|
| 77 | # Currently availyble: MHAlpha, MHTheta <default>
|
|---|
| 78 | # -------------------------------------------------------------------------
|
|---|
| 79 | #MJCut.NameHist: MHAlpha
|
|---|
| 80 |
|
|---|
| 81 | # Be carefull: The spectrum might be done with a different binning
|
|---|
| 82 | # This might change the scale factor
|
|---|
| 83 | #MHThetaSq.NumBinsSignal: 3
|
|---|
| 84 | #MHThetaSq.NumBinsTotal: 75
|
|---|
| 85 |
|
|---|
| 86 | # -------------------------------------------------------------------------
|
|---|
| 87 | # To get a meaningfull false source plot in wobble mode use magic cuts
|
|---|
| 88 | # and the following two lines.
|
|---|
| 89 | # -------------------------------------------------------------------------
|
|---|
| 90 | MJCut.NameHistFS: MHDisp
|
|---|
| 91 |
|
|---|
| 92 | MHDisp.Smearing: 0.06
|
|---|
| 93 | #MHDisp.ScaleMin: 0.325
|
|---|
| 94 | #MHDisp.ScaleMax: 0.475
|
|---|
| 95 |
|
|---|
| 96 | # -------------------------------------------------------------------------
|
|---|
| 97 | # To get use a random source position in the off data with the distribution
|
|---|
| 98 | # in the camera of the on-data
|
|---|
| 99 | # -------------------------------------------------------------------------
|
|---|
| 100 | MJCut.RandomSourcePosition: Yes
|
|---|
| 101 |
|
|---|
| 102 | # -------------------------------------------------------------------------
|
|---|
| 103 | # There are three cuts:
|
|---|
| 104 | # CutQ: Cut after reading file
|
|---|
| 105 | # Cut0: After calculation of src dependent parameters, energy estimation
|
|---|
| 106 | # and disp calculation, before writing summary file
|
|---|
| 107 | # Cut1: After writing summary file before filling false source plots
|
|---|
| 108 | # Cut2: After filling false source plots before filling alpha plots
|
|---|
| 109 | # Cut3: After filling alpha plots (eg. Alpha cuts) for image parameter
|
|---|
| 110 | # display, eg. write only events after alpha cut to ganymed file
|
|---|
| 111 | # -------------------------------------------------------------------------
|
|---|
| 112 | CutQ.Inverted: Yes
|
|---|
| 113 | Cut0.Inverted: Yes
|
|---|
| 114 | Cut1.Inverted: Yes
|
|---|
| 115 | Cut2.Inverted: Yes
|
|---|
| 116 | Cut3.Inverted: Yes
|
|---|
| 117 |
|
|---|
| 118 | # Get rid of triangular events...
|
|---|
| 119 | CutQ.Condition: {0} && {1} && {2} && {3} && ({4} || {5})
|
|---|
| 120 | CutQ.0: MImagePar.fNumIslands<3
|
|---|
| 121 | CutQ.1: MNewImagePar.fNumUsedPixels>5
|
|---|
| 122 | CutQ.2: MNewImagePar.fLeakage1 < 0.3
|
|---|
| 123 | CutQ.3: MImagePar.fNumSatPixelsHG < 1
|
|---|
| 124 | CutQ.4: MNewImagePar.fConcCOG<0.25
|
|---|
| 125 | CutQ.5: log10(MNewImagePar.fConcCOG)<-0.3*(log10(MHillas.fSize)-2.5)-0.24
|
|---|
| 126 |
|
|---|
| 127 | # -------------------------------------------------------------------------
|
|---|
| 128 | # Such a cut could be used to exclude starnge pointing positions. like
|
|---|
| 129 | # the camera center in wobble mode. However this directly influences the
|
|---|
| 130 | # effective observation time, which must then be corrected manually!
|
|---|
| 131 | # It is currently not recommended to use such a cut except for tests!
|
|---|
| 132 | # -------------------------------------------------------------------------
|
|---|
| 133 | #Cut0.Condition: MSrcPosCam.GetDist*MGeomCam.fConvMm2Deg>0.3
|
|---|
| 134 |
|
|---|
| 135 |
|
|---|
| 136 | # ---------- SETUP FOR ONOFF-MODE -----------
|
|---|
| 137 | MAlphaFitter.ScaleMin: 0.137
|
|---|
| 138 | MAlphaFitter.ScaleMax: 0.640
|
|---|
| 139 | MAlphaFitter.BackgroundFitMin: 0.137
|
|---|
| 140 | MAlphaFitter.BackgroundFitMax: 1.000
|
|---|
| 141 | MAlphaFitter.PolynomOrder: 1
|
|---|
| 142 | MAlphaFitter.ScaleMode: OffRegion
|
|---|
| 143 | MAlphaFitter.SignalFunction: ThetaSq
|
|---|
| 144 |
|
|---|
| 145 | # For more details on the magic cuts see its class reference
|
|---|
| 146 | Cut1.Condition: MFMagicCuts
|
|---|
| 147 | Cut1.ThetaCut: None
|
|---|
| 148 | # Parametrization of Disp
|
|---|
| 149 | Cut1.Param0: 1.266195
|
|---|
| 150 | Cut1.Param8: 0.100577
|
|---|
| 151 | Cut1.Param9: 1.80309
|
|---|
| 152 | Cut1.Param10: 2.87177
|
|---|
| 153 | Cut1.Param11: 0.616823
|
|---|
| 154 | # Parametrization for sign of disp (m3long, slope)
|
|---|
| 155 | Cut1.Param5: -0.07
|
|---|
| 156 | Cut1.Param6: 7.2
|
|---|
| 157 | Cut1.Param7: 0.5
|
|---|
| 158 | # ThetaSq-Cut
|
|---|
| 159 | Cut1.Param1: 0.215
|
|---|
| 160 | # Area-Cut
|
|---|
| 161 | Cut1.Param2: 0.215468
|
|---|
| 162 | Cut1.Param3: 5.63973
|
|---|
| 163 | Cut1.Param4: 0.0836169
|
|---|
| 164 |
|
|---|
| 165 | # To switch on random forest uncomment the following lines
|
|---|
| 166 | #Cut1.HadronnessCut: Hadronness
|
|---|
| 167 | #Cut1.Param14: 0.5
|
|---|
| 168 | #Cut1.Param15: 100
|
|---|
| 169 |
|
|---|
| 170 | # To use both, hadronness and area cut do
|
|---|
| 171 | #Cut1.HadronnessCut: All
|
|---|
| 172 |
|
|---|
| 173 | # -------------------------------------------------------------------------
|
|---|
| 174 | # Use this to executa a task (eg to calc hadronness) before Cut1
|
|---|
| 175 | # -------------------------------------------------------------------------
|
|---|
| 176 | # CalcHadronness: MRanForestCalc
|
|---|
| 177 | # CalcHadronness.FileName: rf-separation.root
|
|---|
| 178 | # CalcHadronness.NameOutput: Hadronness
|
|---|
| 179 | # CalcHadronness.Debug: No
|
|---|
| 180 |
|
|---|
| 181 | # -------------------------------------------------------------------------
|
|---|
| 182 | # Use this to executa a task (eg to calc disp) before Cut1
|
|---|
| 183 | # -------------------------------------------------------------------------
|
|---|
| 184 | # CalcDisp: MRanForestCalc
|
|---|
| 185 | # CalcDisp.FileName: rf-disp.root
|
|---|
| 186 | # CalcDisp.NameOutput: Disp
|
|---|
| 187 | # CalcDisp.Debug: No
|
|---|
| 188 | # Cut1.CalcDisp: No
|
|---|
| 189 |
|
|---|
| 190 | # -------------------------------------------------------------------------
|
|---|
| 191 | # Use this to executa a task (eg to estimate energy) after all Cuts
|
|---|
| 192 | # -------------------------------------------------------------------------
|
|---|
| 193 | # EstimateEnergy: MRanForestCalc
|
|---|
| 194 | # EstimateEnergy.NameOutput: MEnergyEst
|
|---|
| 195 | # EstimateEnergy.FileName: ranforest.root
|
|---|
| 196 | # EstimateEnergy.Debug: No
|
|---|
