diff --git a/script/users/TUM/HAT_MultiXPS_2020.py b/script/users/TUM/HAT_MultiXPS_2020.py
new file mode 100644
index 00000000..18fc5e53
--- /dev/null
+++ b/script/users/TUM/HAT_MultiXPS_2020.py
@@ -0,0 +1,60 @@
+import ch.psi.pshell.plot.RangeSelectionPlot.RangeSelection as RangeSelection
+
+ranges = []
+
+#Cu 2p
+rs = RangeSelection(18, 32)
+rs.setVars([0.25, 0.05, 10, 1100.]) #Time, Size, Iteration, photon energy
+ranges.append(rs)
+
+#Na 1s
+rs = RangeSelection(158, 172)
+rs.setVars([0.25, 0.05, 10, 1100.]) #Time, Size, Iteration, photon energy
+ranges.append(rs)
+
+# N 1s
+rs = RangeSelection(143.0, 157.0)
+rs.setVars([0.25, 0.05, 10, 550.]) #Time, Size, Iteration, photon energy
+ranges.append(rs)
+
+# Fermi edge for N 1s
+rs = RangeSelection(540., 550.)
+rs.setVars([0.25, 0.05, 3, 550.]) #Time, Size, Iteration, photon energy
+ranges.append(rs)
+
+# Cl 2p3/2
+rs = RangeSelection(141.2, 155.2)
+rs.setVars([0.25, 0.05, 10, 350.]) #Time, Size, Iteration, photon energy
+ranges.append(rs)
+
+# Fermi edge for Cl 2p3/2
+rs = RangeSelection(340., 350.)
+rs.setVars([0.25, 0.05, 3, 350.]) #Time, Size, Iteration, photon energy
+ranges.append(rs)
+
+
+
+# C 1s
+rs = RangeSelection(151.0, 165.0)
+rs.setVars([0.25, 0.05, 10, 450.]) #Time, Size, Iteration, photon energy
+ranges.append(rs)
+
+# Fermi edge for C 1s
+rs = RangeSelection(440.0, 4450.0)
+rs.setVars([0.25, 0.05, 3, 450.]) #Time, Size, Iteration, photon energy
+ranges.append(rs)
+
+# Valence
+rs = RangeSelection(46.0, 59.0)
+rs.setVars([0.25, 0.05, 30, 60.]) #Time, Size, Iteration, photon energy
+ranges.append(rs)
+
+
+# caution: do not insert multiple run commands in the same script!
+#          earlier data is overwritten!
+
+run("XPSSpectrum", {"save_scienta_image":True, 
+                    "ranges": ranges, 
+                    "ENDSCAN": True, 
+                    "pass_energy":10})
+