- Timestamp:
- 08/24/07 13:44:23 (17 years ago)
- Location:
- trunk/MagicSoft/Mars
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/ganymed.rc
r8668 r8705 128 128 129 129 # ------------------------------------------------------------------------- 130 # Such a cut could be used to exclude starnge pointing positions. like 131 # the camera center in wobble mode. However this directly influences the 132 # effective observation time, which must then be corrected manually! 133 # It is currently not recommended to use such a cut except for tests! 130 # Such a cut could be used to exclude strange pointing positions. like 131 # the camera center in wobble mode. Note, that a cut in the poiting 132 # position directly influences the effective observation time and should 133 # only be applied if the events are lost anyhow. In both cases you must 134 # correct for the loss of observation time calculating your flux! 134 135 # ------------------------------------------------------------------------- 135 136 #Cut0.Condition: MSrcPosCam.GetDist*MGeomCam.fConvMm2Deg>0.3 136 137 137 138 138 # ---------- SETUP FOR ONOFF-MODE ----------- … … 183 183 184 184 # ------------------------------------------------------------------------- 185 # This cut is applied exclusively to the plots in the tab CutT. It is 186 # resonable to show the events within your signal region here to 187 # see if your signal region is clean from strange events (e.g. sparks) 188 # ------------------------------------------------------------------------- 189 CutT.Condition: ThetaSquared.fVal < ThetaSquaredCut.fVal 190 191 # ------------------------------------------------------------------------- 185 192 # Use this to executa a task (eg to calc hadronness) before Cut1 186 193 # ------------------------------------------------------------------------- -
trunk/MagicSoft/Mars/sponde.rc
r8680 r8705 1 EstimateEnergy.Rule: (0.380075+(MPointingPos.fZd*MPointingPos.fZd*0.00109028))*pow(MHillas.fSize,0.892462) 1 # -------------------- Define your energy estimator ------------------- 2 2 3 BinningSize.Raw: 18 49 63600 log 3 # Manually optimized (for example using trainenergy.C) 4 EstimateEnergy.Rule: MHillas.fSize^0.906*pow(10, -0.831*(MHillasSrc.fDist*3.37e-3-0.161)^2*(MHillasSrc.fDist*3.37e-3>0.161) + 2.21*(MHillasSrc.fDist*3.37e-3-0.493)^2*(MHillasSrc.fDist*3.37e-3>0.493) + 8.18e-6*MPointingPos.fZd^2.82 - 0.1*MNewImagePar.fLeakage1^0.2) 5 6 # You can use instead the Random Forest energy estimator 7 # (for example trained by trainenergy.C) 8 #EstimateEnergy: MRanForestCalc 9 #EstimateEnergy.NameOutput: MEnergyEst 10 #EstimateEnergy.FileName: rf-energy.root 11 #EstimateEnergy.Debug: No 12 13 # ------------------------ Define your binnings ----------------------- 14 15 # The energy binning is mandatory 4 16 BinningEnergyEst.Raw: 18 53 35800 log 5 17 6 #MMcSpectrumWeight.NewSlope: -2.26 18 # Binnings used for the Energy- and Disp-tab 19 #BinningSize.Raw: 50 10 100000 log 20 #BinningLeakage.Raw: 50 0 0.3 lin 21 #BinningEnergyResidual.Raw: 51 -1 1 lin 22 #BinningResidualDist.Raw: 51 -1 1 lin 23 24 # A change of the theta binning must be done in ganymed (see below) :( 25 #BinningTheta.Raw: 101 -0.005 0.505 asin 26 27 # --------------- Setup the weighting for your spectrum --------------- 28 # For more details MMcSpectrumWeight::ReadEnv 29 30 # A simple power law with a differential spectral index of -4 31 #MMcSpectrumWeight.NewSlope: -4 32 33 # Give a formula for the spectrum. Capital X is a place-holder 34 # for the energy (eg. the crab-spectrum from our publication) 35 #MMcSpectrumWeight.Formula: pow(X/300, -2.31-0.26*log10(X/300)) 36 37 # ----------------------- Additional setup ---------------------------- 38 39 # In the loop filling the ThetaSq histogram the theta cut is switched off 40 # and the theta-cut is implicitly done in the histogram. Reading the 41 # Monte Carlo for determining the collection area there is no theta- 42 # histogram involved, thus you need the theta cut switched on 43 MJSpectrum.Cut1.ThetaCut: On 44 45 # You may want to change the scale region to a different region than 46 # previously used in your ganymed-file. This might be usefull for 47 # on/off-data. Note, that the result can also depend on what else 48 # was setup for your MAlphaFitter in ganymed. 49 #MAlphaFitter.ScaleMin: 0.137 50 #MAlphaFitter.ScaleMax: 0.500 51 52 # You may use this to change the scale and/or fit region for 53 # determination of the excess vs energy (can be usefull in On/Off mode) 54 #MAlphaFitter.ScaleMode: Fixed 55 56 # You may want to allow scaling of each energy bin individually. (This 57 # can be useful if on- and off- data in wobble mode do not agree well) 58 #MAlphaFitter.ScaleMode: background 59 #MAlphaFitter.BackgroundFitMin: 0.09 60 #MAlphaFitter.BackgroundFitMax: 0.25 61 62 # ------------------------- Additional cut ---------------------------- 63 64 # You can apply a cut in addition to all other cuts to your data and 65 # MOnte Carlos. This can be usefull for example to setup a second 66 # MFMagicCuts with a tighter behaviour assuming that your ganymed 67 # was written with a really loose cut. 68 #CutS.Inverted: Yes 69 #CutS.Condition: MFMagicCuts 70 #CutS.ThetaCut: None 71 #CutS.HadronnesCut: Area 72 # Parametrization of Area 73 #CutS.Param2: 0.215468 74 #CutS.Param3: 5.63973 75 #CutS.Param4: 0.0836169 76 77 # --------------------------- Other binnings -------------------------- 78 79 # Binnings that are used by sponde but can not/shoud not be changed 80 # because they are taken from the ganymed file 81 #BinningTheta 82 #BinningWidth 83 #BinningLength 84 #BinningDist 85 #BinningM3Long 86 #BinningM3Trans 87 #BinningSlope 88 #BinningAsym 89 #BinningConc1 -
trunk/MagicSoft/Mars/sponde_onoff.rc
r8036 r8705 1 # Defaine your energy estimator 2 EstimateEnergy.Rule: 1.40871*pow(MHillas.fSize,0.810253) 1 # -------------------- Define your energy estimator ------------------- 3 2 3 # Manually optimized (for example using trainenergy.C) 4 EstimateEnergy.Rule: MHillas.fSize^0.906*pow(10, -0.831*(MHillasSrc.fDist*3.37e-3-0.161)^2*(MHillasSrc.fDist*3.37e-3>0.161) + 2.21*(MHillasSrc.fDist*3.37e-3-0.493)^2*(MHillasSrc.fDist*3.37e-3>0.493) + 8.18e-6*MPointingPos.fZd^2.82 - 0.1*MNewImagePar.fLeakage1^0.2) 5 6 # You can use instead the Random Forest energy estimator 7 # (for example trained by trainenergy.C) 4 8 #EstimateEnergy: MRanForestCalc 5 9 #EstimateEnergy.NameOutput: MEnergyEst 6 #EstimateEnergy.FileName: r anforest.root10 #EstimateEnergy.FileName: rf-energy.root 7 11 #EstimateEnergy.Debug: No 8 12 9 # Define your binnings 10 BinningSize.Raw: 12 49 6360 log 11 BinningEnergyEst.Raw: 12 53 3580 log 13 # ------------------------ Define your binnings ----------------------- 12 14 13 # Setup the weighting for your spectrum 14 # For more details MMcSpectrumWeight::ReadEnv 15 MMcSpectrumWeight.NewSlope: -2.6 15 # The energy binning is mandatory 16 BinningEnergyEst.Raw: 18 53 35800 log 16 17 17 # Give a formula for the spectrum. Capital X is a place-holder for the energy 18 #MMcSpectrumWeight.Formula: pow(X, -2.6) 18 # Binnings used for the Energy- and Disp-tab 19 #BinningSize.Raw: 50 10 100000 log 20 #BinningLeakage.Raw: 50 0 0.3 lin 21 #BinningEnergyResidual.Raw: 51 -1 1 lin 22 #BinningResidualDist.Raw: 51 -1 1 lin 23 24 # A change of the theta binning must be done in ganymed (see below) :( 25 #BinningTheta.Raw: 101 -0.005 0.505 asin 26 27 # --------------- Setup the weighting for your spectrum --------------- 28 # For more details MMcSpectrumWeight::ReadEnv 29 30 # A simple power law with a differential spectral index of -4 31 #MMcSpectrumWeight.NewSlope: -4 32 33 # Give a formula for the spectrum. Capital X is a place-holder 34 # for the energy (eg. the crab-spectrum from our publication) 35 #MMcSpectrumWeight.Formula: pow(X/300, -2.31-0.26*log10(X/300)) 36 37 # ----------------------- Additional setup ---------------------------- 19 38 20 39 # In the loop filling the ThetaSq histogram the theta cut is switched off 21 # Reading the Monte Carlo for determining the collection area you the 22 # Theta cut switched on 40 # and the theta-cut is implicitly done in the histogram. Reading the 41 # Monte Carlo for determining the collection area there is no theta- 42 # histogram involved, thus you need the theta cut switched on 23 43 MJSpectrum.Cut1.ThetaCut: On 24 44 25 # You may use this to change the scale and/or fit region 26 #MAlphaFitter.BackgroundFitMin: 0.12 27 #MAlphaFitter.BackgroundFitMax: 0.50 45 # You may want to change the scale region to a different region than 46 # previously used in your ganymed-file. This might be usefull for 47 # on/off-data. Note, that the result can also depend on what else 48 # was setup for your MAlphaFitter in ganymed. 49 #MAlphaFitter.ScaleMin: 0.137 50 #MAlphaFitter.ScaleMax: 0.500 28 51 29 # Use this if you like to use the global scale factor for all bins 52 # You may use this to change the scale and/or fit region for 53 # determination of the excess vs energy (can be usefull in On/Off mode) 30 54 #MAlphaFitter.ScaleMode: Fixed 55 56 # You may want to allow scaling of each energy bin individually. (This 57 # can be useful if on- and off- data in wobble mode do not agree well) 58 #MAlphaFitter.ScaleMode: background 59 #MAlphaFitter.BackgroundFitMin: 0.09 60 #MAlphaFitter.BackgroundFitMax: 0.25 61 62 # ------------------------- Additional cut ---------------------------- 63 64 # You can apply a cut in addition to all other cuts to your data and 65 # MOnte Carlos. This can be usefull for example to setup a second 66 # MFMagicCuts with a tighter behaviour assuming that your ganymed 67 # was written with a really loose cut. 68 #CutS.Inverted: Yes 69 #CutS.Condition: MFMagicCuts 70 #CutS.ThetaCut: None 71 #CutS.HadronnesCut: Area 72 # Parametrization of Area 73 #CutS.Param2: 0.215468 74 #CutS.Param3: 5.63973 75 #CutS.Param4: 0.0836169 76 77 # --------------------------- Other binnings -------------------------- 78 79 # Binnings that are used by sponde but can not/shoud not be changed 80 # because they are taken from the ganymed file 81 #BinningTheta 82 #BinningWidth 83 #BinningLength 84 #BinningDist 85 #BinningM3Long 86 #BinningM3Trans 87 #BinningSlope 88 #BinningAsym 89 #BinningConc1
Note:
See TracChangeset
for help on using the changeset viewer.