gstlal_compute_strain: set coherence uncertainty threshold to 0.0025.

lal_insertgap: more informative warning message.

gstlal-calibration/bin/gstlal_compute_strain

@@ -149,7 +149,7 @@ parser.add_option("-v", "--verbose", action = "store_true", help = "Be verbose (
 parser.add_option("--filters-file", metavar="filename", help = "Name of file containing filters (in npz format)")
 parser.add_option("--factors-from-filters-file", action = "store_true", help = "Compute the calibration factors from reference values contained in the filters file instead of from EPICS channels.")
 parser.add_option("--no-coherence", action = "store_true", help = "Gate the calibration factors with a pre-computed coherence channel.")
-parser.add_option("--coherence-uncertainty-threshold", metavar = "float", type = float, default = 0.20, help = "Threshold for the coherence uncertainty for each calibration line. (Default = 0.20)")
+parser.add_option("--coherence-uncertainty-threshold", metavar = "float", type = float, default = 0.0025, help = "Threshold for the coherence uncertainty for each calibration line. (Default = 0.0025)")
 parser.add_option("--coh-unc-sus-line1-channel", metavar="name", default="CAL-CS_TDEP_SUS_LINE1_UNCERTAINTY", help = "Channel name for SUS line 1 coherence uncertainty. (Default=CAL-CS_TDEP_SUS_LINE1_UNCERTAINTY)")
 parser.add_option("--coh-unc-pcaly-line1-channel", metavar="name", default="CAL-CS_TDEP_PCALY_LINE1_UNCERTAINTY", help = "Channel name for PCALY line 1 coherence uncertainty. (Default=CAL-CS_TDEP_PCALY_LINE1_UNCERTAINTY)")
 parser.add_option("--coh-unc-pcaly-line2-channel", metavar="name", default="CAL-CS_TDEP_PCALY_LINE2_UNCERTAINTY", help = "Channel name for PCALY line 2 coherence uncertainty. (Default=CAL-CS_TDEP_PCALY_LINE2_UNCERTAINTY)")
@@ -571,8 +571,8 @@ if not options.no_kappaa or not options.no_kappac or not options.no_fcc or not o
 	exc = calibration_parts.caps_and_progress(pipeline, head_dict["exc"], hoft_caps, "exc")
 
 if not options.no_kappatst or not options.no_kappaa or not options.no_kappac or not options.no_fcc:
-        tstexccaps = "audio/x-raw, format=F64LE, rate=%d" % options.tst_exc_sample_rate
-        tstexc = calibration_parts.caps_and_progress(pipeline, head_dict["tstexc"], tstexccaps, "tstexc")
+	tstexccaps = "audio/x-raw, format=F64LE, rate=%d" % options.tst_exc_sample_rate
+	tstexc = calibration_parts.caps_and_progress(pipeline, head_dict["tstexc"], tstexccaps, "tstexc")
 
 for key in headkeys:
 	if key.startswith("EP"):
@@ -1223,9 +1223,9 @@ if not options.no_fcc:
 mux = pipeparts.mkframecppchannelmux(pipeline, None)
 
 if options.frame_duration is not None:
-        mux.set_property("frame-duration", options.frame_duration)
+	mux.set_property("frame-duration", options.frame_duration)
 if options.frames_per_file is not None:
-        mux.set_property("frames-per-file", options.frames_per_file)
+	mux.set_property("frames-per-file", options.frames_per_file)
 mux.set_property("compression-scheme", options.compression_scheme)
 mux.set_property("compression-level", options.compression_level)
 

gstlal-calibration/gst/lal/gstlal_insertgap.c

@@ -181,8 +181,8 @@ static GstFlowReturn process_inbuf_ ## DTYPE ## COMPLEX(const DTYPE COMPLEX *ind
 	 * First, deal with discontinuity if necessary
 	 */ \
 	if(element->fill_discont && (element->last_sinkbuf_offset_end != 0) && (sinkbuf_pts != element->last_sinkbuf_ets)) { \
-		GST_WARNING_OBJECT(element, "filling in discontinuity of length %" GST_TIME_FORMAT, GST_TIME_ARGS(sinkbuf_pts - element->last_sinkbuf_ets)); \
 		guint64 missing_samples = gst_util_uint64_scale_int_round(sinkbuf_pts - element->last_sinkbuf_ets, element->rate, 1000000000); \
+		GST_WARNING_OBJECT(element, "filling discontinuity of length %f seconds (%lu samples) starting at %f seconds (offset %lu)", (((double) sinkbuf_pts - (double) element->last_sinkbuf_ets) / 1000000000.0), (long unsigned) missing_samples, (double) (element->last_sinkbuf_ets / 1000000000.0), (long unsigned) element->last_sinkbuf_offset_end); \
 		guint standard_blocks = (guint) (missing_samples / max_block_length); \
 		guint64 last_block_length = missing_samples % max_block_length; \
 		DTYPE COMPLEX sample_value; \