From f1972eeaed697c2f98d9d1736dd61c4c6e2d338d Mon Sep 17 00:00:00 2001
From: Kipp Cannon <kipp.cannon@ligo.org>
Date: Tue, 19 Jun 2018 17:02:09 +0900
Subject: [PATCH] lloidparts: enable FIR phase adjustment
- pass the SVD bank's PSD to mkwhitened_multirate_src() for use as a phase reference
---
gstlal-inspiral/python/lloidparts.py | 25 +++++++++++++++++++++----
1 file changed, 21 insertions(+), 4 deletions(-)
diff --git a/gstlal-inspiral/python/lloidparts.py b/gstlal-inspiral/python/lloidparts.py
index c63656b156..bde4baefcf 100644
--- a/gstlal-inspiral/python/lloidparts.py
+++ b/gstlal-inspiral/python/lloidparts.py
@@ -1027,15 +1027,32 @@ def mkLLOIDmulti(pipeline, detectors, banks, psd, psd_fft_length = 32, ht_gate_t
#
# construct dictionaries of whitened, conditioned, down-sampled
- # h(t) streams
+ # h(t) streams. NOTE: we assume all banks for each instrument
+ # were generated with the same processed PSD for that instrument
+ # and just extract the first without checking that this assumption
+ # is correct
#
+ assert psd_fft_length % 4 == 0, "psd_fft_length (= %g) must be multiple of 4" % psd_fft_length
hoftdicts = {}
for instrument in detectors.channel_dict:
- rates = set(rate for bank in banks[instrument] for rate in bank.get_rates())
src, statevector, dqvector = datasource.mkbasicsrc(pipeline, detectors, instrument, verbose)
- assert psd_fft_length % 4 == 0, "psd_fft_length (= %g) must be multiple of 4" % psd_fft_length
- hoftdicts[instrument] = multirate_datasource.mkwhitened_multirate_src(pipeline, src, rates, instrument, psd = psd[instrument], psd_fft_length = psd_fft_length, ht_gate_threshold = ht_gate_threshold, veto_segments = veto_segments[instrument] if veto_segments is not None else None, nxydump_segment = nxydump_segment, track_psd = track_psd, width = 32, statevector = statevector, dqvector = dqvector)
+ hoftdicts[instrument] = multirate_datasource.mkwhitened_multirate_src(
+ pipeline,
+ src = src,
+ rates = set(rate for bank in banks[instrument] for rate in bank.get_rates()),
+ instrument = instrument,
+ psd = psd[instrument],
+ psd_fft_length = psd_fft_length,
+ ht_gate_threshold = ht_gate_threshold,
+ veto_segments = veto_segments[instrument] if veto_segments is not None else None,
+ nxydump_segment = nxydump_segment,
+ track_psd = track_psd,
+ width = 32,
+ statevector = statevector,
+ dqvector = dqvector,
+ fir_whiten_reference_psd = banks[instrument][0].psd[instrument]
+ )
#
# build gate control branches
--
GitLab