Adding an rho integration boundary at rho=semiminor solved the problem with the …

…strange Npred profile.

Adding an rho integration boundary at rho=semiminor solved the problem with the …
…strange Npred profile.
Jürgen Knödlseder
1 parent 89b294ab
Showing 1 changed file with 6 additions and 0 deletions
inst/cta/src/GCTAResponseIrf.cpp
@@ -2409,6 +2409,12 @@ double GCTAResponseIrf::npred_elliptical(const GSource&amp; source,
             bounds.push_back(transition_point);
         }
  
+        // If the integration range includes the semiminor boundary, then
+        // add an integration boundary at that location
+        if (semiminor > rho_min && semiminor < rho_max) {
+            bounds.push_back(semiminor);
+        }
+
         // Integrate kernel
         npred = integral.romberg(bounds, iter_rho);