Commit 8545856a authored by Leo Tsukada's avatar Leo Tsukada

gstlal_inspiral_compute_dtdphideff_cov_matrix : fix the bug to take a list of instrument pairs

parent b803d993
......@@ -52,7 +52,7 @@ psddict = {}
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":
if ifo == "K1":
ifo = "V1"
f = numpy.arange(len(psd[ifo].data.data)) * psd[ifo].deltaF
psddict[ifo_psd] = scipy.interpolate.interp1d(f, psd[ifo].data.data)
......@@ -83,7 +83,8 @@ transdd = {}
# transdd = {frozenset(("H1", "L1")): 1. / numpy.sqrt(sigsqdd["H1"] + sigsqdd["L1"]) , frozenset(("H1", "V1")): 1. / numpy.sqrt(sigsqdd["H1"] + sigsqdd["V1"]), frozenset(("L1", "V1")): 1. / numpy.sqrt(sigsqdd["L1"] + sigsqdd["V1"])}
# transdd = {frozenset(instrument_pair) : 1. / numpy.sqrt(sigsqdd[instrument_pair[0]] + sigsqdd[instrument_pair[1]]) for instrument_pair in TimePhaseSNR.instument_pairs(("H1", "L1", "V1", "K1"))}
for pair in TimePhaseSNR.instument_pairs(("H1", "L1", "V1", "K1")):
pairs = [tuple(sorted(pair)) for pair in itertools.combinations(("H1", "L1", "V1", "K1"), 2)]
for pair in pairs:
a,b = pair
m11 = sigsqtt[a] + sigsqtt[b]
m22 = sigsqpp[a] + sigsqpp[b]
......@@ -107,7 +108,7 @@ for group, mat in zip((h5_transtt, h5_transtp, h5_transpt, h5_transpp, h5_transd
for k,v in mat.items():
group.create_dataset(",".join(sorted(k)), data = float(v))
combos = TimePhaseSNR.instrument_combos(("H1","L1","V1", "K1")) + (("H1",),("L1",),("V1",), ("K1",))
combos = TimePhaseSNR.instrument_combos(("H1","L1","V1", "K1"), min_instruments=1)
norm = dict((frozenset(k), 0.) for k in combos)
norm = {}
h5_norm = f.create_group("norm")
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment