84 lines
2.6 KiB
Python
84 lines
2.6 KiB
Python
import ch.psi.pshell.epics.Positioner as Positioner
|
|
import ch.psi.pshell.epics.ChannelDouble as ChannelDouble
|
|
|
|
dry_run = True
|
|
do_elog = False
|
|
if get_exec_pars().source == CommandSource.ui:
|
|
start = -30.0
|
|
stop = 30.0
|
|
step = 5.0
|
|
nb = 1
|
|
lat = 0.100
|
|
plt = None
|
|
else:
|
|
start = args[0]
|
|
stop = args[1]
|
|
step = args[2]
|
|
nb = int(args[3])
|
|
lat = args[4]
|
|
plt = args[5]
|
|
|
|
if dry_run:
|
|
bph = Positioner("Beam phase", "SINEG01-RSYS:SET-BEAM-PHASE-SIM", "SINEG01-RSYS:SET-BEAM-PHASE-SIM")
|
|
q = ChannelDouble('ICT-Q', "SINEG01-DICT215:B1_CHARGE-SIM")
|
|
q.initialize()
|
|
q.monitored=True
|
|
else:
|
|
bph = Positioner("Beam phase", "SINEG01-RSYS:SET-BEAM-PHASE", "SINEG01-RSYS:SET-BEAM-PHASE")
|
|
st = Stream("ICTstream", dispatcher)
|
|
q = st.addScalar("Charge", "SINEG01-DICT215:B1_CHARGE", 1, 0)
|
|
st.initialize()
|
|
st.start()
|
|
st.waitValueChange(10000)
|
|
|
|
bph.config.minValue =-180.0
|
|
bph.config.maxValue = 360.0
|
|
bph.config.precision = 3
|
|
bph.config.rotation = True
|
|
bph.config.resolution = 1.0
|
|
bph.config.save()
|
|
bph.initialize()
|
|
rph = ChannelDouble('RF phase', "SINEG01-RSYS:SET-VSUM-PHASE-SIM")
|
|
rph.initialize()
|
|
rph.monitored=True
|
|
|
|
rph0 = rph.read()
|
|
|
|
try:
|
|
rph_averager = create_averager(rph, nb, 0.1) # Set polling time to -1 for BS data to get all messages
|
|
q_averager = create_averager(q, nb, 0.1)
|
|
q_averager.monitored=True
|
|
r = lscan(bph, (rph_averager, q_averager), start, stop, step, latency=lat)
|
|
beamphase = r.getPositions(0)
|
|
rfphase = [val.mean for val in r.getReadable(0)]
|
|
rfphaserms = [val.stdev for val in r.getReadable(0)]
|
|
charge = [val.mean for val in r.getReadable(1)]
|
|
chargerms = [val.stdev for val in r.getReadable(1)]
|
|
finally:
|
|
rph.write(rph0)
|
|
bph.close()
|
|
rph.close()
|
|
q.close()
|
|
|
|
#Setting the return value
|
|
index_max = charge.index(max(charge))
|
|
bph_ref = beamphase[index_max] - 80
|
|
rph_ref = rfphase[index_max] - 80
|
|
|
|
print "Beam phase reference = ", bph_ref
|
|
print "RF phase reference = ", rph_ref
|
|
|
|
if do_elog:
|
|
if get_option("Generated data file:\n" + get_exec_pars().path +"\n\n" + "Save to ELOG?", "YesNo") == "Yes":
|
|
log_msg = "Start: %0.1f" % start + "deg\n"
|
|
log_msg = log_msg + "Stop: %0.1f" % stop + "deg\n"
|
|
log_msg = log_msg + "Step: %0.1f" % step + "deg\n"
|
|
log_msg = log_msg + "Nb samples: %0.0f" % nb + "\n"
|
|
log_msg = log_msg + "Latency: %0.3f" % lat + "s\n"
|
|
log_msg = log_msg + "Data file: " + get_exec_pars().path
|
|
elog("Schottky Scan", log_msg, get_plot_snapshots())
|
|
|
|
set_return([bph_ref, rph_ref])
|
|
|
|
|