diff --git a/gstlal-ugly/bin/gstlal_inspiral_metric_overlap b/gstlal-ugly/bin/gstlal_inspiral_metric_overlap index 4728eea54d58a1a754ab5e72a50d072719508f41..c174ef9ce44d8c5fdb50e759d1de5583e5976b7b 100755 --- a/gstlal-ugly/bin/gstlal_inspiral_metric_overlap +++ b/gstlal-ugly/bin/gstlal_inspiral_metric_overlap @@ -49,8 +49,8 @@ output = [] for n, vec1, in enumerate(vec1s): g, det = g_ij(vec1[1:]) def match(vec2, vec1 = vec1, g = g): - return (vec1[0], vec2[0], g_ij.metric_match(g, vec1[1:], vec2[1:])) - thisoutput = [row for row in map(match, vec2s) if row[2] > 0.0] + return (vec1[0], vec2[0], g_ij.pseudo_match(g, vec1[1:], vec2[1:])) + thisoutput = [row for row in map(match, vec2s) if row[2] > 0.01] print n, len(thisoutput) output += thisoutput diff --git a/gstlal-ugly/python/metric.py b/gstlal-ugly/python/metric.py index 0a7891612fa13e0e339f35eb52492696486f7bbf..7d769e6abe910109e568baafb3f67231770cc514 100644 --- a/gstlal-ugly/python/metric.py +++ b/gstlal-ugly/python/metric.py @@ -410,6 +410,11 @@ class Metric(object): else: return 0. + def pseudo_match(self, metric_tensor, c1, c2): + d2 = self.distance(metric_tensor, c1, c2)**2 + d2 = (numpy.arctan(d2**.5 * numpy.pi / 2) / numpy.pi * 2)**2 + return 1. - d2 + def explicit_match(self, c1, c2): def fftmatch(w1, w2):