Inference of Single Waveforms

Gravelamps main function is performing Bayesian parameter estimation and model selection for the variety of implemented gravitational wave lensing models. One such configuration is that for individual gravitational wave events. This regime covers the microlensing and millilensing regimes and benefits from both the wave and geometric optics implementations. To perform the parameter estimation, Gravelamps wraps around the bilby and bilby_pipe frameworks.

Performing inference is handled by the gravelamps_inference entrypoint. This may typically be run as follows:

$ gravelamps_inference -o --option /path/to/ini-file.ini

Commandline Options

-h, --help

Shows the help message and exits.

-v, --verbose

Display the maximum level of information from the logger. Will override the INI setting.

-l, --local

If set, all jobs will be performed locally instead of being built for a scheduler. Will override the INI setting.

Example INI

The following example INI will generate the HTCondor submit files for a job using the isolated point mass model to analyse an unlensed injection.

[output]
outdir = inference-example
label = point-analysis-unlensed-injection
logging_level = INFO

[lens]
module = gravelamps.models.isolated_point
minimum_dimensionless_frequency = 0.01
maximum_dimensionless_frequency = 10000
length_dimensionless_frequency = 20000

minimum_source_position = 0.1
maximum_source_position = 3.0
length_source_position = 1000

[lens.injection]
module = None

[scheduler]
scheduler = condor

[scheduler.condor]
accounting_group = ligo.dev.o4.cbc.lensing.multi

[run]
local = False
includes_injection = True
no_injection_generation = True

[bilby]
detectors = H1
injection = True
zero-noise = True
waveform-generator = gravelamps.waveform.waveform_generator.LensedWaveformGenerator

[prior]
mass_1=36.0
mass_2=29.0
a_1=0.4
a_2=0.3
tilt_1=0.5
tilt_2=1.0
phi_12=1.7
phi_jl=0.3
luminosity_distance=2000.0
theta_jn=0.4
psi=2.659
phase=1.3
geocent_time=1126259642.413
ra=1.375
dec=-1.2108
lens_mass = Uniform(minimum=1, maximum=1000, name='lens_mass')
source_position = Uniform(minimum=0.1, maximum=1.0, name='source_position')
lens_fractional_distance = 0.5 

[injection.parameters]
mass_1=36.0
mass_2=29.0
a_1=0.4
a_2=0.3
tilt_1=0.5
tilt_2=1.0
phi_12=1.7
phi_jl=0.3
luminosity_distance=2000.0
theta_jn=0.4
psi=2.659
phase=1.3
geocent_time=1126259642.413
ra=1.375
dec=-1.2108

For a more detailed look at the options available in the INI file, please see the Gravelamps INI Structure page.