diff --git a/Examples/scatter2d/ConstantBackground.py b/Examples/scatter2d/ConstantBackground.py
index a63dcf6a3cccc86f39880169b755c91647ee9b08..ca0f125643c542ae093b75b3ae9ae5c37cf1aed6 100755
--- a/Examples/scatter2d/ConstantBackground.py
+++ b/Examples/scatter2d/ConstantBackground.py
@@ -7,37 +7,8 @@ from bornagain import deg, nm
 
 
 def get_sample():
-    """
-    Returns a sample with uncorrelated cylinders on a substrate.
-    """
-
-    # Define materials
-    material_Particle = ba.HomogeneousMaterial("Particle", 0.0006, 2e-08)
-    material_Substrate = ba.HomogeneousMaterial("Substrate", 6e-06, 2e-08)
-    material_Vacuum = ba.HomogeneousMaterial("Vacuum", 0, 0)
-
-    # Define form factors
-    ff = ba.FormFactorCylinder(5*nm, 5*nm)
-
-    # Define particles
-    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/Examples/scatter2d/DetectorResolutionFunction.py b/Examples/scatter2d/DetectorResolutionFunction.py
index 1f40c554d49547a57946c21dbee0c8c595800ba0..4c1eeea9e38096390b6f61429a422bd4be330569 100755
--- a/Examples/scatter2d/DetectorResolutionFunction.py
+++ b/Examples/scatter2d/DetectorResolutionFunction.py
@@ -7,37 +7,8 @@ from bornagain import deg, nm
 
 
 def get_sample():
-    """
-    Returns a sample with uncorrelated cylinders on a substrate.
-    """
-
-    # Define materials
-    material_Particle = ba.HomogeneousMaterial("Particle", 0.0006, 2e-08)
-    material_Substrate = ba.HomogeneousMaterial("Substrate", 6e-06, 2e-08)
-    material_Vacuum = ba.HomogeneousMaterial("Vacuum", 0, 0)
-
-    # Define form factors
-    ff = ba.FormFactorCylinder(5*nm, 5*nm)
-
-    # Define particles
-    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):