Accounting for LL0 20-bit DAQs (Quick Fix)

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. I plan on updating this further to revise the gps time to better account for when each DAQ was switched and add the current daq configuration of LHO. However, LHO seem to be experiencing less saturation issues than LLO so I want to get this change in sooner rather than later.

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).")
@@ -38,10 +37,6 @@ def find_lpy(event):
         if first_sat:
             sat_channel, sat_time = first_sat.split(' ', 1)
             
-    if config.IFO == 'L1':
-        twentybit = ['SUS-ETM', 'SUS-ITM', 'SUS-BS', 'SUS-MC1', 'SUS-MC2', 'SUS-MC3', 'SUS-PRM', 'SUS-PR2', 'SUS-PR3', 'SUS-SR2', 'SUS-SRM_M2', 'SUS-SRM_M3']
-    if config.IFO == 'L1' and gps > 1358011999 and any(x in sat_channel for x in twentybit):
-        lpy_lim = 4*config.SATURATION_THRESHOLD*4
 
     # generate the LPY mapping and channels to query from base channel
     lpy_map = channel2lpy_coeffs(sat_channel)
@@ -63,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))