pp/dynesty.md

 # Dynesty
 
-## Run 1
+All runs can be found [on cit](https://ldas-jobs.ligo.caltech.edu/~gregory.ashton/bilby_tests/dynesty/). 
+
+Runs are generated by drawing simulated signals in Gaussian (aLIGO noise curve) data for H1 and L1. The prior distribution from which simulation values are selected is
+
+```python
+chirp_mass = Uniform(name='chirp_mass', minimum=25, maximum=100, unit='$M_{\\odot}$')
+mass_ratio =  Uniform(name='mass_ratio', minimum=0.125, maximum=1)
+a_1 =  0
+a_2 =  0
+tilt_1 =  0
+tilt_2 =  0
+phi_12 =  0
+phi_jl =  0
+dec =  Cosine(name='dec')
+ra =  Uniform(name='ra', minimum=0, maximum=2 * np.pi)
+iota =  Sine(name='iota')
+psi =  Uniform(name='psi', minimum=0, maximum=np.pi)
+phase =  Uniform(name='phase', minimum=0, maximum=2 * np.pi)
+luminosity_distance = PowerLaw(alpha=2, minimum=100, maximum=5000, name='luminosity_distance')
+geocent_time = 0
+```
+
+Then, for the analysis the same prior is used **except** that the `geocent_time` prior is set to
+
+```
+geocent_time = Uniform(minimum=-0.05, maximum=0.05, name='geocent_time', latex_label='$t_c$', unit='$s$')
+```
+
+in all cases we use the `bilby.gw.likelihood.GravitationalWaveTransientLikelihood` with no marginalisation.
+
+## Run 1: `nlive=2000` and `walks=50`
+
+Here is the full set of sampler kwargs:
+
+```python
+{'bound': 'multi', 'sample': 'rwalk', 'verbose': True, 'check_point_delta_t': 600, 'nlive': 2000, 'first_update': None, 'npdim': None, 'rstate': None, 'queue_size': None, 'pool': None, 'use_pool': None, 'live_points': None, 'logl_args': None, 'logl_kwargs': None, 'ptform_args': None, 'ptform_kwargs': None, 'enlarge': None, 'bootstrap': None, 'vol_dec': 0.5, 'vol_check': 2.0, 'facc': 0.5, 'slices': 5, 'walks': 50, 'update_interval': 1200, 'print_func': <bound method Dynesty._print_func of <bilby.core.sampler.dynesty.Dynesty object at 0x2b3c0f61c438>>, 'dlogz': 0.1, 'maxiter': None, 'maxcall': None, 'logl_max': inf, 'add_live': True, 'print_progress': True, 'save_bounds': True}
+```
+
+### PP plot
+
+![image](uploads/354c3d9e27257a476f2f8cee838bb1da/image.png)
+
+### Sampling time
+
+![image](uploads/11b36aeaa61d598be5ff8b122862a36c/image.png)
+
+
+## Old runs
+
+The following are some old results, leaving here in case we want to come back to them. 
+
+### Run 1
 
 Using this kwargs argument for `dynesty`
 
@@ -19,7 +70,7 @@ Produced this PP plot:
 
 ![image](uploads/0cd001c8e5cabf955a2e94a296814e62/image.png)
 
-## Run 2
+### Run 2
 
 Running again with more walks