\(\renewcommand\AA{\unicode{x212B}}\)

SANSWideAngleCorrection v1

../_images/ImageNotFound.png

Enable screenshots using DOCS_SCREENSHOTS in CMake

Summary

Calculate the Wide Angle correction for SANS transmissions.

Properties

Name

Direction

Type

Default

Description

SampleData

Input

MatrixWorkspace

Mandatory

A workspace cropped to the detector to be reduced (the SAME as the input to [[Q1D]]); used to verify the solid angle. The workspace is not modified, just inspected.

TransmissionData

Input

MatrixWorkspace

Mandatory

The transmission data calculated, referred to as <math>T_0</math> in equations in discussion section

OutputWorkspace

Output

MatrixWorkspace

Mandatory

The transmission corrected SANS data, normalised (divided) by <math>T_0</math>, see discussion section

Description

Algorithm Reference Discussion

Looking at Computing guide for Small Angle Scattering Experiments by Ghosh, Egelhaaf & Rennie, we see that for scattering at larger angles the transmission should be modified due to the longer path length after the scattering event.

The longer path length after scattering will also slightly increase the probability of a second scattering event, but this is not dealt with here.

If our on-axis transmission is \(T_0\) through a sample of thickness \(d\), then the transmission at some other thickness \(x\) is \(\exp(-\mu x)\) where attenuation coefficient \(\mu = -\ln( \frac{T_0}{d})\).

If a neutron scatters at angle \(2\theta\) at distance \(x\) into the sample, its total transmission is then:

\(T^{''} = \exp(-\mu x) \exp( \frac{-\mu(d-x)}{\cos(2\theta)})\)

\(T^{''}\) should be integrated and averaged between \(x = 0\) and \(x = d\).

Hammouda, gives an approximate result for the integral, see page 208 of SANS toolbox:

\(T^{'} = \frac{T_0(T_0^A - 1)}{A \ln(T_0)}\)

For:

\(A = \frac{1}{\cos(2\theta)} - 1\)

For example if \(T_0 = 0.2\) and \(2\theta = 40\) then \(T^{'} = 0.158\), a shift of \(~20\)% of the SANS curve. Note that the result is independent of sample thickness.

\(T_0\) is a function of neutron wavelength, whilst \(A\) is a function of detector pixel location.

The output of this algorithm is:

\(OutputWorkspace = \frac{T^\prime}{T_0}\)

Error Propagation

The error propagation follows this formula:

\(OutputWorkspace_{error} = \frac{T_{0E} ^A - 1}{A\ln(T_0E)}\)

Which means, that we do not consider the error in the definition of the \(2\theta\) (the parameter A)

Enabling Wide Angle Correction for Reduction of SANS ISIS

To enable the Wide Angle correction use the User File settings:

SAMPLE/PATH/ON

More information on: SANS_User_File_Commands#SAMPLE

SANS ISIS Reduction

The output of SANSWideAngleCorrection is used as WavePixelAdj parameter at Q1D v2.

Wide Angle Correction and the SANS Reduction

The equation for the reduction is (see Q1D v2)

\(P_I(Q) = \frac{\sum_{\{i, j, n\} \supset \{I\}} S(i,j,n)}{\sum_{\{i, j, n\} \supset \{I\}}M(n)\eta(n)T(n)\Omega_{i j}F_{i j}}\)

But, \(T(n)\) is not really \(T(n)\), because of the wide angles, it is now \(T(n,theta)\) or \(T(n,i,j)\).

So, we decided to have a new factor that changes this equation to:

\(P_I(Q) = \frac{\sum_{\{i, j, n\} \supset \{I\}} S(i,j,n)}{\sum_{\{i, j, n\} \supset \{I\}}M(n)\eta(n)T(n)\Omega_{i j}F_{i j}Corr(i,j,n)}\)

Where Corr (Correction factor) in this case will be:

\(Corr = \frac{T_0^A - 1}{A \ln(T_0)}\)

Which is the OutputWorkspace of SANSWideAngleCorrection.

This parameter enters inside Q1D v2 as WavePixelAdj. But, this is all done for you inside the Reduction Script.

Comparison with Wide Angle Correction at SNS

The transmission correction applied at SNS is described here, and it is applied through the ApplyTransmissionCorrection v1 algorithm. The correction applied there is an approximation for the same equations described here. The picture above compare their results

SNS_ISIS_WideAngleCorrections.png

Note a difference among them is when they are applied. At SNS, the correction is applied before averaging the counters per bin inside Q1D v2 algorithm, while at ISIS, it is used after, inside the Q1D v2 algorithm, for the division of the counters per bin normalized by the transmission counters.

References

Annie Brulet et al. - Improvement of data treatment in small-angle neutron scattering - J. Appl. Cryst. (2007). 40

Ghosh, Egelhaaf & Rennie - Computing guide for Small Angle Scattering Experiments

Usage

Example - Correcting Some Dummy Data

# Create some dummy data, but crop it for quick demonstration purposes.
sample = CreateSimulationWorkspace(Instrument='SANS2D', BinParams=[5,500,100005], UnitX='TOF')
sample = CropWorkspace(sample,StartWorkspaceIndex=0,EndWorkspaceIndex=20)

# Create a dummy transmission workspace.
transmission = CropWorkspace(sample,StartWorkspaceIndex=10,EndWorkspaceIndex=10)
transmission *= 2

corrected_data = SANSWideAngleCorrection(sample, transmission)

print("{:.6f} was corrected to {:.6f}.".format(sample.readY(19)[0], corrected_data.readY(19)[0]))

Output:

1.000000 was corrected to 1.004997.

Categories: AlgorithmIndex | CorrectionFunctions\TransmissionCorrections

Source

Python: SANSWideAngleCorrection.py