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 70f69f8ce267b83bf85a1acfe63a5a3b4eecdafc..539e24317d2ce286038662c226c42fec9051ae4d 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)