sf-daq/dap
0
0
Fork 1
Code Issues Pull Requests 1 Actions Packages Projects Releases Wiki Activity

comments, todo note

Browse Source
This commit is contained in:
augustin_s
2024-08-01 09:16:40 +02:00
parent ac8967611d
commit 3924abcba4
1 changed files with 6 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

9
dap/algos/peakfind.py
View File

@ -13,8 +13,9 @@ def calc_peakfinder_analysis(results, data, pixel_mask_pf):
if k not in results: if k not in results:
return return
x_beam = results["beam_center_x"] - 0.5 # to coordinates where position of first pixel/point is 0.5, 0.5 # to coordinates where position of first pixel/point is 0.5, 0.5
y_beam = results["beam_center_y"] - 0.5 # to coordinates where position of first pixel/point is 0.5, 0.5 x_beam = results["beam_center_x"] - 0.5
y_beam = results["beam_center_y"] - 0.5
hitfinder_min_snr = results["hitfinder_min_snr"] hitfinder_min_snr = results["hitfinder_min_snr"]
hitfinder_min_pix_count = int(results["hitfinder_min_pix_count"]) hitfinder_min_pix_count = int(results["hitfinder_min_pix_count"])
@ -47,9 +48,11 @@ def calc_peakfinder_analysis(results, data, pixel_mask_pf):
number_of_spots = len(peak_list_x) number_of_spots = len(peak_list_x)
results["number_of_spots"] = number_of_spots results["number_of_spots"] = number_of_spots
# to coordinates where position of first pixel/point is 1, 1
results["spot_x"] = [x + 0.5 for x in peak_list_x] results["spot_x"] = [x + 0.5 for x in peak_list_x]
results["spot_y"] = [y + 0.5 for y in peak_list_y] results["spot_y"] = [y + 0.5 for y in peak_list_y]
results["spot_intensity"] = copy(peak_list_value) results["spot_intensity"] = copy(peak_list_value) #TODO: why is this copy needed?
npeaks_threshold_hit = results.get("npeaks_threshold_hit", 15) npeaks_threshold_hit = results.get("npeaks_threshold_hit", 15)