RING background fix in ctskymap (#2839)

Jürgen Knödlseder
1 parent f11482c0
Showing 4 changed files with 8 additions and 2 deletions
ChangeLog
NEWS
doc/source/admin/release_history/1.6.rst
src/ctskymap/ctskymap.cpp
@@ -4,6 +4,7 @@
         ========================
         Log warning if edisp=no (#2842)
+        Fix ring background computation in ctskymap for small roiradius (#2839)
         Add max_iter parameter to ctbutterfly, cterror, ctlike, cttsmap,
           ctulimit and cstsmapsplit
         ctbkgcube now requires an input counts cube (#2763)
@@ -216,7 +216,9 @@ of the input observations (#2054).
 ctskymap - CTA sky mapping
 --------------------------
-Add exclusion_map() methods (#2725)
+Fix an issue with artefacts for small roiradius values (#2839).
+
+Add exclusion_map() methods (#2725).
 Add the "inmap", "iterations" and "threshold" parameters. The hidden
 "inmap" parameter allows specifying a sky map file containing a "COUNTS"
@@ -22,6 +22,8 @@ In particular, this release provides:
 Bug fixes
 ---------
+* [`2839 <https://cta-redmine.irap.omp.eu/issues/2839>`_] -
+  Fix ring background computation in :ref:`ctskymap` for small roiradius
 * [`2721 <https://cta-redmine.irap.omp.eu/issues/2721>`_] -
   Apply event bin weighting in :ref:`ctmodel` for stacked computation
 * [`2694 <https://cta-redmine.irap.omp.eu/issues/2694>`_] -
@@ -1371,7 +1371,8 @@ void ctskymap::compute_ring_values(const int&amp;     ipixel,
                 }
                 // ... otherwise check if pixel is inside source region
-                else if (position.cos_dist(skydir) >= m_cos_roiradius) {
+                else if ((position.cos_dist(skydir) >= m_cos_roiradius) ||
+                         (ipixel == i)) {
                     // Update n_on
                     non += counts(i);