si demonstrator: update HRGNSS SNR settings

siria/si/demonstrator.py

@@ -734,6 +734,11 @@ class DemonstratorConfig(Object):
     create_report = Bool.T(
         default=False,
         help='If ``True``, IDS and Grond reports are generated.')
+    hrgnss_snr_min = Float.T(
+        default=2.,
+        optional=True,
+        help='Minimium absolute Signal-To-Noise amplitude ratio for HRGNSS '
+             'station records.')
 
 
 class SourceInversionDemonstrator(Object):
@@ -954,11 +959,13 @@ class SourceInversionDemonstrator(Object):
                 else:
                     break
 
-            if quantity == 'displacement':
+            if quantity == 'displacement' and \
+                    self.config.hrgnss_snr_min is not None:
+
                 accepted, _ = snr_filter(
                     traces_out,
                     cut_time=self.event.time,
-                    snr_min=3.)
+                    snr_min=self.config.hrgnss_snr_min)
 
                 if not accepted:
                     logger.info(

siria/util.py

@@ -151,8 +151,10 @@ def snr_filter(traces, cut_time, snr_min):
             snr = signal.absmax()[1] / noise.absmax()[1]
             if snr >= snr_min:
                 accepted.append(tr)
+                logger.info('%s accepted with SNR %f', tr.summary, snr)
             else:
                 denied.append(tr)
+                logger.info('%s denied with SNR %f', tr.summary, snr)
         except NoData:
             logger.warning('Did not find sufficient data for %s', tr.summary)