diff --git a/script/users/heinrich/xps-GeTeCu-Ek120.py b/script/users/heinrich/xps-GeTeCu-Ek120.py
new file mode 100644
index 00000000..ba571a12
--- /dev/null
+++ b/script/users/heinrich/xps-GeTeCu-Ek120.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 = []
+
+# Cu 2p
+rs = RangeSelection(111.0,129.0)
+rs.setVars([0.2, 0.1, 1, 199.]) #Time, Size, Iteration, photon energy
+ranges.append(rs)
+
+# Cu 3s (bracketed by Ge 3p)
+rs = RangeSelection(111.0,129.0)
+rs.setVars([0.2, 0.1, 2, 247.]) #Time, Size, Iteration, photon energy
+ranges.append(rs)
+
+# Ge 3s, Te 4s
+rs = RangeSelection(109.0,131.0)
+rs.setVars([0.2, 0.1, 3, 299.]) #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":False, 
+                    "ranges":ranges, 
+                    "ENDSCAN":False, 
+                    "pass_energy":50})