added rounding in cluster finder

CMakeLists.txt

@@ -240,9 +240,7 @@ endif()
 add_library(aare_compiler_flags INTERFACE)
 target_compile_features(aare_compiler_flags INTERFACE cxx_std_17)
 
-if(AARE_PYTHON_BINDINGS)
+
-    add_subdirectory(python)
-endif()
 
 #################
 # MSVC specific #
@@ -308,6 +306,9 @@ target_compile_options(
 
 endif() #GCC/Clang specific
 
+if(AARE_PYTHON_BINDINGS)
+    add_subdirectory(python)
+endif()
 
 if(AARE_ASAN)
     message(STATUS "AddressSanitizer enabled")

RELEASE.md

@@ -16,6 +16,7 @@
 ### Bugfixes: 
 
  - multi threaded cluster finder doesnt drop frames if queues are full 
+ - Round before casting in the cluster finder to avoid biasing clusters by truncating
 
 
 ### 2025.11.21

include/aare/ClusterFinder.hpp

@@ -148,12 +148,27 @@ class ClusterFinder {
                         for (int ic = -dx; ic < dx + has_center_pixel_x; ic++) {
                             if (ix + ic >= 0 && ix + ic < frame.shape(1) &&
                                 iy + ir >= 0 && iy + ir < frame.shape(0)) {
-                                CT tmp =
+
-                                    static_cast<CT>(frame(iy + ir, ix + ic)) -
+                                // If the cluster type is an integral type, and
-                                    static_cast<CT>(
+                                // the pedestal is a floating point type then we
+                                // need to round the value before storing it
+                                if constexpr (std::is_integral_v<CT> /*&&
+                                              std::is_floating_point_v<
+                                                  PEDESTAL_TYPE>*/) {
+                                    auto tmp = std::lround(
+                                        frame(iy + ir, ix + ic) -
                                         m_pedestal.mean(iy + ir, ix + ic));
-                                cluster.data[i] =
+                                    cluster.data[i] = static_cast<CT>(tmp);
-                                    tmp; // Watch for out of bounds access
+                                }
+                                // On the other hand if both are floating point
+                                // or both are integral then we can just static
+                                // cast directly
+                                else {
+                                    auto tmp =
+                                        frame(iy + ir, ix + ic) -
+                                        m_pedestal.mean(iy + ir, ix + ic);
+                                    cluster.data[i] = static_cast<CT>(tmp);
+                                }
                             }
                             i++;
                         }

include/aare/ProducerConsumerQueue.hpp

@@ -28,6 +28,7 @@
 #include <atomic>
 #include <cassert>
 #include <cstdlib>
+#include <cstddef>
 #include <memory>
 #include <stdexcept>
 #include <type_traits>