[ae3ce4e] | 1 | #if !defined(core_shell_cylinder_h) |
---|
| 2 | #define core_shell_cylinder_h |
---|
| 3 | |
---|
[0f5bc9f] | 4 | /** |
---|
| 5 | * Structure definition for core-shell cylinder parameters |
---|
[ae3ce4e] | 6 | */ |
---|
| 7 | //[PYTHONCLASS] = CoreShellCylinderModel |
---|
[0f5bc9f] | 8 | //[DISP_PARAMS] = radius, thickness, length, axis_theta, axis_phi |
---|
[4d3acb6] | 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> |
---|
[96672c0] | 28 | |
---|
[ae3ce4e] | 29 | typedef struct { |
---|
[0f5bc9f] | 30 | /// Scale factor |
---|
[ae3ce4e] | 31 | // [DEFAULT]=scale=1.0 |
---|
| 32 | double scale; |
---|
[0f5bc9f] | 33 | |
---|
[ae3ce4e] | 34 | /// Core radius [A] |
---|
| 35 | // [DEFAULT]=radius=20.0 A |
---|
| 36 | double radius; |
---|
[0f5bc9f] | 37 | |
---|
[ae3ce4e] | 38 | /// Shell thickness [A] |
---|
| 39 | // [DEFAULT]=thickness=10.0 A |
---|
| 40 | double thickness; |
---|
[0f5bc9f] | 41 | |
---|
[ae3ce4e] | 42 | /// Core length [A] |
---|
| 43 | // [DEFAULT]=length=400.0 A |
---|
| 44 | double length; |
---|
[0f5bc9f] | 45 | |
---|
[ae3ce4e] | 46 | /// Core SLD [A-2] |
---|
| 47 | // [DEFAULT]=core_sld=1.0e-6 A-2 |
---|
| 48 | double core_sld; |
---|
[0f5bc9f] | 49 | |
---|
[ae3ce4e] | 50 | /// Shell SLD [A-2] |
---|
| 51 | // [DEFAULT]=shell_sld=4.0e-6 A-2 |
---|
| 52 | double shell_sld; |
---|
[0f5bc9f] | 53 | |
---|
[ae3ce4e] | 54 | /// Solvent SLD [A-2] |
---|
| 55 | // [DEFAULT]=solvent_sld=1.0e-6 A-2 |
---|
| 56 | double solvent_sld; |
---|
[0f5bc9f] | 57 | |
---|
[ae3ce4e] | 58 | /// Incoherent Background [cm-1] |
---|
| 59 | // [DEFAULT]=background=0 cm-1 |
---|
[0f5bc9f] | 60 | double background; |
---|
| 61 | |
---|
[ae3ce4e] | 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; |
---|
[0f5bc9f] | 65 | |
---|
[ae3ce4e] | 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 |
---|
[0f5bc9f] | 68 | double axis_phi; |
---|
| 69 | |
---|
[ae3ce4e] | 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 |
---|