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

GetDetectorOffsets v1

Summary

Creates an OffsetsWorkspace containing offsets for each detector. You can then save these to a .cal file using SaveCalFile.

See Also

AlignComponents

Properties

Name	Direction	Type	Default	Description
InputWorkspace	Input	MatrixWorkspace	Mandatory	A 2D workspace with X values of d-spacing
Step	Input	number	0.001	Step size used to bin d-spacing data
DReference	Input	number	2	Center of reference peak in d-space
XMin	Input	number	0	Minimum of CrossCorrelation data to search for peak, usually negative
XMax	Input	number	0	Maximum of CrossCorrelation data to search for peak, usually positive
GroupingFileName	Input	string		Optional: The name of the output CalFile to save the generated OffsetsWorkspace. Allowed extensions: [‘.cal’]
OutputWorkspace	Output	OffsetsWorkspace	Mandatory	An output workspace containing the offsets.
MaskWorkspace	Output	MaskWorkspace	_empty_	Mask workspace (optional input / output workspace): when specified, if the workspace already exists, any incoming masked detectors will be combined with any additional outgoing masked detectors detected by the algorithm
PeakFunction	Input	string	Gaussian	The function type for fitting the peaks. Allowed values: [‘AsymmetricPearsonVII’, ‘BackToBackExponential’, ‘Bk2BkExpConvPV’, ‘DeltaFunction’, ‘ElasticDiffRotDiscreteCircle’, ‘ElasticDiffSphere’, ‘ElasticIsoRotDiff’, ‘ExamplePeakFunction’, ‘Gaussian’, ‘IkedaCarpenterPV’, ‘Lorentzian’, ‘PearsonIV’, ‘PseudoVoigt’, ‘Voigt’]
EstimateFWHM	Input	boolean	False	Whether to esimate FWHM of peak function when estimating fit parameters
MaxOffset	Input	number	1	Maximum absolute value of offsets; default is 1
OffsetMode	Input	string	Relative	Whether to calculate a relative, absolute, or signed offset. Allowed values: [‘Relative’, ‘Absolute’, ‘Signed’]
DIdeal	Input	number	2	The known peak centre value from the NIST standard information, this is only used in Absolute OffsetMode.

Description

This algorithm requires a workspace that is both in d-spacing, but has also been preprocessed by the CrossCorrelate v1 algorithm. In this first step you select one spectrum to be the reference spectrum and all of the other spectrum are cross correlated against it. Each output spectrum then contains a peak whose location defines the offset from the reference spectrum.

The algorithm iterates over each spectrum in the workspace and fits a PeakFunction (default is a Gaussian function) to the reference peaks. The fit is used to calculate the centre of the fitted peak. Using the fitted peak center, there are 3 options for calculating offset, depending on the choice for the “OffsetMode” input parameter:

Relative (offset relative to reference position):

\(offset = -peakCentre*step/(dreference+PeakCentre*step)\)

Absolute (offset relative to ideal position):

\(offset = -peakCentre*step/(dreference+PeakCentre*step) + (dideal - dreference) / dreference\)

Signed (offset in raw number of bins):

\(offset = -peakCentre\)

This is then written into a .cal file for every detector that contributes to that spectrum. All of the entries in the cal file are initially set to both be included, but also to all group into a single group on DiffractionFocussing v2. The CreateCalFileByNames v1 algorithm can be used to alter the grouping in the cal file.

Usage

import os

# Create a workspace with a Gaussian peak in the centre.
ws = CreateSampleWorkspace(Function='User Defined',UserDefinedFunction='name=Gaussian,Height=1,PeakCentre=10,Sigma=1',XMin=0,XMax=20,BinWidth=0.1)
ws.getAxis(0).setUnit( 'dSpacing' )

# Generate a file path to save the .cal file at.
calFilePath = os.path.expanduser( '~/MantidUsageExample_CalFile.cal' )

# Run the algorithm
msk = GetDetectorOffsets(ws,0.001,10.0,0, 10, calFilePath)

# Read the saved .cal file back in
f = open( calFilePath, 'r' )
file = f.read().split('\n')
f.close()

# Print out first 10 lines of the file
print("{} ...".format(file[0][:55]))
for line in file[1:10]:
    print(line)

Output

# Calibration file for instrument basic_rect written on ...
# Format: number    UDET         offset    select    group
        0            ...     ...       1       1
        1            ...     ...       1       1
        2            ...     ...       1       1
        3            ...     ...       1       1
        4            ...     ...       1       1
        5            ...     ...       1       1
        6            ...     ...       1       1
        7            ...     ...       1       1

Categories: AlgorithmIndex | Diffraction\Calibration

Source

C++ header: GetDetectorOffsets.h

C++ source: GetDetectorOffsets.cpp