Skip to content

Commit

Permalink
Empty tables fixes
Browse files Browse the repository at this point in the history
  • Loading branch information
pannarale committed Jan 30, 2025
1 parent b05201b commit 21c0119
Showing 1 changed file with 18 additions and 8 deletions.
26 changes: 18 additions & 8 deletions bin/pygrb/pycbc_pygrb_page_tables
Original file line number Diff line number Diff line change
Expand Up @@ -497,10 +497,11 @@ if onsource_file:
on_trigs_bestnrs = on_trigs['network/reweighted_snr'][...]

# Gather bestNR index
bestNR_event = np.argmax(on_trigs_bestnrs)
loud_on_bestnr_trigs, loud_on_bestnr = \
(on_trigs['network/event_id'][bestNR_event],
on_trigs_bestnrs[bestNR_event])
if on_trigs_bestnrs.size > 0:
bestNR_event = np.argmax(on_trigs_bestnrs)
loud_on_bestnr_trigs, loud_on_bestnr = \
(on_trigs['network/event_id'][bestNR_event],
on_trigs_bestnrs[bestNR_event])
# If the loudest event has bestnr = 0, there is no event at all!
if loud_on_bestnr == 0:
loud_on_bestnr_trigs = None
Expand Down Expand Up @@ -536,9 +537,6 @@ if onsource_file:
[on_trigs[ifo+'/snr'][ifo_trig_index[ifo]] for ifo in ifos] + \
[loud_on_bestnr]
td.append(d)
else:
td.append(["There are no events"] + [0 for number in range(11)] +
[0 for ifo in ifos] + [0])

# Table header
th = ['p-value', 'GPS time', 'Rec. m1', 'Rec. m2', 'Rec. Mc',
Expand All @@ -547,6 +545,10 @@ if onsource_file:

td = list(zip(*td))

# Handle the case in which there is no data to be placed in the table
if not td:
td = [[] for i in np.arange(len(th))]

# Write to h5 file
logging.info("Writing loudest onsource trigger to h5 file.")
with HFile(lont_h5_outfile, 'w') as lont_h5_fp:
Expand All @@ -562,7 +564,11 @@ if onsource_file:
format_strings.extend(['##.##' for ifo in ifos])
format_strings.extend(['##.##'])

# Table data
# Table data: assemble human readable message when no trigger is recovered
if not loud_on_bestnr_trigs:
td = [["There are no events"] + ["-" for number in range(11)] +
["-" for ifo in ifos] + ["-"]]
td = list(zip(*td))
td = [np.asarray(d) for d in td]
html_table = pycbc.results.html_table(td, th,
format_strings=format_strings,
Expand Down Expand Up @@ -646,6 +652,10 @@ if found_missed_file is not None:
len(td))
td = list(zip(*td))

# Handle the case in which there is no data to be placed in the table
if not td:
td = [[] for i in np.arange(len(th))]

# Write to h5 file
with HFile(qf_h5_outfile, 'w') as qf_h5_fp:
for i, key in enumerate(th):
Expand Down

0 comments on commit 21c0119

Please sign in to comment.