diff --git a/Tests/Examples/CMakeLists.txt b/Tests/Examples/CMakeLists.txt
index d719da12911582ec7027a8e4ff3c82c25f56711b..59fc43eb40f90552aba8de1856d4d5b70fe36e0a 100644
--- a/Tests/Examples/CMakeLists.txt
+++ b/Tests/Examples/CMakeLists.txt
@@ -169,7 +169,6 @@ test_example(scatter2d/Cylinders 2e-10)
 test_example(scatter2d/CylindersAndPrisms 2e-10)
 test_example(scatter2d/CylindersInAverageLayer 2e-10)
 test_example(scatter2d/CylindersInBA 2e-10)
-test_example(scatter2d/DetectorResolutionFunction 2e-10)
 test_example(scatter2d/DodecahedraSAS 2e-10)
 run_example(scatter2d/FindPeaks)
 test_example(scatter2d/HalfSpheresInAverageTopLayer 2e-10)
@@ -195,6 +194,7 @@ test_example(scatter2d/PolarizedSANS 2e-10)
 test_example(scatter2d/PolydisperseCylinders 2e-10)
 run_example(scatter2d/PositionVariance) # TODO -> test
 test_example(scatter2d/RectangularGrating 0.6)
+test_example(scatter2d/Resolution 2e-10)
 test_example(scatter2d/RotatedPyramids 2e-10)
 test_example(scatter2d/SlicedLayer 2e-10)
 test_example(scatter2d/SpheresAtHexLattice 2e-10)
diff --git a/Tests/ReferenceData/MiniExamples/scatter2d/DetectorResolutionFunction.int b/Tests/ReferenceData/MiniExamples/scatter2d/DetectorResolutionFunction.int
deleted file mode 100644
index bbffa0e214de55fc3581d1f1d70ceeb2dcf41cd8..0000000000000000000000000000000000000000
--- a/Tests/ReferenceData/MiniExamples/scatter2d/DetectorResolutionFunction.int
+++ /dev/null
@@ -1,21 +0,0 @@
-# BornAgain Intensity Data
-
-# axis-0
-EquiDivision("q_y (1/nm)", 11, 0, 70)
-
-# axis-1
-EquiDivision("q_z (1/nm)", 11, 0, 70)
-
-# data
-2.065505985016e+09 1.148948748637e+09 6.601376497852e+08 3.721696797312e+08 1.585671987302e+08 4.197134791946e+07 1.289198220746e+07 2.865319996493e+07 4.598194583070e+07 4.447976524630e+07 2.769658991042e+07
-1.217702777557e+09 6.773780363926e+08 3.891965929683e+08 2.194246432607e+08 9.349135440354e+07 2.474733128686e+07 7.600542402230e+06 1.689106040087e+07 2.710531780314e+07 2.621752284391e+07 1.632281530234e+07
-3.484354062550e+08 1.938409367035e+08 1.113768849536e+08 6.279634144931e+07 2.675765093868e+07 7.083265790153e+06 2.174863591412e+06 4.832032667766e+06 7.752474323546e+06 7.495820716745e+06 4.664265227857e+06
-1.369313836173e+07 7.619155457171e+06 4.378471485110e+06 2.469097463956e+06 1.052184043691e+06 2.785002814609e+05 8.542758785081e+04 1.895157621633e+05 3.034095103687e+05 2.923855686831e+05 1.810550917696e+05
-2.308793655674e+07 1.284611930480e+07 7.380370962689e+06 4.161288178566e+06 1.773576663714e+06 4.698070612161e+05 1.442694577481e+05 3.208918888425e+05 5.162123681152e+05 5.010415854464e+05 3.134579577070e+05
-3.394611125214e+07 1.889163608744e+07 1.085479192737e+07 6.121370485228e+06 2.609506407503e+06 6.913676521324e+05 2.121590324579e+05 4.715436883492e+05 7.580457050440e+05 7.348724767930e+05 4.588989247138e+05
-4.995738078622e+06 2.781035111778e+06 1.598260184390e+06 9.016222395545e+05 3.845258576963e+05 1.019290448005e+05 3.126522149421e+04 6.946549503211e+04 1.116608739720e+05 1.081911382341e+05 6.749286786751e+04
-2.811895665040e+06 1.565253872606e+06 8.991080576758e+05 5.069220769891e+05 2.160738430689e+05 5.725536651333e+04 1.754269782408e+04 3.896232985762e+04 6.268754021998e+04 6.085086223777e+04 3.807837182841e+04
-9.075616107462e+06 5.053898388695e+06 2.903787801212e+06 1.637904492382e+06 6.985662785919e+05 1.852413635916e+05 5.673052793807e+04 1.259657877362e+05 2.027136440449e+05 1.967399182308e+05 1.230383402155e+05
-2.580581554380e+06 1.437611885452e+06 8.261999719759e+05 4.662363053822e+05 1.989769785858e+05 5.280792664059e+04 1.616598335301e+04 3.589394104737e+04 5.780331533728e+04 5.612863295175e+04 3.511450074651e+04
-4.709658232181e+05 2.623332798797e+05 1.506339262356e+05 8.491105916443e+04 3.619168102268e+04 9.592884502694e+03 2.928663297868e+03 6.488525652257e+03 1.043903442693e+04 1.013162934169e+04 6.340468232129e+03
-
diff --git a/Tests/ReferenceData/MiniExamples/scatter2d/Resolution.0.int b/Tests/ReferenceData/MiniExamples/scatter2d/Resolution.0.int
new file mode 100644
index 0000000000000000000000000000000000000000..aef7b5e934e860575571a128f52ce8235989ce37
--- /dev/null
+++ b/Tests/ReferenceData/MiniExamples/scatter2d/Resolution.0.int
@@ -0,0 +1,17 @@
+# BornAgain Intensity Data
+
+# axis-0
+EquiDivision("u (mm)", 7, 0, 107)
+
+# axis-1
+EquiDivision("v (mm)", 11, 0, 140)
+
+# data
+7.983153998380e+00 4.224545423750e-02 2.720642290345e-02 2.098662693605e-02 1.552710889535e-02 1.034283868494e-02 5.875571105760e-03 2.631691406116e-03 8.001981567807e-04 1.094220314020e-04 1.646882267356e-07
+1.192629797929e+02 6.265058077342e-01 4.032907865732e-01 3.110887989017e-01 2.303204578129e-01 1.537551678134e-01 8.780205316532e-02 3.977493018966e-02 1.239464737081e-02 1.802356266143e-03 3.012268272074e-06
+3.260799866421e+02 1.708142224692e+00 1.102750183971e+00 8.546304263366e-01 6.368922826433e-01 4.289021818409e-01 2.478025989900e-01 1.140870384229e-01 3.642847205663e-02 5.541658982968e-03 1.157895139460e-05
+2.675244116301e+04 1.400955386371e+02 9.064270626928e+01 7.048779067597e+01 5.277314387636e+01 3.575154665996e+01 2.081186602990e+01 9.673900453759e+00 3.128691099537e+00 4.855574396350e-01 1.114682319120e-03
+3.260799866421e+02 1.708142224692e+00 1.102750183971e+00 8.546304263366e-01 6.368922826433e-01 4.289021818409e-01 2.478025989900e-01 1.140870384229e-01 3.642847205663e-02 5.541658982968e-03 1.157895139460e-05
+1.192629797929e+02 6.265058077342e-01 4.032907865732e-01 3.110887989017e-01 2.303204578129e-01 1.537551678134e-01 8.780205316532e-02 3.977493018966e-02 1.239464737081e-02 1.802356266143e-03 3.012268272074e-06
+7.983153998380e+00 4.224545423750e-02 2.720642290345e-02 2.098662693605e-02 1.552710889535e-02 1.034283868494e-02 5.875571105759e-03 2.631691406116e-03 8.001981567806e-04 1.094220314020e-04 1.646882267356e-07
+
diff --git a/Tests/ReferenceData/MiniExamples/scatter2d/Resolution.1.int b/Tests/ReferenceData/MiniExamples/scatter2d/Resolution.1.int
new file mode 100644
index 0000000000000000000000000000000000000000..15412d72378e81afa2feba8785d5b46ae5fce4f2
--- /dev/null
+++ b/Tests/ReferenceData/MiniExamples/scatter2d/Resolution.1.int
@@ -0,0 +1,17 @@
+# BornAgain Intensity Data
+
+# axis-0
+EquiDivision("u (mm)", 7, 0, 107)
+
+# axis-1
+EquiDivision("v (mm)", 11, 0, 140)
+
+# data
+6.974484032566e+00 3.695794633428e-02 2.377558900386e-02 1.833474147810e-02 7.548969752231e+00 4.862041428170e-02 3.061423693260e-02 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00
+5.263550910517e-01 1.865828894268e-02 1.153317262768e-02 6.876254772596e-03 1.041951268165e+02 5.487718916224e-01 3.528127235842e-01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00
+2.719613729975e-01 2.080400986611e+01 2.422577958646e-01 1.463816900161e-01 7.624151503307e+00 1.511727246888e-01 5.447205768039e-02 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00
+3.552020701315e-02 2.848949429459e+02 1.504418858488e+00 9.693734670764e-01 7.494610377175e-01 1.690859367489e+03 9.226493902518e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00
+5.943576734232e+00 2.515605165543e+01 3.474859506952e+00 2.333414125238e+00 1.368971405679e+00 2.337236728300e+04 1.226163475538e+02 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00
+7.923445881889e+01 6.158821458798e+01 6.670921354781e+01 3.134187628089e+01 1.825153348126e+01 1.698807555606e+03 1.162549826691e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00
+6.178385676186e+00 4.470269050410e+00 2.882146027174e+02 3.753478813397e+00 2.278703536009e+00 1.418377073991e+00 8.289809882657e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00
+
diff --git a/Tests/ReferenceData/MiniExamples/scatter2d/Resolution.2.int b/Tests/ReferenceData/MiniExamples/scatter2d/Resolution.2.int
new file mode 100644
index 0000000000000000000000000000000000000000..edcb68ac0edd575ba2f925e8e4d4b4ac588800b0
--- /dev/null
+++ b/Tests/ReferenceData/MiniExamples/scatter2d/Resolution.2.int
@@ -0,0 +1,17 @@
+# BornAgain Intensity Data
+
+# axis-0
+EquiDivision("u (mm)", 7, 0, 107)
+
+# axis-1
+EquiDivision("v (mm)", 11, 0, 140)
+
+# data
+5.561615318236e+00 7.397072318363e-01 2.553433652307e-02 7.831758660370e-01 6.022460945594e+00 8.079996043494e-01 2.983796324095e-02 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00
+4.213389716563e-01 6.945937287865e-02 1.882142563965e-02 1.058147235560e+01 8.308796803155e+01 1.104789979823e+01 3.438501176333e-01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00
+2.328141362689e+00 1.663060485577e+01 2.319841816403e+00 9.164214327270e-01 6.105792024198e+00 8.997821720978e-01 5.926301467346e-02 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00
+2.894245101939e+01 2.271849071559e+02 3.021132827403e+01 1.134034890376e+00 1.723019908151e+02 1.348433270824e+03 1.789581596581e+02 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00
+7.289674378258e+00 2.100252015999e+01 5.559469178334e+00 3.923829791540e+00 2.373402407000e+03 1.863767618207e+04 2.469777439795e+03 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00
+6.939584898096e+01 6.389269355778e+01 6.259762906006e+01 3.371689190496e+01 1.901428514764e+02 1.356788869878e+03 1.816777303155e+02 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00
+5.396525962190e+00 3.344054627795e+01 2.305091391788e+02 3.247371078293e+01 2.360688543948e+00 2.203136508359e+00 6.750131050555e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00
+
diff --git a/Wrap/Python/ba_plot.py b/Wrap/Python/ba_plot.py
index 932060e0831868aea9be58d3fcf8376f72333bf1..0a5991f99a45d5eae314b06ab9faa8159dfacc18 100644
--- a/Wrap/Python/ba_plot.py
+++ b/Wrap/Python/ba_plot.py
@@ -219,9 +219,12 @@ def plot_specular_curve(result):
 
 def matches_reference(result, fname, subname=""):
     if not tolerance or not reference:
+        print(f"{fname}: no tolerance or no reference")
         return True
     reffile = reference + subname + ".int"
-    return ba.dataMatchesFile(result, reffile, tolerance)
+    ok = ba.dataMatchesFile(result, reffile, tolerance)
+    print(f"{fname} vs {reffile}, tol={tolerance} => ok={ok}")
+    return ok
 
 #  **************************************************************************  #
 #  multiple frames in one plot
diff --git a/auto/Examples/scatter2d/DetectorResolutionFunction.py b/auto/Examples/scatter2d/DetectorResolutionFunction.py
deleted file mode 100755
index 1560c1982bd0df615e2268e8eadb08d067b9e0f5..0000000000000000000000000000000000000000
--- a/auto/Examples/scatter2d/DetectorResolutionFunction.py
+++ /dev/null
@@ -1,28 +0,0 @@
-#!/usr/bin/env python3
-"""
-Cylinder form factor in DWBA with detector resolution function applied
-"""
-import bornagain as ba
-from bornagain import ba_plot as bp, deg, nm
-
-
-def get_sample():
-    from bornagain import std_samples
-    return std_samples.cylinders()
-
-
-def get_simulation(sample):
-    beam = ba.Beam(1e9, 0.1*nm, 0.2*deg)
-    n = 100
-    detector = ba.FlatDetector(n, n, 70, 70, beam, ba.FlatDetector.X, 2000, -35, -35)
-    detector.setResolutionFunction(
-        ba.ResolutionFunction2DGaussian(0.02*deg, 0.02*deg))
-    return ba.ScatteringSimulation(beam, sample, detector)
-
-
-if __name__ == '__main__':
-    sample = get_sample()
-    simulation = get_simulation(sample)
-    result = simulation.simulate()
-    bp.plot_simulation_result(result)
-    bp.show_or_export()
diff --git a/auto/MiniExamples/scatter2d/DetectorResolutionFunction.py b/auto/MiniExamples/scatter2d/DetectorResolutionFunction.py
deleted file mode 100755
index f540d48bc2a64183ede3439d2ce9cc74706024b3..0000000000000000000000000000000000000000
--- a/auto/MiniExamples/scatter2d/DetectorResolutionFunction.py
+++ /dev/null
@@ -1,28 +0,0 @@
-#!/usr/bin/env python3
-"""
-Cylinder form factor in DWBA with detector resolution function applied
-"""
-import bornagain as ba
-from bornagain import ba_plot as bp, deg, nm
-
-
-def get_sample():
-    from bornagain import std_samples
-    return std_samples.cylinders()
-
-
-def get_simulation(sample):
-    beam = ba.Beam(1e9, 0.1*nm, 0.2*deg)
-    n = 11
-    detector = ba.FlatDetector(n, n, 70, 70, beam, ba.FlatDetector.X, 2000, -35, -35)
-    detector.setResolutionFunction(
-        ba.ResolutionFunction2DGaussian(0.02*deg, 0.02*deg))
-    return ba.ScatteringSimulation(beam, sample, detector)
-
-
-if __name__ == '__main__':
-    sample = get_sample()
-    simulation = get_simulation(sample)
-    result = simulation.simulate()
-    bp.plot_simulation_result(result)
-    bp.show_or_export()
diff --git a/rawEx/scatter2d/DetectorResolutionFunction.py b/rawEx/scatter2d/DetectorResolutionFunction.py
deleted file mode 100755
index fdbb269eb201563737e7ef7aca6ebc712153f0e0..0000000000000000000000000000000000000000
--- a/rawEx/scatter2d/DetectorResolutionFunction.py
+++ /dev/null
@@ -1,28 +0,0 @@
-#!/usr/bin/env python3
-"""
-Cylinder form factor in DWBA with detector resolution function applied
-"""
-import bornagain as ba
-from bornagain import ba_plot as bp, deg, nm
-
-
-def get_sample():
-    from bornagain import std_samples
-    return std_samples.cylinders()
-
-
-def get_simulation(sample):
-    beam = ba.Beam(1e9, 0.1*nm, 0.2*deg)
-    n = <%= sm ? 11 : 100 %>
-    detector = ba.FlatDetector(n, n, 70, 70, beam, ba.FlatDetector.X, 2000, -35, -35)
-    detector.setResolutionFunction(
-        ba.ResolutionFunction2DGaussian(0.02*deg, 0.02*deg))
-    return ba.ScatteringSimulation(beam, sample, detector)
-
-
-if __name__ == '__main__':
-    sample = get_sample()
-    simulation = get_simulation(sample)
-    result = simulation.simulate()
-    bp.plot_simulation_result(result)
-    bp.show_or_export()