| ... | @@ -28,3 +28,230 @@ The only change since then is adding parallelization to improve efficiency. |
... | @@ -28,3 +28,230 @@ The only change since then is adding parallelization to improve efficiency. |
|
|
As part of this review, I added some unit tests for the reweighting at the link above.
|
|
As part of this review, I added some unit tests for the reweighting at the link above.
|
|
|
The tests check the case where the likelihoods are the same and for a known difference.
|
|
The tests check the case where the likelihoods are the same and for a known difference.
|
|
|
The reweighting code hasn't changed since the latest release and so the passing unit tests validate at least the latest release.
|
|
The reweighting code hasn't changed since the latest release and so the passing unit tests validate at least the latest release.
|
|
|
|
|
|
|
|
## `bilby_pipe` examples
|
|
|
|
|
|
|
|
Reweight `IMRPhenomXAS` with relative binning to the regular likelihood with `IMRPhenomXHM`.
|
|
|
|
This is a modified version of [this example](https://git.ligo.org/lscsoft/bilby_pipe/-/tree/master/examples/reweighting).
|
|
|
|
The rejection sampling efficiency is ~1 / 3.
|
|
|
|
|
|
|
|
<details><summary>reweighting configuration</summary>
|
|
|
|
|
|
|
|
```json
|
|
|
|
{
|
|
|
|
"likelihood-type": "bilby.gw.likelihood.GravitationalWaveTransient",
|
|
|
|
"waveform_approximant": "IMRPhenomXHM"
|
|
|
|
}
|
|
|
|
```
|
|
|
|
|
|
|
|
</details>
|
|
|
|
|
|
|
|
<details><summary>configuration</summary>
|
|
|
|
|
|
|
|
```python
|
|
|
|
accounting = ligo.dev.o4.cbc.pe.bilby
|
|
|
|
|
|
|
|
label = bbh_injection
|
|
|
|
outdir = outdir_bbh_injection
|
|
|
|
|
|
|
|
detectors = [H1, L1]
|
|
|
|
duration = 64
|
|
|
|
|
|
|
|
sampler = dynesty
|
|
|
|
sampler-kwargs = {'nlive': 100, 'naccept': 20, 'sample': 'acceptance-walk'}
|
|
|
|
|
|
|
|
prior-file = bbh_injection.prior
|
|
|
|
|
|
|
|
injection = True
|
|
|
|
gaussian-noise = True
|
|
|
|
injection-dict = {"mass_1": 41, 'mass_2': 39, 'a_1': 0, 'a_2': 0, 'tilt_1': 0, 'tilt_2': 0,
|
|
|
|
'phi_12': 0, 'phi_jl': 0, 'luminosity_distance': 1000, 'dec': -0.2, 'ra': 1.4,
|
|
|
|
'theta_jn': 0.2, 'psi': 0, 'phase': 0, 'geocent_time': 0}
|
|
|
|
injection-waveform-approximant=IMRPhenomXHM
|
|
|
|
n-simulation = 1
|
|
|
|
n-parallel = 1
|
|
|
|
request-cpus = 4
|
|
|
|
osg = True
|
|
|
|
|
|
|
|
likelihood-type = bilby.gw.likelihood.relative.RelativeBinningGravitationalWaveTransient
|
|
|
|
fiducial_parameters = {"mass_1": 41, 'mass_2': 39, 'a_1': 0, 'a_2': 0, 'tilt_1': 0, 'tilt_2': 0,
|
|
|
|
'phi_12': 0, 'phi_jl': 0, 'luminosity_distance': 1000, 'dec': -0.2, 'ra': 1.4,
|
|
|
|
'theta_jn': 0.2, 'psi': 0, 'phase': 0, 'geocent_time': 0}
|
|
|
|
waveform-approximant = "IMRPhenomXAS"
|
|
|
|
frequency-domain-source-model = bilby.gw.source.lal_binary_black_hole
|
|
|
|
|
|
|
|
reweighting-configuration=reweight.json
|
|
|
|
reweight-nested-samples=False
|
|
|
|
```
|
|
|
|
|
|
|
|
</details>
|
|
|
|
|
|
|
|
<details><summary>Full bilby_pipe log</summary>
|
|
|
|
|
|
|
|
```python
|
|
|
|
12:00 bilby_pipe INFO : Running bilby_pipe version: 1.0.8
|
|
|
|
12:00 bilby_pipe INFO : Running bilby version: 2.0.2
|
|
|
|
12:00 bilby_pipe INFO : Command line arguments: Namespace(ini='outdir_bbh_injection/bbh_injection_config_complete.ini', verbose=False, calibration_model=None, spline_calibration_envelope_dict=None, spline_calibration_nodes=10, spline_calibration_amplitude_uncertainty_dict=None, spline_calibration_phase_uncertainty_dict=None, calibration_prior_boundary='reflective', idx=0, data_dump_file=None, ignore_gwpy_data_quality_check=True, gps_tuple=None, gps_file=None, timeslide_file=None, timeslide_dict=None, trigger_time='0', n_simulation=1, data_dict=None, data_format=None, allow_tape=True, channel_dict=None, gaussian_noise=True, zero_noise=False, coherence_test=False, detectors=['H1', 'L1'], duration=64.0, generation_seed=None, psd_dict=None, psd_fractional_overlap=0.5, post_trigger_duration=2.0, sampling_frequency=4096.0, psd_length=32, psd_maximum_duration=1024, psd_method='median', psd_start_time=None, maximum_frequency=None, minimum_frequency='20', tukey_roll_off=0.4, resampling_method='lal', injection=True, injection_dict='{"mass_1": 41, \'mass_2\': 39, \'a_1\': 0, \'a_2\': 0, \'tilt_1\': 0, \'tilt_2\': 0, \'phi_12\': 0, \'phi_jl\': 0, \'luminosity_distance\': 1000, \'dec\': -0.2, \'ra\': 1.4, \'theta_jn\': 0.2, \'psi\': 0, \'phase\': 0, \'geocent_time\': 0}', injection_file=None, injection_numbers=[None], injection_waveform_approximant='IMRPhenomXHM', injection_waveform_arguments=None, accounting='ligo.dev.o4.cbc.pe.bilby', accounting_user=None, label='bbh_injection_data0_0_generation', local=False, local_generation=False, local_plot=False, outdir='/home/colm.talbot/modules/o4-pe-review/bilby_pipe/examples/reweighting/outdir_bbh_injection', overwrite_outdir=False, periodic_restart_time=28800, request_disk=5.0, request_memory=8.0, request_memory_generation=None, request_cpus=4, conda_env=None, scheduler='condor', scheduler_args=None, scheduler_module=[None], scheduler_env=None, scheduler_analysis_time='7-00:00:00', submit=False, condor_job_priority=0, transfer_files=True, additional_transfer_paths=[None], log_directory=None, osg=True, desired_sites=None, analysis_executable=None, analysis_executable_parser=None, calibration_marginalization=False, distance_marginalization=False, distance_marginalization_lookup_table=None, phase_marginalization=False, time_marginalization=False, jitter_time=True, reference_frame='sky', time_reference='geocent', likelihood_type='bilby.gw.likelihood.relative.RelativeBinningGravitationalWaveTransient', calibration_lookup_table=None, roq_folder=None, roq_linear_matrix=None, roq_quadratic_matrix=None, roq_weights=None, roq_weight_format=None, roq_scale_factor=1.0, fiducial_parameters='{"mass_1": 41, \'mass_2\': 39, \'a_1\': 0, \'a_2\': 0, \'tilt_1\': 0, \'tilt_2\': 0, \'phi_12\': 0, \'phi_jl\': 0, \'luminosity_distance\': 1000, \'dec\': -0.2, \'ra\': 1.4, \'theta_jn\': 0.2, \'psi\': 0, \'phase\': 0, \'geocent_time\': 0}', update_fiducial_parameters=False, epsilon=0.025, extra_likelihood_kwargs=None, plot_trace=False, plot_data=False, plot_injection=False, plot_spectrogram=False, plot_calibration=False, plot_corner=False, plot_marginal=False, plot_skymap=False, plot_waveform=False, plot_format='png', create_summary=False, email=None, notification='Never', queue=None, existing_dir=None, webdir=None, summarypages_arguments=None, result_format='hdf5', final_result=True, final_result_nsamples=20000, default_prior='BBHPriorDict', deltaT=0.2, prior_file='/home/colm.talbot/modules/o4-pe-review/bilby_pipe/examples/reweighting/bbh_injection.prior', prior_dict=None, enforce_signal_duration=True, postprocessing_executable=None, postprocessing_arguments=None, single_postprocessing_executable=None, single_postprocessing_arguments=None, sampler='dynesty', sampling_seed=None, n_parallel=1, sampler_kwargs="{'nlive': 100, 'naccept': 20, 'check_point_plot': True, 'check_point_delta_t': 1800, 'print_method': 'interval-60', 'sample': 'acceptance-walk', 'npool': 4}", reweighting_configuration='/home/colm.talbot/modules/o4-pe-review/bilby_pipe/examples/reweighting/reweight.json', reweight_nested_samples=False, waveform_generator='bilby.gw.waveform_generator.LALCBCWaveformGenerator', reference_frequency=20.0, waveform_approximant='IMRPhenomXAS', catch_waveform_errors=True, pn_spin_order=-1, pn_tidal_order=-1, pn_phase_order=-1, pn_amplitude_order=0, numerical_relativity_file=None, waveform_arguments_dict=None, mode_array=[None], frequency_domain_source_model='bilby.gw.source.lal_binary_black_hole', conversion_function=None, generation_function=None)
|
|
|
|
12:00 bilby_pipe INFO : Unknown command line arguments: []
|
|
|
|
12:00 bilby_pipe INFO : Generation seed set to 56142
|
|
|
|
12:00 bilby_pipe INFO : Requested PSD duration 2048.0=32x64.0 exceeds allowed maximum 1024. Setting psd_duration = 1024
|
|
|
|
12:00 bilby_pipe INFO : Injecting waveform with
|
|
|
|
12:00 bilby_pipe INFO : minimum_frequency = 20.0
|
|
|
|
12:00 bilby_pipe INFO : maximum_frequency = 2048.0
|
|
|
|
12:00 bilby_pipe INFO : trigger_time = 0.0
|
|
|
|
12:00 bilby_pipe INFO : start_time = -62.0
|
|
|
|
12:00 bilby_pipe INFO : duration = 64.0
|
|
|
|
12:00 bilby_pipe INFO : Simulating strain data from psd-colored noise
|
|
|
|
12:00 bilby_pipe INFO : Using waveform arguments: {'reference_frequency': 20.0, 'waveform_approximant': 'IMRPhenomXHM', 'minimum_frequency': 20.0, 'maximum_frequency': 2048.0, 'catch_waveform_errors': True, 'pn_spin_order': -1, 'pn_tidal_order': -1, 'pn_phase_order': -1, 'pn_amplitude_order': 0, 'mode_array': None, 'numerical_relativity_file': None}
|
|
|
|
12:00 bilby INFO : Waveform generator initiated with
|
|
|
|
frequency_domain_source_model: bilby.gw.source.lal_binary_black_hole
|
|
|
|
time_domain_source_model: None
|
|
|
|
parameter_conversion: bilby.gw.conversion.convert_to_lal_binary_black_hole_parameters
|
|
|
|
12:00 bilby INFO : Injected signal in H1:
|
|
|
|
12:00 bilby INFO : optimal SNR = 9.82
|
|
|
|
12:00 bilby INFO : matched filter SNR = 12.78+0.38j
|
|
|
|
12:00 bilby INFO : mass_1 = 41.0
|
|
|
|
12:00 bilby INFO : mass_2 = 39.0
|
|
|
|
12:00 bilby INFO : a_1 = 0.0
|
|
|
|
12:00 bilby INFO : a_2 = 0.0
|
|
|
|
12:00 bilby INFO : tilt_1 = 0.0
|
|
|
|
12:00 bilby INFO : tilt_2 = 0.0
|
|
|
|
12:00 bilby INFO : phi_12 = 0.0
|
|
|
|
12:00 bilby INFO : phi_jl = 0.0
|
|
|
|
12:00 bilby INFO : luminosity_distance = 1000.0
|
|
|
|
12:00 bilby INFO : dec = -0.2
|
|
|
|
12:00 bilby INFO : ra = 1.4
|
|
|
|
12:00 bilby INFO : theta_jn = 0.2
|
|
|
|
12:00 bilby INFO : psi = 0.0
|
|
|
|
12:00 bilby INFO : phase = 0.0
|
|
|
|
12:00 bilby INFO : geocent_time = 0.0
|
|
|
|
12:00 bilby INFO : Injected signal in L1:
|
|
|
|
12:00 bilby INFO : optimal SNR = 11.46
|
|
|
|
12:00 bilby INFO : matched filter SNR = 11.91-0.92j
|
|
|
|
12:00 bilby INFO : mass_1 = 41.0
|
|
|
|
12:00 bilby INFO : mass_2 = 39.0
|
|
|
|
12:00 bilby INFO : a_1 = 0.0
|
|
|
|
12:00 bilby INFO : a_2 = 0.0
|
|
|
|
12:00 bilby INFO : tilt_1 = 0.0
|
|
|
|
12:00 bilby INFO : tilt_2 = 0.0
|
|
|
|
12:00 bilby INFO : phi_12 = 0.0
|
|
|
|
12:00 bilby INFO : phi_jl = 0.0
|
|
|
|
12:00 bilby INFO : luminosity_distance = 1000.0
|
|
|
|
12:00 bilby INFO : dec = -0.2
|
|
|
|
12:00 bilby INFO : ra = 1.4
|
|
|
|
12:00 bilby INFO : theta_jn = 0.2
|
|
|
|
12:00 bilby INFO : psi = 0.0
|
|
|
|
12:00 bilby INFO : phase = 0.0
|
|
|
|
12:00 bilby INFO : geocent_time = 0.0
|
|
|
|
12:00 bilby INFO : Waveform generator initiated with
|
|
|
|
frequency_domain_source_model: bilby.gw.source.lal_binary_black_hole_relative_binning
|
|
|
|
time_domain_source_model: None
|
|
|
|
parameter_conversion: bilby.gw.conversion.convert_to_lal_binary_black_hole_parameters
|
|
|
|
12:00 bilby INFO : Initial fiducial waveforms set up
|
|
|
|
12:00 bilby INFO : Summary Data Obtained
|
|
|
|
12:00 bilby INFO : Fiducial likelihood: 146.46
|
|
|
|
12:00 bilby_pipe INFO : Completed data generation
|
|
|
|
12:00 bilby_pipe INFO : Running bilby_pipe version: 1.0.8
|
|
|
|
12:00 bilby_pipe INFO : Running bilby version: 2.0.2
|
|
|
|
12:00 bilby_pipe INFO : Command line arguments: Namespace(ini='outdir_bbh_injection/bbh_injection_config_complete.ini', verbose=False, calibration_model=None, spline_calibration_envelope_dict=None, spline_calibration_nodes=10, spline_calibration_amplitude_uncertainty_dict=None, spline_calibration_phase_uncertainty_dict=None, calibration_prior_boundary='reflective', idx=0, data_dump_file='outdir_bbh_injection/data/bbh_injection_data0_0_generation_data_dump.pickle', ignore_gwpy_data_quality_check=True, gps_tuple=None, gps_file=None, timeslide_file=None, timeslide_dict=None, trigger_time=None, n_simulation=1, data_dict=None, data_format=None, allow_tape=True, channel_dict=None, gaussian_noise=True, zero_noise=False, coherence_test=False, detectors=['H1', 'L1'], duration=64.0, generation_seed=None, psd_dict=None, psd_fractional_overlap=0.5, post_trigger_duration=2.0, sampling_frequency=4096.0, psd_length=32, psd_maximum_duration=1024, psd_method='median', psd_start_time=None, maximum_frequency=None, minimum_frequency='20', tukey_roll_off=0.4, resampling_method='lal', injection=True, injection_dict='{"mass_1": 41, \'mass_2\': 39, \'a_1\': 0, \'a_2\': 0, \'tilt_1\': 0, \'tilt_2\': 0, \'phi_12\': 0, \'phi_jl\': 0, \'luminosity_distance\': 1000, \'dec\': -0.2, \'ra\': 1.4, \'theta_jn\': 0.2, \'psi\': 0, \'phase\': 0, \'geocent_time\': 0}', injection_file=None, injection_numbers=[None], injection_waveform_approximant='IMRPhenomXHM', injection_waveform_arguments=None, accounting='ligo.dev.o4.cbc.pe.bilby', accounting_user=None, label='bbh_injection_data0_0_analysis_H1L1', local=False, local_generation=False, local_plot=False, outdir='outdir_bbh_injection', overwrite_outdir=False, periodic_restart_time=28800, request_disk=5.0, request_memory=8.0, request_memory_generation=None, request_cpus=4, conda_env=None, scheduler='condor', scheduler_args=None, scheduler_module=[None], scheduler_env=None, scheduler_analysis_time='7-00:00:00', submit=False, condor_job_priority=0, transfer_files=True, additional_transfer_paths=[None], log_directory=None, osg=True, desired_sites=None, analysis_executable=None, analysis_executable_parser=None, calibration_marginalization=False, distance_marginalization=False, distance_marginalization_lookup_table=None, phase_marginalization=False, time_marginalization=False, jitter_time=True, reference_frame='sky', time_reference='geocent', likelihood_type='bilby.gw.likelihood.relative.RelativeBinningGravitationalWaveTransient', calibration_lookup_table=None, roq_folder=None, roq_linear_matrix=None, roq_quadratic_matrix=None, roq_weights=None, roq_weight_format=None, roq_scale_factor=1.0, fiducial_parameters='{"mass_1": 41, \'mass_2\': 39, \'a_1\': 0, \'a_2\': 0, \'tilt_1\': 0, \'tilt_2\': 0, \'phi_12\': 0, \'phi_jl\': 0, \'luminosity_distance\': 1000, \'dec\': -0.2, \'ra\': 1.4, \'theta_jn\': 0.2, \'psi\': 0, \'phase\': 0, \'geocent_time\': 0}', update_fiducial_parameters=False, epsilon=0.025, extra_likelihood_kwargs=None, plot_trace=False, plot_data=False, plot_injection=False, plot_spectrogram=False, plot_calibration=False, plot_corner=False, plot_marginal=False, plot_skymap=False, plot_waveform=False, plot_format='png', create_summary=False, email=None, notification='Never', queue=None, existing_dir=None, webdir=None, summarypages_arguments=None, result_format='hdf5', final_result=True, final_result_nsamples=20000, default_prior='BBHPriorDict', deltaT=0.2, prior_file='/home/colm.talbot/modules/o4-pe-review/bilby_pipe/examples/reweighting/bbh_injection.prior', prior_dict=None, enforce_signal_duration=True, postprocessing_executable=None, postprocessing_arguments=None, single_postprocessing_executable=None, single_postprocessing_arguments=None, sampler='dynesty', sampling_seed=None, n_parallel=1, sampler_kwargs="{'nlive': 100, 'naccept': 20, 'check_point_plot': True, 'check_point_delta_t': 1800, 'print_method': 'interval-60', 'sample': 'acceptance-walk', 'npool': 4}", reweighting_configuration='/home/colm.talbot/modules/o4-pe-review/bilby_pipe/examples/reweighting/reweight.json', reweight_nested_samples=False, waveform_generator='bilby.gw.waveform_generator.LALCBCWaveformGenerator', reference_frequency=20.0, waveform_approximant='IMRPhenomXAS', catch_waveform_errors=True, pn_spin_order=-1, pn_tidal_order=-1, pn_phase_order=-1, pn_amplitude_order=0, numerical_relativity_file=None, waveform_arguments_dict=None, mode_array=[None], frequency_domain_source_model='bilby.gw.source.lal_binary_black_hole', conversion_function=None, generation_function=None)
|
|
|
|
12:00 bilby_pipe INFO : Unknown command line arguments: []
|
|
|
|
12:00 bilby_pipe INFO : Sampling seed set to 677475
|
|
|
|
12:00 bilby_pipe INFO : Found data for detectors = ['H1', 'L1']
|
|
|
|
12:00 bilby_pipe INFO : Using data for detectors = ['H1', 'L1']
|
|
|
|
12:00 bilby_pipe INFO : H1: sampling-frequency=4096.0, segment-start-time=-62.0, duration=64.0
|
|
|
|
12:00 bilby_pipe INFO : L1: sampling-frequency=4096.0, segment-start-time=-62.0, duration=64.0
|
|
|
|
12:00 bilby_pipe INFO : Using conversion_function convert_to_lal_binary_black_hole_parameters
|
|
|
|
12:00 bilby INFO : Waveform generator initiated with
|
|
|
|
frequency_domain_source_model: bilby.gw.source.lal_binary_black_hole_relative_binning
|
|
|
|
time_domain_source_model: None
|
|
|
|
parameter_conversion: bilby.gw.conversion.convert_to_lal_binary_black_hole_parameters
|
|
|
|
12:00 bilby INFO : Initial fiducial waveforms set up
|
|
|
|
12:00 bilby INFO : Summary Data Obtained
|
|
|
|
12:00 bilby INFO : Fiducial likelihood: 146.46
|
|
|
|
12:00 bilby INFO : Running for label 'bbh_injection_data0_0_analysis_H1L1', output will be saved to 'outdir_bbh_injection/result'
|
|
|
|
12:00 bilby WARNING : Parameter fiducial has no default prior and is set to 0, this will not be sampled and may cause an error.
|
|
|
|
12:00 bilby INFO : Using lal version 7.2.4
|
|
|
|
12:00 bilby INFO : Using lal git version Branch: None;Tag: lal-v7.2.4;Id: e5e51f77cbaf264a20a620b4acefbe7ac38c84c6;;Builder: Adam Mercer <adam.mercer@ligo.org>;Repository status: CLEAN: All modifications committed
|
|
|
|
12:00 bilby INFO : Using lalsimulation version 4.0.2
|
|
|
|
12:00 bilby INFO : Using lalsimulation git version Branch: None;Tag: lalsimulation-v4.0.2;Id: 233cc3963d87688c272b7affe7dd0b962e4c11a0;;Builder: Adam Mercer <adam.mercer@ligo.org>;Repository status: CLEAN: All modifications committed
|
|
|
|
12:01 bilby INFO : Analysis priors:
|
|
|
|
12:01 bilby INFO : mass_1=Uniform(minimum=40, maximum=50, name=None, latex_label=None, unit=None, boundary=None)
|
|
|
|
12:01 bilby INFO : mass_2=Uniform(minimum=30, maximum=40, name=None, latex_label=None, unit=None, boundary=None)
|
|
|
|
12:01 bilby INFO : luminosity_distance=bilby.gw.prior.UniformComovingVolume(minimum=100.0, maximum=5000.0, cosmology=FlatLambdaCDM(name="Planck15", H0=67.74 km / (Mpc s), Om0=0.3075, Tcmb0=2.7255 K, Neff=3.046, m_nu=[0. 0. 0.06] eV, Ob0=0.0486), name='luminosity_distance', latex_label='$d_L$', unit=Unit("Mpc"), boundary=None)
|
|
|
|
12:01 bilby INFO : theta_jn=Sine(minimum=0, maximum=3.141592653589793, name='theta_jn', latex_label='$\\theta_{JN}$', unit=None, boundary=None)
|
|
|
|
12:01 bilby INFO : a_1=0.0
|
|
|
|
12:01 bilby INFO : a_2=0.0
|
|
|
|
12:01 bilby INFO : tilt_1=0.0
|
|
|
|
12:01 bilby INFO : tilt_2=0.0
|
|
|
|
12:01 bilby INFO : phi_12=0.0
|
|
|
|
12:01 bilby INFO : phi_jl=0.0
|
|
|
|
12:01 bilby INFO : dec=-0.2
|
|
|
|
12:01 bilby INFO : ra=1.4
|
|
|
|
12:01 bilby INFO : psi=0.0
|
|
|
|
12:01 bilby INFO : phase=0.0
|
|
|
|
12:01 bilby INFO : geocent_time=0.0
|
|
|
|
12:01 bilby INFO : Analysis likelihood class: <class 'bilby.gw.likelihood.relative.RelativeBinningGravitationalWaveTransient'>
|
|
|
|
12:01 bilby INFO : Analysis likelihood noise evidence: -259517.80870582513
|
|
|
|
12:01 bilby INFO : Single likelihood evaluation took 1.982e-03 s
|
|
|
|
12:01 bilby WARNING : Supplied argument 'sampling_seed' not an argument of 'Dynesty', removing.
|
|
|
|
12:01 bilby INFO : Using sampler Dynesty with kwargs {'nlive': 100, 'bound': 'live', 'sample': 'acceptance-walk', 'periodic': None, 'reflective': None, 'update_interval': 600, 'first_update': None, 'npdim': None, 'rstate': None, 'queue_size': 4, 'pool': None, 'use_pool': None, 'live_points': None, 'logl_args': None, 'logl_kwargs': None, 'ptform_args': None, 'ptform_kwargs': None, 'gradient': None, 'grad_args': None, 'grad_kwargs': None, 'compute_jac': False, 'enlarge': None, 'bootstrap': None, 'walks': 100, 'facc': 0.2, 'slices': None, 'fmove': 0.9, 'max_move': 100, 'update_func': None, 'ncdim': None, 'blob': False, 'save_history': False, 'history_filename': None, 'maxiter': None, 'maxcall': None, 'dlogz': 0.1, 'logl_max': inf, 'n_effective': None, 'add_live': True, 'print_progress': True, 'print_func': <bound method Dynesty._print_func of <bilby.core.sampler.dynesty.Dynesty object at 0x7f88a02a7d30>>, 'save_bounds': False, 'checkpoint_file': None, 'checkpoint_every': 60, 'resume': False}
|
|
|
|
12:01 bilby INFO : Checkpoint every check_point_delta_t = 1800s
|
|
|
|
12:01 bilby INFO : Using dynesty version 2.1.1
|
|
|
|
12:01 bilby INFO : Using the bilby-implemented rwalk sampling with an average of 20 accepted steps per MCMC and maximum length 5000
|
|
|
|
12:01 bilby INFO : Setting up multiproccesing pool with 4 processes
|
|
|
|
12:01 bilby INFO : Reading resume file outdir_bbh_injection/result/bbh_injection_data0_0_analysis_H1L1_resume.pickle
|
|
|
|
12:01 bilby WARNING : The resume file outdir_bbh_injection/result/bbh_injection_data0_0_analysis_H1L1_resume.pickle is corrupted or the version of bilby has changed between runs. This resume file will be ignored.
|
|
|
|
12:01 bilby INFO : Generating initial points from the prior
|
|
|
|
986it [0:00:59 bound:46 nc: 57 ncall:3.5e+04 eff:2.8% logz-ratio=135.55+/-0.31 dlogz:4.846>0.1]
|
|
|
|
12:02 bilby INFO : Written checkpoint file outdir_bbh_injection/result/bbh_injection_data0_0_analysis_H1L1_resume.pickle
|
|
|
|
12:02 bilby INFO : Starting to close worker pool.
|
|
|
|
12:02 bilby INFO : Finished closing worker pool.
|
|
|
|
12:02 bilby INFO : Rejection sampling nested samples to obtain 336 posterior samples
|
|
|
|
12:02 bilby INFO : Sampling time: 0:01:43.733058
|
|
|
|
12:02 bilby INFO : Summary of results:
|
|
|
|
nsamples: 336
|
|
|
|
ln_noise_evidence: -259517.809
|
|
|
|
ln_evidence: -259379.101 +/- 0.372
|
|
|
|
ln_bayes_factor: 138.707 +/- 0.372
|
|
|
|
|
|
|
|
12:02 bilby_pipe INFO : Using conversion_function convert_to_lal_binary_black_hole_parameters
|
|
|
|
12:02 bilby INFO : Waveform generator initiated with
|
|
|
|
frequency_domain_source_model: bilby.gw.source.lal_binary_black_hole_relative_binning
|
|
|
|
time_domain_source_model: None
|
|
|
|
parameter_conversion: bilby.gw.conversion.convert_to_lal_binary_black_hole_parameters
|
|
|
|
12:02 bilby INFO : Initial fiducial waveforms set up
|
|
|
|
12:02 bilby INFO : Summary Data Obtained
|
|
|
|
12:02 bilby INFO : Fiducial likelihood: 146.46
|
|
|
|
12:02 bilby_pipe INFO : Setting likelihood_type to bilby.gw.likelihood.GravitationalWaveTransient for reweighting
|
|
|
|
12:02 bilby_pipe INFO : Setting waveform_approximant to IMRPhenomXHM for reweighting
|
|
|
|
12:02 bilby_pipe INFO : Using conversion_function convert_to_lal_binary_black_hole_parameters
|
|
|
|
12:02 bilby INFO : Waveform generator initiated with
|
|
|
|
frequency_domain_source_model: bilby.gw.source.lal_binary_black_hole
|
|
|
|
time_domain_source_model: None
|
|
|
|
parameter_conversion: bilby.gw.conversion.convert_to_lal_binary_black_hole_parameters
|
|
|
|
12:03 bilby_pipe INFO : Using generation_function generate_all_bbh_parameters
|
|
|
|
12:03 bilby INFO : Reweighting posterior with 336 samples
|
|
|
|
Computing likelihoods: 100%|███████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 336/336 [00:00<00:00, 551.57it/s]
|
|
|
|
Computing likelihoods: 100%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 336/336 [00:07<00:00, 42.52it/s]
|
|
|
|
Computing priors: 0%| | 0/336 [00:00<?, ?it/s]12:03 bilby INFO : Checkpointing with 335 samples
|
|
|
|
Computing priors: 100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 336/336 [00:00<00:00, 82743.43it/s]
|
|
|
|
12:03 bilby INFO : Rejection sampling resulted in 97 samples
|
|
|
|
12:03 bilby INFO : Computing per-detector log likelihoods.
|
|
|
|
12:03 bilby INFO : Using a pool with size 4 for nsamples=97
|
|
|
|
100%|█████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 97/97 [00:02<00:00, 41.59it/s]
|
|
|
|
12:03 bilby INFO : Generating sky frame parameters.
|
|
|
|
100%|███████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 97/97 [00:00<00:00, 3983.43it/s]
|
|
|
|
12:03 bilby INFO : Computing SNRs for every sample.
|
|
|
|
12:03 bilby INFO : Using a pool with size 4 for nsamples=97
|
|
|
|
100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 97/97 [00:00<00:00, 12521.47it/s]
|
|
|
|
|
|
|
|
```
|
|
|
|
|
|
|
|
</details>
|
|
|
|
|
|
|
|
For an example with calibration reweighting, see, [this nessai run](https://ldas-jobs.ligo.caltech.edu/~michael.williams/nessai-review/gw/bilby_o4_review/GW150914/outdir_nessai_cal_update_v2/).
|
|
|
|
Here, the rejection sampling efficiency is again ~1 / 3.
|
|
|
|
|
|
|
|
This is the corner plot
|
|
|
|
|
|
|
|
 |
|
|
|
\ No newline at end of file |
