Files
x04sa-es3/script/test/TestNewDiffcalc.py
gac-x04sa 4ca06cc16f Startup
2019-02-11 08:23:01 +01:00

84 lines
2.5 KiB
Python

###################################################################################################\
#Geometry
###################################################################################################
set_geometry("fourcv", True)
###################################################################################################\
#Constraints
###################################################################################################
hkl.con('a_eq_b')
#hkl.con( 'eta', 0)
###################################################################################################\
#Orientation
###################################################################################################
setup_axis(alpha, alpha.getMinValue(), alpha.getMaxValue())
setup_axis(delta, delta.getMinValue(), 90) #delta.getMaxValue())
setup_axis(gamma, gamma.getMinValue(), gamma.getMaxValue())
setup_axis(omegaV, omegaV.getMinValue(), omegaV.getMaxValue())
###################################################################################################\
#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', 11.058, 5.8361, 5.114, 90, 90, 90)
en = 1.5498
#ub.c2th([0, 0, 4], en)
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.6408 , 128.4220, 0.0], en)
#ub.ub()
#ub.setub([[1.22862 ,0.00000,0.00000], [-0.00000,1.07663,0.00000], [-0.00000,-0.00000,0.56820]])
"""
ub.setlat('cubic', 5.114, 5.8361, 11.058, 90, 90, 90)
en = 1.5498
#en=8
#ub.c2th([0, 0, 4], en)
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.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()
###################################################################################################\
#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))[0]
print hkl_to_angles(2, -2, 6.59)[0]