Remnant-fits.md

@@ -636,16 +636,16 @@ $ python force_remnant.py
 ['H1_matched_filter_abs_snr', 'H1_matched_filter_snr_angle', 'H1_optimal_snr', 'H1_spcal_amp_0', 'H1_spcal_amp_1', 'H1_spcal_amp_2', 'H1_spcal_amp_3', 'H1_spcal_amp_4', 'H1_spcal_amp_5', 'H1_spcal_amp_6', 'H1_spcal_amp_7', 'H1_spcal_amp_8', 'H1_spcal_amp_9', 'H1_spcal_phase_0', 'H1_spcal_phase_1', 'H1_spcal_phase_2', 'H1_spcal_phase_3', 'H1_spcal_phase_4', 'H1_spcal_phase_5', 'H1_spcal_phase_6', 'H1_spcal_phase_7', 'H1_spcal_phase_8', 'H1_spcal_phase_9', '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', 'lambda_2', 'log_likelihood', 'loglH1', 'loglL1', 'loglV1', 'logpost', 'log_prior', 'network_matched_filter_snr', 'nullLogL', 'network_optimal_snr', 'phase', 'mass_ratio', 't0', 'temperature', 'geocent_time', 'nLocalTemps', 'randomSeed', 'ra', 'dec', 'luminosity_distance', 'psi', 'chirp_mass', 'a_1', 'a_2', 'alpha', 'theta_jn', 'chain_log_evidence', 'chain_delta_log_evidence', 'chain_log_noise_evidence', 'chain_log_bayes_factor', 'tilt_1', 'tilt_2', 'inverted_mass_ratio', 'mass_1', 'mass_2', 'total_mass', 'symmetric_mass_ratio', 'iota', 'spin_1x', 'spin_1y', 'spin_1z', 'spin_2x', 'spin_2y', 'spin_2z', 'phi_1', 'phi_2', 'chi_eff', 'chi_p', 'final_spin_non_evolved', 'peak_luminosity_non_evolved', 'final_mass_non_evolved', 'cos_tilt_1', 'cos_tilt_2', 'redshift', 'comoving_distance', 'mass_1_source', 'mass_2_source', 'total_mass_source', 'chirp_mass_source', 'final_mass_source_non_evolved', 'radiated_energy_non_evolved', 'H1_time', 'L1_time', 'V1_time', 'H1_matched_filter_snr', 'L1_matched_filter_snr', 'V1_matched_filter_snr', 'cos_theta_jn', 'cos_iota']
 ```
 
-Now we check that spin evolution is performed if forced. We use the following file:
+Now we check that spin evolution is performed if forced. We use the first 10 samples of the following file:
 
 ```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.hdf5 ./precessing_bns_samples.hdf5
+$ 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
 ```
 
 ```python
 from pesummary.io import read
 
-f = read("/home/charlie.hoy/projects/pesummary_review/bbh_fits/precessing_bns_samples.hdf5")
+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, evolve_spins=True)