added basic test to check that containment radius is > 0.0 and < 1e6

nkelhos
1 parent b5029daa
Showing 1 changed file with 7 additions and 0 deletions
inst/cta/test/test_CTA.cpp
@@ -390,6 +390,12 @@ void TestGCTAResponse::test_response_psf(void)
                    gammalib::deg2rad;
         }
         test_value(sum, 1.0, 0.001, "PSF integration for "+eng.print());
+        
+        // test containment radii
+        for ( double f = 0.01 ; f < 1.0 ; f += 0.07 ) {
+            double radius = rsp.psf()->containment_radius( f, eng.log10TeV() ) ;
+            test_assert( radius>0.0 && radius < 1e6 , "PSF containment radius for fraction="+gammalib::str(f) ) ;
+        }
     }
     // Return
@@ -428,6 +434,7 @@ void TestGCTAResponse::test_response_psf_king(void)
                    gammalib::deg2rad;
         }
         test_value(sum, 1.0, 0.001, "PSF integration for "+eng.print());
+        
     }
     // Return