From 37571da2b44b754a26e4692ca1ef9b202bc1e975 Mon Sep 17 00:00:00 2001
From: Leo Tsukada <leo.tsukada@ligo.org>
Date: Wed, 7 Aug 2019 19:07:16 -0700
Subject: [PATCH] gstlal_inspiral_add_dt_dphi_snr_ratio_pdfs : changed to pass
snr dictioary and psd file name to the output .h5 file
---
.../bin/gstlal_inspiral_add_dt_dphi_snr_ratio_pdfs | 7 +++++--
1 file changed, 5 insertions(+), 2 deletions(-)
diff --git a/gstlal-inspiral/bin/gstlal_inspiral_add_dt_dphi_snr_ratio_pdfs b/gstlal-inspiral/bin/gstlal_inspiral_add_dt_dphi_snr_ratio_pdfs
index 70f69f8ce2..539e24317d 100755
--- a/gstlal-inspiral/bin/gstlal_inspiral_add_dt_dphi_snr_ratio_pdfs
+++ b/gstlal-inspiral/bin/gstlal_inspiral_add_dt_dphi_snr_ratio_pdfs
@@ -18,11 +18,14 @@
import sys
import numpy
+import h5py
from gstlal.stats.inspiral_extrinsics import TimePhaseSNR
# Read in and combine all of the input files
-files = sys.argv[1:]
-TPS = TimePhaseSNR.from_hdf5(files[0], files[1:])
+h5_covmat = h5py.File(sys.argv[1])
+kwargs = {"SNR":h5_covmat["SNR"], "psd_fname":h5_covmat["psd"]}
+files = sys.argv[2:]
+TPS = TimePhaseSNR.from_hdf5(files[0], files[1:], **kwargs)
# compute the normalization
time, phase, deff = TimePhaseSNR.tile(NSIDE = 8, NANGLE = 17)
--
GitLab