pshell_config/x03da
0
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Script execution

Browse Source
This commit is contained in:
gac-x03da
2018-09-29 01:05:39 +02:00
parent acb6a72b6c
commit a830fdd450
1 changed files with 198 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

198
script/users/TUM_PhD.py Normal file
View File

@@ -0,0 +1,198 @@
"""
Arguments:
"""
LATENCY = 0.0
MOTORS = [Eph]
#262eV : KE = 70eV for B 1s
#100eV : KE = for Ge 3d
POSITIONS_B = [262]
for nn in range(20):
POSITIONS_B.append(POSITIONS_B[-1]+4)
# BORON 1s
REGION_B = {'name': 'B 1s peak+tail', 'elo': 62., 'ehi': 75., 'estep': 0.10, 'epass': 20, 'tstep': 1, 'iter': 1, 'cis': True}
REGIONS_B = [REGION_B]
# XSP check for beam damage
POSITIONB = [550]
REGION3 = {'name': 'N 1s peak', 'elo': 134., 'ehi': 154, 'estep': 0.05, 'epass': 50, 'tstep': 0.5, 'iter': 1, 'cis': True}
REGIONSB = [REGION3]
CLOSE_SHUTTER_AT_END = True
Scienta.setAcquisitionMode(ch.psi.pshell.epics.Scienta.AcquisitionMode.Swept)
class SpectrumReader(ReadonlyRegisterBase, ReadonlyRegisterArray):
def initialize(self):
#super(SpectrumReader, self).initialize()
self.scan_index = -1
def create_datasets(self):
path = get_exec_pars().scanPath + self.region_name + "/"
self.channel_begin_dataset_name = path + "ScientaChannelBegin"
self.channel_end_dataset_name = path + "ScientaChannelEnd"
self.pass_energy_dataset_name = path + "ScientaPassEnergy"
self.step_energy_dataset_name = path + "ScientaStepEnergy"
self.step_time_dataset_name = path + "ScientaStepTime"
self.iterations_dataset_name = path + "ScientaIterations"
create_dataset(self.channel_begin_dataset_name, 'd')
create_dataset(self.channel_end_dataset_name, 'd')
create_dataset(self.pass_energy_dataset_name, 'd')
create_dataset(self.step_energy_dataset_name, 'd')
create_dataset(self.step_time_dataset_name, 'd')
create_dataset(self.iterations_dataset_name, 'd')
def setup(self):
if self.scan_index != get_exec_pars().index:
self.scan_index = get_exec_pars().index
self.create_datasets()
if self.region_index == 0:
print "scan {0}".format(self.scan_index)
ephot = Eph.read()
try:
if self.region['cis']:
edelta = ephot - self.ephot_start
else:
edelta = 0.0
except AttributeError:
self.ephot_start = ephot
edelta = 0.0
elo = self.region['elo'] + edelta
ehi = self.region['ehi'] + edelta
Scienta.lowEnergy.write(elo)
Scienta.highEnergy.write(ehi)
Scienta.stepSize.write(self.region['estep'])
Scienta.setPassEnergy(self.region['epass'])
Scienta.stepTime.write(self.region['tstep'])
Scienta.setIterations(self.region['iter'])
Scienta.update()
append_dataset(self.channel_begin_dataset_name, elo)
append_dataset(self.channel_end_dataset_name, ehi)
append_dataset(self.pass_energy_dataset_name, self.region['epass'])
append_dataset(self.step_energy_dataset_name, self.region['estep'])
append_dataset(self.step_time_dataset_name, self.region['tstep'])
append_dataset(self.iterations_dataset_name, self.region['iter'])
def read(self):
global current_region_index
current_region_index = self.region_index
self.setup()
print "Photon energy {ephot} eV. Acquiring spectrum {elo}-{ehi} eV.".format(ephot=Eph.read(), elo=Scienta.lowEnergy.value, ehi=Scienta.highEnergy.value)
trig_scienta()
time.sleep(0.5)
sp = Scienta.getSpectrum().read()
return sp
def getSize(self):
nx = int((self.region['ehi'] - self.region['elo']) / self.region['estep']) + 1
return nx
class ImageReader(ReadonlyRegisterBase, ReadonlyRegisterMatrix):
def read(self):
return Scienta.getDataMatrix().read()
def getWidth(self):
nx = int((self.region['ehi'] - self.region['elo']) / self.region['estep']) + 1
return nx
def getHeight(self):
ny = Scienta.slices.read()
return ny
def setup_live_plots(regions):
global live_plots
global current_region_index
names = [region['name'] for region in regions]
live_plots = plot(None, names, title="Live Spectra")
current_region_index = 0
def update_live_plots():
global live_plots
global current_region_index
try:
while get_context().state.running:
y = Scienta.spectrum.take(100)
x = Scienta.spectrumX
try:
series = live_plots[current_region_index].getSeries(0)
series.setData(x, y)
except IndexError:
pass
time.sleep(1.0)
finally:
print "Stopping live spectra"
def do_scan(motors, positions, regions, latency):
global SENSORS
SENSORS = []
for (index, region) in enumerate(regions):
reader = SpectrumReader()
reader.region_index = index
reader.region_name = "region{0}".format(index + 1)
reader.region = region
reader.initialize()
set_device_alias(reader, reader.region_name + "/ScientaSpectrum")
SENSORS.append(reader)
image = ImageReader()
image.region_index = index
image.region = region
image.initialize()
set_device_alias(image, reader.region_name + "/ScientaImage")
SENSORS.append(image)
SENSORS.append(SampleCurrent)
SENSORS.append(RefCurrent)
adjust_sensors()
set_adc_averaging()
vscan(motors, SENSORS, positions, True, latency,False, before_read=wait_beam, after_read = after_readout)
for (index, region) in enumerate(regions):
set_attribute(get_exec_pars().scanPath + "region{0}/ScientaSpectrum".format(index + 1), "RegionName", region['name'])
set_attribute(get_exec_pars().scanPath + "region{0}/ScientaImage".format(index + 1), "RegionName", region['name'])
set_attribute(get_exec_pars().scanPath, "Regions", [region['name'] for region in regions])
try:
setup_live_plots(REGIONSB)
task = fork(update_live_plots)
# theta, x, y, z
sample_pos = [ [-9.1,-0.5, 0.20, 115.0]]
jj = 0
for nn in sample_pos:
ManipulatorX.move(nn[1])
ManipulatorY.move(nn[2])
ManipulatorZ.move(nn[3])
ManipulatorTheta.move(nn[0])
#do_scan(MOTORS, POSITIONB, REGIONSB, LATENCY) #XPS before PhD
do_scan(MOTORS, POSITIONS_B, REGIONS_B, LATENCY) #PhD
#do_scan(MOTORS, POSITIONB, REGIONSB, LATENCY) #XPS after PhD
# if jj == 0:
# do_scan(MOTORS, POSITIONB, REGIONSB, LATENCY) #XPS before PhD
# do_scan(MOTORS, POSITIONS_N, REGIONS_N, LATENCY) #PhD
# do_scan(MOTORS, POSITIONB, REGIONSB, LATENCY) #XPS after PhD
# jj = 1
# elif jj == 1:
# do_scan(MOTORS, POSITIONB, REGIONSB, LATENCY) #XPS before PhD
# do_scan(MOTORS, POSITIONS_C, REGIONS_C, LATENCY) #PhD
# do_scan(MOTORS, POSITIONB, REGIONSB, LATENCY) #XPS after PhD
# jj = 2
# else:
# do_scan(MOTORS, POSITIONB, REGIONSB, LATENCY) #XPS before PhD
# do_scan(MOTORS, POSITIONS_O, REGIONS_O, LATENCY) #PhD
# do_scan(MOTORS, POSITIONB, REGIONSB, LATENCY) #XPS after PhD
# jj = 0
finally:
if CLOSE_SHUTTER_AT_END:
after_scan()