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:
rangeQuantity

The GRB-like-burst sensitive range [Mpc (default)].

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