Skip to content
Snippets Groups Projects
Commit b078e588 authored by Colm Talbot's avatar Colm Talbot
Browse files

Merge branch 'fix_snr_relbin' into 'master'

BUGFIX: Fix SNR calculation with relative binning

See merge request !1198
parents 0651e9fe 1a469132
No related branches found
No related tags found
1 merge request!1198BUGFIX: Fix SNR calculation with relative binning
Pipeline #490092 failed
...@@ -1407,7 +1407,7 @@ def compute_snrs(sample, likelihood, npool=1): ...@@ -1407,7 +1407,7 @@ def compute_snrs(sample, likelihood, npool=1):
if likelihood is not None: if likelihood is not None:
if isinstance(sample, dict): if isinstance(sample, dict):
likelihood.parameters.update(sample) likelihood.parameters.update(sample)
signal_polarizations = likelihood.waveform_generator.frequency_domain_strain(sample) signal_polarizations = likelihood.waveform_generator.frequency_domain_strain(likelihood.parameters.copy())
for ifo in likelihood.interferometers: for ifo in likelihood.interferometers:
per_detector_snr = likelihood.calculate_snrs(signal_polarizations, ifo) per_detector_snr = likelihood.calculate_snrs(signal_polarizations, ifo)
sample['{}_matched_filter_snr'.format(ifo.name)] =\ sample['{}_matched_filter_snr'.format(ifo.name)] =\
...@@ -1461,7 +1461,7 @@ def _compute_snrs(args): ...@@ -1461,7 +1461,7 @@ def _compute_snrs(args):
sample = dict(sample).copy() sample = dict(sample).copy()
likelihood.parameters.update(sample) likelihood.parameters.update(sample)
signal_polarizations = likelihood.waveform_generator.frequency_domain_strain( signal_polarizations = likelihood.waveform_generator.frequency_domain_strain(
sample likelihood.parameters.copy()
) )
snrs = list() snrs = list()
for ifo in likelihood.interferometers: for ifo in likelihood.interferometers:
......
...@@ -105,7 +105,6 @@ for key in [ ...@@ -105,7 +105,6 @@ for key in [
"phase", "phase",
]: ]:
priors[key] = injection_parameters[key] priors[key] = injection_parameters[key]
priors["fiducial"] = 0
# Perform a check that the prior does not extend to a parameter space longer than the data # Perform a check that the prior does not extend to a parameter space longer than the data
priors.validate_prior(duration, minimum_frequency) priors.validate_prior(duration, minimum_frequency)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment