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

TOFSANSResolutionByPixel v1

../_images/ImageNotFound.png

Enable screenshots using DOCS_SCREENSHOTS in CMake

Summary

Calculate the Q resolution for TOF SANS data for each pixel.

See Also

TOFSANSResolution

Properties

Name

Direction

Type

Default

Description

InputWorkspace

Input

MatrixWorkspace

Mandatory

Name the workspace to calculate the resolution for, for each pixel and wavelength

OutputWorkspace

Output

Workspace

Mandatory

Name of the newly created workspace which contains the Q resolution.

DeltaR

Input

number

0

Virtual ring width on the detector (mm).

SampleApertureRadius

Input

number

0

Sample aperture radius, R2 (mm).

SourceApertureRadius

Input

number

0

Source aperture radius, R1 (mm).

SigmaModerator

Input

MatrixWorkspace

Mandatory

Moderator time spread (microseconds) as afunction of wavelength (Angstroms).

CollimationLength

Input

number

0

Collimation length (m)

AccountForGravity

Input

boolean

False

Whether to correct for the effects of gravity

ExtraLength

Input

number

0

Additional length for gravity correction.

Description

Calculates the Q-resolution per pixel according to Mildner and Carpenter equation

\[(\sigma_Q )^2 = \frac{4\pi^2}{12\lambda^2} [ 3(\frac{R_1}{L_1})^2 + 3(\frac{R_2}{L_3})^2 + (\frac{\Delta R}{L_2})^2 ] + Q^2(\frac{\sigma_{\lambda}}{\lambda})^2\]

where \(L1\) and \(L2\) are the collimation length and sample-to-detector distance respectively and

\[\frac{1}{L_3} = \frac{1}{L_1} + \frac{1}{L_2}\]

and

\[(\sigma_{\lambda})^2 = (\Delta \lambda )^2 / 12 + (\sigma_{moderator})^2\]

where \(\sigma_{\lambda}\) is the overall effective standard deviation in wavelength. \(\Delta \lambda\) values are found from the wavelength binning of the InputWorkspace, \(\sigma_{moderator}\) is the moderator time spread (the variation in time for the moderator to emit neutrons of a given wavelength). Note that \(\Delta \lambda\) may be imposed by wavelength steps set elsewhere in Mantid which should be at least as large as the equivalent time bins used in the original histogram data collection. For event mode data \(\Delta \lambda\) is in theory very small, but in practice a histogram in time has to be generated (perhaps using monitor time bins or specifically set event-time-bins), before a rebinning into user provided wavelength steps in InputWorkspace. Again the latter steps should be the largest.

Q values needed here are calculated in the same way as for Q1D, including correction for gravity for which detector coordinates are assumed centred at zero wavelength.

\(\sigma_Q\) is returned as the y-values of the InputWorkspace, and the remaining variables in the main equation above are related to parameters of this algorithm as follows:

  • \(R_1\) equals SourceApertureRadius

  • \(R_2\) equals SampleApertureRadius

  • \(\Delta R\) equals DeltaR

  • \(\sigma_{moderator}\) equals SigmaModerator

  • \(\L_1\) equals CollimationLength

\(\lambda\) in the equation is the midpoint of wavelength histogram bin values of InputWorkspace.

Collimation length \(L_1\) in metres in the equation here is the distance between the first beam defining pinhole (Radius \(R_1\)) and the sample aperture (radius \(R_2\)). (Beware that \(L_1\) is more often the moderator to sample distance.)

For rectangular collimation apertures, size H x W, Mildner & Carpenter say to use \(R = \sqrt{( H^2 +W^2)/6 }\). Note that we are assuming isotropically averaged, scalar \(Q\), and making some small angle approximations. Results on higher angle detectors may not be accurate. For data reduction sliced in different directions on the detector (e.g. GISANS) adjust the calling parameters to suit the collimation in that direction.

Note that \(\Delta\) is the full width of a rectangular distribution in radius or wavelength, for which the standard deviation is \(\sigma=\Delta/\sqrt{12}\). For a Gaussian distribution the FWHM (full width at half maximum) is \(\sqrt{8\ln{2}}\sigma=2.35482\sigma\). For an exponential decay \(e^{-t/\tau}\), the standard deviation (and the mean) is \(\tau\). For non-rectangular distributions these equations allow the equivalent \(\Delta\) to be entered as \(\Delta=\sqrt{12}\sigma\).

This version of the algorithm neglects wavelength-dependent detector detection depth effects.

Categories: AlgorithmIndex | SANS

Source

C++ header: TOFSANSResolutionByPixel.h

C++ source: TOFSANSResolutionByPixel.cpp