DetectorEfficiencyVariation v1

../_images/DetectorEfficiencyVariation-v1_dlg.png

DetectorEfficiencyVariation dialog.

Summary

Compares two white beam vanadium workspaces from the same instrument to find detectors whose efficiencies have changed beyond a threshold.

Properties

Name Direction Type Default Description
WhiteBeamBase Input MatrixWorkspace Mandatory Name of a white beam vanadium workspace
WhiteBeamCompare Input MatrixWorkspace Mandatory Name of a matching second white beam vanadium run from the same instrument
OutputWorkspace Output MatrixWorkspace Mandatory A MaskWorkpace where each spectra that failed the test is masked. Each histogram from the input workspace maps to a histogram in this workspace with one value that indicates if there was a dead detector.
Variation Input number 1.1 Identify histograms whose total number of counts has changed by more than this factor of the median change between the two input workspaces.
StartWorkspaceIndex Input number 0 The index number of the first spectrum to include in the calculation (default: 0)
EndWorkspaceIndex Input number Optional The index number of the last spectrum to include in the calculation (default: the last spectrum in the workspace)
RangeLower Input number Optional No bin with a boundary at an x value less than this will be included in the summation used to decide if a detector is ‘bad’ (default: the start of each histogram)
RangeUpper Input number Optional No bin with a boundary at an x value higher than this value will be included in the summation used to decide if a detector is ‘bad’ (default: the end of each histogram)
NumberOfFailures Output number    

Description

It is intended that the input white beam vanadium workspaces are from the same instrument and were collected before and after an experimental run of interest. First the ratios of the total number of counts in corresponding histograms from each input workspace are calculated and then the median ratio is calculated. Each ratio is compared to the median and a histogram will fail when any of the following conditions are true:

  • (sum1/sum2)/median(sum1/sum2) > Variation
  • (sum1/sum2)/median(sum1/sum2) < 1/Variation

where sum1 is the sum of the counts in a histogram in the workspace WhiteBeamBase and sum2 is the sum of the counts in the equivalent histogram in WhiteBeamCompare. The above equations only make sense for identifying bad detectors if Variation > 1. If a value of less than one is given for Variation then Variation will be set to the reciprocal.

The output workspace contains a MaskWorkspace where those spectra that fail the tests are masked and those that pass them are assigned a single positive value.

Child algorithms used

Uses the Integration v1 algorithm to sum the spectra.

Usage

Example:

import numpy as np

wsBase=CreateSampleWorkspace()
ws=CreateSampleWorkspace()
noisyY =  np.array(ws.readY(0))
noisyY[0]=1e20
ws.setY(50,noisyY)

(wsOut, numFailures) = DetectorEfficiencyVariation(WhiteBeamBase='wsBase', WhiteBeamCompare='ws')

print("{} spectra has been masked in wsOut".format(numFailures))

Output:

1 spectra has been masked in wsOut

Categories: Algorithms | Diagnostics

Source

C++ source: DetectorEfficiencyVariation.cpp (last modified: 2016-10-12)

C++ header: DetectorEfficiencyVariation.h (last modified: 2016-06-07)