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

Script execution

Browse Source
This commit is contained in:
panepucci
2016-02-17 11:59:18 +01:00
parent 141a1eb7c1
commit a73d7a4fe7
1 changed files with 17 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

24
script/local.py
View File

@@ -1,28 +1,38 @@
###################################################################################################
# Deployment specific global definitions - executed after startup.py
###################################################################################################
BEAMLINE_XNAME = "X06DA"
dtor = math.pi /180.0
def A2E(t,h=1.0,k=1.0,l=1.0, deg =True, ln = False):
d0 =2 * 5.43102 * (1.0-lncorr) / math.sqrt(h**2+k**2+l**2)
tt = (t * dtor) if (deg or (t>1.0)) else t
lncorr= 2.e-4 if ln else 0.0
d0=2 * 5.43102 * (1.0-lncorr) / math.sqrt(h**2+k**2+l**2)
tt= (t * dtor) if (deg or (t>1.0)) else t
return 12.39842 / (d0*math.sin(tt))
def angle(e,h=1.0,k=1.0,l=1.0, deg =True, ln = False,bent = False):
lncorr = 2.e-4 if ln else 0.0
d0 =2 * 5.43102 * (1.0-lncorr) / math.sqrt(h**2+k**2+l**2)
a = math.asin(12.39842/d0/e)
if bent:
rho = 2*19.65*8.35/28*math.sin(a)
dt = 0.2e-3/rho*0.279
d0 = 2*5.43102d*(1+dt)/math.sqrt(h^2+k^2+l^2)
a = math.asin(12.39842/d0/E)
a = a/dtor if deg else t
def get_energy(debug_msg = True):
t2 = caget("X06DA-OP-MO1:ROX2.RBV",'d')
e = A2E(abs(t2))
if debug_msg:
print 'Energy [keV]:'+ str(e) + ' Wavelength [A]:' + str(12.39842/e)
return e
run("fit")
run("sete")
run("seteq")
run("rock")
#run("sete")
#run("seteq")