From 544155fecd72923c1142e7e84f6a7f3eec890667 Mon Sep 17 00:00:00 2001
From: Gregory Ashton <gregory.ashton@ligo.org>
Date: Sun, 26 May 2019 20:37:04 -0700
Subject: [PATCH] Extend the light crossing time to the full radius and add
warning message
---
bilby/gw/likelihood.py | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/bilby/gw/likelihood.py b/bilby/gw/likelihood.py
index f98c4b366..7c32e80e0 100644
--- a/bilby/gw/likelihood.py
+++ b/bilby/gw/likelihood.py
@@ -865,6 +865,7 @@ class ROQGravitationalWaveTransient(GravitationalWaveTransient):
indices, in_bounds = self._closest_time_indices(
ifo_time, self.weights['time_samples'])
if not in_bounds:
+ logger.debug("SNR calculation error: requested time at edge of ROQ time samples")
return self._CalculatedSNRs(
d_inner_h=np.nan_to_num(-np.inf), optimal_snr_squared=0,
complex_matched_filter_snr=np.nan_to_num(-np.inf),
@@ -924,7 +925,7 @@ class ROQGravitationalWaveTransient(GravitationalWaveTransient):
"""
- earth_light_crossing_time = radius_of_earth / speed_of_light
+ earth_light_crossing_time = 2 * radius_of_earth / speed_of_light
time_space = self._get_time_resolution()
delta_times = np.arange(
self.priors['geocent_time'].minimum - earth_light_crossing_time,
--
GitLab