diff --git a/script/daq/daq1-setup.py b/script/daq/daq1-setup.py index 6692d12..3cab94c 100644 --- a/script/daq/daq1-setup.py +++ b/script/daq/daq1-setup.py @@ -3,6 +3,12 @@ ################################################################################################### set_geometry("fourcv", True) +###################################################################################################\ +#Constraints +################################################################################################### +hkl.con('a_eq_b') + + ###################################################################################################\ #Orientation ################################################################################################### @@ -42,15 +48,13 @@ ub.newub('test') ub.setub([[ 0.81374, 0.02138, 0.00228], [-0.02199, 0.81543, 0.00123], [ 0.00784, -0.00100, 0.81347]]) ub.checkub() -###################################################################################################\ -#Constraints -################################################################################################### -hkl.con('a_eq_b') print hkl_to_angles(0, 1, 2)[0] print hkl_to_angles(2, 2, 2)[0] print hkl_to_angles(2, 0, 2)[0] -for r in get_ub_matrix(): print r +orig_ub = get_ub_matrix() + + @@ -67,12 +71,24 @@ en = 9.5 ub.addref([0, 1, 2], [ 9.7212 , 9.7488 , 19.5860, 93.1295], en) ub.addref([2, 2, 2], [ 9.8097 , 27.7194 , 20.9065 , 57.5840], en) ub.addref([2, 0, 2], [9.8217 , 19.5132 , 20.2488 , 8.4510], en) -#ub.ub() -#ub.checkub() +ub.ub() +ub.checkub() -for r in get_ub_matrix(): print r +calc_ub = get_ub_matrix() print hkl_to_angles(0, 1, 2)[0] print hkl_to_angles(2, 2, 2)[0] print hkl_to_angles(2, 0, 2)[0] + + +for r in orig_ub: print r +for r in calc_ub: print r +print "---" + +print e +for i in range(3): + for j in range(3): + e[i][j] = abs(orig_ub[i][j] - calc_ub[i][j]) + +for r in e: print r