Smearing Computation

Slit_Smearing

Pinhole_Smearing

2D_Smearing

Slit Smearing

The sit smeared scattering intensity for SAS is defined by

where Norm =

Equation 1

The functions and refer to the slit width weighting function and the slit height weighting determined at the q point, respectively. Here, we assumes that the weighting function is described by a rectangular function, i.e.,

Equation 2

and

Equation 3

so that for and u.

The and stand for the slit height (FWHM/2) and the slit width (FWHM/2) in the q space. Now the integral of Equation 1 is simplified to

Equation 4

Numerical Implementation of Equation 4

Case 1

For = 0 and = constant.

For discrete q values, at the q values from the data points and at the q values extended up to qN= qi + the smeared intensity can be calculated approximately

Equation 5

= 0 for Is in j < i or j > N-1*.

Case 2

For = constant and = 0.

Similarly to Case 1, we get

for qp= qi- and qN= qi+ . = 0 for Is in j < p or j > N-1.

Case 3

For = constant and = constant.

In this case, the best way is to perform the integration, Equation 1, numerically for both slit height and width. However, the numerical integration is not correct enough unless given a large number of iteration, say at least 10000 by 10000 for each element of the matrix, W, which will take minutes and minutes to finish the calculation for a set of typical SAS data. An alternative way which is correct for slit width << slit hight, is used in SasView. This method is a mixed method that combines method 1 with the numerical integration for the slit width.

Equation 7

for qp= qi- and qN= qi+ . = 0 for Is in j < p or j > N-1.

Pinhole Smearing

The pinhole smearing computation is done similar to the case above except that the weight function used is the Gaussian function, so that the Equation 6 for this case becomes

Equation 8

For all the cases above, the weighting matrix W is calculated when the smearing is called at the first time, and it includes the ~ 60 q values (finely binned evenly) below (>0) and above the q range of data in order to cover all data points of the smearing computation for a given model and for a given slit size. The Norm factor is found numerically with the weighting matrix, and considered on Is computation.

2D Smearing

The 2D smearing computation is done similar to the 1D pinhole smearing above except that the weight function used was the 2D elliptical Gaussian function

Equation 9

In Equation 9, x0 = qcos/theta/ and y0 = qsin/theta/, and the primed axes are in the coordinate rotated by an angle /theta/ around the z-axis (below) so that x’0= x0cos/theta/+y0sin/theta/ and y’0= -x0sin/theta/+y0cos/theta/.

Note that the rotation angle is zero for x-y symmetric elliptical Gaussian distribution. The A is a normalization factor.

Now we consider a numerical integration where each bins in /theta/ and R are evenly (this is to simplify the equation below) distributed by /delta//theta/ and /delta/R, respectively, and it is assumed that I(x’, y’) is constant within the bins which in turn becomes

Equation 10

Since we have found the weighting factor on each bin points, it is convenient to transform x’-y’ back to x-y coordinate (rotating it by -/theta/ around z axis). Then, for the polar symmetric smear

Equation 11

where

while for the x-y symmetric smear

Equation 12

where

Here, the current version of the SasView uses Equation 11 for 2D smearing assuming that all the Gaussian weighting functions are aligned in the polar coordinate.

In the control panel, the higher accuracy indicates more and finer binnng points so that it costs more in time.