burst_range¶

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

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

Parameters

psd : FrequencySeries

the instrumental power-spectral-density data

snr : float, optional

the signal-to-noise ratio for which to calculate range, default: 8

energy : float, optional

the relative energy output of the GW burst, defaults to 1e-2 for a GRB-like burst

fmin : float, optional

the lower frequency cutoff of the burst range integral, default: 100 Hz

fmax : float, optional

the upper frequency cutoff of the burst range integral, default: 500 Hz

Returns

range : Quantity

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=30)
>>> print(r)
42.5055584195 Mpc