87 lines
2.5 KiB
Python
87 lines
2.5 KiB
Python
run("diffutils")
|
|
|
|
###################################################################################################\
|
|
#Setup
|
|
###################################################################################################
|
|
|
|
#alpha, delta, gamma, omegaV
|
|
setup_diff(fourcv, energy, ("mu", "delta", "gam", "eta"))
|
|
|
|
|
|
#setup_axis(alpha, alpha.getMinValue(), alpha.getMaxValue())
|
|
#setup_axis(delta, delta.getMinValue(), 90) #delta.getMaxValue())
|
|
#setup_axis(gamma, 0, gamma.getMaxValue())
|
|
#setup_axis(omegaV, omegaV.getMinValue(), omegaV.getMaxValue())
|
|
|
|
|
|
if energy.isSimulated():
|
|
wavelength.write(1.305)
|
|
|
|
|
|
"""
|
|
###################################################################################################\
|
|
#Orientation
|
|
###################################################################################################
|
|
help(ub.ub)
|
|
ub.listub()
|
|
|
|
#alpha delta gamma omegaV
|
|
|
|
# Create a new ub calculation and set lattice parameters
|
|
ub.newub('test')
|
|
|
|
#ub.setlat('cubic', 5.114, 5.8361, 11.058, 90, 90, 90)
|
|
#en = 8
|
|
##ub.c2th([0, 0, 4], en)
|
|
##ub.addref([0, 0, 4]) #From current position and ekergy
|
|
#ub.addref([0, 0, 4], [16.2785, 0.0, 32.5568, 0.0], en)
|
|
##ub.c2th([2, 0, 12], en)
|
|
#ub.addref([2, 0, 12], [71.8285, 37.3082, 138.7440, 0.0], en)
|
|
##ub.c2th([1, -4, 10], en)
|
|
#ub.addref([1, -4, 10], [27.7185, 17.6409 , 128.4220, 0.0], en)
|
|
|
|
#ub.setlat('cubic', 1.0, 1.0, 1.0, 90, 90, 90)
|
|
ub.setlat('cubic', 1.305, 1.305, 1.305, 90, 90, 90)
|
|
|
|
#en = 12.4
|
|
en = 9.5
|
|
#ub.c2th([0, 0, 4], en)
|
|
#ub.addref([0, 0, 4]) #From current position and ekergy
|
|
ub.addref([0, 0, 1], [30.0, 0.0, 60.0, 0.0], en)
|
|
#ub.c2th([2, 0, 12], en)
|
|
ub.addref([1, 0, 1], [20.0, 45.5564,90.000, 44.4437], en)
|
|
#ub.c2th([1, -4, 10], en)
|
|
ub.addref([0, 1, 1], [20.0, 45.5564,90.000, 134.4437], en)
|
|
|
|
ub.ub()
|
|
|
|
#ub.setub([[1.22862,0.00000,0.00000], [-0.00000,1.07663,0.00000], [-0.00000,-0.00000,0.56820]])
|
|
|
|
|
|
# check the state
|
|
|
|
ub.checkub()
|
|
|
|
|
|
|
|
###################################################################################################\
|
|
#Constraints
|
|
###################################################################################################
|
|
help(hkl.con)
|
|
#hkl.con('a_eq_b')
|
|
#hkl.con('eta:0')
|
|
|
|
|
|
hkl.con( 'eta', 0) #OmegaV constant
|
|
#hkl.con( 'mu', 20) #Alpha constant
|
|
|
|
|
|
###################################################################################################\
|
|
#Motion
|
|
###################################################################################################
|
|
|
|
#print angles_to_hkl((16.278, 0.0000, 32.5568, 0.0))
|
|
#print angles_to_hkl((44.3400, 0.0000, 123.7322 , 0.0))
|
|
#print hkl_to_angles(2, -2, 10)
|
|
"""
|