NuRadioReco.modules.LOFAR.stationPulseFinder module

NuRadioReco.modules.LOFAR.stationPulseFinder.find_snr_of_timeseries(timeseries, window_start=0, window_end=-1, noise_start=0, noise_end=-1)

Return the signal-to-noise ratio (SNR) of a given time trace. The additional parameters allow to select a window where the signal is and a separate window where there is only noise.

Parameters:

timeseries: list or np.ndarray

The time trace

window_startint

window_endint

We look for the peak inside the array timeseries[window_start:window_end]

noise_startint

noise_endint

The array timeseries[noise_start:noise_end] is used to calculate the noise level

Returns:

The SNR of the time trace

class NuRadioReco.modules.LOFAR.stationPulseFinder.stationPulseFinder

Bases: object

Look for significant pulses in every station. The module uses beamforming to enhance the sensitivity in direction estimated from the LORA particle data. It also identifies the channels which have an SNR good enough to use for direction fitting later.

Methods

begin([window, noise_window, cr_snr, ...])	Sets the window size to use for pulse finding, as well as the number of samples away from the pulse to use for noise measurements.
run(event, detector)	Run a beamformer on all stations in event to search for significant pulses.

end

begin(window=500, noise_window=10000, cr_snr=3, good_channels=6, logger_level=30)

Sets the window size to use for pulse finding, as well as the number of samples away from the pulse to use for noise measurements. The function also defines what an acceptable SNR is to consider a cosmic-ray signal to be in the trace, as well as the number of good channels a station should have to be kept for further processing.

Parameters:

windowint, default=500

Size of the window to look for pulse

noise_windowint, default=10000

The trace used for noise characterisation goes from sample 0 to the start of the pulse searching window minus this number.

cr_snrfloat, default=3

The minimum SNR a channel should have to be considered having a CR signal.

good_channelsint, default=6

The minimum number of good channels a station should have in order be “triggered”.

logger_levelint, default=logging.WARNING

The logging level to use for the module.

run(event, detector)

Run a beamformer on all stations in event to search for significant pulses.

Parameters:

eventEvent object

The event on which to apply the pulse finding.

detectorDetector object

The detector related to the event.

end()