Gitlab will migrate to a new storage backend starting 0300 UTC on 2020-04-04. We do not anticipate a maintenance window for this migration. Performance may be impacted over the weekend. Thanks for your patience.

Commit f478c017 authored by Aaron Viets's avatar Aaron Viets

calibration_parts.py: remove_lines_with_witnesses(): remove initial samples...

calibration_parts.py:  remove_lines_with_witnesses():  remove initial samples to avoid GSL warnings.
parent ac609c72
Pipeline #50113 passed with stages
in 7 minutes
......@@ -382,10 +382,6 @@ def remove_lines_with_witnesses(pipeline, signal, witnesses, freqs, freq_vars, f
# Find transfer function between witness channel and signal at this frequency
tf_at_f = complex_division(pipeline, line_in_signal, line_in_witness)
# It may be necessary to remove the first few samples since denominator
# samples may have arrived before numerator samples, in which case the
# adder assumes the numerator is one.
tf_at_f = pipeparts.mkgeneric(pipeline, tf_at_f, "lal_insertgap", chop_length = 1000000000) # Removing one second
# Remove worthless data from computation of transfer function if we can
if noisesub_gate_bit is not None:
......@@ -398,10 +394,6 @@ def remove_lines_with_witnesses(pipeline, signal, witnesses, freqs, freq_vars, f
if(i != j):
# Find transfer function between 2 witness channels at this frequency
tf_at_f = complex_division(pipeline, line_in_witnesses[j], line_in_witnesses[i])
# It may be necessary to remove the first few samples since
# denominator samples may have arrived before numerator samples,
# in which case the adder assumes the numerator is one.
tf_at_f = pipeparts.mkgeneric(pipeline, tf_at_f, "lal_insertgap", chop_length = 1000000000)
# Remove worthless data from computation of transfer function if we can
if noisesub_gate_bit is not None:
......@@ -409,6 +401,9 @@ def remove_lines_with_witnesses(pipeline, signal, witnesses, freqs, freq_vars, f
tfs_at_f[(i + 1) * len(witnesses[m]) + j] = pipeparts.mkgeneric(pipeline, tf_at_f, "lal_smoothkappas", default_kappa_re = numpy.random.rand(), default_kappa_im = numpy.random.rand(), array_size = num_median, avg_array_size = num_avg, default_to_median = True, filter_latency = filter_latency)
tfs_at_f = mkinterleave(pipeline, tfs_at_f, complex_data = True)
# It may be necessary to remove data at the beginning so that data that
# is missing do to filtering is not replaced with zeros
tfs_at_f = pipeparts.mkgeneric(pipeline, tfs_at_f, "lal_insertgap", chop_length = 1000000000 * int(1 + (1 - filter_latency) * 21))
tfs_at_f = pipeparts.mkgeneric(pipeline, tfs_at_f, "lal_matrixsolver")
tfs_at_f = mkdeinterleave(pipeline, tfs_at_f, len(witnesses[m]), complex_data = True)
......
......@@ -18,12 +18,10 @@ END = $(shell echo 1234694144 + 4096 + 715 | bc)
#1185771520
SHMRUNTIME = 400
# How much time does the calibration need to settle at the start and end?
PLOT_WARMUP_TIME = 715
#885
PLOT_COOLDOWN_TIME = 715
#4521
PLOT_WARMUP_TIME = 885
PLOT_COOLDOWN_TIME = 4521
GDSCONFIGS = Filters/ER14/GDSFilters/H1GDS_1234630818.ini
GDSCONFIGS = Filters/ER14/GDSFilters/H1GDS_noisesub_test_1232874910.ini
DCSCONFIGS = Filters/ER13/GDSFilters/L1DCS_TEST_1231620000.ini
#../../config_files/O2/H1/tests/H1DCS_AllCorrections_Cleaning.ini
DCSFCCCONFIGS = ../../config_files/O2/H1/tests/H1DCS_FreqIndepAndFccCorrections_Cleaning.ini
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment