\(\renewcommand\AA{\unicode{x212B}}\)
NormaliseToMonitor v1¶
Summary¶
Normalizes a 2D workspace by a specified spectrum, spectrum, described by a monitor ID or spectrun provided in a separate worskspace.
See Also¶
Properties¶
Name |
Direction |
Type |
Default |
Description |
---|---|---|---|---|
InputWorkspace |
Input |
Mandatory |
Name of the input workspace. Must be a non-distribution histogram. |
|
OutputWorkspace |
Output |
Mandatory |
Name to use for the output workspace |
|
MonitorSpectrum |
InOut |
number |
-1 |
The spectrum number within the InputWorkspace you want to normalize by (It can be a monitor spectrum or a spectrum responsible for a group of detectors or monitors) |
MonitorID |
Input |
number |
-1 |
The MonitorID (detector ID), which defines the monitor’s data within the InputWorkspace. Will be overridden by the values correspondent to MonitorSpectrum field if one is provided in the field above. If workspace do not have monitors, the MonitorID can refer to empty data and the field then can accepts any MonitorID within the InputWorkspace. |
MonitorWorkspace |
Input |
A workspace containing one or more spectra to normalize the InputWorkspace by. |
||
MonitorWorkspaceIndex |
InOut |
number |
0 |
The index of the spectrum within the MonitorWorkspace(2 (0<=ind<=nHistograms in MonitorWorkspace) you want to normalize by (usually related to the index, responsible for the monitor’s data but can be any). If no value is provided in this field, ‘’’InputWorkspace’’’ will be normalized by first spectra (with index 0) |
IntegrationRangeMin |
Input |
number |
Optional |
If set, normalization will be by integrated count from this minimum x value |
IntegrationRangeMax |
Input |
number |
Optional |
If set, normalization will be by integrated count up to this maximum x value |
IncludePartialBins |
Input |
boolean |
False |
If true and an integration range is set then partial bins at either end of the integration range are also included |
NormFactorWS |
Output |
Name of the workspace, containing the normalization factor. If this name is empty, normalization workspace is not returned. If the name coincides with the output workspace name, _normFactor suffix is added to this name |
Description¶
Bin-by-bin mode¶
In this, the default scenario, each spectrum in the workspace is normalised on a bin-by-bin basis by the monitor spectrum given. The error on the monitor spectrum is taken into account. The normalisation scheme used is:
\((s_i)_{Norm}=(\frac{s_i}{m_i})*\Delta w_i*\frac{\sum_i{m_i}}{\sum_i(\Delta w_i)}\)
where \(s_i\) is the signal in a bin, \(m_i\) the count in the corresponding monitor bin, \(\Delta w_i\) the bin width, \(\sum_i{m_i}\) the integrated monitor count and \(\sum_i{\Delta w_i}\) the sum of the monitor bin widths. In words, this means that after normalisation each bin is multiplied by the bin width and the total monitor count integrated over the entire frame, and then divided by the total frame width. This leads to a normalised histogram which has unit of counts, as before.
If the workspace does not have common binning, then the monitor spectrum is rebinned internally to match each data spectrum prior to doing the normalisation.
If a bin in the monitor spectrum contains 0 counts there is no special treatment, only a warning is logged. The following options are available to deal with the issue:
Smooth the monitor spectrum before normalisation (using for example FFTSmooth v2, SmoothData v1, SplineSmoothing v1 algorithms).
Replace the infinite numbers and NaNs after normalisation (ReplaceSpecialValues v1 algorithm).
Normalisation by integrated count mode¶
This mode is used if one or both of the relevant ‘IntegrationRange’ optional parameters are set. If either is set to a value outside the workspace range, then it will be reset to the frame minimum or maximum, as appropriate.
The error on the integrated monitor spectrum is taken into account in the normalisation. No adjustment of the overall normalisation takes place, meaning that the output values in the output workspace are technically dimensionless.
Restrictions on the input workspace¶
The data must be histogram, non-distribution data. The exception to the histogram requirement is for workspaces that contain point data with a single count per spectrum. In this case the normalisation is performed by dividing every spectrum by the monitor counts, taking into account the error on the monitor counts.
Detector Scan Workspaces¶
Workspaces that have scanning detectors are supported by this algorithm, both for bin-by-bin mode and normlisation by integrated count. The only option for specifying the monitor is by ‘MonitorID’, attempting to use ‘MonitorSpectrum’ or MonitorWorkspaceIndex’ will throw an error. In this case the ‘NormFactorWS’ output will contain a monitor spectrum for each time index.
Child Algorithms used¶
The ExtractSpectra v1 algorithm is used to pull out the monitor spectrum if it’s part of the InputWorkspace or MonitorWorkspace. For the ‘integrated range’ option, the Integration v1 algorithm is used to integrate the monitor spectrum.
Usage¶
Note
To run these usage examples please first download the usage data, and add these to your path. In Mantid this is done using Manage User Directories.
Example - Normalise to Monitor ID=1
ws =Load('IRS26173.raw')
wsN = NormaliseToMonitor( ws, MonitorID=1 )
print("Without normalisation")
print("Monitor ID=1 {:.3f}, {:.3f}".format(ws.readY(0)[0], ws.readY(0)[1]))
print("Selected data {:.6f}, {:.6f}".format(ws.readY(6)[0], ws.readY(3)[1]))
print("With Normalisation")
print("Monitor ID=1 {:.3f}, {:.3f}".format(wsN.readY(0)[0], wsN.readY(0)[1]))
print("Selected data {:.6f}, {:.6f}".format(wsN.readY(6)[0], wsN.readY(3)[1]))
Output:
Without normalisation
Monitor ID=1 626034.000, 626681.000
Selected data 2.000000, 1.000000
With Normalisation
Monitor ID=1 464872.441, 464872.441
Selected data 1.485135, 0.741801
Categories: AlgorithmIndex | CorrectionFunctions\NormalisationCorrections
Source¶
C++ header: NormaliseToMonitor.h
C++ source: NormaliseToMonitor.cpp