diff --git a/script/staff/Frederik/230506_xps-selPeaks_VB.py b/script/staff/Frederik/230506_xps-selPeaks_VB.py
new file mode 100644
index 00000000..07acd64d
--- /dev/null
+++ b/script/staff/Frederik/230506_xps-selPeaks_VB.py
@@ -0,0 +1,66 @@
+import ch.psi.pshell.plot.RangeSelectionPlot.RangeSelection as RangeSelection
+
+# caution: do not insert multiple runs in the same script!
+#          earlier runs are overwritten!
+
+ranges = []
+
+# overview spec:
+#rs = RangeSelection(100.0,800.0)
+#rs.setVars([0.3, 0.2, 1.0, 800.0]) #Time, Size, Iteration, photon energy
+#ranges.append(rs)
+
+# Au 4f
+#rs = RangeSelection(105.0,114.0) # before 105 to 115 ----   -1eV
+#rs.setVars([0.4, 0.05, 3, 200., "Au 4f"]) #Time, Size, Iteration, photon energy
+#ranges.append(rs)
+
+# C 1s
+#rs = RangeSelection(100.0,108.0) # before 105 to 115 ----   -4eV 
+#rs.setVars([1., 0.04, 3, 394., "C 1s"]) #Time, Size, Iteration, photon energy
+#ranges.append(rs)
+
+# Au 4d5/2
+#rs = RangeSelection(100.0,114.0) # before 100 to 110 ----  +4eV, center was 445eV  
+#rs.setVars([0.8, 0.04, 3, 445., "Au 4d5/2"]) #Time, Size, Iteration, photon energy
+#ranges.append(rs)
+
+# N 1s
+#rs = RangeSelection(100.0,107.0) # before 105 to 115 ----  -5eV
+#rs.setVars([1, 0.04, 3, 507., "N 1s"]) #Time, Size, Iteration, photon energy
+#ranges.append(rs)
+
+# overview
+#rs = RangeSelection(100.0,800.0)
+#rs.setVars([0.3, 0.2, 1, 800., "overview"]) #Time, Size, Iteration, photon energy
+#ranges.append(rs)
+
+# Zn 2p
+#rs = RangeSelection(107.0,112.0) # before 105 to 115 ---  -5eV
+#rs.setVars([0.4, 0.05, 3, 1135., "Zn 2p"]) #Time, Size, Iteration, photon energy
+#ranges.append(rs)
+
+# -------------- VB peaks -----------------------------
+# Zn 3d
+rs = RangeSelection(90.0,106.0)
+rs.setVars([0.4, 0.05, 2, 110.,"Zn2p",50.]) #Time, Size, Iteration, photon energy
+ranges.append(rs)
+
+
+# Zn 3d Homo region near Fermi level
+rs = RangeSelection(103.5,106.)
+rs.setVars([0.5, 0.01, 5, 110.,"Homo", 20.]) #Time, Size, Iteration, photon energy
+ranges.append(rs)
+
+
+# Au 4d3/2
+#rs = RangeSelection(100.0,115.0)
+#rs.setVars([0.4, 0.05, 3, 465., "Au 4d3/2"]) #Time, Size, Iteration, photon energy
+#ranges.append(rs)
+
+
+
+run("XPSSpectrum", {"save_scienta_image":True, 
+                    "ranges":ranges, 
+                    "ENDSCAN":True, 
+                    "pass_energy":20})