I_set = Channel("SINEG01-MSOL130:I-SET")
I_read = Channel("SINEG01-MSOL130:I-READ")
bpm_x = ai1        # Channel("SINEG01-DBPM340:X1")
bpm_y = ai2        # Channel("SINEG01-DBPM340:Y1")

def ccr(mag):
    n = 1
    while n > 0:
        sleep(0.2)
        n = caget(mag + ":I-COMP")

I1 = 2.0
I2 = 15.0
dI = 1.0

mag = [ "SINEG01-MCRX120",
        "SINEG01-MCRY120"]
for m in mag:
    ch = (m + ":I-SET")
    caput(ch, 0.0)
for m in mag:
    ccr(m)
#ccr("SINEG01-MSOL130")

r = lscan(I_set, [I_read, bpm_x, bpm_y], I1, I2, dI)

plot(r.getReadable(2), xdata=r.getReadable(1), title = "Centroid excursion")
set_return(r.print())

msg = str(r)
msg = msg + "\nFile: " + get_context().path
msg = msg + "\n\n" + r.print()
elog("Gun solenoid current scan", msg , get_plot_snapshots())