Skip to content
Snippets Groups Projects
Commit 4d6480d4 authored by Gregory Ashton's avatar Gregory Ashton
Browse files

Replace get event time with GWOSC functionality

parent b8b4a781
No related branches found
No related tags found
No related merge requests found
......@@ -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