Skip to contents

Project GW polarizations onto multiple detectors

proj_network.Rd

Applies antenna pattern projection to gravitational-wave polarizations (hp, hc) for a network of detectors, taking into account relative arrival times due to Earth's rotation.

Usage

proj_network(hp, hc, dets, ra, dec, pol, t_gps, t_ref = "H1")

Arguments

hp

A ts object. Plus polarization waveform.

hc

A ts object. Cross polarization waveform.

dets

A character vector of detector abbreviations (e.g., c("H1", "L1")).

ra

A numeric. Right Ascension of the source (in radians).

dec

A numeric. Declination of the source (in radians).

pol

A numeric. Polarization angle (in radians).

t_gps

A numeric. GPS time of signal arrival at Earth's center.

t_ref

A character. Reference detector for computing relative arrival times (default: "H1"). [Currently unused]

Value

A named list of projected ts objects for each detector, with the following attributes:

"tgps.net"

Vector of GPS arrival times per detector.

"dt_rel"

Matrix of pairwise relative time delays between detectors.

"tgps_remainder"

Rounding residuals for each detector’s delay.

Details

For each detector in dets, the function computes the relative arrival time of the signal based on its sky location. It then applies proj_wave to project the polarizations onto each detector using the corresponding time and antenna pattern.

Note

Requires Python environment with pycbc available via reticulate.

Examples

if (FALSE) { # \dontrun{
hp <- ts(sin(seq(0, 1, length.out = 4096)), start = 1000000000, fs = 4096)
hc <- ts(cos(seq(0, 1, length.out = 4096)), start = 1000000000, fs = 4096)
proj <- proj_network(hp, hc,
    dets = c("H1", "L1"),
    ra = 1.0, dec = -0.5,
    pol = 0.0, t_gps = 1126259462
)
} # }