development of test unit

Luigi Tibaldo
1 parent 63e7f3f8
Showing 4 changed files with 594 additions and 5 deletions
cscripts/csphagen.py
test/data/crab_offaxis.fits → test/data/crab_offaxis1.fits
test/data/crab_offaxis2.fits
test/test_csphagen.py
@@ -92,6 +92,8 @@ class csphagen(ctools.cscript):
             self._excl_reg = gammalib.GSkyRegionMap(
                 self["exclusion"].filename())
             self._has_exclusion = True
+        else:
+            self._has_exclusion = False
  
         # Stacking
         self._stack = self['stack'].boolean()
@@ -240,9 +242,6 @@ class csphagen(ctools.cscript):
         # Run the script
         self.run()
  
-        # Save residual map
-        self.save()
-
         # Return
         return
  
@@ -41,10 +41,11 @@ class Test(test):
         """
         # Call base class constructor
         test.__init__(self)
-        self._myevents = self._datadir + '/crab_offaxis.fits'
+        self._myevents = self._datadir + '/crab_offaxis1.fits'
         self._exclusion = self._datadir + '/crab_exclusion.fits'
         self._nreg_with_excl = 5
-        # number of expected background regions with exclusion
+        self._nreg_wo_excl = 8
+        # number of expected background regions with/wo exclusion
  
         # Return
         return
@@ -108,6 +109,62 @@ class Test(test):
         # Return
         return
  
+    def _test_python(self):
+        """
+        Test csphagen from Python
+        """
+        nbins = 120
+
+        # Same test as from command line
+        phagen = cscripts.csphagen()
+        phagen['inobs'] = self._myevents
+        phagen['caldb'] = self._caldb
+        phagen['irf'] = self._irf
+        phagen['ebinalg'] = 'LOG'
+        phagen['emin'] = 0.1
+        phagen['emax']= 100.
+        phagen['enumbins'] = nbins
+        phagen['coordsys'] = 'CEL'
+        phagen['ra'] = 83.633
+        phagen['dec'] = 22.0145
+        phagen['rad'] = 0.2
+        phagen['stack'] = False
+        phagen['exclusion'] = self._exclusion
+        phagen['outroot'] = 'genpha_py1'
+        phagen['logfile'] = 'csphagen_py1.log'
+        phagen['chatter'] = 1
+
+        # Run script
+        phagen.execute()
+
+        # Check outout
+        self._check_output('genpha_py1', nbins, self._nreg_with_excl)
+
+        # Second test, now without exclusion region
+        phagen = cscripts.csphagen()
+        phagen['inobs'] = self._myevents
+        phagen['caldb'] = self._caldb
+        phagen['irf'] = self._irf
+        phagen['ebinalg'] = 'LOG'
+        phagen['emin'] = 0.1
+        phagen['emax'] = 100.
+        phagen['enumbins'] = nbins
+        phagen['coordsys'] = 'CEL'
+        phagen['ra'] = 83.633
+        phagen['dec'] = 22.0145
+        phagen['rad'] = 0.2
+        phagen['stack'] = False
+        phagen['outroot'] = 'genpha_py2'
+        phagen['logfile'] = 'csphagen_py2.log'
+        phagen['chatter'] = 1
+
+        # Run script
+        phagen.execute()
+
+        # Check outout
+        self._check_output('genpha_py2', nbins, self._nreg_wo_excl)
+
+
     def _check_ebounds(self, table, bins):
         """
         Check EBOUNDS table