Skip to content
Snippets Groups Projects
Commit 0b1b355f authored by Daichi Tsuna's avatar Daichi Tsuna
Browse files

gstlal_inspiral_lvalert_sngls_plotter: back to ligo.lw

parent e4daf190
No related branches found
No related tags found
No related merge requests found
......@@ -50,10 +50,10 @@ import lal
from lal import LIGOTimeGPS
from lal import rate
from lal.utils import CacheEntry
from glue.ligolw import ligolw
from glue.ligolw import table
from glue.ligolw import lsctables
from glue.ligolw import utils as ligolw_utils
from ligo.lw import ligolw
from ligo.lw import table
from ligo.lw import lsctables
from ligo.lw import utils as ligolw_utils
from gstlal import lvalert_helper
from ligo.gracedb import rest as gracedb
......@@ -113,7 +113,7 @@ if sleep > 0:
def parse(c, gps = gps):
xmldoc = ligolw_utils.load_filename(c.path, contenthandler = LIGOLWContentHandler, verbose = True)
sngl_inspiral_table = lsctables.SnglInspiralTable.get_table(xmldoc)
return [row for row in sngl_inspiral_table if row.get_end() < gps + 600 and row.get_end() >= gps - 600]
return [row for row in sngl_inspiral_table if row.end < gps + 600 and row.end >= gps - 600]
def mchirp(m1,m2):
......@@ -160,7 +160,7 @@ maxmchirp, minmchirp = max(mchirps), min(mchirps)
# Get various data by IFO
times, snrs, mchirps, chisqs, ba = {}, {}, {}, {}, {}
for n,ifo in enumerate(["H1", "L1"]):
times[ifo] = numpy.array([float(r.get_end() - gps) for r in sngl if r.ifo == ifo])
times[ifo] = numpy.array([float(r.end - gps) for r in sngl if r.ifo == ifo])
snrs[ifo] = numpy.array([r.snr for r in sngl if r.ifo == ifo])
chisqs[ifo] = numpy.array([r.chisq for r in sngl if r.ifo == ifo])
mchirps[ifo] = numpy.array([mchirp(r.mass1, r.mass2) for r in sngl if r.ifo == ifo])
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment