Skip to content
Snippets Groups Projects
Commit c9deeb82 authored by Colm Talbot's avatar Colm Talbot
Browse files

Merge branch 'update-get-event-time' into 'master'

Replace get event time with GWOSC functionality

See merge request !1051
parents 260e68c3 4d6480d4
No related branches found
No related tags found
1 merge request!1051Replace get event time with GWOSC functionality
Pipeline #343188 passed
......@@ -387,54 +387,31 @@ def zenith_azimuth_to_ra_dec(zenith, azimuth, geocent_time, ifos):
def get_event_time(event):
"""
Get the merger time for known GW events.
See https://www.gw-openscience.org/catalog/GWTC-1-confident/html/
Last update https://arxiv.org/abs/1811.12907:
- GW150914
- GW151012
- GW151226
- GW170104
- GW170608
- GW170729
- GW170809
- GW170814
- GW170817
- GW170818
- GW170823
Get the merger time for known GW events using the gwosc package
Parameters
==========
----------
event: str
Event descriptor, this can deal with some prefixes, e.g.,
'151012', 'GW151012', 'LVT151012'
Event name, e.g. GW150914
Returns
=======
-------
event_time: float
Merger time
"""
event_times = {'150914': 1126259462.4,
'151012': 1128678900.4,
'151226': 1135136350.6,
'170104': 1167559936.6,
'170608': 1180922494.5,
'170729': 1185389807.3,
'170809': 1186302519.8,
'170814': 1186741861.5,
'170817': 1187008882.4,
'170818': 1187058327.1,
'170823': 1187529256.5}
if 'GW' or 'LVT' in event:
event = event[-6:]
Raises
------
ImportError
If the gwosc package is not installed
ValueError
If the event is not in the gwosc dataset
"""
try:
event_time = event_times[event[-6:]]
return event_time
except KeyError:
print('Unknown event {}.'.format(event))
return None
from gwosc import datasets
except ImportError:
raise ImportError("You do not have the gwosc package installed")
return datasets.event_gps(event)
def get_open_strain_data(
......
......@@ -103,21 +103,13 @@ class TestGWUtils(unittest.TestCase):
def test_get_event_time(self):
events = [
"GW150914",
"GW151012",
"GW151226",
"GW170104",
"GW170608",
"GW170729",
"GW170809",
"GW170814",
"GW170817",
"GW170818",
"GW170823",
]
for event in events:
self.assertTrue(isinstance(gwutils.get_event_time(event), float))
self.assertTrue(gwutils.get_event_time("GW010290") is None)
with self.assertRaises(ValueError):
gwutils.get_event_time("GW010290")
def test_read_frame_file(self):
start_time = 0
......
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