Add mass2 and chi_eff to the reweighting script
No longer using my code for the PDB reweighting, just need to add mass2 and chi_eff to the reweighting script since those are being shown in the new figure.
Merge request reports
Activity
requested review from @shanika.galaudage
These are the posteriors I get when I run this script on the PDB results: https://ldas-jobs.ligo.caltech.edu/~sylvia.biscoveanu/nsbh/S230529ay/pdb_reweighted.png
- reweight_pdb.py 0 → 100644
1 import bilby 2 import h5py 3 import pandas as pd 4 import numpy as np 5 from tqdm import tqdm 6 import gwpopulation 7 8 gwpopulation.set_backend("cupy") 9 from gwpopulation.cupy_utils import xp, to_numpy 10 11 xp.random.seed(200105) 12 outdir = "/home/michael.zevin/projects/o4/gw230529/analysis_results/" 13 result = bilby.core.result.read_in_result( 14 outdir + "populations/PowerlawDipBreak/production/with_230529_mass_h_result.json" the weighting method Sylvia is using needs the run with all the events, and then it weights out that event using Eq 46, so the full event result is correct.
changed this line in version 4 of the diff
I made a provisional conda environment where I installed your fork of gwpopulation (from the pending MR that adds the
power_law_dip_break
model) just to make sure this works. But my plan was to just ask to source the conda environment you're using for the full PDB analysis to actually run the "production" version of this script.
Everything looks good @sylvia.biscoveanu!
mentioned in commit 327c5627