#if !defined(binaryHS_h)
#define binaryHS_h

/**
 * Structure definition for binary hard sphere parameters

	[PYTHONCLASS] = BinaryHSModel
	[DISP_PARAMS] = l_radius,s_radius
	[DESCRIPTION] =<text> Model parameters: l_radius : large radius of binary hard sphere
			s_radius : small radius of binary hard sphere
			vol_frac_ls : volume fraction of large spheres
			vol_frac_ss : volume fraction of small spheres
			ls_sld: large sphere  scattering length density
			ss_sld: small sphere scattering length density
			solvent_sld: solvent scattering length density
			background: incoherent background
               </text>
	[FIXED]= l_radius.width;s_radius.width
	[ORIENTATION_PARAMS]= <text> </text>
 */
typedef struct {

	///	large radius of the binary hard sphere [A]
    //  [DEFAULT]=l_radius= 100.0 [A]
    double l_radius;

    ///	small radius of the binary hard sphere [A]
    //  [DEFAULT]=s_radius= 25.0 [A]
    double s_radius;

	///	volume fraction of large spheres
    //  [DEFAULT]=vol_frac_ls= 0.1
    double vol_frac_ls;

	///	volume fraction of small spheres
    //  [DEFAULT]=vol_frac_ss= 0.2
    double vol_frac_ss;

	///	large sphere scattering length density [1/A^(2)]
    //  [DEFAULT]=ls_sld= 3.5e-6 [1/A^(2)]
    double ls_sld;

	///	lsmall sphere scattering length density [1/A^(2)]
    //  [DEFAULT]=ss_sld= 5e-7 [1/A^(2)]
    double ss_sld;

    ///	solvent scattering length density [1/A^(2)]
    //  [DEFAULT]=solvent_sld= 6.36e-6 [1/A^(2)]
    double solvent_sld;

	/// Incoherent Background [1/cm]
	//  [DEFAULT]=background=0.001 [1/cm]
	double background;

} BinaryHSParameters;



/// 1D scattering function
double binaryHS_analytical_1D(BinaryHSParameters *pars, double q);

/// 2D scattering function
double binaryHS_analytical_2D(BinaryHSParameters *pars, double q, double phi);
double binaryHS_analytical_2DXY(BinaryHSParameters *pars, double qx, double qy);

#endif