source: sasview/sansmodels/src/sans/models/c_extensions/ellipsoid.h @ d6da3b1

ESS_GUIESS_GUI_DocsESS_GUI_batch_fittingESS_GUI_bumps_abstractionESS_GUI_iss1116ESS_GUI_iss879ESS_GUI_iss959ESS_GUI_openclESS_GUI_orderingESS_GUI_sync_sascalccostrafo411magnetic_scattrelease-4.1.1release-4.1.2release-4.2.2release_4.0.1ticket-1009ticket-1094-headlessticket-1242-2d-resolutionticket-1243ticket-1249ticket885unittest-saveload
Last change on this file since d6da3b1 was 4628e31, checked in by Jae Cho <jhjcho@…>, 14 years ago

changed the unit of angles into degrees

  • Property mode set to 100644
File size: 2.4 KB
Line 
1#if !defined(ellipsoid_h)
2#define ellipsoid_h
3
4/**
5 * Structure definition for ellipsoid parameters
6 * The ellipsoid has axes radius_b, radius_b, radius_a.
7 * Ref: Jan Skov Pedersen, Advances in Colloid and Interface Science, 70 (1997) 171-210
8 */
9 //[PYTHONCLASS] = EllipsoidModel
10 //[DISP_PARAMS] = radius_a, radius_b, axis_theta, axis_phi
11 //[DESCRIPTION] = <text>"P(q.alpha)= scale*f(q)^(2)+ bkg, where f(q)= 3*(sld_ell
12 //             - sld_solvent)*V*[sin(q*r(Ra,Rb,alpha))
13 //             -q*r*cos(qr(Ra,Rb,alpha))]
14 //             /[qr(Ra,Rb,alpha)]^(3)"
15 //
16 //     r(Ra,Rb,alpha)= [Rb^(2)*(sin(alpha))^(2)
17 //     + Ra^(2)*(cos(alpha))^(2)]^(1/2)
18 //
19 //             scatter_sld: SLD of the scatter
20 //             solvent_sld: SLD of the solvent
21 //     sldEll: SLD of ellipsoid
22 //             sldSolv: SLD of solvent
23 //             V: volune of the Eliipsoid
24 //             Ra: radius along the rotation axis
25 //             of the Ellipsoid
26 //             Rb: radius perpendicular to the
27 //             rotation axis of the ellipsoid
28 //             </text>
29 //[FIXED]= <text> axis_phi.width; axis_theta.width;radius_a.width;
30 //radius_b.width; length.width; r_minor.width;
31 //r_ratio.width</text>
32 //[ORIENTATION_PARAMS]=  axis_phi.width; axis_theta.width;axis_phi; axis_theta
33
34
35typedef struct {
36    /// Scale factor
37    //  [DEFAULT]=scale=1.0
38    double scale;
39
40    /// Rotation axis radius_a [A]
41    //  [DEFAULT]=radius_a=20.0 [A]
42    double radius_a;
43
44    /// Radius_b [A]
45    //  [DEFAULT]=radius_b=400 [A]
46    double radius_b;
47
48    /// sldEll [1/A^(2)]
49    //  [DEFAULT]=sldEll=4.0e-6 [1/A^(2)]
50    double sldEll;
51
52    /// sld of solvent [1/A^(2)]
53    //  [DEFAULT]=sldSolv=1.0e-6 [1/A^(2)]
54    double sldSolv;
55
56        /// Incoherent Background [1/cm]
57        //  [DEFAULT]=background=0 [1/cm]
58        double background;
59
60    /// Orientation of the long axis of the ellipsoid w/respect incoming beam [deg]
61    //  [DEFAULT]=axis_theta=90.0 [deg]
62    double axis_theta;
63    /// Orientation of the long axis of the ellipsoid in the plane of the detector [deg]
64    //  [DEFAULT]=axis_phi=0.0 [deg]
65    double axis_phi;
66} EllipsoidParameters;
67
68
69
70/// 1D scattering function
71double ellipsoid_analytical_1D(EllipsoidParameters *pars, double q);
72
73/// 2D scattering function
74double ellipsoid_analytical_2D(EllipsoidParameters *pars, double q, double phi);
75double ellipsoid_analytical_2DXY(EllipsoidParameters *pars, double qx, double qy);
76double ellipsoid_analytical_2D_scaled(EllipsoidParameters *pars, double q, double q_x, double q_y);
77
78#endif
Note: See TracBrowser for help on using the repository browser.