diff --git a/tests/test_SV2o2.py b/tests/test_SV2o2.py
new file mode 100644
index 0000000..18e804d
--- /dev/null
+++ b/tests/test_SV2o2.py
@@ -0,0 +1,39 @@
+import os.path
+import numpy as np
+from numpy.testing import assert_allclose
+import pytest
+
+from optifik.scheludko import thickness_from_scheludko
+from optifik.io import load_spectrum
+from optifik.analysis import smooth_intensities
+
+import yaml
+
+
+def load():
+    FOLDER = os.path.join('tests', 'spectraVictor2', 'order2')
+
+    yaml_file = os.path.join(FOLDER, 'known_value.yaml')
+    with open(yaml_file, "r") as yaml_file:
+        thickness_dict = yaml.safe_load(yaml_file)
+    data = [(os.path.join(FOLDER, fn), val) for fn, val in thickness_dict.items()]
+    return data
+
+
+@pytest.mark.parametrize("spectrum_path, expected", load())
+def test_SV2o2(spectrum_path, expected):
+    lambdas, raw_intensities = load_spectrum(spectrum_path, lambda_min=450)
+    smoothed_intensities = smooth_intensities(raw_intensities)
+
+    refractive_index =  1.324188 + 3102.060378 / (lambdas**2)
+    prominence = 0.01
+
+
+    thickness_scheludko = thickness_from_scheludko(lambdas, smoothed_intensities,
+                                                 refractive_index=refractive_index,
+                                                 min_peak_prominence=prominence,
+                                                 plot=False)
+    result = thickness_scheludko.thickness
+
+    assert_allclose(result, expected, rtol=1e-1)
+
diff --git a/tests/test_SV2o3.py b/tests/test_SV2o3.py
new file mode 100644
index 0000000..478b22d
--- /dev/null
+++ b/tests/test_SV2o3.py
@@ -0,0 +1,39 @@
+import os.path
+import numpy as np
+from numpy.testing import assert_allclose
+import pytest
+
+from optifik.scheludko import thickness_from_scheludko
+from optifik.io import load_spectrum
+from optifik.analysis import smooth_intensities
+
+import yaml
+
+
+def load():
+    FOLDER = os.path.join('tests', 'spectraVictor2', 'order3')
+
+    yaml_file = os.path.join(FOLDER, 'known_value.yaml')
+    with open(yaml_file, "r") as yaml_file:
+        thickness_dict = yaml.safe_load(yaml_file)
+    data = [(os.path.join(FOLDER, fn), val) for fn, val in thickness_dict.items()]
+    return data
+
+
+@pytest.mark.parametrize("spectrum_path, expected", load())
+def test_SV2o3(spectrum_path, expected):
+    lambdas, raw_intensities = load_spectrum(spectrum_path, lambda_min=450)
+    smoothed_intensities = smooth_intensities(raw_intensities)
+
+    refractive_index =  1.324188 + 3102.060378 / (lambdas**2)
+    prominence = 0.02
+
+
+    thickness_scheludko = thickness_from_scheludko(lambdas, smoothed_intensities,
+                                                 refractive_index=refractive_index,
+                                                 min_peak_prominence=prominence,
+                                                 plot=False)
+    result = thickness_scheludko.thickness
+
+    assert_allclose(result, expected, rtol=1e-1)
+
diff --git a/tests/test_SV2o4.py b/tests/test_SV2o4.py
new file mode 100644
index 0000000..caac1bc
--- /dev/null
+++ b/tests/test_SV2o4.py
@@ -0,0 +1,39 @@
+import os.path
+import numpy as np
+from numpy.testing import assert_allclose
+import pytest
+
+from optifik.scheludko import thickness_from_scheludko
+from optifik.io import load_spectrum
+from optifik.analysis import smooth_intensities
+
+import yaml
+
+
+def load():
+    FOLDER = os.path.join('tests', 'spectraVictor2', 'order4')
+
+    yaml_file = os.path.join(FOLDER, 'known_value.yaml')
+    with open(yaml_file, "r") as yaml_file:
+        thickness_dict = yaml.safe_load(yaml_file)
+    data = [(os.path.join(FOLDER, fn), val) for fn, val in thickness_dict.items()]
+    return data
+
+
+@pytest.mark.parametrize("spectrum_path, expected", load())
+def test_SV2o4(spectrum_path, expected):
+    lambdas, raw_intensities = load_spectrum(spectrum_path, lambda_min=450)
+    smoothed_intensities = smooth_intensities(raw_intensities)
+
+    refractive_index =  1.324188 + 3102.060378 / (lambdas**2)
+    prominence = 0.02
+
+
+    thickness_scheludko = thickness_from_scheludko(lambdas, smoothed_intensities,
+                                                 refractive_index=refractive_index,
+                                                 min_peak_prominence=prominence,
+                                                 plot=False)
+    result = thickness_scheludko.thickness
+
+    assert_allclose(result, expected, rtol=1e-1)
+