From 2742f513532b586dc9792e628f160490a0320087 Mon Sep 17 00:00:00 2001
From: panepucci <ezequiel.panepucci@psi.ch>
Date: Fri, 18 Dec 2015 09:39:48 +0100
Subject: [PATCH] Startup

---
 script/test2.py | 42 ++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 42 insertions(+)
 create mode 100644 script/test2.py

diff --git a/script/test2.py b/script/test2.py
new file mode 100644
index 0000000..d3cec62
--- /dev/null
+++ b/script/test2.py
@@ -0,0 +1,42 @@
+
+
+if !noref:
+    caput('X06DA-OP-MO:MONOBEAM', -1 )
+e = get_energy()
+if dx is None:
+    dx = 7.5e-2/e
+
+# reinforce theta1 and 2 setting (setting could be RDBD 0.0001 and RTRY 1 when seteq does not finish nicely)
+# put 'Retry deadband' to 0.00004 and 'Retries Max' to 5
+time.sleep(0.2)
+caput('X06DA-OP-MO1:ROX1.RDBD',0.00004)
+caput('X06DA-OP-MO1:ROX2.RDBD',0.00004)
+caput('X06DA-OP-MO1:ROX1.RTRY',5)
+caput('X06DA-OP-MO1:ROX2.RTRY',5)
+
+#time stamp
+print time.strftime('%X %x')
+
+#MVR,axis,-dx
+#SCAN,axis,0,2*dx,20,/fit,/centre,SILENT = silent,FFAIL = ffail,time=tt,/deriv,$
+#data=d
+axis_pos = axis.read()
+result = lscan(axis, beam_ref, -dx, dx, 20, latency = 0.3, relative = True, context = None, before_read = None, after_read = None)
+(ydata, xdata) = (result.getReadable(0), result.getPositions(0))
+(norm, mean, sigma) = fit(ydata, xdata) 
+if (mean is not None) and (mean <= (axis_po + dx))  and  (mean >= (axis_po - dx)):
+    caput(axis, mean)
+else:
+    print 'fit failed - centering on maximum'
+    max_x= xdata[ydata.index(max(ydata))]
+    caput(axis, max_x)    
+
+if seti0 and not ffail:
+    time.sleep(2)
+    run("setI0")    
+#add after_rock BPM1:SUM as reference for monitoring monochromator thermal drift
+if (noref==False) and not ffail:
+    time.sleep(3) # wait a bit until mono theta finish moving        
+    br = beam_ref.read()
+    caput('X06DA-OP-MO:MONOBEAM', br)
+    print 'rock_success new mono beam at ' + str(br) + ' at' + str(e) + ' keV'