|
24 | 24 | import MDAnalysis as mda |
25 | 25 | from MDAnalysis.analysis.rdf import InterRDF |
26 | 26 |
|
27 | | -import maicos |
28 | | -from maicos.lib.math import atomic_form_factor, rdf_structure_factor |
| 27 | +from scatterkit.lib.math import atomic_form_factor, rdf_structure_factor |
| 28 | +import scatterkit |
29 | 29 |
|
30 | 30 | # %% |
31 | 31 | # The `water` system consists of 510 water molecules in the liquid state. The |
|
47 | 47 | # Extract small angle x-ray scattering (SAXS) intensities |
48 | 48 | # ------------------------------------------------------- |
49 | 49 | # |
50 | | -# Let us use the :class:`maicos.Saxs` class of MAICoS and apply it to all atoms in the |
| 50 | +# Let us use the :class:`scatterkit.Saxs` class of MAICoS and apply it to all atoms in the |
51 | 51 | # system: |
52 | 52 |
|
53 | | -saxs = maicos.Saxs(u.atoms).run(stop=30) |
| 53 | +saxs = scatterkit.Saxs(u.atoms).run(stop=30) |
54 | 54 |
|
55 | 55 | # %% |
56 | 56 | # .. Note:: |
|
109 | 109 | # contributions individually. Let us calculate both oxygen and hydrogen contributions, |
110 | 110 | # respectively: |
111 | 111 |
|
112 | | -saxs_O = maicos.Saxs(group_O).run(stop=30) |
113 | | -saxs_H = maicos.Saxs(group_H).run(stop=30) |
| 112 | +saxs_O = scatterkit.Saxs(group_O).run(stop=30) |
| 113 | +saxs_H = scatterkit.Saxs(group_H).run(stop=30) |
114 | 114 |
|
115 | 115 | # %% |
116 | 116 | # Let us plot the results for the structure factor, the squared atomic form factor as |
117 | 117 | # well scattering intensities together. For computing the atomic form factor we will use |
118 | | -# :func:`maicos.lib.math.atomic_form_factor`. Note that for the ``structure_factors`` |
| 118 | +# :func:`scatterkit.lib.math.atomic_form_factor`. Note that for the ``structure_factors`` |
119 | 119 | # and the ``scattering_intensities` we access the results directly from the ``results`` |
120 | 120 | # attribute without storing them in individual variables as before: |
121 | 121 |
|
|
192 | 192 | # %% |
193 | 193 | # We use ``exclude_same="residue"`` to exclude atomic self contributions resulting in a |
194 | 194 | # large peak at 0. Next, we convert the RDF into a structure factor using |
195 | | -# :func:`maicos.lib.math.rdf_structure_factor` and the number density of the |
| 195 | +# :func:`scatterkit.lib.math.rdf_structure_factor` and the number density of the |
196 | 196 | # oxygens. |
197 | 197 |
|
198 | 198 | density = group_O.n_atoms / u.trajectory.ts.volume |
|
0 commit comments