[503a972] | 1 | #if !defined(o_h) |
---|
| 2 | #define sphere_sld_h |
---|
| 3 | #include "parameters.hh" |
---|
| 4 | |
---|
| 5 | /** |
---|
| 6 | * Structure definition for sphere parameters |
---|
| 7 | */ |
---|
| 8 | //[PYTHONCLASS] = SphereSLDModel |
---|
| 9 | //[DISP_PARAMS] = rad_core0, thick_inter0 |
---|
| 10 | //[DESCRIPTION] =<text>Calculate neutron reflectivity using the Parratt iterative formula |
---|
| 11 | // Parameters: |
---|
| 12 | // background:background |
---|
| 13 | // scale: scale factor |
---|
| 14 | // sld_core0: the SLD of the substrate |
---|
| 15 | // sld_solv: the SLD of the incident medium |
---|
| 16 | // or superstrate |
---|
| 17 | // sld_flatN: the SLD of the flat region of |
---|
| 18 | // the N'th layer |
---|
| 19 | // thick_flatN: the thickness of the flat |
---|
| 20 | // region of the N'th layer |
---|
| 21 | // func_interN: the function used to describe |
---|
| 22 | // the interface of the N'th layer |
---|
| 23 | // nu_interN: the coefficient for the func_interN |
---|
| 24 | // thick_interN: the thickness of the interface |
---|
| 25 | // of the N'th layer |
---|
| 26 | // Note: the layer number starts to increase |
---|
| 27 | // from the bottom (substrate) to the top. |
---|
| 28 | // </text> |
---|
| 29 | //[FIXED]= <text>rad_core0.width; thick_inter0.width</text> |
---|
| 30 | //[NON_FITTABLE_PARAMS]= <text>n_shells;func_inter0;func_inter1;func_inter2;func_inter3;func_inter4;func_inter5;func_inter5;func_inter7;func_inter8;func_inter9;func_inter10 </text> |
---|
| 31 | //[ORIENTATION_PARAMS]= <text> </text> |
---|
| 32 | |
---|
| 33 | class SphereSLDModel{ |
---|
| 34 | public: |
---|
| 35 | // Model parameters |
---|
| 36 | /// number of shells |
---|
| 37 | // [DEFAULT]=n_shells=1 |
---|
| 38 | Parameter n_shells; |
---|
| 39 | /// Scale factor |
---|
| 40 | // [DEFAULT]=scale= 1.0 |
---|
| 41 | Parameter scale; |
---|
| 42 | /// thick_inter0 [A] |
---|
| 43 | // [DEFAULT]=thick_inter0=50.0 [A] |
---|
| 44 | Parameter thick_inter0; |
---|
| 45 | /// func_inter0 |
---|
| 46 | // [DEFAULT]=func_inter0= 0 |
---|
| 47 | Parameter func_inter0; |
---|
| 48 | /// sld_core0 [1/A^(2)] |
---|
| 49 | // [DEFAULT]=sld_core0= 2.07e-6 [1/A^(2)] |
---|
| 50 | Parameter sld_core0; |
---|
| 51 | /// sld_solv [1/A^(2)] |
---|
| 52 | // [DEFAULT]=sld_solv= 1.0e-6 [1/A^(2)] |
---|
| 53 | Parameter sld_solv; |
---|
| 54 | /// Background |
---|
| 55 | // [DEFAULT]=background=0 |
---|
| 56 | Parameter background; |
---|
| 57 | |
---|
| 58 | // [DEFAULT]=sld_flat1=4.0e-06 [1/A^(2)] |
---|
| 59 | Parameter sld_flat1; |
---|
| 60 | // [DEFAULT]=sld_flat2=3.5e-06 [1/A^(2)] |
---|
| 61 | Parameter sld_flat2; |
---|
| 62 | // [DEFAULT]=sld_flat3=4.0e-06 [1/A^(2)] |
---|
| 63 | Parameter sld_flat3; |
---|
| 64 | // [DEFAULT]=sld_flat4=3.5e-06 [1/A^(2)] |
---|
| 65 | Parameter sld_flat4; |
---|
| 66 | // [DEFAULT]=sld_flat5=4.0e-06 [1/A^(2)] |
---|
| 67 | Parameter sld_flat5; |
---|
| 68 | // [DEFAULT]=sld_flat6=3.5e-06 [1/A^(2)] |
---|
| 69 | Parameter sld_flat6; |
---|
| 70 | // [DEFAULT]=sld_flat7=4.0e-06 [1/A^(2)] |
---|
| 71 | Parameter sld_flat7; |
---|
| 72 | // [DEFAULT]=sld_flat8=3.5e-06 [1/A^(2)] |
---|
| 73 | Parameter sld_flat8; |
---|
| 74 | // [DEFAULT]=sld_flat9=4.0e-06 [1/A^(2)] |
---|
| 75 | Parameter sld_flat9; |
---|
| 76 | // [DEFAULT]=sld_flat10=3.5e-06 [1/A^(2)] |
---|
| 77 | Parameter sld_flat10; |
---|
| 78 | |
---|
| 79 | // [DEFAULT]=thick_inter1=50.0 [A] |
---|
| 80 | Parameter thick_inter1; |
---|
| 81 | // [DEFAULT]=thick_inter2=50.0 [A] |
---|
| 82 | Parameter thick_inter2; |
---|
| 83 | // [DEFAULT]=thick_inter3=50.0 [A] |
---|
| 84 | Parameter thick_inter3; |
---|
| 85 | // [DEFAULT]=thick_inter4=50.0 [A] |
---|
| 86 | Parameter thick_inter4; |
---|
| 87 | // [DEFAULT]=thick_inter5=50.0 [A] |
---|
| 88 | Parameter thick_inter5; |
---|
| 89 | // [DEFAULT]=thick_inter6=50.0 [A] |
---|
| 90 | Parameter thick_inter6; |
---|
| 91 | // [DEFAULT]=thick_inter7=50.0 [A] |
---|
| 92 | Parameter thick_inter7; |
---|
| 93 | // [DEFAULT]=thick_inter8=50.0 [A] |
---|
| 94 | Parameter thick_inter8; |
---|
| 95 | // [DEFAULT]=thick_inter9=50.0 [A] |
---|
| 96 | Parameter thick_inter9; |
---|
| 97 | // [DEFAULT]=thick_inter10=50.0 [A] |
---|
| 98 | Parameter thick_inter10; |
---|
| 99 | |
---|
| 100 | // [DEFAULT]=thick_flat1=100 [A] |
---|
| 101 | Parameter thick_flat1; |
---|
| 102 | // [DEFAULT]=thick_flat2=100 [A] |
---|
| 103 | Parameter thick_flat2; |
---|
| 104 | // [DEFAULT]=thick_flat3=100 [A] |
---|
| 105 | Parameter thick_flat3; |
---|
| 106 | // [DEFAULT]=thick_flat4=100 [A] |
---|
| 107 | Parameter thick_flat4; |
---|
| 108 | // [DEFAULT]=thick_flat5=100 [A] |
---|
| 109 | Parameter thick_flat5; |
---|
| 110 | // [DEFAULT]=thick_flat6=100 [A] |
---|
| 111 | Parameter thick_flat6; |
---|
| 112 | // [DEFAULT]=thick_flat7=100 [A] |
---|
| 113 | Parameter thick_flat7; |
---|
| 114 | // [DEFAULT]=thick_flat8=100 [A] |
---|
| 115 | Parameter thick_flat8; |
---|
| 116 | // [DEFAULT]=thick_flat9=100 [A] |
---|
| 117 | Parameter thick_flat9; |
---|
| 118 | // [DEFAULT]=thick_flat10=100 [A] |
---|
| 119 | Parameter thick_flat10; |
---|
| 120 | |
---|
| 121 | // [DEFAULT]=func_inter1=0 |
---|
| 122 | Parameter func_inter1; |
---|
| 123 | // [DEFAULT]=func_inter2=0 |
---|
| 124 | Parameter func_inter2; |
---|
| 125 | // [DEFAULT]=func_inter3=0 |
---|
| 126 | Parameter func_inter3; |
---|
| 127 | // [DEFAULT]=func_inter4=0 |
---|
| 128 | Parameter func_inter4; |
---|
| 129 | // [DEFAULT]=func_inter5=0 |
---|
| 130 | Parameter func_inter5; |
---|
| 131 | // [DEFAULT]=func_inter6=0 |
---|
| 132 | Parameter func_inter6; |
---|
| 133 | // [DEFAULT]=func_inter7=0 |
---|
| 134 | Parameter func_inter7; |
---|
| 135 | // [DEFAULT]=func_inter8=0 |
---|
| 136 | Parameter func_inter8; |
---|
| 137 | // [DEFAULT]=func_inter9=0 |
---|
| 138 | Parameter func_inter9; |
---|
| 139 | // [DEFAULT]=func_inter10=0 |
---|
| 140 | Parameter func_inter10; |
---|
| 141 | |
---|
| 142 | // [DEFAULT]=nu_inter1=2.5 |
---|
| 143 | Parameter nu_inter1; |
---|
| 144 | // [DEFAULT]=nu_inter2=2.5 |
---|
| 145 | Parameter nu_inter2; |
---|
| 146 | // [DEFAULT]=nu_inter3=2.5 |
---|
| 147 | Parameter nu_inter3; |
---|
| 148 | // [DEFAULT]=nu_inter4=2.5 |
---|
| 149 | Parameter nu_inter4; |
---|
| 150 | // [DEFAULT]=nu_inter5=2.5 |
---|
| 151 | Parameter nu_inter5; |
---|
| 152 | // [DEFAULT]=nu_inter6=2.5 |
---|
| 153 | Parameter nu_inter6; |
---|
| 154 | // [DEFAULT]=nu_inter7=2.5 |
---|
| 155 | Parameter nu_inter7; |
---|
| 156 | // [DEFAULT]=nu_inter8=2.5 |
---|
| 157 | Parameter nu_inter8; |
---|
| 158 | // [DEFAULT]=nu_inter9=2.5 |
---|
| 159 | Parameter nu_inter9; |
---|
| 160 | // [DEFAULT]=nu_inter10=2.5 |
---|
| 161 | Parameter nu_inter10; |
---|
| 162 | |
---|
| 163 | // [DEFAULT]=npts_inter=35.0 |
---|
| 164 | Parameter npts_inter; |
---|
| 165 | // [DEFAULT]=nu_inter0=2.5 |
---|
| 166 | Parameter nu_inter0; |
---|
| 167 | // [DEFAULT]=rad_core0=50.0 [A] |
---|
| 168 | Parameter rad_core0; |
---|
| 169 | |
---|
| 170 | // Constructor |
---|
| 171 | SphereSLDModel(); |
---|
| 172 | |
---|
| 173 | // Operators to get I(Q) |
---|
| 174 | double operator()(double q); |
---|
| 175 | double operator()(double qx, double qy); |
---|
| 176 | double calculate_ER(); |
---|
| 177 | double evaluate_rphi(double q, double phi); |
---|
| 178 | }; |
---|
| 179 | |
---|
| 180 | |
---|
| 181 | #endif |
---|