1 | #if !defined(core_shell_cylinder_h) |
---|
2 | #define core_shell_cylinder_h |
---|
3 | |
---|
4 | /** |
---|
5 | * Structure definition for core-shell cylinder parameters |
---|
6 | */ |
---|
7 | //[PYTHONCLASS] = CoreShellCylinderModel |
---|
8 | //[DISP_PARAMS] = radius, thickness, length, axis_theta, axis_phi |
---|
9 | //[DESCRIPTION] = <text>P(q,alpha)= scale/Vs*f(q)^(2) + bkg Where:\n\ |
---|
10 | // f(q)= 2(core_sld- solvant_sld)* Vc*sin[qLcos(alpha/2)]/\n\ |
---|
11 | // [qLcos(alpha/2)]*J1(qRsin(alpha))/[qRsin(alpha)] +\n 2(shell_sld-solvent_sld)*Vs |
---|
12 | // *sin[q(L+T)cos(alpha/2)]/[[q(L+T)cos(alpha/2)] |
---|
13 | // *J1(q(R+T)sin(alpha))/q(R+T)sin(alpha)] |
---|
14 | // alpha:is the angle between the axis of the cylinder and the q-vector |
---|
15 | // Vs: the volume of the outer shell |
---|
16 | // Vc: the volume of the core |
---|
17 | // L: the length of the core |
---|
18 | // shell_sld: the scattering length density of the shell |
---|
19 | // solvent_sld: the scattering length density of the solvent |
---|
20 | // bkg: the background |
---|
21 | // T: the thickness |
---|
22 | // R+T: is the outer radius |
---|
23 | // L+2T: The total length of the outershell |
---|
24 | // J1: the first order Bessel function |
---|
25 | // theta: axis_theta of the cylinder |
---|
26 | // phi: the axis_phi of the cylinder |
---|
27 | // </text> |
---|
28 | |
---|
29 | typedef struct { |
---|
30 | /// Scale factor |
---|
31 | // [DEFAULT]=scale=1.0 |
---|
32 | double scale; |
---|
33 | |
---|
34 | /// Core radius [A] |
---|
35 | // [DEFAULT]=radius=20.0 A |
---|
36 | double radius; |
---|
37 | |
---|
38 | /// Shell thickness [A] |
---|
39 | // [DEFAULT]=thickness=10.0 A |
---|
40 | double thickness; |
---|
41 | |
---|
42 | /// Core length [A] |
---|
43 | // [DEFAULT]=length=400.0 A |
---|
44 | double length; |
---|
45 | |
---|
46 | /// Core SLD [A-2] |
---|
47 | // [DEFAULT]=core_sld=1.0e-6 A-2 |
---|
48 | double core_sld; |
---|
49 | |
---|
50 | /// Shell SLD [A-2] |
---|
51 | // [DEFAULT]=shell_sld=4.0e-6 A-2 |
---|
52 | double shell_sld; |
---|
53 | |
---|
54 | /// Solvent SLD [A-2] |
---|
55 | // [DEFAULT]=solvent_sld=1.0e-6 A-2 |
---|
56 | double solvent_sld; |
---|
57 | |
---|
58 | /// Incoherent Background [cm-1] |
---|
59 | // [DEFAULT]=background=0 cm-1 |
---|
60 | double background; |
---|
61 | |
---|
62 | /// Orientation of the long axis of the core-shell cylinder w/respect incoming beam [rad] |
---|
63 | // [DEFAULT]=axis_theta=1.57 rad |
---|
64 | double axis_theta; |
---|
65 | |
---|
66 | /// Orientation of the long axis of the core-shell cylinder in the plane of the detector [rad] |
---|
67 | // [DEFAULT]=axis_phi=0.0 rad |
---|
68 | double axis_phi; |
---|
69 | |
---|
70 | } CoreShellCylinderParameters; |
---|
71 | |
---|
72 | |
---|
73 | |
---|
74 | /// 1D scattering function |
---|
75 | double core_shell_cylinder_analytical_1D(CoreShellCylinderParameters *pars, double q); |
---|
76 | |
---|
77 | /// 2D scattering function |
---|
78 | double core_shell_cylinder_analytical_2D(CoreShellCylinderParameters *pars, double q, double phi); |
---|
79 | double core_shell_cylinder_analytical_2DXY(CoreShellCylinderParameters *pars, double qx, double qy); |
---|
80 | double core_shell_cylinder_analytical_2D_scaled(CoreShellCylinderParameters *pars, double q, double q_x, double q_y); |
---|
81 | |
---|
82 | #endif |
---|