refactoring

instruments/AD_utilites.py

+"""
+General functions for dealing with Area Detectors and Cameras at 29ID
+
+work in progress need to redo
+"""
+##############################################################################################################
+##############################     General Area Detector            ##############################
+##############################################################################################################
+import datetime
+import re
+from os import listdir,mkdir,chown,system,chmod
+from os.path import join, isfile, exists, dirname
+from time import sleep
+
+from epics import caget, caput
+from .IEX_endstations import *
+from files_and_folders import get_next_fileNumber
+from scanRecord import *
+
+def AD_CurrentDirectory(ADplugin):
+    """
+    returns the current directory for area detector SavePlugin
+    handles both Winodws and linux IOCs
+    ADplugin = "29idc_ps1:TIFF1:"; $(P)$(SavePlugin)
+    """
+    SubDir=caget(ADplugin+"FilePath",as_string=True)
+    if SubDir[0] == 'X':
+        Dir='/net/s29data/export/data_29idb/'
+        SubDir=SubDir.split('\\')[1:]
+    elif SubDir[0] == 'Y':
+        Dir='/net/s29data/export/data_29idc/'
+        SubDir=SubDir.split('\\')[1:]
+    elif SubDir[0] == 'Z':
+        Dir='/net/s29data/export/data_29idd/'    
+        SubDir=SubDir.split('\\')[1:]
+    else: 
+        Dir = SubDir
+        SubDir=[]
+    FilePath=join(Dir,*SubDir,'')
+    return FilePath 
+
+def AD_prefix(ADplugin):
+    """
+    returns the prefix for AreaDetector plugin based on ADplugin 
+    """
+    prefix = caget(ADplugin+"FileName_RBV",as_string=True)
+    return prefix
+
+def AD_EnableStats(ADplugin):
+    """
+    Enabling the statistics in an AreaDector
+    ADplugin = "29idc_ps1:Stats1:"; (ADplugin=$(P)$(StatsPlugin))
+    """
+    caput(ADplugin+"EnableCallbacks","Enable")
+    caput(ADplugin+"ComputeStatistics","Yes")
+    caput(ADplugin+"ComputeCentroid","Yes")
+    
+
+def AD_SaveFileSetup(ADplugin,mda,**kwargs):
+    """    
+    ADplugin = "29id_ps1:TIFF1:" which IOC and which filesaving plugin
+            (ADplugin=$(P)$(SavePlugin))
+    uses to get the current MDA directory and then set the path to one up + /dtype
+    MDA_CurrentDirectory(scanIOC=None)
+    
+    **kwargs (defaults)
+        scanIOC = BL_ioc()
+        userpath = extracted from ScanRecord unless specified
+        subfolder = taken from ADplugin unless specified
+            filepath = userpath/subfolder
+        
+        prefix = default same as subfolder
+        ext = file extension is extracted for ADplugin unless specified
+            (TIFF -> tif, HDF -> h5, ...)
+        FileTemplate="%s%s_%4.4d."+ext; format for filename first %s = filepath, second %s = prefix
+        
+    """
+    kwargs.setdefault("userpath",dirname(dirname(mda.filepath)))
+    kwargs.setdefault("subfolder",ADplugin.split(":")[-2][:-1])
+
+    kwargs.setdefault("prefix",ADplugin.split(":")[-2][:-1])
+    extDict={"TIFF":".tif","HDF":"h5"}
+    kwargs.setdefault("ext",ADplugin.split(":")[-2][:-1])
+    kwargs.setdefault("FileTemplate","%s%s_%4.4d."+kwargs["ext"])
+    
+    kwargs.setdefault("debug",False)
+    
+    if kwargs['debug']:
+        print("kwargs: ",kwargs)
+
+    fpath=join(kwargs['userpath'],kwargs['subfolder'],'')
+    print("\nFolder: " + fpath)
+    if not (exists(fpath)):
+        fileNumber=1
+    else:
+        fileNumber=get_next_fileNumber(fpath,kwargs["prefix"])
+    print("NextFile: "+str(fileNumber))
+    caput(ADplugin+"CreateDirectory",-1) #allows IOC to create directories    
+    caput(ADplugin+"FilePath",fpath)
+    caput(ADplugin+"FileName",kwargs["prefix"])
+    caput(ADplugin+"FileNumber",fileNumber)
+    
+    #setup AD
+    caput(ADplugin+"FileTemplate",kwargs["FileTemplate"])
+    caput(ADplugin+"AutoIncrement","Yes")
+    caput(ADplugin+"AutoSave","Yes")
+    
+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)
+    return Prefix
+   
+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)
+    current_run=re.findall("\d\d\d\d_\d", fpath)[0]
+    return current_run
+   
+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]
+    folder_name[:-1]
+    SubDir=caget(ADplugin+":FilePath",as_string=True)
+    if SubDir[0] == 'X':
+        current_user='Staff'
+    elif SubDir[0] == 'Y':
+        m=SubDir.find(AD_CurrentRun(ADplugin))
+        n=SubDir.find(folder_name)
+        current_user=SubDir[m+7:n]
+    else: current_user=None
+    return current_user
+
+def AD_DoneSingleSave(ADplugin,**kwargs):
+    """
+    sets and AD up ready to save images
+        Acquire -> Done
+        ImageMode -> Single
+        Save -> Enable
+        e.g. ADplugin="29id_ps2:TIFF1:"
+        
+        **kwargs:
+        P=first part of ADplugin if not specified
+        R="cam1:"; other AD have "det1:"
+    """
+    kwargs.setdefault("P",ADplugin.split(":")[0]+":")
+    kwargs.setdefault("R","cam1:")
+    caput(kwargs["P"]+kwargs["R"]+"Acquire","Done",wait=True,timeout=5*60);sleep(.1)
+    caput(kwargs["P"]+kwargs["R"]+"ImageMode","Single");sleep(.5)
+    caput(ADplugin+"EnableCallbacks","Enable");sleep(.1)
+        
+def AD_FreeRun(ADplugin,**kwargs):
+    """
+    sets and AD to disable saving and free run
+        Saving -> Disable
+        Acquire -> Done
+        ImageMode -> Single
+
+        e.g. ADplugin="29id_ps2:TIFF1:"
+        
+        **kwargs:
+        P=first part of ADplugin if not specified
+        R="cam1:"; other AD have "det1:"
+    """
+    kwargs.setdefault("P",ADplugin.split(":")[0]+":")
+    kwargs.setdefault("R","cam1:")
+    caput(kwargs["P"]+kwargs["R"]+"Acquire","Done",wait=True,timeout=5*60)
+    caput(ADplugin+"EnableCallbacks","Disable");sleep(.1)
+    caput(kwargs["P"]+kwargs["R"]+"ImageMode","Continuous");sleep(.1)
+    caput(kwargs["P"]+kwargs["R"]+"Acquire","Acquire");sleep(.1)
+        
+
+def AD_snap(ADplugin,**kwargs):
+    """
+    takes an image and save the image for ADplugin
+        e.g. ADplugin="29id_ps2:TIFF1:"
+        
+    use AD_SaveFileSetup to set filepath, prefix etc.
+    use AD_CurrentDirectory to see current directory
+    use AD_prefix to see current prefix
+    
+    **kwargs:
+        P=first part of ADplugin if not specified
+        R="cam1:"; other AD have "det1:"
+        
+        ExposureTime: changes both the exposure time and the acquire time for the snapshot
+                      resets after acquisition
+        FreeRun: True => disable setting and go back to continuous acquision 
+                 False => leave saving enabled and camera in single acquision
+    """
+    kwargs.setdefault("P",ADplugin.split(":")[0]+":")
+    kwargs.setdefault("R","cam1:")
+    kwargs.setdefault("FreeRun",True)
+
+    expT=caget(kwargs["P"]+kwargs["R"]+"AcquireTime_RBV")
+    acqT=caget(kwargs["P"]+kwargs["R"]+"AcquirePeriod_RBV")
+    
+    AD_DoneSingleSave(ADplugin,**kwargs)
+    
+    if "ExposureTime" in kwargs:
+        caput(kwargs["P"]+kwargs["R"]+"AcquireTime",kwargs["ExposureTime"])
+        caput(kwargs["P"]+kwargs["R"]+"AcquireTime",kwargs["ExposureTime"]+.01)
+        
+    caput(kwargs["P"]+kwargs["R"]+"Acquire","Acquire",wait=True,timeout=5*60)
+    
+    if "ExposureTime" in kwargs:
+        caput(kwargs["P"]+kwargs["R"]+"AcquireTime",expT)
+        caput(kwargs["P"]+kwargs["R"]+"AcquireTime",acqT) 
+
+    if kwargs["FreeRun"]:
+        sleep(.1)
+        AD_FreeRun(ADplugin,**kwargs)
+        
+
+
+def AD_ScanTrigger(ADplugin,mda,**kwargs):
+    """
+    Add Triggering of AreaDetector to scanIOC
+    ADplugin = "29idc_ps1:TIFF1:" (ADplugin=$(P)$(SavePlugin))
+    
+    **kwargs
+        scanIOC = "29id"+BL_ioc() if not specified
+        scanDIM = 1
+        P=first part of ADplugin if not specified
+        R="cam1:"; other AD have "det1"
+        detTrig = 2; detectorTrigger number
+    """
+    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"])    
+    trigger=".T"+str(kwargs["detTrig"])+"PV"
+    caput(scanPV+trigger,kwargs["P"]+kwargs["R"]+"Acquire",wait=True,timeout=5*60)
+    
+def ADplugin_ScanSetup(ADplugin,mda, **kwargs):
+    """
+    stop the acquisition, puts in ImageMode=Single
+    enables saving
+    add to detector trigger
+    Does not press go
+    
+    ADplugin = "29idc_ps1:TIFF1:"; (ADplugin=$(P)$(SavePlugin))
+    **kwargs
+        # AD_ScanTrigger
+        scanIOC = "29id"+BL_ioc() if not specified
+        scanDIM = 1
+        P=first part of ADplugin if not specified
+        R="cam1:"; other AD have "det1:"
+        detTrig = 2; detectorTrigger number
+        
+        # AD_SaveFileSetup
+        filepath=userpath (from BL_ioc scanRecord)+"/dtype"
+         (e.g. filepath="/net/s29data/export/data_29id"+folder+"/"+run+"/"+userName+"/"+df)
+        dtype = taken from ADplugin
+        FileTemplate="%s%s_%4.4d."+dtype; format for filename first %s = filepath, second %s = prefix
+        prefix = dtype by default
+    
+    """
+    #from AD_ScanTrigger
+    kwargs.setdefault("P",ADplugin.split(":")[0]+":")
+    kwargs.setdefault("R","cam1:")
+
+    AD_DoneSingleSave(ADplugin,**kwargs)
+    
+    AD_SaveFileSetup(ADplugin,**kwargs)
+    AD_ScanTrigger(ADplugin, **kwargs)
+    trigger=".T"+str(kwargs["detTrig"])+"PV"
+    scanPV=mda.ioc+"scan"+str(kwargs["scanDIM"])  
+    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+",'')")
+
+def AD_ROI_setup(AD,ROInum,xcenter=500,ycenter=500,xsize=50,ysize=50,binX=1,binY=1):  
+    """
+    AD = "29id_ps4"
+    AD = "29iddMPA"
+    """
+    # roiNUM=1  MPA_ROI_SetUp(535,539,50,50)  center of MCP
+    
+    ADplugin=AD+':ROI'+str(ROInum)+':'
+    xstart=xcenter-xsize/2.0
+    ystart=ycenter-ysize/2.0
+    caput(ADplugin+'MinX',xstart)
+    caput(ADplugin+'MinY',ystart)
+    caput(ADplugin+'SizeX',xsize)
+    caput(ADplugin+'SizeY',ysize)
+    caput(ADplugin+'BinX',binX)
+    caput(ADplugin+'BinY',binY)
+    caput(ADplugin+'EnableCallbacks','Enable')
+    print(ADplugin+' - '+caget(ADplugin+'EnableCallbacks_RBV',as_string=True))
+    #MPA_ROI_Stats(roiNUM)
+    
+def AD_OVER_SetUp(AD,ROInum,OVERnum,linewidth=5,shape='Rectangle'):
+    """
+    AD = "29id_ps4"
+    AD = "29iddMPA"
+    shape= 'Cross', 'Rectangle', 'Ellipse','Text'
+    """
+    OVER1=AD+":Over1:"+str(OVERnum)+":"
+    ROI=AD+":ROI"+str(ROInum)+":"
+    
+    caput(ROI+'EnableCallbacks','Enable')
+    caput(OVER1+"Name","ROI"+str(ROInum))
+    caput(OVER1+"Shape",shape)
+    caput(OVER1+"Red",0)
+    caput(OVER1+"Green",255)
+    caput(OVER1+"Blue",0)
+    caput(OVER1+'WidthX',linewidth)
+    caput(OVER1+'WidthY',linewidth)
+    
+    caput(OVER1+"PositionXLink.DOL",ROI+"MinX_RBV CP")
+    caput(OVER1+"SizeXLink.DOL",ROI+"SizeX_RBV CP")
+    caput(OVER1+"PositionYLink.DOL",ROI+"MinY_RBV CP")
+    caput(OVER1+"SizeYLink.DOL",ROI+"SizeY_RBV CP")
+   
+    caput(OVER1+"Use","Yes")
+
+    
+def AD_OverLayCenter(x,y,AD,OverLay=1,Num=1):
+    """
+    Sets CenterX and CenterY for  AD:Over(OverLay):Num:
+    eg. 29id_ps2:Over1:1: AD='29id_ps2',Overlay=1, Num=1
+    """
+    caput(AD+":Over"+str(OverLay)+":"+str(Num)+":CenterX",x)
+    caput(AD+":Over"+str(OverLay)+":"+str(Num)+":CenterY",y)
+    
+def AD_OverLayCenter_get(AD,OverLay=1,Num=1):
+    """
+    Sets CenterX and CenterY for  AD:Over(OverLay):Num:
+    eg. 29id_ps2:Over1:1: AD='29id_ps2',Overlay=1, Num=1
+    """
+    print('x = '+str(caget(AD+":Over"+str(OverLay)+":"+str(Num)+":CenterX")))
+    print('y = '+str(caget(AD+":Over"+str(OverLay)+":"+str(Num)+":CenterY")))
+
+def Cam_ScanSetup(mda,scanDIM,camNUM):        #individual files saving
+    pvCam="29id_ps"+str(camNUM)+":"
+    pvIOC=mda.ioc
+    #beforscan
+    Cam_SaveStrSeq(camNUM)
+    caput(pvIOC+"scan"+str(scanDIM)+".BSPV",pvIOC+"userStringSeq2.PROC")
+    #scan record (filename and trigger)
+    nextfile=str(caget(pvIOC+"saveData_baseName"))+str(caget(pvIOC+"saveData_scanNumber"))
+    filepath=caget(pvIOC+"saveData_fileSystem")
+    filepath="/net"+filepath[1:len(filepath)]+"/tif"
+    caput(pvCam+"TIFF1:FilePath",filepath)
+    caput(pvCam+"TIFF1:FileName",nextfile)
+    caput(pvCam+"TIFF1:FileWriteMode","Single")
+    caput(pvIOC+"scan1.T2PV",pvCam+"cam1:Acquire")
+    #afterscan
+    Cam_FreeStrSeq(camNUM)
+    caput(pvIOC+"scan"+str(scanDIM)+".ASPV",pvIOC+"userStringSeq1.PROC")
+    caput(pvIOC+"scan1.ASCD",1)
+    print("DON'T FORGET TO CLEAR SCAN RECORD AT THE END OF THE SCRIPT!!!!! ")
+    print("Use Script: Cam_ScanClear(scanIOC,scanDIM)")
+
+
+
+def Cam_FreeRun(camNUM):
+    camNUM=str(camNUM)
+    pv="29id_ps"+camNUM+":"
+    caput(pv+"TIFF1:AutoSave",0)
+    caput(pv+"TIFF1:EnableCallbacks",0)
+    sleep(0.5)
+    caput(pv+"cam1:ImageMode",2)
+    caput(pv+"cam1:Acquire",1)
+    caput(pv+"cam1:AcquireTime",0.015)
+    caput(pv+"cam1:AcquirePeriod",0.030)
+    print("Cam"+str(camNUM)+": Free run mode")
+
+
+
+
+
+def Cam_SaveMode(camNUM):
+    camNUM=str(camNUM)
+    pv="29id_ps"+camNUM+":"
+    caput(pv+"cam1:Acquire",0)
+    sleep(0.5)
+    caput(pv+"cam1:ImageMode",0)
+    caput(pv+"TIFF1:AutoSave",1)
+    caput(pv+"TIFF1:EnableCallbacks",1)
+    caput(pv+"TIFF1:FileNumber",1)
+    caput(pv+"cam1:AcquireTime",0.015)
+    caput(pv+"cam1:AcquirePeriod",0.030)
+    print("Cam"+str(camNUM)+": Saving mode")
+
+def Cam_Start(camNUM):
+    Cam_FreeRun(camNUM)
+
+def Cam_Stop(camNUM):
+    camNUM=str(camNUM)
+    pv="29id_ps"+camNUM+":"
+    caput(pv+"cam1:Acquire",0)
+
+
+
+def Cam_ROI_SetUp(xstart,ystart,xsize,ysize,camNUM,roiNUM=1,binX=1,binY=1):
+    pv="29id_ps"+str(camNUM)+":ROI"+str(roiNUM)+':'
+    caput(pv+'MinX',xstart)
+    caput(pv+'MinY',ystart)
+    caput(pv+'SizeX',xsize)
+    caput(pv+'SizeY',ysize)
+    caput(pv+'BinX',binX)
+    caput(pv+'BinY',binY)
+    caput(pv+'EnableCallbacks','Enable')
+    print(('ROI'+str(roiNUM)+' - '+caget(pv+'EnableCallbacks_RBV',as_string=True)))
+
+def Cam_ROI_Stats(xstart,ystart,xsize,ysize,camNUM,roiNUM=1,binX=1,binY=1):
+    pv="29id_ps"+str(camNUM)+":Stats1:"
+    Cam_ROI_SetUp(xstart,ystart,xsize,ysize,camNUM,roiNUM,binX,binY)
+    caput(pv+'EnableCallbacks','Enable')
+    roi=caget(pv+'NDArrayPort_RBV',as_string=True)
+    print((roi+' Stats => '+pv+'Total_RBV'))
+    print('To set-up as detector use:')
+    print('Cam_ROI_Det(detNUM,scanIOC,camNUM='+str(camNUM)+',roiNUM='+str(roiNUM)+')')
+    
+def Cam_ROI_Det(detNUM,scanIOC,camNUM,roiNUM,scanDIM=1):
+    pvdet='29id'+scanIOC+':scan'+str(scanDIM)+'.D'+str(detNUM)+'PV'
+    pvroi="29id_ps"+str(camNUM)+":Stats1:Total_RBV"
+    caput(pvdet,pvroi)
+    print('ROI stats set up as detector D'+str(detNUM)+' in '+scanIOC+' scan'+str(scanDIM))
+
+### Scan Record Set-Up:
+
+def Scan_Cam_Go(mda,VAL,RBV,scanDIM,start,stop,step,camNUM):
+    mda.fillin(scanDIM,VAL,RBV,start,stop,step)
+    Cam_ScanSetup(scanDIM,camNUM)
+
+    caput("29id"+scanIOC+":scan"+str(scanDIM)+".PASM","STAY")
+
+    # Scan_Go() without Before_After_Scan()
+    BL_mode=BL_Mode_Read()[0]
+    if BL_mode != 2:
+        Check_MainShutter()
+    FileName = caget("29id"+scanIOC+":saveData_baseName",as_string=True)
+    FileNum  = caget("29id"+scanIOC+":saveData_scanNumber")
+    print(FileName+": #"+str(FileNum)+" started at ", dateandtime())
+    caput("29id"+scanIOC+":scan"+str(scanDIM)+".EXSC",1,wait=True,timeout=900000)  #pushes scan button
+    print(FileName+": #"+str(FileNum)+" finished at ", dateandtime())
+
+    Cam_ScanClear(scanIOC,scanDIM)
+    caput("29id"+scanIOC+":scan"+str(scanDIM)+".PASM","PRIOR POS")
+
+def Scan_Cam_Pos2_Go(VAL1,RBV1,VAL2,RBV2,scanDIM,start1,stop1,step1,start2,stop2,camNUM):
+    scanIOC=BL_ioc()
+    pvCam="29id_ps"+str(camNUM)+":"
+    Cam_ScanSetup(scanDIM,camNUM)
+    Scan_FillIn(VAL1,RBV1,scanIOC,scanDIM,start1,stop1,step1)
+    Scan_FillIn_Pos2(VAL2,RBV2,scanIOC,scanDIM,start2,stop2)
+    caput("29id"+scanIOC+":scan"+str(scanDIM)+".PASM","STAY")
+    mda.goScan_Go(scanIOC,scanDIM=1 )
+    Cam_ScanClear(scanIOC,scanDIM)
+    caput("29id"+scanIOC+":scan"+str(scanDIM)+".PASM","PRIOR POS")
+
+def Cam_ScanClear(scanIOC,scanDIM):
+    caput("29id"+scanIOC+":scan"+str(scanDIM)+".BSPV","")
+    caput("29id"+scanIOC+":scan"+str(scanDIM)+".ASPV","")
+    caput("29id"+scanIOC+":scan"+str(scanDIM)+".T2PV","")
+    caput("29id"+scanIOC+":scan"+str(scanDIM)+".DDLY",0.5)
+    #print "Scan Record cleared from Camera Settings"
+
+
+
+### Image Acquisition:
+
+def TakeImageSetFolder(camNUM,NewFolder=""):
+    scanIOC=BL.ioc
+    pvCam="29id_ps"+str(camNUM)+":"
+    pvIOC="29id"+scanIOC+":"
+    #beforscan
+    Cam_SaveMode(camNUM)
+    #FilePath
+    filepath=caget(pvIOC+"saveData_fileSystem")+"/tif"
+    filepath="/net"+filepath[1:len(filepath)]+"/"+NewFolder
+    caput(pvCam+"TIFF1:FilePath",filepath)
+    print("WARNING: Make sure NewFolder exists !!!")
+
+def TakeImageSetFileNum(camNUM,n):
+    scanIOC=BL.mda.ioc
+    pvCam="29id_ps"+str(camNUM)+":"
+    pvIOC="29id"+scanIOC+":"
+    #beforscan
+#    Cam_SaveMode(camNUM)
+    #Filename
+    if n<10:
+        nextfile="29id"+scanIOC+"_000"+str(n)
+    elif n<100:
+        nextfile="29id"+scanIOC+"_00"+str(n)
+    elif n<1000:
+        nextfile="29id"+scanIOC+"_0"+str(n)
+    elif n<10000:
+        nextfile="29id"+scanIOC+"_"+str(n)
+    caput(pvCam+"TIFF1:FileName",nextfile)
+    caput(pvCam+"TIFF1:AutoIncrement",1)
+    caput(pvCam+"TIFF1:FileNumber",2)
+
+def TakeImageSetup(camNUM):
+    scanIOC=BL.mda.ioc
+    pvCam="29id_ps"+str(camNUM)+":"
+    pvIOC="29id"+scanIOC+":"
+    #beforscan
+    Cam_SaveMode(camNUM)
+    #Filename
+    nextfile=str(caget(pvIOC+"saveData_baseName"))+str(caget(pvIOC+"saveData_scanNumber")-1)
+    filepath=caget(pvIOC+"saveData_fileSystem")+"/tif"
+    if scanIOC == "Test":
+        filepath="/"+filepath[1:len(filepath)]        # for 29idTest
+    else:
+        filepath="/net"+filepath[1:len(filepath)]    # for 29idb/c/d
+
+    caput(pvCam+"TIFF1:FilePath",filepath)
+    caput(pvCam+"TIFF1:FileName",nextfile)
+    #afterscan
+
+def TakeImage(camNUM,AcqTime):
+    AcqPeriode=AcqTime+0.040
+    pvCam="29id_ps"+str(camNUM)+":"
+    caput(pvCam+"cam1:AcquireTime",AcqTime)
+    caput(pvCam+"cam1:AcquirePeriod",AcqPeriode)
+    sleep(1)
+    caput(pvCam+"cam1:Acquire",1,wait=True,timeout=500)
+    sleep(1)
+    TiffNum =caget(pvCam+"TIFF1:FileNumber_RBV")-1
+    FileName= caget(pvCam+"TIFF1:FileName_RBV",as_string=True)
+    print("\n"+FileName+": #"+str(TiffNum)+" started at ", dateandtime())
+    print("\n================================================\n")
+
+
+
+
+    
+    
+
+
+

instruments/FMB_mirrors.py

+from time import sleep
+
+from epics import caget, caput
+
+from .IEX_endstations import *
+from .utilities import read_dict, print_warning_message
+from .m3r import m3r_branch
+
+M0M1_fpath="/home/beams22/29IDUSER/Documents/User_Macros/Macros_29id/IEX_Dictionaries/Dict_IDCal.txt"
+
+
+def FMB_mirror_ioc(mirror_num):
+    """
+    returns the ioc name for the given mirror number: 
+    miror_num = 0 / 1 / 3  for M0 / M1 /M3R respectively
+    """
+    ioc = '29id_'+['M0','M1','','M3R'][mirror_num]+":"
+    return ioc
+
+def FMB_mirror_status(mirror_num):
+    """
+    returns the ioc name for the given mirror number: 
+    miror_num = 0 / 1 / 3  for M0 / M1 /M3R respectively
+
+    status =1 when positioned
+    """
+    ioc = '29id_'+['M0','M1','','M3R'][mirror_num]+":"
+    pv = FMB_mirror_ioc(mirror_num)
+    status=caget(pv+'SYSTEM_STS')
+
+    return status
+
+def FMB_mirror_axis_position(mirror_num,axis, verbose=True):
+    """
+    returns the readback and set point of an FMB mirror
+    """
+    pv = FMB_mirror_ioc(mirror_num)
+  
+    rbv = round(caget(pv+axis+"_MON"),3)
+    sp = round(caget(pv+axis+"_POS_SP"),3)
+
+    if verbose:
+        print (pv+axis+": "+str(rbv))
+    
+    return rbv,sp
+
+def FMB_mirror_get(mirror_num,verbose=True):
+    """
+    get and returns the current mirror position 
+    """
+    axis_labels=['Tx','Ty','Tz','Rx','Ry','Rz']
+    vals=[]
+    message = "\nM"+str(mirror_num)+" @ "
+    for axis in axis_labels:
+        rbv,sp = FMB_mirror_axis_position(mirror_num,axis, verbose=False)
+        vals.append(rbv)
+        message =+ "%.3f"+"/" % rbv 
+    if verbose:
+        print(message)
+        if mirror_num == 3:
+            mirror_branch = m3r_branch()
+            print(" => In "+mirror_branch+" branch")
+    return vals
+
+def FMB_mirror_move(mirror_num,axis,val,verbose=True):
+    """
+    miror_num = 0 / 1 / 3  for M0 / M1 /M3R respectively
+    and axis:
+        TX = lateral                 RX = Roll
+        TY = vertical                RY = Pitch
+        TZ = longitudinal            RZ = Yaw
+
+    "Previously: Move_M0M1
+    """  
+    pv = FMB_mirror_ioc(mirror_num)
+    axes = ['TX','TY','TZ','RX','RY','RZ']
+    if axis in axes:
+        caput(pv+axis+"_SP.PREC",3)  
+        caput(pv+axis+"_POS_SP")
+        caput(pv+"MOVE_CMD.PROC",1,wait=True,timeout=18000)
+        while True:
+            if FMB_mirror_status != 1:
+                sleep(.5)
+            else:
+                break
+        FMB_mirror_axis_position(mirror_num,axis,verbose=verbose)
+            
+    else:
+        print_warning_message(axis+' is not a valid axis chose '+axes)
+
+def FMB_mirror_move_all(mirror_num,position_list,verbose=True):
+    """
+    miror_num = 0 / 1 / 3  for M0 / M1 /M3R respectively
+    position_list = [TX,TY,TZ,RX,RY,RZ]
+
+    "Previously: Move_M0M1
+    """  
+    pv = FMB_mirror_ioc(mirror_num)
+    for axis in ['TX','TY','TZ','RX','RY','RZ']:
+        caput(pv+axis+"_SP.PREC",3)  
+        caput(pv+axis+"_POS_SP")
+    caput(pv+"MOVE_CMD.PROC",1,wait=True,timeout=18000)
+    while True:
+        if FMB_mirror_status != 1:
+            sleep(.5)
+        else:
+            break
+    FMB_mirror_axis_position(mirror_num,axis,verbose=verbose)
+    
+
+def FMB_mirror_tweak(mirror_num,axis,val,verbose=False):
+    """
+    miror_num = 0 / 1 / 3  for M0 / M1 /M3R respectively
+    and axis:
+        TX = lateral                 RX = Roll
+        TY = vertical                RY = Pitch
+        TZ = longitudinal            RZ = Yaw
+
+    "Previously: Move_M0M1
+    """
+    pv = FMB_mirror_ioc(mirror_num)
+    previous_position = FMB_mirror_axis_position(mirror_num,axis,verbose=False)
+    new_position = previous_position[0] + val
+    FMB_mirror_move(mirror_num,axis,new_position,verbose=False)
+
+    if verbose:
+     print(pv+" "+str(previous_position[0])+" -> "+str(new_position[0]))
+
+
+def FMB_mirror_scan(mirror_num,start,stop,step):
+    """
+    e.g. Scan_FMB_mirror("m1:TX",-5,5,.25)
+
+        TX = lateral                 RX = Yaw
+        TY = vertical                RY = Pitch
+        TZ = longitudinal            RZ = Roll
+    """
+    pv = FMB_mirror_ioc(mirror_num)
+    
+    # database sets .PREC==0.  We want more digits than that.
+    caput(pv+"_SP.PREC",3)
+    mda.fillin(pv+"_SP",pv+"_MON",start,stop,step)
+
+
+
+
+def M0M1_table(run,mirror):
+    """
+    Prints the positions TX / TY / TZ / RX / RY / RZ for either Mirror = 0 or 1 (M0 or M1) for the specified Run
+    Run='default' give a reasonable starting position after homing
+    M0M1_SP() will put those values as the set points, you will need to push the Move button
+    """
+    M0M1_Pos=read_dict(M0M1_fpath)
+    return M0M1_Pos[run][mirror[-1]]
+
+def M0M1_SP(run,mirror,Go=False):
+    """
+    Gets the values from the mirror position from a previous run and
+    put values as defined by M0M1_Table as the set points
+    Go = True / False (default); True to moves to that position
+    """
+    mirror_pos=M0M1_table(run,mirror).split('/')
+    motor=['TX','TY','TZ','RX','RY','RZ']
+    mirror=mirror[-1]
+    for i in range(len(motor)):
+        PV="29id_m"+str(mirror)+":"+motor[i]+"_POS_SP"
+        Val=mirror_pos[i] #float(MirrorPos[i])
+        print(PV+" = "+Val)
+        caput(PV,Val)
+    sleep(0.5)
+    if Go:
+        caput('29id_m'+str(mirror)+':MOVE_CMD.PROC',0,wait=True,timeout=18000)
+    else:
+        print(" caput(\'29id_m"+str(mirror)+":MOVE_CMD.PROC\',0)")

instruments/IEX_VPU.py

+
+
+from math import *
+from time import sleep
+import numpy.polynomial.polynomial as poly
+
+from epics import caget, caput
+
+from .utilities import dateandtime, print_warning_message, read_dict
+from .shutters import main_shutter_check_open
+from .VLS_PGM import mono_grating_get
+from .userCalcs import userCalcOut_clear
+
+IDcal_fpath="/home/beams22/29IDUSER/Documents/User_Macros/Macros_29id/IEX_Dictionaries/Dict_IDCal.txt"
+
+##############################################################################################################
+################################             ID limits and calibration             ##############################
+##############################################################################################################
+
+def ID_calc_SP(mono_grating,ID_mode,hv_eV):    # Mode = state (0=RCP,1=LCP,2=V,3=H)
+    """Calculate the ID SP for a given polarization mode and energy;
+    with Mode = 0 (RCP),1 (LCP), 2 (V), 3 (H)
+    
+    Previously: ID_Calc
+    """
+    
+    
+    if type(ID_mode)== str:
+        ID_state=ID_state_mode()[ID_mode]
+    try:
+        K=ID_Coef(mono_grating,ID_mode,hv_eV)
+        ID=poly.polyval(hv_eV,K)
+    except KeyError:
+        message_string='Not a valid ID mode!'+"\nValid Modes: "+str(ID_mode_list)
+        print_warning_message(message_string)
+        ID=caget("ID29:EnergySeteV")
+    return round(ID,1)
+
+##############################################################################################################
+################################             ID Functions             ##############################
+##############################################################################################################
+
+
+def ID_wait_for_permission():
+    """
+    Monitors the ID permissions and waits for the ID to be in User Mode and then breaks
+    Checks the status every 30 seconds
+
+    Previously: WaitForPermission
+    """
+    while True:
+        ID_Access=caget("ID29:AccessSecurity.VAL")
+        if (ID_Access!=0):
+            print("Checking ID permission, please wait..."+dateandtime())
+            sleep(30)
+        else:
+            print("ID now in user mode -"+dateandtime())
+            break
+
+##############################################################################################################
+def ID_get_all(verbose=False):
+    """ 
+    returns dictionart with: ID_Mode, ID_QP_ratio, ID_SP, ID_RBV
+
+    Previously: Get_energy
+    """
+    vals={    
+        "ID_mode":ID_mode_get(verbose=False),
+        "ID_QP_ratio":ID_QP_ratio_get(verbose=False),
+        "ID_sp":ID_SP_get_eV(verbose=False),
+        "ID_rbv":ID_rbv_get_eV(verbose=False)
+    }
+
+    if verbose:
+        print(" ID SP  : "+"%.2f" % vals["ID_sp"] , "eV    ID mode : "+vals["ID_mode"])
+        print(" ID RBV : "+"%.2f" % vals['ID_rbv'], "eV    QP mode : "+str(vals['ID_QP_ratio']) +" %")
+    return vals
+
+def ID_energy_range(ID_mode=None):  
+    """
+    Returns the ID_min_SP, ID_max_SP for a given ID mode
+    
+    Previously: ID_Range
+    """
+    if ID_mode == None:
+        ID_state=caget("ID29:ActualMode")
+    else:
+        ID_state = ID_state_mode(ID_mode)
+    #       RCP,LCP, V , H , HN
+    ID_min_SP = [400,400,440,250,250]
+    ID_max_SP = 3800
+    #hv_min=[390,390,430,245,245]
+
+    return ID_min_SP[ID_state],ID_max_SP
+
+def ID_mode_list():
+    """
+    returns the ID_mode_List
+    current mode = 
+    """
+    return ["RCP", "LCP", "V", "H", "HN"]
+
+
+def ID_state_mode(ID_mode):
+    """ 
+    Returns the state number if mode is a string 
+    Returns the mode string if mode is an int or float
+
+    Previously: ID_State2Mode(which,mode)
+    """
+    try:
+        ID_mode = ID_mode.upper()
+        ID_state =ID_mode_list()[ID_mode].index(ID_mode)
+        return ID_state
+
+    except:
+        message_string='Not a valid ID mode!'+"\nValid Modes: "+str(ID_mode_list)
+        print_warning_message(message_string)
+
+def ID_state_get():
+    """
+    Returns the current ID state
+    """
+    ID_state = caget("ID29:ActualMode")
+
+def ID_mode_get(verbose=True):
+    """
+    Returns the current ID mode
+    """
+    ID_state = caget("ID29:ActualMode")
+    ID_mode = ID_mode_list()[ID_state]
+
+    if verbose:
+        print('ID mode: '+ID_mode)
+    return ID_mode
+
+def ID_mode_put(ID_mode):
+    """
+    writes the desired mode to the correct ID pv
+    """
+    caput("ID29:DesiredMode.VAL",ID_mode,wait=True,timeout=18000)     # RCP
+
+
+def ID_ready(verbose=True):
+    """
+    check both the read and busy pv every 2 seconds
+
+    Previously: ID_Ready
+    """
+    while True:
+        RBV=caget("ID29:Energy.VAL")
+        checkready=caget("ID29:feedback.VAL")
+        checkbusy=caget("ID29:BusyRecord")
+        if (checkready!="Ready") or (RBV < 3.7):
+            sleep(2)
+        elif ((checkready=="Ready") and (RBV > 3.7)) and (checkbusy=="Busy"):
+            caput("ID29:Busy.VAL",0)
+        else:
+            break
+    if verbose:
+        print("ID Ready")       
+
+def ID_stop(verbose=True):
+    """
+    waits for the ID to be in user mode and then turns it off
+
+    Previously: ID_Stop
+    """
+    ID_wait_for_permission()
+    caput("ID29:Main_on_off.VAL",1,wait=True,timeout=18000)
+
+    sleep(5)
+    if verbose:
+        print("ID is now off")
+
+
+def ID_off():
+    """
+    calls ID_stop
+    """
+    ID_stop()
+
+
+def ID_start(ID_mode='RCP'):
+    """
+    waits for ID permission and then 
+    starts ID with a specific polarization
+    mode = \"H\", \"V\", \"RCP\" or \"LCP\"
+
+    """
+    
+    ID_wait_for_permission()
+    #turning on the ID; default mode = RCP
+    print("Starting ID  -  "+dateandtime())
+    caput("ID29:EnergySet.VAL",3.8)
+    caput("ID29:Main_on_off.VAL",0,wait=True,timeout=18000)
+    ID_ready()
+    ID_mode_put(ID_mode)
+   
+    ID_wait_for_permission()
+    main_shutter_check_open()
+
+    print('ID is now on, please set your energy')
+
+def ID_switch_mode(ID_mode):
+    """
+    Change ID polarization; which = \"H\", \"V\", \"RCP\" or \"LCP\"
+        if ID_mode = current mode then does nothing
+
+    WARNING: Does not set the energy
+
+    Previously Switch_IDMode
+    """
+    ID_wait_for_permission()
+    main_shutter_check_open()
+    ID_state = ID_state_get()
+
+    try:
+        if ID_state_mode(ID_mode) != ID_state:
+            print("Turning ID off...")
+            ID_stop(verbose=True)
+            sleep(10)
+            
+            print("Switching ID mode, please wait...")
+            ID_mode_put(ID_mode) 
+            ID_ready()
+        print("ID Mode:",ID_mode)
+    except:
+        print_warning_message("Not a valid ID mode")
+
+def ID_SP_get(verbose=False):
+    """
+    returns the ID setpoint in keV
+    """
+    ID_SP = caget("ID29:EnergyScanSet.VAL")
+    if verbose:
+        print("ID_SP: ", ID_SP)
+    return ID_SP
+
+def ID_SP_get_eV(verbose=False):
+    """
+    returns the ID setpoint in eV
+    """
+    ID_SP = ID_SP_get(verbose=False)
+    if verbose:
+        print("ID_SP: ", ID_SP)
+    return ID_SP
+
+def ID_rbv_get(verbose=False):
+    """
+    returns the readback value for the 
+    """
+    ID_RBV = caget("ID29:EnergyRBV")
+    if verbose:
+        print("ID_RBV: ", ID_RBV)
+    return ID_RBV  
+
+def ID_rbv_get_eV(verbose=False):
+    """
+    returns the readback value for the 
+    """
+    ID_RBV = ID_rbv_get(verbose=False)
+    if verbose:
+        print("ID_RBV: ", ID_RBV)
+    return ID_RBV 
+
+def ID_restart():
+    """
+    turns off the ID and turns it back on with the same set point
+    
+    Previously: ID_Restart
+    """
+    ID_mode=ID_mode_list[ID_state_get()]
+    hv=ID_SP_get_eV(verbose=False)
+
+    print("Restarting ID", dateandtime())
+    ID_stop(verbose=False)
+    ID_ready(verbose=False)
+  
+    ID_start(verbose=False)
+    
+    print("ID is back ON", dateandtime())
+    ID_SP_put(hv)
+
+def ID_SP_put(hv_eV,verbose=True):
+    """
+    "Sets the ID set point to a specific value (hv(eV)) which will not likely be optimum"
+
+    Previously: SetID_Raw
+    """
+    ID_wait_for_permission()
+    main_shutter_check_open()
+
+    ID_max,ID_min = ID_energy_range()
+    ID_SP=min(max(hv_eV,ID_min),ID_max)*1.0
+    
+    if hv_eV < ID_min or hv_eV > ID_max:
+        message_string="Set point out of BL energy range \nPlease select a different energy."
+        print_warning_message(message_string)
+    else:
+        ID_SP_RBV=round(caget("ID29:EnergySet.VAL"),3)*1000
+        if ID_SP == ID_SP_RBV:                # checks if ID is already at the SP energy
+            ID_RBV=caget("ID29:EnergyRBV")
+            if verbose:
+                print("ID SET : "+"%.1f" % ID_SP, "eV")
+                print("ID RBV : "+"%.1f" % ID_RBV, "eV")
+                print(caget('ID29:TableDirection',as_string=True))
+        else:
+            caput("ID29:EnergyScanSet.VAL",ID_SP/1000,wait=True,timeout=18000)
+            sleep(0.5)
+            caput("ID29:EnergyScanSet.VAL",(ID_SP+0.001)/1000,wait=True,timeout=18000)
+            caput('ID29:StartRamp.VAL',1)
+            
+            sleep(5)
+            ID_ready(verbose=False)
+            
+            ID_RBV=caget("ID29:EnergyRBV")
+            print("\nID SET : "+"%.1f" % ID_SP, "eV") 
+            print("ID RBV : "+"%.1f" % ID_RBV, "eV")
+            print(caget('ID29:TableDirection',as_string=True))
+
+            ID_diff = abs(ID_RBV-ID_SP)
+            ID_bw = ID_SP*0.095
+            if ID_diff > ID_bw:
+                sleep(20)
+                ID_RBV=caget("ID29:EnergyRBV")
+                ID_diff = abs(ID_RBV-ID_SP)
+                print("\nID RBV : "+"%.1f" % ID_RBV, "eV")
+                if ID_diff > ID_bw:
+                    ID_restart
+
+def ID_energy_set(hv_eV):
+    """
+    Sets optimum ID set point for hv(eV) (max intensity)
+    and opens the main shutter
+
+    Note that QP is generally not calibrated
+
+    Previously: SetID
+    """
+    ID_mode = ID_mode_list[ID_state_get()]
+    mono_grating = mono_grating_get()
+    ID_SP = ID_calc_SP(mono_grating,ID_mode,hv_eV)
+   
+    ID_SP(hv_eV)
+
+    
+def ID_QP_ratio_get(verbose=True):
+    """
+    gets the read back for the QP ratio
+    calculate the QP ratio
+
+    
+    """
+    Byq=caget("ID29:ByqRdbk")
+    Vcoil=caget("ID29:ByRdbk.VAL")
+    ratio_calc=Byq/Vcoil
+    ratio_RBV=caget("ID29:QuasiRatio.RVAL")
+    if verbose:
+        print("QP ratio =", round(ratio_RBV,3)*100,"%")
+        print("QP RBV   =", ratio_RBV,"%")
+
+    if abs(ratio_RBV-ratio_calc)>1:
+        message_string="QP RBV and QP calc do not agree \nCheck Interlocks"
+        print_warning_message(message_string)
+    return ratio_RBV, ratio_calc
+
+def ID_QP_mode_set(ID_mode,QP_ratio):
+    """
+    switched to QP mode, if not currently in
+    sets the QP ratio (QP ration min is 70) and polarization (ID mode)
+
+    does not set the energy, need to use ID_SP
+
+    Previously: Switch_IDQP
+    """
+    QP_min = 70
+    if QP_ratio < QP_min:
+        message_string="QP ratio is too small, setting it to minimum allowed value ("+str(QP_min)+")"
+        print_warning_message(message_string)
+    QP_ratio=max(70,QP_ratio)
+    
+    ratio_RBV, ratio_calc=ID_QP_ratio_get()
+    if ratio_RBV != QP_ratio:
+        ID_stop()
+        caput("ID29:QuasiRatioIn.C",QP_ratio)
+    
+    ID_start(ID_mode)
+    sleep(15)
+
+    Byq=caget("ID29:ByqRdbk")
+    Vcoil=caget("ID29:ByRdbk.VAL")
+    ratio=Byq/Vcoil
+    ratio_RBV=caget("ID29:QuasiRatio.RVAL")
+    print("QP ratio =", round(ratio,3)*100,"%")
+    print("QP RBV   =", ratio_RBV,"%")
+    sleep(15)
+
+def ID_Coef(grt,ID_mode,hv_eV):    # Mode = state (0=RCP,1=LCP,2=V,3=H); 
+    
+    """Return the ID coeff for a given polarization mode and energy;
+    with Mode = 0 (RCP),1 (LCP), 2 (V), 3 (H).
+    Current coefficient dictionary:
+        /home/beams22/29IDUSER/Documents/User_Macros/Macros_29id/IEX_Dictionaries/Dict_IDCal.txt
+    
+    Previously: ID_Coef
+    """
+
+    def ListRange(grt,ID_mode,IDdict):  # extract the list of break pts for a given mode/grt 
+        tmp_list=[]
+        for item in (IDdict[grt][ID_mode]):
+            tmp_list.append(item[0])  
+        return tmp_list
+
+
+    def FindRange(hv_eV,range_list):         # returns the index for the corresponding range
+        B = [x - hv_eV for x in range_list]
+        #print(B)
+        index = [i for (i, x) in enumerate(B) if x > 0]
+        #print(index)
+        return(index[0])
+    
+    try:
+        ID_function=read_dict(IDcal_fpath)
+    
+    except KeyError:
+        print("Unable to read dictionary") 
+        
+    try:   
+        Lrange = ListRange(grt,ID_mode,ID_function)
+        Erange = FindRange(hv_eV,Lrange)
+        K = ID_function[grt][ID_mode][Erange][1]
+        return K
+        
+    except KeyError:
+        print("WARNING: PLease select one of the following:")
+        print("        mode 0 = RCP")
+        print("        mode 1 = LCP")
+        print("        mode 2 = V")
+        print("        mode 3 = H")
+
+def ID_scan_pvs():
+    """ 
+    returns the rbv and val for scanning 
+    """
+    val_pv="ID29:EnergyScanSeteV"
+    rbv_pv="ID29:EnergySetRBV"
+    return val_pv, rbv_pv
+
+def ID_scan_fillin(mda,scan_dim,start,stop,step,**kwargs):
+    """
+    fills in the scanRecord for scanning the ID set point
+
+    **kwargs => scanRecord.fillin kwargs
+    """
+    #Setting up the ScanRecord for ID in Table mode
+    val_pv, rbv_pv = ID_scan_pvs()
+    mda.fillin(scan_dim,val_pv,rbv_pv,start,stop,step,**kwargs)
+
+
+def ID_scan_fillin_table(mda,scan_dim,ID_array,**kwargs):
+    """
+    fills in the scanRecord for scanning the ID set point
+
+    **kwargs => scanRecord.fillin kwargs
+    """
+    #Setting up the ScanRecord for ID in Table mode
+    val_pv, rbv_pv = ID_scan_pvs()
+    mda.fillin.table(scan_dim,val_pv,rbv_pv,ID_array,**kwargs)
+
+##############################################################################################################
+##############################             ID direction table        ##############################
+##############################################################################################################
+
+def ID_Table():
+    """
+
+    Previously: ID_Table
+    """
+    table = caget("ID29:TableDirection")    # up = 1 , down = 0
+    By = caget("ID29:ByPolaritySet")        # pos = 1, neg = 0
+
+    mode = ID_mode_get()
+
+    if By > 0:
+        print("\nBy > 0")
+        if table == 1:            # By=1, table = 1     => -1    => A=B
+            print("table = up")
+            ID_direction = -1
+        elif table ==  0:        # By=1, table = 0     => +1    => A#B
+            print("table = down")
+            ID_direction = 1
+    elif By <= 0:
+        print("\nBy < 0")
+        if table == 1:            # By=0, table = 1     => +1    => A=B
+            print("table = up")
+            ID_direction = 1
+        elif table ==  0:        # By=0, table = 0     => -1    => A#B
+            print("table = down")
+            ID_direction = -1
+
+
+    if mode == "H" and mode == "RCP":
+        if By > 0 and table == 0:
+            print_warning_message("will do a long hysteresis if decreasing energy !!!")
+#    if Mode == "HN" and Mode == "LCP":
+#        if By = 0 and table == 1:
+#            print "WARNING: will do a long hysteresis if decreasing energy !!!"
+    print("ID direction", ID_direction)
+    return ID_direction
+
+def ID_table_userCalcOut(): 
+    """
+       # Work in progress
+
+    Previously:ID_Table_CalcOut
+    """
+    n=4
+    userCalcOut_clear("b",n)
+    pvstr="29idb:userCalcOut"+str(n)
+    caput(pvstr+".DESC","ID_Table")
+    table="ID29:TableDirection"    # up = 1 , down = 0
+    By="ID29:ByPolaritySet"        # pos = 1, neg = 0
+    caput(pvstr+".INPA",table+" CP NMS")
+    caput(pvstr+".INPB",By+" CP NMS")
+    caput(pvstr+".CALC$","A#B")
+    caput(pvstr+".OOPT","On Change")
+    caput(pvstr+".DOPT","Use CALC")
\ No newline at end of file

instruments/IEX_dictionaries/Dict_Slit.txt

+
+======= 20201231: 
+{'HEG':{'S1H':0,'S1V':0.25,'S2H':0,'S2V':-0.75},'MEG':{'S1H':0.25,'S1V':0.125,'S2H':-0.5,'S2V':0.5}}
+
+======= 20210101: 
+{'HEG': {'S1H':0,'S1V':0.25,'S2H':0,'S2V':-0.75}, 'MEG': {'S1H':0.25,'S1V':0.125,'S2H':-0.5,'S2V':0.5}}
+
+======= 20210126: 
+{'HEG': {'S1H': 0, 'S1V': 0.25, 'S2H': 0, 'S2V': -0.75}, 'MEG': {'S1H': 0.0, 'S1V': 0.0, 'S2H': 0.25, 'S2V': 0.0}}
+
+======= 20210127: 
+{'HEG': {'S1H': 0.0, 'S1V': 0.0, 'S2H': 0.0, 'S2V': 0.0}, 'MEG': {'S1H': 0.0, 'S1V': 0.0, 'S2H': 0.25, 'S2V': 0.0}}
+
+======= 20210209: 
+{'HEG': {'S1H': 0.0, 'S1V': 0.0, 'S2H': 0.0, 'S2V': 0.0}, 'MEG': {'S1H': 0.0, 'S1V': 0.0, 'S2H': 0.0, 'S2V': 0.0}}
+
+======= 20210316: 
+{'HEG': {'S1H': 0.0, 'S1V': 0.0, 'S2H': 0.0, 'S2V': 0.0}, 'MEG': {'S1H': -0.125, 'S1V': 0.0, 'S2H': -0.325, 'S2V': -0.25}}
+
+======= 20210323: 
+{'HEG': {'S1H': 0.0, 'S1V': 0.0, 'S2H': 0.0, 'S2V': 0.0}, 'MEG': {'S1H': -0.5, 'S1V': 0.0, 'S2H': -0.25, 'S2V': 0.0}}
+
+======= 20210330: 
+{'HEG': {'S1H': 0.0, 'S1V': 0.0, 'S2H': 0.0, 'S2V': 0.0}, 'MEG': {'S1H': -0.125, 'S1V': 0.0, 'S2H': -0.25, 'S2V': 0.0}}
+
+======= 20210406: 
+{'HEG': {'S1H': 0.0, 'S1V': 0.0, 'S2H': 0.0, 'S2V': 0.0}, 'MEG': {'S1H': 0.0, 'S1V': 0.0, 'S2H': 0.0, 'S2V': 0.0}}
+
+======= 20210408: 
+{'HEG': {'S1H': 0.0, 'S1V': 0.0, 'S2H': 0.0, 'S2V': 0.0}, 'MEG': {'S1H': -0.125, 'S1V': -0.125, 'S2H': 0.0, 'S2V': -0.25}}
+
+======= 20210413: 
+{'HEG': {'S1H': 0.0, 'S1V': 0.0, 'S2H': 0.0, 'S2V': 0.0}, 'MEG': {'S1H': -0.125, 'S1V': -0.1, 'S2H': -0.25, 'S2V': 0.0}}
+
+======= 20210413: 
+{'HEG': {'S1H': 0.0, 'S1V': 0.0, 'S2H': 0.0, 'S2V': 0.0}, 'MEG': {'S1H': 0.125, 'S1V': -0.1, 'S2H': 0.25, 'S2V': -0.25}}
+
+======= 20210603: 
+{'HEG': {'S1H': 0.0, 'S1V': 0.0, 'S2H': 0.0, 'S2V': 0.0}, 'MEG': {'S1H': 0.125, 'S1V': 0.0, 'S2H': 0.5, 'S2V': 0.0}}
+
+======= 20210603: 
+{'HEG': {'S1H': 0.125, 'S1V': 0.0, 'S2H': 0.75, 'S2V': -1.0}, 'MEG': {'S1H': 0.125, 'S1V': 0.0, 'S2H': 0.5, 'S2V': 0.0}}
+
+======= 20210603: 
+{'HEG': {'S1H': 0.125, 'S1V': -0.125, 'S2H': 0.75, 'S2V': -1.0}, 'MEG': {'S1H': 0.125, 'S1V': 0.0, 'S2H': 0.5, 'S2V': 0.0}}
+
+======= 20210622: 
+{'HEG': {'S1H': 0.125, 'S1V': -0.125, 'S2H': 0.75, 'S2V': -1.0}, 'MEG': {'S1H': 0.0, 'S1V': -0.125, 'S2H': -0.25, 'S2V': -0.25}}
+
+======= 20210629: 
+{'HEG': {'S1H': 0.125, 'S1V': -0.125, 'S2H': 0.75, 'S2V': -1.0}, 'MEG': {'S1H': 0.0, 'S1V': 0.0, 'S2H': -0.5, 'S2V': 0.25}}
+
+======= 20210708: 
+{'HEG': {'S1H': 0.125, 'S1V': -0.125, 'S2H': 0.75, 'S2V': -1.0}, 'MEG': {'S1H': 0.125, 'S1V': -0.125, 'S2H': -0.25, 'S2V': -0.25}}
+
+======= 20210714: 
+{'HEG': {'S1H': 0.125, 'S1V': -0.125, 'S2H': 0.75, 'S2V': -1.0}, 'MEG': {'S1H': 0.125, 'S1V': 0.0, 'S2H': -0.75, 'S2V': 0.0}}
+
+======= 20210727: 
+{'HEG': {'S1H': 0.125, 'S1V': -0.125, 'S2H': 0.75, 'S2V': -1.0}, 'MEG': {'S1H': 0.125, 'S1V': -0.125, 'S2H': -0.25, 'S2V': -0.25}}
+
+======= 20210914: 
+{'HEG': {'S1H': 0.25, 'S1V': 0.0, 'S2H': -0.25, 'S2V': 0.0}, 'MEG': {'S1H': 0.125, 'S1V': -0.125, 'S2H': -0.25, 'S2V': -0.25}}
+
+======= 20210914: 
+{'HEG': {'S1H': 0.25, 'S1V': 0.0, 'S2H': -0.25, 'S2V': 0.0}, 'MEG': {'S1H': 0.125, 'S1V': 0.0, 'S2H': -0.25, 'S2V': 0.0}}
+
+======= 20211004: 
+{'HEG': {'S1H': 0.25, 'S1V': 0.0, 'S2H': -0.25, 'S2V': 0.0}, 'MEG': {'S1H': 0.125, 'S1V': 0.0, 'S2H': -0.5, 'S2V': 0.0}}
+
+======= 20211019: 
+{'HEG': {'S1H': 0.25, 'S1V': 0.0, 'S2H': -0.25, 'S2V': 0.0}, 'MEG': {'S1H': 0.125, 'S1V': 0.0, 'S2H': -0.75, 'S2V': 0.0}}
+
+======= 20211103: 
+{'HEG': {'S1H': 0.25, 'S1V': 0.0, 'S2H': -0.25, 'S2V': 0.0}, 'MEG': {'S1H': 0.0, 'S1V': -0.125, 'S2H': -0.75, 'S2V': 0.0}}
+
+======= 20211109: 
+{'HEG': {'S1H': -0.125, 'S1V': 0.0, 'S2H': -0.25, 'S2V': -0.25}, 'MEG': {'S1H': 0.0, 'S1V': -0.125, 'S2H': -0.75, 'S2V': 0.0}}
+
+======= 20211117: 
+{'HEG': {'S1H': -0.125, 'S1V': 0.0, 'S2H': -0.25, 'S2V': -0.25}, 'MEG': {'S1H': 0.125, 'S1V': 0.0, 'S2H': -0.25, 'S2V': 0.0}}
+
+======= 20211117: 
+{'HEG': {'S1H': -0.125, 'S1V': 0.0, 'S2H': -0.25, 'S2V': -0.25}, 'MEG': {'S1H': 0.125, 'S1V': 0.0, 'S2H': -0.25, 'S2V': 0.0}}
+
+======= 20211214: 
+{'HEG': {'S1H': -0.125, 'S1V': 0.0, 'S2H': -0.25, 'S2V': -0.25}, 'MEG': {'S1H': 0.25, 'S1V': 0.0, 'S2H': -0.75, 'S2V': 0.0}}
+
+======= 20220202: 
+{'HEG': {'S1H': 0.25, 'S1V': 0.0, 'S2H': 0.25, 'S2V': 0.25}, 'MEG': {'S1H': 0.25, 'S1V': 0.0, 'S2H': -0.75, 'S2V': 0.0}}
+
+======= 20220202: 
+{'HEG': {'S1H': 0.25, 'S1V': 0.0, 'S2H': 0.25, 'S2V': 0.25}, 'MEG': {'S1H': 0.0, 'S1V': 0.0, 'S2H': 0.0, 'S2V': 0.0}}
+
+======= 20220209: 
+{'HEG': {'S1H': 0.25, 'S1V': 0.0, 'S2H': 0.25, 'S2V': 0.25}, 'MEG': {'S1H': 0.125, 'S1V': -0.125, 'S2H': -0.25, 'S2V': 0.0}}
+
+======= 20220215: 
+{'HEG': {'S1H': 0.25, 'S1V': 0.0, 'S2H': 0.25, 'S2V': 0.25}, 'MEG': {'S1H': 0.125, 'S1V': 0.0, 'S2H': -0.25, 'S2V': 0.125}}
+
+======= 20220215: 
+{'HEG': {'S1H': 0.25, 'S1V': 0.0, 'S2H': 0.25, 'S2V': 0.25}, 'MEG': {'S1H': 0.125, 'S1V': 0.0, 'S2H': -0.25, 'S2V': 0.0}}
+
+======= 20220222: 
+{'HEG': {'S1H': 0.25, 'S1V': 0.0, 'S2H': 0.25, 'S2V': 0.25}, 'MEG': {'S1H': 0.1, 'S1V': 0.0, 'S2H': 0.22, 'S2V': 0.0}}
+
+======= 20220228: 
+{'HEG': {'S2V': 0.25, 'S2H': 0.25, 'S1V': 0.0,'S1H': 0.25 }, 'MEG': {'S2V': 0.0, 'S2H': 0.22, 'S1V': 0.0, 'S1H': 0.1}}
+
+======= 20220302: 
+{'HEG': {'S2V': 0.25, 'S2H': 0.25, 'S1V': 0.0, 'S1H': 0.25}, 'MEG': {'S2V': 0.0, 'S2H': 0.352, 'S1V': -0.132, 'S1H': 0.12}}
+
+======= 20220302: 
+{'HEG': {'S2V': -0.3, 'S2H': 0.2, 'S1V': 0.0, 'S1H': 0.0}, 'MEG': {'S2V': 0.0, 'S2H': 0.352, 'S1V': -0.132, 'S1H': 0.12}}
+
+======= 20220308: 
+{'HEG': {'S2V': -0.3, 'S2H': 0.2, 'S1V': 0.0, 'S1H': 0.0}, 'MEG': {'S2V': 0.2, 'S2H': 0.2, 'S1V': -0.08, 'S1H': 0.16}}
+
+======= 20220310: 
+{'HEG': {'S2V': -0.3, 'S2H': 0.2, 'S1V': 0.0, 'S1H': 0.0}, 'MEG': {'S2V': 0.133, 'S2H': 0.245, 'S1V': -0.105, 'S1H': 0.17}}
+
+======= 20220310: 
+{'HEG': {'S2V': -0.3, 'S2H': 0.2, 'S1V': 0.0, 'S1H': 0.0}, 'MEG': {'S2V': 0.107, 'S2H': 0.109, 'S1V': -0.105, 'S1H': 0.166}}
+
+======= 20220316: 
+{'HEG': {'S2V': -0.3, 'S2H': 0.2, 'S1V': 0.0, 'S1H': 0.0}, 'MEG': {'S2V': 0.0, 'S2H': 0.284, 'S1V': -0.055, 'S1H': 0.156}}
+
+======= 20220322: 
+{'HEG': {'S2V': -0.244, 'S2H': -0.125, 'S1V': 0.0, 'S1H': 0.09}, 'MEG': {'S2V': 0.0, 'S2H': 0.284, 'S1V': -0.055, 'S1H': 0.156}}
+
+======= 20220322: 
+{'HEG': {'S2V': -0.244, 'S2H': -0.125, 'S1V': 0.0, 'S1H': 0.09}, 'MEG': {'S2V': 0.0976, 'S2H': -0.002, 'S1V': -0.0353, 'S1H': 0.1558}}
+
+======= 20220329: 
+{'HEG': {'S2V': -0.244, 'S2H': -0.125, 'S1V': 0.0, 'S1H': 0.09}, 'MEG': {'S2V': 0.05, 'S2H': -0.02, 'S1V': -0.05, 'S1H': 0.1}}

instruments/IEX_dictionaries/Dict_TempDiode.txt

+{"SI435":[0.00000,1.66751,1.64531,1.61210,1.57373,1.53247,1.49094,1.45196,1.41723,1.38705,1.36089,1.33785,1.31699,1.29756,1.27912,1.26154,1.24489,1.22917,1.21406,1.19855,1.18193,1.16246,1.13841,1.12425,1.11828,1.11480,1.11217,1.10996,1.10828,1.10643,1.10465,1.10295,1.10124,1.09952,1.09791,1.09629,1.09468,1.09306,1.09145,1.08983,1.08821,1.08659,1.08497,1.08334,1.08171,1.08008,1.07844,1.07681,1.07517,1.07353,1.07188,1.07023,1.06858,1.06693,1.06528,1.06362,1.06196,1.06029,1.05863,1.05696,1.05528,1.05360,1.05192,1.05024,1.04856,1.04687,1.04518,1.04349,1.04179,1.04010,1.03839,1.03669,1.03498,1.03327,1.03156,1.02985,1.02814,1.02642,1.02471,1.02299,1.02127,1.01957,1.01785,1.01612,1.01439,1.01267,1.01093,1.00918,1.00744,1.00569,1.00393,1.00218,1.00042,0.99865,0.99687,0.99510,0.99332,0.99153,0.98975,0.98795,0.98615,0.98434,0.98254,0.98073,0.97891,0.97710,0.97527,0.97344,0.97161,0.96978,0.96794,0.96609,0.96424,0.96239,0.96053,0.95867,0.95680,0.95492,0.95304,0.95116,0.94927,0.94738,0.94549,0.94358,0.94167,0.93976,0.93785,0.93592,0.93398,0.93203,0.93008,0.92812,0.92616,0.92418,0.92221,0.92022,0.91823,0.91623,0.91423,0.91222,0.91021,0.90819,0.90617,0.90414,0.90212,0.90008,0.89805,0.89601,0.89397,0.89293,0.88988,0.88783,0.88578,0.88372,0.88166,0.87959,0.87752,0.87545,0.87337,0.87129,0.86921,0.86712,0.86503,0.86294,0.86084,0.85874,0.85664,0.85453,0.85242,0.85030,0.84818,0.84606,0.84393,0.84180,0.83967,0.83754,0.83540,0.83325,0.83111,0.82896,0.82680,0.82465,0.82249,0.82032,0.81816,0.81599,0.81381,0.81163,0.80945,0.80727,0.80508,0.80290,0.80071,0.79851,0.79632,0.79412,0.79192,0.78972,0.78752,0.78532,0.78311,0.78090,0.77869,0.77648,0.77426,0.77205,0.76983,0.76761,0.76539,0.76317,0.76094,0.75871,0.75648,0.75425,0.75202,0.74978,0.74755,0.74532,0.74308,0.74085,0.73861,0.73638,0.73414,0.73191,0.72967,0.72743,0.72520,0.72296,0.72072,0.71848,0.71624,0.71400,0.71176,0.70951,0.70727,0.70503,0.70278,0.70054,0.69829,0.69604,0.69379,0.69154,0.68929,0.68704,0.68479,0.68253,0.68028,0.67802,0.67576,0.67351,0.67124,0.66898,0.66672,0.66445,0.66219,0.65992,0.65765,0.65538,0.65310,0.65083,0.64855,0.64628,0.64400,0.64172,0.63944,0.63716,0.63487,0.63259,0.63030,0.62802,0.62573,0.62344,0.62115,0.61885,0.61656,0.61427,0.61197,0.60968,0.60738,0.60509,0.60279,0.60050,0.59820,0.59590,0.59360,0.59131,0.58901,0.58671,0.58441,0.58211,0.57980,0.57750,0.57520,0.57289,0.57059,0.56828,0.56598,0.56367,0.56136,0.55905,0.55674,0.55443,0.55211,0.54980,0.54748,0.54516,0.54285,0.54053,0.53820,0.53588,0.53356,0.53123,0.52891,0.52658,0.52425,0.52192,0.51958,0.51725,0.51492,0.51258,0.51024,0.50790,0.50556,0.50322,0.50088,0.49854,0.49620,0.49385,0.49151,0.48916,0.48681,0.48446,0.48212,0.47977,0.47741,0.47506,0.47271,0.47036,0.46801,0.46565,0.46330,0.46095,0.45860,0.45624,0.45389,0.45154,0.44918,0.44683,0.44447,0.44212,0.43976,0.43741,0.43505,0.43270,0.43034,0.42799,0.42563,0.42327,0.42092,0.41856,0.41620,0.41384,0.41149,0.40913,0.40677,0.40442,0.40206,0.39970,0.39735,0.39499,0.39263,0.39027,0.38792,0.38556,0.38320,0.38085,0.37849,0.37613,0.37378,0.37142,0.36906,0.36670,0.36435,0.36199,0.35963,0.35728,0.35492,0.35256,0.35021,0.34785,0.34549,0.34313,0.34078,0.33842,0.33606,0.33371,0.33135,0.32899,0.32667,0.32428,0.32192],"SI410":[0.0000,1.7191,1.7086,1.6852,1.6530,1.6124,1.5659,1.5179,1.4723,1.4309,1.3956,1.3656,1.3385,1.3142,1.2918,1.2712,1.2517,1.2333,1.2151,1.1963,1.1759,1.1524,1.1293,1.1192,1.1146,1.1114,1.1090,1.1069,1.1049,1.1031,1.1014,1.0997,1.0980,1.0964,1.0949,1.0933,1.0917,1.0902,1.0886,1.0871,1.0855,1.0839,1.0824,1.0808,1.0792,1.0776,1.0760,1.0744,1.0728,1.0712,1.0696,1.0679,1.0663,1.0646,1.0630,1.0613,1.0597,1.0580,1.0563,1.0547,1.0530,1.0513,1.0497,1.0480,1.0463,1.0446,1.0429,1.0412,1.0395,1.0378,1.0361,1.0344,1.0327,1.0310,1.0293,1.0276,1.0259,1.0242,1.0224,1.0207,1.0190,1.0172,1.0155,1.0137,1.0120,1.0102,1.0085,1.0067,1.0049,1.0032,1.0014,0.9996,0.9978,0.9960,0.9942,0.9924,0.9905,0.9887,0.9869,0.9851,0.9832,0.9814,0.9795,0.9777,0.9758,0.9740,0.9721,0.9703,0.9684,0.9665,0.9646,0.9628,0.9609,0.9590,0.9571,0.9552,0.9533,0.9514,0.9495,0.9476,0.9457,0.9437,0.9418,0.9398,0.9379,0.9359,0.9340,0.9320,0.9300,0.9281,0.9261,0.9241,0.9222,0.9202,0.9182,0.9162,0.9142,0.9122,0.9102,0.9082,0.9062,0.9042,0.9022,0.9002,0.8982,0.8962,0.8942,0.8921,0.8901,0.8881,0.8860,0.8840,0.8820,0.8799,0.8779,0.8758,0.8738,0.8717,0.8696,0.8676,0.8655,0.8634,0.8613,0.8593,0.8572,0.8551,0.8530,0.8509,0.8488,0.8467,0.8446,0.8425,0.8404,0.8383,0.8362,0.8341,0.8320,0.8299,0.8278,0.8257,0.8235,0.8214,0.8193,0.8171,0.8150,0.8129,0.8107,0.8086,0.8064,0.8043,0.8021,0.8000,0.7979,0.7957,0.7935,0.7914,0.7892,0.7871,0.7849,0.7828,0.7806,0.7784,0.7763,0.7741,0.7719,0.7698,0.7676,0.7654,0.7632,0.7610,0.7589,0.7567,0.7545,0.7523,0.7501,0.7479,0.7457,0.7435,0.7413,0.7391,0.7369,0.7347,0.7325,0.7303,0.7281,0.7259,0.7237,0.7215,0.7193,0.7170,0.7148,0.7126,0.7104,0.7082,0.7060,0.7037,0.7015,0.6993,0.6971,0.6948,0.6926,0.6904,0.6881,0.6859,0.6837,0.6814,0.6792,0.6770,0.6747,0.6725,0.6702,0.6680,0.6657,0.6635,0.6612,0.6590,0.6567,0.6545,0.6522,0.6500,0.6477,0.6455,0.6432,0.6410,0.6387,0.6365,0.6342,0.6319,0.6297,0.6274,0.6251,0.6229,0.6206,0.6183,0.6161,0.6138,0.6115,0.6092,0.6070,0.6047,0.6024,0.6001,0.5979,0.5956,0.5933,0.5910,0.5887,0.5865,0.5842,0.5819,0.5796,0.5773,0.5750,0.5727,0.5705,0.5682,0.5659,0.5636,0.5613,0.5590,0.5567,0.5544,0.5521,0.5498,0.5475,0.5452,0.5429,0.5406,0.5383,0.5360,0.5337,0.5314,0.5291,0.5268,0.5245,0.5222,0.5199,0.5176,0.5153,0.5130,0.5107,0.5084,0.5061,0.5038,0.5015,0.4992,0.4968,0.4945,0.4922,0.4899,0.4876,0.4853,0.4830,0.4806,0.4783,0.4760,0.4737,0.4714,0.4690,0.4667,0.4644,0.4621,0.4597,0.4574,0.4551,0.4528,0.4504,0.4481,0.4458,0.4435,0.4411,0.4388,0.4365,0.4341,0.4318,0.4295,0.4271,0.4248,0.4225,0.4201,0.4178,0.4154,0.4131,0.4108,0.4084,0.4061,0.4037,0.4014,0.3991,0.3967,0.3944,0.4154,0.4131,0.4108,0.4084,0.4061,0.4037,0.4014,0.3991,0.3967,0.3709,0.3685,0.3662,0.3638,0.3615,0.3591,0.3567,0.3544,0.3520,0.3497,0.3473,0.3449,0.3426,0.3402,0.3379,0.3355,0.3331,0.3308,0.3284,0.3260,0.3237,0.3213,0.3189,0.3165,0.3142,0.3118,0.3094,0.3071,0.3047,0.3023,0.2999,0.2976,0.2952,0.2928,0.2904,0.2880,0.2857,0.2833,0.2809,0.2785,0.2761,0.2738,0.2714,0.2690,0.2666,0.2642,0.2618,0.2594,0.2570,0.2547,0.2523,0.2499,0.2475,0.2451,0.2427,0.2403,0.2379,0.2355,0.2331,0.2307,0.2283,0.2259,0.2235,0.2211,0.2187,0.2163,0.2139,0.2115,0.2091,0.2067,0.2043]}

instruments/IEX_dictionaries/DiodeResponsivityCurve

+Responsivity Energy
+0.109341 10.0752
+0.1075 10.168
+0.105579 10.2777
+0.103737 10.3979
+0.101656 10.5007
+0.099895 10.614
+0.0978137 10.719
+0.0959725 10.8153
+0.0943715 10.9124
+0.0926103 11.0203
+0.0910093 11.1193
+0.0896484 11.1792
+0.0880474 11.2695
+0.0864464 11.3809
+0.0847653 11.514
+0.0830042 11.6591
+0.0812431 11.8378
+0.0798822 11.9762
+0.0790816 12.2033
+0.0788414 12.457
+0.0790015 12.7502
+0.079882 12.992
+0.081403 13.2029
+0.083084 13.4172
+0.0842847 13.5257
+0.0866862 13.733
+0.0882072 13.844
+0.0898882 14.006
+0.0918894 14.1572
+0.0939708 14.3613
+0.0965323 14.5553
+0.0983735 14.6862
+0.100215 14.8049
+0.101976 14.9113
+0.103817 15.0453
+0.106058 15.2077
+0.107899 15.3169
+0.109981 15.4546
+0.112062 15.5796
+0.113983 15.7196
+0.115824 15.8325
+0.117585 15.9748
+0.119747 16.1328
+0.121588 16.2342
+0.123349 16.3654
+0.12535 16.5421
+0.127832 16.7207
+0.130233 16.9466
+0.132419 17.1223
+0.135033 17.342
+0.137321 17.5348
+0.139608 17.7901
+0.141896 18.0644
+0.144837 18.3742
+0.14756 18.6576
+0.149957 18.881
+0.152462 19.1722
+0.154532 19.4348
+0.157364 19.7345
+0.159325 20.0048
+0.161395 20.2615
+0.163356 20.5216
+0.165752 20.7673
+0.167822 21.016
+0.170545 21.3039
+0.172942 21.5957
+0.175774 21.9287
+0.178389 22.248
+0.180131 22.4379
+0.182419 22.7259
+0.184707 22.9981
+0.186667 23.2932
+0.188193 23.6323
+0.189718 24.0376
+0.191352 24.4498
+0.193203 24.975
+0.194837 25.4466
+0.196036 25.8609
+0.197016 26.3044
+0.198105 26.9839
+0.199086 27.5401
+0.19963 27.9886
+0.200284 28.5655
+0.201062 29.1788
+0.202043 29.8563
+0.202874 30.4717
+0.203931 31.2588
+0.204762 31.8625
+0.205385 32.3423
+0.205728 32.6468
+0.2059 33.0831
+0.205685 33.4729
+0.205685 33.9202
+0.205814 34.4273
+0.205942 34.8329
+0.206423 35.5101
+0.206929 36.335
+0.207582 37.2736
+0.208236 38.2689
+0.20878 39.2241
+0.209543 40.3401
+0.210196 41.4878
+0.211068 43.0324
+0.211721 44.1815
+0.212375 45.5159
+0.213029 46.8506
+0.213847 48.4147
+0.21499 51.095
+0.217278 56.3781
+0.220137 63.5332
+0.222901 71.5962
+0.225284 79.93
+0.227667 88.1947
+0.23024 98.0003
+0.232909 110.697
+0.235768 124.745
+0.238437 137.322
+0.240915 151.521
+0.24406 168.367
+0.24692 187.965
+0.25016 211.325
+0.252924 234.82
+0.256451 265.864
+0.259215 294.732
+0.26217 329.039
+0.263599 344.016
+0.26522 365.622
+0.266649 386.769
+0.267602 405.322
+0.268365 421.79
+0.269127 442.023
+0.270176 464.314
+0.271034 482.048
+0.271863 502.915
+0.272661 524.768
+0.273353 544.89
+0.273888 560.168
+0.274296 577.813
+0.274432 596.754
+0.274432 613.265
+0.274387 630.233
+0.274387 655.755
+0.274577 677.334
+0.274896 698.718
+0.275163 720.777
+0.275588 742.318
+0.275801 762.005
+0.275908 786.062
+0.275961 813.536
+0.275961 836.479
+0.275961 860.069
+0.275961 881.435
+0.275854 906.293
+0.275804 934.582
+0.275686 968.151
+0.27545 999.14
+0.275214 1031.12
+0.275095 1056.11
+0.274918 1083.07
+0.274682 1109.32
+0.274446 1141.94
+0.27421 1174.05
+0.273856 1210.1
+0.273561 1248.84
+0.273207 1287.19
+0.272819 1322.48
+0.272499 1357.16
+0.272086 1386.55
+0.271732 1413.01
+0.271437 1436.35
+0.271161 1464.13
+0.270905 1482.33
+0.270669 1501.13
+0.270374 1520.17
+0.269961 1543.33
+0.26943 1568.82
+0.268401 1596.21
+0.26707 1629.27
+0.265713 1664.55
+0.264474 1696.32
+0.26347 1724.34
+0.26228 1755.08
+0.261215 1786.94
+0.259937 1822.37
+0.259032 1849.4
+0.258155 1871.57
+0.257397 1894.63
+0.256576 1915.74
+0.255377 1947.28
+0.275058 1957.05
+0.275058 2160.44
+0.275058 2389.74
+0.275058 2507.09
+0.275058 2604.07
+0.275058 2715.63
+0.275058 2962.14
+0.275058 3276.53
+0.275058 3649.72
+0.275058 3981.03

instruments/IEX_dictionaries/DiodeResponsivityCurve_2

+Responsivity	Energy
+0.109341	10.0752
+0.1075	10.168
+0.105579	10.2777
+0.103737	10.3979
+0.101656	10.5007
+0.099895	10.614
+0.0978137	10.719
+0.0959725	10.8153
+0.0943715	10.9124
+0.0926103	11.0203
+0.0910093	11.1193
+0.0896484	11.1792
+0.0880474	11.2695
+0.0864464	11.3809
+0.0847653	11.514
+0.0830042	11.6591
+0.0812431	11.8378
+0.0798822	11.9762
+0.0790816	12.2033
+0.0788414	12.457
+0.0790015	12.7502
+0.079882	12.992
+0.081403	13.2029
+0.083084	13.4172
+0.0842847	13.5257
+0.0866862	13.733
+0.0882072	13.844
+0.0898882	14.006
+0.0918894	14.1572
+0.0939708	14.3613
+0.0965323	14.5553
+0.0983735	14.6862
+0.100215	14.8049
+0.101976	14.9113
+0.103817	15.0453
+0.106058	15.2077
+0.107899	15.3169
+0.109981	15.4546
+0.112062	15.5796
+0.113983	15.7196
+0.115824	15.8325
+0.117585	15.9748
+0.119747	16.1328
+0.121588	16.2342
+0.123349	16.3654
+0.12535	16.5421
+0.127832	16.7207
+0.130233	16.9466
+0.132419	17.1223
+0.135033	17.342
+0.137321	17.5348
+0.139608	17.7901
+0.141896	18.0644
+0.144837	18.3742
+0.14756	18.6576
+0.149957	18.881
+0.152462	19.1722
+0.154532	19.4348
+0.157364	19.7345
+0.159325	20.0048
+0.161395	20.2615
+0.163356	20.5216
+0.165752	20.7673
+0.167822	21.016
+0.170545	21.3039
+0.172942	21.5957
+0.175774	21.9287
+0.178389	22.248
+0.180131	22.4379
+0.182419	22.7259
+0.184707	22.9981
+0.186667	23.2932
+0.188193	23.6323
+0.189718	24.0376
+0.191352	24.4498
+0.193203	24.975
+0.194837	25.4466
+0.196036	25.8609
+0.197016	26.3044
+0.198105	26.9839
+0.199086	27.5401
+0.19963	27.9886
+0.200284	28.5655
+0.201062	29.1788
+0.202043	29.8563
+0.202874	30.4717
+0.203931	31.2588
+0.204762	31.8625
+0.205385	32.3423
+0.205728	32.6468
+0.2059	33.0831
+0.205685	33.4729
+0.205685	33.9202
+0.205814	34.4273
+0.205942	34.8329
+0.206423	35.5101
+0.206929	36.335
+0.207582	37.2736
+0.208236	38.2689
+0.20878	39.2241
+0.209543	40.3401
+0.210196	41.4878
+0.211068	43.0324
+0.211721	44.1815
+0.212375	45.5159
+0.213029	46.8506
+0.213847	48.4147
+0.21499	51.095
+0.217278	56.3781
+0.220137	63.5332
+0.222901	71.5962
+0.225284	79.93
+0.227667	88.1947
+0.23024	98.0003
+0.232909	110.697
+0.235768	124.745
+0.238437	137.322
+0.240915	151.521
+0.24406	168.367
+0.24692	187.965
+0.25016	211.325
+0.252924	234.82
+0.256451	265.864
+0.259215	294.732
+0.26217	329.039
+0.263599	344.016
+0.26522	365.622
+0.266649	386.769
+0.267602	405.322
+0.268365	421.79
+0.269127	442.023
+0.270176	464.314
+0.271034	482.048
+0.271863	502.915
+0.272661	524.768
+0.273353	544.89
+0.273888	560.168
+0.274296	577.813
+0.274432	596.754
+0.274432	613.265
+0.274387	630.233
+0.274387	655.755
+0.274577	677.334
+0.274896	698.718
+0.275163	720.777
+0.275588	742.318
+0.275801	762.005
+0.275908	786.062
+0.275961	813.536
+0.275961	836.479
+0.275961	860.069
+0.275961	881.435
+0.275854	906.293
+0.275804	934.582
+0.275686	968.151
+0.27545	999.14
+0.275214	1031.12
+0.275095	1056.11
+0.274918	1083.07
+0.274682	1109.32
+0.274446	1141.94
+0.27421	1174.05
+0.273856	1210.1
+0.273561	1248.84
+0.273207	1287.19
+0.272819	1322.48
+0.272499	1357.16
+0.272086	1386.55
+0.271732	1413.01
+0.271437	1436.35
+0.271161	1464.13
+0.270905	1482.33
+0.270669	1501.13
+0.270374	1520.17
+0.269961	1543.33
+0.26943	1568.82
+0.268401	1596.21
+0.26707	1629.27
+0.265713	1664.55
+0.264474	1696.32
+0.26347	1724.34
+0.26228	1755.08
+0.261215	1786.94
+0.259937	1822.37
+0.259032	1849.4
+0.258155	1871.57
+0.257397	1894.63
+0.256576	1915.74
+0.255377	1947.28
+0.275058	1957.05
+0.275058	2160.44
+0.275058	2389.74
+0.275058	2507.09
+0.275058	2604.07
+0.275058	2715.63
+0.275058	2962.14
+0.275058	3276.53
+0.275058	3649.72
+0.275058	3981.03

instruments/IEX_dictionaries/M0M1_positions.py

+M0M1_Pos['default']={0:'-0.400/-22.5/0/0/0.000/0.0000',1:'0.400/-21.5/0/0/8.400/2.800'}
+M0M1_Pos['2019_1']= {0:'0.3010/-22.5/0/0/-0.151/0.0393',1:'1.449/-21.5/0/0/8.339/2.700'}
+M0M1_Pos['2019_2']= {0:'-0.400/-22.5/0/0/ 0.000/0.000',1:'0.400 /-21.5/0/0/8.436/3.000'}
\ No newline at end of file

instruments/IEX_dictionaries/calibrationToFix.txt

+
+======= 20211108: 
+{'HEG': {0: [[600, [-6635.813431810828, 53.075898114105726, -0.15253171108859886, 0.00019792480095656869, -9.598209896301075e-08]], [3000, [46.600135556920925, 0.8412063733137993, 0.00023492720237763438, -1.2593652551917328e-07, 2.5435326557098538e-11]]], 3: [[400, [1051.452570904546, -13.099314837713008, 0.06939879665006896, -0.00014784669111041307, 1.155031075284652e-07]], [600, [2782.786054275039, -21.145755715684032, 0.06554921272398154, -8.534237586539917e-05, 4.134061079155964e-08]], [3000, [87.23060103156234, 0.7171860430279529, 0.0003793918782124725, -1.9605428816497513e-07, 3.6770897051532086e-11]]], 2: [[590, [-12422.076427583952, 100.54771064816772, -0.29845648627641835, 0.00039688772576782855, -1.97439876026206e-07]], [3000, [-33.02431482418943, 1.1302724323594304, -0.00013479679549416608, 7.22272393581498e-08, -1.373032035520199e-11]]], 1: [[600.0, [-4765.159827798981, 39.08508762376964, -0.11372474852884872, 0.00015063687981483728, -7.456557379431018e-08]], [3000, [84.53488099145231, 0.7066601865786498, 0.00043283422433330405, -2.2608518373215027e-07, 4.202604094132355e-11]]]},
+{'MEG': {
+0: [[600.0, [-5940.176848239725, 48.30636368582396, -0.14075467796987542, 0.00018579017325356665, -9.17653113254405e-08]], [2200.0, [20.344164869431648, 0.942351240250947, 0.00010121556917763097, -4.949527030585497e-08, 9.675448511744878e-12]], [3000.0, [365297.2627682823, -552.812440238765, 0.3141479258407871, -7.90104068928604e-05, 7.434922126836773e-09]]], 
+
+3: [[600.0, [-634.0261758446517, 7.262520018224911, -0.022419323576083252, 3.49681129089603e-05, -1.999932008653819e-08]], [2200.0, [58.93716149881088, 0.8198484295498225, 0.0002490011121555607, -1.2688431897796926e-07, 2.37752539333468e-11]], [3000.0, [112089.95145715227, -159.53263761672966, 0.08554887203588527, -2.0076278347329386e-05, 1.749084702165079e-09]]], 
+
+2: [[600.0, [-604.5006436697821, 3.0530278050672597, 0.0006285185258970196, -7.876285868395855e-06, 6.618688516815772e-09]], [2200.0, [34.67137980324179, 0.8848938354598829, 0.00017775770898386393, -9.469862981232054e-08, 1.8324158207531115e-11]], [3000.0, [409195.45122128195, -612.7658620540126, 0.3445915452908224, -8.582156181499232e-05, 8.000918394998031e-09]]], 
+
+1: [[600.0, [-3677.005333394253, 29.772255818860682, -0.0843978889369737, 0.00011023738262993015, -5.399238835077749e-08]], [2200.0, [46.243321909129264, 0.8454613578261618, 0.00024604317670495594, -1.1956387793694198e-07, 2.0350719268411266e-11]], [3000.0, [-99353.40775512981, 160.7507495025834, -0.0962536588691068, 2.5766240641565706e-05, -2.583875704169572e-09]]]}
\ No newline at end of file

instruments/IEX_dictionaries/mono_parameters.txt

+   MonoParm["20160822"]={
+        'Au':[0.3237,-12.04],'Si':[0.3237,-2.04],'Carbon':[0.3237,7.96],
+        'MEG_Imp':[0.7290,4.872,1200.0,-7.34e-05],'HEG_JY':[0.7302,69.0,2400.0,-5.58e-05],'MEG_JY':[0.8,131.624,1200.0,0.0]}
+    MonoParm["20190719_4"]={
+        'Au':[0.3266,-12.04],'Si':[0.3237,-2.04],'Carbon':[0.3237,7.96],
+        'MEG_Imp':[0.7306,4.872,1200.0,-7.15e-05],'HEG_JY':[0.7278,69.0,2400.0,-5.535e-05],'MEG_JY':[0.8,131.624,1200.0,0.0]}
+    MonoParm['20191018']= {
+        'Au':[0.3203,-12.04],'Si':[0.3237,-2.04],'Carbon':[0.3237,7.96],
+        'MEG_Imp':[0.7306,4.872,1200.0,-7.15e-05],'HEG_JY':[0.7123,69.0,2400.0,-5.535e-05],'MEG_JY':[0.7306,131.624,1200.0,0.0],}
+   
\ No newline at end of file

instruments/IEX_endstations.py

+from math import floor
+import time
+
+from epics import caput
+"""
+Endstation class is used to contain information about ioc, branch, mode
+
+it makes the default stuff work
+
+this will prompt for the endstation and will set the default parameters, look at the 
+init kwargs to see which defaults you can change.
+"""
+
+class Endstation:
+    """
+    used for defining which endstation in which you are running
+    for short names and ioc info
+
+    BL = endstation('ARPES')
+        BL.ioc = "ARPES" / "Kappa" / "Test"
+        BL.mode = "User" / "Staff" / "No Xray"
+        BL.branch = 'c' / 'd'
+        BL.mda => mda scanRecord
+
+    """
+
+    def __init__(self,endstation_name,scan_ioc,xrays,BL_mode,mda_scanRecord):
+        """
+        intializes the several beamline variables 
+
+        endstation_name = 'ARPES' / 'Kappa'
+        
+        BL = Endstation()
+            BL.endstation => endstation_name
+            BL.xrays => True/False 
+            BL.mode => previously: BL_Mode_Set
+            BL.folder => 'b','c','d'
+            BL.prefix => 'ARPES_','Kappa_'
+            BL.ioc => previously: BL_ioc()
+            BL.mda_filepath
+            
+        """
+        global BL
+        endstations_list = ['ARPES','Kappa']   
+        BL_mode_list = ['user','staff']
+
+        if self.endstation_name in endstations_list:
+            self.endstation=endstation_name
+        else:
+            print('Not a valid Endstation choice')
+            print('Endstations: '+endstations_list)
+            return
+
+        if BL_mode in BL_mode_list:
+            self.mode = BL_mode
+            if BL_mode.lower == 'user':
+                if endstation_name == 'ARPES':
+                    self.folder = 'c'
+                    self.prefix = 'ARPES_'
+                elif endstation_name == 'Kappa':
+                    self.folder = 'd'
+                    self.prefix = 'Kappa_'
+                #elif endstation_name == 'Octupole':
+                    #self.folder = 'e'
+            elif BL_mode.lower == 'staff':
+                self.folder = 'b' #overwrite folder for staff mode
+        else:
+            print('Not a valid BL_mode choice')
+            print('BL_modes: '+BL_mode_list) 
+            return  
+
+        self.xrays = xrays
+        self.ioc = scan_ioc
+        self.mda = mda_scanRecord
+
+        global mda
+        mda = mda_scanRecord
+
+        global Motors
+
+
+
+    def set_logfile_path():
+        """
+        sets the default logfile path
+        """
+
+##############################################################################################################
+##############################               BL commands            ##############################
+##############################################################################################################
+def BL_ioc():
+    """
+    returns the branch from the Endstation instance
+    """
+    return BL.ioc
+
+def BL_mode():
+    """
+    returns the beamline mode, User / Staff / No_Xray
+    """
+    return BL.ioc
+
+def BL_mda_prefix():
+    """
+    returns the mda file prefix
+    """
+    return BL.prefix
+
+def BL_mda_filepath():
+    """
+    returns the mda file prefix
+    """
+    return BL.filepath
+
+def scalar_cts(self,integration_time=0.1,verbose=True,**kwargs):
+    """
+    Sets the integration time for the scalers
+    kwargs:
+        mcp = True/False to sum mpa
+    Previously: cts, Kappa counts 
+    """
+
+    if BL.endstation == 'ARPES':
+        pass
+    elif BL.endstation == 'Kappa':
+        Kappa_scalar_pv = "29idMZ0:scaler1.TP"
+        mpa_Proc1_pv = "29iddMPA:Proc1:"
+
+        caput(Kappa_scalar_pv,integration_time)
+
+        if kwargs["mpa"]:
+            caput(mpa_Proc1_pv+'NumFilter',floor(time))
+
+        if verbose: 
+            print("Integration time set to:", str(time))
\ No newline at end of file

instruments/Kappa_Euler.py

+from epics import caget, caput
+
+from userCalcs import userCalcOut_clear
+
+#### Obsolete?
+
+def Bragg_Index_CalcOut(d,eV,th):
+    n=8
+    userCalcOut_pv = userCalcOut_clear(mda,n)
+    
+    h=4.135667516e-15
+    c=299792458
+    f=h*c*1e9*10
+
+    caput(userCalcOut_pv+".DESC","Bragg_Index")
+    caput(userCalcOut_pv+".A",d)
+    caput(userCalcOut_pv+".B",th)
+    caput(userCalcOut_pv+".C",eV)
+    caput(userCalcOut_pv+".D",pi)
+    caput(userCalcOut_pv+".E",f)
+    caput(userCalcOut_pv+".CALC$","SIN(B*D/"+str(180)+")*"+str(2)+"*A*C/E")
+
+
+
+def KtoE_th_CalcOut():
+    n=1
+    userCalcOut_pv = userCalcOut_clear(mda,n)
+
+    caput(userCalcOut_pv+".DESC","KtoE_th")
+    caput(userCalcOut_pv+".INPA","29idKappa:m8.RBV CP NMS")    #A=kth
+    caput(userCalcOut_pv+".INPB","29idKappa:m7.RBV CP NMS")    #B=kap
+    caput(userCalcOut_pv+".INPC","29idKappa:m1.RBV CP NMS")    #C=kphi
+    caput(userCalcOut_pv+".D",3.141592653589793)
+    caput(userCalcOut_pv+".E",180)
+    caput(userCalcOut_pv+".F",50)
+    caput(userCalcOut_pv+".H",57.045)
+    caput(userCalcOut_pv+".CALC$","((A-(G-H))*D/E-ATAN(TAN(B*D/E/2.0)*COS(F*D/E)))*E/D")
+    caput(userCalcOut_pv+".DOPT",0)           # Use CAL
+
+def KtoE_chi_CalcOut():
+    n=2
+    userCalcOut_pv = userCalcOut_clear(mda,n)
+
+    caput(userCalcOut_pv+".DESC","KtoE_chi")
+    caput(userCalcOut_pv+".INPA","29idKappa:m8.RBV CP NMS")    #A=kth
+    caput(userCalcOut_pv+".INPB","29idKappa:m7.RBV CP NMS")    #B=kap
+    caput(userCalcOut_pv+".INPC","29idKappa:m1.RBV CP NMS")    #C=kphi
+    caput(userCalcOut_pv+".D",3.141592653589793)
+    caput(userCalcOut_pv+".E",180)
+    caput(userCalcOut_pv+".F",50)
+    caput(userCalcOut_pv+".CALC$","(2*ASIN(SIN(B*D/E/2.0) * SIN(F*D/E)))*E/D")
+    caput(userCalcOut_pv+".DOPT",0)           # Use CAL
+
+def KtoE_phi_CalcOut():
+    n=3
+    userCalcOut_pv = userCalcOut_clear(mda,n)
+    
+    caput(userCalcOut_pv+".DESC","KtoE_phi")
+    caput(userCalcOut_pv+".INPA","29idKappa:m8.RBV CP NMS")    #A=kth
+    caput(userCalcOut_pv+".INPB","29idKappa:m7.RBV CP NMS")    #B=kap
+    caput(userCalcOut_pv+".INPC","29idKappa:m1.RBV CP NMS")    #C=kphi
+    caput(userCalcOut_pv+".D",3.141592653589793)
+    caput(userCalcOut_pv+".E",180)
+    caput(userCalcOut_pv+".F",50)
+    caput(userCalcOut_pv+".CALC$","(C*D/E-ATAN(TAN(B*D/E/2.0)*COS(F*D/E)))*E/D")
+    caput(userCalcOut_pv+".DOPT",0)           # Use CAL
+
+
+
+
+def EtoK_kth_CalcOut():
+    n=4
+    userCalcOut_pv = userCalcOut_clear(mda,n)
+
+    caput(userCalcOut_pv+".DESC","EtoK_kth")
+    caput(userCalcOut_pv+".INPA","29idKappa:userCalcOut1.VAL CP NMS")    #A=th
+    caput(userCalcOut_pv+".INPB","29idKappa:userCalcOut2.VAL CP NMS")    #B=chi
+    caput(userCalcOut_pv+".INPC","29idKappa:userCalcOut3.VAL CP NMS")    #C=phi
+    caput(userCalcOut_pv+".D",3.141592653589793)
+    caput(userCalcOut_pv+".E",180)
+    caput(userCalcOut_pv+".F",50)
+    caput(userCalcOut_pv+".CALC$","(A*D/E-ASIN(-TAN(B*D/E/2.0)/TAN(F*D/E)))*E/D")
+    caput(userCalcOut_pv+".DOPT",0)           # Use CAL
+
+def EtoK_kap_CalcOut():
+    n=5
+    scanIOC='Kappa'
+    ClearCalcOut(scanIOC,n)
+    pvcal="29id"+scanIOC+":userCalcOut"+str(n)
+    caput(pvcal+".DESC","EtoK_kap")
+    caput(pvcal+".INPA","29idKappa:userCalcOut1.VAL CP NMS")    #A=th
+    caput(pvcal+".INPB","29idKappa:userCalcOut2.VAL CP NMS")    #B=chi
+    caput(pvcal+".INPC","29idKappa:userCalcOut3.VAL CP NMS")    #C=phi
+    caput(pvcal+".D",3.141592653589793)
+    caput(pvcal+".E",180)
+    caput(pvcal+".F",50)
+    caput(pvcal+".CALC$","((2*ASIN(SIN(B*D/E/2.0) / SIN(F*D/E))))*E/D")
+    caput(pvcal+".DOPT",0)           # Use CAL
+
+def EtoK_kphi_CalcOut():
+    n=6
+    userCalcOut_pv = userCalcOut_clear(mda,n)
+    
+    caput(userCalcOut_pv+".DESC","EtoK_kphi")
+    caput(userCalcOut_pv+".INPA","29idKappa:userCalcOut1.VAL CP NMS")    #A=th
+    caput(userCalcOut_pv+".INPB","29idKappa:userCalcOut2.VAL CP NMS")    #B=chi
+    caput(userCalcOut_pv+".INPC","29idKappa:userCalcOut3.VAL CP NMS")    #C=phi
+    caput(userCalcOut_pv+".D",3.141592653589793)
+    caput(userCalcOut_pv+".E",180)
+    caput(userCalcOut_pv+".F",50)
+    caput(userCalcOut_pv+".CALC$","(C*D/E-ASIN(-TAN(B*D/E/2.0)/TAN(F*D/E)))*E/D")
+    caput(userCalcOut_pv+".DOPT",0)           # Use CAL
+
+
+
+
+
+

instruments/Kappa_det.py

+from time import sleep
+from epics import caget, caput
+ 
+from .IEX_endstations import Motors
+
+##############################################################################################################
+################################            Kappa detector class             ##############################
+##############################################################################################################
+class Kappa_Detector:
+    """
+    class for Kappa detector
+    """
+ 
+    def __init__(self):
+        self.get()
+
+    def get(self):
+        """
+        gets the current det_name and position
+
+        det_name, tth_val
+        """
+        #det_name = caget('29idKappa:userStringSeq6.STR1')
+        det_name = caget(det_set_pv)
+        self.name = det_name
+        tth_val = Motors.get('tth')
+        return self.name, tth_val
+
+    def set(self,det_name,move=True):
+        """
+        sets the tth detector 
+        and moves the motor to tth value for the new detecotr
+        
+        det_names
+            d3: large diode
+            d4: small diode
+            mcp: for the mpa
+            apd: avalanche photodiode (not currently installed) 
+            yag: yag 'fluorescence screen'
+        """
+        #get current value for tth
+        tth_val = Motors.get('tth')
+
+        #change det 
+        if det_name in det_list:
+            caput(det_set_pv,det_name)
+            if move:
+                Motors.move('tth',tth_val,wait=True,verbose=False)
+
+    def tth0_set(move):
+        """
+        resetting the tth motor zero to correspond to direct beam
+
+        only works with d4
+        """
+        current_det=caget(det_set_pv,as_string=True)
+        tth_pv = Motors._motor_dictionary['th'][3]
+
+        if current_det != 'd4':
+            print('tth0 can only be redefined with d4')
+        else:
+            foo=input('Are you sure you want to reset tth0 (y or n)? >')
+            if foo == 'Y' or foo == 'y' or foo == 'yes'or foo == 'YES':
+                caput(tth_pv+'.SET','Set')
+                sleep(0.5)
+                caput(tth_pv+'.VAL',0)
+                caput(tth_pv+'.DVAL',0)
+                sleep(0.5)
+                caput(tth_pv+'.SET','Use')
+                print("tth position reset to 0")
+            else:
+                print("That's ok, everybody can get cold feet in tough situation...")
+

instruments/Kappa_optimization.py

+import matplotlib.pyplot as plt
+
+from epics import caget,caput
+from .Kappa import *
+from .Motors import *
+
+from IEX_plotting_and_analysis.IEX_plotting_and_analysis.ScanFunctions_plot import fit_mda, mda_1D
+
+##### Kappa Optimization scans #####
+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 to d4')
+        tthdet.set('d4',move=False)
+        mvtth(0)
+    if moveZ is not None:
+        mvz(z-1000)
+    for i in range(iteration):
+        scantth(start,stop,step,'relative')
+        scannum = mda.lastFileNum()
+        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,move=False)
+    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=mda.fileNum()
+        z0=fit_mda(scannum,det,z_w,'erf',graph=graph)
+        mvz(z0)
+    scanz(-700,700,50,'relative')
+    scannum = mda.lastFileNum()
+    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,move=False)
+    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(mda.lastFileNum(),det,0.1,'gauss',graph=graph)
+    mvkth(new_kth)
+    scankth(-0.2*i,0.2*i,0.02*i,'relative')
+    new_kth=fit_mda(mda.lastFileNum(),det,0.1,'gauss',graph=graph)
+    mvkth(new_kth)
+    scannum=mda.lastFileNum
+    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)