"""
parallel photon energy and manipulator scan
"""
import math
global VECTOR, SENSORS

MODE = "swept"
LATENCY = 0.0
SENSORS = [Counts, Scienta.dataMatrix, SampleCurrent, RefCurrent, MachineCurrent]

Scienta.setAcquisitionMode(ch.psi.pshell.epics.Scienta.AcquisitionMode.Swept)
#writables = (Eph, Scienta.lowEnergy, Scienta.highEnergy, ManipulatorZ)
writables = (Eph, Scienta.lowEnergy, Scienta.highEnergy, ManipulatorZ)

adjust_sensors()
set_adc_averaging()
set_preference(Preference.PLOT_TYPES, {'Scienta spectrum':1})

n = 66
Estep=1.0
Zstart=115.5
Zstep=-0.05
Estart=135.0
Ewidth=19.5

try:
    eph = [Estart + i * Estep for i in range(n)]
    elo = [Estart-Ewidth + i * Estep for i in range(n)]
    ehi = [Estart-0.5 + i * Estep for i in range(n)]
    z =   [Zstart   + i * Zstep for i in range(n)]
    VECTOR = [[eph[i], elo[i], ehi[i], z[i]] for i in range(len(eph))]
    vscan(writables, SENSORS, VECTOR, True, LATENCY, False, before_read=before_readout, after_read = after_readout)
finally:
    after_scan()