Transfer function

In this example we demonstrate how to calculate the transfer function between two TimeSeries signals.

All ground-based gravitational wave observatories would be unable to operate if they did not employ sophisticated ground-motion suppression technology to prevent vibrations from the local (or remote) environment from transferring through to optical components.

The impact of the seismic isolation system can be seen by calculating the transfer function between the ground motion at the laboratory and that of the optical suspension points.

from gwpy.time import tconvert
from gwpy.timeseries import TimeSeriesDict
from gwpy.plot import BodePlot

For this example we will use data from the Auxiliary Channel Three Hour Release:

start = tconvert("August 14 2017 10:25")
end = start + 1800
gndchannel = "L1:ISI-GND_STS_ITMY_Y_DQ"
suschannel = "L1:ISI-ITMY_SUSPOINT_ITMY_EUL_L_DQ"

We can call the get() method of the TimeSeriesDict to retrieve all data in a single operation:

data = TimeSeriesDict.get(
    [gndchannel, suschannel],
    start,
    end,
    verbose=True,
    host="nds.gwosc.org",
)
gnd = data[gndchannel]
sus = data[suschannel]

- Attempting data access from files
/home/duncan.macleod/gwpy-nightly-build/conda/envs/gwpy-nightly-3.11/lib/python3.11/site-packages/igwn_auth_utils/requests.py:46: DeprecationWarning: Support for identity-based X.509 credentials for LIGO.ORG is being dropped.
Calls to this utility will stop working on/around 20 May 2025.

For details on this change please see

https://computing.docs.ligo.org/guide/compsoft/roadmap/LVK/x509_retirement/

If you have questions about this message, or its implications, please
consider opening an IGWN Computing Help Desk ticket:

https://git.ligo.org/computing/helpdesk/-/issues/new
  return func(*args, **kwargs)
Cannot locate the following channel(s) in any known frametype at GPS=1186741518:
    L1:ISI-GND_STS_ITMY_Y_DQ
    L1:ISI-ITMY_SUSPOINT_ITMY_EUL_L_DQ
[files on tape have not been checked, use allow_tape=True for a complete search]
Data access from files failed
- Attempting data access from NDS2

The transfer function between time series is easily computed with the transfer_function() method:

tf = gnd.transfer_function(sus, fftlength=128, overlap=64)

The BodePlot knows how to separate a complex-valued FrequencySeries into magnitude and phase:

plot = BodePlot(tf)
plot.maxes.set_title(
    r"L1 ITMY ground $\rightarrow$ SUS transfer function",
)
plot.maxes.set_xlim(5e-2, 30)
plot.show()

# This example demonstrates the impressive noise suppression of the LIGO
# seismic isolation system. For more details, please see
# https://www.ligo.caltech.edu/page/vibration-isolation.