176 lines
6.7 KiB
Python
176 lines
6.7 KiB
Python
#Parameters
|
|
'''
|
|
FIELD = "Hx"
|
|
RANGES = [(-0.2, -0.1, 1.5),(-0.1, 0.0, 1.0),(0.0, 0.1, 1.0),(0.1, 0.2, 1.5)]
|
|
ENERGIES = (707.90, 703.90)
|
|
ENERGY_CHANGE_SLEEP = 0.5
|
|
MODE = 'CIRC +'
|
|
OFFSET = -1.0
|
|
'''
|
|
"""
|
|
FIELD = field_x
|
|
RANGES =[(B1, B2, rr1),(B2,-B2,rr2),(-B2,-B1,rr1)]
|
|
ENERGIES = [E1,E2]
|
|
ENERGY_CHANGE_SLEEP = 0.5
|
|
MODE =None
|
|
OFFSET=None
|
|
"""
|
|
|
|
FIELD_PRECISION = 0.01
|
|
if FIELD == field_z: FIELD = "Hz"
|
|
field = field_z if FIELD == "Hz" else field_x
|
|
field_done = field_z_done if FIELD == "Hz" else field_x_done
|
|
ramp_rate = ramp_rate_z if FIELD == "Hz" else ramp_rate_x
|
|
|
|
#Pre-actions
|
|
|
|
if NO_BEAM_CHECK == False:
|
|
wait_beam()
|
|
if MODE is not None:
|
|
pol_mode.write(MODE)
|
|
if OFFSET is not None:
|
|
pol_offset.write(OFFSET) #caput('X07MA-ID:ENERGY-OFFS', OFFSET)
|
|
#caputq('X07MA-PC:CSCALER.INPB', '1')
|
|
|
|
if len(ENERGIES) ==2:
|
|
dif_series = plot([],"Dif", title="Dif")[0].getSeries(0)
|
|
else:
|
|
dif_series = None
|
|
|
|
set_preference(Preference.ENABLED_PLOTS, ['field', 'tey_norm', 'trans_norm'])
|
|
set_preference(Preference.PLOT_TYPES, {'tey_norm':1, 'trans_norm':1})
|
|
|
|
field_var = 'fieldAnalogZ' if (FIELD == "Hz") else 'fieldAnalogX'
|
|
scan = ManualScan(['field', 'Energy'], ['TEY', 'I0', 'trans', 'polarization', 'polAngle', 'temperature', 'RingCurrent', field_var, 'tey_norm','trans_norm'], [0.0, ENERGIES[0]], [0.0, ENERGIES[-1]], [0, len(ENERGIES)-1])
|
|
scan.start()
|
|
|
|
write_logs()
|
|
|
|
for (START_FIELD, END_FIELD, RAMP_RATE) in RANGES:
|
|
#ramp_rate.write(RAMP_RATE)
|
|
if abs(field.readback.read() - START_FIELD) > FIELD_PRECISION:
|
|
ramp_rate.write(2.0)
|
|
field.write(START_FIELD)
|
|
#if iPS == True:
|
|
#caput("X07MA-PC-MAG:X:DMD",float(START_FIELD))
|
|
|
|
if field_done is None:
|
|
#ramp_rate.write(2.0)
|
|
ramp_rate.waitInPosition(2.0,10000)
|
|
time.sleep(1.0)
|
|
caput("X07MA-ES1-IPS:ACTION", 1)
|
|
time.sleep(3.0)
|
|
wait_device_in_range(field, START_FIELD, FIELD_PRECISION)
|
|
print "Waiting field done..."
|
|
while caget("X07MA-ES1-IPS:OUTPUT_RBV", "i") != 0 :
|
|
time.sleep(0.1)
|
|
print "Done waiting"
|
|
else:
|
|
# OLD AND PRONE TO HANG
|
|
#time.sleep(3.0)
|
|
#wait_device(field_done, 1) #ramp_done.wait_for_value(1.0)
|
|
|
|
# NEW, copied from hystscan_esr, april 2021
|
|
time.sleep(3.0)
|
|
# wait_device(field_done, 1) #ramp_done.wait_for_value(1.0)
|
|
|
|
### this waiting procedure was introduced because sometimes the field does not move and the whole script gets stuck then.
|
|
difference = field.readback.read() - START_FIELD
|
|
difference0 = difference # difference0 is the field offset before entering the stabilization loop
|
|
last_difference = difference # last_difference is updated during the stabilization loop
|
|
while abs(difference) > FIELD_PRECISION:
|
|
time.sleep(5.0)
|
|
difference = field.readback.read() - START_FIELD
|
|
# resend a start ramp in case nothing is happening, ie when it is stuck in this waitfield
|
|
if abs(difference-difference0) < FIELD_PRECISION:
|
|
caput('X07MA-PC-MAG:STARTRAMP.PROC',1)
|
|
print('Field setting stuck: re-starting ramp')
|
|
# also resend a start ramp in case nothing is happening wrt to the last 5 seconds
|
|
elif abs(difference-last_difference) < FIELD_PRECISION:
|
|
caput('X07MA-PC-MAG:STARTRAMP.PROC',1)
|
|
print('Field setting stuck: re-starting ramp')
|
|
last_difference = difference
|
|
|
|
|
|
|
|
|
|
|
|
|
|
#open_valve(2.0)
|
|
|
|
if RAMP_RATE is not None:
|
|
ramp_rate.write(RAMP_RATE)
|
|
ramp_rate.waitInPosition(RAMP_RATE,10000)
|
|
|
|
print "Set end field: ", END_FIELD
|
|
field.write(END_FIELD) #caputq('X07MA-PC-PS2:M:GO.A', END_FIELD)
|
|
|
|
index = 0
|
|
while(True):
|
|
for setpoint2 in ENERGIES:
|
|
#Energy.put(setpoint2, timeout=None) # TODO: Set appropriate timeout
|
|
energy.write(setpoint2)
|
|
sleep( ENERGY_CHANGE_SLEEP ) # Settling time
|
|
|
|
#TODO: change to Controlled variable?
|
|
readback2 = energy_readback.read() #EnergyReadback.get()
|
|
#if abs(readback2 - setpoint2) > 0.5 : # TODO: Check accuracy
|
|
# raise Exception('Energy could not be set to the value ' + str(setpoint2))
|
|
|
|
#while abs(readback2 - setpoint2) > 0.1 : # TODO: Check accuracy
|
|
#print "Energy could not be set to the value ", setpoint2, "Readback:", readback2
|
|
#time.sleep(2.0)
|
|
#energy.write(setpoint2)
|
|
#sleep( ENERGY_CHANGE_SLEEP ) # Settling time
|
|
#readback2 = energy_readback.read()
|
|
|
|
|
|
#readback2 = energy_readback.read() #EnergyReadback.get()
|
|
if abs(readback2 - setpoint2) > 0.1 : # TODO: Check accuracy
|
|
#raise Exception('Energy could not be set to the value ' + str(setpoint2))
|
|
#print 'Energy could not be set to the value, try again. ' + str(setpoint2)
|
|
energy.write(setpoint2)
|
|
sleep( ENERGY_CHANGE_SLEEP ) # Settling time
|
|
|
|
#Detector field readback
|
|
field_readback = field.readback.read() #fieldReadback.get()
|
|
|
|
detector1 = signal_tey.read() #TEY.get()
|
|
detector2 = signal_i0.read() #I0.get()
|
|
detector3 = signal_trans.read() #trans.get()
|
|
detector4 = float(pol_mode.getPositions().index(pol_mode.readback.read())) #polarization.get()
|
|
detector5 = pol_angle.read() #polAngle.get()
|
|
detector6 = temperature.readback.read() #temperature.get()
|
|
detector7 = current.read()
|
|
detector8 = signal_field_analog_x.read() # fieldAnalogX.get()
|
|
|
|
tey_norm = detector1/detector2
|
|
trans_norm = detector3/detector2
|
|
|
|
if dif_series is not None:
|
|
if setpoint2 == ENERGIES[0]:
|
|
first = tey_norm
|
|
else:
|
|
dif = tey_norm-first
|
|
dif_series.appendData(field_readback,dif)
|
|
scan.append ([field_readback, setpoint2], [field_readback, readback2], [detector1, detector2, detector3, detector4, detector5, detector6, detector7, detector8, tey_norm, trans_norm])
|
|
|
|
if field_done.read() == 1: #If not ramping #ramp_done.get() == 1:
|
|
print "Reached end field. Readback=", field_readback
|
|
break
|
|
index = index+1
|
|
|
|
scan.end()
|
|
|
|
log_scan_filename()
|
|
|
|
if RAMP_RATE is not None:
|
|
if FIELD == "Hz":
|
|
print "Reset ramp rate-z"
|
|
ramp_rate.write(0.35)
|
|
ramp_rate.waitInPosition(0.35,10000)
|
|
else:
|
|
print "Reset ramp rate-x"
|
|
ramp_rate.write(2.0)
|
|
ramp_rate.waitInPosition(2.0,10000)
|