O4-Review.md

@@ -221,7 +221,85 @@ Ran the above [Fiducial 4s BBH Inj](#fiducial-bbh-inj) with nCPU 32, 64, 128.
 
 ## 4s BBH PP Test
 
+<details><summary>Configs Contents</summary>
+
+[Link to config](https://git.ligo.org/avi.vajpeyi/pbilby_o4_review/-/blob/main/pp_tests/bbh_4s/4s_bbh_0.ini)
+
+```ini
+prior-file = 4s_bbh.prior
+injection-file = injections.dat
+gaussian-noise = True
+maximum-frequency = 1024
+minimum-frequency = 20
+trigger-time = 0
+injection = True
+detectors = [H1, L1]
+psd_dict = None
+duration = 4
+distance-marginalization = True
+phase-marginalization = False
+time-marginalization = True
+jitter-time = True
+reference-frame = H1L1
+time-reference = geocent
+waveform-approximant = IMRPhenomXPHM
+catch-waveform-errors = True
+enforce-signal-duration = False
+nlive = 1000
+nact = 50
 
+```
+
+Prior
+```python
+chirp_mass = bilby.gw.prior.UniformInComponentsChirpMass(minimum=25, maximum=50, name='chirp_mass', latex_label='$\\mathcal{M}$', unit=None, boundary=None)
+mass_ratio = bilby.gw.prior.UniformInComponentsMassRatio(minimum=0.05, maximum=1, name='mass_ratio', latex_label='$q$', unit=None, boundary=None)
+mass_1 = Constraint(minimum=1, maximum=1000, name=None, latex_label=None, unit=None)
+mass_2 = Constraint(minimum=1, maximum=1000, name=None, latex_label=None, unit=None)
+# total_mass = Constraint(minimum=150, maximum=1000)
+a_1 = Uniform(minimum=0, maximum=0.99, name=None, latex_label=None, unit=None, boundary=None)
+a_2 = Uniform(minimum=0, maximum=0.99, name=None, latex_label=None, 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 = bilby.gw.prior.UniformSourceFrame(minimum=1.0, maximum=1000.0, name='luminosity_distance', latex_label='$d_L$', unit="Mpc", boundary=None)
+dec = Cosine(minimum=-1.5707963267948966, maximum=1.5707963267948966, name='dec', latex_label='$\\mathrm{DEC}$', unit=None, boundary=None)
+ra = Uniform(minimum=0, maximum=6.283185307179586, name='ra', latex_label='$\\mathrm{RA}$', unit=None, boundary='periodic')
+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')
+geocent_time = Uniform(minimum=-0.1, maximum=0.1, name=None, latex_label=None, unit=None, boundary=None)
+
+```
+
+</details>
+
+
+
+<details><summary>KS-test p-value</summary>
+
+```
+tilt_1: 0.2110473235496638
+ra: 0.4553802594360755
+a_1: 0.16749876930052376
+phi_12: 0.9312673273519888
+tilt_2: 0.6359606676409679
+chirp_mass: 0.5084285364847729
+phi_jl: 0.7645417984623673
+phase: 0.3009719585091924
+geocent_time: 0.5782924423635027
+theta_jn: 0.3038257852818535
+mass_ratio: 0.8624046637946167
+psi: 0.8546023833772283
+dec: 0.957635050095067
+a_2: 0.4764282474559757
+luminosity_distance: 0.5814225969546702
+Combined p-value: 0.9074045216633773
+```
+
+
+</details>