pshell_config/sf-op
0
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Closedown

Browse Source
This commit is contained in:
sfop
2017-06-12 11:46:48 +02:00
parent 4c95a51af6
commit 7ae0abd3a0
7 changed files with 59 additions and 56 deletions
Download Patch File Download Diff File Expand all files Collapse all files
script/RFscan/GunScan.py
+1 -1
View File
@@ -45,7 +45,7 @@ try:
r = lscan(phase, [xb, dxb], start, stop, step , latency=lat)
rf_phase = r.getPositions(0)
E = [energy0 * (1 + val.mean / 1e6 / disp) for val in r.getReadable(0)]
dE = [energy0 * ( val.mean / 1e6 / abs(disp)) for val in r.getReadable(1)]
dE = [abs(energy0 * (val.mean / 1e6 / disp)) for val in r.getReadable(1)]
finally:
phase.write(phase0)
#phase.close()
script/RFscan/SchottkyScan.py
+40 -42
View File
@@ -1,63 +1,61 @@
import ch.psi.pshell.epics.ControlledVariable as ControlledVariable
if get_exec_pars().source == CommandSource.ui:
start = -10.0
stop = 180.0
start = -30.0
stop = 30.0
step = 5.0
nb = 3
lat = 0.300
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]
bphase = ControlledVariable("Beam phase", "SINEG01-RSYS:SET-BEAM-PHASE", "SINEG01-RSYS:GET-BEAM-PHASE")
bphase.config.minValue =-180.0
bphase.config.maxValue = 360.0
bphase.config.precision = 3
bphase.config.resolution = 1.0
bphase.config.save()
bphase.initialize()
st = Stream("Schottky", dispatcher)
q = st.addScalar("Charge", "SINEG01-DICT215:B1_CHARGE", 1, 0)
phase = ControlledVariable("Phase", "SINEG01-RSYS:SET-VSUM-PHASE", "SINEG01-RSYS:GET-VSUM-PHASE")
phase.config.minValue =-180.0
phase.config.maxValue = 360.0
phase.config.precision = 3
phase.config.resolution = 1.0
phase.config.save()
phase.initialize()
st = Stream("ICTstream", dispatcher)
q = st.addScalar("Charge", "SINEG01-DICT215:B1_CHARGE", 1, 0)
st.initialize()
st.start()
rphase = Channel("SINEG01-RSYS:GET-VSUM-PHASE", type = 'd', alias = 'RF phase')
bphase0 = bphase.read()
st.waitValueChange(10000)
phase0 = phase.read()
try:
q_averager = create_averager(q, nb, 0.100)
rphase_averager = create_averager(rphase, nb, 0.100)
r = lscan(bphase, (q_averager, rphase_averager), start, stop, step, latency=lat)
beamphase = r.getPositions(0)
charge = [val.mean for val in r.getReadable(0)]
chargerms = [val.stdev for val in r.getReadable(0)]
rfphase = [val.mean for val in r.getReadable(1)]
rfphaserms = [val.stdev for val in r.getReadable(1)]
qb = create_averager(q, nb, 0.100)
r = lscan(phase, (qb), start, stop, step, latency=lat)
rf_phase = r.getPositions(0)
charge = [val.mean for val in r.getReadable(0)]
chargerms = [val.stdev for val in r.getReadable(0)]
finally:
bphase.write(bphase0)
bphase.close()
phase.write(phase0)
phase.close()
st.close()
q.close()
rphase.close()
#Setting the return value
print "MAX CHARGE = " , max(charge)
index_max = charge.index(max(charge))
print "INDEX = " , index_max
rphase_ref = rfphase[index_max] - 80
bphase_ref = bphase[index_max] - 80
set_return([rphase_ref, bphase_ref])
x = rf_phase
y = charge
index_max = y.index(max(y))
phase_ref = x[index_max] - 80
"""
# save the entry in the logbook
#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" % nb + "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())
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(phase_ref)