diff --git a/Examples/fit/scatter2d/Honeycomb_fit.py b/Examples/fit/scatter2d/Honeycomb_fit.py
index 3bc0f5961d4333c8e8c019bfd974dde29182ca7f..3f01274449a186a73a51e7949eb33847bcfd36fc 100755
--- a/Examples/fit/scatter2d/Honeycomb_fit.py
+++ b/Examples/fit/scatter2d/Honeycomb_fit.py
@@ -12,7 +12,7 @@ import numpy
 import matplotlib.pyplot as plt
 from scipy.optimize import differential_evolution
 import bornagain as ba
-from bornagain import angstrom
+from bornagain import angstrom, sample_tools as st
 
 # number of points on which the computed result is plotted
 scan_size = 1500
@@ -158,12 +158,12 @@ def plot_sld_profile(fitParams, filename):
     plt.figure()
     parameters = dict(fitParams, **fixedParams)
 
-    z_300_p, sld_300_p = ba.materialProfile(get_sample(parameters, 1))
-    z_300_m, sld_300_m = ba.materialProfile(get_sample(parameters, -1))
+    z_300_p, sld_300_p = st.materialProfile(get_sample(parameters, 1))
+    z_300_m, sld_300_m = st.materialProfile(get_sample(parameters, -1))
 
-    z_150_p, sld_150_p = ba.materialProfile(
+    z_150_p, sld_150_p = st.materialProfile(
         get_sample(parameters, 1, ms150=parameters["ms150"]))
-    z_150_m, sld_150_m = ba.materialProfile(
+    z_150_m, sld_150_m = st.materialProfile(
         get_sample(parameters, -1, ms150=parameters["ms150"]))
 
     plt.figure()
@@ -179,7 +179,7 @@ def plot_sld_profile(fitParams, filename):
     plt.legend()
     plt.tight_layout()
     plt.savefig(filename)
-    plt.close()
+#    plt.close()
 
 
 ####################################################################