Refining multi-tile detector metrology with DIALS

Warning

This tutorial is out of date! Default parameters and output have changed.

Introduction

At the end of the Processing in Detail, we showed plots from dials.analyse_output, including images of the positional residuals as a function of position on the detector face. There were clearly some systematic effects, suggesting whole-tile shifts or rotations. Generally, these are small (typically much less than the pixel size) and are expected. DECTRIS provides calibration tables from factory metrology experiments, which XDS can use. Currently we are not applying the correction tables to reflection positions in DIALS. In some situations we may have to work with a much less well calibrated or characterised detector than the Pilatus instruments installed at Diamond. Nevertheless, the multi-panel detector models in DIALS allow us to discover the appropriate corrections assuming we have good diffraction data to refine against.

Here we have access to a thaumatin dataset collected with low dose to avoid radiation damage and with the detector 400 mm from the sample to ensure the coverage of reflections extends right out towards the corners of the images.

Preparing for multi-tile refinement

Usually multi-panel detectors like the Pilatus P6M are treated as if they were single panels (with dead regions) in DIALS. Clearly, if we are to make corrections to the relative positions and orientations of the panels, we need a more detailed detector model in which the panels are treated separately. Currently there isn’t really a neat way to switch between the single and multiple panel models, but until we have a native mechanism for doing this in DIALS there is a hidden workaround that does this. We just need to set an environment variable:

export P6M_60_PANEL=1

With that in place, we start as usual by importing the dataset:

dials.import Thaum_M10S15_3_*.cbf

The fact that the multiple panel detector model was used is clear from the output of:

dials.show datablock.json

We can now inspect images from the dataset with:

dials.image_viewer datablock.json

However, the dataset is at such low dose that it can be difficult to find the spots. It might be helpful to us a more mature image viewer such as ADXV or Albula for this.

We have to run spot finding throughout the dataset. It is a 360 degree sweep so this will take a few minutes. We can use more processes to move a little quicker:

dials.find_spots datablock.json nproc=8 min_spot_size=3

Here we chose to set min_spot_size=3, which overrides the default of 6 used for this detector model. We did this because otherwise this weakly diffracting dataset produces rather few strong spots, whereas we want as much coverage as we can get, right out to the corners

Now to index the data. Although it is a well-diffracting crystal, running indexing with defaults finds an approximate supercell with the a and b lengths slightly more than doubled. Inspection of the found spots with dials.reciprocal_lattice_viewer shows curved lines of spots rather than a regular grid, indicating a poor starting model. There are plenty of things that can be done to get indexing to find the correct cell here. If you have time you might like to try a few options out. However, if not, we won’t worry too much about it, because if we pass in the known unit cell for thaumatin then it works just fine. We also chose to apply tetragonal symmetry immediately:

dials.index datablock.json strong.pickle space_group="P 4" unit_cell="58 58 150 90 90 90"

The output of refinement in the highest resolution macrocycle is as follows:

################################################################################
Starting refinement (macro-cycle 5)
################################################################################


  Summary statistics for observations matched to predictions:
  --------------------------------------------------------------------------
  |                   | Min     | Q1       | Med        | Q3      | Max    |
  --------------------------------------------------------------------------
  | Xc - Xo (mm)      | -2.492  | -0.04261 | -0.0004832 | 0.0424  | 2.131  |
  | Yc - Yo (mm)      | -1.697  | -0.03481 | 0.001733   | 0.03682 | 1.656  |
  | Phic - Phio (deg) | -0.9645 | -0.06454 | -0.0007655 | 0.06379 | 0.8353 |
  | X weights         | 80.55   | 127.8    | 131.8      | 133.9   | 135.2  |
  | Y weights         | 77.02   | 127.6    | 131.8      | 134     | 135.2  |
  | Phi weights       | 26.58   | 32.59    | 32.65      | 32.65   | 32.65  |
  --------------------------------------------------------------------------


  Refinement steps:
  -------------------------------------------------
  | Step | Nref  | RMSD_X   | RMSD_Y   | RMSD_Phi |
  |      |       | (mm)     | (mm)     | (deg)    |
  -------------------------------------------------
  | 0    | 18025 | 0.067273 | 0.063827 | 0.086195 |
  | 1    | 18025 | 0.067369 | 0.063706 | 0.086197 |
  | 2    | 18025 | 0.067439 | 0.063629 | 0.086155 |
  | 3    | 18025 | 0.067531 | 0.063529 | 0.086115 |
  | 4    | 18025 | 0.067587 | 0.063456 | 0.086122 |
  | 5    | 18025 | 0.067607 | 0.063419 | 0.086155 |
  | 6    | 18025 | 0.067617 | 0.063401 | 0.086176 |
  | 7    | 18025 | 0.06762  | 0.063395 | 0.086181 |
  -------------------------------------------------
  RMSD no longer decreasing

  RMSDs by experiment:
  ----------------------------------------------
  | Exp | Nref  | RMSD_X  | RMSD_Y  | RMSD_Z   |
  |     |       | (px)    | (px)    | (images) |
  ----------------------------------------------
  | 0   | 18025 | 0.39314 | 0.36858 | 0.24623  |
  ----------------------------------------------

  RMSDs by panel:
  -----------------------------------------------
  | Panel | Nref | RMSD_X  | RMSD_Y  | RMSD_Z   |
  |       |      | (px)    | (px)    | (images) |
  -----------------------------------------------
  | 0     | 10   | 0.53322 | 0.82846 | 0.14379  |
  | 1     | 140  | 0.37822 | 0.43153 | 0.26887  |
  | 2     | 297  | 0.41585 | 0.40848 | 0.26118  |
  | 3     | 146  | 0.5563  | 0.40867 | 0.27493  |
  | 4     | 17   | 0.49115 | 0.43483 | 0.21656  |
  | 5     | 44   | 0.41713 | 0.46071 | 0.27804  |
  | 6     | 405  | 0.44932 | 0.51832 | 0.26569  |
  | 7     | 758  | 0.33153 | 0.28766 | 0.25461  |
  | 8     | 478  | 0.38555 | 0.40525 | 0.26244  |
  | 9     | 96   | 0.47338 | 0.45054 | 0.24652  |
  | 10    | 152  | 0.40191 | 0.959   | 0.2976   |
  | 11    | 701  | 0.31004 | 0.3035  | 0.24597  |
  | 12    | 1128 | 0.37366 | 0.24303 | 0.24016  |
  | 13    | 802  | 0.33404 | 0.32094 | 0.24614  |
  | 14    | 201  | 0.50746 | 0.47124 | 0.27169  |
  | 15    | 231  | 0.40185 | 0.56122 | 0.3029   |
  | 16    | 745  | 0.21926 | 0.30075 | 0.21426  |
  | 17    | 831  | 0.17864 | 0.17706 | 0.2176   |
  | 18    | 778  | 0.36795 | 0.21167 | 0.20696  |
  | 19    | 269  | 0.44757 | 0.40028 | 0.23809  |
  | 20    | 205  | 0.48078 | 0.65283 | 0.34241  |
  | 21    | 467  | 0.37577 | 0.35789 | 0.20992  |
  | 22    | 370  | 0.22329 | 0.29697 | 0.20943  |
  | 23    | 447  | 0.36708 | 0.22722 | 0.20176  |
  | 24    | 231  | 0.28466 | 0.52795 | 0.31729  |
  | 25    | 3    | 0.36377 | 1.1547  | 0.70452  |
  | 26    | 28   | 0.14842 | 0.29869 | 0.23437  |
  | 27    | 22   | 0.13375 | 0.2666  | 0.17175  |
  | 28    | 17   | 0.1222  | 0.14108 | 0.17261  |
  | 29    | 4    | 0.40725 | 0.85708 | 0.42044  |
  | 30    | 145  | 0.60805 | 0.57195 | 0.43388  |
  | 31    | 281  | 0.18353 | 0.21351 | 0.2061   |
  | 32    | 250  | 0.19984 | 0.15417 | 0.1993   |
  | 33    | 294  | 0.15006 | 0.27665 | 0.21823  |
  | 34    | 158  | 0.44084 | 0.50436 | 0.38365  |
  | 35    | 203  | 0.86526 | 0.47248 | 0.32568  |
  | 36    | 525  | 0.54992 | 0.31614 | 0.19036  |
  | 37    | 521  | 0.22611 | 0.12494 | 0.21894  |
  | 38    | 624  | 0.23105 | 0.20793 | 0.22246  |
  | 39    | 259  | 0.43408 | 0.43312 | 0.2724   |
  | 40    | 134  | 0.55289 | 0.3935  | 0.26933  |
  | 41    | 557  | 0.29543 | 0.24331 | 0.23127  |
  | 42    | 771  | 0.46252 | 0.19655 | 0.2348   |
  | 43    | 658  | 0.24901 | 0.31232 | 0.24217  |
  | 44    | 199  | 0.44994 | 0.43547 | 0.28695  |
  | 45    | 54   | 0.5679  | 0.44719 | 0.23983  |
  | 46    | 333  | 0.44117 | 0.44019 | 0.24586  |
  | 47    | 592  | 0.3056  | 0.44653 | 0.24795  |
  | 48    | 411  | 0.67885 | 0.39219 | 0.2588   |
  | 49    | 93   | 0.45355 | 0.5659  | 0.2611   |
  | 50    | 5    | 1.0009  | 0.46232 | 0.21423  |
  | 51    | 162  | 0.72162 | 0.36619 | 0.23884  |
  | 52    | 324  | 0.51574 | 0.48569 | 0.2488   |
  | 53    | 224  | 0.41145 | 0.40151 | 0.2671   |
  | 54    | 16   | 0.5362  | 0.50568 | 0.21454  |
  | 55    | 1    | 0.85351 | 0.57076 | 0.24389  |
  | 56    | 37   | 0.72428 | 0.38373 | 0.27121  |
  | 57    | 106  | 0.80208 | 0.98469 | 0.27939  |
  | 58    | 61   | 0.5144  | 1.3281  | 0.30346  |
  | 59    | 4    | 0.33371 | 0.50058 | 0.24425  |
  -----------------------------------------------
  Final refined crystal models:
  model 1 (192715 reflections):
  Crystal:
      Unit cell: (57.834, 57.834, 150.022, 90.000, 90.000, 90.000)
      Space group: P 4
      U matrix:  {{ 0.4122, -0.9018,  0.1299},
                  { 0.2361, -0.0320, -0.9712},
                  { 0.8800,  0.4310,  0.1997}}
      B matrix:  {{ 0.0173,  0.0000,  0.0000},
                  { 0.0000,  0.0173,  0.0000},
                  { 0.0000,  0.0000,  0.0067}}
      A = UB:    {{ 0.0071, -0.0156,  0.0009},
                  { 0.0041, -0.0006, -0.0065},
                  { 0.0152,  0.0075,  0.0013}}

This refinement was performed moving all the panels as a rigid block, as usual. With overall positional RMSDs within 40% of the pixel size and a quarter of the image width in \phi we can see straight away that we are dealing with a fairly good quality dataset. There are a few outliers of well over 1 mm on the detector surface and nearly 1 degree in \phi though, which we would prefer not to include in refinement. The outliers are not as bad if we had kept min_spot_size=6, but the detector coverage is worse in that case. Although from the indexing results it seems that coverage of reflections on the outer panels is rather low, so far we let refinement take a random subset of the data in order to index quicker, so there’s no need to worry about that yet.

Now we will refine the detector as a rigid block again, turning on outlier rejection and requesting to use all reflections to get the best we can out of the dataset. We will also keep the refined reflections file for analysis. The final parameter here, close_to_spindle_cutoff=0.01 allows reflections closer to the spindle to be included in refinement (default value is 0.05, and if set to 0.0 no reflections will be rejected for being too close). Without this option, the central panels are very sparse:

dials.refine indexed.pickle experiments.json \
 outlier.algorithm=tukey use_all_reflections=true close_to_spindle_cutoff=0.01 \
 output.reflections=refined_reflections_lev0.pickle \
 output.experiments=refined_experiments_lev0.json

Here is the output:

The following parameters have been modified:

output {
  reflections = refined_reflections_lev0.pickle
  experiments = refined_experiments_lev0.json
}
refinement {
  reflections {
    use_all_reflections = true
    close_to_spindle_cutoff = 0.01
    outlier {
      algorithm = null *tukey
    }
  }
}
input {
  experiments = experiments.json
  reflections = indexed.pickle
}

Configuring refiner

Summary statistics for observations matched to predictions:
-----------------------------------------------------------------------
|                   | Min    | Q1       | Med       | Q3      | Max   |
-----------------------------------------------------------------------
| Xc - Xo (mm)      | -2.487 | -0.04174 | 0.000276  | 0.0424  | 2.128 |
| Yc - Yo (mm)      | -1.75  | -0.03468 | 0.00168   | 0.03656 | 1.654 |
| Phic - Phio (deg) | -5.549 | -0.06627 | -0.002108 | 0.063   | 3.27  |
| X weights         | 80.55  | 127.9    | 131.9     | 134     | 135.2 |
| Y weights         | 77.02  | 127.7    | 131.9     | 134     | 135.2 |
| Phi weights       | 25.97  | 32.57    | 32.65     | 32.65   | 32.65 |
-----------------------------------------------------------------------

6375 reflections have been rejected as outliers

Summary statistics for observations matched to predictions:
--------------------------------------------------------------------------
|                   | Min     | Q1       | Med        | Q3      | Max    |
--------------------------------------------------------------------------
| Xc - Xo (mm)      | -0.2933 | -0.04128 | -3.167e-05 | 0.04112 | 0.3725 |
| Yc - Yo (mm)      | -0.4707 | -0.03376 | 0.001821   | 0.03588 | 0.4945 |
| Phic - Phio (deg) | -0.6965 | -0.06574 | -0.002213  | 0.06229 | 0.7919 |
| X weights         | 80.55   | 128      | 132        | 134     | 135.2  |
| Y weights         | 77.02   | 127.8    | 131.9      | 134.1   | 135.2  |
| Phi weights       | 26.58   | 32.57    | 32.65      | 32.65   | 32.65  |
--------------------------------------------------------------------------

Performing refinement...

Refinement steps:
--------------------------------------------------
| Step | Nref   | RMSD_X   | RMSD_Y   | RMSD_Phi |
|      |        | (mm)     | (mm)     | (deg)    |
--------------------------------------------------
| 0    | 186203 | 0.064091 | 0.057786 | 0.086136 |
| 1    | 186203 | 0.064045 | 0.057829 | 0.08608  |
| 2    | 186203 | 0.063949 | 0.05791  | 0.086068 |
| 3    | 186203 | 0.063825 | 0.058023 | 0.086026 |
| 4    | 186203 | 0.063734 | 0.058114 | 0.085958 |
| 5    | 186203 | 0.063682 | 0.058167 | 0.085909 |
| 6    | 186203 | 0.063654 | 0.058198 | 0.085887 |
| 7    | 186203 | 0.063645 | 0.058208 | 0.085882 |
| 8    | 186203 | 0.063644 | 0.05821  | 0.085882 |
--------------------------------------------------
RMSD no longer decreasing

RMSDs by experiment:
-----------------------------------------------
| Exp | Nref   | RMSD_X  | RMSD_Y  | RMSD_Z   |
|     |        | (px)    | (px)    | (images) |
-----------------------------------------------
| 0   | 186203 | 0.37002 | 0.33843 | 0.24538  |
-----------------------------------------------

RMSDs by panel:
-----------------------------------------------
| Panel | Nref | RMSD_X  | RMSD_Y  | RMSD_Z   |
|       |      | (px)    | (px)    | (images) |
-----------------------------------------------
| 0     | 63   | 0.47216 | 0.59331 | 0.22352  |
| 1     | 1363 | 0.35019 | 0.46301 | 0.26549  |
| 2     | 3076 | 0.4128  | 0.40658 | 0.26659  |
| 3     | 1630 | 0.46232 | 0.39692 | 0.25547  |
| 4     | 112  | 0.39711 | 0.42575 | 0.23841  |
| 5     | 394  | 0.35145 | 0.51263 | 0.25474  |
| 6     | 3489 | 0.3238  | 0.48667 | 0.25694  |
| 7     | 6172 | 0.30869 | 0.27402 | 0.25323  |
| 8     | 3908 | 0.37309 | 0.38322 | 0.25577  |
| 9     | 680  | 0.46698 | 0.44533 | 0.24812  |
| 10    | 1178 | 0.33544 | 0.51129 | 0.26289  |
| 11    | 5690 | 0.28409 | 0.29009 | 0.24376  |
| 12    | 9385 | 0.36731 | 0.23495 | 0.23873  |
| 13    | 6878 | 0.31962 | 0.31088 | 0.23691  |
| 14    | 1749 | 0.43068 | 0.47078 | 0.25696  |
| 15    | 1870 | 0.29306 | 0.46181 | 0.26494  |
| 16    | 6328 | 0.19137 | 0.29002 | 0.21465  |
| 17    | 7599 | 0.17418 | 0.1739  | 0.22093  |
| 18    | 7072 | 0.36776 | 0.18662 | 0.20821  |
| 19    | 2609 | 0.37377 | 0.42663 | 0.24206  |
| 20    | 1806 | 0.40342 | 0.57026 | 0.32007  |
| 21    | 4247 | 0.35488 | 0.34512 | 0.19181  |
| 22    | 3558 | 0.21403 | 0.30996 | 0.20695  |
| 23    | 4103 | 0.36559 | 0.19154 | 0.17851  |
| 24    | 2270 | 0.25744 | 0.48077 | 0.26616  |
| 25    | 480  | 0.38505 | 1.0035  | 0.73284  |
| 26    | 930  | 0.15033 | 0.26661 | 0.34611  |
| 27    | 971  | 0.11836 | 0.19105 | 0.18709  |
| 28    | 975  | 0.11907 | 0.25901 | 0.30046  |
| 29    | 549  | 0.41024 | 0.95027 | 0.72705  |
| 30    | 1478 | 0.56625 | 0.3835  | 0.34664  |
| 31    | 3216 | 0.16373 | 0.19856 | 0.19257  |
| 32    | 2826 | 0.20149 | 0.14689 | 0.19485  |
| 33    | 3285 | 0.12942 | 0.30274 | 0.20848  |
| 34    | 1965 | 0.38874 | 0.45564 | 0.38875  |
| 35    | 1934 | 0.61707 | 0.35339 | 0.25067  |
| 36    | 5875 | 0.55413 | 0.29556 | 0.19246  |
| 37    | 5992 | 0.2137  | 0.12029 | 0.21502  |
| 38    | 6589 | 0.2254  | 0.21047 | 0.21193  |
| 39    | 2763 | 0.39851 | 0.41673 | 0.27272  |
| 40    | 1512 | 0.50654 | 0.36893 | 0.25051  |
| 41    | 6476 | 0.28308 | 0.21635 | 0.22668  |
| 42    | 9212 | 0.40577 | 0.18283 | 0.2252   |
| 43    | 7521 | 0.22936 | 0.29881 | 0.23864  |
| 44    | 2392 | 0.40816 | 0.42881 | 0.2751   |
| 45    | 611  | 0.56917 | 0.48084 | 0.26191  |
| 46    | 4379 | 0.42672 | 0.42447 | 0.25285  |
| 47    | 7749 | 0.29309 | 0.4378  | 0.24631  |
| 48    | 5470 | 0.67537 | 0.37528 | 0.25187  |
| 49    | 1210 | 0.42778 | 0.53938 | 0.27266  |
| 50    | 127  | 0.92194 | 0.45501 | 0.24629  |
| 51    | 2235 | 0.71785 | 0.33424 | 0.25781  |
| 52    | 4379 | 0.41055 | 0.38155 | 0.25735  |
| 53    | 2757 | 0.44268 | 0.40149 | 0.26646  |
| 54    | 327  | 0.52856 | 0.51306 | 0.28205  |
| 56    | 457  | 0.60967 | 0.429   | 0.26482  |
| 57    | 1507 | 0.70752 | 0.67612 | 0.25844  |
| 58    | 786  | 0.53669 | 0.50635 | 0.27127  |
| 59    | 39   | 0.43205 | 0.65796 | 0.25585  |
-----------------------------------------------
Saving refined experiments to refined_experiments_lev0.json
Saving refined reflections to refined_reflections_lev0.pickle

Outlier rejection has cleaned up the positional residuals so now the greatest deviation is within 0.4 mm of the predicted position. The angular extreme is now just over 0.4 degrees. Coverage of the outer and central panels (where reflections are in the backstop shadow or thrown away for being too close to the spindle) is still a little low. Notably, panel 55 (a corner panel) is completely missing. If we had more datasets recorded at the same detector distance (and more time to process them) we could combine them in a multi-crystal joint refinement job to increase the coverage of panels further. However, for the purposes of this tutorial we will see what we can get with this single dataset.

Before moving on to the multi-panel refinement job we will take a look at the refined reflections file:

dials.analyse_output refined_reflections_lev0.pickle grid_size=5,12

Here we had to tell dials.analyse_output about the arrangement of the panels, as it does not use the refined_experiments_lev0.json file so cannot figure this out itself.

Here are the positional residual plots for X and Y, analysis/centroid/centroid_diff_x.png and analysis/centroid/centroid_diff_y.png. The multi-panel versions of these plots are not as compact as the single tile version presented in the Processing in Detail. However, careful comparison of the plots is enough to show that the same pattern of shifts is present.

../../_images/centroid_diff_x_multi_panel_lev0.png ../../_images/centroid_diff_y_multi_panel_lev0.png

Multi-tile refinement

Now we repeat refinement, but we allow the panels to move independently. In DIALS multi-panel detectors are represented by a hierarchical model. The highest level hieararchy_level=0 means to treat the whole detector unit as a rigid block. Some detectors, notably the CS-PAD used at LCLS beamlines, have a real hierarchy of a few levels deep. The Pilatus P6M has a very simple hierarchy, with a single lower level, hieararchy_level=1, in which every panel is treated separately. We now start from the previous refinement run specifying this hierarchy level:

dials.refine indexed.pickle refined_experiments_lev0.json outlier.algorithm=tukey \
 use_all_reflections=true output.reflections=refined_reflections_lev1.pickle \
 close_to_spindle_cutoff=0.01 bin_size_fraction=0 hierarchy_level=1 \
 output.experiments=refined_experiments_lev1.json

You may have noticed that apart from hierarchy_level=1 there was an additional parameter added to this command compared to the previous refinement run, namely bin_size_fraction=0. This sets the RMSD target for refinement to zero, so that refinement will never terminate because the RMSDs are ‘good enough’, only if they converge so that their rate of decrease on subsequent steps falls to zero. This is necessary because the extra freedom allowed by parameterising each panel individually allows the RMSDs to fall lower than the default target. There are 366 parameters in total for this refinement run. This can be seen by checking the file dials.refine.debug.log once refinement is underway.

Warning

This job took 17 minutes to run on a Linux desktop with a Core i7 CPU running at 3.07GHz, and uses about 4 GB of RAM.

Refinement is single-process at the moment, unfortunately, so we can’t yet make use of parallelism here to speed the job up. The output is as follows:

The following parameters have been modified:

output {
  experiments = refined_experiments_lev1.json
  reflections = refined_reflections_lev1.pickle
}
refinement {
  parameterisation {
    detector {
      hierarchy_level = 1
    }
  }
  target {
    bin_size_fraction = 0
  }
  reflections {
    use_all_reflections = true
    close_to_spindle_cutoff = 0.01
    outlier {
      algorithm = null *tukey
    }
  }
}
input {
  experiments = refined_experiments.json
  reflections = indexed.pickle
}

Configuring refiner

Summary statistics for observations matched to predictions:
------------------------------------------------------------------------
|                   | Min    | Q1       | Med        | Q3      | Max   |
------------------------------------------------------------------------
| Xc - Xo (mm)      | -2.496 | -0.04178 | 0.0004518  | 0.04222 | 2.133 |
| Yc - Yo (mm)      | -1.903 | -0.03577 | 0.0006705  | 0.03588 | 1.656 |
| Phic - Phio (deg) | -5.576 | -0.06467 | -0.0007391 | 0.06414 | 3.292 |
| X weights         | 80.55  | 127.9    | 131.9      | 134     | 135.2 |
| Y weights         | 77.02  | 127.7    | 131.9      | 134     | 135.2 |
| Phi weights       | 25.97  | 32.57    | 32.65      | 32.65   | 32.65 |
------------------------------------------------------------------------

6433 reflections have been rejected as outliers

Summary statistics for observations matched to predictions:
--------------------------------------------------------------------------
|                   | Min     | Q1       | Med        | Q3      | Max    |
--------------------------------------------------------------------------
| Xc - Xo (mm)      | -0.2916 | -0.04138 | 0.0001204  | 0.041   | 0.3683 |
| Yc - Yo (mm)      | -0.4838 | -0.03481 | 0.0008164  | 0.03518 | 0.4917 |
| Phic - Phio (deg) | -0.6969 | -0.06416 | -0.0008647 | 0.06351 | 0.7644 |
| X weights         | 80.55   | 128      | 132        | 134     | 135.2  |
| Y weights         | 77.02   | 127.8    | 131.9      | 134.1   | 135.2  |
| Phi weights       | 26.58   | 32.57    | 32.65      | 32.65   | 32.65  |
--------------------------------------------------------------------------

Performing refinement...

Refinement steps:
--------------------------------------------------
| Step | Nref   | RMSD_X   | RMSD_Y   | RMSD_Phi |
|      |        | (mm)     | (mm)     | (deg)    |
--------------------------------------------------
| 0    | 186145 | 0.063617 | 0.058127 | 0.085801 |
| 1    | 186145 | 0.056976 | 0.05538  | 0.085719 |
| 2    | 186145 | 0.049808 | 0.052619 | 0.085597 |
| 3    | 186145 | 0.04634  | 0.051408 | 0.085475 |
| 4    | 186145 | 0.045568 | 0.051142 | 0.085391 |
| 5    | 186145 | 0.04538  | 0.051014 | 0.085337 |
| 6    | 186145 | 0.045228 | 0.050729 | 0.085285 |
| 7    | 186145 | 0.045054 | 0.050076 | 0.085215 |
| 8    | 186145 | 0.044868 | 0.049034 | 0.085115 |
| 9    | 186145 | 0.044746 | 0.048167 | 0.085025 |
| 10   | 186145 | 0.044708 | 0.047833 | 0.08498  |
| 11   | 186145 | 0.044695 | 0.047759 | 0.084969 |
| 12   | 186145 | 0.04468  | 0.047738 | 0.084967 |
| 13   | 186145 | 0.044668 | 0.047726 | 0.084967 |
| 14   | 186145 | 0.044662 | 0.047721 | 0.084967 |
| 15   | 186145 | 0.04466  | 0.047719 | 0.084967 |
--------------------------------------------------
RMSD no longer decreasing

RMSDs by experiment:
-----------------------------------------------
| Exp | Nref   | RMSD_X  | RMSD_Y  | RMSD_Z   |
|     |        | (px)    | (px)    | (images) |
-----------------------------------------------
| 0   | 186145 | 0.25965 | 0.27743 | 0.24276  |
-----------------------------------------------

RMSDs by panel:
------------------------------------------------
| Panel | Nref | RMSD_X   | RMSD_Y  | RMSD_Z   |
|       |      | (px)     | (px)    | (images) |
------------------------------------------------
| 0     | 64   | 0.32036  | 0.57819 | 0.22482  |
| 1     | 1361 | 0.34304  | 0.4111  | 0.26504  |
| 2     | 3079 | 0.3561   | 0.36483 | 0.26682  |
| 3     | 1632 | 0.37843  | 0.37838 | 0.25453  |
| 4     | 112  | 0.39065  | 0.39675 | 0.23706  |
| 5     | 394  | 0.30539  | 0.48132 | 0.25442  |
| 6     | 3490 | 0.29682  | 0.35249 | 0.25661  |
| 7     | 6178 | 0.28525  | 0.27486 | 0.25318  |
| 8     | 3909 | 0.31728  | 0.30551 | 0.25565  |
| 9     | 675  | 0.39465  | 0.42802 | 0.24507  |
| 10    | 1175 | 0.30057  | 0.46365 | 0.26215  |
| 11    | 5685 | 0.23485  | 0.27143 | 0.24326  |
| 12    | 9386 | 0.20171  | 0.18887 | 0.23872  |
| 13    | 6880 | 0.23494  | 0.23081 | 0.23689  |
| 14    | 1746 | 0.31911  | 0.40137 | 0.25357  |
| 15    | 1870 | 0.27575  | 0.44911 | 0.26334  |
| 16    | 6321 | 0.18338  | 0.21691 | 0.21382  |
| 17    | 7608 | 0.14673  | 0.13459 | 0.22118  |
| 18    | 7075 | 0.15323  | 0.1772  | 0.20816  |
| 19    | 2608 | 0.27977  | 0.40106 | 0.23663  |
| 20    | 1803 | 0.24748  | 0.54944 | 0.31271  |
| 21    | 4249 | 0.15151  | 0.21436 | 0.18971  |
| 22    | 3560 | 0.11336  | 0.12431 | 0.20676  |
| 23    | 4103 | 0.10087  | 0.16239 | 0.17534  |
| 24    | 2271 | 0.20061  | 0.46907 | 0.24929  |
| 25    | 476  | 0.26774  | 0.94809 | 0.69743  |
| 26    | 927  | 0.10167  | 0.22193 | 0.33398  |
| 27    | 974  | 0.097735 | 0.13559 | 0.1876   |
| 28    | 977  | 0.069491 | 0.1907  | 0.29589  |
| 29    | 549  | 0.19349  | 0.91362 | 0.67233  |
| 30    | 1474 | 0.31601  | 0.36268 | 0.32879  |
| 31    | 3209 | 0.13728  | 0.15038 | 0.18643  |
| 32    | 2829 | 0.10162  | 0.11078 | 0.19453  |
| 33    | 3276 | 0.10363  | 0.17357 | 0.2043   |
| 34    | 1951 | 0.20521  | 0.4038  | 0.36983  |
| 35    | 1931 | 0.37408  | 0.28915 | 0.24379  |
| 36    | 5869 | 0.18862  | 0.14476 | 0.1912   |
| 37    | 5991 | 0.13887  | 0.11716 | 0.21506  |
| 38    | 6591 | 0.15001  | 0.19861 | 0.21179  |
| 39    | 2763 | 0.25238  | 0.36458 | 0.26476  |
| 40    | 1511 | 0.47841  | 0.27181 | 0.24826  |
| 41    | 6461 | 0.27368  | 0.18468 | 0.22587  |
| 42    | 9213 | 0.19848  | 0.15776 | 0.22531  |
| 43    | 7521 | 0.21923  | 0.24706 | 0.2383   |
| 44    | 2390 | 0.32645  | 0.40052 | 0.27085  |
| 45    | 611  | 0.55212  | 0.29673 | 0.26008  |
| 46    | 4378 | 0.38894  | 0.24988 | 0.25198  |
| 47    | 7747 | 0.2824   | 0.2289  | 0.24614  |
| 48    | 5470 | 0.29905  | 0.31716 | 0.25135  |
| 49    | 1210 | 0.41311  | 0.47572 | 0.2691   |
| 50    | 127  | 0.69191  | 0.39228 | 0.2447   |
| 51    | 2231 | 0.47797  | 0.32129 | 0.25717  |
| 52    | 4378 | 0.38899  | 0.31331 | 0.25676  |
| 53    | 2758 | 0.39246  | 0.40955 | 0.26575  |
| 54    | 327  | 0.48186  | 0.52671 | 0.2798   |
| 56    | 458  | 0.57234  | 0.42803 | 0.2643   |
| 57    | 1507 | 0.494    | 0.41177 | 0.25741  |
| 58    | 786  | 0.48028  | 0.48906 | 0.26898  |
| 59    | 40   | 0.39499  | 0.38962 | 0.2466   |
------------------------------------------------
Saving refined experiments to refined_experiments_lev1.json
Saving refined reflections to refined_reflections_lev1.pickle

Following refinement, we repeat the analysis of positional residuals:

mv analysis analysis_lev0
dials.analyse_output refined_reflections.pickle grid_size=5,12
mv analysis analysis_lev1

The positional residual plots for X and Y, analysis_lev1/centroid/centroid_diff_x.png and analysis_lev1/centroid/centroid_diff_y.png make it clear that despite poor coverage on some panels, the systematic shifts have been cleaned up by the refinement job.

../../_images/centroid_diff_x_multi_panel_lev1.png ../../_images/centroid_diff_y_multi_panel_lev1.png

Applying the corrected metrology

Although we have corrected for small shifts and rotations apparently present from spot positions recorded in one dataset, what we would really like to do is apply these corrections to a different dataset. Short of rewriting the dxtbx Format object for the detector to incorporate the corrections, we could try to ‘copy and paste’ the detector from one dataset to the other.

We will choose the standard tutorial data to try this, from the Processing in Detail. First we have to process that data using the multi-panel version of the Pilatus P6M detector model. Assuming the environment variable P6M_60_PANEL=1 is set in this terminal we just need to repeat the commands from the tutorial:

mkdir tutorial_data
cd !$
dials.import /path/to/th_8_2*cbf
dials.find_spots datablock.json nproc=4
dials.index datablock.json strong.pickle space_group="P4"
dials.refine experiments.json indexed.pickle outlier.algorithm=tukey use_all_reflections=true bin_size_fraction=0.0

Note these are the overall RMSDs (comparable to the results from the Processing in Detail, as we’d expect):

RMSDs by experiment:
----------------------------------------------
| Exp | Nref  | RMSD_X  | RMSD_Y  | RMSD_Z   |
|     |       | (px)    | (px)    | (images) |
----------------------------------------------
| 0   | 53916 | 0.26341 | 0.21858 | 0.106    |
----------------------------------------------

Now we do the scan-varying refinement and integrate:

dials.refine refined_experiments.json indexed.pickle outlier.algorithm=tukey use_all_reflections=true bin_size_fraction=0.0 scan_varying=true output.experiments=sv_refined_experiments.json
dials.integrate sv_refined_experiments.json indexed.pickle outlier.algorithm=null nproc=4
dials.export integrated.pickle sv_refined_experiments.json mtz.hklout=integrated.mtz ignore_panels=true
dials.analyse_output integrated.pickle grid_size=5,12

From the end of dials.integrate.log:

Summary of integration results binned by resolution
----------------------------------------------------------------------------------------------------------
d min |  d max | # full | # part | # over | # ice | # sum | # prf | <Ibg> | <I/sigI> | <I/sigI> | <CC prf>
      |        |        |        |        |       |       |       |       |    (sum) |    (prf) |
----------------------------------------------------------------------------------------------------------
 1.17 |   1.19 |    324 |      3 |      0 |     0 |   327 |   260 |  0.04 |     0.37 |     0.54 |     0.11
 1.19 |   1.21 |   1093 |      7 |      0 |     0 |  1100 |   984 |  0.04 |     0.43 |     0.54 |     0.09
 1.21 |   1.23 |   2317 |     14 |      0 |     0 |  2331 |  2138 |  0.05 |     0.48 |     0.57 |     0.10
 1.23 |   1.26 |   3776 |     27 |      0 |     0 |  3803 |  3581 |  0.05 |     0.53 |     0.66 |     0.12
 1.26 |   1.28 |   5416 |     34 |      0 |     0 |  5450 |  5100 |  0.05 |     0.59 |     0.74 |     0.14
 1.28 |   1.31 |   7235 |     54 |      0 |     0 |  7289 |  6806 |  0.06 |     0.64 |     0.81 |     0.16
 1.31 |   1.35 |   9504 |     63 |      0 |     0 |  9567 |  8952 |  0.06 |     0.77 |     0.95 |     0.19
 1.35 |   1.38 |  12495 |     87 |      0 |     0 | 12582 | 11826 |  0.07 |     0.90 |     1.11 |     0.21
 1.38 |   1.42 |  16967 |    111 |      0 |     0 | 17078 | 16125 |  0.08 |     0.98 |     1.20 |     0.23
 1.42 |   1.47 |  20187 |    153 |      0 |     0 | 20340 | 19437 |  0.08 |     1.20 |     1.43 |     0.27
 1.47 |   1.52 |  23623 |    258 |      0 |     0 | 23881 | 22761 |  0.09 |     1.47 |     1.71 |     0.30
 1.52 |   1.58 |  24126 |    271 |      0 |     0 | 24397 | 23289 |  0.10 |     1.75 |     2.03 |     0.35
 1.58 |   1.66 |  25506 |    281 |      0 |     0 | 25787 | 24863 |  0.11 |     2.17 |     2.46 |     0.41
 1.66 |   1.74 |  24269 |    241 |      0 |     0 | 24510 | 23389 |  0.12 |     2.71 |     3.01 |     0.46
 1.74 |   1.85 |  24780 |    238 |      0 |     0 | 25018 | 24019 |  0.15 |     3.50 |     3.82 |     0.52
 1.85 |   2.00 |  25710 |    278 |      0 |     0 | 25988 | 25059 |  0.19 |     4.87 |     5.20 |     0.59
 2.00 |   2.20 |  24757 |    274 |      0 |     0 | 25031 | 24000 |  0.26 |     6.58 |     6.94 |     0.66
 2.20 |   2.51 |  25673 |    252 |      0 |     0 | 25925 | 25119 |  0.29 |     8.84 |     9.15 |     0.69
 2.51 |   3.17 |  25214 |    297 |      0 |     0 | 25511 | 24607 |  0.37 |    12.78 |    13.01 |     0.73
 3.17 | 151.26 |  25766 |    303 |      0 |     0 | 26069 | 25226 |  0.45 |    25.48 |    25.36 |     0.74
----------------------------------------------------------------------------------------------------------

Summary of integration results for the whole dataset
----------------------------------------------
Number fully recorded                 | 338720
Number partially recorded             | 4747
Number with overloaded pixels         | 0
Number in powder rings                | 0
Number processed with summation       | 331984
Number processed with profile fitting | 317541
<Ibg>                                 | 0.18
<I/sigI> (summation)                  | 5.63
<I/sigI> (profile fitting)            | 5.92
<CC prf>                              | 0.43
----------------------------------------------

The plot of correlations between strong (reference) reflections and the profiles used to fit them shows up the systematic effects of uncorrected tile shifts and misorientations:

../../_images/reference_corr_vs_xy_multipanel_uncorrected.png

Now, how can we apply the metrology? Here we will use dials.combine_experiments with the reference_from_experiment.detector option to overwrite the detector model from our experiments.json. We don’t really want the combined experiments file, only this side-effect, so we immediately split it again:

dials.combine_experiments experiments=../refined_experiments_lev1.json experiments=refined_experiments.json reflections=../refined_reflections_lev1.pickle reflections=indexed.pickle reference_from_experiment.detector=0
dials.split_experiments combined_experiments.json combined_reflections.pickle

This results in a few files, of which experiments_1.json is interesting. It contains the updated detector - but beware the detector distance is now completely wrong! It is at about 400 mm rather than 265 mm. We could correct that by editing experiments_1.json directly, but actually there is no need. dials.refine is extremely forgiving of bad starting geometry, though we should remember to fix the beam and crystal models:

dials.refine experiments_1.json indexed.pickle output.experiments=corrected_refined_experiments.json beam.fix=all crystal.fix=all

A snippet from the log file shows that the detector distance offset was largely corrected in a single step:

Refinement steps:
------------------------------------------------
| Step | Nref | RMSD_X   | RMSD_Y   | RMSD_Phi |
|      |      | (mm)     | (mm)     | (deg)    |
------------------------------------------------
| 0    | 4049 | 40.622   | 40.688   | 0.019202 |
| 1    | 4049 | 0.28874  | 0.30251  | 0.019202 |
| 2    | 4049 | 0.045937 | 0.045205 | 0.019202 |
------------------------------------------------

Now we’ll let the crystal and beam refine along with the new detector to RMSD convergence:

dials.refine corrected_refined_experiments.json indexed.pickle outlier.algorithm=tukey use_all_reflections=true bin_size_fraction=0.0 output.experiments=corrected_refined_experiments.json

Here is the output:

The following parameters have been modified:

output {
  experiments = corrected_refined_experiments.json
}
refinement {
  target {
    bin_size_fraction = 0.0
  }
  reflections {
    use_all_reflections = true
    outlier {
      algorithm = null *tukey
    }
  }
}
input {
  experiments = corrected_refined_experiments.json
  reflections = indexed.pickle
}

Configuring refiner

Summary statistics for observations matched to predictions:
-------------------------------------------------------------------------
|                   | Min     | Q1       | Med       | Q3      | Max    |
-------------------------------------------------------------------------
| Xc - Xo (mm)      | -0.4701 | -0.02988 | -0.01025  | 0.02159 | 0.6902 |
| Yc - Yo (mm)      | -0.7267 | -0.02846 | -0.005739 | 0.02339 | 1.247  |
| Phic - Phio (deg) | -1.409  | -0.01035 | 5.562e-05 | 0.01104 | 0.4029 |
| X weights         | 113.5   | 133.7    | 134.5     | 135     | 135.2  |
| Y weights         | 109.2   | 133.4    | 134.5     | 135     | 135.2  |
| Phi weights       | 155.5   | 175.8    | 177       | 177.5   | 177.8  |
-------------------------------------------------------------------------

2382 reflections have been rejected as outliers

Summary statistics for observations matched to predictions:
--------------------------------------------------------------------------
|                   | Min     | Q1        | Med       | Q3      | Max    |
--------------------------------------------------------------------------
| Xc - Xo (mm)      | -0.1148 | -0.02997  | -0.01092  | 0.01966 | 0.331  |
| Yc - Yo (mm)      | -0.1397 | -0.02791  | -0.005742 | 0.02263 | 0.2531 |
| Phic - Phio (deg) | -0.1881 | -0.009998 | 0.0001419 | 0.0109  | 0.1176 |
| X weights         | 113.5   | 133.7     | 134.6     | 135     | 135.2  |
| Y weights         | 109.9   | 133.5     | 134.5     | 135     | 135.2  |
| Phi weights       | 158.2   | 175.9     | 177       | 177.5   | 177.8  |
--------------------------------------------------------------------------

Performing refinement...

Refinement steps:
-------------------------------------------------
| Step | Nref  | RMSD_X   | RMSD_Y   | RMSD_Phi |
|      |       | (mm)     | (mm)     | (deg)    |
-------------------------------------------------
| 0    | 53738 | 0.04524  | 0.04251  | 0.015875 |
| 1    | 53738 | 0.043495 | 0.041514 | 0.015395 |
| 2    | 53738 | 0.042657 | 0.040762 | 0.01535  |
| 3    | 53738 | 0.040455 | 0.038826 | 0.015287 |
| 4    | 53738 | 0.0356   | 0.034575 | 0.015251 |
| 5    | 53738 | 0.028146 | 0.027998 | 0.015246 |
| 6    | 53738 | 0.021919 | 0.022913 | 0.015245 |
| 7    | 53738 | 0.019867 | 0.021931 | 0.015256 |
| 8    | 53738 | 0.019579 | 0.02202  | 0.015267 |
| 9    | 53738 | 0.019553 | 0.02204  | 0.015269 |
| 10   | 53738 | 0.019552 | 0.022041 | 0.015269 |
-------------------------------------------------
RMSD no longer decreasing

RMSDs by experiment:
----------------------------------------------
| Exp | Nref  | RMSD_X  | RMSD_Y  | RMSD_Z   |
|     |       | (px)    | (px)    | (images) |
----------------------------------------------
| 0   | 53738 | 0.11368 | 0.12815 | 0.10179  |
----------------------------------------------

RMSDs by panel:
-------------------------------------------------
| Panel | Nref | RMSD_X   | RMSD_Y   | RMSD_Z   |
|       |      | (px)     | (px)     | (images) |
-------------------------------------------------
| 2     | 18   | 0.18212  | 0.28586  | 0.12145  |
| 6     | 125  | 0.15295  | 0.24972  | 0.12934  |
| 7     | 609  | 0.1286   | 0.21324  | 0.12291  |
| 8     | 265  | 0.1243   | 0.21443  | 0.10878  |
| 11    | 1189 | 0.14339  | 0.19084  | 0.10244  |
| 12    | 2389 | 0.097393 | 0.16185  | 0.1081   |
| 13    | 1550 | 0.084417 | 0.14168  | 0.098364 |
| 14    | 148  | 0.14653  | 0.17253  | 0.11116  |
| 15    | 254  | 0.21387  | 0.15463  | 0.10746  |
| 16    | 2728 | 0.11345  | 0.15427  | 0.081722 |
| 17    | 4265 | 0.060489 | 0.1104   | 0.096919 |
| 18    | 2952 | 0.07167  | 0.092076 | 0.093342 |
| 19    | 565  | 0.1626   | 0.11391  | 0.11975  |
| 20    | 534  | 0.18454  | 0.16728  | 0.16016  |
| 21    | 2772 | 0.092555 | 0.12647  | 0.076534 |
| 22    | 2873 | 0.047736 | 0.10398  | 0.090391 |
| 23    | 2686 | 0.063075 | 0.074814 | 0.075038 |
| 24    | 729  | 0.16869  | 0.11128  | 0.14306  |
| 25    | 119  | 0.81495  | 0.29087  | 0.44711  |
| 26    | 406  | 0.080441 | 0.12121  | 0.098888 |
| 27    | 425  | 0.052654 | 0.15365  | 0.0754   |
| 28    | 364  | 0.06931  | 0.066238 | 0.073655 |
| 29    | 100  | 0.85448  | 0.22713  | 0.24889  |
| 30    | 579  | 0.21735  | 0.11564  | 0.18439  |
| 31    | 1994 | 0.061756 | 0.09509  | 0.068866 |
| 32    | 1924 | 0.091952 | 0.11146  | 0.083513 |
| 33    | 2166 | 0.088252 | 0.10009  | 0.085247 |
| 34    | 519  | 0.13737  | 0.16646  | 0.1872   |
| 35    | 634  | 0.16379  | 0.093519 | 0.12164  |
| 36    | 2835 | 0.067371 | 0.080697 | 0.081958 |
| 37    | 3737 | 0.094318 | 0.086918 | 0.10408  |
| 38    | 2806 | 0.098988 | 0.13001  | 0.082836 |
| 39    | 417  | 0.17099  | 0.18814  | 0.12774  |
| 40    | 227  | 0.20834  | 0.14785  | 0.10324  |
| 41    | 1817 | 0.082715 | 0.1174   | 0.10219  |
| 42    | 2385 | 0.1096   | 0.12543  | 0.10965  |
| 43    | 1589 | 0.12167  | 0.17861  | 0.099929 |
| 44    | 88   | 0.21088  | 0.2111   | 0.10221  |
| 45    | 28   | 0.19222  | 0.16491  | 0.16288  |
| 46    | 472  | 0.094714 | 0.15589  | 0.1232   |
| 47    | 1047 | 0.11775  | 0.17961  | 0.12234  |
| 48    | 307  | 0.12735  | 0.22443  | 0.11915  |
| 51    | 22   | 0.11972  | 0.18238  | 0.11037  |
| 52    | 80   | 0.12845  | 0.24552  | 0.1396   |
-------------------------------------------------
Saving refined experiments to corrected_refined_experiments.json

As a reminder, before metrology correction we had these refined RMSDs from scan- static refinement:

RMSDs by experiment:
----------------------------------------------
| Exp | Nref  | RMSD_X  | RMSD_Y  | RMSD_Z   |
|     |       | (px)    | (px)    | (images) |
----------------------------------------------
| 0   | 53916 | 0.26341 | 0.21858 | 0.106    |
----------------------------------------------

After correction they are as follows:

RMSDs by experiment:
----------------------------------------------
| Exp | Nref  | RMSD_X  | RMSD_Y  | RMSD_Z   |
|     |       | (px)    | (px)    | (images) |
----------------------------------------------
| 0   | 53738 | 0.11368 | 0.12815 | 0.10179  |
----------------------------------------------

Let’s now do scan-varying refinement then integrate the dataset with corrected metrology:

dials.refine corrected_refined_experiments.json indexed.pickle outlier.algorithm=tukey use_all_reflections=true bin_size_fraction=0.0 scan_varying=true output.experiments=corrected_sv_refined_experiments.json
dials.integrate corrected_sv_refined_experiments.json indexed.pickle outlier.algorithm=null nproc=4 output.reflections=corrected_integrated.pickle
dials.export corrected_integrated.pickle corrected_sv_refined_experiments.json mtz.hklout=corrected_integrated.mtz ignore_panels=true
dials.analyse_output corrected_integrated.pickle grid_size=5,12

From the integration log:

Summary of integration results binned by resolution
----------------------------------------------------------------------------------------------------------
d min |  d max | # full | # part | # over | # ice | # sum | # prf | <Ibg> | <I/sigI> | <I/sigI> | <CC prf>
      |        |        |        |        |       |       |       |       |    (sum) |    (prf) |
----------------------------------------------------------------------------------------------------------
 1.17 |   1.19 |    340 |      3 |      0 |     0 |   343 |   278 |  0.04 |     0.28 |     0.41 |     0.06
 1.19 |   1.21 |   1125 |      6 |      0 |     0 |  1131 |  1011 |  0.05 |     0.37 |     0.42 |     0.06
 1.21 |   1.23 |   2346 |     13 |      0 |     0 |  2359 |  2176 |  0.05 |     0.47 |     0.49 |     0.07
 1.23 |   1.26 |   3776 |     24 |      0 |     0 |  3800 |  3565 |  0.05 |     0.50 |     0.61 |     0.10
 1.26 |   1.29 |   5448 |     37 |      0 |     0 |  5485 |  5123 |  0.05 |     0.58 |     0.71 |     0.12
 1.29 |   1.32 |   7241 |     52 |      0 |     0 |  7293 |  6804 |  0.06 |     0.66 |     0.79 |     0.15
 1.32 |   1.35 |   9522 |     71 |      0 |     0 |  9593 |  8986 |  0.06 |     0.77 |     0.94 |     0.19
 1.35 |   1.38 |  12528 |     84 |      0 |     0 | 12612 | 11851 |  0.07 |     0.92 |     1.11 |     0.22
 1.38 |   1.43 |  16954 |    113 |      0 |     0 | 17067 | 16130 |  0.08 |     1.01 |     1.22 |     0.24
 1.43 |   1.47 |  20247 |    154 |      0 |     0 | 20401 | 19499 |  0.08 |     1.22 |     1.45 |     0.29
 1.47 |   1.52 |  23553 |    267 |      0 |     0 | 23820 | 22740 |  0.09 |     1.50 |     1.74 |     0.33
 1.52 |   1.58 |  24022 |    267 |      0 |     0 | 24289 | 23215 |  0.10 |     1.78 |     2.06 |     0.38
 1.58 |   1.66 |  25479 |    286 |      0 |     0 | 25765 | 24860 |  0.11 |     2.19 |     2.49 |     0.44
 1.66 |   1.74 |  24305 |    241 |      0 |     0 | 24546 | 23395 |  0.12 |     2.74 |     3.05 |     0.51
 1.74 |   1.85 |  24742 |    236 |      0 |     0 | 24978 | 23965 |  0.15 |     3.54 |     3.87 |     0.58
 1.85 |   2.00 |  25657 |    281 |      0 |     0 | 25938 | 24999 |  0.19 |     4.90 |     5.24 |     0.65
 2.00 |   2.20 |  24755 |    268 |      0 |     0 | 25023 | 24005 |  0.26 |     6.61 |     6.97 |     0.72
 2.20 |   2.52 |  25596 |    250 |      0 |     0 | 25846 | 25058 |  0.29 |     8.88 |     9.20 |     0.76
 2.52 |   3.17 |  25140 |    298 |      0 |     0 | 25438 | 24540 |  0.37 |    12.81 |    13.04 |     0.78
 3.17 | 151.25 |  25720 |    303 |      0 |     0 | 26023 | 25193 |  0.44 |    25.52 |    25.39 |     0.79
----------------------------------------------------------------------------------------------------------

Summary of integration results for the whole dataset
----------------------------------------------
Number fully recorded                 | 338446
Number partially recorded             | 4836
Number with overloaded pixels         | 0
Number in powder rings                | 0
Number processed with summation       | 331750
Number processed with profile fitting | 317393
<Ibg>                                 | 0.18
<I/sigI> (summation)                  | 5.65
<I/sigI> (profile fitting)            | 5.94
<CC prf>                              | 0.47
----------------------------------------------

We now have a nicer plot of correlations between profiles and the strong reflections(notwithstanding a few plotting artefacts):

../../_images/reference_corr_vs_xy_multipanel_corrected.png

By comparison with the previous results we can see that correcting the panel shifts and mis-orientations has improved the overall profile fitting mean \frac{I}{\sigma_I} and correlation coefficients. However, comparing the tables of results binned by resolution we see that the improvements are limited to the low resolution, and the results are worse at high resolution. For example, we can make a quick plot of the profile fitting correlation coeffients.

../../_images/CC_metrology_comparison.png

From this we see the crossover point at about \frac{1}{d_\textrm{mid}^2} = 0.55 or a resolution of about 1.35 Angstroms.

So, finally, we conclude that although we have the means to do metrology analysis and corrections of modular detectors, we have learned that it is critical to use data that gives coverage right out to the corners of the detector so that we can apply the corrected metrology to other datasets effectively. With dials.refine this could be achieved by joint refinement of data from multiple crystals, but it will take a significant time to run and use a lot of memory. It is best to do this right, and exactly once!

What to do next

  • Panel coverage was rather low in the corners. We could try combining more metrology datasets together to improve coverage and make an effort at factoring out crystal-dependent effects.
  • Once we are satisfied with the corrected detector model we can write the changes directly back into the dxtbx Format object, so that they are available immediately for any dataset, without having to perform the ‘copy-and-paste’ operation presented here.

Acknowledgements

Dave Hall (Diamond Light Source) for collecting the data.