def inspiral_coinc_compare(a, offseta, b, offsetb, light_travel_time, e_thinca_parameter):
"""
Returns False (a & b are coincident) if they pass the ellipsoidal
thinca test. Otherwise return True.
light_travel_time if given in units of seconds.
"""
if offseta: a.set_end(a.get_end() + offseta)
if offsetb: b.set_end(b.get_end() + offsetb)
try:
# FIXME: should it be "<" or "<="?
coincident = xlaltools.XLALCalculateEThincaParameter(a, b) <= e_thinca_parameter
except ValueError:
# ethinca test failed to converge == events are not
# coincident
coincident = False
if offseta: a.set_end(a.get_end() - offseta)
if offsetb: b.set_end(b.get_end() - offsetb)
return not coincident
def inspiral_coinc_compare(a, offseta, b, offsetb, light_travel_time,
e_thinca_parameter):
"""
Returns False (a & b are coincident) if they pass the ellipsoidal
thinca test.
"""
if offseta: a.end += offseta
if offsetb: b.end += offsetb
try:
# FIXME: should it be "<" or "<="?
coincident = xlaltools.XLALCalculateEThincaParameter(
a, b) <= e_thinca_parameter
except ValueError:
# ethinca test failed to converge == events are not
# coincident
coincident = False
if offseta: a.end -= offseta
if offsetb: b.end -= offsetb
return not coincident