\(\renewcommand\AA{\unicode{x212B}}\)
CreatePSDBleedMask dialog.
Table of Contents
Runs a diagnostic test for saturation of PSD tubes and creates a MaskWorkspace marking the failed tube spectra.
| Name | Direction | Type | Default | Description |
|---|---|---|---|---|
| InputWorkspace | Input | MatrixWorkspace | Mandatory | The name of the input workspace. |
| OutputWorkspace | Output | MatrixWorkspace | Mandatory | The name of the output MaskWorkspace which will contain the result masks. |
| MaxTubeFramerate | Input | number | Mandatory | The maximum rate allowed for a tube in counts/us/frame. |
| NIgnoredCentralPixels | Input | number | 80 | The number of pixels about the centre to ignore. |
| NumberOfFailures | Output | number | An output property containing the number of masked tubes |
The algorithm identifies tubes, which become saturated due to high neutron flux from intense Bragg peaks affecting these tubes, so these tubes are not counting neutrons but show constant high counts over the tube in the elastic line(s) energy region(s). After that, the algorithm generates masking workspace to eliminate such tubes from the final results.
The bleeding effect occurs due to some old data acquisition electronics and observed as homogeneous high counts reading, symmetric over the whole tube length. The bleeding occurs at or close to incident energy or specific energies in multirep mode, where bright Bragg reflections hit the tube in positions, symmetric with respect to the tube centre.
The example of the workspace, affected by the bleed current and requesting bleeding corrections is presented on the following picture:
First of all, the algorithm attempts to find all tubes within the instrument attached to the workspace. If successful, each tube is tested for saturation above the level defined by the production of MaxTubeFramerate property by the goodfrm log value, retrieved from the workspace.
If the total signal, summed over all tubes pixels, excluding the number of central pixels specified by the NIgnoredCentralPixels property, exceeds the threshold, specified by the MaxTubeFramerate multiplied by the goodfrm value, the tube becomes masked.
The following image shows two instrument views obtained for two converted to energy transfer workspaces, recorded for two incident energies in multirep mode. Second energy range suffers from bleeding signal as on the picture above, so proper bleeding corrections are calculated. First energy range does not show bleeding signal, so no bleeding corrections applied to it.
Example:
import numpy as np
ws=CreateSampleWorkspace()
AddSampleLog(ws,"goodfrm","10","Number")
noisyY = np.array(ws.readY(0))
noisyY[0]=1e20
ws.setY(50,noisyY)
(wsOut, numFailures) = CreatePSDBleedMask(ws,MaxTubeFramerate=10, NIgnoredCentralPixels=2)
print("{} spectra have been masked in wsOut".format(numFailures))
Output:
10 spectra have been masked in wsOut
Categories: AlgorithmIndex | Diagnostics