/*
TODO: Add 2D model
*/
#if !defined(stacked_disks_h)
#define stacked_disks_h
/** Structure definition for stacked disks parameters
* [PYTHONCLASS] = StackedDisksModel
* [DISP_PARAMS] = length, radius, axis_theta, axis_phi
[DESCRIPTION] =
[FIXED]= length.width; radius.width; axis_theta.width; axis_phi.width
[ORIENTATION_PARAMS]= axis_phi; axis_theta; axis_phi.width; axis_theta.width
**/
typedef struct {
/// Scale factor
// [DEFAULT]=scale=0.01
double scale;
/// radius of the staked disk [A]
// [DEFAULT]=radius=3000 [A]
double radius;
/// Length of the staked disk [A]
// [DEFAULT]=length=10.0 [A]
double length;
/// Thickness of the staked disk [A]
// [DEFAULT]=thickness=15.0 [A]
double thickness;
/// Core scattering length density[1/AČ]
// [DEFAULT]=core_sld=4e-6 [1/AČ]
double core_sld;
/// layer scattering length density[1/AČ]
// [DEFAULT]=layer_sld=-4e-7 [1/AČ]
double layer_sld;
/// solvent scattering length density[1/AČ]
// [DEFAULT]=solvent_sld=5.0e-6 [1/AČ]
double solvent_sld;
/// number of layers
// [DEFAULT]=nlayers=1
double nlayers;
/// GSD of disks spacing
// [DEFAULT]=spacing=0
double spacing;
/// Incoherent Background [1/cm]
// [DEFAULT]=background=0.001 [1/cm]
double background;
/// Orientation of the staked disk axis w/respect incoming beam [rad]
// [DEFAULT]=axis_theta=1.0 [rad]
double axis_theta;
/// Orientation of the staked disk in the plane of the detector [rad]
// [DEFAULT]=axis_phi=1.0 [rad]
double axis_phi;
} StackedDisksParameters;
/// 1D scattering function
double stacked_disks_analytical_1D(StackedDisksParameters *pars, double q);
/// 2D scattering function
double stacked_disks_analytical_2D(StackedDisksParameters *pars, double q, double phi);
double stacked_disks_analytical_2DXY(StackedDisksParameters *pars, double qx, double qy);
double stacked_disks_analytical_2D_scaled(StackedDisksParameters *pars, double q, double q_x, double q_y);
#endif