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

IdentifyNoisyDetectors v1

../_images/ImageNotFound.png

Enable screenshots using DOCS_SCREENSHOTS in CMake

Summary

This algorithm creates a single-column workspace where the Y values are populated withs 1s and 0s, 0 signifying that the detector is to be considered ‘bad’ based on the method described below.

See Also

CreatePSDBleedMask

Properties

Name

Direction

Type

Default

Description

InputWorkspace

Input

MatrixWorkspace

Mandatory

OutputWorkspace

Output

MatrixWorkspace

Mandatory

RangeLower

Input

number

2000

The lower integration range

RangeUpper

Input

number

19000

The upper integration range

Description

The process for this algorithm is:

  • The standard deviation for each pixel within the pre-set range is calculated.

  • The mean value and standard deviation of these standard deviation values is calculated.

  • Any detector/pixel for which it’s standard deviation value satisfied the following conditions is considered to be “noisy”:

    • sdev(pixel) < mean(sdevs) - 3 * sdev(sdevs)

    • sdev(pixel) > mean(sdevs) + 3 * sdev(sdevs)

    • sdev(pixel) < mean(sdevs) * 0.0001

This is repeated three times from the second step.

This uses the Integration v1, Power v1 and Divide v1 algorithms for the first step.

The lower bound for the integration is currently fixed to 2000.

The upper bound for the integration is currently fixed to 19000.

Usage

Example: A clean run no noisy detectors

ws = CreateSampleWorkspace()
wsOut = IdentifyNoisyDetectors(ws)
print('The output workspace is filled with {:.0f} for good detectors and 0 for noisy ones.'.format( wsOut.readY(0)[0]))
sum = SumSpectra(wsOut)
print("{:.0f} good spectra left".format(sum.readY(0)[0]))

Output:

The output workspace is filled with 1 for good detectors and 0 for noisy ones.
200 good spectra left

Example: With lots of noisy detectors

ws = CreateSampleWorkspace()
#make the spectra vary
for i in range(ws.getNumberHistograms()):
    ws.setY(i, ws.readY(i)*2**i)

wsOut = IdentifyNoisyDetectors(ws)
sum = SumSpectra(wsOut)
print('{:.0f} good spectra left from an original {:.0f}.'.format(sum.readY(0)[0], wsOut.getNumberHistograms()))

Output:

15 good spectra left from an original 200.

Categories: AlgorithmIndex | Diagnostics

Source

C++ header: IdentifyNoisyDetectors.h

C++ source: IdentifyNoisyDetectors.cpp