Fix exception when calling ctools clear() methods (#1835)

Jürgen Knödlseder
1 parent 04b1443d
Showing 24 changed files with 319 additions and 86 deletions
ChangeLog
NEWS
README.md
src/ctbin/ctbin.cpp
src/ctbkgcube/ctbkgcube.cpp
src/ctbutterfly/ctbutterfly.cpp
src/ctcubemask/ctcubemask.cpp
src/ctedispcube/ctedispcube.cpp
src/cterror/cterror.cpp
src/ctexpcube/ctexpcube.cpp
src/ctlike/ctlike.cpp
src/ctmapcube/ctmapcube.cpp
src/ctmodel/ctmodel.cpp
src/ctobssim/ctobssim.cpp
src/ctpsfcube/ctpsfcube.cpp
src/ctselect/ctselect.cpp
src/ctskymap/ctskymap.cpp
src/cttsmap/cttsmap.cpp
src/ctulimit/ctulimit.cpp
src/support/ctool.cpp
-2016-07-31
+2016-08-04
         * Version 1.2.0 released
         ========================
+        Fix exception when calling ctools clear() methods (#1835)
         Add support modules to cscripts (#1830)
         Fix lacking coverage measure for cscripts execution (#1827)
         Set "addbounds" parameters by default to "no" (#1819)
 New Features and Important Changes in ctools 1.2.0
-31 July 2016
+4 August 2016
 Introduction
@@ -27,82 +27,106 @@ None
 ctobssim - CTA observation simulation
 -------------------------------------
-None
+The clear() method now calls GApplication::clear() to conserve the properties
+of the tool (#1835).
 ctselect - CTA event selection
 ------------------------------
-None
+The clear() method now calls GApplication::clear() to conserve the properties
+of the tool (#1835).
 ctbin - CTA event binning
 -------------------------
-None
+The clear() method now calls GApplication::clear() to conserve the properties
+of the tool (#1835).
 ctlike - Maximum likelihood fitting
 -----------------------------------
-None
+The clear() method now calls GApplication::clear() to conserve the properties
+of the tool (#1835).
 cttsmap - Generation of Test Statistic map
 ------------------------------------------
-None
+The clear() method now calls GApplication::clear() to conserve the properties
+of the tool (#1835).
 ctskymap - CTA sky mapping
 --------------------------
-None
+The clear() method now calls GApplication::clear() to conserve the properties
+of the tool (#1835).
 ctmodel - CTA counts cube modelling
 -----------------------------------
-None
+The clear() method now calls GApplication::clear() to conserve the properties
+of the tool (#1835).
 ctexpcube - CTA exposure cube generation
 ----------------------------------------
 Set hidden "addbounds" parameter by default to "no" (#1819).
+The clear() method now calls GApplication::clear() to conserve the properties
+of the tool (#1835).
+
 ctpsfcube - CTA point spread function cube generation
 -----------------------------------------------------
 Set hidden "addbounds" parameter by default to "no" (#1819).
+The clear() method now calls GApplication::clear() to conserve the properties
+of the tool (#1835).
+
 ctedispcube - CTA energy dispersion cube generation
 ---------------------------------------------------
 Set hidden "addbounds" parameter by default to "no" (#1819).
+The clear() method now calls GApplication::clear() to conserve the properties
+of the tool (#1835).
+
 ctbkgcube - CTA background cube generation
 ------------------------------------------
 Set hidden "addbounds" parameter by default to "no" (#1819).
+The clear() method now calls GApplication::clear() to conserve the properties
+of the tool (#1835).
+
 ctmapcube - CTA map cube generation
 -----------------------------------
-None
+The clear() method now calls GApplication::clear() to conserve the properties
+of the tool (#1835).
 ctcubemask - CTA cube filtering
 -------------------------------
-None
+The clear() method now calls GApplication::clear() to conserve the properties
+of the tool (#1835).
 ctbutterfly - Generate butterfly diagrams
 -----------------------------------------
-None
+The clear() method now calls GApplication::clear() to conserve the properties
+of the tool (#1835).
 ctulimit - Compute upper limit
 ------------------------------
-None
+The clear() method now calls GApplication::clear() to conserve the properties
+of the tool (#1835).
 cterror - Likelihood profile error computation
 ----------------------------------------------
-None
+The clear() method now calls GApplication::clear() to conserve the properties
+of the tool (#1835).
 cscripts
 ctools information
 ==================
-* Version:             1.2.0.dev1 (31 July 2016)
+* Version:             1.2.0.dev1 (4 August 2016)
 * Author:              Juergen Knoedlseder (jurgen.knodlseder@irap.omp.eu)
 * GammaLib dependency: 1.2.0.dev1
@@ -195,13 +195,17 @@ void ctbin::clear(void)
     // Free members
     free_members();
     this->ctool::free_members();
-    this->GApplication::free_members();
+
+    // Clear base class (needed to conserve tool name and version)
+    this->GApplication::clear();
     // Initialise members
-    this->GApplication::init_members();
     this->ctool::init_members();
     init_members();
+    // Write header into logger
+    log_header();
+
     // Return
     return;
 }
@@ -175,20 +175,26 @@ ctbkgcube&amp; ctbkgcube::operator=(const ctbkgcube&amp; app)
  ==========================================================================*/
 /***********************************************************************//**
- * @brief Clear instance
+ * @brief Clear ctbkgcube tool
+ *
+ * Clears ctbkgcube tool.
  ***************************************************************************/
 void ctbkgcube::clear(void)
 {
     // Free members
     free_members();
     this->ctool::free_members();
-    this->GApplication::free_members();
+
+    // Clear base class (needed to conserve tool name and version)
+    this->GApplication::clear();
     // Initialise members
-    this->GApplication::init_members();
     this->ctool::init_members();
     init_members();
+    // Write header into logger
+    log_header();
+
     // Return
     return;
 }
@@ -176,20 +176,26 @@ ctbutterfly&amp; ctbutterfly::operator=(const ctbutterfly&amp; app)
  ==========================================================================*/
 /***********************************************************************//**
- * @brief Clear instance
+ * @brief Clear ctbutterfly tool
+ *
+ * Clears ctbutterfly tool.
  ***************************************************************************/
 void ctbutterfly::clear(void)
 {
     // Free members
     free_members();
     this->ctool::free_members();
-    this->GApplication::free_members();
+
+    // Clear base class (needed to conserve tool name and version)
+    this->GApplication::clear();
     // Initialise members
-    this->GApplication::init_members();
     this->ctool::init_members();
     init_members();
+    // Write header into logger
+    log_header();
+
     // Return
     return;
 }
@@ -313,10 +319,6 @@ void ctbutterfly::run(void)
     double cov_pg = m_covariance(inx_prefactor,inx_index);
     double cov_gg = m_covariance(inx_index,inx_index);
-    // Initialise min and max intensities
-    m_min_intensities.assign(m_ebounds.size(), 1e30);
-    m_max_intensities.assign(m_ebounds.size(), 0.0);
-
     // Compute eigenvectors and eigenvalues
     double  lambda1;
     double  lambda2;
@@ -346,6 +348,12 @@ void ctbutterfly::run(void)
     // Write header into logger
     log_header1(TERSE, "Generate butterfly");
+    // Initialise result arrays
+    m_energies.assign(m_ebounds.size(), 0.0);
+    m_intensities.assign(m_ebounds.size(), 0.0);
+    m_min_intensities.assign(m_ebounds.size(), 1e30);
+    m_max_intensities.assign(m_ebounds.size(), 0.0);
+
     // Walk around the error ellipse
     int    nt = 360;
     double dt = gammalib::twopi/double(nt);
@@ -382,7 +390,7 @@ void ctbutterfly::run(void)
             // Evaluate power law
             double intensity = plaw.eval(energy, time);
-            
+
             // Get extremes
             if (intensity < m_min_intensities[k]) {
                 m_min_intensities[k] = intensity;
@@ -395,13 +403,10 @@ void ctbutterfly::run(void)
     } // endfor: looped over all angles
-    // Write header into logger
-    log_header1(TERSE, "Results");
-
-    // Setup power law model
+    // Setup power law model for mean prefactor and index
     GModelSpectralPlaw plaw(prefactor_mean, index_mean, pivot);
-    // Show results
+    // Compute energy and intensity vectors for mean power law
     for (int k = 0; k < m_ebounds.size(); ++k) {
         // Get the energy of current bin
@@ -412,16 +417,24 @@ void ctbutterfly::run(void)
         double intensity = plaw.eval(energy, time);
         // Store results
-        m_energies.push_back(energy.MeV());
-        m_intensities.push_back(intensity);
+        m_energies[k]    = energy.MeV();
+        m_intensities[k] = intensity;
+
+    } // endfor: looped over energies
+
+    // Write header into logger
+    log_header1(TERSE, "Results");
+
+    // Show results
+    for (int k = 0; k < m_ebounds.size(); ++k) {
         // Write results into logger
-        log_value(NORMAL, "Intensity at "+energy.print(),
-                  gammalib::str(intensity)+" ("+
+        log_value(NORMAL, "Intensity at "+gammalib::str(m_energies[k])+" MeV",
+                  gammalib::str(m_intensities[k])+" ("+
                   gammalib::str(m_min_intensities[k])+", "+
                   gammalib::str(m_max_intensities[k])+") ph/cm2/s/MeV");
-    } // endfor
+    } // endfor: looped over energies
     // Restore energy dispersion flag for all CTA observations
     for (int i = 0; i < m_obs.size(); ++i) {
@@ -177,20 +177,26 @@ ctcubemask&amp; ctcubemask::operator=(const ctcubemask&amp; app)
  ==========================================================================*/
 /***********************************************************************//**
- * @brief Clear instance
+ * @brief Clear ctcubemask tool
+ *
+ * Clears ctcubemask tool.
  ***************************************************************************/
 void ctcubemask::clear(void)
 {
     // Free members
     free_members();
     this->ctool::free_members();
-    this->GApplication::free_members();
+
+    // Clear base class (needed to conserve tool name and version)
+    this->GApplication::clear();
     // Initialise members
-    this->GApplication::init_members();
     this->ctool::init_members();
     init_members();
+    // Write header into logger
+    log_header();
+
     // Return
     return;
 }
@@ -195,13 +195,17 @@ void ctedispcube::clear(void)
     // Free members
     free_members();
     this->ctool::free_members();
-    this->GApplication::free_members();
+
+    // Clear base class (needed to conserve tool name and version)
+    this->GApplication::clear();
     // Initialise members
-    this->GApplication::init_members();
     this->ctool::init_members();
     init_members();
+    // Write header into logger
+    log_header();
+
     // Return
     return;
 }
@@ -188,22 +188,26 @@ cterror&amp; cterror::operator=(const cterror&amp; app)
  ==========================================================================*/
 /***********************************************************************//**
- * @brief Clear instance
+ * @brief Clear cterror tool
  *
- * Resets the cterror tool to a clean initial state.
+ * Clears cterror tool.
  ***************************************************************************/
 void cterror::clear(void)
 {
     // Free members
     free_members();
     this->ctool::free_members();
-    this->GApplication::free_members();
+
+    // Clear base class (needed to conserve tool name and version)
+    this->GApplication::clear();
     // Initialise members
-    this->GApplication::init_members();
     this->ctool::init_members();
     init_members();
+    // Write header into logger
+    log_header();
+
     // Return
     return;
 }
@@ -175,20 +175,26 @@ ctexpcube&amp; ctexpcube::operator=(const ctexpcube&amp; app)
  ==========================================================================*/
 /***********************************************************************//**
- * @brief Clear instance
+ * @brief Clear ctexpcube tool
+ *
+ * Clears ctexpcube tool.
  ***************************************************************************/
 void ctexpcube::clear(void)
 {
     // Free members
     free_members();
     this->ctool::free_members();
-    this->GApplication::free_members();
+
+    // Clear base class (needed to conserve tool name and version)
+    this->GApplication::clear();
     // Initialise members
-    this->GApplication::init_members();
     this->ctool::init_members();
     init_members();
+    // Write header into logger
+    log_header();
+
     // Return
     return;
 }
@@ -171,20 +171,26 @@ ctlike&amp; ctlike::operator=(const ctlike&amp; app)
  ==========================================================================*/
 /***********************************************************************//**
- * @brief Clear instance
+ * @brief Clear ctlike tool
+ *
+ * Clears ctlike tool.
  ***************************************************************************/
 void ctlike::clear(void)
 {
     // Free members
     free_members();
     this->ctool::free_members();
-    this->GApplication::free_members();
+
+    // Clear base class (needed to conserve tool name and version)
+    this->GApplication::clear();
     // Initialise members
-    this->GApplication::init_members();
     this->ctool::init_members();
     init_members();
+    // Write header into logger
+    log_header();
+
     // Return
     return;
 }
@@ -151,20 +151,26 @@ ctmapcube&amp; ctmapcube::operator=(const ctmapcube&amp; app)
  ==========================================================================*/
 /***********************************************************************//**
- * @brief Clear tool
+ * @brief Clear ctmapcube tool
+ *
+ * Clears ctmapcube tool.
  ***************************************************************************/
 void ctmapcube::clear(void)
 {
     // Free members
     free_members();
     this->ctool::free_members();
-    this->GApplication::free_members();
+
+    // Clear base class (needed to conserve tool name and version)
+    this->GApplication::clear();
     // Initialise members
-    this->GApplication::init_members();
     this->ctool::init_members();
     init_members();
+    // Write header into logger
+    log_header();
+
     // Return
     return;
 }
@@ -175,20 +175,26 @@ ctmodel&amp; ctmodel::operator=(const ctmodel&amp; app)
  ==========================================================================*/
 /***********************************************************************//**
- * @brief Clear instance
+ * @brief Clear ctmodel tool
+ *
+ * Clears ctmodel tool.
  ***************************************************************************/
 void ctmodel::clear(void)
 {
     // Free members
     free_members();
     this->ctool::free_members();
-    this->GApplication::free_members();
+
+    // Clear base class (needed to conserve tool name and version)
+    this->GApplication::clear();
     // Initialise members
-    this->GApplication::init_members();
     this->ctool::init_members();
     init_members();
+    // Write header into logger
+    log_header();
+
     // Return
     return;
 }
@@ -203,13 +203,17 @@ void ctobssim::clear(void)
     // Free members
     free_members();
     this->ctool::free_members();
-    this->GApplication::free_members();
+
+    // Clear base class (needed to conserve tool name and version)
+    this->GApplication::clear();
     // Initialise members
-    this->GApplication::init_members();
     this->ctool::init_members();
     init_members();
+    // Write header into logger
+    log_header();
+
     // Return
     return;
 }
@@ -174,20 +174,26 @@ ctpsfcube&amp; ctpsfcube::operator=(const ctpsfcube&amp; app)
  ==========================================================================*/
 /***********************************************************************//**
- * @brief Clear instance
+ * @brief Clear ctpsfcube tool
+ *
+ * Clears ctpsfcube tool.
  ***************************************************************************/
 void ctpsfcube::clear(void)
 {
     // Free members
     free_members();
     this->ctool::free_members();
-    this->GApplication::free_members();
+
+    // Clear base class (needed to conserve tool name and version)
+    this->GApplication::clear();
     // Initialise members
-    this->GApplication::init_members();
     this->ctool::init_members();
     init_members();
+    // Write header into logger
+    log_header();
+
     // Return
     return;
 }
@@ -179,20 +179,26 @@ ctselect&amp; ctselect::operator=(const ctselect&amp; app)
  ==========================================================================*/
 /***********************************************************************//**
- * @brief Clear instance
+ * @brief Clear ctselect tool
+ *
+ * Clears ctselect tool.
  ***************************************************************************/
 void ctselect::clear(void)
 {
     // Free members
     free_members();
     this->ctool::free_members();
-    this->GApplication::free_members();
+
+    // Clear base class (needed to conserve tool name and version)
+    this->GApplication::clear();
     // Initialise members
-    this->GApplication::init_members();
     this->ctool::init_members();
     init_members();
+    // Write header into logger
+    log_header();
+
     // Return
     return;
 }
@@ -175,20 +175,26 @@ ctskymap&amp; ctskymap::operator=(const ctskymap&amp; app)
  ==========================================================================*/
 /***********************************************************************//**
- * @brief Clear instance
+ * @brief Clear ctskymap tool
+ *
+ * Clears ctskymap tool.
  ***************************************************************************/
 void ctskymap::clear(void)
 {
     // Free members
     free_members();
     this->ctool::free_members();
-    this->GApplication::free_members();
+
+    // Clear base class (needed to conserve tool name and version)
+    this->GApplication::clear();
     // Initialise members
-    this->GApplication::init_members();
     this->ctool::init_members();
     init_members();
+    // Write header into logger
+    log_header();
+
     // Return
     return;
 }
@@ -173,20 +173,26 @@ cttsmap&amp; cttsmap::operator=(const cttsmap&amp; app)
  ==========================================================================*/
 /***********************************************************************//**
- * @brief Clear instance
+ * @brief Clear cttsmap tool
+ *
+ * Clears cttsmap tool.
  ***************************************************************************/
 void cttsmap::clear(void)
 {
     // Free members
     free_members();
     this->ctool::free_members();
-    this->GApplication::free_members();
+
+    // Clear base class (needed to conserve tool name and version)
+    this->GApplication::clear();
     // Initialise members
-    this->GApplication::init_members();
     this->ctool::init_members();
     init_members();
+    // Write header into logger
+    log_header();
+
     // Return
     return;
 }
@@ -198,13 +198,17 @@ void ctulimit::clear(void)
     // Free members
     free_members();
     this->ctool::free_members();
-    this->GApplication::free_members();
+
+    // Clear base class (needed to conserve tool name and version)
+    this->GApplication::clear();
     // Initialise members
-    this->GApplication::init_members();
     this->ctool::init_members();
     init_members();
+    // Write header into logger
+    log_header();
+
     // Return
     return;
 }
@@ -276,9 +276,7 @@ void ctool::copy_members(const ctool&amp; app)
 void ctool::free_members(void)
 {
     // Write separator into logger
-    if (logTerse()) {
-        log << std::endl;
-    }
+    log_string(TERSE, "\n");
     // Return
     return;
@@ -179,14 +179,14 @@ class Test(test):
         self._check_cube(evt, 5542)
         # Now clear ctbin tool
-        #bin.clear()         # Segfaults, see issue #1835
+        bin.clear()
         # Check that cleared ctbin has an empty observation container and
         # counts cube
-        #self.test_value(bin.obs().size(), 0,
-        #     'Check that empty ctbin has an empty observation container')
-        #self.test_value(bin.cube().size(), 0,
-        #     'Check that empty ctbin has an empty counts cube')
+        self.test_value(bin.obs().size(), 0,
+             'Check that empty ctbin has an empty observation container')
+        self.test_value(bin.cube().size(), 0,
+             'Check that empty ctbin has an empty counts cube')
         # Prepare observation container for stacked analysis
         cta = gammalib.GCTAObservation(self._events)
@@ -116,8 +116,6 @@ class Test(test):
         bkgcube['inobs']    = self._events
         bkgcube['inmodel']  = self._model
         bkgcube['incube']   = 'NONE'
-        bkgcube['outcube']  = 'ctbkgcube_py1.fits'
-        bkgcube['outmodel'] = 'ctbkgcube_py1.xml'
         bkgcube['caldb']    = self._caldb
         bkgcube['irf']      = self._irf
         bkgcube['ebinalg']  = 'LOG'
@@ -131,6 +129,8 @@ class Test(test):
         bkgcube['proj']     = 'CAR'
         bkgcube['xref']     = 83.63
         bkgcube['yref']     = 22.01
+        bkgcube['outcube']  = 'ctbkgcube_py1.fits'
+        bkgcube['outmodel'] = 'ctbkgcube_py1.xml'
         bkgcube['logfile']  = 'ctbkgcube_py1.log'
         bkgcube['chatter']  = 2
@@ -142,6 +142,87 @@ class Test(test):
         # Check result file
         self._check_result_file('ctbkgcube_py1.fits')
+        # Set-up ctbkgcube from counts cube
+        bkgcube = ctools.ctbkgcube()
+        bkgcube['inobs']     = self._events
+        bkgcube['inmodel']   = self._model
+        bkgcube['incube']    = self._cntcube
+        bkgcube['caldb']     = self._caldb
+        bkgcube['irf']       = self._irf
+        bkgcube['outcube']   = 'ctbkgcube_py2.fits'
+        bkgcube['outmodel']  = 'ctbkgcube_py2.xml'
+        bkgcube['logfile']   = 'ctbkgcube_py2.log'
+        bkgcube['chatter']   = 3
+        bkgcube['publish']   = True
+        bkgcube['addbounds'] = True
+
+        # Execute ctbkgcube tool
+        bkgcube.logFileOpen()   # Make sure we get a log file
+        bkgcube.execute()
+
+        # Check result file
+        self._check_result_file('ctbkgcube_py2.fits')
+
+        # Copy ctbkgcube tool and execute copied tool
+        cpy_bkgcube = bkgcube
+        cpy_bkgcube['outcube']  = 'ctbkgcube_py3.fits'
+        cpy_bkgcube['outmodel'] = 'ctbkgcube_py3.xml'
+        cpy_bkgcube['logfile']  = 'ctbkgcube_py3.log'
+        cpy_bkgcube['chatter']  = 4
+        cpy_bkgcube.execute()
+
+        # Check result file
+        self._check_result_file('ctbkgcube_py3.fits')
+
+        # Clear ctbkgcube tool
+        bkgcube.clear()
+
+        # TODO: Do some test after clearing of ctbkgcube tool
+
+        # Prepare observation container
+        obs1 = gammalib.GCTAObservation(self._events)
+        obs2 = gammalib.GCTAObservation(self._cntcube)
+        obs3 = gammalib.GLATObservation()
+        obs1.id('0001')
+        obs2.id('0002')
+        obs3.id('0001')
+        obs3.events(gammalib.GLATEventList())
+        obs = gammalib.GObservations()
+        obs.append(obs1)
+        obs.append(obs2)
+        obs.append(obs3)
+        obs.models(gammalib.GModels(self._model))
+
+        # Set-up ctbkgcube from observation container
+        bkgcube = ctools.ctbkgcube(obs)
+        #bkgcube['inobs']    = self._events
+        #bkgcube['inmodel']  = self._model
+        bkgcube['incube']   = 'NONE'
+        bkgcube['caldb']    = self._caldb
+        bkgcube['irf']      = self._irf
+        bkgcube['ebinalg']  = 'LOG'
+        bkgcube['emin']     = 0.1
+        bkgcube['emax']     = 100
+        bkgcube['enumbins'] = 20
+        bkgcube['nxpix']    = 10
+        bkgcube['nypix']    = 10
+        bkgcube['binsz']    = 0.4
+        bkgcube['coordsys'] = 'CEL'
+        bkgcube['proj']     = 'CAR'
+        bkgcube['xref']     = 83.63
+        bkgcube['yref']     = 22.01
+        bkgcube['outcube']  = 'ctbkgcube_py4.fits'
+        bkgcube['outmodel'] = 'ctbkgcube_py4.xml'
+        bkgcube['logfile']  = 'ctbkgcube_py4.log'
+        bkgcube['chatter']  = 4
+
+        # Execute ctbkgcube tool
+        bkgcube.logFileOpen()   # Make sure we get a log file
+        bkgcube.execute()
+
+        # Check result file
+        self._check_result_file('ctbkgcube_py4.fits')
+
         # Return
         return
@@ -138,13 +138,33 @@ class Test(test):
         butterfly['logfile'] = 'ctbutterfly_py2.log'
         butterfly['chatter'] = 3
-        # Run ctbutterfly tool
+        # Execute ctbutterfly tool
         butterfly.logFileOpen()   # Make sure we get a log file
         butterfly.execute()
         # Check result file
         self._check_result_file('ctbutterfly_py2.dat')
+        # Recover observation container
+        obs = butterfly.obs()
+
+        # TODO: Do someting test on observation container
+
+        # Copy ctbkgcube tool and execute copy
+        cpy_butterfly = butterfly
+        cpy_butterfly['outfile'] = 'ctbutterfly_py3.dat'
+        cpy_butterfly['logfile'] = 'ctbutterfly_py3.log'
+        cpy_butterfly['chatter'] = 4
+        cpy_butterfly.execute()
+
+        # Check result file
+        self._check_result_file('ctbutterfly_py3.dat')
+
+        # Clear ctbkgcube tool
+        butterfly.clear()
+
+        # TODO: Do some test after clearing
+
         # Return
         return
@@ -139,6 +139,7 @@ class Test(test):
         mask['outcube'] = 'ctcubemask_py2.fits'
         mask['logfile'] = 'ctcubemask_py2.log'
         mask['chatter'] = 3
+        mask['publish'] = True
         # Execute ctcubemask tool
         mask.logFileOpen()   # Make sure we get a log file
@@ -147,6 +148,21 @@ class Test(test):
         # Check result file
         self._check_result_file('ctcubemask_py2.fits', events=5542)
+        # Copy ctcubemask tool and execute copy
+        cpy_mask = mask
+        cpy_mask['outcube'] = 'ctcubemask_py3.fits'
+        cpy_mask['logfile'] = 'ctcubemask_py3.log'
+        cpy_mask['chatter']  = 4
+        cpy_mask.execute()
+
+        # Check result file
+        self._check_result_file('ctcubemask_py3.fits', events=5542)
+
+        # Clear ctcubemask tool
+        mask.clear()
+
+        # TODO: Do some test after clearing
+
         # Return
         return