mask1 = flex.bool(flex.grid(int(ysize / div), int(xsize / div)))
mask2 = flex.bool(flex.grid(int(ysize / div), int(xsize / div)))
n = 5
dm = n * 0.157 * pi / 180.0
# dm = n * 1.0 * pi / 180.0
# dm /= pi / 2
print(dm)
xyz = detector.get_pixel_lab_coord((0, 1200))
s1 = matrix.col(xyz).normalize() * s0.length()
phi = 0.0
cs = CoordinateSystem(m2, s0, s1, phi)
print("Limits: ", cs.limits())
# from_angle = FromRotationAngleAccurate(cs)
# phi_list = []
# e3_list = []
# for t in range(-100, 100):
# p = phi - dm * t / 10
# print p, from_angle(p)
# phi_list.append(p)
# e3_list.append(from_angle(p))
# from matplotlib import pylab
# pylab.plot(phi_list, e3_list)
# pylab.axvline(x=phi)
# pylab.show()
mask1 = flex.bool(flex.grid(int(ysize / div), int(xsize / div))) mask2 = flex.bool(flex.grid(int(ysize / div), int(xsize / div))) n = 5 dm = n * 0.157 * pi / 180.0 #dm = n * 1.0 * pi / 180.0 #dm /= pi / 2 print dm xyz = detector.get_pixel_lab_coord((0, 1200)) s1 = matrix.col(xyz).normalize() * s0.length() phi = 0.0 cs = CoordinateSystem(m2, s0, s1, phi) print "Limits: ", cs.limits() #from_angle = FromRotationAngleAccurate(cs) #phi_list = [] #e3_list = [] #for t in range(-100, 100): # p = phi - dm * t / 10 # print p, from_angle(p) # phi_list.append(p) # e3_list.append(from_angle(p)) #from matplotlib import pylab #pylab.plot(phi_list, e3_list) #pylab.axvline(x=phi) #pylab.show()
