Add show_residuals.py unit test

Jürgen Knödlseder
1 parent 6e6c2ca9
Showing 5 changed files with 23 additions and 24 deletions
cscripts/Makefile.am
cscripts/cslightcrv.py
examples/show_residuals.py
test/data/residual.fits
test/test_examples.py.in
@@ -195,6 +195,7 @@ testdata = $(top_srcdir)/test/data/crab.xml \
            $(top_srcdir)/test/data/sensitivity_integral.dat \
            $(top_srcdir)/test/data/skymap.fits \
            $(top_srcdir)/test/data/spectrum.fits \
+           $(top_srcdir)/test/data/residual.fits \
            $(top_srcdir)/test/data/lightcurve.fits \
            $(top_srcdir)/test/data/ts.dat \
            $(top_srcdir)/test/data/models_spatial.xml \
@@ -558,12 +558,12 @@ class cslightcrv(ctools.csobservation):
             # fill results table with zeros
             else:
                 self._log_value(gammalib.TERSE, 'Warning',
-                                'No observations available in this time bin',
-                                'Skip bin.')
+                                'No observations available in this time bin, '
+                                'skip bin.')
  
                 # Set all results to 0
                 for par in pars:
-                    result['values'][par] = 0.0
+                    result['values'][par]        = 0.0
                     result['values']['e_' + par] = 0.0
  
                 # Append result
@@ -286,7 +286,7 @@ def show_residuals():
     Show residuals
     """
     # Set usage string
-    usage = 'show_spectrum.py [-p plotfile] [-h hdu] [file]'
+    usage = 'show_residual.py [-p plotfile] [-h hdu] [file]'
  
     # Set default options
     options = [{'option': '-p', 'value': ''},
@@ -297,7 +297,7 @@ def show_residuals():
  
     # Extract script parameters from options
     plotfile = options[0]['value']
-    hdu = options[1]['value']
+    hdu      = options[1]['value']
  
     # Show residuals
     plot_residuals(args[0], plotfile, hdu)
@@ -2,7 +2,7 @@
 # ==========================================================================
 # Performs unit tests for the ctools example scripts
 #
-# Copyright (C) 2016-2017 Juergen Knoedlseder
+# Copyright (C) 2016-2018 Juergen Knoedlseder
 #
 # This program is free software: you can redistribute it and/or modify
 # it under the terms of the GNU General Public License as published by
@@ -89,26 +89,8 @@ class Test(gammalib.GPythonTestSuite):
         self.name('examples')
  
         # Append tests
-        #self.append(self.test_generate_prod3_irfs,
-        #            'Test generate_prod3_irfs.py')
         self.append(self.test_make_pointings,
                     'Test make_pointings.py')
-        #self.append(self.test_make_spectrum,
-        #            'Test make_spectrum.py')
-        #self.append(self.test_make_ts_distributions,
-        #            'Test make_ts_distributions.py')
-        #self.append(self.test_pipeline_binned_disk,
-        #            'Test pipeline_binned_disk.py')
-        #self.append(self.test_pipeline_binned_mem,
-        #            'Test pipeline_binned_mem.py')
-        #self.append(self.test_pipeline_stacked_disk,
-        #            'Test pipeline_stacked_disk.py')
-        #self.append(self.test_pipeline_stacked_mem,
-        #            'Test pipeline_stacked_mem.py')
-        #self.append(self.test_pipeline_unbinned_disk,
-        #            'Test pipeline_unbinned_disk.py')
-        #self.append(self.test_pipeline_unbinned_mem,
-        #            'Test pipeline_unbinned_mem.py')
         self.append(self.test_show_butterfly,
                     'Test show_butterfly.py')
         self.append(self.test_show_irf,
@@ -123,6 +105,8 @@ class Test(gammalib.GPythonTestSuite):
                     'Test show_pull_evolution.py')
         self.append(self.test_show_pull_histogram,
                     'Test show_pull_histogram.py')
+        self.append(self.test_show_residuals,
+                    'Test show_residuals.py')
         self.append(self.test_show_response,
                     'Test show_response.py')
         self.append(self.test_show_sensitivity,
@@ -384,6 +368,20 @@ class Test(gammalib.GPythonTestSuite):
         # Return
         return
  
+    # Test show_residuals
+    def test_show_residuals(self):
+        """
+        Test show_residuals
+        """
+        # Set arguments
+        args = '-p example_show_residuals.eps %s/residual.fits' % self._datadir
+
+        # Execute script
+        self._execute_python('show_residuals', args=args)
+
+        # Return
+        return
+
     # Test show_response
     def test_show_response(self):
         """