Error injecting signal in zero-noise
I wrote an ini
to inject a signal in zero noise + run PE on the data. However, I am being requested to provide a channel-dict when I run the below's generation step.
I might be misunderstanding the instructions to inject a signal in zero noise... It might be helpful to add a zero-noise
injection ini example.
ini
label="fast_injection"
outdir="outdir_fast_injection"
submit=True
accounting = ligo.dev.o3.cbc.pe.lalinference
################################################################################
## Data generation arguments
################################################################################
trigger-time=0
gaussian-noise=False
zero-noise=True
################################################################################
## Detector arguments
################################################################################
coherence-test=False
detectors=['H1']
duration=4.0
generation-seed=0
maximum-frequency=1024
minimum-frequency=20
################################################################################
## Injection arguments
################################################################################
injection=True
injection-dict={'chirp_mass': 17.0, 'mass_ratio': 0.3, 'a_1': 0.3, 'a_2': 0.2, 'tilt_1': 1.0, 'tilt_2': 2.0, 'phi_12': 5.0, 'phi_jl': 2.5, 'luminosity_distance': 500.0, 'dec': 0.2, 'ra': 3.9, 'theta_jn': 1.9, 'psi': 2.7, 'phase': 3.6, 'geocent_time': 0.0}
################################################################################
## Likelihood arguments
################################################################################
distance-marginalization=True
phase-marginalization=True
time-marginalization=True
jitter-time=True
################################################################################
## Prior arguments
################################################################################
deltaT=0.2
prior-dict = {
mass_1 = Constraint(name='mass_1', minimum=10, maximum=80),
mass_2 = Constraint(name='mass_2', minimum=10, maximum=80),
mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1, latex_label="$q$"),
chirp_mass = Uniform(name='chirp_mass', minimum=25, maximum=40, latex_label="$M_{c}$"),
luminosity_distance = PowerLaw(alpha=2, name='luminosity_distance', minimum=50, maximum=2000, unit='Mpc', latex_label='$d_L$'),
theta_jn = Sine(name='theta_jn'),
a_1 = 0.3,
a_2 = 0.2,
tilt_1 = 1.0,
tilt_2 = 2.0,
phi_12 = 5.0,
phi_jl = 2.5,
dec = 0.2,
ra = 3.9,
psi = 2.7,
phase = 3.6,
geocent_time = 0.0,
}
################################################################################
## Waveform arguments
################################################################################
waveform-approximant=IMRPhenomPv2
frequency-domain-source-model=lal_binary_black_hole
################################################################################
## Sampler arguments
################################################################################
sampler = dynesty
sampler-kwargs = {'nlive': 1500, 'walks':100, 'nact': 10}
sampling-seed = 0
logs
13:46 bilby_pipe INFO : Running bilby_pipe version: 1.0.2: (CLEAN) 68e07af 2020-10-09 09:52:42 -0500
13:46 bilby_pipe INFO : Running bilby: 1.0.3: (CLEAN) f15c154a 2020-11-12 11:17:04 +1100
13:46 bilby_pipe INFO : Command line arguments: Namespace(accounting='ligo.dev.o3.cbc.pe.lalinference', analysis_executable='/cvmfs/oasis.opensciencegrid.org/ligo/sw/conda/envs/igwn-py38/bin/bilby_pipe_analysis', calibration_model=None, catch_waveform_errors=False, channel_dict=None, coherence_test=False, condor_job_priority=0, convert_to_flat_in_component_mass=False, create_plots=False, create_summary=False, data_dict=None, data_dump_file=None, data_format=None, default_prior='BBHPriorDict', deltaT=0.2, detectors=["'H1'"], distance_marginalization=True, distance_marginalization_lookup_table=None, duration=4.0, email=None, existing_dir=None, extra_likelihood_kwargs=None, frequency_domain_source_model='lal_binary_black_hole', gaussian_noise=False, generation_seed=0, gps_file=None, gps_tuple=None, idx=0, ignore_gwpy_data_quality_check=True, ini='fast_bbh_injection.ini', injection=True, injection_dict="{'chirp_mass': 17.0, 'mass_ratio': 0.3, 'a_1': 0.3, 'a_2': 0.2, 'tilt_1': 1.0, 'tilt_2': 2.0, 'phi_12': 5.0, 'phi_jl': 2.5, 'luminosity_distance': 500.0, 'dec': 0.2, 'ra': 3.9, 'theta_jn': 1.9, 'psi': 2.7, 'phase': 3.6, 'geocent_time': 0.0}", injection_file=None, injection_numbers=None, injection_waveform_approximant=None, jitter_time=True, label='regular_rwalk_fast_bbh_injection_data0_0-0_generation', 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_regular_rwalk_fast_bbh_injection', periodic_restart_time=28800, phase_marginalization=True, 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='{mass-1: Constraint(name=\'mass_1\', minimum=10, maximum=80), mass_2 = Constraint(name=\'mass_2\', minimum=10, maximum=80), mass_ratio = Uniform(name=\'mass_ratio\', minimum=0.125, maximum=1, latex_label="$q$"), chirp_mass = Uniform(name=\'chirp_mass\', minimum=25, maximum=40, latex_label="$M_{c}$"), luminosity_distance = PowerLaw(alpha=2, name=\'luminosity_distance\', minimum=50, maximum=2000, unit=\'Mpc\', latex_label=\'$d_L$\'), theta_jn = Sine(name=\'theta_jn\'), a_1 = 0.3, a_2 = 0.2, tilt_1 = 1.0, tilt_2 = 2.0, phi_12 = 5.0, phi_jl = 2.5, dec = 0.2, ra = 3.9, psi = 2.7, phase = 3.6, geocent_time = 0.0, }', prior_file=None, 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, request_cpus=1, request_memory=4.0, request_memory_generation=None, resampling_method='lal', roq_folder=None, roq_scale_factor=1, roq_weights=None, sampler='dynesty', sampler_kwargs="{'nlive': 1500, 'walks':100, 'nact': 10}", sampling_frequency=4096, sampling_seed=0, 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=True, summarypages_arguments=None, time_marginalization=True, time_reference='geocent', timeslide_dict=None, timeslide_file=None, transfer_files=True, trigger_time='0', 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=True)
13:46 bilby_pipe INFO : Unknown command line arguments: []
13:46 bilby_pipe INFO : Generation seed set to 0
13:46 bilby_pipe INFO : Setting trigger time 0.0
13:46 bilby_pipe INFO : Setting prior-file to None
13:46 bilby_pipe INFO : Setting segment duration 4.0s
13:46 bilby_pipe INFO : PSD duration set to 128.0s, 32x the duration 4.0s
13:46 bilby_pipe INFO : No calibration_model model provided, calibration marginalization will not be used
13:46 bilby_pipe INFO : Setting PSD for H1 from data
13:46 bilby_pipe INFO : Using default PSD start time -128.0 relative to start time
13:46 bilby_pipe INFO : Getting psd-segment data for H1
Traceback (most recent call last):
...
data = DataGenerationInput(args, unknown_args)
...
psd_data = self.__get_psd_data(det)
File "../site-packages/bilby_pipe/data_generation.py", line 611, in __get_psd_data
det, self.get_channel_type(det), actual_psd_start_time, actual_psd_end_time
File "../site-packages/bilby_pipe/data_generation.py", line 380, in get_channel_type
raise BilbyPipeError("No channel-dict argument provided")
bilby_pipe.utils.BilbyPipeError: No channel-dict argument provided