Skip to content

Changes

Page history
Update O4 Review authored by Avi Vajpeyi's avatar Avi Vajpeyi
Show whitespace changes
Inline Side-by-side
O4-Review.md
View page @ 29e3848a
...@@ -180,7 +180,236 @@ Note: the JS-Divergence values are median +- std values obtained with bootstrapp ...@@ -180,7 +180,236 @@ Note: the JS-Divergence values are median +- std values obtained with bootstrapp
## GW191109 ## GW191109
<details><summary>Configs</summary>
https://git.ligo.org/avi.vajpeyi/pbilby_o4_review/-/tree/main/GW191109
```python
chirp_mass = bilby.gw.prior.UniformInComponentsChirpMass(minimum=26.700411297080883, maximum=96.67976382167427, name='chirp_mass', latex_label='$\\mathcal{M}$', unit='$M_{\\odot}$', boundary=None)
mass_ratio = bilby.gw.prior.UniformInComponentsMassRatio(minimum=0.05, maximum=1.0, name='mass_ratio', latex_label='$q$', unit=None, boundary=None)
mass_1 = Constraint(minimum=1, maximum=1000, name='mass_1', latex_label='$m_1$', unit=None)
mass_2 = Constraint(minimum=1, maximum=1000, name='mass_2', latex_label='$m_2$', unit=None)
a_1 = Uniform(minimum=0, maximum=0.99, name='a_1', latex_label='$a_1$', unit=None, boundary=None)
a_2 = Uniform(minimum=0, maximum=0.99, name='a_2', latex_label='$a_2$', unit=None, boundary=None)
tilt_1 = Sine(minimum=0, maximum=3.141592653589793, name='tilt_1', latex_label='$\\theta_1$', unit=None, boundary=None)
tilt_2 = Sine(minimum=0, maximum=3.141592653589793, name='tilt_2', latex_label='$\\theta_2$', unit=None, boundary=None)
phi_12 = Uniform(minimum=0, maximum=6.283185307179586, name='phi_12', latex_label='$\\Delta\\phi$', unit=None, boundary='periodic')
phi_jl = Uniform(minimum=0, maximum=6.283185307179586, name='phi_jl', latex_label='$\\phi_{JL}$', unit=None, boundary='periodic')
luminosity_distance = PowerLaw(alpha=2, minimum=100, maximum=10000, name='luminosity_distance', latex_label='$d_L$', unit='Mpc', boundary=None)
theta_jn = Sine(minimum=0, maximum=3.141592653589793, name='theta_jn', latex_label='$\\theta_{JN}$', unit=None, boundary=None)
psi = Uniform(minimum=0, maximum=3.141592653589793, name='psi', latex_label='$\\psi$', unit=None, boundary='periodic')
phase = Uniform(minimum=0, maximum=6.283185307179586, name='phase', latex_label='$\\phi$', unit=None, boundary='periodic')
azimuth = Uniform(minimum=0, maximum=6.283185307179586, name=None, latex_label='$\\epsilon$', unit=None, boundary='periodic')
zenith = Sine(minimum=0, maximum=3.141592653589793, name=None, latex_label='$\\kappa$', unit=None, boundary=None)
geocent_time = Uniform(minimum=1257296855.1164582, maximum=1257296855.316458, name='geocent_time', latex_label='$t_c$', unit='$s$', boundary=None)
recalib_H1_amplitude_0 = Gaussian(mu=-0.009695320217683778, sigma=0.028499697803266222, name='recalib_H1_amplitude_0', latex_label='$A^H1_0$', unit=None, boundary='reflective')
...
recalib_L1_frequency_9 = DeltaFunction(peak=447.9999999999999, name='recalib_L1_frequency_9', latex_label='$f^L1_9$', unit=None)
time_jitter = Uniform(minimum=-0.0009765625, maximum=0.0009765625, name=None, latex_label=None, unit=None, boundary='periodic')
```
```
################################################################################
## Calibration arguments
################################################################################
calibration-model=CubicSpline
spline-calibration-envelope-dict={H1:calibration/H1.dat,L1:calibration/L1.dat}
spline-calibration-nodes=10
spline-calibration-amplitude-uncertainty-dict=None
spline-calibration-phase-uncertainty-dict=None
calibration-prior-boundary=None
################################################################################
## Data generation arguments
################################################################################
ignore-gwpy-data-quality-check=True
gps-tuple=None
gps-file=None
timeslide-file=None
timeslide-dict=None
trigger-time=1257296855.216458
n-simulation=0
data-dict=None
data-format=None
allow-tape=True
channel-dict={H1:DCS-CALIB_STRAIN_CLEAN_SUB60HZ_C01_T1700406_v4, L1:DCS-CALIB_STRAIN_CLEAN_SUB60HZ_C01_T1700406_v4}
gaussian-noise=False
zero-noise=False
################################################################################
## Detector arguments
################################################################################
coherence-test=False
detectors=['H1', 'L1']
duration=4.0
generation-seed=None
psd-dict={H1:psds/H1-psd.dat, L1:psds/L1-psd.dat}
psd-fractional-overlap=0.5
post-trigger-duration=2.0
sampling-frequency=1024.0
psd-length=4
psd-maximum-duration=1024
psd-method=median
psd-start-time=None
maximum-frequency={ 'H1': 448, 'L1': 448, }
minimum-frequency={ 'H1': 20, 'L1': 20, }
tukey-roll-off=0.4
resampling-method=lal
################################################################################
## Injection arguments
################################################################################
injection=False
injection-dict=None
injection-file=None
injection-numbers=[None]
injection-waveform-approximant=None
injection-waveform-arguments=None
################################################################################
## Job submission arguments
################################################################################
accounting=ligo.dev.o4.cbc.pe.bilby
accounting-user=None
label=GW191109_dynesty
local=False
local-generation=True
local-plot=False
outdir=outdir_dynesty_A
overwrite-outdir=False
periodic-restart-time=28800
request-disk=5.0
request-memory=8.0
request-memory-generation=None
request-cpus=16
conda-env=None
scheduler=condor
scheduler-args=None
scheduler-module=[None]
scheduler-env=None
scheduler-analysis-time=7-00:00:00
submit=False
condor-job-priority=0
transfer-files=False
additional-transfer-paths=[None]
log-directory=None
osg=False
desired-sites=None
analysis-executable=None
analysis-executable-parser=None
################################################################################
## Likelihood arguments
################################################################################
calibration-marginalization=False
distance-marginalization=True
distance-marginalization-lookup-table=/home/gregory.ashton/public_html/bilby_review/bilby_o4_review_230228/GW191109/TD.npz
phase-marginalization=False
time-marginalization=True
jitter-time=True
reference-frame=H1L1
time-reference=geocent
likelihood-type=GravitationalWaveTransient
calibration-lookup-table=None
roq-folder=None
roq-linear-matrix=None
roq-quadratic-matrix=None
roq-weights=None
roq-weight-format=None
roq-scale-factor=1.0
fiducial-parameters=None
update-fiducial-parameters=False
epsilon=0.025
extra-likelihood-kwargs=None
################################################################################
## Output arguments
################################################################################
plot-trace=False
plot-data=False
plot-injection=False
plot-spectrogram=False
plot-calibration=False
plot-corner=False
plot-marginal=False
plot-skymap=False
plot-waveform=False
plot-format=png
create-summary=False
email=None
notification=Never
queue=None
existing-dir=None
webdir=None
summarypages-arguments=None
result-format=hdf5
final-result=True
final-result-nsamples=20000
################################################################################
## Prior arguments
################################################################################
default-prior=BBHPriorDict
deltaT=0.2
prior-file=/home/gregory.ashton/public_html/bilby_review/bilby_o4_review_230228/GW191109/Prod5.prior
prior-dict=None
enforce-signal-duration=True
################################################################################
## Post processing arguments
################################################################################
postprocessing-executable=None
postprocessing-arguments=None
single-postprocessing-executable=None
single-postprocessing-arguments=None
################################################################################
## Sampler arguments
################################################################################
sampler=bilby_mcmc
sampling-seed=None
n-parallel=4
sampler-kwargs={'nlive': 1600, 'sample': 'rwalk', 'walks': 100, 'nact': 50, 'check_point_delta_t':1800, 'check_point_plot':True}
reweighting-configuration=None
reweight-nested-samples=True
################################################################################
## Waveform arguments
################################################################################
waveform-generator=bilby.gw.waveform_generator.WaveformGenerator
reference-frequency=20.0
waveform-approximant=IMRPhenomXPHM
catch-waveform-errors=True
pn-spin-order=-1
pn-tidal-order=-1
pn-phase-order=-1
pn-amplitude-order=4
numerical-relativity-file=None
waveform-arguments-dict=None
mode-array=[None]
frequency-domain-source-model=lal_binary_black_hole
conversion-function=None
generation-function=None
```
</details>
<details><summary>Results</summary> <details><summary>Results</summary>
| | JS Divergence | | | JS Divergence |
|:--------------------|:--------------------| |:--------------------|:--------------------|
| mass_2 | 0.00911 +/- 0.00030 | | mass_2 | 0.00911 +/- 0.00030 |
......