Data Accompanying: "Deep follow-up of GW151226: ordinary binary or low-mass-ratio system?"

Data accompanying the arxiv.org/abs/2203.13406 paper

This folder contain all the results for the Vajpeyi et al (2022) paper.
The structure of this repository is as follows:

├── GW150914
│   ├── literature_samples
│   │   └── mateu_qxeff_samples.npz
│   ├── results
│   │   ├── ptA_qxeff_deep_followup_result.json
│   │   └── ptB_qxeff_deep_followup_result.json
│   └── run_configs
│       ├── GWTC1_H1_PSD.dat
│       ├── GWTC1_L1_PSD.dat
│       ├── H1_data.gwf
│       ├── L1_data.gwf
│       ├── full.ini
│       ├── get_data.py
│       ├── ptA.ini
│       └── ptB.ini
├── GW151226
│   ├── literature_samples
│   │   ├── chia_qxeff_samples.npz
│   │   ├── mateu_qxeff_samples.npz
│   │   └── nitz_qxeff_samples.npz
│   ├── results
│   │   ├── chia_qxeff_deep_followup_result.json
│   │   ├── mateu_qxeff_deep_followup_result.json
│   │   └── nitz_qxeff_deep_followup_result.json
│   └── run_configs
│       ├── GWTC1_H1_PSD.dat
│       ├── GWTC1_L1_PSD.dat
│       ├── H1_data.gwf
│       ├── L1_data.gwf
│       ├── full.ini
│       ├── get_data.py
│       ├── ptChia.ini
│       ├── ptMateu.ini
│       └── ptNitz.ini
├── Injection
│   ├── results
│   │   ├── full_pe_on_injection_result.json
│   │   ├── ptA_qxeff_deep_followup_result.json
│   │   └── ptB_qxeff_deep_followup_result.json
│   └── run_configs
│       ├── full_pe.ini
│       ├── full_pe.prior
│       ├── ptA.ini
│       ├── ptA.prior
│       ├── ptB.ini
│       └── ptB.prior
├── README.md
├── print_odds.py
└── requirements.txt

The literature_samples folder contains the q-xeff samples from the literature.
The results folder contains the deep followup results for the various case studies.
The run_configs folder contains the configuration files analyses.

To reproduce the values from the Vajpeyi et al (2022) paper tables, run the following commands:

pip install -r requirements.txt
python print_odds.py

This will print the following: 

-----------------------GW151226 Odds-----------------------
-----------------------  -------------  -------------  ------------  ----------------
Pts                      KDE Post Odds  Prior Odds     Bayes Factor  'Deep' Post Odds
KDE-Mateu Mateu VS Chia  3.45+/-0.14    1.677+/-0.013  1.00+/-0.21   1.7+/-0.4
KDE-Chia Mateu VS Chia   0.328+/-0.004  1.676+/-0.012  1.00+/-0.21   1.7+/-0.4
KDE-Nitz Mateu VS Chia   1.715+/-0.014  1.676+/-0.012  1.00+/-0.21   1.7+/-0.4
-----------------------  -------------  -------------  ------------  ----------------
-----------------------  -------------  -------------  ------------  ----------------
Pts                      KDE Post Odds  Prior Odds     Bayes Factor  'Deep' Post Odds
KDE-Mateu Mateu VS Nitz  1.46+/-0.04    1.076+/-0.011  1.31+/-0.21   1.42+/-0.29
KDE-Chia Mateu VS Nitz   0.927+/-0.011  1.076+/-0.010  1.31+/-0.21   1.42+/-0.29
KDE-Nitz Mateu VS Nitz   0.870+/-0.006  1.076+/-0.012  1.31+/-0.21   1.41+/-0.29
-----------------------  -------------  -------------  ------------  ----------------
----------------------  -------------  -------------  ------------  ----------------
Pts                     KDE Post Odds  Prior Odds     Bayes Factor  'Deep' Post Odds
KDE-Mateu Nitz VS Chia  2.38+/-0.10    1.560+/-0.012  0.76+/-0.21   1.19+/-0.25
KDE-Chia Nitz VS Chia   0.355+/-0.004  1.556+/-0.011  0.76+/-0.21   1.18+/-0.25
KDE-Nitz Nitz VS Chia   1.973+/-0.014  1.556+/-0.011  0.76+/-0.21   1.18+/-0.25
----------------------  -------------  -------------  ------------  ----------------



------------------------GW150914 A VS B------------------------
PtA q-xeff: (0.7, 0.02)
PtB q-xeff: (0.78, -0.04)
---------------  -------------  -------------  ------------  ----------------
Pts              KDE Post Odds  Prior Odds     Bayes Factor  'Deep' Post Odds
GW150914 A VS B  0.248+/-0.016  1.056+/-0.013  0.2+/-0.4     0.24+/-0.10
---------------  -------------  -------------  ------------  ----------------



------------------------Injection A VS B------------------------
PtA q-xeff: (0.3605278878376452, 0.1729095104595008)
PtB q-xeff: (0.3923, 0.2461)
----------------  -------------  -------------  ------------  ----------------
Pts               KDE Post Odds  Prior Odds     Bayes Factor  'Deep' Post Odds
Injection A VS B  0.474+/-0.018  1.344+/-0.010  0.35+/-0.24   0.47+/-0.11
----------------  -------------  -------------  ------------  ----------------

To rerun the deep followup analyses, run the following commands:

deep_followup_setup <ini_file> 
parallel_bilby_analysis <args printed by deep_followup_setup>

External Data

This work utilised some external data. The data is described below:


1. q-xeff Posteriors
* The q-xeff samples labelled as 'Mateu' are obtained from Mateu-Lucena et al (2022)
* The q-xeff samples labelled as 'Nitz' are obtained from Nitz et al (2021)
* The q-xeff samples labelled as 'Chia' are obtained from the authors of Chia et al (2022), via email correspondence

2. Gravitational Wave Strain Data
* The raw strain data from the LIGO Hanford and Livingston detectors are obtained from GWOSC. The get_data.py script in the run_configs folder demonstrates how the data is downloaded.

3. Gravitational Wave PSD Data
* The PSDs for the GW150914 and GW151226 analyses are obtained from the LIGO dcc

References

* Nitz et al (2021): https://iopscience.iop.org/article/10.3847/1538-4357/ac1c03
* Mateu-Lucena et al (2022): https://doi.org/10.1093/mnras/stac2724
* Chia et al (2022): https://journals.aps.org/prd/abstract/10.1103/PhysRevD.106.024009
* LIGO dcc: https://dcc-llo.ligo.org/LIGO-P1900011/public