import  ch.psi.pshell.epics.ChannelDouble as ChannelDouble
A1 = ChannelDouble("x",            "SARUN13-DBPM070:X-REF-FB")
A2 = ChannelDouble("y",            "SARUN13-DBPM070:Y-REF-FB")
camera_name  = "SARBD01-DSCR110"
cam_server.start(camera_name+"_sp", True)    #True for shared pipeline 
wait_cam_server_message()
S1   = cam_server.stream.getChild("intensity")
A1.initialize()
A2.initialize()
S1.initialize()
A1_init = A1.read()
A2_init = A2.read()
A1i    = A1_init-0.1
A1f    = A1_init+0.1
A2i    = A2_init-0.1
A2f    = A2_init+0.1
nstep  = 11
lat    = 1.0
nav    = 10
plt    = plot(None, title="Output")[0]
plt.clear()
plt.setStyle(plt.Style.ErrorY)
plt.addSeries(LinePlotErrorSeries("Sensor1", Color.red))
def after_sample(record, scan):
    plt.getSeries(0).appendData(record.positions[0], record.readables[0].mean, record.readables[0].stdev)
try:
    S1_averager = create_averager(S1, nav, 0.1)
    r = ascan((A1,A2), (S1_averager), (A1i,A2i), (A1f,A2f), (nstep,nstep), latency=lat, after_read = after_sample)
    Act1     = r.getPositions(0)
    S1mean   = [val.mean for val in r.getReadable(0)]
    S1rms    = [val.stdev for val in r.getReadable(0)]
finally:
    A1.write(A1_init)
    A1.close()
    S1.close()
    S2.close()