Skip to content
Snippets Groups Projects
New look: Off

'chirp_mass' not found when `get_initial_point_from_prior` is run.

    • Closed Issue created by Avi Vajpeyi @avi.vajpeyi 
    trigger-time = 1126259462.3910  # Trigger time from G184098

label = GW150914
outdir = outdir

detectors = [H1, L1]
duration = 4
channel-dict = {H1:DCS-CALIB_STRAIN_C02, L1:DCS-CALIB_STRAIN_C02}

maximum-frequency=1024
minimum-frequency=20
sampling-frequency=2048
reference-frequency = 20
waveform-approximant = IMRPhenomPV2

deltaT = 0.2
time-marginalization=False
distance-marginalization=False
phase-marginalization=False

prior-file = 4s
distance-marginalization-lookup-table=lookup.npz

nodes=1
ntasks-per-node=1
time=24:00:00
nlive=1500
nact=10
n-parallel = 1

    11:50 bilby INFO    : Setting sampling seed = 545228
11:50 bilby INFO    : sampling_keys=['chirp_mass', 'mass_ratio', 'a_1', 'a_2', 'tilt_1', 'tilt_2', 'phi_12', 'phi_jl', 'luminosity_distance', 'dec', 'ra', 'theta_jn', 'psi', 'phase', 'geocent_time']
11:50 bilby INFO    : Periodic keys: ['phi_12', 'phi_jl', 'ra', 'psi', 'phase']
11:50 bilby INFO    : Reflective keys: []
11:50 bilby INFO    : Using priors:
11:50 bilby INFO    : chirp_mass: Uniform(minimum=12.299703, maximum=45, name='chirp_mass', latex_label='$\\mathcal{M}$', unit='$M_{\\odot}$', boundary=None)
11:50 bilby INFO    : mass_ratio: Uniform(minimum=0.125, maximum=1, name='mass_ratio', latex_label='$q$', unit=None, boundary=None)
11:50 bilby INFO    : mass_1: Constraint(minimum=1.001398, maximum=1000, name='mass_1', latex_label='$m_1$', unit=None)
11:50 bilby INFO    : mass_2: Constraint(minimum=1.001398, maximum=1000, name='mass_2', latex_label='$m_2$', unit=None)
11:50 bilby INFO    : a_1: Uniform(minimum=0, maximum=0.88, name='a_1', latex_label='$a_1$', unit=None, boundary=None)
11:50 bilby INFO    : a_2: Uniform(minimum=0, maximum=0.88, name='a_2', latex_label='$a_2$', unit=None, boundary=None)
11:50 bilby INFO    : tilt_1: Sine(minimum=0, maximum=3.141592653589793, name='tilt_1', latex_label='$\\theta_1$', unit=None, boundary=None)
11:50 bilby INFO    : tilt_2: Sine(minimum=0, maximum=3.141592653589793, name='tilt_2', latex_label='$\\theta_2$', unit=None, boundary=None)
11:50 bilby INFO    : phi_12: Uniform(minimum=0, maximum=6.283185307179586, name='phi_12', latex_label='$\\Delta\\phi$', unit=None, boundary='periodic')
11:50 bilby INFO    : phi_jl: Uniform(minimum=0, maximum=6.283185307179586, name='phi_jl', latex_label='$\\phi_{JL}$', unit=None, boundary='periodic')
11:50 bilby INFO    : luminosity_distance: UniformSourceFrame(minimum=100.0, maximum=5000.0, cosmology=FlatLambdaCDM(name="Planck15", H0=67.7 km / (Mpc s), Om0=0.307, Tcmb0=2.725 K, Neff=3.05, m_nu=[0.   0.   0.06] eV, Ob0=0.0486), name='luminosity_distance', latex_label='$d_L$', unit=Unit("Mpc"), boundary=None)
11:50 bilby INFO    : dec: Cosine(minimum=-1.5707963267948966, maximum=1.5707963267948966, name='dec', latex_label='$\\mathrm{DEC}$', unit=None, boundary=None)
11:50 bilby INFO    : ra: Uniform(minimum=0, maximum=6.283185307179586, name='ra', latex_label='$\\mathrm{RA}$', unit=None, boundary='periodic')
11:50 bilby INFO    : theta_jn: Sine(minimum=0, maximum=3.141592653589793, name='theta_jn', latex_label='$\\theta_{JN}$', unit=None, boundary=None)
11:50 bilby INFO    : psi: Uniform(minimum=0, maximum=3.141592653589793, name='psi', latex_label='$\\psi$', unit=None, boundary='periodic')
11:50 bilby INFO    : phase: Uniform(minimum=0, maximum=6.283185307179586, name='phase', latex_label='$\\phi$', unit=None, boundary='periodic')
11:50 bilby INFO    : geocent_time: Uniform(minimum=1126259462.2910001, maximum=1126259462.491, name='geocent_time', latex_label='$t_c$', unit='$s$', boundary=None)
11:50 bilby INFO    : Resume file /home/avi.vajpeyi/projects/IGWN_tests/run_a/outdir/result/GW150914_0_checkpoint_resume.pickle does not exist.
11:50 bilby INFO    : read_saved_state: 0:00:00.000319
11:50 bilby INFO    : Initializing sampling points with pool size=39
===========================worker with rank 1 crashed===========================
Traceback (most recent call last):
  File "/cvmfs/oasis.opensciencegrid.org/ligo/sw/conda/envs/igwn-py37-testing/lib/python3.7/site-packages/parallel_bilby/schwimmbad_fast.py", line 96, in __init__
    self.wait()
  File "/cvmfs/oasis.opensciencegrid.org/ligo/sw/conda/envs/igwn-py37-testing/lib/python3.7/site-packages/parallel_bilby/schwimmbad_fast.py", line 196, in wait
    result = func(arg)
  File "/cvmfs/oasis.opensciencegrid.org/ligo/sw/conda/envs/igwn-py37-testing/lib/python3.7/site-packages/parallel_bilby/utils.py", line 67, in get_initial_point_from_prior
    theta = prior_transform_function(unit)
  File "/cvmfs/oasis.opensciencegrid.org/ligo/sw/conda/envs/igwn-py37-testing/lib/python3.7/site-packages/parallel_bilby/analysis.py", line 375, in prior_transform_function
    return priors.rescale(sampling_keys, u_array)
  File "/cvmfs/oasis.opensciencegrid.org/ligo/sw/conda/envs/igwn-py37-testing/lib/python3.7/site-packages/bilby/core/prior/dict.py", line 733, in rescale
    return [result[key] for key in keys]
  File "/cvmfs/oasis.opensciencegrid.org/ligo/sw/conda/envs/igwn-py37-testing/lib/python3.7/site-packages/bilby/core/prior/dict.py", line 733, in <listcomp>
    return [result[key] for key in keys]
KeyError: 'chirp_mass'
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 1 in communicator MPI_COMM_WORLD
with errorcode 0.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------

    Package Versions

    bilby 1.1.2
bilby_pipe 1.0.4
parallel_bilby 1.0.1
    Edited

    Designs

      Child items ...

      2. Activity

      • Avi Vajpeyi added Bug label

        added Bug label

      • Avi Vajpeyi changed the description

        changed the description

      • Avi Vajpeyi
        Avi Vajpeyi @avi.vajpeyi ·
        Author Maintainer

        Testing same ini with bilby_pipe 

        15:19 bilby_pipe INFO    : Running bilby_pipe version: 1.0.4: (CLEAN) a0f8deb 2021-05-14 15:08:45 +0000
15:19 bilby_pipe INFO    : Running bilby: 1.1.2: release
15:19 bilby_pipe INFO    : Command line arguments: Namespace(accounting='ligo.dev.o3.cbc.pe.lalinference', allow_tape=True, analysis_executable=None, calibration_model=None, catch_waveform_errors=False, channel_dict='{H1:DCS-CALIB_STRAIN_C02, L1:DCS-CALIB_STRAIN_C02}', coherence_test=False, condor_job_priority=0, conversion_function=None, convert_to_flat_in_component_mass=False, create_plots=False, create_summary=False, data_dict=None, data_dump_file='outdir/data/GW150914_data0_1126259462-391_generation_data_dump.pickle', data_format=None, default_prior='BBHPriorDict', deltaT=0.2, detectors=['H1', 'L1'], distance_marginalization=False, distance_marginalization_lookup_table='lookup.npz', duration=4.0, email=None, existing_dir=None, extra_likelihood_kwargs=None, frequency_domain_source_model='lal_binary_black_hole', gaussian_noise=False, generation_function=None, generation_seed=None, gps_file=None, gps_tuple=None, idx=0, ignore_gwpy_data_quality_check=True, ini='outdir/GW150914_config_complete.ini', injection=False, injection_dict=None, injection_file=None, injection_numbers=['None'], injection_waveform_approximant=None, jitter_time=True, label='GW150914_data0_1126259462-391_analysis_H1L1_dynesty', likelihood_type='GravitationalWaveTransient', local=False, local_generation=False, local_plot=False, log_directory=None, maximum_frequency='1024', minimum_frequency='20', mode_array=[None], n_parallel=1, n_simulation=0, notification='Never', numerical_relativity_file=None, online_pe=False, osg=False, outdir='outdir', overwrite_outdir=False, periodic_restart_time=28800, phase_marginalization=False, plot_calibration=False, plot_corner=False, plot_format='png', plot_marginal=False, plot_skymap=False, plot_waveform=False, pn_amplitude_order=0, pn_phase_order=-1, pn_spin_order=-1, pn_tidal_order=-1, post_trigger_duration=2.0, postprocessing_arguments=None, postprocessing_executable=None, prior_dict=None, prior_file='4s', psd_dict=None, psd_fractional_overlap=0.5, psd_length=32, psd_maximum_duration=1024, psd_method='median', psd_start_time=None, reference_frame='sky', reference_frequency=20.0, request_cpus=1, request_memory=4.0, request_memory_generation=None, resampling_method='lal', result_format='json', roq_folder=None, roq_scale_factor=1.0, roq_weights=None, sampler='dynesty', sampler_kwargs="{'nlive': 1000, 'walks': 50, 'check_point_plot': True, 'n_check_point': 10000}", sampling_frequency=2048.0, sampling_seed=None, scheduler='condor', scheduler_analysis_time='7-00:00:00', scheduler_args=None, scheduler_env=None, scheduler_module=[None], single_postprocessing_arguments=None, single_postprocessing_executable=None, singularity_image=None, spline_calibration_amplitude_uncertainty_dict=None, spline_calibration_envelope_dict=None, spline_calibration_nodes=5, spline_calibration_phase_uncertainty_dict=None, submit=False, summarypages_arguments=None, time_marginalization=False, time_reference='geocent', timeslide_dict=None, timeslide_file=None, transfer_files=True, trigger_time='1126259462.3910', tukey_roll_off=0.4, verbose=False, waveform_approximant='IMRPhenomPV2', waveform_arguments_dict=None, waveform_generator='bilby.gw.waveform_generator.WaveformGenerator', webdir=None, zero_noise=False)
15:19 bilby_pipe INFO    : Sampling seed set to 751567
15:19 bilby_pipe INFO    : No calibration_model model provided, calibration marginalization will not be used
15:19 bilby_pipe INFO    : Found data for detectors = ['H1', 'L1']
15:19 bilby_pipe INFO    : Using data for detectors = ['H1', 'L1']
15:19 bilby_pipe INFO    : H1: sampling-frequency=2048.0, segment-start-time=1126259460.3909912, duration=4.0
15:19 bilby_pipe INFO    : L1: sampling-frequency=2048.0, segment-start-time=1126259460.3909912, duration=4.0
15:19 bilby_pipe INFO    : Using the lal_binary_black_hole source model
15:19 bilby_pipe INFO    : Using conversion_function convert_to_lal_binary_black_hole_parameters
15:19 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
15:19 bilby_pipe INFO    : Using generation_function generate_all_bbh_parameters
15:19 bilby INFO    : Running for label 'GW150914_data0_1126259462-391_analysis_H1L1_dynesty', output will be saved to 'outdir/result'
15:19 bilby INFO    : Using lal version 7.1.2
15:19 bilby INFO    : Using lal git version Branch: None;Tag: lal-v7.1.2;Id: c21902de2364e911fab0c5fa47ea5b954fbced73;;Builder: Adam Mercer <adam.mercer@ligo.org>;Repository status: CLEAN: All modifications committed
15:19 bilby INFO    : Using lalsimulation version 2.5.1
15:19 bilby INFO    : Using lalsimulation git version Branch: None;Tag: lalsimulation-v2.5.1;Id: bdbbfabd0f14ce7fcb8d574f1a07e911e62c85b5;;Builder: Adam Mercer <adam.mercer@ligo.org>;Repository status: CLEAN: All modifications committed
15:19 bilby INFO    : Search parameters:
15:19 bilby INFO    :   chirp_mass = Uniform(minimum=12.299703, maximum=45, name='chirp_mass', latex_label='$\\mathcal{M}$', unit='$M_{\\odot}$', boundary=None)
15:19 bilby INFO    :   mass_ratio = Uniform(minimum=0.125, maximum=1, name='mass_ratio', latex_label='$q$', unit=None, boundary=None)
15:19 bilby INFO    :   a_1 = Uniform(minimum=0, maximum=0.88, name='a_1', latex_label='$a_1$', unit=None, boundary=None)
15:19 bilby INFO    :   a_2 = Uniform(minimum=0, maximum=0.88, name='a_2', latex_label='$a_2$', unit=None, boundary=None)
15:19 bilby INFO    :   tilt_1 = Sine(minimum=0, maximum=3.141592653589793, name='tilt_1', latex_label='$\\theta_1$', unit=None, boundary=None)
15:19 bilby INFO    :   tilt_2 = Sine(minimum=0, maximum=3.141592653589793, name='tilt_2', latex_label='$\\theta_2$', unit=None, boundary=None)
15:19 bilby INFO    :   phi_12 = Uniform(minimum=0, maximum=6.283185307179586, name='phi_12', latex_label='$\\Delta\\phi$', unit=None, boundary='periodic')
15:19 bilby INFO    :   phi_jl = Uniform(minimum=0, maximum=6.283185307179586, name='phi_jl', latex_label='$\\phi_{JL}$', unit=None, boundary='periodic')
15:19 bilby INFO    :   luminosity_distance = UniformSourceFrame(minimum=100.0, maximum=5000.0, cosmology=FlatLambdaCDM(name="Planck15", H0=67.7 km / (Mpc s), Om0=0.307, Tcmb0=2.725 K, Neff=3.05, m_nu=[0.   0.   0.06] eV, Ob0=0.0486), name='luminosity_distance', latex_label='$d_L$', unit=Unit("Mpc"), boundary=None)
15:19 bilby INFO    :   dec = Cosine(minimum=-1.5707963267948966, maximum=1.5707963267948966, name='dec', latex_label='$\\mathrm{DEC}$', unit=None, boundary=None)
15:19 bilby INFO    :   ra = Uniform(minimum=0, maximum=6.283185307179586, name='ra', latex_label='$\\mathrm{RA}$', unit=None, boundary='periodic')
15:19 bilby INFO    :   theta_jn = Sine(minimum=0, maximum=3.141592653589793, name='theta_jn', latex_label='$\\theta_{JN}$', unit=None, boundary=None)
15:19 bilby INFO    :   psi = Uniform(minimum=0, maximum=3.141592653589793, name='psi', latex_label='$\\psi$', unit=None, boundary='periodic')
15:19 bilby INFO    :   phase = Uniform(minimum=0, maximum=6.283185307179586, name='phase', latex_label='$\\phi$', unit=None, boundary='periodic')
15:19 bilby INFO    :   geocent_time = Uniform(minimum=1126259462.2910001, maximum=1126259462.491, name='geocent_time', latex_label='$t_c$', unit='$s$', boundary=None)
15:19 bilby INFO    :   mass_1 = Constraint(minimum=1.001398, maximum=1000, name='mass_1', latex_label='$m_1$', unit=None)
15:19 bilby INFO    :   mass_2 = Constraint(minimum=1.001398, maximum=1000, name='mass_2', latex_label='$m_2$', unit=None)
15:19 bilby INFO    : Generating frequency domain strain from given time domain strain.
15:19 bilby INFO    : Applying a tukey window with alpha=0.2, roll off=0.4
15:19 bilby INFO    : Generating frequency domain strain from given time domain strain.
15:19 bilby INFO    : Applying a tukey window with alpha=0.2, roll off=0.4
15:19 bilby INFO    : Single likelihood evaluation took 9.262e-03 s
0it [00:00, ?it/s]15:19 bilby INFO    : Using sampler Dynesty with kwargs {'bound': 'multi', 'sample': 'rwalk', 'verbose': True, 'periodic': None, 'reflective': None, 'check_point_delta_t': 1800, 'nlive': 1000, 'first_update': None, 'walks': 50, 'npdim': None, 'rstate': None, 'queue_size': 1, 'pool': None, 'use_pool': None, 'live_points': None, 'logl_args': None, 'logl_kwargs': None, 'ptform_args': None, 'ptform_kwargs': None, 'enlarge': 1.5, 'bootstrap': None, 'vol_dec': 0.5, 'vol_check': 8.0, 'facc': 0.2, 'slices': 5, 'update_interval': 600, 'print_func': <bound method Dynesty._print_func of <bilby.core.sampler.dynesty.Dynesty object at 0x7fc967954090>>, 'dlogz': 0.1, 'maxiter': None, 'maxcall': None, 'logl_max': inf, 'add_live': True, 'print_progress': True, 'save_bounds': False, 'n_effective': None, 'maxmcmc': 5000, 'nact': 5}
15:19 bilby INFO    : Checkpoint every check_point_delta_t = 600s
15:19 bilby INFO    : Using dynesty version 1.0.1
15:19 bilby INFO    : Using the bilby-implemented rwalk sample method with ACT estimated walks
15:19 bilby INFO    : Resume file outdir/result/GW150914_data0_1126259462-391_analysis_H1L1_dynesty_resume.pickle does not exist.
15:19 bilby INFO    : Generating initial points from the prior
207it [00:19, 43.53it/s, bound:0 nc:  1 ncall:1.2e+03 eff:16.8% logz-ratio=-11.23+/-0.08 dlogz:52.850>0.1]^C15:19 bilby INFO    : Run interrupted by signal 2: checkpoint and exit on 77
15:19 bilby INFO    : Written checkpoint file outdir/result/GW150914_data0_1126259462-391_analysis_H1L1_dynesty_resume.pickle

        So it is able to get past the "Generating initial points from the prior" step

        Edited by Avi Vajpeyi
      • Avi Vajpeyi mentioned in merge request !109 (merged)

        mentioned in merge request !109 (merged)

      • Avi Vajpeyi mentioned in commit 463059d3

        mentioned in commit 463059d3

      • Avi Vajpeyi closed with merge request !109 (merged)

        closed with merge request !109 (merged)

      Please register or sign in to reply