adjusted weighs filtering further\

swiss-seismological-service · Dec 13, 2024 · 68cc2e6 · 68cc2e6
1 parent 1d2a777
commit 68cc2e6
Showing 1 changed file with 5 additions and 4 deletions.
diff --git a/seismostats/analysis/bvalue/base.py b/seismostats/analysis/bvalue/base.py
@@ -93,10 +93,11 @@ def _filtering(self):
         '''
         Filter out magnitudes below the completeness magnitude.
         '''
-        self.magnitudes = self.magnitudes[self.magnitudes
-                                          >= self.mc - self.delta_m / 2]
-        self.weights = self.weights[self.magnitudes
-                                    >= self.mc - self.delta_m / 2]
+        idx = self.magnitudes >= self.mc - self.delta_m / 2
+        self.magnitudes = self.magnitudes[idx]
+
+        if self.weights is not None:
+            self.weights = self.weights[idx]
 
     def _sanity_checks(self):
         '''