pp/dynesty.md

@@ -55,3 +55,56 @@ phase =  Uniform(name='phase', minimum=0, maximum=2 * np.pi, periodic_boundary=T
 ![image](uploads/c86ad17839756e84b804f98db07c836e/image.png)
 ![image](uploads/c7ded58e8db76251cce48bae297c192a/image.png)
 ![image](uploads/4ab487489ea384d4f80fdcdaf85f2f3b/image.png)
+
+### runD
+
+```python
+# config.ini
+label = dynesty_runA
+outdir = .
+accounting = ligo.dev.o3.cbc.pe.lalinference
+
+detectors = [H1, L1]
+duration = 16
+
+trigger-time = 0
+deltaT = 0.2
+prior-file = chirp_mass_mass_ratio_spins.prior
+
+sampler = dynesty
+sampler-kwargs = {nlive: 1000, walks: 100, n_check_point: 5000}
+
+create-plots = False
+injection = True
+n-injection = 100
+
+time-marginalization = True
+distance-marginalization = True
+phase-marginalization = True
+```
+
+```
+# chirp_mass_mass_ratio_spins.prior
+chirp_mass = Uniform(name='chirp_mass', minimum=6, maximum=100, unit='$M_{\\odot}$')
+mass_ratio =  Uniform(name='mass_ratio', minimum=0.125, maximum=1)
+mass_1 = Constraint(name='mass_1', minimum=5, maximum=100)
+mass_2 = Constraint(name='mass_2', minimum=5, maximum=100)
+a_1 =  Uniform(name='a_1', minimum=0, maximum=0.99)
+a_2 =  Uniform(name='a_2', minimum=0, maximum=0.99)
+tilt_1 =  Sine(name='tilt_1')
+tilt_2 =  Sine(name='tilt_2')
+phi_12 =  Uniform(name='phi_12', minimum=0, maximum=2 * np.pi, periodic_boundary=True)
+phi_jl =  Uniform(name='phi_jl', minimum=0, maximum=2 * np.pi, periodic_boundary=True)
+ra =  Uniform(name='ra', minimum=0, maximum=2 * np.pi, periodic_boundary=True)
+dec =  Cosine(name='dec')
+luminosity_distance = bilby.gw.prior.UniformComovingVolume(name='luminosity_distance', minimum=1e2, maximum=5e3, unit='Mpc', periodic_boundary=False)
+theta_jn = Sine(name='iota')
+geocent_time = Uniform(name='geocent_time', minimum=-0.1, maximum=0.1)
+psi =  Uniform(name='psi', minimum=0, maximum=np.pi, periodic_boundary=True)
+phase =  Uniform(name='phase', minimum=0, maximum=2 * np.pi, periodic_boundary=True)
+```
+
+![image](uploads/b685bad143d84d89c66a86884ad6ad01/image.png)
+![image](uploads/bea266d1eece5cc4084c53ce1298b243/image.png)
+
+I haven't added the sampling times as they where skewed by evictions and the fact that we don't track sampling times across evictions.
\ No newline at end of file