Changeset eb69cce in sasmodels for sasmodels/models/cylinder.py
- Timestamp:
- Nov 30, 2015 7:18:41 PM (8 years ago)
- Branches:
- master, core_shell_microgels, costrafo411, magnetic_model, release_v0.94, release_v0.95, ticket-1257-vesicle-product, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests
- Children:
- d18f8a8
- Parents:
- d138d43
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
sasmodels/models/cylinder.py
rd138d43 reb69cce 12 12 .. math:: 13 13 14 P( Q,\alpha) = {\text{scale} \over V} F^2(Q) + \text{background}14 P(q,\alpha) = \frac{\text{scale}}{V} F^2(q) + \text{background} 15 15 16 16 where … … 18 18 .. math:: 19 19 20 F( Q) = 2 (\Delta \rho) V21 {\sin \left(Q\tfrac12 L\cos\alpha \right)22 \over Q\tfrac12 L \cos \alpha}23 {J_1 \left(Q R \sin \alpha\right) \over QR \sin \alpha}20 F(q) = 2 (\Delta \rho) V 21 \frac{\sin \left(q\tfrac12 L\cos\alpha \right)} 22 {q\tfrac12 L \cos \alpha} 23 \frac{J_1 \left(q R \sin \alpha\right)}{q R \sin \alpha} 24 24 25 25 and $\alpha$ is the angle between the axis of the cylinder and $\vec q$, $V$ … … 44 44 45 45 NB: The 2nd virial coefficient of the cylinder is calculated based on the 46 radius and length values, and used as the effective radius for $S( Q)$47 when $P( Q) \cdot S(Q)$ is applied.46 radius and length values, and used as the effective radius for $S(q)$ 47 when $P(q) \cdot S(q)$ is applied. 48 48 49 49 The output of the 1D scattering intensity function for randomly oriented … … 52 52 .. math:: 53 53 54 P( Q) = {\text{scale} \overV}55 \int_0^{\pi/2} F^2( Q,\alpha) \sin \alpha\ d\alpha + \text{background}54 P(q) = \frac{\text{scale}}{V} 55 \int_0^{\pi/2} F^2(q,\alpha) \sin \alpha\ d\alpha + \text{background} 56 56 57 The *theta* and *phi* parameters are not used for the 1D output. Our 58 implementation of the scattering kernel and the 1D scattering intensity 59 use the c-library from NIST. 57 The $\theta$ and $\phi$ parameters are not used for the 1D output. 60 58 61 59 Validation … … 82 80 .. math:: 83 81 84 P( Q) = \int_0^{\pi/2} d\phi85 \int_0^\pi p(\theta, \phi) P_0( Q,\alpha) \sin \theta\ d\theta82 P(q) = \int_0^{\pi/2} d\phi 83 \int_0^\pi p(\theta, \phi) P_0(q,\alpha) \sin \theta\ d\theta 86 84 87 85 88 86 where $p(\theta,\phi)$ is the probability distribution for the orientation 89 and $P_0( Q,\alpha)$ is the scattering intensity for the fully oriented87 and $P_0(q,\alpha)$ is the scattering intensity for the fully oriented 90 88 system. Since we have no other software to compare the implementation of 91 89 the intensity for fully oriented cylinders, we can compare the result of … … 129 127 130 128 # [ "name", "units", default, [lower, upper], "type", "description"], 131 parameters = [["sld", " 1e-6/Ang^2", 4, [-inf, inf], "",129 parameters = [["sld", "4e-6/Ang^2", 4, [-inf, inf], "", 132 130 "Cylinder scattering length density"], 133 131 ["solvent_sld", "1e-6/Ang^2", 1, [-inf, inf], "",
Note: See TracChangeset
for help on using the changeset viewer.