Do no longer use a finer energy binning for the background cube in obsutils.get_stacked_response().

The interpolation of the stacked background cube is problematic when combining various energy thresholds. The background cube should *not* be interpolated but should be identical to the counts cube. Both are defined in the same data space.

Do no longer use a finer energy binning for the background cube in obsutils.get_stacked_response().
The interpolation of the stacked background cube is problematic when combining various energy thresholds. The background cube should *not* be interpolated but should be identical to the counts cube. Both are defined in the same data space.
Jürgen Knödlseder
1 parent 83fd1279
Showing 1 changed file with 4 additions and 2 deletions
cscripts/obsutils.py
@@ -633,7 +633,9 @@ def get_stacked_response(obs, xref, yref, binsz=0.05, nxpix=200, nypix=200,
         psfcube.logFileOpen()
     psfcube.run()
  
-    # Create background cube
+    # Create background cube. Note that we use the same energy binning as
+    # for the counts cube since no interpolation should be actually done
+    # for the background cube.
     bkgcube = ctools.ctbkgcube(obs)
     bkgcube['incube']    = 'NONE'
     bkgcube['usepnt']    = usepnt
@@ -641,7 +643,7 @@ def get_stacked_response(obs, xref, yref, binsz=0.05, nxpix=200, nypix=200,
     bkgcube['binsz']     = binsz
     bkgcube['nxpix']     = nxpix
     bkgcube['nypix']     = nypix
-    bkgcube['enumbins']  = _enumbins
+    bkgcube['enumbins']  = enumbins
     bkgcube['emin']      = emin
     bkgcube['emax']      = emax
     bkgcube['coordsys']  = coordsys