Gitlab will migrate to a new storage backend starting 0300 UTC on 2020-04-04. We do not anticipate a maintenance window for this migration. Performance may be impacted over the weekend. Thanks for your patience.

Commit 1a4fa6a5 authored by Sylvia Biscoveanu's avatar Sylvia Biscoveanu Committed by Paul Lasky

Load data by channel name example

parent a407eb3c
#!/usr/bin/env python
"""
Tutorial to demonstrates reading and plotting the GraceDB data for GW170608,
the lightest GWTC-1 BBH event.
This tutorial must be run on a LIGO cluster.
"""
import bilby
label = 'GW170608'
outdir = label + '_out'
# List of detectors involved in the event. GW170608 is pre-August 2017 so only
# H1 and L1 were involved
detectors = ['H1', 'L1']
# Names of the channels and query types to use, and the event GraceDB ID - ref.
# https://ldas-jobs.ligo.caltech.edu/~eve.chase/monitor/online_pe/C02_clean/
# C02_HL_Pv2/lalinferencenest/IMRPhenomPv2pseudoFourPN/config.ini
channel_names = ['H1:DCH-CLEAN_STRAIN_C02', 'L1:DCH-CLEAN_STRAIN_C02']
gracedb = 'G288686'
# Duration of data around the event to use
duration = 16
# Duration of PSD data
psd_duration = 1024
# Minimum frequency and reference frequency
minimum_frequency = 10 # Hz
sampling_frequency = 2048.
# Get trigger time
candidate = bilby.gw.utils.gracedb_to_json(gracedb, outdir=outdir)
trigger_time = candidate['gpstime']
gps_start_time = trigger_time + 2. - duration
# Load the PSD data starting after the segment you want to analyze
psd_start_time = gps_start_time + duration
# Set up interferometer objects from the cache files
interferometers = bilby.gw.detector.InterferometerList([])
for channel_name in channel_names:
interferometer = bilby.gw.detector.load_data_by_channel_name(
start_time=gps_start_time, segment_duration=duration,
psd_duration=psd_duration, psd_start_time=psd_start_time,
channel_name=channel_name, sampling_frequency=sampling_frequency, outdir=outdir)
interferometer.minimum_frequency = minimum_frequency
interferometers.append(interferometer)
# Plot the strain data
interferometers.plot_data(outdir=outdir, label=label)
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment