diff --git a/script/local.py b/script/local.py index 7f06f16..6dec498 100644 --- a/script/local.py +++ b/script/local.py @@ -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): 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 + 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("rock") -run("sete") -run("seteq") +#run("sete") +#run("seteq")