V271 bugfixes (#4972)

* Fix a couple of output/upload value issues in pycbc live (#4953)

* get singles back to being reduced chisquared in uploads

* Save ranking statistic information in the "likelihood" part of the coinc table

* Not all CandidateForGraceDB will have stat in the coinc_results

* Fix bug where single detector event fit_threshold was being read from wrong location (#4971)

* Minor increment to version number

---------

Co-authored-by: Gareth S Cabourn Davies <gareth.cabourndavies@ligo.org>

examples/live/make_singles_significance_fits.py

@@ -21,7 +21,6 @@
 # Set attributes to be something appropriate
 f.attrs['start_date'] = 0
 f.attrs['end_date'] = n_days
-f.attrs['fit_threshold'] = 6.5
 
 f['bins_edges'] = np.logspace(np.log10(min_duration),
                               np.log10(max_duration),
@@ -35,6 +34,7 @@
     ifo_group.attrs['live_time'] = np.round(n_days * 0.7 * 86400.)
     ifo_group.attrs['mean_alpha'] = alpha
     ifo_group.attrs['total_counts'] = daily_counts_per_bin * n_days * n_bins
+    ifo_group.attrs['fit_threshold'] = 6.5
 
     # Make daily fits datasets, this doesn't atually matter for the test, but
     # means that the file format matches in case of future requirements

pycbc/events/single.py

@@ -303,7 +303,8 @@ def check(self, trigs, data_reader):
 
         # fill in a new candidate event
         candidate = {
-            f'foreground/{self.ifo}/{k}': cutall_trigs[k][i] for k in trigs
+            f'foreground/{self.ifo}/{k}': cut_trigs[k][sngl_idx][i]
+            for k in trigs
         }
         candidate['foreground/stat'] = rank[i]
         candidate['foreground/ifar'] = ifar
@@ -319,7 +320,7 @@ def calculate_ifar(self, sngl_ranking, duration):
             with HFile(self.fit_file, 'r') as fit_file:
                 bin_edges = fit_file['bins_edges'][:]
                 live_time = fit_file[self.ifo].attrs['live_time']
-                thresh = fit_file.attrs['fit_threshold']
+                thresh = fit_file[self.ifo].attrs['fit_threshold']
                 dist_grp = fit_file[self.ifo][self.sngl_ifar_est_dist]
                 rates = dist_grp['counts'][:] / live_time
                 coeffs = dist_grp['fit_coeff'][:]

pycbc/io/gracedb.py

@@ -136,7 +136,10 @@ def __init__(self, coinc_ifos, ifos, coinc_results, **kwargs):
         coinc_event_row.time_slide_id = lsctables.TimeSlideID(0)
         coinc_event_row.process_id = proc_id
         coinc_event_row.coinc_event_id = coinc_id
-        coinc_event_row.likelihood = 0.
+        if 'foreground/stat' in coinc_results:
+            coinc_event_row.likelihood = coinc_results['foreground/stat']
+        else:
+            coinc_event_row.likelihood = 0.
         coinc_event_table.append(coinc_event_row)
         outdoc.childNodes[0].appendChild(coinc_event_table)
 

setup.py

@@ -95,7 +95,7 @@ def __getattr__(self, attr):
         vinfo = _version_helper.generate_git_version_info()
     except:
         vinfo = vdummy()
-        vinfo.version = '2.7.0'
+        vinfo.version = '2.7.1'
         vinfo.release = 'True'
 
     version_script = f"""# coding: utf-8