PyCBC 2.8.1__tar.gz → 2.8.3__tar.gz

{PyCBC-2.8.1 → PyCBC-2.8.3}/PKG-INFO

@@ -1,9 +1,9 @@
 Metadata-Version: 2.1
 Name: PyCBC
-Version: 2.8.1
+Version: 2.8.3
 Summary: Core library to analyze gravitational-wave data, find signals, and study their parameters.
 Home-page: http://pycbc.org/
-Download-URL: https://github.com/gwastro/pycbc/tarball/v2.8.1
+Download-URL: https://github.com/gwastro/pycbc/tarball/v2.8.3
 Author: The PyCBC team
 Author-email: alex.nitz@gmail.org
 Keywords: ligo,physics,gravity,signal processing,gravitational waves

{PyCBC-2.8.1 → PyCBC-2.8.3}/PyCBC.egg-info/PKG-INFO

@@ -1,9 +1,9 @@
 Metadata-Version: 2.1
 Name: PyCBC
-Version: 2.8.1
+Version: 2.8.3
 Summary: Core library to analyze gravitational-wave data, find signals, and study their parameters.
 Home-page: http://pycbc.org/
-Download-URL: https://github.com/gwastro/pycbc/tarball/v2.8.1
+Download-URL: https://github.com/gwastro/pycbc/tarball/v2.8.3
 Author: The PyCBC team
 Author-email: alex.nitz@gmail.org
 Keywords: ligo,physics,gravity,signal processing,gravitational waves

{PyCBC-2.8.1 → PyCBC-2.8.3}/PyCBC.egg-info/SOURCES.txt

@@ -318,7 +318,9 @@ examples/multi_inspiral/check_faceon_faceaway_trigs.py
 examples/multi_inspiral/check_gw170817_trigs.py
 examples/multi_inspiral/clean.sh
 examples/multi_inspiral/faceon_faceaway.sh
-examples/multi_inspiral/run.sh
+examples/multi_inspiral/gw170817_h.sh
+examples/multi_inspiral/gw170817_hl.sh
+examples/multi_inspiral/gw170817_hlv.sh
 examples/noise/frequency.py
 examples/noise/timeseries.py
 examples/psd/analytic.py

{PyCBC-2.8.1 → PyCBC-2.8.3}/bin/pycbc_multi_inspiral

@@ -624,7 +624,7 @@ with ctx:
                 if not inj_filter_rejector[ifo].template_segment_checker(
                     bank, t_num, stilde[ifo]
                 ):
-                    logging.info(
+                    logging.debug(
                         "Skipping segment %d/%d with template %d/%d as no "
                         "detectable injection is present",
                         s_num + 1,
@@ -648,7 +648,11 @@ with ctx:
             if not analyse_segment:
                 continue
             logging.info(
-                "Analyzing segment %d/%d", s_num + 1, len(segments[ifo])
+                "Analyzing template %d/%d, segment %d/%d",
+                t_num + 1,
+                n_bank,
+                s_num + 1,
+                len(segments[ifo])
             )
             # The following dicts with IFOs as keys are created to store
             # copies of the matched filtering results computed below.
@@ -664,7 +668,7 @@ with ctx:
             # - The list of normalized SNR values at the trigger locations
             snr = dict.fromkeys(args.instruments)
             for ifo in args.instruments:
-                logging.info(
+                logging.debug(
                     "Filtering template %d/%d, ifo %s", t_num + 1, n_bank, ifo
                 )
                 # The following lines unpack and store copies of the matched
@@ -696,10 +700,10 @@ with ctx:
 
             # Loop over (short) time-slides, starting with the zero-lag
             for slide in range(num_slides):
-                logging.info("Analyzing slide %d/%d", slide, num_slides)
+                logging.debug("Analyzing slide %d/%d", slide, num_slides)
                 # Loop over sky positions
                 for position_index in range(len(sky_positions)):
-                    logging.info(
+                    logging.debug(
                         "Analyzing sky position %d/%d",
                         position_index + 1,
                         len(sky_positions),
@@ -741,7 +745,7 @@ with ctx:
                                 args.instruments[0]
                             ]
                         )
-                    logging.info(
+                    logging.debug(
                         "Found %d coincident triggers", len(coinc_idx)
                     )
                     # Time delay is applied to indices to have them at the IFOs
@@ -753,9 +757,14 @@ with ctx:
                         % len(snr_dict[ifo])
                         for ifo in args.instruments
                     }
-                    # Calculate the coincident and coherent SNR.
-                    # First check there is enough data to compute the SNRs.
-                    if len(coinc_idx) != 0 and nifo > 1:
+                    # Calculate the quadrature sum of individual detector SNRs,
+                    # regardless of how many IFOs are available.
+                    # This accounts for a coherent SNR component and an
+                    # incoherent one.  It is referred to as coincident SNR
+                    # throughout the code, even in the single IFO case.
+                    # The calculation is performed only if triggers were
+                    # gathered so far.
+                    if len(coinc_idx) != 0:
                         # Find coinc SNR at trigger times and apply coinc SNR
                         # threshold (which depopulates coinc_idx accordingly)
                         (
@@ -768,32 +777,24 @@ with ctx:
                             args.coinc_threshold,
                             time_delay_idx[slide][position_index],
                         )
-                        logging.info(
+                        logging.debug(
                             "%d triggers above coincident SNR threshold",
                             len(coinc_idx),
                         )
                         if len(coinc_idx) != 0:
-                            logging.info(
+                            logging.debug(
                                 "With max coincident SNR = %.2f",
                                 max(rho_coinc),
                             )
-                    # If there is only one IFO, just take its triggers
-                    # and their SNRs
-                    elif len(coinc_idx) != 0 and nifo == 1:
-                        coinc_triggers = {
-                            args.instruments[0]: snr_dict[args.instruments[0]][
-                                coinc_idx_det_frame[args.instruments[0]]
-                            ]
-                        }
                     else:
                         coinc_triggers = {}
-                        logging.info("No coincident triggers were found")
+                        logging.debug("No coincident triggers were found")
                     # If there are triggers above coinc threshold and more
-                    # than 2 IFOs, then calculate the coherent statistics for
+                    # than 1 IFO, then calculate the coherent statistics for
                     # them and apply the cut on coherent SNR (with threshold
                     # equal to the coinc SNR one)
-                    if len(coinc_idx) != 0 and nifo > 2:
-                        logging.info("Calculating their coherent statistics")
+                    if len(coinc_idx) != 0 and nifo > 1:
+                        logging.debug("Calculating their coherent statistics")
                         # Plus and cross antenna pattern dictionaries
                         fp = {
                             ifo: antenna_pattern[ifo][position_index][0]
@@ -893,12 +894,12 @@ with ctx:
                                 project,
                                 rho_coinc,
                             )
-                        logging.info(
+                        logging.debug(
                             "%d triggers above coherent SNR threshold",
                             len(rho_coh),
                         )
                         if len(coinc_idx) != 0:
-                            logging.info(
+                            logging.debug(
                                 "With max coherent SNR = %.2f", max(rho_coh)
                             )
                             # Calculate the null SNR and apply the null SNR cut
@@ -918,11 +919,11 @@ with ctx:
                                 snrv=coinc_triggers,
                                 index=coinc_idx,
                             )
-                            logging.info(
+                            logging.debug(
                                 "%d triggers above null threshold", len(null)
                             )
                             if len(coinc_idx) != 0:
-                                logging.info(
+                                logging.debug(
                                     "With max null SNR = %.2f", max(null)
                                 )
                     # Now calculate the individual detector chi2 values
@@ -967,7 +968,7 @@ with ctx:
                             chisq, chisq_dof, coherent_ifo_trigs
                         )
                         # Calculate chisq reweighted SNR
-                        if nifo > 2:
+                        if nifo > 1:
                             reweighted_snr = ranking.newsnr(
                                 rho_coh,
                                 network_chisq_dict,
@@ -983,13 +984,6 @@ with ctx:
                                 null_grad=args.null_grad,
                                 null_step=args.null_step,
                             )
-                        elif nifo == 2:
-                            reweighted_snr = ranking.newsnr(
-                                rho_coinc,
-                                network_chisq_dict,
-                                q=args.chisq_index,
-                                n=args.chisq_nhigh,
-                            )
                         else:
                             rho_sngl = abs(
                                 snr_dict[args.instruments[0]][
@@ -1050,11 +1044,9 @@ with ctx:
                                 ifo_names,
                                 [ifo_out_vals[n] for n in ifo_names],
                             )
-                        if nifo > 2:
+                        if nifo > 1:
                             network_out_vals['coherent_snr'] = rho_coh
                             network_out_vals['null_snr'] = null
-                        elif nifo == 2:
-                            network_out_vals['coherent_snr'] = rho_coinc
                         else:
                             network_out_vals['coherent_snr'] = abs(
                                 snr_dict[args.instruments[0]][
@@ -1089,7 +1081,7 @@ with ctx:
                 cluster_window = int(template.chirp_length * sample_rate)
             # Cluster template events by slide
             for slide in range(num_slides):
-                logging.info("Clustering slide %d", slide)
+                logging.debug("Clustering slide %d", slide)
                 event_mgr.cluster_template_network_events(
                     'time_index', 'reweighted_snr', cluster_window, slide=slide
                 )

{PyCBC-2.8.1 → PyCBC-2.8.3}/bin/pygrb/pycbc_grb_trig_combiner

@@ -137,7 +137,7 @@ def merge_hdf5_files(inputfiles, outputfile, verbose=False, **compression_kw):
     all_event_id_names = sngl_event_id_names + network_event_id_names
 
     # handle search datasets as a special case
-    # (they will the same in all files)
+    # (they will be the same in all files)
     search_datasets = set(filter(lambda x: "/search/" in x, dataset_names))
     gating_datasets = set(filter(lambda x: "/gating/" in x, dataset_names))
     once_only_datasets = search_datasets.union(gating_datasets)

{PyCBC-2.8.1 → PyCBC-2.8.3}/bin/pygrb/pycbc_pygrb_efficiency

@@ -564,6 +564,12 @@ if len(found_after_vetoes['network/template_id']):
 # ==========
 logging.info("Plotting.")
 
+# Calculate distances (as means) for the horizontal axis, only if there are
+# found injections
+if len(found_after_vetoes['network/template_id']):
+    dist_plot_vals = [np.asarray(dist_bin).mean()
+                      for dist_bin in dist_bins]
+
 # Plot efficiency using loudest background
 if opts.background_output_file:
     fig = plt.figure()
@@ -571,10 +577,6 @@ if opts.background_output_file:
 
     # Without found injections, do not determine these quantities to plot
     if len(found_after_vetoes['network/template_id']):
-        # Calculate distances (horizontal axis) as means
-        dist_plot_vals = [np.asarray(dist_bin).mean()
-                          for dist_bin in dist_bins]
-
         # Calculate error bars for efficiency/distance plots and datafiles
         # using max BestNR of background
         yerr_low_mc, yerr_high_mc, fraction_mc = efficiency_with_errs(

{PyCBC-2.8.1 → PyCBC-2.8.3}/bin/pygrb/pycbc_pygrb_plot_injs_results

@@ -474,12 +474,16 @@ ax.grid()
 # Handle axis limits when plotting spins
 if "spin" in x_qty and missed_inj['spin1_a'].size:
     max_missed_inj = missed_inj['spin1_a'].max()
-    ax.set_xlim([0, np.ceil(10 * max(max_missed_inj,
-                                     found_inj[x_qty].max())) / 10])
+    ax.set_xlim([0, np.ceil(10 * max_missed_inj) / 10])
+    if found_inj[x_qty].size:
+        ax.set_xlim([0, np.ceil(10 * max(max_missed_inj,
+                                         found_inj[x_qty].max())) / 10])
 if "spin" in y_qty and missed_inj['spin2_a'].size:
     max_missed_inj = missed_inj['spin2_a'].max()
-    ax.set_ylim([0, np.ceil(10 * max(max_missed_inj,
-                                     found_inj[y_qty].max())) / 10])
+    ax.set_ylim([0, np.ceil(10 * max_missed_inj) / 10])
+    if found_inj[y_qty].size:
+        ax.set_ylim([0, np.ceil(10 * max(max_missed_inj,
+                                         found_inj[y_qty].max())) / 10])
 
 # Handle axis limits when plotting inclination
 if "incl" in x_qty or "incl" in y_qty:

PyCBC-2.8.3/examples/multi_inspiral/check_gw170817_trigs.py

@@ -0,0 +1,46 @@
+#!/usr/bin/env python
+
+# Read a pycbc_multi_inspiral HDF5 trigger file and check that it contains
+# triggers compatible with GW170817
+# 2022 Andrew Williamson, Tito Dal Canton
+
+import sys
+import logging
+import h5py
+
+
+with h5py.File('GW170817_test_output.hdf', 'r') as f:
+    end_time = f['network/end_time_gc'][:]
+    coh_snr = f['network/coherent_snr'][:]
+    rw_snr = f['network/reweighted_snr'][:]
+    slide_id = f['network/slide_id'][:]
+    present_detectors = ''.join(
+        d for d in 'HLV' if f'{d}1' in f.keys()
+    )
+
+# basic sanity checks
+assert (end_time > 0).all()
+assert (coh_snr > 0).all()
+assert (rw_snr > 0).all()
+
+# search for trigs compatible with GW170817
+expected_time = 1187008882.43
+required_net_snr = {'H': 17, 'HL': 28, 'HLV': 28}
+mask = (
+    (abs(expected_time - end_time) < 0.1)
+    & (coh_snr > required_net_snr[present_detectors])
+    & (rw_snr > required_net_snr[present_detectors])
+    & (slide_id == 0)
+)
+n = mask.sum()
+if n > 0:
+    print(
+        f'PASS: GW170817 found with max coherent SNR {max(coh_snr[mask]):.2f}, '
+        f'max reweighted SNR {max(rw_snr[mask]):.2f}'
+    )
+    status = 0
+else:
+    print('FAIL: GW170817 not found')
+    status = 1
+
+sys.exit(status)

PyCBC-2.8.3/examples/multi_inspiral/clean.sh

@@ -0,0 +1 @@
+rm -f bank_veto_bank.xml *.hdf

{PyCBC-2.8.1 → PyCBC-2.8.3}/examples/multi_inspiral/faceon_faceaway.sh

@@ -32,8 +32,8 @@ PAD=8
 BANK_FILE=gw170817_single_template.hdf
 BANK_VETO_FILE=bank_veto_bank.xml
 CHANNEL=SIMULATED_GW170817
-RA=5.016076270234897
-DEC=-0.408407044967
+RA="5.016076270234897 rad"
+DEC="-0.408407044967 rad"
 
 echo -e "\\n\\n>> [`date`] Getting template bank"
 wget -nv -nc ${CONFIG_URL}/${BANK_FILE}

PyCBC-2.8.3/examples/multi_inspiral/gw170817_h.sh

@@ -0,0 +1,62 @@
+#!/bin/bash
+
+set -e
+
+CONFIG_URL=https://github.com/gwastro/pycbc-config/raw/master/test/multi_inspiral
+BANK_FILE=gw170817_single_template.hdf
+echo -e "\\n\\n>> [`date`] Getting template bank"
+wget -nv -nc ${CONFIG_URL}/${BANK_FILE}
+
+EVENT=1187008882
+PAD=8
+START_PAD=111
+END_PAD=17
+GPS_START=$((EVENT - 192 - PAD))
+GPS_END=$((EVENT + 192 + PAD))
+TRIG_START=$((GPS_START + START_PAD))
+TRIG_END=$((GPS_END - END_PAD))
+OUTPUT=GW170817_test_output.hdf
+
+echo -e "\\n\\n>> [`date`] Running pycbc_multi_inspiral on GW170817 data"
+pycbc_multi_inspiral \
+    --verbose \
+    --projection left+right \
+    --instruments H1 \
+    --channel-name H1:GWOSC-16KHZ_R1_STRAIN \
+    --frame-type H1:GWOSC \
+    --trigger-time ${EVENT} \
+    --gps-start-time ${GPS_START} \
+    --gps-end-time ${GPS_END} \
+    --trig-start-time ${TRIG_START} \
+    --trig-end-time ${TRIG_END} \
+    --ra '3.44527994344 rad' \
+    --dec '-0.408407044967 rad' \
+    --bank-file ${BANK_FILE} \
+    --approximant IMRPhenomD \
+    --order -1 \
+    --low-frequency-cutoff 30 \
+    --sngl-snr-threshold 3.0 \
+    --chisq-bins "0.9*get_freq('fSEOBNRv4Peak',params.mass1,params.mass2,params.spin1z,params.spin2z)**(2./3.)" \
+    --pad-data 8 \
+    --strain-high-pass 25 \
+    --sample-rate 4096 \
+    --autogating-threshold 100 \
+    --autogating-cluster 0.5 \
+    --autogating-width 0.25 \
+    --autogating-taper 0.25 \
+    --autogating-pad 0 \
+    --cluster-method window \
+    --cluster-window 0.1 \
+    --segment-length 256 \
+    --segment-start-pad ${START_PAD} \
+    --segment-end-pad ${END_PAD} \
+    --psd-estimation median \
+    --psd-segment-length 32 \
+    --psd-segment-stride 8 \
+    --psd-num-segments 29 \
+    --do-shortslides \
+    --slide-shift 1 \
+    --output ${OUTPUT}
+
+echo -e "\\n\\n>> [`date`] Checking output files"
+python check_gw170817_trigs.py

PyCBC-2.8.3/examples/multi_inspiral/gw170817_hl.sh

@@ -0,0 +1,64 @@
+#!/bin/bash
+
+set -e
+
+CONFIG_URL=https://github.com/gwastro/pycbc-config/raw/master/test/multi_inspiral
+BANK_FILE=gw170817_single_template.hdf
+echo -e "\\n\\n>> [`date`] Getting template bank"
+wget -nv -nc ${CONFIG_URL}/${BANK_FILE}
+
+EVENT=1187008882
+PAD=8
+START_PAD=111
+END_PAD=17
+GPS_START=$((EVENT - 192 - PAD))
+GPS_END=$((EVENT + 192 + PAD))
+TRIG_START=$((GPS_START + START_PAD))
+TRIG_END=$((GPS_END - END_PAD))
+OUTPUT=GW170817_test_output.hdf
+
+echo -e "\\n\\n>> [`date`] Running pycbc_multi_inspiral on GW170817 data"
+pycbc_multi_inspiral \
+    --verbose \
+    --projection left+right \
+    --instruments H1 L1 \
+    --channel-name \
+        H1:GWOSC-16KHZ_R1_STRAIN \
+        L1:GWOSC-16KHZ_R1_STRAIN \
+    --frame-type H1:GWOSC L1:GWOSC \
+    --trigger-time ${EVENT} \
+    --gps-start-time ${GPS_START} \
+    --gps-end-time ${GPS_END} \
+    --trig-start-time ${TRIG_START} \
+    --trig-end-time ${TRIG_END} \
+    --ra '3.44527994344 rad' \
+    --dec '-0.408407044967 rad' \
+    --bank-file ${BANK_FILE} \
+    --approximant IMRPhenomD \
+    --order -1 \
+    --low-frequency-cutoff 30 \
+    --sngl-snr-threshold 3.0 \
+    --chisq-bins "0.9*get_freq('fSEOBNRv4Peak',params.mass1,params.mass2,params.spin1z,params.spin2z)**(2./3.)" \
+    --pad-data 8 \
+    --strain-high-pass 25 \
+    --sample-rate 4096 \
+    --autogating-threshold 100 \
+    --autogating-cluster 0.5 \
+    --autogating-width 0.25 \
+    --autogating-taper 0.25 \
+    --autogating-pad 0 \
+    --cluster-method window \
+    --cluster-window 0.1 \
+    --segment-length 256 \
+    --segment-start-pad ${START_PAD} \
+    --segment-end-pad ${END_PAD} \
+    --psd-estimation median \
+    --psd-segment-length 32 \
+    --psd-segment-stride 8 \
+    --psd-num-segments 29 \
+    --do-shortslides \
+    --slide-shift 1 \
+    --output ${OUTPUT}
+
+echo -e "\\n\\n>> [`date`] Checking output files"
+python check_gw170817_trigs.py

PyCBC-2.8.1/examples/multi_inspiral/run.sh → PyCBC-2.8.3/examples/multi_inspiral/gw170817_hlv.sh

@@ -1,25 +1,11 @@
-#!/bin/bash -e
+#!/bin/bash
+
+set -e
 
 CONFIG_URL=https://github.com/gwastro/pycbc-config/raw/master/test/multi_inspiral
 BANK_FILE=gw170817_single_template.hdf
-BANK_VETO_FILE=bank_veto_bank.xml
-H1_FRAME=https://www.gwosc.org/eventapi/html/GWTC-1-confident/GW170817/v3/H-H1_GWOSC_4KHZ_R1-1187006835-4096.gwf
-H1_CHANNEL=GWOSC-4KHZ_R1_STRAIN
-L1_FRAME=https://dcc.ligo.org/public/0144/T1700406/003/L-L1_CLEANED_HOFT_C02_T1700406_v3-1187008667-4096.gwf
-L1_CHANNEL=DCH-CLEAN_STRAIN_C02_T1700406_v3
-V1_FRAME=https://www.gwosc.org/eventapi/html/GWTC-1-confident/GW170817/v3/V-V1_GWOSC_4KHZ_R1-1187006835-4096.gwf
-V1_CHANNEL=GWOSC-4KHZ_R1_STRAIN
-
-for IFO in H1 L1 V1; do
-    echo -e "\\n\\n>> [`date`] Getting ${IFO} frame"
-    FRAME=${IFO}_FRAME
-    wget -nv -nc ${!FRAME}
-done
-
 echo -e "\\n\\n>> [`date`] Getting template bank"
 wget -nv -nc ${CONFIG_URL}/${BANK_FILE}
-echo -e "\\n\\n>> [`date`] Bank veto bank"
-wget -nv -nc ${CONFIG_URL}/${BANK_VETO_FILE}
 
 EVENT=1187008882
 PAD=8
@@ -36,13 +22,18 @@ pycbc_multi_inspiral \
     --verbose \
     --projection left+right \
     --instruments H1 L1 V1 \
+    --channel-name \
+        H1:GWOSC-16KHZ_R1_STRAIN \
+        L1:GWOSC-16KHZ_R1_STRAIN \
+        V1:GWOSC-16KHZ_R1_STRAIN \
+    --frame-type H1:GWOSC L1:GWOSC V1:GWOSC \
     --trigger-time ${EVENT} \
     --gps-start-time ${GPS_START} \
     --gps-end-time ${GPS_END} \
     --trig-start-time ${TRIG_START} \
     --trig-end-time ${TRIG_END} \
-    --ra 3.44527994344 \
-    --dec -0.408407044967 \
+    --ra '3.44527994344 rad' \
+    --dec '-0.408407044967 rad' \
     --bank-file ${BANK_FILE} \
     --approximant IMRPhenomD \
     --order -1 \
@@ -52,11 +43,11 @@ pycbc_multi_inspiral \
     --pad-data 8 \
     --strain-high-pass 25 \
     --sample-rate 4096 \
-    --channel-name H1:${H1_CHANNEL} L1:${L1_CHANNEL} V1:${V1_CHANNEL} \
-    --frame-files \
-        H1:`basename ${H1_FRAME}` \
-        L1:`basename ${L1_FRAME}` \
-        V1:`basename ${V1_FRAME}` \
+    --autogating-threshold 100 \
+    --autogating-cluster 0.5 \
+    --autogating-width 0.25 \
+    --autogating-taper 0.25 \
+    --autogating-pad 0 \
     --cluster-method window \
     --cluster-window 0.1 \
     --segment-length 256 \

{PyCBC-2.8.1 → PyCBC-2.8.3}/pycbc/filter/resample.py

@@ -128,7 +128,7 @@ def lfilter(coefficients, timeseries):
     else:
         # recursively perform which saves a bit on memory usage
         # but must keep within recursion limit
-        chunksize = max(fillen * 5, len(timeseries) // 128)
+        chunksize = max(fillen * 5, len(timeseries) // 2)
         part1 = lfilter(coefficients, timeseries[0:chunksize])
         part2 = lfilter(coefficients, timeseries[chunksize - fillen:])
         out = timeseries.copy()

{PyCBC-2.8.1 → PyCBC-2.8.3}/pycbc/pnutils.py

@@ -1025,6 +1025,8 @@ def jframe_to_l0frame(mass1, mass2, f_ref, phiref=0., thetajn=0., phijl=0.,
         binary (in solar masses)
     f_ref : float
         The reference frequency.
+    phiref : float
+        The reference phase.
     thetajn : float
         Angle between the line of sight and the total angular momentume J.
     phijl : float

{PyCBC-2.8.1 → PyCBC-2.8.3}/pycbc/results/pygrb_postprocessing_utils.py

@@ -345,7 +345,7 @@ def load_data(input_file, ifos, rw_snr_threshold=None, data_tag=None,
             # The dataset contains search information or missed injections
             # information, not properties of triggers or found injections:
             # just copy it
-            if 'search' in path or 'missed' in path:
+            if 'search' in path or 'missed' in path or 'gating' in path:
                 trigs_dict[path] = dset[:]
             # The dataset is trig/inj info at an IFO:
             # cut with the correct index
@@ -401,9 +401,12 @@ def apply_vetoes_to_found_injs(found_missed_file, found_injs, ifos,
     keep_keys = keys if keys else found_injs.keys()
 
     if not found_missed_file or ifos[0]+'/end_time' not in found_injs.keys():
+        t_id_key = 'network/template_id'
+        if t_id_key not in keep_keys:
+            keep_keys = list(keep_keys+[t_id_key])
         empty_dict = dict.fromkeys(keep_keys, numpy.array([]))
-        empty_dict['network/template_id'] = \
-            empty_dict['network/template_id'].astype(dtype=numpy.int64)
+        empty_dict[t_id_key] = \
+            empty_dict[t_id_key].astype(dtype=numpy.int64)
         return (empty_dict, empty_dict, None, None)
 
     found_idx = numpy.arange(len(found_injs[ifos[0]+'/end_time'][:]))

{PyCBC-2.8.1 → PyCBC-2.8.3}/pycbc/results/scatter_histograms.py

@@ -385,6 +385,7 @@ def create_marginalized_hist(ax, values, label, percentiles=None,
     """
     if fillcolor is None:
         htype = 'step'
+        fillcolor = 'none'
     else:
         htype = 'stepfilled'
     if rotated:

{PyCBC-2.8.1 → PyCBC-2.8.3}/pycbc/version.py

@@ -1,27 +1,27 @@
 # coding: utf-8
-# Generated by setup.py for PyCBC on 2025-04-02 08:22:20 +0000.
+# Generated by setup.py for PyCBC on 2025-05-13 08:01:03 +0000.
 
 # general info
-version = '2.8.1'
-date = '2025-04-01 23:56:55 +0000'
+version = '2.8.3'
+date = '2025-05-13 06:43:32 +0000'
 release = True
-last_release = '2.8.1'
+last_release = '2.8.3'
 
 # git info
-git_hash = '2daede81e4763d03c723e9f32492fa602b896914'
+git_hash = 'd5f6fe7b8ff2f8dc98ff10d4924ac814b62dfdf8'
 git_branch = 'None'
-git_tag = 'v2.8.1'
+git_tag = 'v2.8.3'
 git_author = 'Tito Dal Canton <tito.dalcanton@ijclab.in2p3.fr>'
 git_committer = 'GitHub <noreply@github.com>'
 git_status = 'CLEAN: All modifications committed'
 git_builder = 'Unknown User <>'
-git_build_date = '2025-04-02 08:22:20 +0000'
-git_verbose_msg = """Version: 2.8.1
+git_build_date = '2025-05-13 08:01:03 +0000'
+git_verbose_msg = """Version: 2.8.3
 Branch: None
-Tag: v2.8.1
-Id: 2daede81e4763d03c723e9f32492fa602b896914
+Tag: v2.8.3
+Id: d5f6fe7b8ff2f8dc98ff10d4924ac814b62dfdf8
 Builder: Unknown User <>
-Build date: 2025-04-02 08:22:20 +0000
+Build date: 2025-05-13 08:01:03 +0000
 Repository status is CLEAN: All modifications committed"""
 
 from pycbc._version import *