Skip to content
Snippets Groups Projects

Compare revisions

Changes are shown as if the source revision was being merged into the target revision. Learn more about comparing revisions.

Source

Select target project
No results found

Target

Select target project
  • 29id/iexcode
1 result
Show changes
Commits on Source (27)
Hide whitespace changes
Inline Side-by-side
Showing
with 563 additions and 552 deletions
iexcode/APRES_cheatsheet.md
View file @ acc5252a
......@@ -49,12 +49,12 @@ ___
### **Motor Moves and Scans**
-------------------------------------------------------------
mprint() #prints current position
mv_sample(position_list) #moves to position
mvsample(position_list) #moves to position
-------------------------------------------------------------
mvx(val) / mvy(val) / mvz(val) / mvth(val)
scanx(start,stop,step) / dscanx(delta_start,delta_stop,step) ...
-------------------------------------------------------------
ARPES_map_sample(y_start,y_stop,y_step,z_start,z_stop,z_step)
sample_map(y_start,y_stop,y_step,z_start,z_stop,z_step)
-------------------------------------------------------------
mvfocus(x_val)
mvy_fixed_focus(y_val)
......
iexcode/instruments/AD_utilities.py
View file @ acc5252a
......@@ -21,7 +21,13 @@ from iexcode.instruments.scanRecord import mda_filepath, mda_ioc
def AD_exposure_time(ADplugin,exp_time,**kwargs):
"""
sets the aquiretime and the exposure time
overwrites what is in kwargs
ADplugin = $(P)$(SavePlugin) (e.g. 29id_ps1:TIFF1:); only needs P with the :
exp_time: exposure time in seconds
**kwargs:
P = camera ioc (default is to split ADplugin)
R = "cam1:" by default
"""
kwargs.setdefault("P",ADplugin.split(":")[0]+":")
kwargs.setdefault("R","cam1:")
......@@ -36,7 +42,8 @@ def AD_filepath(ADplugin):
"""
returns the current directory for area detector SavePlugin
handles both Winodws and linux IOCs
ADplugin = "29idc_ps1:TIFF1:"; $(P)$(SavePlugin)
ADplugin = $(P)$(SavePlugin) (e.g. 29id_ps1:TIFF1:)
"""
SubDir=caget(ADplugin+"FilePath",as_string=True)
if SubDir[0] == 'X':
......@@ -57,6 +64,8 @@ def AD_filepath(ADplugin):
def AD_prefix(ADplugin):
"""
returns the prefix for AreaDetector plugin based on ADplugin
ADplugin = $(P)$(SavePlugin) (e.g. 29id_ps1:TIFF1:)
"""
prefix = caget(ADplugin+"FileName_RBV",as_string=True)
return prefix
......@@ -64,6 +73,7 @@ def AD_prefix(ADplugin):
def AD_EnableStats(ADplugin):
"""
Enabling the statistics in an AreaDector
ADplugin = "29idc_ps1:Stats1:"; (ADplugin=$(P)$(StatsPlugin))
"""
caput(ADplugin+"EnableCallbacks","Enable")
......@@ -73,6 +83,8 @@ def AD_EnableStats(ADplugin):
def AD_last_filename(ADplugin,verbose=True):
"""
returns the full path for the last saved image
ADplugin = $(P)$(SavePlugin) (e.g. 29id_ps1:TIFF1:)
"""
fpath = caget(ADplugin+"FullFileName_RBV",as_string=True)
if verbose:
......@@ -128,6 +140,7 @@ def AD_SaveFileSetup(ADplugin,**kwargs):
def AD_CurrentPrefix(ADplugin):
"""
returns the prefix (without _) for area detector SavePlugin
ADplugin = "29id_ps1:TIFF1:"; $(P)$(SavePlugin)
"""
Prefix=caget(ADplugin+'FileName',as_string=True)
......@@ -136,6 +149,7 @@ def AD_CurrentPrefix(ADplugin):
def AD_CurrentRun(ADplugin):
"""
returns the curent run specified in the filepath for area detector SavePlugin
ADplugin = "29id_ps1:TIFF1:"; $(P)$(SavePlugin)
"""
fpath=caget(ADplugin+"FilePath",as_string=True)
......@@ -145,6 +159,7 @@ def AD_CurrentRun(ADplugin):
def AD_CurrentUser(ADplugin):
"""
returns the curent user specified in the filepath for area detector SavePlugin
ADplugin = "29id_ps1:TIFF1:"; $(P)$(SavePlugin)
"""
folder_name = ADplugin.split(":")[1]
......@@ -165,11 +180,11 @@ def AD_DoneSingleSave(ADplugin,**kwargs):
Acquire -> Done
ImageMode -> Single
Save -> Enable
e.g. ADplugin="29id_ps2:TIFF1:"
e.g. ADplugin="29id_ps2:TIFF1:"; only needs ioc with :
**kwargs:
P=first part of ADplugin if not specified
R="cam1:"; other AD have "det1:"
P = cam ioc, first part of ADplugin if not specified
R = "cam1:"; other AD have "det1:"
"""
kwargs.setdefault("P",ADplugin.split(":")[0]+":")
kwargs.setdefault("R","cam1:")
......@@ -184,10 +199,9 @@ def AD_FreeRun(ADplugin,**kwargs):
Acquire -> Done
ImageMode -> Single
e.g. ADplugin="29id_ps2:TIFF1:"
e.g. ADplugin="29id_ps2:TIFF1:" only needs ioc with :
**kwargs:
P=first part of ADplugin if not specified
R="cam1:"; other AD have "det1:"
"""
kwargs.setdefault("P",ADplugin.split(":")[0]+":")
......@@ -200,6 +214,9 @@ def AD_FreeRun(ADplugin,**kwargs):
def AD_Done(ADplugin,**kwargs):
"""
sets acquisition to Done
ADplugin="29id_ps2:TIFF1:" only needs ioc with :
**kwargs:
P=first part of ADplugin if not specified
R="cam1:"; other AD have "det1:"
......@@ -254,21 +271,23 @@ def AD_snap(ADplugin,**kwargs):
def AD_ScanTrigger(ADplugin,**kwargs):
"""
Add Triggering of AreaDetector to scanIOC
ADplugin = "29idc_ps1:TIFF1:" (ADplugin=$(P)$(SavePlugin))
**kwargs
scanIOC = "29id"+BL_ioc() if not specified
scanIOC = default is mda_ioc()
scanDIM = 1
P=first part of ADplugin if not specified
R="cam1:"; other AD have "det1"
detTrig = 2; detectorTrigger number
"""
kwargs.setdefault("scanIOC",mda_ioc())
kwargs.setdefault("scanDIM",1)
kwargs.setdefault("P",ADplugin.split(":")[0]+":")
kwargs.setdefault("R","cam1:")
kwargs.setdefault("detTrig",2)
scanPV=mda_ioc()+"scan"+str(kwargs["scanDIM"])
scanPV = kwargs["scanIOC"]+"scan"+str(kwargs["scanDIM"])
trigger=".T"+str(kwargs["detTrig"])+"PV"
caput(scanPV+trigger,kwargs["P"]+kwargs["R"]+"Acquire",wait=True,timeout=5*60)
......@@ -297,8 +316,11 @@ def ADplugin_ScanSetup(ADplugin, **kwargs):
"""
#from AD_ScanTrigger
kwargs.setdefault("scanIOC",mda_ioc())
kwargs.setdefault("scanDIM",1)
kwargs.setdefault("P",ADplugin.split(":")[0]+":")
kwargs.setdefault("R","cam1:")
kwargs.setdefault('detTrig',2)
AD_DoneSingleSave(ADplugin,**kwargs)
......@@ -306,11 +328,22 @@ def ADplugin_ScanSetup(ADplugin, **kwargs):
AD_exposure_time(ADplugin,kwargs["ExposureTime"])
AD_SaveFileSetup(ADplugin,**kwargs)
AD_ScanTrigger(ADplugin, **kwargs)
trigger = ".T"+str(kwargs["detTrig"])+"PV"
scanPV = mda_ioc()+"scan"+str(kwargs["scanDIM"])
AD_ScanTrigger(ADplugin, **kwargs)
print("WARNING: you need to need to disable saving and clear the trigger by hand after the scan")
print("\tAD_FreeRun("+ADplugin+"); caput("+scanPV+trigger+",'')")
print("\tAD_FreeRun("+ADplugin+"); scan_triggers_clear()")
def AD_attributes_set(ADplugin, filepath):
"""
filepath is the full path to the xml file
filepath = '/xorApps/epics/synApps_6_1/ioc/29idcScienta/iocBoot/ioc29idcScienta/HDF5Attributes.xml'
filepath = '/home/beams/29IDUSER/Documents/User_Macros/Macros_29id/AD_Attributes/DiaGon_Attributes.xml'
"""
current_xml = caget(ADplugin+"NDAttributesFile",as_string=True)
print ('Changed attributes file from ',current_xml)
caput(ADplugin+"NDAttributesFile",filepath)
sleep(1)
print('To: ',caget(ADplugin+"NDAttributesStatus",as_string=True))
def AD_ROI_setup(AD,ROInum,xcenter=500,ycenter=500,xsize=50,ysize=50,binX=1,binY=1):
"""
......
iexcode/instruments/ARPES.py
View file @ acc5252a
......@@ -82,6 +82,8 @@ def ARPES_init(*userName,**kwargs):
user_name = mda_user_name()
else:
user_name = userName[0]
#if you give a userName then you set the folder to that user
kwargs['set_folders']=True
if user_name == 'staff':
user_name = 'Staff'
kwargs['mode'] = 'staff'
......@@ -580,7 +582,7 @@ def _ARPES_DefaultPosition(destination):
DefaultPosition={
'measure':(0.0,0.0,-13.5,0.0,0.0,0.0),
'LEED':(-0.0, 0.0, -141.5, 89.5, 0.0, 0.0),
'transfer':(-0,0.0,-141.5,180.0,0.0,0.0),
'transfer':(-2.0,0.0,-141.5,180.0,0.0,0.0),
}
if destination in DefaultPosition:
pos=DefaultPosition[destination]
......
iexcode/instruments/IEX_BL_config.py
View file @ acc5252a
......@@ -88,8 +88,8 @@ def get_all(verbose=True):
"""
vals={}
vals.update(iex.BL.endstation_get(verbose))
return vals
if 'verbose'==False:
return vals
##############################################################################################################
......
iexcode/instruments/IEX_VPU.py
View file @ acc5252a
......@@ -413,7 +413,10 @@ def ID_energy_set(keV,verbose=True):
if QP_ratio[1] == None:
#ID_off() #JM look here
#ID_on()
print_warning_message('check QP')
if ID_get() > 3.7:
message = 'ID off'
else:
print_warning_message('check QP')
eV = keV*1000.0
keV = ID_calc_eV(mono_grating,ID_mode,eV,QP_ratio[0])/1000.0
ID_set(keV,verbose=verbose)
......
iexcode/instruments/Kappa.py
View file @ acc5252a
......@@ -38,6 +38,8 @@ from iexcode.instruments.spec_stuff import folders_spec
default_ioc = '29idKappa:'
kth_offset_pv = '29idKappa:userCalcOut1.G'
physical_motors = ['x','y','z','tth','kth','kap','kphi']
pseudo_motors = ['th','chi','phi']
#############################################################################
def kappa_init(*userName,**kwargs):
"""
......@@ -66,8 +68,6 @@ def kappa_init(*userName,**kwargs):
#motors
physical_motors = ['x','y','z','tth','kth','kap','kphi']
pseudo_motors = ['th','chi','phi']
global kappa_Motors
kappa_Motors = Motors('kappa',_kappa_motor_dictionary(),physical_motors,pseudo_motors)
......@@ -86,6 +86,8 @@ def kappa_init(*userName,**kwargs):
user_name = mda_user_name()
else:
user_name = userName[0]
#if you give a userName then you set the folder to that user
kwargs['set_folders']=True
if user_name == 'staff':
user_name = 'Staff'
kwargs['mode'] = 'staff'
......@@ -99,23 +101,22 @@ def kappa_init(*userName,**kwargs):
folders_kappa(user_name,**kwargs)
#MPA
global MPA
iex.mpa = MPA()
MPA=iex.mpa
global mpa
mpa = MPA()
#global detectors
global tth_pv,tthdet
tthdet = Kappa_Detector(kappa_Motors)
tthdet = Kappa_Detector(kappa_Motors,tth_detector_dictionary())
tth_pv = _kappa_motor_dictionary()['tth'][3]
iex.tth = tthdet
#global detectors
global tey,d3,d4,mesh
tey = SRS("29idMZ0:scaler1.S2", '29idd:A1')
d3 = SRS("29idMZ0:scaler1.S3", '29idd:A2')
d4 = SRS("29idMZ0:scaler1.S4", '29idd:A3')
mesh = SRS("29idMZ0:scaler1.S14", '29idd:A4')
tthdet
#logging (logging is dependent on global detectors)
if kwargs['set_folders']:
......@@ -124,9 +125,6 @@ def kappa_init(*userName,**kwargs):
log_init('kappa',user_name,_kappa_log_dictionary,set_file=False)
#default scan settings
kappa_default_scan_settings()
......@@ -143,15 +141,15 @@ def kappa_init(*userName,**kwargs):
#stores detecotrs
iex.kappa_detectors_dictionary={
kappa_global_detectors = {
'tey':tey,
'd3':d3,
'd4':d4,
'mesh':mesh,
'tthdet':tthdet
'tthdet':tthdet,
'mpa':mpa,
}
iex.kappa_detectors_dictionary = kappa_global_detectors
#return any detectors or other parameters that you want access to from jupyter
return tey,d3,d4,mesh,tth_pv,tthdet
......@@ -198,14 +196,20 @@ def kappa_motor_detectors():
57:"29idKappa:Euler_PhiRBV",
}
return motors
def tth_detector_dictionary():
detNums={
'tthdet':15,
'd3':33,
'd4':34,
'mpa':30
}
"""
Used by Kappa_det to set detectors
"""
detNums = {}
detNums['tth_detNum']=15
detNums['d3']="29idMZ0:scaler1.S3" #33
detNums['d4']="29idMZ0:scaler1.S4" #34
detNums['mpa']="29iddMPA:det1:TotalRate_RBV" #30
detNums['apd']=""
detNums['yag']=""
return detNums
def _kappa_detector_dictionary(**kwargs):
"""
......@@ -266,6 +270,7 @@ def _kappa_detector_dictionary(**kwargs):
det_dict.update(scalers)
det_dict.update(motors)
det_dict.update(tthdet_detector)
det_dict.update(m3r)
if kwargs['add_vortex']:
det_dict.update(vortex)
if kwargs['add_MPA']:
......@@ -409,7 +414,7 @@ def kappa_reset(**kwargs):
#resetting the MPA ROIs
if kwargs['reset_MPA_ROIs']:
MPA.mpa_ROI_reset()
mpa.mpa_ROI_reset()
#resetting mono and other beamline stuff
if BLconfig_xrays():
......@@ -760,6 +765,13 @@ def PI_motors_sync():
caput(val_pv,current_rbv)
print('PI motors VAL synced to RBV')
def PI_motor_servo(Enable_Disable='Enable'):
motor_dictionary = _kappa_motor_dictionary()
for motor in _PI_motor_list():
servo_pv = motor_dictionary[motor][3]+'.CNEN'
caput(servo_pv,Enable_Disable)
print('PI motors servos '+Enable_Disable)
def Euler_motors_sync():
"""
Syncs the Euler motores
......@@ -962,18 +974,24 @@ def scanth2th_table(tth_table,th_offset, **kwargs):
#############################################################################################################
############################## Preset Positions ##############################
##############################################################################################################
def _kappaTransfer_StrSeq():
def _kappaTransfer_StrSeq(**kwargs):
kwargs.setdefault('seq_num',4)
seq_num=kwargs['seq_num']
ioc = mda_ioc()
#User= [ DESC, x, y, z, tth, kth, kap, kphi]
User = ["kappa Transfer",0, -2650, -650, 0, 57, 0, -88]
n=4
#_kappaPreset_StrSeq(n,User)
pv,proc = userStringSeq_pvs(ioc, seq_num)
return proc
def _kappaGrazing_StrSeq(): #Need to determine positions and then add to the kappa graphic
def _kappaGrazing_StrSeq(**kwargs): #Need to determine positions and then add to the kappa graphic
kwargs.setdefault('seq_num',3)
seq_num=kwargs['seq_num']
ioc = mda_ioc()
#Dial= [ DESC, x, y, z, tth, kth, kap, kphi]
User = ["kappa Grazing",0, 0, 0, 0, 57.045, 134.76,57.045]
n=3
#_kappaPreset_StrSeq(n,User)
pv,proc = userStringSeq_pvs(ioc, seq_num)
return proc
def kappa_ResetPreset():
_kappaGrazing_StrSeq()
......@@ -1007,7 +1025,7 @@ def kappa_light(ON_OFF):
##############################################################################################################
######################## Scan Temp and Pressure ##############################
##############################################################################################################
def kappa_temperature_pressure_scan(scan_dim=1):
def kappa_temperature_pressure_scan(scan_dim=1,**kwargs):
"""
starts scan to monitor temperature and pressure in the kappa chamber
......@@ -1022,7 +1040,8 @@ def kappa_temperature_pressure_scan(scan_dim=1):
caput(pv+".D02PV","29idd:tc1:getVal_B.VAL")
caput(pv+".D03PV","29idb:VS11D.VAL")
#time scan
iex.BL.mda.time_go()
duration_min=10000
scan_time(duration_min,step_min=1,**kwargs)
......@@ -1043,7 +1062,7 @@ def kappa_safe_state(**kwargs):
if kwargs["mpa_off"]:
try:
MPA.HV_off()
mpa.HV_off()
except:
print('MPA is not running')
branch = 'd'
......
iexcode/instruments/Kappa_det.py
View file @ acc5252a
......@@ -2,20 +2,14 @@ from time import sleep
from epics import caget, caput
import iexcode.instruments.cfg as iex
#from iexcode.instruments.Kappa import tth_detector_dictionary
from iexcode.instruments.scanRecord import scan_detectors_set
##############################################################################################################
################################ Kappa detector class ##############################
##############################################################################################################
det_set_pv = '29idKappa:det:set'
det_list = ['d3', 'd4', 'mcp', 'apd', 'yag']
def tth_detectors(name):
detNums={
'd3':33,
'd4':34,
'mpa':30
}
#JM finish here
class Kappa_Detector:
......@@ -23,9 +17,11 @@ class Kappa_Detector:
class for Kappa detector
"""
def __init__(self,kappa_Motors):
def __init__(self,kappa_Motors,tth_detector_dictionary):
self.motors = kappa_Motors
self.detector_dictionary = tth_detector_dictionary
self.get()
def get(self):
"""
......@@ -59,6 +55,16 @@ class Kappa_Detector:
caput(det_set_pv,det_name)
if move:
self.motors.move('tth',tth_val,wait=True,verbose=False)
#replace det: 15 with SRS for current detector
tth_detNum = self.detector_dictionary['tth_detNum']
det_pv = self.detector_dictionary[self.name]
try:
scan_detectors_set(detector_dictionary={tth_detNum:det_pv})
except:
scan_detectors_set(detector_dictionary={tth_detNum:""})
print_warning("default tth_detNum("+str(tth_detNum)+") not set")
sleep(2)
def tth0_set(self):
"""
......
iexcode/instruments/Lakeshore_335.py
View file @ acc5252a
......@@ -10,17 +10,182 @@ from epics import caget, caput
from iexcode.instruments.IEX_BL_config import *
from iexcode.instruments.utilities import read_dict
def Lakeshore_reset(pv,d):
class Lakeshore_335:
"""
resets the lake short to the default paramters defind in the dictionary d
Lakeshore 355 data class
"""
for key in d.keys():
caput(pv+key,d[key])
def __init__(self,ioc,**kwargs):
"""
ioc = "29idARPES:"
**kwargs
Q = "TC1:"
defaults => dictionary of default parameter used by reset
"""
kwargs.setdefault('Q','TC1:')
self.pv = ioc+"LS335:"+kwargs['Q']
self.reset(**kwargs)
self.get()
pass
def reset(self,**kwargs):
"""
resets the lake short to the default paramters defind in the dictionary defaults
**kwargs
defaults = {"read.SCAN":".5 second","OUT1:Cntrl":"A","OUT2:Cntrl":"B","OUT1:Mode":"Closed Loop"}
"""
kwargs.setdefault("defaults",{"read.SCAN":".5 second","OUT1:Cntrl":"A",
"OUT2:Cntrl":"B","OUT1:Mode":"Closed Loop","P1":1,"I1":5,"D1":1,"RampR1":5,"OnRamp1":"ON"})
d=kwargs['defaults']
for key in d.keys():
caput(self.pv+key,d[key])
def heater_range_list(self):
return ["off","low","medium","high"]
def pv_dictionary(self):
d={
"INA":None,
"INB":None,
"IN1":None,
"IN2":None,
"IN1:Name_RBV":"",
"IN2:Name_RBV":"",
"OUT1:SP":None,
"OUT1:SP_RBV":None,
"HTR1":None,
"HTR1:Range":"",
"HTR1:Range_RBV":"",
"OUT1:Cntrl_RBV":"",
"P1":None,
"P1_RBV":None,
"I1":None,
"I1_RBV":None,
"D1":None,
"D1_RBV":None,
"RampR1":None,
"RampR1_RBV":None,
"OnRamp1":"",
"OnRamp1_RBV":"",
}
def get(self):
"""gets the current valuse
"""
d = self.pv_dictionary()
caput(self.pv+"read.PROC",1)
sleep(1)
caput(self.pv+"read.PROC",0)
for key in d.keys():
if d[key] != None:
caget(self.pv+key,as_string=True)
else:
caget(self.pv+key)
return d
def set(self,TA_sp,bandwidth=5,seconds=5*60,heater_range=None,out_num=1):
"""
TA_SP = temperature set point
bandwidth = the temperature range acceptable
seconds = time to wait for stablized
heater_range = "off","low","medium","high"
= None then uses current setting
out_num = control loop number (1 or 2)
"""
d = self.pv_dictionary
TA = caget(self.pv+"INA")
while True:
#sepoint with +/- current temperature
if TA - bandwidth <= TA_sp > TA + bandwidth:
#regulate off sample
self.regulate_TA(TA_sp,heater_range,out_num)
if self.TA_stable(TA_sp,bandwidth,seconds) == True:
break
else:
self.regulate_TB(TA_sp,heater_range,out_num)
caget()
def TA_stable(self,T_sp,bandwidth=5,seconds=5*60,out_num=1):
"""
returns True or False is the temperature is with in the bandwidth over a given
time interval in seconds
T_sp = temperature set point
bandwidth = +/- Temperature to be stablized
seconds = time to wait to stablize
out_num = control loop number (1 or 2)
"""
read_interval = 0.5
name_rbv = caget(self.pv+"IN"+str(out_num)+":Name_RBV", as_string=True)
for t in range(0,seconds/read_interval): #sampling every 0.5 seconds
T = caget(self.pv+"IN"+str(out_num))
if T - bandwidth <= T_sp > T + bandwidth:
return False
else:
sleep(read_interval)
return True
def regulate_TA(self,T_sp,heater_range=None,out_num=1):
"""
"""
self.Cntrl(1,'A')
self.Cntrl(2,'B')
sleep(0.5)
self.set_point(1,T_sp,heater_range)
self.set_point(2,0,heater_range='off')
def regulate_TB(self,T_sp,heater_range=None,out_num=1):
self.Cntrl(1,'B')
self.Cntrl(2,'A')
sleep(0.5)
self.set_point(1,T_sp,heater_range)
self.set_point(2,0,heater_range='off')
def Cntrl(self,out_num,sensor):
"""
out_num = control loop number (1 or 2)
sensor = 'A' or 'B'
"""
caput(self.pv+"OUT"+str(out_num)+":Cntrl",sensor)
sleep(.5)
def set_point(self,out_num,T_sp,heater_range=None):
caput(self.pv+"OUT"+str(out_num)+":SP")
self.set_heater_range(out_num,heater_range)
def set_heater_range(self,out_num,heater_range):
"""
sets the heater range
out_num = control loop number (1 or 2)
"""
if heater_range == None:
heater_range = caget(self.pv+"OUT"+str(out_num)+"HTR1:Range_RBV",as_string=True)
else:
heater_range = heater_range.upper()
if heater_range not in self.heater_range_list():
message = "heater_range = "+heater_range+" not a valid choice"
message += "\nchoose one of "+str(self.heater_range_list)
print_warning_message(message)
return
else:
caput(self.pv+"OUT"+str(out_num)+"HTR1:Range",heater_range)
sleep(.5)
return caget(self.pv+"OUT"+str(out_num)+"HTR1:Range_RBV",as_string=True)
#############################################################################################################
############################## Lakeshore 335 Diode Curves ##############################
##############################################################################################################
#############################################################################################################
#---Settings as of 12/04/2018---
#DiodeCurve_Write(22,"SI435")
#DiodeCurve_SetInput("A",22,400)
......@@ -145,7 +310,7 @@ def DiodeCurve_SetInput(Channel,Curve,Tmax):
caput(PV,cmd,wait=True,timeout=1800)
##############################################################################################################
############################## PID Settings ##############################
############################## PID Settings for open cycle ##############################
##############################################################################################################
......@@ -159,157 +324,7 @@ def DiodeCurve_SetInput(Channel,Curve,Tmax):
#Lakeshore 355 I = time in seconds for one oscilation / 1000; D=50-100 for starters
def PID_dict(which,T): #Dictionary of PID setting to try for the ARPES chamber
"""Dictionary for common PID setting for the ARPES system.
Since the PID is dependent on the cooling power, this will change with the cryogen and the flow rate:
- Flow = "LHe,LN" => LHe = needle valve = "on" and flow = 70
"""
PID={}
# PID[which,T]=[P,I,D,Range] Range=HIGH,MEDIUM,LOW,OFF
PID['RT', 378.0]=[200.0,0.0,0.0,'MEDIUM']
PID['RT', 375.0]=[200.0,0.0,0.0,'MEDIUM']
PID['RT', 338.0]=[100.0,0.0,0.0,'MEDIUM']
PID['RT', 298.0]=[20.0,0.0,0.0,'LOW']
PID['GHe', 378.0]=[210.0,0.0,0.0,'MEDIUM']
PID['GHe', 375.0]=[210.0,0.0,0.0,'MEDIUM']
PID['GHe', 338.0]=[130.0,0.0,0.0,'MEDIUM']
PID['GHe', 298.0]=[60.0,0.0,0.0,'MEDIUM']
PID['GHe', 250.0]=[10.0,0.0,0.0,'LOW']
PID['LHe',298.0]=[20.0,5.0,1.0,'HIGH'] # +/- 1.1 deg, over 10 min, heater power 53%
PID['LHe',200.0]=[10.0,5.0,1.0,'HIGH'] # +/- 1.1 deg, over 10 min, heater power 53%
PID['LHe',150.0]=[200.0,0.1,0,'HIGH'] # Stablized at 153, heater power 43%
PID['LHe',180.0]=[200.0,0.1,0,'HIGH']
PID['LHe',230.0]=[200.0,0.1,0,'HIGH']
PID['LHe',300.0]=[200.0,0.1,0,'HIGH']
PID['LHe',40.0]=[200.0,0.1,0,'HIGH'] #
PID['LN2',230]=[10.0,0.6,100,'HIGH'] #stable 237.83, needle valve=ON; flow 6 half turns from max)
PID['LN2',180]=[50.0,4,100,'MEDIUM']
return PID[which,T]
def PID_set(which,T,heater='ON'):
"""
Uses preset PID settings as defined in PID_dict()
To cool down: heater='OFF'
"""
P,I,D,Range=PID_dict(which,T)
caput("29idARPES:LS335:TC1:P1",P)
caput("29idARPES:LS335:TC1:I1",I)
caput("29idARPES:LS335:TC1:D1",D)
if heater == 'ON':
caput("29idARPES:LS335:TC1:HTR1:Range",Range)
elif heater == 'OFF':
caput("29idARPES:LS335:TC1:HTR1:Range",'OFF')
print('\nP = ',PID_dict(which,T)[0])
print('I = ',PID_dict(which,T)[1])
print('D = ',PID_dict(which,T)[2])
print('Range = ',PID_dict(which,T)[3])
def SetT_Sample(which,T,precision=5,minutes=15,offset=6):
"""
Sets PID settings for a given set point as defined in PID_dict().
Available set points:
which = 'RT' : T = 298, 338, 375
which = 'LHe': T = 150, 200
which = 'GHe': T = 250, 298, 338, 370
which = 'LN2': T = 230, 180
precision is temperature in K
"""
current_T=caget("29idARPES:LS335:TC1:IN1")
print('\nSet T to '+ str(T)+' K @ '+ dateandtime())
if T>current_T:
PID_set(which,T,'ON')
Range=PID_dict(which,T)[3]
SetT_Up(T,offset,Range,precision,minutes)
else:
PID_set(which,T,'OFF')
Range=PID_dict(which,T)[3]
SetT_Down(T,offset,Range,precision,minutes)
caput("29idARPES:LS335:TC1:OUT1:SP",T)
def SetT_Up(T,offset,Range,precision=5,minutes=15):
if Range=='LOW': Range=1
if Range=='MEDIUM': Range=2
if Range=='HIGH': Range=3
t=15
u=0
caput('29idARPES:LS335:TC1:OUT1:SP',T) # set the set point to T
while True: # while TA < T
TA=caget('29idARPES:LS335:TC1:IN1') # get current temperature at the sample
TB=caget('29idARPES:LS335:TC1:IN2') # get current temperature (B=cold finger)
if TA<T-precision: # if it hasn't reach SP
caput("29idARPES:LS335:TC1:HTR1:Range",min(Range,3)) # increase heater range to Range +1
while True: # while TB < T+offset:
TB=caget('29idARPES:LS335:TC1:IN2') # get current temperature at the cold finger
if (t%120)==0:
print('\nTA = ',TA,' TB = ',TB, ' @ ',dateandtime())
if TB<(T+offset) and t<=minutes*60: #if it hasn't reach the SP+offser
sleep(15)
t+=15
#print t, TA, TB
elif TB<(T+offset) and t>minutes*60:
heater_power=caget('29idARPES:LS335:TC1:HTR1')
heater_range=caget('29idARPES:LS335:TC1:HTR1:Range')
if heater_power > 90:
caput("29idARPES:LS335:TC1:HTR1:Range",min(heater_range+1,3))
print('Change Range to',caget('29idARPES:LS335:TC1:HTR1:Range'),' @ ',dateandtime())
elif heater_range==3 or heater_power<=90:
P=caget("29idARPES:LS335:TC1:P1")
caput("29idARPES:LS335:TC1:P1",P*1.5)
print('Change P to',caget("29idARPES:LS335:TC1:P1"),' @ ',dateandtime())
t=0
else: #else
break # break
caput("29idARPES:LS335:TC1:HTR1:Range",'OFF') # turn off the heater
elif TA>T-precision: # if it has reach the set point
break # break
print('TA = ',TA,' TB = ',TB, ' @ ',dateandtime()) # print temperatures
caput("29idARPES:LS335:TC1:HTR1:Range",Range) # set the heater range to preset value
def SetT_Down(T,offset,Range,precision=5,minutes=15):
t=0
caput('29idARPES:LS335:TC1:OUT1:SP',T) # set the set point to T
while True: # while TA < T
TA=caget('29idARPES:LS335:TC1:IN1')
TB=caget('29idARPES:LS335:TC1:IN2')
if (t%120)==0:
print('\nTA = ',TA,' TB = ',TB, ' @ ',dateandtime())
if TA>T+precision:
sleep(15)
t+=15
elif t>minutes*60:
P=caget("29idARPES:LS335:TC1:P1")
caput("29idARPES:LS335:TC1:P1",P/1.5)
t=0
else:
break
caput("29idARPES:LS335:TC1:HTR1:Range",Range)
print('TA = ',TA,' TB = ',TB, ' @ ',dateandtime())
def Get_PID(which='LHe'):
T=caget("29idARPES:LS335:TC1:IN1")
SP=caget("29idARPES:LS335:TC1:OUT1:SP")
P=caget("29idARPES:LS335:TC1:P1")
I=caget("29idARPES:LS335:TC1:I1")
D=caget("29idARPES:LS335:TC1:D1")
Range=caget("29idARPES:LS335:TC1:HTR1:Range",as_string=True)
# print SP,P,I,D,Range
print("Current T:", T,'K')
print("PID[\'"+which+"\',"+str(SP)+"]=["+str(P)+","+str(I)+","+str(D)+",\'"+Range+"\']")
#functions use from ScanFunctions_IEX import BL_ioc
def PVLakeshores(which):
"""
which = RSXS / ARPES /Hydra / DetPool
......@@ -324,98 +339,5 @@ def PVLakeshores(which):
}
return d[which]
def dateandtime():
return strftime("%a %d %b %Y %H:%M:%S",localtime())
def tempr(which=None):
"""
reads the temperature in the current branch unless otherwise specificed
"""
if which == None:
which=BL_ioc()
PVs=PVLakeshores(which)
sample=PVs[0]+PVs[1]
print(which+" Sample Temperature: {} K".format(caget(sample)))
def temps(tempset,which=None):
"""
Set sthe sample temperature set point and waits for it to stablelize for the current branch
unless specified by which = Kappa / ARPES /Hydra set by PVLakeshores
"""
if which == None:
which=BL_ioc()
P, Sample, ColdFinger, SetPoint = PVLakeshores(which)
print("Initial Sample Temperature:",caget(P+Sample)," K",dateandtime())
caput(P+SetPoint,1.0*tempset)
stop_var=0
b=0
sleep(0)
while b == 0:
delta=abs(caget(P+Sample)-1.0*tempset)
if abs(delta) < 1:
c=0
while c < 10:
sleep(1)
delta=delta+(caget(P+Sample)-1.0*tempset)
c=c+1
if abs(delta/10) < 1:
print("Stable Sample Temperature:",caget(P+Sample)," K",dateandtime())
b=1
else:
temp1=caget(P+Sample)
sleep(10)
temp2=caget(P+Sample)
if abs(temp1-temp2) < 0.1:
sleep(30)
temp2=caget(P+Sample)
if abs(temp1-temp2) < .5:
print("UNABLE TO STABLIZE TEMPERATURE! Stopped at T=:",caget(P+Sample)," K",dateandtime())
b=1
def ARPES_warming(Tmax,PowerMax=50,Tdiff=10):
"""
Ramps the temperature up for ARPES using open loop and manual power at PowerMax with
the heater on High until Tmax - Tdif is reached and then switches to
closed loop with the heater on Medium
PowerMax=50,Tdiff=10
PowerMax=55,Tdiff=25
"""
P="29idARPES:LS335:TC1:"
#Needs modified in controller changed
Tsample=P+"IN1"
Tcontrol=P+"IN2"
Tsetpoint=P+"OUT1:SP"
HeaterRange=P+"HTR1:Range"
ControlMode=P+"OUT1:Mode"
ManualPower=P+"OUT1:MOUT"
caput(Tsetpoint,Tmax)
caput(ManualPower,PowerMax)
caput(ControlMode,"Open Loop")
caput(HeaterRange,'HIGH')
print("T", caget(Tsample))
print("Started warming:", dateandtime())
while True:
TA=caget(Tsample)
TB=caget(Tsample)
if TB >= Tmax:
print ("Control is to high",TB)
break
if TA < Tmax - Tdiff:
sleep(60)
else:
break
print("T", caget(Tsample))
print("Switch to closed loop", dateandtime())
caput(HeaterRange,'MEDIUM');sleep(1)
caput(ControlMode,"Closed Loop");sleep(1)
caput(HeaterRange,'OFF');sleep(1)
caput(HeaterRange,'MEDIUM');sleep(1)
temps(Tmax,which='ARPES')
iexcode/instruments/Motors.py
View file @ acc5252a
......@@ -264,6 +264,12 @@ class Motors:
kwargs.setdefault('cts',0.1)
kwargs.setdefault('mpa',False)
if kwargs['debug']:
print('kwargs:', kwargs)
print('inner_loop_list = ',inner_loop_list)
print('outer_loop_list = ',outer_loop_list)
rbv_pv,val_pv,sgm_pv,pv = self._motor_dictionary[inner_loop_list[0]]
inner_loop_list.pop(0)
inner_loop_list.insert(0,rbv_pv)
......
iexcode/instruments/Octupole.py
View file @ acc5252a
......@@ -80,6 +80,8 @@ def Octupole_init(*userName,**kwargs):
user_name = mda_user_name()
else:
user_name = userName[0]
#if you give a userName then you set the folder to that user
kwargs['set_folders']=True
if user_name == 'staff':
user_name = 'Staff'
kwargs['mode'] = 'staff'
......
iexcode/instruments/Scienta.py
View file @ acc5252a
......@@ -436,13 +436,23 @@ class Scienta:
caput(self.PHV+"sweepStopEnergy.VAL",KEstop);time.sleep(0.25) #added 20220317 - sweep calcs need time to process or you get weird stuff
caput(self.PHV+"sweepStepEnergy.VAL",KEstep)
time.sleep(1)
#Re-Calc
#print('Yes Re-Calc')
#caput(self.PHV+"SweepReCalc.PROC",1)
#time.sleep(1)
#Load spectra parameters
caput(self.PHV+"SweepLoadFanout.PROC",1)
time.sleep(1)
#check that there are no zeros in the KE_array
KE_array = caget(self.PHV+"KE_array")
if np.min(KE_array)<1:
print("KE_array had zeros")
caput(self.PHV+"SweepReCalc.PROC",1)
time.sleep(2)
caput(self.PHV+"SweepLoadFanout.PROC",1)
time.sleep(1)
if len(KE_array) != caget(self.PHV+"scan2.NPTS"):
print("number of scan point and KE_array miss match")
caput(self.PHV+"SweepReCalc.PROC",1)
time.sleep(2)
caput(self.PHV+"SweepLoadFanout.PROC",1)
time.sleep(1)
def _spectra_BabySwept(self, KEcenter):
......@@ -520,6 +530,9 @@ class Scienta:
# writing scan parameters
else:
#put the analzer at inital value so calcuations are correct
self.put(KElist[0],PassEnergy,"Angular")
time.sleep(2)
if len(EAlist)==7: #Swept Mode
scanType="Sweep"
self._spectra_Swept(KElist[0],KElist[1],KElist[2],PassEnergy)
......
iexcode/instruments/VLS_PGM.py
View file @ acc5252a
......@@ -225,14 +225,15 @@ def mono_energy_set(hv_eV,verbose=True):
caput(pv,hv_eV,wait=True,timeout=60)
time.sleep(2.5)
mono_status = mono_status_get()
while True:
if mono_status > 1:
if mono_status_get() > 1:
mono_resest_pitch()
caput(pv,hv_eV,wait=True,timeout=60)
time.sleep(2.5)
else:
break
if verbose:
print("Mono set to",str(round(hv_eV,3)),"eV")
......@@ -247,8 +248,9 @@ def mono_scan_pvs():
def mono_status_get():
"""
returns the status of the mirror (GRT * MIR)
1 = ready
2 = moving
1 = positioned/ready
0 = somebody moving
1 > somebody not happy
Previously Mono_Status
......@@ -256,8 +258,8 @@ def mono_status_get():
pvs = _mono_pvs()
mir_P_status = int(caget(pvs['mir_P_status']))
grt_P_status = int(caget(pvs['grt_P_status']))
mir_X_status = int(caget(pvs['mir_P_status']))
grt_X_status = int(caget(pvs['grt_P_status']))
mir_X_status = int(caget(pvs['mir_X_status']))
grt_X_status = int(caget(pvs['grt_X_status']))
mirror_status = mir_P_status * mir_X_status
grating_status = grt_P_status * grt_X_status
return mirror_status*grating_status
......@@ -315,13 +317,15 @@ def _mono_optic_change(optic,name, verbose=True):
if type_get != name:
caput(type_sp,num,wait=True,timeout=18000)
caput(type_move,1,wait=True,timeout=18000)
time.sleep(1) #needed
while True:
if mono_status_get() > 1:
if mono_status_get() != 1:
time.sleep(5)
else:
break
time.sleep(1)
mono_energy_set(hv_eV)
time.sleep(1)
else:
if verbose:
......
iexcode/instruments/Vortex.py
View file @ acc5252a
......@@ -32,7 +32,7 @@ class Vortex:
#set hi and low
caput(self._vortex_pv+'mca1.R0LO',round(center-width/2))
caput(self._vortex_pv+'mca1.R0HI',round(center-width/2))
caput(self._vortex_pv+'mca1.R0HI',round(center+width/2))
#pv to store energy set point
caput(self._calc_pv,energy)
......
iexcode/instruments/cfg.py
View file @ acc5252a
......@@ -3,5 +3,4 @@ BL=None
#Kappa global stuff
mpa=None
det=None
kappa_detectors_dictionary={}
\ No newline at end of file
iexcode/instruments/electron_analyzer.py
View file @ acc5252a
......@@ -686,15 +686,16 @@ def scanEA_motor(EAlist, motor,start,stop,step,mode='absolute',**kwargs):
EAparms=scanEA(EAlist,**EAkwargs)
#Setting up the ScanRecord for motor scans
execute = kwargs['execute']
kwargs.update({'execute':False})
scan_dim=2 #hard coded
kwargs.update({'scan_dim':scan_dim})
EA.Motors.scan(motor,start,stop,step,**kwargs)
if kwargs['debug']:
print("ScanGo scan_dim = ",scan_dim)
if kwargs["execute"]==True:
if execute == True:
#Scanning
EA.put(EAlist[1],EAparms['PassEnergy'],LensMode="Angular",Frames=EAparms['Frames'],**kwargs)
scan_go(**kwargs)
......@@ -724,6 +725,7 @@ def scanEA_map_sample(EAlist,y_list,z_list,**kwargs):
"""
kwargs.setdefault("execute",True)
kwargs.setdefault("debug",False)
kwargs.setdefault('outer_scan_dim',3)
#setting up EA
EAkwargs={
......@@ -733,18 +735,18 @@ def scanEA_map_sample(EAlist,y_list,z_list,**kwargs):
EAparms=scanEA(EAlist,**EAkwargs)
#Setting up the ScanRecord for motor scans
outer_scan_dim=3 #hard coded
inner_loop_list = y_list.insert(0,"y")
outer_loop_list = z_list.insert(0,"z")
EA.Motors.scan_2D(inner_loop_list,outer_loop_list,outer_scan_dim,**kwargs)
execute = kwargs['execute']
kwargs.update({'execute':False})
EA.Motors.scan_2D(['y',y_list[0],y_list[1],y_list[2]],['z',z_list[0],z_list[1],z_list[2]],**kwargs)
if kwargs['debug']:
print("ScanGo scan_dim = ",outer_scan_dim)
print("ScanGo scan_dim = ",kwargs['outer_scan_dim'])
if kwargs['execute']:
if execute:
#Scanning
EA.put(EAlist[1],EAparms['PassEnergy'],LensMode="Angular",Frames=EAparms['Frames'],**kwargs)
scan_go(scan_dim=outer_scan_dim)
scan_go(scan_dim=kwargs['outer_scan_dim'])
EAlog_update()
#After scan
......
iexcode/instruments/m3r_align.py
View file @ acc5252a
......@@ -146,24 +146,25 @@ def m3r_align(target_pixel=None, precision=0.5, num_tries=10, verbose=True,debug
if target_pixel == None:
target_pixel = caget( m3r_align_pv_dictionary()['desired_pixel'])
if iex.BL.branch =='d':
# coarse aligment
pixel_delta = m3r_align_stats_centroid_posx()-target_pixel
m3r_align_move_pitch_piezo(pixel_delta)
try_num = 1
while True:
if hasattr(iex.BL, 'branch'):
if iex.BL.branch =='d':
# coarse aligment
pixel_delta = m3r_align_stats_centroid_posx()-target_pixel
if debug:
print(try_num,pixel_delta)
if abs(pixel_delta) <= precision:
break
elif try_num > num_tries:
break
else:
#fine alignment
m3r_align_move_pitch_piezo(pixel_delta)
try_num = 1
while True:
pixel_delta = m3r_align_stats_centroid_posx()-target_pixel
if debug:
print('Fine align')
m3r_align_move_pitch_piezo(pixel_delta)
print(try_num,pixel_delta)
if abs(pixel_delta) <= precision:
break
elif try_num > num_tries:
break
else:
#fine alignment
if debug:
print('Fine align')
m3r_align_move_pitch_piezo(pixel_delta)
if verbose:
......
iexcode/instruments/scanRecord.py
View file @ acc5252a
......@@ -56,13 +56,27 @@ def scan_positioner_after_set(**kwargs):
def scan_positioner_settling_time_set(**kwargs):
iex.BL.mda.positioner_settling_time_set(**kwargs)
def scan_triggers_set(**kwargs):
iex.BL.mda.triggers_set(**kwargs)
def scan_triggers_set(reset=False,**kwargs):
"""
sets the trigger1 to be the lower dim execute, does not change the other triggers
does nothing if scan_dim = 1
**kwargs:
scan_dim: 1 (default)
"""
iex.BL.mda.triggers_set(reset,**kwargs)
def scan_triggers_clear(**kwargs):
"""
Clears all scan detectors triggers
**kwargs:
scan_dim: 1 (default)
"""
iex.BL.mda.triggers_clear(**kwargs)
def scan_detectors_set(**kwargs):
"""
detector_dictionary = {num:pv} #of detectors to change
"""
iex.BL.mda.detectors_set(**kwargs)
def scan_reset_all(**kwargs):
......@@ -516,7 +530,7 @@ class ScanRecord:
#detector triggers
self.triggers_clear(**kwargs)
self.triggers_set(**kwargs)
self.triggers_set(reset=True,**kwargs)
self.before_scan_set(**kwargs)
self.after_scan_set(**kwargs)
......@@ -781,7 +795,7 @@ class ScanRecord:
Clear all scan detectors triggers
**kwargs:
scan_dim: 1 (default)
verbose: Falase (default)
verbose: False (default)
Previously: Clear_Scan_Triggers
"""
......@@ -790,8 +804,15 @@ class ScanRecord:
scan_pv = self.ioc+"scan"+str(kwargs['scan_dim'])
for trig_num in range(1,5):
caput(scan_pv+'.T'+str(trig_num)+'PV',"")
#clear trigger 2-4
for trig_num in range(2,5):
caput(scan_pv+'.T'+str(trig_num)+'PV',"")
#set trigger 1
if kwargs['scan_dim']==1:
caput(scan_pv+'.T'+str(trig_num)+'PV',"")
else:
scan_trigger = self.ioc+"scan"+str(kwargs['scan_dim']-1) + ".EXSC"
caput(scan_pv+'.T'+str(trig_num)+'PV',scan_trigger)
def triggers_get(self,**kwargs):
"""
......@@ -805,31 +826,28 @@ class ScanRecord:
for tigger_num in range(1,5):
scan_trigger.update({tigger_num:caget(scan_pv+".T"+str(tigger_num)+"PV")})
def triggers_set(self,**kwargs):
def triggers_set(self,reset=False,**kwargs):
"""
sets the trigger1 to be the lower dim execute, does not change the other triggers
does nothing if scan_dim = 1
does nothing if scan_dim = None
**kwargs:
scan_dim: 1 (default)
verbose: False (default)
"""
kwargs.setdefault('scan_dim',1)
kwargs.setdefault('trigger_dictionary',self.trigger_dictionary)
kwargs.setdefault('verbose',False)
kwargs.setdefault('scan_dim',1)
scan_pv = self.ioc+"scan"+str(kwargs['scan_dim'])
#clear all triggers for scan_dim
self.triggers_clear(**kwargs)
#reset scan_dim trigger
inner_scan_trigger = self.ioc+"scan"+str(kwargs['scan_dim']-1)+".EXSC"
if kwargs['scan_dim'] > 1:
caput(scan_pv+".T1PV",inner_scan_trigger) #Resets the higher dimensional scans to trigger the inner scan
if reset:
if kwargs['scan_dim'] == 1:
kwargs.setdefault('trigger_dictionary',self.trigger_dictionary)
else:
kwargs.setdefault('trigger_dictionary',{1:inner_scan_trigger,2:"",3:"",4:""})
else:
kwargs.setdefault('trigger_dictionary',{})
#set detector triggers to dictionary values
trigger_dictionary = kwargs['trigger_dictionary']
if len(trigger_dictionary.keys())>0:
for tigger_num in trigger_dictionary.keys():
caput(scan_pv+".T"+str(tigger_num)+"PV",trigger_dictionary[tigger_num])
......@@ -870,7 +888,7 @@ class ScanRecord:
det_num=(str(i).zfill(2))
det_pv=caget(scan_pv+".D"+det_num+"PV")
if det_pv !='':
det=PV(det_pv); time.sleep(0.1)#smallest sleep to allow for PV traffic
det=PV(det_pv); time.sleep(0.25)#sleep to allow for PV traffic
if not det.connected:
print("Detector "+det_num+" has a bad PV: "+det.pvname+" not connected")
det_ok = False
......@@ -881,7 +899,7 @@ class ScanRecord:
pos_num=str(i)
pos_pv=caget(scan_pv+".P"+pos_num+"PV")
if pos_pv != '':
pos=PV(pos_pv); time.sleep(0.1)#smallest sleep to allow for PV traffic
pos=PV(pos_pv); time.sleep(0.25)#increased from .1
if not pos.connected:
print("Positioner "+pos_num+" has a BAD PV: "+pos.pvname+" not connected")
pos_ok = False
......
iexcode/instruments/xrays.py
View file @ acc5252a
......@@ -294,28 +294,34 @@ def apertures_set(c_2B=1,c_1A=1,verbose=True):
slit2B_set(Size2B[0],Size2B[1],center=(H2_center,V2_center),verbose=verbose)
def mvID(val):
def mvID(energy_eV):
"""
Sets the ID absolute set point (not optimized , mono & apertures stay fixed).
Sets the ID absolute set point in eV
Note: flux is not optimized and mono & apertures stay fixed
to go to most flux use energy
"""
ID_set_eV(val)
ID_set_eV(energy_eV)
def mvID_calibrated(val):
def mvID_calibrated(energy_eV):
"""
Sets the ID peak to val (calibrated).
Sets the ID peak to eV (calibrated).
Note: mono & apertures stay fixed
"""
ID_energy_set_eV(val)
ID_energy_set_eV(energy_eV)
def scan_ID(ID_sp_start,ID_sp_stop,ID_sp_step,**kwargs):
def scanID(ID_sp_start,ID_sp_stop,ID_sp_step,**kwargs):
"""
scan the ID set point
"""
kwargs.setdefault('execute',True)
kwargs.setdefault('scan_dim',1)
val_pv,rbv_pv = ID_scan_pvs()
scan_fillin(val_pv,rbv_pv,ID_sp_start,ID_sp_stop,ID_sp_step,**kwargs)
#Scanning
if kwargs['execute']:
scan_go(scan_dim=kwargs['scan_dim'])
def mvmono(val):
......
iexcode/macros/BL_shutdown.py
View file @ acc5252a
from epics import caget,caput
from ..instruments.electron_analyzer import EA
from ..instruments.shutters import main_shutter_close,branch_shutter_close
from ..instruments.diagnostics import diagnostics_all_out
from iexcode.instruments.electron_analyzer import EA
from iexcode.instruments.shutters import main_shutter_close,branch_shutter_close
from iexcode.instruments.diagnostics import diagnostics_all_out
from iexcode.instruments.xrays import grating, mvmono
from iexcode.instruments.ARPES import ARPES_safestate
from iexcode.instruments.kappa import kappa_safestate
from iexcode.instruments.Octupole import Octupole_safestate
##############################################################################################################
############################################ Shut down #################################################
......@@ -10,14 +14,20 @@ from ..instruments.diagnostics import diagnostics_all_out
def BL_Shutdown():
BL_CloseAllShutters()
print('\nClosed all beamline shutters')
BL_Valve_All(state="CLOSE")
print('\nClosed all beamline valves')
diagnostics_all_out()
print('\nAll diagnostics out')
try:
EA.off()
ARPES_safestate()
kappa_safestate()
Octupole_safestate()
except:
print('EA is not running, visually confirm HV is off')
caput("29iddau1:dau1:011:DAC_Set",0) #RSXS HV = "OFF"
print('Check that endstations are in a safe state')
grating('MEG')
mvmono(500)
def BL_CloseAllShutters():
main_shutter_close()
......
iexcode/macros/Kappa_optimization.py
View file @ acc5252a
from time import sleep
import iexcode.instruments.cfg as iex
from iexcode.instruments.utilities import *
from iexcode.instruments.utilities import play_sound,print_warning_message,play_sound
from iexcode.instruments.Kappa import uan
from iexcode.instruments.Motors import *
from iexcode.instruments.Motors import mvz,mvtth,dscanz,dscantth,dscanth
from iexcode.instruments.scanRecord import last_mda
from iexplot.mda_quick_plot import *
from iexplot.mda_quick_plot import fit_mda
try:
tey = iex.kappa_detectors_dictionary['tey']
d3 = iex.kappa_detectors_dictionary['d3']
d4 = iex.kappa_detectors_dictionary['d4']
mesh = iex.kappa_detectors_dictionary['mesh']
tthdet = iex.kappa_detectors_dictionary['tthdet']
detector_dictionary = iex.kappa_detectors_dictionary
tey = detector_dictionary['tey']
d3 = detector_dictionary['d3']
d4 = detector_dictionary['d4']
mesh = detector_dictionary['mesh']
tthdet = detector_dictionary['tthdet']
except:
print('you must kappa_init prior to importing')
##### Kappa Optimization scans - 2022 versions #####
def Find_kth_zero(th_0,th_180):
"""
th_0 is the motor position for specular near zero
th_180 is motor position for spectular near 180
th_0 is the motor position for specular near th=zero
th_180 is motor position for spectular near th=180
"""
Offset =0.5*(180 - th_180 - th_0)
print("with motor at 0, the actual value is ",Offset)
def align_d4(z0,th,detNum=34,plot=False):
''' th = 0 or 180
'''
def align_d4(th,**kwargs):
"""
sets tth.det to 'd4'; moves the sample out of the way
z= z0-1500
th = 0 or 180
scans tth around zero to set the user dial on d4
**kwargs:
detNum: default => 15
plot: True/False to show the plot (default => False)
sound: True/False (default=False)
"""
kwargs.setdefault('detNum',15)
kwargs.setdefault('plot',False)
kwargs.setdefault('sound',False)
#set the detector
tthdet.set('d4')
sleep(1)
d4.get_all(verbose=False)
d4gain = d4.gain
d4units = d4.current_unit
#moving sample out of the way
z0=iex.BL.motors.get('z')
uan(0,th)
mvz(z0-1500)
#setting gain for direct beam and scanning
d4.setgain(10,'uA')
dscantth(-0.25,0.25,0.01)
sleep(1)
d4.setgain(1,'nA')
#setting thing back
d4.setgain(d4gain,d4units)
mvz(z0)
#fitting
scanNum=last_mda()
tth0 =fit_mda(scanNum,detNum,'poly',plot=plot)
tth0 =fit_mda(scanNum,kwargs['detNum'],'gauss',plot=kwargs['plot'])
print('tth0 = '+str(round(tth0,3)))
play_sound()
if kwargs['sound']:
play_sound()
return tth0
def align_z0_chi0(question=True,detNum=33):
if question:
foo=input_d('Are you at th=90/270, chi=0, phi=90 (Y or N)? >')
if foo == 'Y'.lower() or foo == 'yes'.lower():
tthdet.set('d3')
sleep(1)
mvtth(0)
dscanz(-750,750,50)
sleep(1)
scanNum=last_mda()
z0=fit_mda(n,detNum,0,'erf',graph=False)
print('z0 = '+str(round(z0,0)))
play_sound()
return z0
else:
print('Canceling')
return
else:
tthdet.set('d3')
sleep(1)
def align_z0_chi0(**kwargs):
"""
sets tth.det to 'd3', and tth=0
scanz z and fits erf
**kwargs:
detNum: default => 15
plot: True/False to show the plot (default => False)
sound: True/False (default=False)
"""
kwargs.setdefault('detNum',15)
kwargs.setdefault('plot',False)
kwargs.setdefault('sound',False)
if int(iex.BL.motors.get('th')) == 0 or int(iex.BL.motors.get('th')) == 270:
tthdet.set('d3',move=False)
mvtth(0)
#scans z
dscanz(-750,750,50)
sleep(1)
n=last_mda()
fit_mda(n,25,0,'poly',graph=False)
z0=fit_z(n)
scanNum=last_mda()
z0=fit_mda(scanNum,kwargs['detNum'],'erf',plot=kwargs['plot'])
print('z0 = '+str(round(z0,0)))
play_sound()
if kwargs['sound']:
play_sound()
return z0
else:
print('th = ',round(iex.BL.motors.get('th'),0))
print_warning_message('th is not at 0 or 270')
def align_z0(omega_off):
def align_z0(omega_off,**kwargs):
"""
finds the correct z value
omega = th_offset
......@@ -87,50 +107,90 @@ def align_z0(omega_off):
- goes to grazing incidence
- scans z +/- 750 of current positions
**kwargs:
detNum: default => 15
plot: True/False to show the plot (default => False)
sound: True/False (default=False)
"""
tthdet.set('d3')
sleep(1)
kwargs.setdefault('detNum',15)
kwargs.setdefault('plot',False)
kwargs.setdefault('sound',False)
#setting the detecotr
tthdet.set('d3',move=False)
uan(0,omega_off)
#scanning
dscanz(-750,750,50)
sleep(1)
n=last_mda()
fit_mda(n,25,0,'poly',graph=False)
#z0=fit_z(n)
#print('z0 = '+str(round(z0,0)))
play_sound()
#return z0
scanNum = last_mda()
z0 = fit_mda(scanNum,kwargs['detNum'],'erf',plot=kwargs['plot'])
if kwargs['sound']:
play_sound()
return z0
def align_z180(omega_off):
tthdet.set('d3')
sleep(1)
def align_z180(omega_off,**kwargs):
""" finds the correct z value for tth = 180
omega = th_offset
Used to find the sample z0 (splitting the beam) by doing the following
- sets to d3
- goes to grazing incidence
- scans z +/- 750 of current positions
**kwargs:
detNum: default => 15
plot: True/False to show the plot (default => False)
sound: True/False (default=False)
"""
kwargs.setdefault('detNum',15)
kwargs.setdefault('plot',False)
kwargs.setdefault('sound',False)
tthdet.set('d3',move=False)
uan(0,180+omega_off)
dscanz(-750,750,50)
sleep(1)
n=last_mda()
#fit_mda(n,25,0,'poly',graph=False)
#z180=fit_z(n)
#print('z180 = '+str(round(z180,0)))
play_sound()
#return z180
scanNum=last_mda()
z180 = fit_mda(scanNum,kwargs['detNum'],'erf',plot=kwargs['plot'])
if kwargs['sound']:
play_sound()
return z180
def align_th0(omega_off,z):
def align_th0(omega_off,z,**kwargs):
"""
sets d4
goes to specular tth = 30 (default see kwargs); th = tth/2+omega_off
does rocking curve
**kwargs:
detNum: default => 15
plot: True/False to show the plot (default => False)
sound: True/False (default=False)
tth: default =>30
"""
kwargs.setdefault('detNum',15)
kwargs.setdefault('plot',False)
kwargs.setdefault('sound',False)
kwargs.setdefault('tth',30)
#set the detector
tthdet.set('d4')
sleep(1)
#uan(20,10+omega)
uan(30,15+omega_off)
d4.get_all(verbose=False)
d4gain = d4.gain
d4units = d4.current_unit
#moving sample
uan(kwargs['tth'],kwargs['tth']/2+omega_off)
mvz(z)
sleep(30)
d4.setgain(1,'nA')
#scanning
dscanth(-0.25,0.25,0.01)
n=last_mda()
fit_mda(n,25,0,'poly',graph=False)
th_spec=fit_d4(n)
#fitting
scanNum=last_mda()
th_spec = fit_mda(scanNum,kwargs['detNum'],'gauss',plot=kwargs['plot'])
print('th specular = '+str(round(th_spec,3)))
tth=round(caget('29idKappa:m9.RBV'),3)
tth=round(iex.BL.motors.get('tth'),3)
omega_off=th_spec-tth/2
print('omega = th-tth/2 = '+str(round(omega_off,3)))
play_sound()
if kwargs['sound']:
play_sound()
return omega_off
......@@ -143,10 +203,9 @@ def align_th180(omega_off,z):
d4.setgain(1,'nA')
dscanth(-0.25,0.25,0.01)
n=last_mda()
fit_mda(n,25,0,'poly',graph=False)
th_spec=fit_d4(n)
th_spec = fit_mda(n,25,0,'poly',graph=False)
print('th specular = '+str(round(th_spec,3)))
tth=round(caget('29idKappa:m9.RBV'),3)
tth=round(iex.BL.motors.get('tth'),3)
omega180=th_spec-tth/2-180
print('omega180 = th-tth/2-180 = '+str(round(omega180,3)))
play_sound()
......@@ -154,97 +213,3 @@ def align_th180(omega_off,z):
##### Kappa Optimization scans - old versions #####
def Opt_d4(iteration,moveZ,graph='y',srs=None,start=-0.5,stop=0.5,step=0.04):
tth_w=0.1
if srs==None: det=21
else: det=34
z=caget('29idKappa:m4.RBV')
current_det=caget('29idKappa:userStringSeq6.STR1',as_string=True)
if current_det != 'd4':
print('Switching detector omega')
if moveZ is not None:
mvz(z-1000)
for i in range(iteration):
scantth(start,stop,step,'relative')
scannum=last_mda()
tth0=fit_mda(scannum,det,tth_w,'gauss',graph=graph)
mvtth(tth0)
mvz(z)
return tth0,scannum
def Opt_z0(iteration,movetth,movekth,det='d4',srs=None,graph='y'):
z_w=400
current_det=caget('29idKappa:userStringSeq6.STR1',as_string=True)
if current_det != det:
print('Switching detector to '+det)
tthdet.set(det)
if det=='d3' and srs==None: det=20
if det=='d3' and srs!=None: det=33
if det=='d4' and srs==None: det=21
if det=='d4' and srs!=None: det=34
if movetth is not None:
mvtth(movetth)
if movekth is not None:
mvkth(movekth)
if iteration>1:
scanz(-2000,2000,250,'relative')
scannum=last_mda()
z0=fit_mda(scannum,det,z_w,'erf',graph=graph)
mvz(z0)
scanz(-700,700,50,'relative')
scannum=last_mda()
z0=fit_mda(scannum,det,z_w,'erf',graph=graph)
mvz(z0)
return z0,scannum
def Opt_kth(kth_0,theta,det='d4',srs=None,graph='y'):
current_det=caget('29idKappa:userStringSeq6.STR1',as_string=True)
if current_det != det:
print('Switching detector to '+det)
tthdet.set(det)
if det == 'd4':
if srs==None: det=21
else: det=34
i=1
elif det == 'd3':
if srs==None: det=20
else: det=33
i=5
mvtth(theta*2)
mvkth(kth_0+theta)
scankth(-0.5*i,0.5*i,0.05*i,'relative')
new_kth=fit_mda(last_mda(),det,0.1,'gauss',graph=graph)
mvkth(new_kth)
scankth(-0.2*i,0.2*i,0.02*i,'relative')
new_kth=fit_mda(last_mda(),det,0.1,'gauss',graph=graph)
mvkth(new_kth)
scannum=last_mda()
kth0=round(new_kth-theta,3)
print('\nkth0 = ',kth0)
print('To plot use:')
print(' fit_mda('+str(scannum)+',21,0.2,"gauss",mytitle=\'theta ='+str(theta)+'\')')
print('If you are happy with that, you can set this value as kth0 using:')
print(' kth0_set('+str(kth0)+')')
return kth0,scannum
def plot_opt(opt_scan_list,energy_list,det,mytitle=''):
fig,(ax1)=plt.subplots(1,1)
fig.set_size_inches(5,4)
for i,j in zip(opt_scan_list,energy_list):
x,y,x_name,y_name=mda_1D(i,det,1,0)
xdata = np.array(x)
ydata = np.array(y)
Imax=np.max(ydata)
ax1.plot(x,y/Imax,label=str(j)+" eV")
ax1.set(ylabel='Intensity')
ax1.set(xlabel=x_name)
ax1.ticklabel_format(style='sci', axis='y', scilimits=(0,0))
ax1.grid(color='lightgray', linestyle='-', linewidth=0.5)
ax1.legend(bbox_to_anchor=[1.2, 1],ncol=2,shadow=True, title=mytitle, fancybox=True)