Skip to content

Changes

Page history
Update Remnant fits authored by Charlie Hoy's avatar Charlie Hoy
Show whitespace changes
Inline Side-by-side
Remnant-fits.md
View page @ 17757cae
...@@ -639,23 +639,25 @@ $ python force_remnant.py ...@@ -639,23 +639,25 @@ $ python force_remnant.py
Now we check that spin evolution is performed if forced. We use the first 10 samples of the following file: Now we check that spin evolution is performed if forced. We use the first 10 samples of the following file:
```bash ```bash
$ rsync -avz --rsh="gsissh -p 2222" ldas-pcdev1.ligo.caltech.edu:/home/gregory.ashton/public_html/O3/pe_O3_S190425z_git_repo/pe_samples/PROD0_posterior_samples.dat ./precessing_bns_samples.dat $ rsync -avz --rsh="gsissh -p 2222" ldas-pcdev1.ligo.caltech.edu:/home/charlie.hoy/projects/pesummary_review/remnant_fits/parallel_evolved/PhenomPv3HM_samples.hdf5 .
$ head -n 11 precessing_bns_samples.dat > precessing_bns_samples_first10.dat
``` ```
and run the following: and run the following:
```python ```python
from pesummary.io import read from pesummary.io import read, write
f = read("/home/charlie.hoy/projects/pesummary_review/bbh_fits/precessing_bns_samples_first10.dat") f = read("/home/charlie.hoy/projects/pesummary_review/bbh_fits/precessing_bns_samples.hdf5")
f.generate_all_posterior_samples(force_remnant_computation=True, evolve_spins='ISCO', approximant="IMRPhenomPv2") write(f.parameters, f.samples[:10], file_format="dat", filename="precessing_bns_samples.dat")
f = read("/home/charlie.hoy/projects/pesummary_review/bbh_fits/precessing_bns_samples.dat")
assert "final_mass" not in f.parameters
print(f.parameters)
f.generate_all_posterior_samples(force_remnant_computation=True, f_low=20., evolve_spins='ISCO', approximant="IMRPhenomPv2")
assert "final_mass" in f.parameters assert "final_mass" in f.parameters
print(f.parameters) print(f.parameters)
``` ```
```bash ```bash
$ python force_remnant.py
2020-05-22 17:38:16 PESummary WARNING : Posterior samples for tidal deformability found in the posterior table. Applying BBH remnant fits to this system. This may not give sensible results. 2020-05-22 17:38:16 PESummary WARNING : Posterior samples for tidal deformability found in the posterior table. Applying BBH remnant fits to this system. This may not give sensible results.
2020-05-22 17:43:36 PESummary INFO : Evolving spins | 100% | 10/10 | 05:19 2020-05-22 17:43:36 PESummary INFO : Evolving spins | 100% | 10/10 | 05:19
['L1_matched_filter_abs_snr', 'L1_matched_filter_snr_angle', 'L1_optimal_snr', 'l1_spcal_amp_0', 'l1_spcal_amp_1', 'l1_spcal_amp_2', 'l1_spcal_amp_3', 'l1_spcal_amp_4', 'l1_spcal_amp_5', 'l1_spcal_amp_6', 'l1_spcal_amp_7', 'l1_spcal_amp_8', 'l1_spcal_amp_9', 'l1_spcal_phase_0', 'l1_spcal_phase_1', 'l1_spcal_phase_2', 'l1_spcal_phase_3', 'l1_spcal_phase_4', 'l1_spcal_phase_5', 'l1_spcal_phase_6', 'l1_spcal_phase_7', 'l1_spcal_phase_8', 'l1_spcal_phase_9', 'V1_matched_filter_abs_snr', 'V1_matched_filter_snr_angle', 'V1_optimal_snr', 'v1_spcal_amp_0', 'v1_spcal_amp_1', 'v1_spcal_amp_2', 'v1_spcal_amp_3', 'v1_spcal_amp_4', 'v1_spcal_amp_5', 'v1_spcal_amp_6', 'v1_spcal_amp_7', 'v1_spcal_amp_8', 'v1_spcal_amp_9', 'v1_spcal_phase_0', 'v1_spcal_phase_1', 'v1_spcal_phase_2', 'v1_spcal_phase_3', 'v1_spcal_phase_4', 'v1_spcal_phase_5', 'v1_spcal_phase_6', 'v1_spcal_phase_7', 'v1_spcal_phase_8', 'v1_spcal_phase_9', 'azimuth', 'deltalogl', 'deltalogll1', 'deltaloglv1', 'lambda_1', 'lambda_2', 'log_likelihood', 'log_prior', 'logw', 'network_matched_filter_snr', 'network_optimal_snr', 'phase', 'phi_12', 'phi_jl', 'mass_ratio', 't0', 'geocent_time', 'ra', 'dec', 'luminosity_distance', 'psi', 'chirp_mass', 'a_1', 'a_2', 'tilt_1', 'tilt_2', 'lal_amporder', 'lal_approximant', 'lal_pnorder', 'sky_frame', 'f_ref', 'flow', 'psdscaleflag', 'signalmodelflag', 'spcal_active', 'spcal_npts', 'alpha', 'theta_jn', 'l1_spcal_freq_0', 'l1_spcal_freq_1', 'l1_spcal_freq_2', 'l1_spcal_freq_3', 'l1_spcal_freq_4', 'l1_spcal_freq_5', 'l1_spcal_freq_6', 'l1_spcal_freq_7', 'l1_spcal_freq_8', 'l1_spcal_freq_9', 'v1_spcal_freq_0', 'v1_spcal_freq_1', 'v1_spcal_freq_2', 'v1_spcal_freq_3', 'v1_spcal_freq_4', 'v1_spcal_freq_5', 'v1_spcal_freq_6', 'v1_spcal_freq_7', 'v1_spcal_freq_8', 'v1_spcal_freq_9', 'mass_1', 'mass_2', 'symmetric_mass_ratio', 'total_mass', 'l1v1_relative_phase', 'H1_time', 'L1_time', 'V1_time', 'h1l1_delay', 'h1v1_delay', 'l1h1_delay', 'l1v1_delay', 'v1h1_delay', 'v1l1_delay', 'spin_1z', 'spin_2z', 'chi_eff', 'chi_tot', 'chi_p', 'redshift', 'mass_1_source', 'mass_2_source', 'total_mass_source', 'chirp_mass_source', 'lam_tilde', 'dlam_tilde', 'iota', 'theta1', 'phi_1', 'theta2', 'phi_2', 'beta', 'final_spin', 'final_mass', 'mf_source', 'e_rad', 'peak_luminosity', 'cosiota', 'cos_tilt_1', 'cos_tilt_2', 'cos_theta_jn', 'cosbeta', 'inverted_mass_ratio', 'spin_1x', 'spin_1y', 'spin_2x', 'spin_2y', 'lambda_tilde', 'delta_lambda', 'comoving_distance', 'final_mass_source', 'radiated_energy', 'L1_matched_filter_snr', 'V1_matched_filter_snr', 'cos_iota'] ['L1_matched_filter_abs_snr', 'L1_matched_filter_snr_angle', 'L1_optimal_snr', 'l1_spcal_amp_0', 'l1_spcal_amp_1', 'l1_spcal_amp_2', 'l1_spcal_amp_3', 'l1_spcal_amp_4', 'l1_spcal_amp_5', 'l1_spcal_amp_6', 'l1_spcal_amp_7', 'l1_spcal_amp_8', 'l1_spcal_amp_9', 'l1_spcal_phase_0', 'l1_spcal_phase_1', 'l1_spcal_phase_2', 'l1_spcal_phase_3', 'l1_spcal_phase_4', 'l1_spcal_phase_5', 'l1_spcal_phase_6', 'l1_spcal_phase_7', 'l1_spcal_phase_8', 'l1_spcal_phase_9', 'V1_matched_filter_abs_snr', 'V1_matched_filter_snr_angle', 'V1_optimal_snr', 'v1_spcal_amp_0', 'v1_spcal_amp_1', 'v1_spcal_amp_2', 'v1_spcal_amp_3', 'v1_spcal_amp_4', 'v1_spcal_amp_5', 'v1_spcal_amp_6', 'v1_spcal_amp_7', 'v1_spcal_amp_8', 'v1_spcal_amp_9', 'v1_spcal_phase_0', 'v1_spcal_phase_1', 'v1_spcal_phase_2', 'v1_spcal_phase_3', 'v1_spcal_phase_4', 'v1_spcal_phase_5', 'v1_spcal_phase_6', 'v1_spcal_phase_7', 'v1_spcal_phase_8', 'v1_spcal_phase_9', 'azimuth', 'deltalogl', 'deltalogll1', 'deltaloglv1', 'lambda_1', 'lambda_2', 'log_likelihood', 'log_prior', 'logw', 'network_matched_filter_snr', 'network_optimal_snr', 'phase', 'phi_12', 'phi_jl', 'mass_ratio', 't0', 'geocent_time', 'ra', 'dec', 'luminosity_distance', 'psi', 'chirp_mass', 'a_1', 'a_2', 'tilt_1', 'tilt_2', 'lal_amporder', 'lal_approximant', 'lal_pnorder', 'sky_frame', 'f_ref', 'flow', 'psdscaleflag', 'signalmodelflag', 'spcal_active', 'spcal_npts', 'alpha', 'theta_jn', 'l1_spcal_freq_0', 'l1_spcal_freq_1', 'l1_spcal_freq_2', 'l1_spcal_freq_3', 'l1_spcal_freq_4', 'l1_spcal_freq_5', 'l1_spcal_freq_6', 'l1_spcal_freq_7', 'l1_spcal_freq_8', 'l1_spcal_freq_9', 'v1_spcal_freq_0', 'v1_spcal_freq_1', 'v1_spcal_freq_2', 'v1_spcal_freq_3', 'v1_spcal_freq_4', 'v1_spcal_freq_5', 'v1_spcal_freq_6', 'v1_spcal_freq_7', 'v1_spcal_freq_8', 'v1_spcal_freq_9', 'mass_1', 'mass_2', 'symmetric_mass_ratio', 'total_mass', 'l1v1_relative_phase', 'H1_time', 'L1_time', 'V1_time', 'h1l1_delay', 'h1v1_delay', 'l1h1_delay', 'l1v1_delay', 'v1h1_delay', 'v1l1_delay', 'spin_1z', 'spin_2z', 'chi_eff', 'chi_tot', 'chi_p', 'redshift', 'mass_1_source', 'mass_2_source', 'total_mass_source', 'chirp_mass_source', 'lam_tilde', 'dlam_tilde', 'iota', 'theta1', 'phi_1', 'theta2', 'phi_2', 'beta', 'final_spin', 'final_mass', 'mf_source', 'e_rad', 'peak_luminosity', 'cosiota', 'cos_tilt_1', 'cos_tilt_2', 'cos_theta_jn', 'cosbeta', 'inverted_mass_ratio', 'spin_1x', 'spin_1y', 'spin_2x', 'spin_2y', 'lambda_tilde', 'delta_lambda', 'comoving_distance', 'final_mass_source', 'radiated_energy', 'L1_matched_filter_snr', 'V1_matched_filter_snr', 'cos_iota']
......