Justine Devin / gammalib

Commit 8b05373833c6e97160ee79e85dca669f06fac7c9

Authored by Jürgen Knödlseder
1 parent a70c7a78

Set integration precision for diffuse Npred computation to 1e-5.

Inline Side-by-side
Showing 2 changed files with 4 additions and 2 deletions
inst/cta/src/GCTAResponse.cpp
  View file @8b05373
... ... @@ -1817,7 +1817,8 @@ double GCTAResponse::npred_diffuse(const GSource& source,
1817 1817  
1818 1818 // Integrate over theta
1819 1819 GIntegral integral(&integrand);
1820   - integral.eps(1.0e-4);
  1820 + //integral.eps(1.0e-4);
  1821 + integral.eps(1.0e-5);
1821 1822 npred = integral.romb(0.0, roi_psf_radius);
1822 1823  
1823 1824 // Apply deadtime correction
... ...
inst/cta/src/GCTAResponse_helpers.cpp
  View file @8b05373
... ... @@ -966,7 +966,8 @@ double cta_npred_diffuse_kern_theta::eval(const double& theta)
966 966  
967 967 // Integrate over phi
968 968 GIntegral integral(&integrand);
969   - integral.eps(1.0e-4);
  969 + //integral.eps(1.0e-4);
  970 + integral.eps(1.0e-5);
970 971 npred = integral.romb(0.0, gammalib::twopi) * sin_theta;
971 972  
972 973 // Debug: Check for NaN
... ...