Page History
Content
Table of Contents |
---|
Problem
Meng sees a problem with interpretation of data. She observes in data large sample of events which have 2-peaks in the qH around two most probable positions 1.1 and 1.7 1/Å. In contrast, our model of lattice does not show these peaks simultaneously.
...
A few crosschecks are presented below in order to address this issue.
Lattice
Triclinic crystal cell primitive vectors
We use a few methods from pyimgalgos.FiberIndexing
to generate lattice:
...
Code Block |
---|
# 3-d space primitive vectors: # a1 = ( 18.550000, 0.000000, 0.000000) # a2 = ( -5.435624, 26.645524, 0.000000) # a3 = ( 0.000000, 0.000000, 4.860100) |
Reciprocal space primitive vectors
b1, b2, b3 = reciprocal_from_bravias(a1, a2, a3)
- converts 3-d primitive vectors to reciprocal space primitive vectors.
...
Code Block |
---|
# reciprocal space primitive vectors: # b1 = ( 0.053908, 0.010997, -0.000000) # b2 = ( 0.000000, 0.037530, 0.000000) # b3 = ( 0.000000, -0.000000, 0.205757) |
Lattice generator and test
Then, call to method
x, y, z, r, h, k, l = lattice(b1, b2, b3, hmax, kmax, lmax, cdtype)
...
Code Block | ||
---|---|---|
| ||
( h, k) R(h,k)[1/A] ( 0, 0) 0.0000 ( 0, 1) 0.0375 ( 1, 0) 0.0550 (-1, 1) 0.0601 ( 1, 1) 0.0725 ( 0, 2) 0.0751 (-1, 2) 0.0837 ( 1, 2) 0.1015 (-2, 1) 0.1089 ( 2, 0) 0.1100 ( 0, 3) 0.1126 (-1, 3) 0.1150 (-2, 2) 0.1202 ( 2, 1) 0.1232 ( 1, 3) 0.1348 (-2, 3) 0.1408 ( 2, 2) 0.1451 (-1, 4) 0.1492 ( 0, 4) 0.1501 (-3, 1) 0.1618 ( 3, 0) 0.1651 (-3, 2) 0.1671 (-2, 4) 0.1675 ( 1, 4) 0.1699 ( 2, 3) 0.1724 ( 3, 1) 0.1764 (-3, 3) 0.1803 (-1, 5) 0.1847 ( 0, 5) 0.1876 ( 3, 2) 0.1945 (-2, 5) 0.1977 (-3, 4) 0.1997 ( 2, 4) 0.2031 ( 1, 5) 0.2058 (-4, 1) 0.2157 ( 3, 3) 0.2176 (-4, 2) 0.2179 ( 4, 0) 0.2201 (-1, 6) 0.2209 (-3, 5) 0.2238 ( 0, 6) 0.2252 (-4, 3) 0.2263 (-2, 6) 0.2300 ( 4, 1) 0.2305 ( 2, 5) 0.2357 (-4, 4) 0.2403 ( 1, 6) 0.2423 ( 3, 4) 0.2443 ( 4, 2) 0.2463 (-3, 6) 0.2512 (-4, 5) 0.2591 ( 4, 3) 0.2665 ( 2, 6) 0.2697 ( 3, 5) 0.2736 (-4, 6) 0.2817 ( 4, 4) 0.2901 ( 3, 6) 0.3046 ( 4, 5) 0.3165 ( 4, 6) 0.3449 |
Rotation of the crystal lattice
Both in test_lattice
and make_lookup_table
methods crystal lattice in reciprocal space is rotated around Z and Y axes by omega_deg and beta_deg, respectively using calls:
...
Movies with frame delay 50, 100, 200, 400ms:
Look-up table
Look-up table is generated in method make_lookup_table
_v2 called as
...
Code Block |
---|
# file name: ./v02-lut-cxif5315-r0169-2016-02-03T15:10:48.txt # photon energy = 6003.1936 eV # wavelength = 2.0653 A # wave number/Evald radius k = 1/lambda = 0.484187 1/A # sigma_ql = 0.001453 1/A (approximately = k * <pixel size>/ # sigma_qt = 0.000484 1/A (approximately = k * <pixel size>/<sample-to-detector distance> = k*100um/100mm) # 3*sigma_ql = 0.004358 1/A # 3*sigma_qt = 0.001453 1/A # Triclinic crystal cell parameters: # a = 18.55 A # b = 27.19 A # c = 4.86 A # alpha = 90.00 deg # beta = 90.00 deg # gamma = 78.47 deg # 3-d space primitive vectors: # a1 = ( 18.550000, 0.000000, 0.000000) # a2 = ( -5.435624, 26.645524, 0.000000) # a3 = ( 0.000000, 0.000000, 4.860100) # reciprocal space primitive vectors: # b1 = ( 0.053908, 0.010997, -0.000000) # b2 = ( 0.000000, 0.037530, 0.000000) # b3 = ( 0.000000, -0.000000, 0.205757) # _________________________________________________________________________________________ # beta 0.00 omega 0.00 degree # index beta omega h k l dr[1/A] R(h,k,l) qv[1/A] qh[1/A] P(omega) 1 0.00 0.00 0 -1 0 0.001452 0.037530 0.000000 -0.037446 0.606641 1 0.00 0.00 0 1 0 0.001452 0.037530 0.000000 0.037446 0.606641 1 0.00 0.00 -1 6 0 -0.003549 0.220861 0.000000 0.221647 0.050590 # beta 0.00 omega 0.50 degree # index beta omega h k l dr[1/A] R(h,k,l) qv[1/A] qh[1/A] P(omega) 2 0.00 0.50 0 -1 0 0.001779 0.037530 0.000000 -0.037419 0.472496 2 0.00 0.50 0 1 0 0.001126 0.037530 0.000000 0.037469 0.740615 ... |
Check for look-up table
Method cxif5315/plot-lattice-from-data.py generates two plots from look-up table;
- intensities for omega vs. qh
- reconstructed lattice from measured omega and qh
- the same image convolved with 2-d gaussian weight matrix arr_2d_gauss(2, 1.5)
Lattice reconstructed by Meng
Original and convoluted images
...