He3TubeEfficiency v1

../_images/He3TubeEfficiency-v1_dlg.png

He3TubeEfficiency dialog.

Summary

He3 tube efficiency correction.

Properties

Name Direction Type Default Description
InputWorkspace Input MatrixWorkspace Mandatory Name of the input workspace
OutputWorkspace Output MatrixWorkspace Mandatory Name of the output workspace, can be the same as the input
ScaleFactor Input number 1 Constant factor with which to scale the calculateddetector efficiency. Same factor applies to all efficiencies.
TubePressure Input dbl list   Provide overriding the default tube pressure. The pressure must be specified in atm.
TubeThickness Input dbl list   Provide overriding the default tube thickness. The thickness must be specified in metres.
TubeTemperature Input dbl list   Provide overriding the default tube temperature. The temperature must be specified in Kelvin.

Description

This algorithm corrects the detection efficiency of He3 tubes using an exponential function and certain detector properties. The correction scheme is given by the following:

\epsilon = \frac{A}{1-e^{\frac{-\alpha P (L - 2W) \lambda}{T sin(\theta)}}}

where A is a dimensionless scaling factor, \alpha is a constant with units (Kelvin / (metres\: \AA\: atm)), P is pressure in units of atm, L is the tube diameter in units of metres, W is the tube thickness in units of metres, T is the temperature in units of Kelvin, sin(\theta) is the angle of the neutron trajectory with respect to the long axis of the He3 tube and \lambda is in units of \AA.

The Optional properties that are of num list type can be used in the following manner. If no input value is given, the detector parameter is pulled from the detector itself. If a single value is used as input, that value is applied to all detectors. If an array of values is used, that array must be the same size as the number of spectra in the workspace. If it is not, the spectra indicies that do not have values will be zeroed in the output workspace.

Restrictions on Input Workspace

The input workspace must be in units of wavelength.

Usage

The correction applies equally well to histogram and event workspaces. An example for a histogram workspace will be given.

# create workspace with correct units
ws = CreateSampleWorkspace(Function='Flat background',NumBanks='1',BankPixelWidth='2',
                           XUnit='Wavelength',XMax='10',BinWidth='0.01')

# workspace does not have correct tube parameters, so we need to provide some
# parameters taken from the ARCS instrument at the SNS
ws1 = He3TubeEfficiency(ws, TubePressure='10',TubeThickness='0.0008',TubeTemperature='290')

# show the before and after for the first detector for the first 5 bins
print "Original counts:", str(ws.readY(0)[:5])
print "Corrected counts:", str(ws1.readY(0)[:5])

Output:

Original counts: [ 1.  1.  1.  1.  1.]
Corrected counts: [ 417.07353051  139.35837696   83.81566631   60.01187609   46.787726  ]

Categories: Algorithms | CorrectionFunctions | EfficiencyCorrections