diff --git a/inst/cta/src/GCTAResponse.cpp b/inst/cta/src/GCTAResponse.cpp
index 65fcfe1..05865cd 100644
--- a/inst/cta/src/GCTAResponse.cpp
+++ b/inst/cta/src/GCTAResponse.cpp
@@ -1817,7 +1817,8 @@ double GCTAResponse::npred_diffuse(const GSource& source,
 
             // Integrate over theta
             GIntegral integral(&integrand);
-            integral.eps(1.0e-4);
+            //integral.eps(1.0e-4);
+            integral.eps(1.0e-5);
             npred = integral.romb(0.0, roi_psf_radius);
 
             // Apply deadtime correction
diff --git a/inst/cta/src/GCTAResponse_helpers.cpp b/inst/cta/src/GCTAResponse_helpers.cpp
index a963ac1..913fc88 100644
--- a/inst/cta/src/GCTAResponse_helpers.cpp
+++ b/inst/cta/src/GCTAResponse_helpers.cpp
@@ -966,7 +966,8 @@ double cta_npred_diffuse_kern_theta::eval(const double& theta)
 
         // Integrate over phi
         GIntegral integral(&integrand);
-        integral.eps(1.0e-4);
+        //integral.eps(1.0e-4);
+        integral.eps(1.0e-5);
         npred = integral.romb(0.0, gammalib::twopi) * sin_theta;
 
         // Debug: Check for NaN