diff --git a/Examples/scatter2d/Cylinders.py b/Examples/scatter2d/Cylinders.py
index 720e63951b676f852181d18d6bfe1475c8543099..f698de647b603fac8c34c1a8899024c85b6a2ec2 100755
--- a/Examples/scatter2d/Cylinders.py
+++ b/Examples/scatter2d/Cylinders.py
@@ -8,31 +8,8 @@ from bornagain import deg, nm
 
 
 def get_sample():
-    # Define materials
-    material_Particle = ba.HomogeneousMaterial("Particle", 6e-4, 2e-08)
-    material_Substrate = ba.HomogeneousMaterial("Substrate", 6e-6, 2e-08)
-    material_Vacuum = ba.HomogeneousMaterial("Vacuum", 0, 0)
-
-    # Define particles
-    ff = ba.FormFactorCylinder(5*nm, 5*nm)
-    particle = ba.Particle(material_Particle, ff)
-
-    # Define particle layouts
-    layout = ba.ParticleLayout()
-    layout.addParticle(particle)
-    layout.setTotalParticleSurfaceDensity(0.01)
-
-    # Define layers
-    layer_1 = ba.Layer(material_Vacuum)
-    layer_1.addLayout(layout)
-    layer_2 = ba.Layer(material_Substrate)
-
-    # Define sample
-    sample = ba.MultiLayer()
-    sample.addLayer(layer_1)
-    sample.addLayer(layer_2)
-
-    return sample
+    from bornagain import std_samples
+    return std_samples.cylinders()
 
 
 def get_simulation(sample):
diff --git a/Wrap/Python/std_samples.py b/Wrap/Python/std_samples.py
index 68a005db8986fc512a988f8920e725c328de816a..61b135c848acf2e6cb2676e6342b4e3f1bdd3810 100644
--- a/Wrap/Python/std_samples.py
+++ b/Wrap/Python/std_samples.py
@@ -2,11 +2,11 @@
 BornAgain collection of standard sample.
 """
 import bornagain as ba
-from bornagain import deg, angstrom
+from bornagain import deg, nm
 
 def alternating_layers():
     """
-    Consists of 20 alternating Ti and Ni layers.
+    Returns multilayer sample consisting of 20 alternating Ti and Ni layers.
     """
 
     # Define materials
@@ -17,8 +17,8 @@ def alternating_layers():
 
     # Define layers
     ambient_layer = ba.Layer(m_ambient)
-    ti_layer = ba.Layer(m_ti, 30*angstrom)
-    ni_layer = ba.Layer(m_ni, 70*angstrom)
+    ti_layer = ba.Layer(m_ti, 3*nm)
+    ni_layer = ba.Layer(m_ni, 7*nm)
     substrate_layer = ba.Layer(m_substrate)
 
     # Define sample
@@ -30,3 +30,34 @@ def alternating_layers():
     sample.addLayer(substrate_layer)
 
     return sample
+
+
+def cylinders():
+    """
+    Returns sample consisting of dilute cylinders on substrate.
+    """
+    # Define materials
+    material_Particle = ba.HomogeneousMaterial("Particle", 6e-4, 2e-08)
+    material_Substrate = ba.HomogeneousMaterial("Substrate", 6e-6, 2e-08)
+    material_Vacuum = ba.HomogeneousMaterial("Vacuum", 0, 0)
+
+    # Define particles
+    ff = ba.FormFactorCylinder(5*nm, 5*nm)
+    particle = ba.Particle(material_Particle, ff)
+
+    # Define particle layouts
+    layout = ba.ParticleLayout()
+    layout.addParticle(particle)
+    layout.setTotalParticleSurfaceDensity(0.01)
+
+    # Define layers
+    layer_1 = ba.Layer(material_Vacuum)
+    layer_1.addLayout(layout)
+    layer_2 = ba.Layer(material_Substrate)
+
+    # Define sample
+    sample = ba.MultiLayer()
+    sample.addLayer(layer_1)
+    sample.addLayer(layer_2)
+
+    return sample
diff --git a/auto/Wrap/libBornAgainSim_wrap.cpp b/auto/Wrap/libBornAgainSim_wrap.cpp
index d56e0967d649eb352ad5a0280031e46e93880da6..0632dedf53d40c114639909125c22e28e1d8e202 100644
--- a/auto/Wrap/libBornAgainSim_wrap.cpp
+++ b/auto/Wrap/libBornAgainSim_wrap.cpp
@@ -5292,7 +5292,7 @@ SWIG_AsVal_std_complex_Sl_double_Sg_  (PyObject *o, std::complex<double>* val)
 
 
 SWIGINTERNINLINE PyObject*
-SWIG_From_std_complex_Sl_double_Sg_  (/*@SWIG:/usr/local/share/swig/4.0.2/typemaps/swigmacros.swg,104,%ifcplusplus@*/
+SWIG_From_std_complex_Sl_double_Sg_  (/*@SWIG:/usr/share/swig4.0/typemaps/swigmacros.swg,104,%ifcplusplus@*/
 
 const std::complex<double>&
 
diff --git a/hugo/content/py/simulation/gisas/_index.md b/hugo/content/py/simulation/gisas/_index.md
index 34ece8da5cb47e29e8d00eec18bc3820a336f2a3..dafbca65890c0e76589a49cc8322330e7c798fb8 100644
--- a/hugo/content/py/simulation/gisas/_index.md
+++ b/hugo/content/py/simulation/gisas/_index.md
@@ -21,7 +21,7 @@ we consider a dilute random assembly of monodisperse cylindrical disks on a subs
 {{< highlightfile file="Examples/scatter2d/Cylinders.py" >}}
 <p>
 
-### Explanation
+### TODO MOVE Explanation
 
 The ***Define materials*** stance defines three materials in terms
 of their refractive index.