From 98ba59bb4d6a62c97d6fc416ded1b49ab9e7e5c5 Mon Sep 17 00:00:00 2001
From: Chad Hanna <chad.hanna@comp-hd-002.gwave.ics.psu.edu>
Date: Wed, 5 Jun 2019 12:54:30 -0400
Subject: [PATCH] inspiral_extrinsics: update dd covariance matrix term
---
gstlal-inspiral/python/stats/inspiral_extrinsics.py | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/gstlal-inspiral/python/stats/inspiral_extrinsics.py b/gstlal-inspiral/python/stats/inspiral_extrinsics.py
index daf6e0513c..9b875d7290 100644
--- a/gstlal-inspiral/python/stats/inspiral_extrinsics.py
+++ b/gstlal-inspiral/python/stats/inspiral_extrinsics.py
@@ -1152,7 +1152,7 @@ class TimePhaseSNR(object):
transtp = {frozenset(['V1', 'H1']): -1.6205914017601712, frozenset(['H1', 'L1']): -3.0734607125393643, frozenset(['V1', 'L1']): -1.7020167980031775}
transpt = {frozenset(['V1', 'H1']): -0, frozenset(['H1', 'L1']): -0, frozenset(['V1', 'L1']): -0}
transpp = {frozenset(['V1', 'H1']): 1.2422379329336195, frozenset(['H1', 'L1']): 2.6540061786001892, frozenset(['V1', 'L1']): 1.2984050700516516}
- transdd = {frozenset(['V1', 'H1']): 5.0, frozenset(['H1', 'L1']): 5.0, frozenset(['V1', 'L1']): 5.0}
+ transdd = {frozenset(['V1', 'H1']): 4.20997920997921, frozenset(['H1', 'L1']): 16.554054054054053, frozenset(['V1', 'L1']): 4.5884393063583815}
def __init__(self, tree_data = None, margsky = None, verbose = False, margstart = 0, margstop = None):
"""
@@ -1392,7 +1392,7 @@ class TimePhaseSNR(object):
coordtransmat = numpy.array([[self.transtt[frozenset((ifo1, ifo2))], self.transtp[frozenset((ifo1, ifo2))]],[self.transpt[frozenset((ifo1, ifo2))], self.transpp[frozenset((ifo1, ifo2))]]])
out[:,slc[0]], out[:,slc[1]] = numpy.dot(coordtransmat, dtdphivec)
# FIXME should this be the ratio - 1 or without the -1 ???
- out[:,slc[2]] = (deff[ifo1] / deff[ifo2] - 1) * self.transdd[frozenset((ifo1, ifo2))]
+ out[:,slc[2]] = numpy.log(deff[ifo1] / deff[ifo2]) * self.transdd[frozenset((ifo1, ifo2))]
return out
--
GitLab