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

Use dset.omega_for_bins for columnfiles rather than dset.omega #439

Merged
merged 10 commits into from
Mar 12, 2025

Conversation

jadball
Copy link
Contributor

@jadball jadball commented Mar 6, 2025

Also added a hopefully useful crystal reference frames explanation notebook that derives some functions to convert between common crystallographic frames.

Copy link

Check out this pull request on  ReviewNB

See visual diffs & provide feedback on Jupyter Notebooks.


Powered by ReviewNB

@@ -0,0 +1,1112 @@
{
Copy link
Member

@jonwright jonwright Mar 7, 2025

Choose a reason for hiding this comment

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

WIthin a cif file people will find atom_site_fract_[xyz] andatom_site_cartn_[xyz]

People use 'xyz' for these numbers (e.g. EXYZ shelx instruction).

I don't remember reading [uvw] for co-ordinates of atoms.


Reply via ReviewNB

Copy link
Contributor Author

Choose a reason for hiding this comment

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

[uvw] is commonly used to refer to a point in a lattice via the real-space basis vectors. https://www.doitpoms.ac.uk/tlplib/miller_indices/printall.php

atom_site_fract_[xyz] also use these same real-space basis vectors: https://mmcif.wwpdb.org/dictionaries/mmcif_ma.dic/Items/_atom_site.fract_x.html

atom_site_cartn_[xyz] use a cartesian frame which is our (xyz) frame in this document.

@@ -0,0 +1,1112 @@
{
Copy link
Member

@jonwright jonwright Mar 7, 2025

Choose a reason for hiding this comment

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

These are the coordinates you find in a PDB file. Also atom_site_cartn_[xyz]


Reply via ReviewNB

Copy link
Contributor Author

Choose a reason for hiding this comment

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

Yes, agreed! Will add to usage

@@ -0,0 +1,1112 @@
{
Copy link
Member

@jonwright jonwright Mar 7, 2025

Choose a reason for hiding this comment

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

Using G for this vector in documentation for ImageD11 will be confusing. There is only one kind of G vector in ImageD11/polyxsim/grainspotter and that one is the scattering vectors in sample co-ordinates.

In ImageD11 you have 'k-vectors' that are the scattering vectors measured in the laboratory frame, and then 'g' vectors that have been rotated to the frame of the sample. Both are in the frame of the sample. This one looks like the hkl space where we find out how close is an observed spot position compared to integer hkl. It is in the frame of a crystal/unit cell


Reply via ReviewNB

Copy link
Member

Choose a reason for hiding this comment

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

Or maybe not? G is in orthogonal xyz and hkl is in crystal space. So I don't get the hkl subscript. Is this G = (1,1,1) or G = (0.123,0.234,0.346) ?

Copy link
Contributor Author

Choose a reason for hiding this comment

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

This is G = (1, 1, 1) - it's in the reciprocal space basis. This agrees with the FABLE Geometry definition:

Let G be represented in the reciprocal lattice system by the integer Miller indices Ghkl = (h, k, l)t.

G in orthogonal xyz space is referred to as G_c (cartesian, w.r.t the crystal), or G_l in the orthogonal xyz lab frame.

@@ -0,0 +1,1112 @@
{
Copy link
Member

@jonwright jonwright Mar 7, 2025

Choose a reason for hiding this comment

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

F in ImageD11 is the deformation gradient tensor. It is not B transpose.


Reply via ReviewNB

Copy link
Contributor Author

Choose a reason for hiding this comment

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

The problem seems to be a lack of letters in the alphabet!

@jadball jadball changed the title Guidance on reference frame conversions Use dset.omega_for_bins for columnfiles rather than dset.omega Mar 12, 2025
@jadball
Copy link
Contributor Author

jadball commented Mar 12, 2025

@jonwright commits landed here for omega_for_bins because of the open PR. End-to-end testing seems to have worked on my side. Happy to merge?

@jonwright jonwright merged commit 9797d3f into FABLE-3DXRD:master Mar 12, 2025
6 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.

2 participants