# Create a CoincEvent table.
coinc_table = lsctables.New(lsctables.CoincTable)
out_xmldoc.childNodes[0].appendChild(coinc_table)
# Create a CoincInspiral table.
coinc_inspiral_table = lsctables.New(lsctables.CoincInspiralTable)
out_xmldoc.childNodes[0].appendChild(coinc_inspiral_table)
# Precompute values that are common to all simulations.
detectors = [
lalsimulation.DetectorPrefixToLALDetector(ifo) for ifo in opts.detector
]
responses = [det.response for det in detectors]
locations = [det.location for det in detectors]
for sim_inspiral in progress.iterate(sim_inspiral_table):
# Unpack some values from the row in the table.
m1 = sim_inspiral.mass1
m2 = sim_inspiral.mass2
f_low = sim_inspiral.f_lower if opts.f_low is None else opts.f_low
DL = sim_inspiral.distance
ra = sim_inspiral.longitude
dec = sim_inspiral.latitude
inc = sim_inspiral.inclination
phi = sim_inspiral.coa_phase
psi = sim_inspiral.polarization
epoch = lal.LIGOTimeGPS(sim_inspiral.geocent_end_time,
sim_inspiral.geocent_end_time_ns)
gmst = lal.GreenwichMeanSiderealTime(epoch)
waveform = (sim_inspiral.waveform
# Create a CoincEvent table.
coinc_table = lsctables.New(lsctables.CoincTable)
out_xmldoc.childNodes[0].appendChild(coinc_table)
# Create a CoincInspiral table.
coinc_inspiral_table = lsctables.New(lsctables.CoincInspiralTable)
out_xmldoc.childNodes[0].appendChild(coinc_inspiral_table)
# Precompute values that are common to all simulations.
detectors = [lalsimulation.DetectorPrefixToLALDetector(ifo)
for ifo in opts.detector]
responses = [det.response for det in detectors]
locations = [det.location for det in detectors]
for sim_inspiral in progress.iterate(sim_inspiral_table):
# Unpack some values from the row in the table.
m1 = sim_inspiral.mass1
m2 = sim_inspiral.mass2
f_low = sim_inspiral.f_lower if opts.f_low is None else opts.f_low
DL = sim_inspiral.distance
ra = sim_inspiral.longitude
dec = sim_inspiral.latitude
inc = sim_inspiral.inclination
phi = sim_inspiral.coa_phase
psi = sim_inspiral.polarization
epoch = lal.LIGOTimeGPS(
sim_inspiral.geocent_end_time, sim_inspiral.geocent_end_time_ns)
gmst = lal.GreenwichMeanSiderealTime(epoch)
waveform = sim_inspiral.waveform if opts.waveform is None else opts.waveform
