""" __author__ = 'Duncan Macleod <*****@*****.**>' # .. currentmodule:: gwpy.segments # # Getting the segments # -------------------- # # First, we need to fetch the Open Data timeline segments from LOSC. # To do that we can call the :meth:`DataQualityFlag.fetch_open_data` method # using ``'H1_DATA'`` as the flag (for an explanation of what this means, # read up on `The S6 Data Release <https://losc.ligo.org/S6/>`__). from gwpy.segments import DataQualityFlag h1segs = DataQualityFlag.fetch_open_data('H1_DATA', 'Sep 16 2010', 'Sep 17 2010') # For sanity, lets plot these segments: splot = h1segs.plot(figsize=[12, 3]) splot.show() splot.close() # hide # We see that the LIGO Hanford Observatory detector was operating for the # majority of the day, with a few outages of ~30 minutes or so. # We can use the :func:`abs` function to display the total amount of time # spent taking data: print(abs(h1segs.active))
from gwpy.segments import DataQualityFlag h1segs = DataQualityFlag.fetch_open_data('H1_DATA', 'Sep 12 2015', 'Jan 19 2016')
""" __author__ = 'Duncan Macleod <*****@*****.**>' # .. currentmodule:: gwpy.segments # # Getting the segments # -------------------- # # First, we need to fetch the Open Data timeline segments from LOSC. # To do that we can call the :meth:`DataQualityFlag.fetch_open_data` method # using ``'H1_DATA'`` as the flag (for an explanation of what this means, # read up on `The S6 Data Release <https://losc.ligo.org/S6/>`__). from gwpy.segments import DataQualityFlag h1segs = DataQualityFlag.fetch_open_data('H1_DATA', 'Sep 16 2010', 'Sep 17 2010') # For sanity, lets plot these segments: splot = h1segs.plot(figsize=[12, 3]) splot.show() splot.close() # hide # We see that the LIGO Hanford Observatory detector was operating for the # majority of the day, with a few outages of ~30 minutes or so. # We can use the :func:`abs` function to display the total amount of time # spent taking data: print(abs(h1segs.active))
`Observing Run 1 (O1) <https://www.gw-openscience.org/O1/>`__ have been released by |GWOSC|_. This example demonstrates how to download segment information into a :class:`~gwpy.segments.DataQualityFlag`, and then plot them. """ __author__ = 'Duncan Macleod <*****@*****.**>' __currentmodule__ = 'gwpy.segments' # All we need to do is import the `DataQualityFlag` object, and then call # the :meth:`DataQualityFlag.fetch_open_data` method to query for, and download # the segments for all of O1: from gwpy.segments import DataQualityFlag h1segs = DataQualityFlag.fetch_open_data('H1_DATA', 'Sep 12 2015', 'Jan 19 2016') # We can then generate a plot of the times when LIGO-Hanford was operating: plot = h1segs.plot(color='gwpy:ligo-hanford') plot.show() # That's a lot of segments. We can pare-down the list a little to display # only the segments from the first month of the run: h1month1 = DataQualityFlag.fetch_open_data('H1_DATA', 'Sep 12 2015', 'Oct 12 2015') # We can also download the LIGO-Livingston segments from the same period # and display them alongside, as well as those segments during which both # interferometers were operating at the same time
The data from the full Observing Run 1 (O1) have been released through the `LIGO Open Science Center <https://losc.ligo.org/detector_status/>`__ (LOSC). This example demonstrates how to download segment information into a :class:`~gwpy.segments.DataQualityFlag`, and then plot them. """ __author__ = 'Duncan Macleod <*****@*****.**>' __currentmodule__ = 'gwpy.segments' # All we need to do is import the `DataQualityFlag` object, and then call # the :meth:`DataQualityFlag.fetch_open_data` method to query for, and download # the segments for all of O1: from gwpy.segments import DataQualityFlag h1segs = DataQualityFlag.fetch_open_data('H1_DATA', 'Sep 12 2015', 'Jan 19 2016') # We can then generate a plot of the times when LIGO-Hanford was operating: plot = h1segs.plot(color='gwpy:ligo-hanford') plot.show() # That's a lot of segments. We can pare-down the list a little to display # only the segments from the first month of the run: h1month1 = DataQualityFlag.fetch_open_data('H1_DATA', 'Sep 12 2015', 'Oct 12 2015') # We can also download the LIGO-Livingston segments from the same period # and display them alongside, as well as those segments during which both # interferometers were operating at the same time