swh:1:snp:3a699297f000109a1bc833f294a54171df990207
Tip revision: d79adac1f3fd34f1c6bb602c3d87cd79e39023eb authored by Alex Nitz on 25 August 2021, 11:40:32 UTC
Update setup.py (#3772)
Update setup.py (#3772)
Tip revision: d79adac
frame.rst
###################################################
Reading Gravitational-wave Frames
###################################################
============
Introduction
============
All data generated and recorded by the current generation of ground-based laser-interferometer gravitational-wave detectors are recorded in gravitational-wave frame (GWF) files. These files typically contain data from a number of sources bundled into a single, time-stamped set, along with the metadata for each channel.
=====================
Querying a LDR server
=====================
The LIGO Data Replicator (LDR) is a tool for replicating data sets to the different data grids. If you have access to a LDR server you can read GWF files using ``pycbc.frame`` module as follows::
>>> from pycbc import frame
>>> tseries = frame.query_and_read_frame("G1_RDS_C01_L3", "G1:DER_DATA_H", 1049587200, 1049587200 + 60)
This returns a ``TimeSeries`` instance of the data. Note if you do not have access to frames through an LDR server then you will need to copy the frames to your run location.
Alternatively, if you just want the location of the frame files, you can do::
>>> from pycbc import frame
>>> frame_files = frame.frame_paths("G1_RDS_C01_L3", 1049587200, 1049587200 + 60)
This will return a ``list`` of the frame files' paths.
=====================
Reading a frame file
=====================
The ``pycbc.frame`` module provides methods for reading these files into ``TimeSeries`` objects as follows::
>>> from pycbc import frame
>>> data = frame.read_frame('G-G1_RDS_C01_L3-1049587200-60.gwf', 'G1:DER_DATA_H')
Here the first argument is the path to the frame file of interest, while the second lists the `data channel` of interest whose data exist within the file.
=====================
Writing a frame file
=====================
The ``pycbc.frame`` modules provides a method for writing ``TimeSeries`` instances to GWF as follows::
>>> from pycbc import frame
>>> from pycbc import types
>>> data = types.TimeSeries(numpy.ones(16384 * 16), delta_t=1.0 / 16384)
>>> frame.write_frame("./test.gwf", "H1:TEST_UNITY", data)
Here the first argument is the path to the output frame file, the second is the name of the channel, and the last is the ``TimeSeries`` instance to be written to the frame.
====================
Method documentation
====================
.. automodule:: pycbc.frame
:noindex:
:members: