Update VItamin and Bilby with Different Noises authored by Wei Changfeng's avatar Wei Changfeng
...@@ -89,7 +89,21 @@ Then, in response to the suspicion that there may be errors in my previous scrip ...@@ -89,7 +89,21 @@ Then, in response to the suspicion that there may be errors in my previous scrip
bilby.core.utils.setup_logger(outdir=outdir, label=label) bilby.core.utils.setup_logger(outdir=outdir, label=label)
``` ```
However this doesn't matter to the results. However this doesn't matter to the results.
(ii) Priors are used differently. Priors used by tutorial are:(only different parts are shown) (ii) Setting of start time. In the tutorial, duration is 4 second and start time is set as following:
```
ifos.set_strain_data_from_power_spectral_densities(
sampling_frequency=sampling_frequency, duration=duration,
start_time=ref_geocent_time-3)
```
Mine:
```
ifos.set_strain_data_from_power_spectral_densities(
sampling_frequency=sampling_frequency, duration=duration,
start_time=ref_geocent_time-duration/2.0)
```
I think, my way of setting start time is Ok.
(iii) Priors are used differently. Priors used by tutorial are:(only different parts are shown)
``` ```
priors['geocent_time'] = bilby.core.prior.Uniform( priors['geocent_time'] = bilby.core.prior.Uniform(
minimum=ref_geocent_time - 0.1, minimum=ref_geocent_time - 0.1,
...@@ -108,7 +122,15 @@ priors['geocent_time'] = bilby.core.prior.Uniform( ...@@ -108,7 +122,15 @@ priors['geocent_time'] = bilby.core.prior.Uniform(
maximum=ref_geocent_time + 1, maximum=ref_geocent_time + 1,
name='geocent_time', latex_label='$t_c$', unit='$s$') name='geocent_time', latex_label='$t_c$', unit='$s$')
``` ```
(iii) Result function. In this part, the additional specified priors are the correct format from Daniel. However, when I output priors, their distance priors are different:
>
'luminosity_distance': Uniform(minimum=1000, maximum=5000, name='luminosity_distance', latex_label='$d_L$', unit='$Mpc_{\\odot}$', boundary=None)
>
'luminosity_distance': UniformSourceFrame(minimum=100.0, maximum=5000.0, cosmology=FlatLambdaCDM(name="Planck15", H0=67.7 km / (Mpc s), Om0=0.307, Tcmb0=2.725 K, Neff=3.05, m_nu=[0. 0. 0.06] eV, Ob0=0.0486), name='luminosity_distance', latex_label='$d_L$', unit=Unit("Mpc"), boundary=None)
(iv) Result function.
Tutorial: Tutorial:
``` ```
result = bilby.run_sampler( result = bilby.run_sampler(
... ...
......