Changeset 17574ae in sasview for sansmodels/src/sans/models/media/smear_computation.html
- Timestamp:
- Aug 14, 2012 6:00:12 PM (13 years ago)
- Branches:
- master, ESS_GUI, ESS_GUI_Docs, ESS_GUI_batch_fitting, ESS_GUI_bumps_abstraction, ESS_GUI_iss1116, ESS_GUI_iss879, ESS_GUI_iss959, ESS_GUI_opencl, ESS_GUI_ordering, ESS_GUI_sync_sascalc, costrafo411, magnetic_scatt, release-4.1.1, release-4.1.2, release-4.2.2, release_4.0.1, ticket-1009, ticket-1094-headless, ticket-1242-2d-resolution, ticket-1243, ticket-1249, ticket885, unittest-saveload
- Children:
- 5a313a3
- Parents:
- 96277b9
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
sansmodels/src/sans/models/media/smear_computation.html
r76b5220 r17574ae 4 4 <meta http-equiv=Content-Type content="text/html; charset=windows-1252"> 5 5 <meta name=Generator content="Microsoft Word 12 (filtered)"> 6 <style>7 <!--8 /* Font Definitions */9 @font-face10 {font-family:Wingdings;11 panose-1:5 0 0 0 0 0 0 0 0 0;}12 @font-face13 {font-family:"Cambria Math";14 panose-1:2 4 5 3 5 4 6 3 2 4;}15 @font-face16 {font-family:Calibri;17 panose-1:2 15 5 2 2 2 4 3 2 4;}18 @font-face19 {font-family:Tahoma;20 panose-1:2 11 6 4 3 5 4 4 2 4;}21 /* Style Definitions */22 p.MsoNormal, li.MsoNormal, div.MsoNormal23 {margin-top:0in;24 margin-right:0in;25 margin-bottom:10.0pt;26 margin-left:0in;27 line-height:115%;28 font-size:11.0pt;29 font-family:"Calibri","sans-serif";}30 p.MsoAcetate, li.MsoAcetate, div.MsoAcetate31 {mso-style-link:"Balloon Text Char";32 margin:0in;33 margin-bottom:.0001pt;34 font-size:8.0pt;35 font-family:"Tahoma","sans-serif";}36 p.MsoListParagraph, li.MsoListParagraph, div.MsoListParagraph37 {margin-top:0in;38 margin-right:0in;39 margin-bottom:10.0pt;40 margin-left:.5in;41 line-height:115%;42 font-size:11.0pt;43 font-family:"Calibri","sans-serif";}44 p.MsoListParagraphCxSpFirst, li.MsoListParagraphCxSpFirst, div.MsoListParagraphCxSpFirst45 {margin-top:0in;46 margin-right:0in;47 margin-bottom:0in;48 margin-left:.5in;49 margin-bottom:.0001pt;50 line-height:115%;51 font-size:11.0pt;52 font-family:"Calibri","sans-serif";}53 p.MsoListParagraphCxSpMiddle, li.MsoListParagraphCxSpMiddle, div.MsoListParagraphCxSpMiddle54 {margin-top:0in;55 margin-right:0in;56 margin-bottom:0in;57 margin-left:.5in;58 margin-bottom:.0001pt;59 line-height:115%;60 font-size:11.0pt;61 font-family:"Calibri","sans-serif";}62 p.MsoListParagraphCxSpLast, li.MsoListParagraphCxSpLast, div.MsoListParagraphCxSpLast63 {margin-top:0in;64 margin-right:0in;65 margin-bottom:10.0pt;66 margin-left:.5in;67 line-height:115%;68 font-size:11.0pt;69 font-family:"Calibri","sans-serif";}70 span.BalloonTextChar71 {mso-style-name:"Balloon Text Char";72 mso-style-link:"Balloon Text";73 font-family:"Tahoma","sans-serif";}74 .MsoPapDefault75 {margin-bottom:10.0pt;76 line-height:115%;}77 @page WordSection178 {size:8.5in 11.0in;79 margin:1.0in 1.0in 1.0in 1.0in;}80 div.WordSection181 {page:WordSection1;}82 /* List Definitions */83 ol84 {margin-bottom:0in;}85 ul86 {margin-bottom:0in;}87 -->88 </style>89 6 90 7 </head> … … 96 13 <p class=MsoNormal><span style='font-size:16.0pt;line-height:115%;font-family: 97 14 "Times New Roman","serif"'><h4>Smear Computation </h4></span></p> 98 <p class=MsoNormal> </p> 15 99 16 100 17 <ul style='margin-top:0in' type=disc> … … 106 23 </li> 107 24 </ul> 108 <p class=MsoNormal> </p> 109 <p class=MsoNormal> </p> 25 110 26 <p class=MsoListParagraph><span style='font-size:14.0pt;line-height:115%; 111 27 font-family:"Times New Roman","serif"'><h5><a name="Slit Smear">Slit Smear</a></h5></span></p> … … 115 31 116 32 <p class=MsoNormal><img width=349 height=49 117 src=" sm_image002.gif" align=left hspace=12></p>118 119 <p class=MsoNormal><span style='font-family:"Times New Roman","serif"'> 120 33 src="./img/sm_image002.gif" align=left hspace=12></p> 34 35 <p class=MsoNormal><span style='font-family:"Times New Roman","serif"'> 36 1)</span><br clear=all> 121 37 <span style='font-family:"Times New Roman","serif"'>where Norm = <span 122 38 style='position:relative;top:15.0pt'><img width=137 height=49 123 src=" sm_image003.gif"></span>.</span></p>39 src="./img/sm_image003.gif"></span>.</span></p> 124 40 125 41 <p class=MsoNormal><span style='font-family:"Times New Roman","serif"'>The 126 42 functions <span style='position:relative;top:6.0pt'><img width=43 height=25 127 src=" sm_image004.gif"></span>and <span style='position:43 src="./img/sm_image004.gif"></span>and <span style='position: 128 44 relative;top:6.0pt'><img width=43 height=25 129 src=" sm_image005.gif"></span>refer to the slit width weighting45 src="./img/sm_image005.gif"></span>refer to the slit width weighting 130 46 function and the slit height weighting determined at the q point, respectively. 131 47 Here, we assumes that the weighting function is described by a rectangular 132 48 function, i.e.,</span></p> 133 49 134 50 <p class=MsoNormal><span style='position:relative;top:7.0pt'><img width=134 135 height=26 src=" sm_image006.gif">136 51 height=26 src="./img/sm_image006.gif"> 52 </span><span style='font-family:"Times New Roman","serif";position:relative; 137 53 top:7.0pt'>2)</span></p> 138 54 … … 140 56 141 57 <p class=MsoNormal><span style='position:relative;top:7.0pt'><img width=136 142 height=26 src=" sm_image007.gif"></span>,143 58 height=26 src="./img/sm_image007.gif"></span>, 59 <span style='font-family:"Times New Roman","serif"'>3)</span></p> 144 60 145 61 <p class=MsoNormal><span style='font-family:"Times New Roman","serif"'>so that </span><span 146 62 style='position:relative;top:6.0pt'><img width=58 height=23 147 src=" sm_image008.gif"></span><span style='position:relative;148 top:16.0pt'><img width=76 height=51 src=" sm_image009.gif"></span><span149 style='font-family:"Times New Roman","serif"'>for</span> 63 src="./img/sm_image008.gif"></span> <span style='position:relative; 64 top:16.0pt'><img width=76 height=51 src="./img/sm_image009.gif"></span> <span 65 style='font-family:"Times New Roman","serif"'>for</span> <span 150 66 style='position:relative;top:3.0pt'><img width=40 height=15 151 src=" sm_image010.gif"></span><span style='font-family:67 src="./img/sm_image010.gif"></span> <span style='font-family: 152 68 "Times New Roman","serif"'>and <i>u</i>. The </span><span style='position:relative; 153 top:6.0pt'><img width=28 height=24 src=" sm_image011.gif"></span><span69 top:6.0pt'><img width=28 height=24 src="./img/sm_image011.gif"></span> <span 154 70 style='font-family:"Times New Roman","serif"'>and </span><span 155 71 style='position:relative;top:6.0pt'><img width=28 height=24 156 src=" sm_image012.gif"> </span><span style='font-family:72 src="./img/sm_image012.gif"> </span><span style='font-family: 157 73 "Times New Roman","serif"'>stand for the slit height (FWHM/2) and the slit 158 74 width (FWHM/2) in the q space. Now the integral of Eq. (1) is simplified to</span></p> 159 75 160 76 <p class=MsoNormal><img width=283 height=52 161 src=" sm_image013.gif" align=left hspace=12><span162 style='font-family:"Times New Roman","serif"'> 163 77 src="./img/sm_image013.gif" align=left hspace=12><span 78 style='font-family:"Times New Roman","serif"'> 79 4)</span></p> 164 80 165 81 <p class=MsoNormal><span style='font-family:"Times New Roman","serif"; … … 174 90 </span></span><span style='font-family:"Times New Roman","serif"'>For </span><span 175 91 style='position:relative;top:6.0pt'><img width=28 height=24 176 src=" sm_image014.gif"></span>= 0<span style='font-family:92 src="./img/sm_image014.gif"></span>= 0 <span style='font-family: 177 93 "Times New Roman","serif"'>and </span><span style='position:relative; 178 top:6.0pt'><img width=28 height=24 src=" sm_image015.gif"></span>=94 top:6.0pt'><img width=28 height=24 src="./img/sm_image015.gif"></span> = 179 95 <span style='font-family:"Times New Roman","serif"'>constant:</span></p> 180 96 181 97 <p class=MsoListParagraphCxSpMiddle style='margin-left:.25in'> 182 <img 183 src="sm_image016.gif"></p> 98 <img src="./img/sm_image016.gif"></p> 184 99 185 100 <p class=MsoListParagraphCxSpMiddle style='margin-left:.25in'><span 186 101 style='font-family:"Times New Roman","serif"'>For discrete q values, at the q 187 values from the data points and at the q values extended up to 102 values from the data points and at the q values extended up to q<sub>N</sub>= 188 103 q<sub>i</sub> + </span><span style='position:relative;top:6.0pt'><img width=28 189 height=24 src=" sm_image011.gif"></span><span104 height=24 src="./img/sm_image011.gif"></span><span 190 105 style='font-family:"Times New Roman","serif"'>, the smeared intensity can be 191 106 calculated approximately,</span></p> 192 107 193 108 <p class=MsoListParagraphCxSpMiddle style='margin-left:.25in'><img 194 src=" sm_image017.gif">. <span195 style='font-family:"Times New Roman","serif"'>5)</span></p>109 src="./img/sm_image017.gif">. 110 <span style='font-family:"Times New Roman","serif"'>5)</span></p> 196 111 197 112 <p class=MsoListParagraphCxSpMiddle style='margin-left:.25in'><span 198 113 style='position:relative;top:7.0pt'><img width=23 height=25 199 src=" sm_image018.gif"></span><span style='font-family:114 src="./img/sm_image018.gif"></span> <span style='font-family: 200 115 "Times New Roman","serif"'>= 0 for <i>I<sub>s</sub></i> in</span> <i><span 201 116 style='font-family:"Times New Roman","serif"'>j < i</span></i><span … … 207 122 <p class=MsoListParagraphCxSpMiddle style='margin-left:.25in;text-indent:-.25in'><span 208 123 style='font-family:"Times New Roman","serif"'>2)<span style='font:7.0pt "Times New Roman"'> 209 </span></span><span style='font-family:"Times New Roman","serif"'>For 124 </span></span><span style='font-family:"Times New Roman","serif"'>For </span><span 210 125 style='position:relative;top:6.0pt'><img width=28 height=24 211 src=" sm_image014.gif"></span>= <span style='font-family:212 "Times New Roman","serif"'>constant </span> 126 src="./img/sm_image014.gif"></span>= <span style='font-family: 127 "Times New Roman","serif"'>constant </span> <span style='font-family:"Times New Roman","serif"'>and 213 128 </span><span style='position:relative;top:6.0pt'><img width=28 height=24 214 src=" sm_image015.gif"></span>= <span style='font-family:129 src="./img/sm_image015.gif"></span> = <span style='font-family: 215 130 "Times New Roman","serif"'>0:</span></p> 216 131 … … 219 134 220 135 <p class=MsoListParagraphCxSpMiddle style='margin-left:.25in'> 221 <img 222 src="sm_image019.gif"> 223 <span 224 style='font-family:"Times New Roman","serif"'>6)</span></p> 225 226 <p class=MsoListParagraphCxSpMiddle style='margin-left:.25in'><span 227 style='font-family:"Times New Roman","serif"'>for q<sub>p</sub> = q<sub>i</sub> 136 <img src="./img/sm_image019.gif"> 137 <span style='font-family:"Times New Roman","serif"'>6)</span></p> 138 139 <p class=MsoListParagraphCxSpMiddle style='margin-left:.25in'><span 140 style='font-family:"Times New Roman","serif"'>for q<sub>p</sub> = q<sub>i</sub> 228 141 - </span><span style='position:relative;top:6.0pt'><img width=28 height=24 229 src=" sm_image012.gif"></span><span style='font-family:230 "Times New Roman","serif"'> 142 src="./img/sm_image012.gif"></span><span style='font-family: 143 "Times New Roman","serif"'> and</span> <span style='font-family:"Times New Roman","serif"'>q<sub>N</sub> 231 144 = q<sub>i</sub> + </span><span style='position:relative;top:6.0pt'><img 232 width=28 height=24 src=" sm_image012.gif"></span>.<span145 width=28 height=24 src="./img/sm_image012.gif"></span>. <span 233 146 style='position:relative;top:7.0pt'><img width=23 height=25 234 src=" sm_image018.gif"></span><span style='font-family:147 src="./img/sm_image018.gif"></span> <span style='font-family: 235 148 "Times New Roman","serif"'>= 0 for <i>I<sub>s</sub></i> in</span> <i><span 236 149 style='font-family:"Times New Roman","serif"'>j < p</span></i><span … … 241 154 <p class=MsoListParagraphCxSpMiddle style='margin-left:.25in;text-indent:-.25in'><span 242 155 style='font-family:"Times New Roman","serif"'>3)<span style='font:7.0pt "Times New Roman"'> 243 </span></span><span style='font-family:"Times New Roman","serif"'>For 156 </span></span><span style='font-family:"Times New Roman","serif"'>For </span><span 244 157 style='position:relative;top:6.0pt'><img width=28 height=24 245 src=" sm_image014.gif"></span>= <span style='font-family:246 "Times New Roman","serif"'>constant </span> 158 src="./img/sm_image014.gif"></span>= <span style='font-family: 159 "Times New Roman","serif"'>constant </span> <span style='font-family:"Times New Roman","serif"'>and 247 160 </span><span style='position:relative;top:6.0pt'><img width=28 height=24 248 src=" sm_image015.gif"></span>= <span style='font-family:161 src="./img/sm_image015.gif"></span> = <span style='font-family: 249 162 "Times New Roman","serif"'>constant:</span></p> 250 163 … … 256 169 W, which will take minutes and minutes to finish the calculation for a set of 257 170 typical SANS data. An alternative way which is correct for slit width << 258 slit hight, is used in the SANSView: 171 slit hight, is used in the SANSView: This method is a mixed method that 259 172 combines the method 1) with the numerical integration for the slit width.</span></p> 260 173 … … 263 176 264 177 <p class=MsoListParagraphCxSpMiddle style='margin-left:.25in'> 265 <img 266 src="sm_image020.gif"> <span style='font-family: 178 <img src="./img/sm_image020.gif"> <span style='font-family: 267 179 "Times New Roman","serif"'>(7)</span></p> 268 180 269 181 <p class=MsoListParagraphCxSpMiddle style='margin-left:.25in'><span 270 style='font-family:"Times New Roman","serif"'>for 182 style='font-family:"Times New Roman","serif"'>for q<sub>p</sub> = q<sub>i</sub> 271 183 - </span><span style='position:relative;top:6.0pt'><img width=28 height=24 272 src=" sm_image012.gif"></span><span style='font-family:273 "Times New Roman","serif"'> 184 src="./img/sm_image012.gif"></span><span style='font-family: 185 "Times New Roman","serif"'> and</span> <span style='font-family:"Times New Roman","serif"'>q<sub>N</sub> 274 186 = q<sub>i</sub> + </span><span style='position:relative;top:6.0pt'><img 275 width=28 height=24 src=" sm_image012.gif"></span>.<span187 width=28 height=24 src="./img/sm_image012.gif"></span>. <span 276 188 style='position:relative;top:7.0pt'><img width=23 height=25 277 src=" sm_image018.gif"></span><span style='font-family:189 src="./img/sm_image018.gif"></span> <span style='font-family: 278 190 "Times New Roman","serif"'>= 0 for <i>I<sub>s</sub></i> in</span> <i><span 279 191 style='font-family:"Times New Roman","serif"'>j < p</span></i><span … … 294 206 case becomes</span></p> 295 207 296 <p class=MsoNormal><img 297 src="sm_image021.gif"><span 298 style='font-family:"Times New Roman","serif"'> (8)</span></p> 208 <p class=MsoNormal><img src="./img/sm_image021.gif"><span 209 style='font-family:"Times New Roman","serif"'> (8)</span></p> 299 210 300 211 <p class=MsoNormal><span style='font-family:"Times New Roman","serif"'>For all … … 303 214 evenly) below (>0) and above the q range of data in order to cover all data 304 215 points of the smearing computation for a given model and for a given slit size. 305 216 The <i>Norm</i> factor is found numerically with the weighting matrix, and 306 217 considered on <i>I<sub>s</sub></i> computation.</span></p> 307 218 … … 316 227 except that the weight function used was the 2D elliptical Gaussian function</span></p> 317 228 318 <p class=MsoNormal><img 319 src="sm_image022.gif"><span 320 style='font-family:"Times New Roman","serif"'> (9)</span></p> 229 <p class=MsoNormal><img src="./img/sm_image022.gif"><span 230 style='font-family:"Times New Roman","serif"'> (9)</span></p> 321 231 322 232 <p class=MsoNormal><span style='font-family:"Times New Roman","serif"'>In Eq … … 326 236 , and the primed axes are in the coordinate rotated by an angle </span><span 327 237 style='font-family:Symbol'>theta</span><span style='font-family:"Times New Roman","serif"'> 328 around z-axis (below) so that x<sub>0</sub> = 238 around z-axis (below) so that x<sub>0</sub> = x<sub>0</sub>cos</span><span 329 239 style='font-family:Symbol'>(theta) + </span><span style='font-family:"Times New Roman","serif"'>y<sub>0</sub> 330 sin</span><span style='font-family:Symbol'>(theta) 331 "Times New Roman","serif"'>and y<sub>0</sub> = 240 sin</span><span style='font-family:Symbol'>(theta) </span><span style='font-family: 241 "Times New Roman","serif"'>and y<sub>0</sub> = -x<sub>0</sub>sin</span><span 332 242 style='font-family:Symbol'>(theta) + </span><span style='font-family:"Times New Roman","serif"'>y<sub>0</sub> 333 243 cos</span><span style='font-family:Symbol'>(theta) .</span><span style='font-family: 334 244 "Times New Roman","serif"'> Note that the rotation angle is zero for x-y 335 245 symmetric elliptical Gaussian distribution</span><span style='font-family:Symbol'>. 336 </span><span style='font-family:"Times New Roman","serif"'>The 246 </span><span style='font-family:"Times New Roman","serif"'>The A is a 337 247 normalization factor.</span></p> 338 248 339 249 <p class=MsoNormal align=center style='text-align:center'><span 340 250 style='font-family:"Times New Roman","serif"'><img width=439 height=376 341 id="Object 1" src=" sm_image023.gif"></span></p>251 id="Object 1" src="./img/sm_image023.gif"></span></p> 342 252 343 253 <p class=MsoNormal><span style='font-family:"Times New Roman","serif"'> </span></p> … … 352 262 assumed that I(x, y) is constant within the bins which in turn becomes</span></p> 353 263 354 <p class=MsoNormal><img 355 src="sm_image024.gif"></p> 356 357 <p class=MsoNormal> <span 264 <p class=MsoNormal><img src="./img/sm_image024.gif"></p> 265 266 <p class=MsoNormal> <span 358 267 style='font-family:"Times New Roman","serif"'>(10)</span></p> 359 268 … … 362 271 transform x-y back to x-y coordinate (rotating it by -</span><span 363 272 style='font-family:Symbol'>(theta)</span><span style='font-family:"Times New Roman","serif"'> 364 around z axis). Then, for the polar symmetric smear,</span></p> 365 366 <p class=MsoNormal><img 367 src="sm_image025.gif"><span 368 style='position:relative;top:35.0pt'> </span>(11)</p> 273 around z axis). Then, for the polar symmetric smear,</span></p> 274 275 <p class=MsoNormal><img src="./img/sm_image025.gif"><span 276 style='position:relative;top:35.0pt'> </span>(11)</p> 369 277 370 278 <p class=MsoNormal><span style='font-family:"Times New Roman","serif"'>where,</span></p> 371 279 372 <p class=MsoNormal><img 373 src="sm_image026.gif"></p> 280 <p class=MsoNormal><img src="./img/sm_image026.gif"></p> 374 281 375 282 <p class=MsoNormal><span style='font-family:"Times New Roman","serif"'>while 376 283 for the x-y symmetric smear,</span></p> 377 284 378 <p class=MsoNormal><img 379 src="sm_image027.gif"><span 380 style='font-family:"Times New Roman","serif"'> (12)</span></p> 285 <p class=MsoNormal><img src="./img/sm_image027.gif"><span 286 style='font-family:"Times New Roman","serif"'> (12)</span></p> 381 287 382 288 <p class=MsoNormal><span style='font-family:"Times New Roman","serif"'>where,</span></p> 383 289 384 <p class=MsoNormal><img 385 src="sm_image028.gif"></p> 290 <p class=MsoNormal><img src="./img/sm_image028.gif"></p> 386 291 387 292 <p class=MsoNormal><span style='font-family:"Times New Roman","serif"'>Here, the
Note: See TracChangeset
for help on using the changeset viewer.