Accounting for current DAQ configuration (LLO and LHO) after Jan 17 2023

This code takes into account the 20-bit daq configuration for LLO as of the writing of this alog https://alog.ligo-la.caltech.edu/aLOG/index.php?callRep=63133 and the changes made here https://alog.ligo-la.caltech.edu/aLOG/index.php?callRep=63808. The code also accounts for the LHO daq configuration in this alog https://alog.ligo-wa.caltech.edu/aLOG/index.php?callRep=66836. This checks saturation thresholds in saturations.py and lpy.py through a function in config.py. This function needs to be updated further so the gps time better accounts for when each DAQ was switched and it needs to be updated in the event of any future daq changes.

locklost/plugins/lpy.py

@@ -27,7 +27,6 @@ def find_lpy(event):
     infile_csv = 'saturations.csv'
     inpath_csv = event.path(infile_csv)
     sat_channel = ''
-    lpy_lim = 4*config.SATURATION_THRESHOLD
 
     if not os.path.exists(inpath_csv):
         logging.info("LPY plot bypassed (no saturating suspensions).")
@@ -37,6 +36,7 @@ def find_lpy(event):
         first_sat = f.readline()
         if first_sat:
             sat_channel, sat_time = first_sat.split(' ', 1)
+            
 
     # generate the LPY mapping and channels to query from base channel
     lpy_map = channel2lpy_coeffs(sat_channel)
@@ -58,7 +58,8 @@ def find_lpy(event):
         for title in titles:
             for corner, coeff in lpy_map[title].items():
                 idx = channels.index(get_sus_channel(base_channel, corner))
-                lpy[title] += (bufs[idx].data * coeff)/lpy_lim
+                saturation_threshold = config.get_saturation_threshold(get_sus_channel(base_channel, corner), gps, config.IFO)
+                lpy[title] += (bufs[idx].data * coeff)/ (4 * saturation_threshold)
 
         colors = ['#2c7fb8', '#e66101', '#5e3c99']
         fig, axs = plt.subplots(3, sharex=True, figsize=(27, 16))