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

defined values outside interpolation range

Browse Source
This commit is contained in:
stahn
2024-09-23 17:27:00 +02:00
parent 6c2f9dc28b
commit f7bbce9d50
1 changed files with 14 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

15
libeos/file_reader.py
View File

@@ -224,11 +224,24 @@ class AmorData:
# remove 'partially filled' pulses
self.pulseTimeS = self.pulseTimeS[1:-1]
def associate_pulse_with_current(self):
if self.monitorType == 'protonCharge':
self.currentTime -= self.seriesStartTime
currentInterpolator = sp.interpolate.interp1d(self.currentTime, self.current, kind='previous', bounds_error=False)
charge = np.array(currentInterpolator(self.pulseTimeS) * 2*self.tau *1e-3, dtype=float)
# filter low-current pulses
charge = np.where(charge > 2*self.tau * 1e-1, charge, 0)
#self.event
# TODO: activate the following filter AND remove the respective events :
# remove events (wallTime, tof and pixelID) from stream for pulses with too low monitor signal
# probably using np.isin()
def define_monitor(self):
if self.monitorType == 'protonCharge':
# associate each pulse with a proton current
self.currentTime -= self.seriesStartTime
currentInterpolator = sp.interpolate.interp1d(self.currentTime, self.current, kind='previous', fill_value=0)
currentInterpolator = sp.interpolate.interp1d(self.currentTime, self.current, kind='previous', bounds_error=False, fill_value=0)
charge = np.array(currentInterpolator(self.pulseTimeS) * 2*self.tau *1e-3, dtype=float)
# TODO: activate the following filter AND remove the respective events :
# remove events (wallTime, tof and pixelID) from stream for pulses with too low monitor signal