detectors/aare
25
0
Fork 0
mirror of https://github.com/slsdetectorgroup/aare.git synced 2026-06-04 12:18:54 +02:00
Code Issues Actions Packages Releases Activity

with processing

Browse Source
This commit is contained in:
Erik Fröjdh
2026-06-03 15:21:30 +02:00
parent 9113d097c3
commit aeee9b3415
1 changed files with 27 additions and 29 deletions
Download Patch File Download Diff File Expand all files Collapse all files
src/hist/PedestalTrackingPixelHistogram.cpp
+27 -29
View File
@@ -239,35 +239,33 @@ void PedestalTrackingPixelHistogram::worker_loop(int thread_id) {
// never fires, recovering plain histogram-only behaviour
// (modulo the per-pixel gate evaluation). The [xmin, xmax)
// histogram gate lives inside PixelHistogramImpl::fill.
// auto &my_std = partial_std_[thread_id];
// const auto n_sigma = n_sigma_.load(std::memory_order_relaxed);
// for (const auto &frame : *images) {
// const auto cols = frame.shape(1);
// for (int local_row = 0; local_row < local_rows; ++local_row)
// {
// const auto row =
// static_cast<ssize_t>(first_row + local_row);
// for (ssize_t col = 0; col < cols; ++col) {
// const FrameType raw = frame(row, col);
// const AxisType val =
// static_cast<AxisType>(raw) -
// static_cast<AxisType>(my_pedestal.mean(
// static_cast<uint32_t>(local_row),
// static_cast<uint32_t>(col)));
// my_hist.fill_unchecked(local_row,
// static_cast<int>(col),
// val);
// const AxisType sigma = my_std(local_row, col);
// if (sigma > AxisType{0.0} &&
// std::abs(static_cast<AxisType>(val)) <
// n_sigma * sigma) {
// my_pedestal.template push<FrameType>(
// static_cast<uint32_t>(local_row),
// static_cast<uint32_t>(col), raw);
// }
// }
// }
// }
auto &my_std = partial_std_[thread_id];
const auto n_sigma = n_sigma_.load(std::memory_order_relaxed);
for (const auto &frame : *images) {
const auto cols = frame.shape(1);
for (int local_row = 0; local_row < local_rows; ++local_row) {
const auto row =
static_cast<ssize_t>(first_row + local_row);
for (ssize_t col = 0; col < cols; ++col) {
const FrameType raw = frame(row, col);
const AxisType val =
static_cast<AxisType>(raw) -
static_cast<AxisType>(my_pedestal.mean(
static_cast<uint32_t>(local_row),
static_cast<uint32_t>(col)));
my_hist.fill_unchecked(local_row, static_cast<int>(col),
val);
const AxisType sigma = my_std(local_row, col);
if (sigma > AxisType{0.0} &&
std::abs(static_cast<AxisType>(val)) <
n_sigma * sigma) {
my_pedestal.template push<FrameType>(
static_cast<uint32_t>(local_row),
static_cast<uint32_t>(col), raw);
}
}
}
}
break;
}
}