Commit 55f9af023a891010d24f07b686d52095c861f40f

Authored by Jürgen Knödlseder
1 parent 6f6cc745

Set default number of angular bins for PSF cube to 200 and reduce spatial resolu…

…tion by a factor of 20 w/r to the counts cube.

Also document the new anumbins parameter.
cscripts/csspec.py
... ... @@ -122,7 +122,7 @@ class csspec(ctools.cscript):
122 122 pars.append(gammalib.GApplicationPar("nxpix","i","a","200","","","Size of the X axis in pixels"))
123 123 pars.append(gammalib.GApplicationPar("nypix","i","a","200","","","Size of the Y axis in pixels"))
124 124 pars.append(gammalib.GApplicationPar("binsz","r","a","0.02","","","Pixel size (deg/pixel)"))
125   - pars.append(gammalib.GApplicationPar("anumbins","i","h","20","","","Angular separation bins for PSF cube in binned mode"))
  125 + pars.append(gammalib.GApplicationPar("anumbins","i","h","200","","","Angular separation bins for PSF cube in binned mode"))
126 126 pars.append(gammalib.GApplicationPar("calc_ts","b","h","yes","yes|no","","Compute TS value in each bin"))
127 127 pars.append(gammalib.GApplicationPar("calc_ulim","b","h","yes","yes|no","","Compute upper limit in each bin"))
128 128 pars.append(gammalib.GApplicationPar("fix_srcs","b","h","yes","yes|no","","Fix other skymodel parameters"))
... ... @@ -404,9 +404,9 @@ class csspec(ctools.cscript):
404 404 psfcube["edisp"] = self.m_edisp
405 405 psfcube["xref"] = self.m_xref
406 406 psfcube["yref"] = self.m_yref
407   - psfcube["binsz"] = self.m_binsz
408   - psfcube["nxpix"] = self.m_nxpix
409   - psfcube["nypix"] = self.m_nypix
  407 + psfcube["binsz"] = self.m_binsz*20.0 # Slow Psf variation
  408 + psfcube["nxpix"] = self.m_nxpix/20
  409 + psfcube["nypix"] = self.m_nypix/20
410 410 psfcube["enumbins"] = self.m_ebins
411 411 psfcube["anumbins"] = self.m_anumbins
412 412 psfcube["emin"] = emin.TeV()
... ...
doc/source/reference_manual/csspec.rst
... ... @@ -81,6 +81,9 @@ General parameters
81 81  
82 82 ``binsz [real]``
83 83 Pixel size (in degrees/pixel).
  84 +
  85 +``(anumbins = 200) [integer]``
  86 + Number of angular separation bins.
84 87  
85 88 ``(calc_ts = yes) [boolean]``
86 89 Compute TS for each spectral point?
... ...