gwastro/pycbc-tutorials

Learn how to use PyCBC to analyze gravitational-wave data and do parameter inference.

PyCBC: Python Software to Study Gravitational Waves

PyCBC is software developed by a collaboration of LIGO, Virgo, and independent scientists. It is open source and freely available. We use PyCBC in the detection of gravitational waves from binary mergers such as GW150914. These examples explore how to analyze gravitational wave data, how we find potential signals, and learn about them.

These notebooks are available to view, download, or run in interactive sessions.

Run tutorials from your browser!

Gravitational-wave Data Analysis

Tutorial 1: Accessing Gravitational-wave data

Tutorial 2: Data visualization and basic signal processing

Tutorial 3: Matched filtering to identify signals

Tutorial 4: Signal Consistency and Basic Significance Testing

Gravitational-wave Inference

Many of these tutorials will require you to make edits to config files as part of their exercises. At the moment this isn't easy to do on services like google colab. However, you can do them either on your local machine or by using services such as mybinder or sciserver which allow a full juypterhub experience with the ability to view and edit text files. Below we give links to the tutorials which should directly work in google colab, however.

Tutorial 0: Overview (this is a summary of inference tutorials 1-6; meant for a shorter tutorial session)

See video presentation of Tutorial 0 by Collin Capano

Tutorial 1: Models (likelihood function you are trying to sample + priors) and Samplers by Hand

Tutorial 5: Results files and Plotting

Tutorial 7: Adding Custom Waveforms to PyCBC

Tutorial 8: Black hole ringdown analysis with PyCBC

Tutorial 9: Adding custom models to PyCBC (e.g. alternate likelihoods, GRB / Kilonova models, etc.)

Some things that you may learn

• How to access LIGO data
• How to do some basic signal processing
• Data visualization of LIGO data in time-frequency plots
• Matched filtering to extract a known signal

Other ways to run in the browser

pycbc-tutorials's Issues

gw151226_snr.ipynb: A few minor corrections

Some minor corrections of typos that are not essential but may confuse the beginners:

In PyCBC-Tutorials/examples/gw151226_snr.ipynb:

In [3]: # Generate a waveform similar to LVT151012: ------------------------> GW151226 is correct

In [8]: # Get the data for H1 from the LOSC site: ------------------------------> L1 is correct
# Read the Hanford data and remove low-frequency content: -----> Livingston is correct
# Let's look closer around the H1 peak: ---------------------------------> L1 is correct

In PyCBC-Tutorials/tutorial/4_ChisqSignificance.ipynb:
The block In[2] is repeated.

Thank you for the wonderful tutorials.

Kamran Derakhshani

Files including read_frame() do not work

All PyCBC-Tutorial files that include read_frame() do NOT work, with such an error:

File "/home/master/miniconda3/envs/py27/lib/python2.7/site-packages/pycbc/frame/frame.py", line 135, in locations_to_cache
cache = lalframe.FrOpen(str(dir_name), str(file_name)).cache

RuntimeError: RETURN: untrapped XLAL error

My system:
Ubuntu 14.04 LTS
Python 2.7 (under Spyder in miniconda3)
PyCBC 1.14.4
LALFrame 1.4.5

Thank you.

Inference1: Issue with the SingleTemplate model

It seems that the models.SingleTemplate(...) functions have a problem. It returns the following error:
TypeError: init() got multiple values for keyword argument 'psds'

Issue in Tutorial 0 section.2 read_frame()

read_frame() gives out an error:

module 'lalframe' has no attribute 'FrSetMode'

I have run the code on the same colab notebook.

Frame file binder links not working

I'm trying to run your notebook on colab.research.google.com. I run into trouble where the binder frame file links are not working. E.g., in gw150914_audio.ipynb, the notebook tries downloading

https://raw.githubusercontent.com/ligo-cbc/binder/master/H-H1_LOSC_4_V2-1126259446-32.gwf

which returns a 15-byte file that can't be opened. But if I replace the URL with

https://losc.ligo.org/s/events/GW150914/H-H1_LOSC_4_V2-1126259446-32.gwf

then the notebook works as expected.

Perhaps the frame URLs need to be updated?

Make tutorial 3 more stand-a-lone

The matched filtering tutorials is not quite stand-a-lone for class / workshop use because it doesn't explain how to read frame files to do the challenge sets. It would be good to add an optional refresher at the beginning to demonstrate how to do this.

HDF5 file

Hello. I have data for neutron star mergers in hdf5 file format saved in my computer. How can I upload it so that I can visualize the data and follow tutorial 2 for my own data?

Posterior Layout Adjustment

I would like to know how to modify the posterior layout like color/line style/ Axes line size...etc. I have found documentation of pycbc_inference_plot_posterior dated back to 2016, which does not seem to apply to the pycbc version used here. Thanks!