burst_range¶

gwpy.astro.burst_range(psd: FrequencySeries, snr: float = 8, energy: float = 0.01, fmin: float = 100, fmax: float = 500) → Quantity[source]¶

Calculate the integrated GRB-like GW burst range from a strain PSD.

Parameters:¶

psdFrequencySeries: The instrumental power-spectral-density data.
snrfloat, optional: The signal-to-noise ratio for which to calculate range.
energyfloat, optional: The relative energy output of the GW burst.
fminfloat, optional: The lower frequency cutoff of the burst range integral.
fmaxfloat, optional: The upper frequency cutoff of the burst range integral.

Returns:¶

Examples

Grab some data for LIGO-Livingston around GW150914 and generate a PSD

>>> from gwpy.timeseries import TimeSeries
>>> hoft = TimeSeries.fetch_open_data("H1", 1126259446, 1126259478)
>>> hoff = hoft.psd(fftlength=4)

Now we can calculate the burst_range():

>>> from gwpy.astro import burst_range
>>> r = burst_range(hoff, fmin: float = 30)
>>> print(r)
42.5055584195 Mpc