diff --git a/script/staff/xps-multiple_peaks.py b/script/staff/xps-multiple_peaks.py
new file mode 100644
index 00000000..c3d80c79
--- /dev/null
+++ b/script/staff/xps-multiple_peaks.py
@@ -0,0 +1,36 @@
+import ch.psi.pshell.plot.RangeSelectionPlot.RangeSelection as RangeSelection
+
+# caution: do not insert multiple runs in the same script!
+#          earlier runs are overwritten!
+
+ranges = []
+
+# Zn 2p
+rs = RangeSelection(105.0,115.0)
+rs.setVars([0.4, 0.05, 3, 1135.]) #Time, Size, Iteration, photon energy
+ranges.append(rs)
+
+# N 1s
+rs = RangeSelection(105.0,115.0)
+rs.setVars([0.4, 0.05, 3, 515.]) #Time, Size, Iteration, photon energy
+ranges.append(rs)
+
+# C 1s
+rs = RangeSelection(105.0,115.0)
+rs.setVars([0.4, 0.05, 3, 400.]) #Time, Size, Iteration, photon energy
+ranges.append(rs)
+
+# C 1s
+#rs = RangeSelection(114.0,126.0)
+#rs.setVars([0.2, 0.1, 5, 409.]) #Time, Size, Iteration, photon energy
+#ranges.append(rs)
+
+# O 1s
+#rs = RangeSelection(115.0,125.0)
+#rs.setVars([0.2, 0.1, 5, 655.5]) #Time, Size, Iteration, photon energy
+#ranges.append(rs)
+
+run("XPSSpectrum", {"save_scienta_image":True, 
+                    "ranges":ranges, 
+                    "ENDSCAN":True, 
+                    "pass_energy":50})