From 88bd7df797bed9450313b7a4cb7547c3636f81d5 Mon Sep 17 00:00:00 2001
From: Laurent Mercadier <laurent.mercadier@xfel.eu>
Date: Wed, 19 Feb 2020 13:30:27 +0100
Subject: [PATCH] Find trace using threshold value instead of pulse pattern
---
xgm.py | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/xgm.py b/xgm.py
index d38165e..40d61ea 100644
--- a/xgm.py
+++ b/xgm.py
@@ -882,7 +882,8 @@ def autoFindFastAdcPeaks(data, channel=5, threshold=35000, display=False, plot=F
key = f'FastADC{channel}raw'
if key not in data:
raise ValueError(f'{key} not found in data set')
- trace = data[key].where(data['npulses_sase3']>0, drop=True).isel(trainId=0).values
+ tid = data[key].where(data[key]>threshold, drop=True).trainId[0]
+ trace = data[key].sel(trainId=tid)
centers, peaks = find_peaks(trace, height=threshold, width=(4, None))
c = centers[0]
w = np.average(peaks['widths']).astype(int)
--
GitLab