Skip to content
Snippets Groups Projects

Add mass2 and chi_eff to the reweighting script

Merged Sylvia Biscoveanu requested to merge postprocessing into master
2 unresolved threads
+ 2
0
@@ -24,9 +24,11 @@ data = h5py.File(filename, "r")
df = {}
df["mass_1"] = data["combined_imrphm_high_spin"]["posterior_samples"]["mass_1_source"]
df["mass_ratio"] = data["combined_imrphm_high_spin"]["posterior_samples"]["mass_ratio"]
df["mass_2"] = data["combined_imrphm_high_spin"]["posterior_samples"]["mass_2_source"]
df["a_1"] = data["combined_imrphm_high_spin"]["posterior_samples"]["a_1"]
df["a_2"] = data["combined_imrphm_high_spin"]["posterior_samples"]["a_2"]
df["spin_1z"] = data["combined_imrphm_high_spin"]["posterior_samples"]["spin_1z"]
df["chi_eff"] = data["combined_imrphm_high_spin"]["posterior_samples"]["chi_eff"]
df["redshift"] = data["combined_imrphm_high_spin"]["posterior_samples"]["redshift"]
df["prior"] = df["mass_1"]*(1 + data["combined_imrphm_high_spin"]["posterior_samples"]["redshift"])**2*\
redshift_prior.prob(data["combined_imrphm_high_spin"]["posterior_samples"]["redshift"])
Loading