diff --git a/cscripts/csbkgmodel.py b/cscripts/csbkgmodel.py
index e5009a6..2874724 100755
--- a/cscripts/csbkgmodel.py
+++ b/cscripts/csbkgmodel.py
@@ -355,6 +355,9 @@ class csbkgmodel(ctools.csobservation):
         like = ctools.ctlike(obs)
         like.run()
 
+        # Extract optimiser
+        opt = like.opt()
+
         # Extract fitted model
         model = like.obs().models()[0].copy()
 
@@ -376,6 +379,9 @@ class csbkgmodel(ctools.csobservation):
             like = ctools.ctlike(obs)
             like.run()
 
+            # Extract optimiser
+            opt = like.opt()
+
             # Extract fitted model
             model = like.obs().models()[0].copy()
 
@@ -388,6 +394,9 @@ class csbkgmodel(ctools.csobservation):
                 if spectral[iint].error() == 0.0:
                     spectral.remove(nodes - i - 1)
 
+        # Write optimizer
+        self._log_string(gammalib.EXPLICIT, str(opt))
+
         # Return model
         return model