diff --git a/Core/Export/SampleToPython.cpp b/Core/Export/SampleToPython.cpp
index 07301b1735112f5822230d9b284648c32cbee560..3a3b544cdf005206be7fb396940b6f1431e41e36 100644
--- a/Core/Export/SampleToPython.cpp
+++ b/Core/Export/SampleToPython.cpp
@@ -155,10 +155,10 @@ const std::map<MATERIAL_TYPES, std::string> factory_names{
 std::string SampleToPython::defineMaterials() const {
     const auto themap = m_materials->materialMap();
     if (themap.empty())
-        return "# No Materials.\n\n";
+        return "";
     std::ostringstream result;
     result << std::setprecision(12);
-    result << indent() << "# Define Materials\n";
+    result << indent() << "# Define materials\n";
     std::set<std::string> visitedMaterials;
     for (auto it : themap) {
         const std::string& key = it.first;
diff --git a/Examples/Python/sim01_Particles/CylindersAndPrisms.py b/Examples/Python/sim01_Particles/CylindersAndPrisms.py
index d3e490a7b510606263511ef41f6a19a12da80012..cbab1b3b6d47c381c6f2adc6dc66d69f608dbf28 100644
--- a/Examples/Python/sim01_Particles/CylindersAndPrisms.py
+++ b/Examples/Python/sim01_Particles/CylindersAndPrisms.py
@@ -10,7 +10,7 @@ def get_sample():
     Returns a sample with uncorrelated cylinders and prisms on a substrate.
     """
 
-    # Define Materials
+    # 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, 0.0)
diff --git a/Examples/Python/sim01_Particles/CylindersInBA.py b/Examples/Python/sim01_Particles/CylindersInBA.py
index 7cd32576df32694e30353c7cdb45927b7b7cea62..26855e44a93ba0f63acf6bb3e266b453afc3050d 100644
--- a/Examples/Python/sim01_Particles/CylindersInBA.py
+++ b/Examples/Python/sim01_Particles/CylindersInBA.py
@@ -11,7 +11,7 @@ def get_sample():
     implying a simulation in plain Born approximation.
     """
 
-    # Define Materials
+    # Define materials
     material_Particle = ba.HomogeneousMaterial("Particle", 0.0006, 2e-08)
     material_Vacuum = ba.HomogeneousMaterial("Vacuum", 0.0, 0.0)
 
diff --git a/Examples/Python/sim01_Particles/CylindersInDWBA.py b/Examples/Python/sim01_Particles/CylindersInDWBA.py
index c4091f92df2d0e58cb31be1a69795526915e0f66..28d884bf62df1e399244fdd48aba86bc76b148d0 100644
--- a/Examples/Python/sim01_Particles/CylindersInDWBA.py
+++ b/Examples/Python/sim01_Particles/CylindersInDWBA.py
@@ -10,7 +10,7 @@ def get_sample():
     Returns a sample with cylinders on a substrate.
     """
 
-    # Define Materials
+    # 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, 0.0)
diff --git a/Examples/Python/sim01_Particles/CylindersWithSizeDistribution.py b/Examples/Python/sim01_Particles/CylindersWithSizeDistribution.py
index 61919fbc83ec29ce4e77be80a26beba590ca7852..adc1fd1decbcae25800c687551080346903e08d2 100644
--- a/Examples/Python/sim01_Particles/CylindersWithSizeDistribution.py
+++ b/Examples/Python/sim01_Particles/CylindersWithSizeDistribution.py
@@ -11,7 +11,7 @@ def get_sample():
     The cylinders have a Gaussian size distribution.
     """
 
-    # Define Materials
+    # Define materials
     material_Particle = ba.HomogeneousMaterial("Particle", 0.0006, 2e-08)
     material_Vacuum = ba.HomogeneousMaterial("Vacuum", 0.0, 0.0)
 
diff --git a/Examples/Python/sim01_Particles/RotatedPyramids.py b/Examples/Python/sim01_Particles/RotatedPyramids.py
index 01b6ceea543caf71b93b165f8bc214138a04cb06..0077becb7bfe0fbb2de1ba8ed717227810b81960 100644
--- a/Examples/Python/sim01_Particles/RotatedPyramids.py
+++ b/Examples/Python/sim01_Particles/RotatedPyramids.py
@@ -10,7 +10,7 @@ def get_sample():
     Returns a sample with rotated pyramids on top of a substrate.
     """
 
-    # Define Materials
+    # 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, 0.0)
diff --git a/Examples/Python/sim01_Particles/TwoTypesOfCylindersWithSizeDistribution.py b/Examples/Python/sim01_Particles/TwoTypesOfCylindersWithSizeDistribution.py
index 90f1a3135a3c0e0120c2fa30e622db60221b136f..4609d50824a4e035e5a8a960b09a00ec1564754b 100644
--- a/Examples/Python/sim01_Particles/TwoTypesOfCylindersWithSizeDistribution.py
+++ b/Examples/Python/sim01_Particles/TwoTypesOfCylindersWithSizeDistribution.py
@@ -11,7 +11,7 @@ def get_sample():
     The cylinders are a 95:5 mixture of two different size distributions.
     """
 
-    # Define Materials
+    # Define materials
     material_Particle = ba.HomogeneousMaterial("Particle", 0.0006, 2e-08)
     material_Vacuum = ba.HomogeneousMaterial("Vacuum", 0.0, 0.0)
 
diff --git a/Examples/Python/sim02_Complexes/BiMaterialCylinders.py b/Examples/Python/sim02_Complexes/BiMaterialCylinders.py
index 74b0d0abf673c87ef7f3a0c3f71aede1c3177c92..d2c4b45a9bb0b4c96391021b1a34f63c5fec72c6 100644
--- a/Examples/Python/sim02_Complexes/BiMaterialCylinders.py
+++ b/Examples/Python/sim02_Complexes/BiMaterialCylinders.py
@@ -35,7 +35,7 @@ def get_sample():
     Particle shifted down to cross interface.
     """
 
-    # Define Materials
+    # Define materials
     material_Ag = ba.HomogeneousMaterial("Ag", 1.245e-05, 5.419e-07)
     material_Substrate = ba.HomogeneousMaterial("Substrate", 3.212e-06,
                                                 3.244e-08)
diff --git a/Examples/Python/sim02_Complexes/CoreShellNanoparticles.py b/Examples/Python/sim02_Complexes/CoreShellNanoparticles.py
index 4465546e4784ebc39a788d5a473622b12e9ad089..0ff5326cf1f6ba7d591fdfa20063ec2540d3fde2 100644
--- a/Examples/Python/sim02_Complexes/CoreShellNanoparticles.py
+++ b/Examples/Python/sim02_Complexes/CoreShellNanoparticles.py
@@ -10,7 +10,7 @@ def get_sample():
     Returns a sample with box-shaped core-shell particles on a substrate.
     """
 
-    # Define Materials
+    # Define materials
     material_Core = ba.HomogeneousMaterial("Core", 6e-05, 2e-08)
     material_Shell = ba.HomogeneousMaterial("Shell", 0.0001, 2e-08)
     material_Vacuum = ba.HomogeneousMaterial("Vacuum", 0.0, 0.0)
diff --git a/Examples/Python/sim02_Complexes/HexagonalLatticesWithBasis.py b/Examples/Python/sim02_Complexes/HexagonalLatticesWithBasis.py
index f4c6f774a26adce960ef431c41edc51146473e30..860732d97d96c21cdd8a05304edf2b4bc84b9855 100644
--- a/Examples/Python/sim02_Complexes/HexagonalLatticesWithBasis.py
+++ b/Examples/Python/sim02_Complexes/HexagonalLatticesWithBasis.py
@@ -12,7 +12,7 @@ def get_sample():
     forming two hexagonal close packed layers.
     """
 
-    # Define Materials
+    # 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, 0.0)
diff --git a/Examples/Python/sim02_Complexes/MesoCrystal.py b/Examples/Python/sim02_Complexes/MesoCrystal.py
index 6dfb587d18048df529c6ed1bd6c736952554abdd..3bd203daadd04ab6d752e726b4b8dc616ddf2e31 100644
--- a/Examples/Python/sim02_Complexes/MesoCrystal.py
+++ b/Examples/Python/sim02_Complexes/MesoCrystal.py
@@ -10,7 +10,7 @@ def get_sample():
     Returns a sample with a cylindrically shaped mesocrystal on a substrate.
     """
 
-    # Define Materials
+    # 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, 0.0)
diff --git a/Examples/Python/sim02_Complexes/ParticlesCrossingInterface.py b/Examples/Python/sim02_Complexes/ParticlesCrossingInterface.py
index 314f0754b04dad36d6d817035268fb40e0e971ac..4065dbd16b0bebe16f00d52bd210fcf489f372e0 100644
--- a/Examples/Python/sim02_Complexes/ParticlesCrossingInterface.py
+++ b/Examples/Python/sim02_Complexes/ParticlesCrossingInterface.py
@@ -28,7 +28,7 @@ def get_sample():
     Returns a sample with uncorrelated cylinders and prisms on a substrate.
     """
 
-    # Define Materials
+    # 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, 0.0)
diff --git a/Examples/Python/sim03_Structures/ApproximationDA.py b/Examples/Python/sim03_Structures/ApproximationDA.py
index 28f80aa0a355b8248d4fa994471773080e3dd3e6..c072968b804ba59f5effe416b447e620af14c038 100644
--- a/Examples/Python/sim03_Structures/ApproximationDA.py
+++ b/Examples/Python/sim03_Structures/ApproximationDA.py
@@ -11,7 +11,7 @@ def get_sample():
     The cylinder positions are modelled in Decoupling Approximation.
     """
 
-    # Define Materials
+    # 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, 0.0)
diff --git a/Examples/Python/sim03_Structures/ApproximationLMA.py b/Examples/Python/sim03_Structures/ApproximationLMA.py
index c0903baf4cdca78d9f6916b5f5d1d31e0574df84..73c8eebb9b3947f23b9da1de47a4608d8b17e872 100644
--- a/Examples/Python/sim03_Structures/ApproximationLMA.py
+++ b/Examples/Python/sim03_Structures/ApproximationLMA.py
@@ -11,7 +11,7 @@ def get_sample():
     The cylinder positions are modelled in Local Monodisperse Approximation.
     """
 
-    # Define Materials
+    # 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, 0.0)
diff --git a/Examples/Python/sim03_Structures/ApproximationSSCA.py b/Examples/Python/sim03_Structures/ApproximationSSCA.py
index 94467eda54cb32a58de67c43adedfc968ca7d535..abc4a8cb83ba6dd9afedf41a54b0aa6ac54e53cf 100644
--- a/Examples/Python/sim03_Structures/ApproximationSSCA.py
+++ b/Examples/Python/sim03_Structures/ApproximationSSCA.py
@@ -11,7 +11,7 @@ def get_sample():
     The cylinder positions are modelled in Size-Spacing Coupling  Approximation.
     """
 
-    # Define Materials
+    # 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, 0.0)
diff --git a/Examples/Python/sim03_Structures/CosineRipplesAtRectLattice.py b/Examples/Python/sim03_Structures/CosineRipplesAtRectLattice.py
index eaab091abfd82c38665d93668c3a26604a2fd80f..b4ae226890224b83dc62439efe2f945e434d0c20 100644
--- a/Examples/Python/sim03_Structures/CosineRipplesAtRectLattice.py
+++ b/Examples/Python/sim03_Structures/CosineRipplesAtRectLattice.py
@@ -12,7 +12,7 @@ def get_sample():
     The structure is modelled as a 2D Lattice.
     """
 
-    # Define Materials
+    # 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, 0.0)
diff --git a/Examples/Python/sim03_Structures/Interference1DLattice.py b/Examples/Python/sim03_Structures/Interference1DLattice.py
index b8d65c90852da57bd5bad12a57e0ca5c00d61152..a23eb272e885653ac054d3856755eda4b986c84e 100644
--- a/Examples/Python/sim03_Structures/Interference1DLattice.py
+++ b/Examples/Python/sim03_Structures/Interference1DLattice.py
@@ -13,7 +13,7 @@ def get_sample():
     modelled by very long boxes forming a 1D lattice with Cauchy correlations.
     """
 
-    # Define Materials
+    # 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, 0.0)
diff --git a/Examples/Python/sim03_Structures/Interference1DRadialParaCrystal.py b/Examples/Python/sim03_Structures/Interference1DRadialParaCrystal.py
index 717b727f63c0473d7347a467ff84465f558460fe..12e550a5a155149517eeca0327417769bcc95125 100644
--- a/Examples/Python/sim03_Structures/Interference1DRadialParaCrystal.py
+++ b/Examples/Python/sim03_Structures/Interference1DRadialParaCrystal.py
@@ -13,7 +13,7 @@ def get_sample():
     Returns a sample with cylinders on a substrate that form a radial paracrystal.
     """
 
-    # Define Materials
+    # 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, 0.0)
diff --git a/Examples/Python/sim03_Structures/Interference2DCenteredSquareLattice.py b/Examples/Python/sim03_Structures/Interference2DCenteredSquareLattice.py
index 9e8b5117265bd8a3ad359b3ec32b4a0910779fe9..03049b3fb7b280223d910daf2450e1c90088194a 100644
--- a/Examples/Python/sim03_Structures/Interference2DCenteredSquareLattice.py
+++ b/Examples/Python/sim03_Structures/Interference2DCenteredSquareLattice.py
@@ -12,7 +12,7 @@ def get_sample():
     forming a 2D centered square lattice
     """
 
-    # Define Materials
+    # 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, 0.0)
diff --git a/Examples/Python/sim03_Structures/Interference2DParaCrystal.py b/Examples/Python/sim03_Structures/Interference2DParaCrystal.py
index f92d6a9492492be566b942ce3326f32490b82fe4..dfae0d5c546f0ba8c4e51376d145db064bdf33d3 100644
--- a/Examples/Python/sim03_Structures/Interference2DParaCrystal.py
+++ b/Examples/Python/sim03_Structures/Interference2DParaCrystal.py
@@ -10,7 +10,7 @@ def get_sample():
     Returns a sample with cylinders on a substrate, forming a 2D paracrystal
     """
 
-    # Define Materials
+    # 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, 0.0)
diff --git a/Examples/Python/sim03_Structures/Interference2DRotatedSquareLattice.py b/Examples/Python/sim03_Structures/Interference2DRotatedSquareLattice.py
index d6d588cf6be7359ae43612e3aba0274b7131f301..33ec3f79809db7ed9ef3666054966b6723863b70 100644
--- a/Examples/Python/sim03_Structures/Interference2DRotatedSquareLattice.py
+++ b/Examples/Python/sim03_Structures/Interference2DRotatedSquareLattice.py
@@ -11,7 +11,7 @@ def get_sample():
     Returns a sample with cylinders on a substrate, forming a rotated 2D lattice
     """
 
-    # Define Materials
+    # 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, 0.0)
diff --git a/Examples/Python/sim03_Structures/Interference2DSquareFiniteLattice.py b/Examples/Python/sim03_Structures/Interference2DSquareFiniteLattice.py
index 5360333a20dcd29a0517e149ddddacafafa69114..e018b063c317bb9935640c8dce3f10b1fd4ffff5 100644
--- a/Examples/Python/sim03_Structures/Interference2DSquareFiniteLattice.py
+++ b/Examples/Python/sim03_Structures/Interference2DSquareFiniteLattice.py
@@ -11,7 +11,7 @@ def get_sample():
     Returns a sample with cylinders on a substrate, forming a 2D square lattice.
     """
 
-    # Define Materials
+    # 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, 0.0)
diff --git a/Examples/Python/sim03_Structures/Interference2DSquareLattice.py b/Examples/Python/sim03_Structures/Interference2DSquareLattice.py
index 633c66ec50e0c97284f1879b0c6a0bc6e9879686..affd7715ec59e66c594c4acf5a22332d0da8d990 100644
--- a/Examples/Python/sim03_Structures/Interference2DSquareLattice.py
+++ b/Examples/Python/sim03_Structures/Interference2DSquareLattice.py
@@ -11,7 +11,7 @@ def get_sample():
     Returns a sample with cylinders on a substrate, forming a 2D square lattice.
     """
 
-    # Define Materials
+    # 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, 0.0)
diff --git a/Examples/Python/sim03_Structures/SpheresAtHexLattice.py b/Examples/Python/sim03_Structures/SpheresAtHexLattice.py
index 1c83edd25c117414bcee19727bdc1000488a69d6..7f3c49d042b4f7e05fa5b942c17900b3816cb40e 100644
--- a/Examples/Python/sim03_Structures/SpheresAtHexLattice.py
+++ b/Examples/Python/sim03_Structures/SpheresAtHexLattice.py
@@ -11,7 +11,7 @@ def get_sample():
     forming a hexagonal 2D lattice.
     """
 
-    # Define Materials
+    # 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, 0.0)
diff --git a/Examples/Python/sim03_Structures/TriangularRipple.py b/Examples/Python/sim03_Structures/TriangularRipple.py
index 580f2961dd61d56877f53b3e055cd46e36497cbf..2815287042a10a22f4e7cb04849b1132650fd8b9 100644
--- a/Examples/Python/sim03_Structures/TriangularRipple.py
+++ b/Examples/Python/sim03_Structures/TriangularRipple.py
@@ -12,7 +12,7 @@ def get_sample():
     forming a 1D Paracrystal.
     """
 
-    # Define Materials
+    # 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, 0.0)
diff --git a/Examples/Python/sim04_Multilayers/BuriedParticles.py b/Examples/Python/sim04_Multilayers/BuriedParticles.py
index 2042a3f8c757d56245a60f971f6a74d34b65baa1..d543571da20b84d28b480c3c3dcc2cec353c3d3c 100644
--- a/Examples/Python/sim04_Multilayers/BuriedParticles.py
+++ b/Examples/Python/sim04_Multilayers/BuriedParticles.py
@@ -11,7 +11,7 @@ def get_sample():
     between vacuum and substrate.
     """
 
-    # Define Materials
+    # Define materials
     material_IntermLayer = ba.HomogeneousMaterial("IntermLayer", 3.45e-06,
                                                   5.24e-09)
     material_Particle = ba.HomogeneousMaterial("Particle", 0.0, 0.0)
diff --git a/Examples/Python/sim04_Multilayers/HalfSpheresInAverageTopLayer.py b/Examples/Python/sim04_Multilayers/HalfSpheresInAverageTopLayer.py
index c108ed596f9ae21d80b9b95090fcc7399dcea355..a691fbab52d14a78ea9c97b1a90129aade4781e1 100644
--- a/Examples/Python/sim04_Multilayers/HalfSpheresInAverageTopLayer.py
+++ b/Examples/Python/sim04_Multilayers/HalfSpheresInAverageTopLayer.py
@@ -14,7 +14,7 @@ def get_sample():
     Returns a sample with cylinders on a substrate.
     """
 
-    # Define Materials
+    # Define materials
     material_Particle = ba.HomogeneousMaterial("Particle", 3e-05, 2e-08)
     material_Substrate = ba.HomogeneousMaterial("Substrate", 6e-06, 2e-08)
     material_Vacuum = ba.HomogeneousMaterial("Vacuum", 0.0, 0.0)
diff --git a/Examples/Python/sim05_Magnetism/MagneticSpheres.py b/Examples/Python/sim05_Magnetism/MagneticSpheres.py
index 598b2a2342913fd7042de05f42f7c648f188dc1f..1b6277ea0a801dfb4a98908c95612aabf6fd4a44 100644
--- a/Examples/Python/sim05_Magnetism/MagneticSpheres.py
+++ b/Examples/Python/sim05_Magnetism/MagneticSpheres.py
@@ -13,7 +13,7 @@ def get_sample():
     Returns a sample with magnetic spheres in the substrate.
     """
 
-    # Define Materials
+    # Define materials
     magnetic_field = kvector_t(0, 0, 10000000)
     material_Particle = ba.HomogeneousMaterial("Particle", 2e-05, 4e-07,
                                                magnetic_field)
diff --git a/Examples/Python/sim11_Device/AxesInDifferentUnits.py b/Examples/Python/sim11_Device/AxesInDifferentUnits.py
index 45050234bd891f5c5cdbc9edae939a8e36d95629..95eb46e47e4608d0c88c4b765ea8a801be03adf9 100644
--- a/Examples/Python/sim11_Device/AxesInDifferentUnits.py
+++ b/Examples/Python/sim11_Device/AxesInDifferentUnits.py
@@ -13,7 +13,7 @@ def get_sample():
     Returns a sample with uncorrelated cylinders on a substrate.
     """
 
-    # Define Materials
+    # 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, 0.0)
diff --git a/Examples/Python/sim11_Device/BeamDivergence.py b/Examples/Python/sim11_Device/BeamDivergence.py
index ab5b733e9a231f9f7049558245cf21fe032f22c5..9fdc00c06f875932142dd7e433719193725f30e0 100644
--- a/Examples/Python/sim11_Device/BeamDivergence.py
+++ b/Examples/Python/sim11_Device/BeamDivergence.py
@@ -10,7 +10,7 @@ def get_sample():
     Returns a sample with uncorrelated cylinders on a substrate.
     """
 
-    # Define Materials
+    # 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, 0.0)
diff --git a/Examples/Python/sim11_Device/ConstantBackground.py b/Examples/Python/sim11_Device/ConstantBackground.py
index 31a9948b0644630c87a450aefd046049b228b81e..75935b4b9c9b04a5bcd5f0a08c49f783e4219943 100644
--- a/Examples/Python/sim11_Device/ConstantBackground.py
+++ b/Examples/Python/sim11_Device/ConstantBackground.py
@@ -10,7 +10,7 @@ def get_sample():
     Returns a sample with uncorrelated cylinders on a substrate.
     """
 
-    # Define Materials
+    # 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, 0.0)
diff --git a/Examples/Python/sim11_Device/DetectorResolutionFunction.py b/Examples/Python/sim11_Device/DetectorResolutionFunction.py
index d8289bfe055af8cdba75e5f18cb0a22cf8198e56..5add835f7998d4aad0f4c544aee6734061adfb37 100644
--- a/Examples/Python/sim11_Device/DetectorResolutionFunction.py
+++ b/Examples/Python/sim11_Device/DetectorResolutionFunction.py
@@ -10,7 +10,7 @@ def get_sample():
     Returns a sample with uncorrelated cylinders on a substrate.
     """
 
-    # Define Materials
+    # 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, 0.0)
diff --git a/Examples/Python/sim11_Device/OffSpecularSimulation.py b/Examples/Python/sim11_Device/OffSpecularSimulation.py
index 8ad41f0794d08bf10a081a5bc7fdd875747a84e7..fae853ac2cf200101b68f245ea4b10bfbfcc603c 100644
--- a/Examples/Python/sim11_Device/OffSpecularSimulation.py
+++ b/Examples/Python/sim11_Device/OffSpecularSimulation.py
@@ -15,7 +15,7 @@ def get_sample():
     modelled by infinitely long boxes forming a 1D lattice.
     """
 
-    # Define Materials
+    # 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, 0.0)
diff --git a/Examples/Python/sim11_Device/RectangularDetector.py b/Examples/Python/sim11_Device/RectangularDetector.py
index 2e3ec81a23f3c23a9176b99dd2f546e11e276b46..cabd270157fa3e12a2b23ae0b7e9920dbe2968df 100644
--- a/Examples/Python/sim11_Device/RectangularDetector.py
+++ b/Examples/Python/sim11_Device/RectangularDetector.py
@@ -18,7 +18,7 @@ def get_sample():
     Returns a sample with cylindrical particles on a substrate.
     """
 
-    # Define Materials
+    # 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, 0.0)
diff --git a/Examples/Python/sim21_Reflectometry/BasicPolarizedReflectometry.py b/Examples/Python/sim21_Reflectometry/BasicPolarizedReflectometry.py
index 3ab08f8de5bef01547828eea51b1a2c6e4a678aa..a61b8f9da614369882137b1e201b4eba94cbb7da 100644
--- a/Examples/Python/sim21_Reflectometry/BasicPolarizedReflectometry.py
+++ b/Examples/Python/sim21_Reflectometry/BasicPolarizedReflectometry.py
@@ -14,7 +14,7 @@ def get_sample():
     Defines sample and returns it
     """
 
-    # Define Materials
+    # Define materials
     material_Ambient = ba.MaterialBySLD("Ambient", 0.0, 0.0)
     magnetic_field = kvector_t(0, 100000000, 0)
     material_Layer = ba.MaterialBySLD("Layer", 0.0001, 1e-08, magnetic_field)
diff --git a/Examples/Python/sim21_Reflectometry/PolarizedNoAnalyzer.py b/Examples/Python/sim21_Reflectometry/PolarizedNoAnalyzer.py
index b561f1fe641803792433d689b1937d280558c51e..3be7294ae9257c39797f09005a7d29121cd65613 100644
--- a/Examples/Python/sim21_Reflectometry/PolarizedNoAnalyzer.py
+++ b/Examples/Python/sim21_Reflectometry/PolarizedNoAnalyzer.py
@@ -14,7 +14,7 @@ def get_sample():
     Defines sample and returns it
     """
 
-    # Define Materials
+    # Define materials
     material_Ambient = ba.MaterialBySLD("Ambient", 0.0, 0.0)
     magnetic_field = kvector_t(50000000, 86602540.3784, 0)
     material_Layer = ba.MaterialBySLD("Layer", 0.0001, 1e-08, magnetic_field)
diff --git a/Examples/Python/sim21_Reflectometry/PolarizedSpinFlip.py b/Examples/Python/sim21_Reflectometry/PolarizedSpinFlip.py
index 4a6eceaedea789f02a46ad18faffdc2212b7610f..901087ebb03e3a7cdba54d300c2bdb7811e9a2e2 100644
--- a/Examples/Python/sim21_Reflectometry/PolarizedSpinFlip.py
+++ b/Examples/Python/sim21_Reflectometry/PolarizedSpinFlip.py
@@ -15,7 +15,7 @@ def get_sample():
     Defines sample and returns it
     """
 
-    # Define Materials
+    # Define materials
     material_Ambient = ba.MaterialBySLD("Ambient", 0.0, 0.0)
     magnetic_field = kvector_t(50000000, 86602540.3784, 0)
     material_Layer = ba.MaterialBySLD("Layer", 0.0001, 1e-08, magnetic_field)
diff --git a/Examples/Python/sim22_OffSpecular/BoxesWithSpecularPeak.py b/Examples/Python/sim22_OffSpecular/BoxesWithSpecularPeak.py
index d7c474fa36e58d2032d8789edaf752f84163b379..e5d58ff02734e888aacc75f51f14792d2d29278d 100644
--- a/Examples/Python/sim22_OffSpecular/BoxesWithSpecularPeak.py
+++ b/Examples/Python/sim22_OffSpecular/BoxesWithSpecularPeak.py
@@ -10,7 +10,7 @@ def get_sample():
     Returns a sample with boxes on a substrate.
     """
 
-    # Define Materials
+    # Define materials
     material_Particle = ba.HomogeneousMaterial("Particle", 3e-05, 2e-08)
     material_Substrate = ba.HomogeneousMaterial("Substrate", 6e-06, 2e-08)
     material_Vacuum = ba.HomogeneousMaterial("Vacuum", 0.0, 0.0)
diff --git a/Examples/Python/sim23_SAS/PolarizedSANS.py b/Examples/Python/sim23_SAS/PolarizedSANS.py
index 21d6c46a5c116ab170fd9286a2b491045d9ec70b..eff2352864e67fb3f72c99bb0de61adebcfe60f9 100644
--- a/Examples/Python/sim23_SAS/PolarizedSANS.py
+++ b/Examples/Python/sim23_SAS/PolarizedSANS.py
@@ -15,7 +15,7 @@ def get_sample():
     Returns a sample with a magnetic core-shell particle in a solvent.
     """
 
-    # Define Materials
+    # Define materials
     magnetic_field = kvector_t(0, 0, 10000000)
     material_Core = ba.HomogeneousMaterial("Core", 6e-06, 2e-08, magnetic_field)
     material_Shell = ba.HomogeneousMaterial("Shell", 1e-07, 2e-08)
diff --git a/Examples/Python/sim29_DepthProbe/DepthProbe.py b/Examples/Python/sim29_DepthProbe/DepthProbe.py
index c8c767e9db63a2aed1c5772f9b0d8fecb26e4989..22de45ab8dae4c9be84d3f70d54bbcf8cb0b22af 100644
--- a/Examples/Python/sim29_DepthProbe/DepthProbe.py
+++ b/Examples/Python/sim29_DepthProbe/DepthProbe.py
@@ -55,7 +55,7 @@ def get_sample():
     Constructs a sample with one resonating Ti/Pt layer
     """
 
-    # Define Materials
+    # Define materials
     material_D2O = ba.HomogeneousMaterial("D2O", 0.00010116, 1.809e-12)
     material_Pt = ba.HomogeneousMaterial("Pt", 0.00010117, 3.01822e-08)
     material_Si = ba.HomogeneousMaterial("Si", 3.3009e-05, 0.0)
diff --git a/Examples/Python/sim31_Parameterization/AccessingSimulationResults.py b/Examples/Python/sim31_Parameterization/AccessingSimulationResults.py
index 78d1c876d7eba49d564dd5ccd368c02c5d59e0fe..5dd69b16e2eebbd937a9ad7d7c6944ce0eb5346e 100644
--- a/Examples/Python/sim31_Parameterization/AccessingSimulationResults.py
+++ b/Examples/Python/sim31_Parameterization/AccessingSimulationResults.py
@@ -15,7 +15,7 @@ def get_sample():
     Returns a sample with uncorrelated cylinders on a substrate.
     """
 
-    # Define Materials
+    # 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, 0.0)
diff --git a/Examples/Python/sim31_Parameterization/SimulationParameters.py b/Examples/Python/sim31_Parameterization/SimulationParameters.py
index 407583f79a66d76413a14a0f54f56750342eeeb1..195210a3e01c8bc8fc9c285a56b841766fbd7b7e 100644
--- a/Examples/Python/sim31_Parameterization/SimulationParameters.py
+++ b/Examples/Python/sim31_Parameterization/SimulationParameters.py
@@ -16,7 +16,7 @@ def get_sample():
     Parameter set is fixed.
     """
 
-    # Define Materials
+    # 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, 0.0)