From 25e3606ec2f21b25c9347211a01122e4fb57a61c Mon Sep 17 00:00:00 2001
From: Leo Tsukada <leo.tsukada@ligo.org>
Date: Fri, 12 Jul 2019 11:37:51 -0700
Subject: [PATCH] gstlal_inspiral_compute_dtdphideff_cov_matrix : use Virgo psd
for KAGRA
---
.../bin/gstlal_inspiral_compute_dtdphideff_cov_matrix | 8 ++++++--
1 file changed, 6 insertions(+), 2 deletions(-)
diff --git a/gstlal-inspiral/bin/gstlal_inspiral_compute_dtdphideff_cov_matrix b/gstlal-inspiral/bin/gstlal_inspiral_compute_dtdphideff_cov_matrix
index 7f3ff95475..8eb5bfb7e7 100755
--- a/gstlal-inspiral/bin/gstlal_inspiral_compute_dtdphideff_cov_matrix
+++ b/gstlal-inspiral/bin/gstlal_inspiral_compute_dtdphideff_cov_matrix
@@ -49,9 +49,13 @@ rho = {"L1": args.L_snr, "H1": args.H_snr, "V1": args.V_snr, "K1": args.K_snr}
psd = series.read_psd_xmldoc(ligolw_utils.load_filename(refpsd, verbose = True, contenthandler = series.PSDContentHandler))
psddict = {}
-for ifo in psd:
+for ifo in rho:
+ ifo_psd = ifo
+ # FIXME use Virgo PSD for KAGRA temporarily. This has to be fixed once KAGRA PSD is measured.
+ if rho == "K1":
+ ifo = "V1"
f = numpy.arange(len(psd[ifo].data.data)) * psd[ifo].deltaF
- psddict[ifo] = scipy.interpolate.interp1d(f, psd[ifo].data.data)
+ psddict[ifo_psd] = scipy.interpolate.interp1d(f, psd[ifo].data.data)
def moment(fmin, fmax, n, m1, m2, s1, s2, ifo, psddict = psddict, delta_f = 0.25):
farr = numpy.linspace(0, fmax, fmax / delta_f + delta_f)
--
GitLab