Changeset 8705 for trunk


Ignore:
Timestamp:
08/24/07 13:44:23 (17 years ago)
Author:
tbretz
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/ganymed.rc

    r8668 r8705  
    128128
    129129# -------------------------------------------------------------------------
    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!
    134135# -------------------------------------------------------------------------
    135136#Cut0.Condition: MSrcPosCam.GetDist*MGeomCam.fConvMm2Deg>0.3
    136 
    137137
    138138# ---------- SETUP FOR ONOFF-MODE -----------
     
    183183
    184184# -------------------------------------------------------------------------
     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# -------------------------------------------------------------------------
     189CutT.Condition: ThetaSquared.fVal < ThetaSquaredCut.fVal
     190
     191# -------------------------------------------------------------------------
    185192# Use this to executa a task (eg to calc hadronness) before Cut1
    186193# -------------------------------------------------------------------------
  • 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 -------------------
    22
    3 BinningSize.Raw:        18   49   63600 log
     3# Manually optimized (for example using trainenergy.C)
     4EstimateEnergy.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
    416BinningEnergyEst.Raw:   18   53   35800 log
    517
    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
     43MJSpectrum.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 -------------------
    32
     3# Manually optimized (for example using trainenergy.C)
     4EstimateEnergy.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)
    48#EstimateEnergy: MRanForestCalc
    59#EstimateEnergy.NameOutput: MEnergyEst
    6 #EstimateEnergy.FileName: ranforest.root
     10#EstimateEnergy.FileName: rf-energy.root
    711#EstimateEnergy.Debug: No
    812
    9 # Define your binnings
    10 BinningSize.Raw:        12   49   6360 log
    11 BinningEnergyEst.Raw:   12   53   3580 log
     13# ------------------------ Define your binnings -----------------------
    1214
    13 # Setup the weighting for your spectrum
    14 # For more details MMcSpectrumWeight::ReadEnv
    15 MMcSpectrumWeight.NewSlope: -2.6
     15# The energy binning is mandatory
     16BinningEnergyEst.Raw:   18   53   35800 log
    1617
    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 ----------------------------
    1938
    2039# 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
    2343MJSpectrum.Cut1.ThetaCut: On
    2444
    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
    2851
    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)
    3054#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.