\(\renewcommand\AA{\unicode{x212B}}\)
Table of Contents
Filters out events associated with pulses that happen when proton charge is lower than a given percentage of the average.
Name | Direction | Type | Default | Description |
---|---|---|---|---|
InputWorkspace | Input | EventWorkspace | Mandatory | An event workspace |
OutputWorkspace | Output | EventWorkspace | Mandatory | The name to use for the output workspace |
LowerCutoff | Input | number | 95 | The percentage of the average to use as the lower bound |
The algorithm looks at sample logs (“proton_charge”), finds the mean, and rejects any events that occurred during a pulse that was below a certain percentage of that mean. This effectively removes neutrons from the background that were measured while the accelerator was not actually producing neutrons, reducing background noise.
Example - Using a simple proton charge log
ws = CreateSampleWorkspace("Event",BankPixelWidth=1)
AddTimeSeriesLog(ws, Name="proton_charge", Time="2010-01-01T00:00:00", Value=100)
AddTimeSeriesLog(ws, Name="proton_charge", Time="2010-01-01T00:10:00", Value=100)
AddTimeSeriesLog(ws, Name="proton_charge", Time="2010-01-01T00:20:00", Value=100)
AddTimeSeriesLog(ws, Name="proton_charge", Time="2010-01-01T00:30:00", Value=100)
AddTimeSeriesLog(ws, Name="proton_charge", Time="2010-01-01T00:40:00", Value=15)
AddTimeSeriesLog(ws, Name="proton_charge", Time="2010-01-01T00:50:00", Value=100)
AddSampleLog(ws,"gd_prtn_chrg", "1e6", "Number")
wsFiltered = FilterBadPulses(ws)
print("The number of events that remain: %i" % wsFiltered.getNumberEvents())
print("compared to the number in the unfiltered workspace: %i" % ws.getNumberEvents())
Output:
The number of events that remain: 950
compared to the number in the unfiltered workspace: 1900
Categories: AlgorithmIndex | Events\EventFiltering
C++ header: FilterBadPulses.h (last modified: 2021-03-31)
C++ source: FilterBadPulses.cpp (last modified: 2021-03-31)