Make sure that GCTAEdisp2D energy dispersion is never negative

Jürgen Knödlseder
1 parent ed7efb06
Showing 1 changed file with 5 additions and 0 deletions
inst/cta/src/GCTAEdisp2D.cpp
@@ -206,6 +206,11 @@ double GCTAEdisp2D::operator()(const double&amp; logEobs,
     // Compute edisp
     edisp = m_edisp(m_inx_matrix, arg[0], arg[1], arg[2]);
  
+    // Make sure that energy dispersion is non-negative
+    if (edisp < 0.0) {
+        edisp = 0.0;
+    }
+
     // Return
     return edisp;
 }