diff --git a/pyfstat/mcmc_based_searches.py b/pyfstat/mcmc_based_searches.py
index 78a1af122445ed56720a66769c7d88d1d1d8cba8..83b8a281868012b1844570b23709a87047d0f9fd 100644
--- a/pyfstat/mcmc_based_searches.py
+++ b/pyfstat/mcmc_based_searches.py
@@ -1732,7 +1732,7 @@ class MCMCSearch(core.BaseSearchClass):
             logging.info("theta0 index: {}".format(self.theta0_idx))
         logging.info("Max twoF: {} with parameters:".format(max_twoF))
         for k in np.sort(list(max_twoFd.keys())):
-            print("  {:10s} = {:1.9e}".format(k, max_twoFd[k]))
+            logging.info("  {:10s} = {:1.9e}".format(k, max_twoFd[k]))
         logging.info("Median +/- std for production values")
         for k in np.sort(list(median_std_d.keys())):
             if "std" not in k: