From 396742e19e7815bfeb743ea0db2e61b7139616b0 Mon Sep 17 00:00:00 2001
From: Marko Mecina <marko.mecina@univie.ac.at>
Date: Wed, 13 Dec 2023 17:22:40 +0100
Subject: [PATCH] integrate events in evtmap
---
Ccs/tools/dataprocessing/smile_raw_ce_converter.py | 5 +++--
1 file changed, 3 insertions(+), 2 deletions(-)
diff --git a/Ccs/tools/dataprocessing/smile_raw_ce_converter.py b/Ccs/tools/dataprocessing/smile_raw_ce_converter.py
index 7f41a3b..9eccf34 100755
--- a/Ccs/tools/dataprocessing/smile_raw_ce_converter.py
+++ b/Ccs/tools/dataprocessing/smile_raw_ce_converter.py
@@ -92,14 +92,15 @@ def mk_ed(data):
# also add an HDU with event map
nodes = np.zeros((2, 2, NROWS_FT, NCOLS_FT))
- nodes[:] = np.nan
for _, _, ccd, col, row, node, fx in bindata:
try:
- nodes[ccd, node, row - 2:row + 3, col - 2:col + 3] = fx.reshape(5, 5)
+ nodes[ccd, node, row - 2:row + 3, col - 2:col + 3] += fx.reshape(5, 5)
except:
print(col, row, 'FAILED')
+ nodes[:] = np.nan
+
ed_img = _assemble_ft_frames_to_fp_view(nodes)
hdl.append(fits.ImageHDU(data=ed_img, name='EVTMAP'))
--
GitLab