FindSXPeaks v1

../_images/FindSXPeaks-v1_dlg.png

FindSXPeaks dialog.

Summary

Takes a 2D workspace as input and finds the highest intensity point in each 1D spectrum. This is used in particular for single crystal as a quick way to find strong peaks.

Properties

Name Direction Type Default Description
InputWorkspace Input MatrixWorkspace Mandatory The name of the Workspace2D to take as input
RangeLower Input number Optional The X value to search from (default 0)
RangeUpper Input number Optional The X value to search to (default total number of bins)
StartWorkspaceIndex Input number 0 Start workspace index (default 0)
EndWorkspaceIndex Input number Optional End workspace index (default to total number of histograms)
PeakFindingStrategy Input string StrongestPeakOnly Different options for peak finding.1. StrongestPeakOnly: Looks only for the the strongest peak in each spectrum (provided there is one). This options is more performant than the AllPeaks option. 2. AllPeaks: This strategy will find all peaks in each spectrum. This is slower than StrongestPeakOnly. Note that the recommended ResolutionStrategy in this mode is AbsoluteResolution. Allowed values: [‘StrongestPeakOnly’, ‘AllPeaks’]
SignalBackground Input number 10 Multiplication factor for the signal background. Peaks which are below the estimated background are discarded. The background is estimated to be an average of the first and the last signal and multiplied by the SignalBackground property.
AbsoluteBackground Input number 30 Peaks which are below the specified absolute background are discarded. The background is gloabally specified for all spectra. Inspect your data in the InstrumentView to get a good feeling for the background threshold. Background thresholds which are too low will mistake noise for peaks.
ResolutionStrategy Input string RelativeResolution Different options for the resolution.1. RelativeResolution: This defines a relative tolerance needed to avoid peak duplication in number of pixels. This selection will enable the Resolution property and disable the XResolution, PhiResolution, ThetaResolution. 1. AbsoluteResolution: This defines an absolute tolerance needed to avoid peak duplication in number of pixels. This selection will disable the Resolution property and enable the XResolution, PhiResolution, ThetaResolution. Allowed values: [‘RelativeResolution’, ‘AbsoluteResolution’]
Resolution Input number 0.01 Tolerance needed to avoid peak duplication in number of pixels
XResolution Input number 0 Absolute tolerance in time-of-flight or d-spacing needed to avoid peak duplication in number of pixels. The values are specified in either microseconds or angstroms.
PhiResolution Input number 1 Absolute tolerance in the phi coordinate needed to avoid peak duplication in number of pixels. The values are specified in degrees.
TwoThetaResolution Input number 1 Absolute tolerance of two theta value needed to avoid peak duplication in number of pixels. The values are specified in degrees.
OutputWorkspace Output PeaksWorkspace   The name of the PeaksWorkspace in which to store the list of peaks found

Description

Detector-space, single crystal peak finding. The algorithm can be used to search through each spectra and find either all peaks or only the strongest peak within a given x range. The x units of each spectrum can either be TOF or d-spacing.

Strongest peak mode:

  • The highest intensity bin is taken to be the peak, so the algorithm only finds one peak per spectra. In this case the background is calculated as the average of the start and end intensity multiplied by the provided SignalBackground parameter. Peaks that are not above the background are culled.
  • This peak finding mode is mainly used with the ResolutionStrategy set to RelativeResolution. In this case, the resolution parameter defines a tolerance which is compared to the absolute difference between the parameters \phi, 2\theta, and t of two found peaks. If the absolute difference between any of the parameters for two peaks is greater than the product of the tolerance value and the parameter value then the two peaks are classed as not the same. i.e. if |\phi_1 - \phi_2| > tolerance * \phi_1 then peaks 1 & 2 are not the same (as well as similar definitions for 2\theta and t).

All peaks mode:

  • All peaks in each spectra above a certain threshold are detected as peaks (see below for the threshold).
  • In this case the user provides an absolute, global background/threshold value which is set with the AbsoluteBackground parameter. Peaks that are not above the absolute background are culled. A good way of identifying the background is to inspect your data set in pick mode within the InstrumentView. Ensure that your absolute background has been set high enough, else the algorithm will pick up noise as peaks.
  • This peak finding mode provides the best results when used with the ResolutionStrategy set to AbsoluteResolution, which allows the user to set absolute resolution values for the XUnitResolution (either in units of TOF or d-spacing), PhiResolution and the TwoThetaResolution parameters. These resolution parameters define tolerances which are compared to the absolute difference between the parameters \phi, 2\theta, and t of two found peaks. If the absolute difference between any of the parameters for two peaks is greater than the absolute tolerance then the two peaks are classed as not the same. i.e. if |\phi_1 - \phi_2| >  PhiTolerance then peaks 1 & 2 are not the same (as well as similar definitions for 2\theta and t).

General points:

  • Calculated Qlab follows the Busy, Levy 1967 convention.

Usage

Example

# create histogram workspace
ws=CreateSampleWorkspace()

wsPeaks = FindSXPeaks(ws)

print "Peaks found: " + str(wsPeaks.getNumberPeaks())

Output:

Peaks found: 174

Categories: Algorithms | Crystal\Peaks | Optimization\PeakFinding

Source

C++ source: FindSXPeaks.cpp

C++ header: FindSXPeaks.h