Skip to content

implement a new sampling parameters: mu_2

Soichiro Morisaki requested to merge soichiro/bilby:implement_mu2 into master

This MR implements a new sampling parameter \mu_2, a combination of Post-Newtonian coefficients introduced by https://arxiv.org/abs/2203.05216. This can be used instead of \chi_1.

In this implementation, the prior is decomposed as \pi(\mathcal{M}, q, \mu_2, \chi_2) = \pi (\mu_2|\mathcal{M}, q, \chi_2) \pi(\mathcal{M}, q, \chi_2), and \pi (\mu_2|\mathcal{M}, q, \chi_2) is implemented by using ConditionalPrior. Since \mu_2 is a linear combination of \chi_1, it is reduced to \chi_1 prior with linear transformation from \mu_2 to \chi_1.

This MR does the following:

  • bilby/gw/conversion.py: implement conversion function from \mu_2 to \chi_1
  • bilby/gw/prior.py: implement conditional prior of \mu_2, and conditional prior of \chi_{1,\perp} conditioned on \mu_2

Here are auto-correlation plots from bilby-mcmc: without mu_2 and with mu_2, which were generated by this script.

This will be prototype of !1162, as this does not implement \mu_1.

Edited by Soichiro Morisaki

Merge request reports