Skip to content
Snippets Groups Projects
Commit 25e3606e authored by Leo Tsukada's avatar Leo Tsukada
Browse files

gstlal_inspiral_compute_dtdphideff_cov_matrix : use Virgo psd for KAGRA

parent 2c973a6f
No related branches found
No related tags found
No related merge requests found
......@@ -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)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment