Reduce memory usage by disposing the event lists after usage.

Jürgen Knödlseder
1 parent 6c66c03a
Showing 5 changed files with 18 additions and 3 deletions
NEWS
src/ctbin/ctbin.cpp
src/ctbkgcube/ctbkgcube.cpp
src/ctmodel/ctmodel.cpp
src/ctobssim/ctobssim.cpp
@@ -85,6 +85,9 @@ is set to the map centre. This is only identical to the pointing in
 case of binning a single event list and using the usepnt=yes option
 (or specifying the pointing manually as the xref and yref parameters).
+Memory usage has been reduced by disposing the event lists after their
+usage.
+
 8. ctlike - Maximum likelihood fitting
 --------------------------------------
@@ -141,6 +144,9 @@ the use for 2D response functions in FITS files.
 ctmodel now also supports the handling of energy dispersion. Energy
 dispersion can be switched on or off using the hidden edisp parameter.
+Memory usage has been reduced by disposing the event lists after their
+usage.
+
 12. ctexpcube - CTA exposure cube generation
 --------------------------------------------
@@ -165,6 +171,9 @@ ctbkgcube computes the background cube for an observation definition
 file. The cube definition may be taken from a counts cube or may be
 specified by user parameters.
+Memory usage has been reduced by disposing the event lists after their
+usage.
+
 15. ctcubemask - CTA cube filtering
 -----------------------------------
@@ -263,6 +263,9 @@ void ctbin::run(void)
             // Fill the cube
             fill_cube(obs);
+            // Dispose events to free memory
+            obs->dispose_events();
+
         } // endif: CTA observation found
     } // endfor: looped over observations
@@ -322,6 +322,9 @@ void ctbkgcube::run(void)
             // Fill the cube
             fill_cube(obs);
+            // Dispose events to free memory
+            obs->dispose_events();
+
         } // endif: CTA observation found
     } // endfor: looped over observations
@@ -273,6 +273,9 @@ void ctmodel::run(void)
             // Fill the cube
             fill_cube(obs);
+            // Dispose events to free memory
+            obs->dispose_events();
+
         } // endif: CTA observation found
     } // endfor: looped over observations
@@ -481,8 +484,6 @@ void ctmodel::get_parameters(void)
         // Read observation definition filename and response parameters
         m_obsfile = (*this)["obsfile"].filename();
-        m_caldb   = (*this)["caldb"].string();
-        m_irf     = (*this)["irf"].string();
         // If no observation definition file has been specified then read all
         // parameters that are necessary to create an observation from scratch
@@ -323,7 +323,6 @@ void ctobssim::run(void)
         // is enabled, this loop will be parallelized.
         #pragma omp for
         for (int i = 0; i < m_obs.size(); ++i) {
-std::cout << "Observation " << i << ": " << get_current_rss() << std::endl;
             // Get pointer on CTA observation
             GCTAObservation* obs = dynamic_cast<GCTAObservation*>(m_obs[i]);