sinq-reflectometry/eos
0
1
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases 5 Wiki Activity

division by zero in ts mode fixed

Browse Source
This commit is contained in:
stahn
2025-09-05 10:42:31 +02:00
parent f5b4f8bba2
commit 03539eaea6
2 changed files with 5 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files
libeos/file_reader.py
+4 -1
View File
@@ -333,7 +333,10 @@ class AmorData:
self.wallTime_e = self.wallTime_e[filter_e]
logging.info(f' low-beam (<{self.config.lowCurrentThreshold} mC) rejected pulses: {np.shape(self.monitorPerPulse)[0]-1-np.shape(goodTimeS)[0]} out of {np.shape(self.monitorPerPulse)[0]-1}')
logging.info(f' with {np.shape(filter_e)[0]-np.shape(self.tof_e)[0]} events')
logging.info(f' average counts per pulse = {np.shape(self.tof_e)[0] / np.shape(goodTimeS[goodTimeS!=0])[0]:7.1f}')
if np.shape(goodTimeS[goodTimeS!=0])[0]:
logging.info(f' average counts per pulse = {np.shape(self.tof_e)[0] / np.shape(goodTimeS[goodTimeS!=0])[0]:7.1f}')
else:
logging.info(f' average counts per pulse = undefined')
def filter_strange_times(self):
# 'strange' tof times are those with t > 2 tau (originating from the efu)
libeos/reduction.py
+1 -1
View File
@@ -465,7 +465,7 @@ class AmorReduction:
if self.monitor > 1e-6 :
ref_lz *= self.normMonitor / self.monitor
else:
logging.info(' too small monitor value for normalisation -> ignoring monitors')
logging.info(' low monitor -> zero output')
err_lz = ref_lz * np.sqrt( 1/(int_lz+.1) + 1/norm_lz )
# TODO: allow for non-ideal Delta lambda / lambda (rather than 2.2%)