#if !defined(o_h)
#define onion_h

/**
 * Structure definition for sphere parameters
 */
 //[PYTHONCLASS] = OnionModel
 //[DISP_PARAMS] = rad_core0, thick_shell1,thick_shell2,thick_shell3,thick_shell4, thick_shell5,thick_shell6,thick_shell7,thick_shell8,thick_shell9,thick_shell10
 //[DESCRIPTION] =<text>Form factor of mutishells normalized by the volume. Here each shell is described
 //						 by an exponential function;
 //							I)
 //							For A_shell != 0,
 //							f(r) = B*exp(A_shell*(r-r_in)/thick_shell)+C
 //                         where
 //								B=(sld_out-sld_in)/(exp(A_shell)-1)
 //								C=sld_in-B.
 //							Note that in the above case,
 //								the function becomes a linear function
 //								as A_shell --> 0+ or 0-.
 //							II)
 //							For the exact point of A_shell == 0,
 //							f(r) = sld_in ,i.e., it crosses over flat function
 //							Note that the 'sld_out' becaomes NULL in this case.
 //
 //				background:background,
 //				rad_core0: radius of sphere(core)
 //				thick_shell#:the thickness of the shell#
 //				sld_core0: the SLD of the sphere
 //				sld_solv: the SLD of the solvent
 //				sld_shell: the SLD of the shell#
 //				A_shell#: the coefficient in the exponential function
 //		</text>
 //[FIXED]=  <text>rad_core0.width;thick_shell1.width;thick_shell2.width;thick_shell3.width;thick_shell4.width;thick_shell5.width;thick_shell6.width;thick_shell7.width;thick_shell8.width;thick_shell9.width;thick_shell10.width</text>
 //[ORIENTATION_PARAMS]= <text> </text>

typedef struct {
	/// number of shells
	//  [DEFAULT]=n_shells=1
	int n_shells;
    /// Scale factor
    //  [DEFAULT]=scale= 1.0
	double scale;
    ///	Radius of sphere [A]
    //  [DEFAULT]=rad_core0=200.0 [A]
	double rad_core0;
	///	sld_core [1/A^(2)]
	//  [DEFAULT]=sld_core0= 1.0e-6 [1/A^(2)]
	double sld_core0;
	///	sld_solv [1/A^(2)]
	//  [DEFAULT]=sld_solv= 6.4e-6 [1/A^(2)]
	double sld_solv;
	/// Incoherent Background [1/cm]
	//  [DEFAULT]=background=0 [1/cm]
	double background;

    //  [DEFAULT]=sld_out_shell1=2.0e-06 [1/A^(2)]
    double sld_out_shell1;
    //  [DEFAULT]=sld_out_shell2=2.5e-06 [1/A^(2)]
    double sld_out_shell2;
    //  [DEFAULT]=sld_out_shell3=3.0e-06 [1/A^(2)]
    double sld_out_shell3;
    //  [DEFAULT]=sld_out_shell4=3.5e-06 [1/A^(2)]
    double sld_out_shell4;
    //  [DEFAULT]=sld_out_shell5=4.0e-06 [1/A^(2)]
    double sld_out_shell5;
    //  [DEFAULT]=sld_out_shell6=4.5e-06 [1/A^(2)]
    double sld_out_shell6;
    //  [DEFAULT]=sld_out_shell7=5.0e-06 [1/A^(2)]
    double sld_out_shell7;
    //  [DEFAULT]=sld_out_shell8=5.5e-06 [1/A^(2)]
    double sld_out_shell8;
    //  [DEFAULT]=sld_out_shell9=6.0e-06 [1/A^(2)]
    double sld_out_shell9;
    //  [DEFAULT]=sld_out_shell10=6.2e-06 [1/A^(2)]
    double sld_out_shell10;

    //  [DEFAULT]=sld_in_shell1=1.7e-06 [1/A^(2)]
    double sld_in_shell1;
    //  [DEFAULT]=sld_in_shell2=2.2e-06 [1/A^(2)]
    double sld_in_shell2;
    //  [DEFAULT]=sld_in_shell3=2.7e-06 [1/A^(2)]
    double sld_in_shell3;
    //  [DEFAULT]=sld_in_shell4=3.2e-06 [1/A^(2)]
    double sld_in_shell4;
    //  [DEFAULT]=sld_in_shell5=3.7e-06 [1/A^(2)]
    double sld_in_shell5;
    //  [DEFAULT]=sld_in_shell6=4.2e-06 [1/A^(2)]
    double sld_in_shell6;
    //  [DEFAULT]=sld_in_shell7=4.7e-06 [1/A^(2)]
    double sld_in_shell7;
    //  [DEFAULT]=sld_in_shell8=5.2e-06 [1/A^(2)]
    double sld_in_shell8;
    //  [DEFAULT]=sld_in_shell9=5.7e-06 [1/A^(2)]
    double sld_in_shell9;
    //  [DEFAULT]=sld_in_shell10=6.0e-06 [1/A^(2)]
    double sld_in_shell10;

    //  [DEFAULT]=A_shell1=1.0
    double A_shell1;
    //  [DEFAULT]=A_shell2=1.0
    double A_shell2;
    //  [DEFAULT]=A_shell3=1.0
    double A_shell3;
    //  [DEFAULT]=A_shell4=1.0
    double A_shell4;
    //  [DEFAULT]=A_shell5=1.0
    double A_shell5;
    //  [DEFAULT]=A_shell6=1.0
    double A_shell6;
    //  [DEFAULT]=A_shell7=1.0
    double A_shell7;
    //  [DEFAULT]=A_shell8=1.0
    double A_shell8;
    //  [DEFAULT]=A_shell9=1.0
    double A_shell9;
    //  [DEFAULT]=A_shell10=1.0
    double A_shell10;

    //  [DEFAULT]=thick_shell1=50.0 [A]
    double thick_shell1;
    //  [DEFAULT]=thick_shell2=50.0 [A]
    double thick_shell2;
    //  [DEFAULT]=thick_shell3=50.0 [A]
    double thick_shell3;
    //  [DEFAULT]=thick_shell4=50.0 [A]
    double thick_shell4;
    //  [DEFAULT]=thick_shell5=50.0 [A]
    double thick_shell5;
    //  [DEFAULT]=thick_shell6=50.0 [A]
    double thick_shell6;
    //  [DEFAULT]=thick_shell7=50.0 [A]
    double thick_shell7;
    //  [DEFAULT]=thick_shell8=50.0 [A]
    double thick_shell8;
    //  [DEFAULT]=thick_shell9=50.0 [A]
    double thick_shell9;
    //  [DEFAULT]=thick_shell10=50.0 [A]
    double thick_shell10;

    //  [DEFAULT]=func_shell1=2
    int func_shell1;
    //  [DEFAULT]=func_shell2=2
    int func_shell2;
    //  [DEFAULT]=func_shell3=2
    int func_shell3;
    //  [DEFAULT]=func_shell4=2
    int func_shell4;
    //  [DEFAULT]=func_shell5=2
    int func_shell5;
    //  [DEFAULT]=func_shell6=2
    int func_shell6;
    //  [DEFAULT]=func_shell7=2
    int func_shell7;
    //  [DEFAULT]=func_shell8=2
    int func_shell8;
    //  [DEFAULT]=func_shell9=2
    int func_shell9;
    //  [DEFAULT]=func_shell10=2
    int func_shell10;
} OnionParameters;

double so_kernel(double dq[], double q);

/// 1D scattering function
double onion_analytical_1D(OnionParameters *pars, double q);

/// 2D scattering function
double onion_analytical_2D(OnionParameters *pars, double q, double phi);
double onion_analytical_2DXY(OnionParameters *pars, double qx, double qy);

#endif