diff --git a/script/jan/hyst_Dy_singleenergy.py b/script/jan/hyst_Dy_singleenergy.py
new file mode 100644
index 0000000..fb214e5
--- /dev/null
+++ b/script/jan/hyst_Dy_singleenergy.py
@@ -0,0 +1,53 @@
+off_cplus = -6.5
+off_cminus = -2.5
+rrate = 2.0        # ramping speed of mag field
+e1 = 1290.2        # energy 1 (usually on edge)
+e2 = 1283.0        # energy 2 (pre edge)
+
+B1 = 4.0            # field span
+B2 = -B1
+nr_loops = 1        # number of loops to acquire
+
+
+open_valve()
+
+loop_count = 0
+while loop_count < nr_loops:
+
+    set_pol_cplus(offset = off_cplus)
+    set_file('Dy_hyst_plus')
+    hyst_cont(field_x, B1, B2, rrate, [e1])
+    
+    set_pol_cminus(offset = off_cminus)
+    set_file('Dy_hyst_minus')
+    hyst_cont(field_x, B2, B1, rrate, [e1])
+    
+    set_pol_cminus(offset = off_cminus)
+    set_file('Dy_hyst_minus')
+    hyst_cont(field_x, B1, B2, rrate, [e1])
+    
+    set_pol_cplus(offset = off_cplus)
+    set_file('Dy_hyst_plus')
+    hyst_cont(field_x, B2, B1, rrate, [e1])
+
+    set_pol_cplus(offset = off_cplus)
+    set_file('Dy_hyst_plus')
+    hyst_cont(field_x, B1, B2, rrate, [e2])
+    
+    set_pol_cminus(offset = off_cminus)
+    set_file('Dy_hyst_minus')
+    hyst_cont(field_x, B2, B1, rrate, [e2])
+    
+    set_pol_cminus(offset = off_cminus)
+    set_file('Dy_hyst_minus')
+    hyst_cont(field_x, B1, B2, rrate, [e2])
+    
+    set_pol_cplus(offset = off_cplus)
+    set_file('Dy_hyst_plus')
+    hyst_cont(field_x, B2, B1, rrate, [e2])
+
+    loop_count = loop_count + 1
+    print "Loop # ",loop_count," finished."
+
+close_valve()
+