added fix_spat_for_ts parameter to ctlike

Michael Mayer
1 parent e7db1c12
Showing 3 changed files with 24 additions and 0 deletions
src/ctlike/ctlike.cpp
src/ctlike/ctlike.hpp
src/ctlike/ctlike.par
@@ -262,6 +262,26 @@ void ctlike::run(void)
         double  logL_src = m_logL;
         GModels models   = m_obs.models();
  
+        // Fix spatial parameters if requested
+        if (m_fix_spat_for_ts) {
+            for(int i = 0; i < models.size(); ++i) {
+
+                // Continue only if model is skymodel
+                GModelSky* sky= dynamic_cast<GModelSky*>(models[i]);
+                if (sky != NULL) {
+
+                    // Fix spatial parameters
+                    GModelSpatial* spat = sky->spatial();
+                    for (int j = 0; j < spat->size(); j++) {
+                        (*spat)[j].fix();
+                    } // endfor: loop over spatial parameters
+
+                } // endif: there was a sky model
+
+            } // endfor: Loop over models
+
+        } // endif: spatial parameter should be fixed
+
         // Loop over stored models, remove source and refit
         for (int i = 0; i < ts_srcs.size(); ++i) {
             models.remove(ts_srcs[i]);
@@ -394,6 +414,7 @@ void ctlike::get_parameters(void)
     // Get other parameters
     m_refit       = (*this)["refit"].boolean();
     m_apply_edisp = (*this)["edisp"].boolean();
+    m_fix_spat_for_ts = (*this)["fix_spat_for_ts"].boolean();
  
     // Optionally read ahead parameters so that they get correctly
     // dumped into the log file
@@ -76,6 +76,7 @@ protected:
     bool          m_refit;       //!< Refitting
     std::string   m_outmodel;    //!< Source model output XML file
     bool          m_apply_edisp; //!< Apply energy dispersion?
+    bool          m_fix_spat_for_ts; //!< Fix spatial parameters for TS computation?
  
     // Members
     GObservations m_obs;        //!< Observations
@@ -28,6 +28,8 @@ outmodel, f, a, &quot;crab_results.xml&quot;,,, &quot;Source model output file&quot;
 stat,     s, h, "POISSON",,, "Optimization statistics"
 edisp,    b, h, no,,, "Apply energy dispersion?"
 refit,    b, h, no,,, "Do refitting?"
+fix_spat_for_ts,b,h,no,,,"Fix spatial parameters for TS computation?"
+
 #
 # Standard parameters
 #====================