Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Implement reindexing against a reference file (pdb, mtz, cif) in dials.reindex #2420

Merged

Conversation

jbeilstenedmands
Copy link
Contributor

Due to the series of cell transformations in cosym, when a reference dataset is given, the indexing mode can change, meaning that the final data are not reindexed into the correct mode for some use cases. This PR introduces a final test on the reindexed data, to calculate the CC with the reference dataset (in its true symmetry, without any cell reduction), and applying a final reindexing to the whole dataset to be in the most consistent indexing mode.

A few logging changes are required in dials.reindex as a side effect.

Copy link
Contributor

@rjgildea rjgildea left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Thanks for adding this, lgtm!

src/dials/command_line/reindex.py Outdated Show resolved Hide resolved
jbeilstenedmands and others added 3 commits May 31, 2023 14:36
Co-authored-by: Richard Gildea <rjgildea@users.noreply.github.com>
Add missing wavelength in cosym reference calculation, and ksol, bsol when d_min given
@jbeilstenedmands
Copy link
Contributor Author

After further thought and testing, I'm not sure if what I was trying to do is particularly robust. As the data are unscaled, it does not become reliable to calculate a cc against the reference (I was trying to do this on the expt and refl file data to avoid the symmetry transformations, but then the intensities haven't even had the normalisation applied). I think we can achieve what we want by a simple call to dials.reindex using a reference after scaling, so I think I will back out some of the changes to cosym and keep the recent commits to add the functionality to dials.reindex

@jbeilstenedmands jbeilstenedmands changed the title Add extra final reindexing test when using dials.cosym with a reference dataset. Implement reindexing against a reference file (pdb, mtz, cif) in dials.reindex Jun 1, 2023
@jbeilstenedmands jbeilstenedmands merged commit ea9b6ea into dials:main Jun 1, 2023
10 checks passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

3 participants