/*
	TODO: Add 2D model
*/

#if !defined(stacked_disks_h)
#define stacked_disks_h
/** Structure definition for stacked disks parameters
 * [PYTHONCLASS] = StackedDisksModel
 * [DISP_PARAMS] = core_thick, layer_thick, radius, axis_theta, axis_phi
   [DESCRIPTION] = <text> One layer of disk consists of a core, a top layer, and a bottom layer.
		radius =  the radius of the disk
		core_thick = thickness of the core
		layer_thick = thickness of a layer
		core_sld = the SLD of the core
		layer_sld = the SLD of the layers
		n_stacking = the number of the disks
		sigma_d =  Gaussian STD of d-spacing
		solvent_sld = the SLD of the solvent
		</text>
	[FIXED]= <text>core_thick.width;layer_thick.width; radius.width; axis_theta.width; axis_phi.width</text>
	[ORIENTATION_PARAMS]= <text>axis_phi; axis_theta; axis_phi.width; axis_theta.width</text>


 **/
typedef struct {
    /// Scale factor
    //  [DEFAULT]=scale=0.01
    double scale;
	/// radius of the staked disk [A]
    //  [DEFAULT]=radius=3000 [A]
    double radius;
    /// Thickness of the core disk [A]
    //  [DEFAULT]=core_thick=10.0 [A]
    double core_thick;
	/// Thickness of the staked disk [A]
    //  [DEFAULT]=layer_thick=15.0 [A]
    double layer_thick;
	/// Core scattering length density[1/A^(2)]
    //  [DEFAULT]=core_sld=4e-6 [1/A^(2)]
    double core_sld;
	/// layer scattering length density[1/A^(2)]
    //  [DEFAULT]=layer_sld=-4e-7 [1/A^(2)]
    double layer_sld;
	/// solvent scattering length density[1/A^(2)]
    //  [DEFAULT]=solvent_sld=5.0e-6 [1/A^(2)]
    double solvent_sld;
    /// number of stacking
    //  [DEFAULT]=n_stacking=1
	double n_stacking;
    /// GSD of disks sigma_d
    //  [DEFAULT]=sigma_d=0
    double sigma_d;
	/// 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=0.0 [rad]
    double axis_theta;
    /// Orientation of the  staked disk in the plane of the detector [rad]
    //  [DEFAULT]=axis_phi=0.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