Files
x09lb/script/templates/Regine.py
gac-x09lb efdbfe9940 Closedown
2022-12-03 18:56:07 +01:00

93 lines
3.0 KiB
Python

"""
Arguments:
DRY_RUN (BOOL)
CENTER_X, CENTER_Y (FLOAT)
STEPS_X, STEPS_Y (INT)
STEP_SIZE_X, STEP_SIZE_Y (FLOAT)
EXPOSURES (list or None)
SETTLING_TIME (double)
ZIGZAG (BOOLEAN)
COMPRESSION (BOOLEAN)
NOISE (FLOAT)
"""
import random
position_plot = None if "POSITION_PLOT" not in globals() else POSITION_PLOT
#Debugging
if (get_exec_pars().args is None) and (get_exec_pars().script=="Regine"):
DRY_RUN = False
CENTER_X = 11458.0
CENTER_Y = 10398.0
EXPOSURES = [1] #unit:s
STEPS_X = 30
STEPS_Y= 0
STEP_SIZE_X = 2.5
STEP_SIZE_Y = 5
NOISE=0.0
SETTLING_TIME = 1.0
ZIGZAG = True
COMPRESSION = False
position_plot=None
if not position_plot: position_plot=plot(None, title="Motor Positions")[0]
#Constants
SENSORS = [current, ccd.dataMatrix, interf_0, interf_1]
POSITIONERS = [dummy_x, dummy_y] if DRY_RUN else [sample_x, sample_y]
ENABLED_PLOTS = [ccd.dataMatrix]
CUSTOM_PLOT_TYPES = {ccd.dataMatrix:"ch.psi.pshell.plot.MatrixPlotRenderer"}
RANGE_X=[CENTER_X - STEP_SIZE_X*STEPS_X, CENTER_X + STEP_SIZE_X*STEPS_X]
RANGE_Y=[CENTER_Y - STEP_SIZE_Y*STEPS_Y, CENTER_Y + STEP_SIZE_Y*STEPS_Y]
position_plot.clear()
position_plot.addSeries(LinePlotSeries("positions"))
position_plot.getAxis(AxisId.X).label = POSITIONERS[0].name
position_plot.getAxis(AxisId.X).setRange(RANGE_X[0]-STEP_SIZE_X, RANGE_X[1]+STEP_SIZE_X)
position_plot.getAxis(AxisId.Y).label = POSITIONERS[1].name
position_plot.getAxis(AxisId.Y).setRange(RANGE_Y[0]-STEP_SIZE_Y, RANGE_Y[1]+STEP_SIZE_Y)
def after_read(record, scan):
after_readout(record, scan)
position_plot.getSeries(0).appendData(record[POSITIONERS[0]], record[POSITIONERS[1]])
if EXPOSURES:
POSITIONERS = POSITIONERS + [exposure_index()]
SENSORS = SENSORS + [exposure()]
RANGE_E=[0, len(EXPOSURES)-1]
def gen():
x_index = y_index = e_index = 0
for x_step in range(-STEPS_X, STEPS_X+1):
range_y= range(STEPS_Y, -STEPS_Y-1, -1) if (ZIGZAG and(x_index%2 ==1)) else range(-STEPS_Y, STEPS_Y+1)
for y_step in range_y:
xpos = CENTER_X + x_step*STEP_SIZE_X + (random.random()-0.5)*NOISE*2*STEP_SIZE_X
ypos = CENTER_Y + y_step*STEP_SIZE_Y + (random.random()-0.5)*NOISE*2*STEP_SIZE_Y
if EXPOSURES:
range_e= range(len(EXPOSURES)-1,-1, -1) if (ZIGZAG and(y_index%2 ==1)) else range(len(EXPOSURES))
for e in range_e:
yield([xpos,ypos,e])
else:
yield([xpos,ypos])
y_index=y_index+1
x_index=x_index+1
try:
r = vscan(POSITIONERS, SENSORS , gen(), False,\
SETTLING_TIME, relative=False, zigzag = ZIGZAG, initial_move=False, \
before_read=before_readout, after_read = after_read, \
manual_range=RANGE_X, manual_range_y=RANGE_Y, \
compression = COMPRESSION, enabled_plots=ENABLED_PLOTS, \
keep=False, check_positions=False, plot_types=CUSTOM_PLOT_TYPES)
set_return(r)
finally:
after_scan()