Skip to content
Snippets Groups Projects

Enable posterior plots of added cbc parameters

Merged Paul Lasky requested to merge enable_posterior_plots_of_added_CBC_parameters into master
1 file
+ 13
0
Compare changes
  • Side-by-side
  • Inline
+ 13
0
@@ -370,16 +370,29 @@ class Result(dict):
""" Construct widely used derived parameters of CBCs """
self.posterior['mass_chirp'] = (self.posterior.mass_1 * self.posterior.mass_2) ** 0.6 / (
self.posterior.mass_1 + self.posterior.mass_2) ** 0.2
self.search_parameter_keys.append('mass_chirp')
self.parameter_labels.append('$\mathcal{M}$')
self.posterior['q'] = self.posterior.mass_2 / self.posterior.mass_1
self.search_parameter_keys.append('q')
self.parameter_labels.append('$q$')
self.posterior['eta'] = (self.posterior.mass_1 * self.posterior.mass_2) / (
self.posterior.mass_1 + self.posterior.mass_2) ** 2
self.search_parameter_keys.append('eta')
self.parameter_labels.append('$\eta$')
self.posterior['chi_eff'] = (self.posterior.a_1 * np.cos(self.posterior.tilt_1)
+ self.posterior.q * self.posterior.a_2 * np.cos(self.posterior.tilt_2)) / (
1 + self.posterior.q)
self.search_parameter_keys.append('chi_eff')
self.parameter_labels.append('$\chi_{\mathrm eff}$')
self.posterior['chi_p'] = np.maximum(self.posterior.a_1 * np.sin(self.posterior.tilt_1),
(4 * self.posterior.q + 3) / (3 * self.posterior.q + 4) * self.posterior.q
* self.posterior.a_2 * np.sin(self.posterior.tilt_2))
self.search_parameter_keys.append('chi_p')
self.parameter_labels.append('$\chi_{\mathrm p}$')
def check_attribute_match_to_other_object(self, name, other_object):
""" Check attribute name exists in other_object and is the same
Loading