From ca874e05a54cd295df18ea9e0436f9eabb1d2b94 Mon Sep 17 00:00:00 2001
From: Michael Mayer <michael.mayer@desy.de>
Date: Thu, 5 Nov 2015 16:23:03 +0100
Subject: [PATCH] added fix_spat_for_ts parameter to ctlike

---
 src/ctlike/ctlike.cpp | 21 +++++++++++++++++++++
 src/ctlike/ctlike.hpp |  1 +
 src/ctlike/ctlike.par |  2 ++
 3 files changed, 24 insertions(+), 0 deletions(-)

diff --git a/src/ctlike/ctlike.cpp b/src/ctlike/ctlike.cpp
index d34e80f..2766b9e 100644
--- a/src/ctlike/ctlike.cpp
+++ b/src/ctlike/ctlike.cpp
@@ -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
diff --git a/src/ctlike/ctlike.hpp b/src/ctlike/ctlike.hpp
index 88a7fb0..ee40cc7 100644
--- a/src/ctlike/ctlike.hpp
+++ b/src/ctlike/ctlike.hpp
@@ -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
diff --git a/src/ctlike/ctlike.par b/src/ctlike/ctlike.par
index 5ae14f0..9e869e1 100644
--- a/src/ctlike/ctlike.par
+++ b/src/ctlike/ctlike.par
@@ -28,6 +28,8 @@ outmodel, f, a, "crab_results.xml",,, "Source model output file"
 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
 #====================
--
libgit2 0.21.4