axis = th1
ffail = False
tt = 0.2
seti0 = False
dx = None
noref = False


if noref == False:
    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, mono_beam, -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)):
    axis.move(mean)
else:
    print 'fit failed - centering on maximum'
    max_x= xdata[ydata.index(max(ydata))]
    axis.move(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 = mono_beam.read()
    caput('X06DA-OP-MO:MONOBEAM', br)
    print 'rock_success new mono beam at ' + str(br) + ' at' + str(e) + ' keV'