Merge pull request #51 from LIBRA-project/patch-window-avg

Bugfix wrong array for avg count rate calculation

libra_toolbox/neutron_detection/diamond/process_data.py

@@ -113,10 +113,10 @@ def get_avg_rate(self, t_min: float, t_max: float, energy_window: tuple = None):
             time_values = time_values[peak_mask]
 
         # Create mask to only count pulses of any energy in section time window
-        idx = np.logical_and(self.time_values > t_min, self.time_values < t_max)
+        idx = np.logical_and(time_values > t_min, time_values < t_max)
 
-        counts = len(self.time_values[idx])
-        error = np.sqrt(len(self.time_values[idx]))
+        counts = len(time_values[idx])
+        error = np.sqrt(len(time_values[idx]))
         delta_t = t_max - t_min
         count_rate = counts / delta_t
         count_rate_err = error / delta_t