Observing Scenarios Injection Tool (bayestar-inject)

Rough-cut injection tool.

The idea is to efficiently sample events, uniformly in “sensitive volume” (differential comoving volume divided by 1 + z), and from a distribution of masses and spins, such that later detection cuts will not reject an excessive number of events.

This occurs in two steps. First, we divide the intrinsic parameter space into a very coarse 10x10x10x10 grid and calculate the maximum horizon distance in each grid cell. Second, we directly sample injections jointly from the mass and spin distribution and a uniform and isotropic spatial distribution with a redshift cutoff that is piecewise constant in the masses and spins.

usage: bayestar-inject [-h] [--seed SEED] [--version]
                       [-l CRITICAL|ERROR|WARNING|INFO|DEBUG|NOTSET]
                       [--cosmology {WMAP1,WMAP3,WMAP5,WMAP7,WMAP9,Planck13,Planck15,Planck18}]
                       (--distribution {bns_astro,bns_broad,nsbh_astro,nsbh_broad,bbh_astro,bbh_broad} | --distribution-samples DISTRIBUTION_SAMPLES)
                       --reference-psd PSD.xml[.gz] [--f-low F_LOW]
                       [--snr-threshold SNR_THRESHOLD]
                       [--min-triggers MIN_TRIGGERS] [--min-snr MIN_SNR]
                       [--max-distance Mpc] [--waveform WAVEFORM]
                       [--nsamples NSAMPLES] [-o INJ.xml[.gz]] [-j [JOBS]]

Named Arguments

--version

show program’s version number and exit

-l, --loglevel

Default: INFO

--cosmology

Possible choices: WMAP1, WMAP3, WMAP5, WMAP7, WMAP9, Planck13, Planck15, Planck18

Cosmological model

Default: “Planck15”

--distribution

Possible choices: bns_astro, bns_broad, nsbh_astro, nsbh_broad, bbh_astro, bbh_broad

Use a preset distribution

--distribution-samples

Load samples of the intrinsic mass and spin distribution from any file that can be read as an Astropy table. The table columns should be mass1, mass2, spin1z, and spin2z.

--reference-psd

PSD file

--f-low

Low frequency cutoff in Hz

Default: 25.0

--snr-threshold

Single-detector SNR threshold

Default: 4.0

--min-triggers

Emit coincidences only when at least this many triggers are found

Default: 2

--min-snr

Minimum decisive SNR of injections given the reference PSDs. Deprecated; use the synonymous –snr-threshold option instead.

--max-distance

Maximum luminosity distance for injections

--waveform

Waveform approximant

Default: “o2-uberbank”

--nsamples

Output this many injections

Default: 100000

-o, --output

Output file, optionally gzip-compressed

Default: -

-j, --jobs

Number of threads

Default: 1

random number generator options

Options that affect the Numpy pseudo-random number genrator

--seed

Pseudo-random number generator seed [default: initialized from /dev/urandom or clock]

Python helper functions

class ligo.skymap.tool.bayestar_inject.GWCosmo(cosmology)

Evaluate GW distance figures of merit for a given cosmology.

Parameters:

cosmoastropy.cosmology.FLRW

The cosmological model.