Several changes to fix line subtraction when lines turn on and off

@louis.dartez pointed out an issue on September 5 when the calibration lines were present in the NOLINES channel during a time when the calibration lines were turned off. This was caused by the fact that, due to the line subtraction gating, the subtraction signal was still being subtracted when the lines were turned off. Here is a spectrogram showing what happened:

In order to fix this, it was necessary to implement a method that distinguishes between zeros and missing data in line witness channels, so that the subtraction can continue to function during brief data dropouts, but actual zeros will cause the subtraction to stop (i.e., to subtract zeros). During these times, it is still important to gate the line subtraction TFs, so that they are not updated with incoherent data that corrupts the subtraction once the lines turn back on. Here is a spectrogram showing the improvement after the first fix:

Although this fixes the issue during the time when the lines are off, there is still a problem visible starting just after the 6-hour mark in the above plot. I discovered that this is because the line subtraction was not updating as intended due to secular changes in the TFs. (When significant secular change is detected in the line subtracion TFs, the algorithm is supposed to switch to a shorter running median so that the subtraction can adjust faster--see details in G2400875.) The reason, as I found, was that the GstController was not working in this case. The fix was to implement another way to control the timing of the line subtraction updates. The the resulting improvement in the spectrongram can be seen here:

Finally, I analyzed another stretch of data in response to this aLOG that revealed another minor issue. In this spectrogram, the problem can be seen from around 19 - 22 UTC. The problem was greatly reduced by fixing the fix to the line subtraction updating due to secular change in the TFs noted above, but there was one remaining issue: the gating was not triggered by the calibration measurements being done around that time because the lock state during calibration measurements (700) is greater than the lock state during low noise (600). So, I included an option in gstlal_compute_strain to require the lock state to be equal to a particular value rather than greater than that value. Here is a spectrogram of that time after the final fix:

tests/check_calibration/Makefile

@@ -10,7 +10,7 @@ OBSRUN = O3
 # Determines whether to find calibration filters in the old SVN or the git repo.  Starting with PreER15 and O4, we switched to the git repo.  Set to 'svn' or 'git'.
 FILTERSRC = svn
 
-START = $(shell echo 1410782418 | bc)
+START = $(shell echo 1409576018 | bc)
 # 1409576018 or 1410782418 gstlal-calibration-1.5.8 linesub issue
 # 1399057926 line sub not working
 # 1371862818 oversubtraction study
@@ -21,7 +21,7 @@ START = $(shell echo 1410782418 | bc)
 # 1238288418
 # 1269090018
 # 1269114645-21449
-END = $(shell echo 1410839118 | bc)
+END = $(shell echo 1409604100 | bc)
 # 1409604100 or 1410839118 gstlal-calibration-1.5.8 linesub issue
 # 1399086062 line sub not working
 # 1371911734 oversubtraction study
@@ -423,7 +423,7 @@ noise_subtraction_ASD_GDS: $(IFO)1_hoft_GDS_frames.cache $(IFO)1_hoft_GDS_frames
 	python3 plot_ASD.py --gps-start-time $(PLOT_START) --gps-end-time $(PLOT_END) --ifo $(IFO)1 --frame-cache-list $(IFO)1_hoft_GDS_frames.cache,$(IFO)1_hoft_GDS_frames_shortTF.cache,$(IFO)1_hoft_GDS_frames_longTF.cache,$(IFO)1_hoft_GDS_frames.cache --channel-list GDS-CALIB_STRAIN,GDS-CALIB_STRAIN_NOLINES,GDS-CALIB_STRAIN_NOLINES,GDS-CALIB_STRAIN_NOLINES --sample-rate 16384 --fft-time 512 --fft-spacing 256 --window=0 --freq-res 0.05 --frequency-min 1 --frequency-max=8192 --ASD-max=1e-17 --labels 'strain,128s\ estimate,4096s\ estimate,switch'
 
 noise_subtraction_spectrogram_GDS: $(IFO)1_hoft_GDS_frames.cache
-	python3 plot_spectrogram.py --gps-start-time $(PLOT_START) --gps-end-time $(PLOT_END) --ifo $(IFO)1 --frame-cache $(IFO)1_hoft_GDS_frames.cache --channel-name GDS-CALIB_STRAIN_NOLINES --sample-rate 16384 --fft-time 1 --fft-spacing 4 --window=3 --frequency-min 10 --frequency-max=1000 --update-time 16 --amplitude-min 0.01 --amplitude-max 100
+	python3 plot_spectrogram.py --gps-start-time $(PLOT_START) --gps-end-time $(PLOT_END) --ifo $(IFO)1 --frame-cache $(IFO)1_hoft_GDS_frames.cache --channel-name GDS-CALIB_STRAIN_NOLINES --sample-rate 16384 --window=3 --frequency-min 10 --frequency-max=1000 --amplitude-min 0.01 --amplitude-max 100
 
 noise_subtraction_ASD_GDS1: $(IFO)1_hoft_GDS_frames.cache $(IFO)1_hoft_GDS_frames_shortTF.cache $(IFO)1_hoft_GDS_frames_longTF.cache
 	python3 plot_ASD.py --gps-start-time $(PLOT_START) --gps-end-time $(PLOT_END) --ifo $(IFO)1 --frame-cache-list $(IFO)1_hoft_GDS_frames.cache,$(IFO)1_hoft_GDS_frames_shortTF.cache,$(IFO)1_hoft_GDS_frames_longTF.cache,$(IFO)1_hoft_GDS_frames.cache --channel-list GDS-CALIB_STRAIN,GDS-CALIB_STRAIN_NOLINES,GDS-CALIB_STRAIN_NOLINES,GDS-CALIB_STRAIN_NOLINES --sample-rate 16384 --fft-time 4096 --fft-spacing 2048 --freq-res 0.001 --frequency-min 410.27 --frequency-max=410.33 --ASD-max=1e-19 --labels 'strain,128s\ estimate,4096s\ estimate,switch'