mirror of
https://github.com/slsdetectorgroup/aare.git
synced 2026-06-05 10:38:41 +02:00
87 lines
3.2 KiB
C++
87 lines
3.2 KiB
C++
// SPDX-License-Identifier: MPL-2.0
|
|
#include "aare/ClusterCollector.hpp"
|
|
#include "aare/ClusterFileSink.hpp"
|
|
#include "aare/ClusterFinder.hpp"
|
|
#include "aare/ClusterFinderMT.hpp"
|
|
#include "aare/ClusterVector.hpp"
|
|
#include "aare/NDView.hpp"
|
|
#include "aare/Pedestal.hpp"
|
|
#include "np_helper.hpp"
|
|
|
|
#include <cstdint>
|
|
#include <filesystem>
|
|
#include <pybind11/pybind11.h>
|
|
#include <pybind11/stl.h>
|
|
#include <pybind11/stl_bind.h>
|
|
|
|
namespace py = pybind11;
|
|
using pd_type = double;
|
|
|
|
using namespace aare;
|
|
|
|
#pragma GCC diagnostic push
|
|
#pragma GCC diagnostic ignored "-Wunused-parameter"
|
|
|
|
template <typename T, uint8_t ClusterSizeX, uint8_t ClusterSizeY,
|
|
typename CoordType = uint16_t>
|
|
void define_ClusterFinder(py::module &m, const std::string &typestr) {
|
|
auto class_name = fmt::format("ClusterFinder_{}", typestr);
|
|
|
|
using ClusterType = Cluster<T, ClusterSizeX, ClusterSizeY, CoordType>;
|
|
|
|
py::class_<ClusterFinder<ClusterType, uint16_t, pd_type>>(
|
|
m, class_name.c_str())
|
|
.def(py::init<Shape<2>, pd_type, size_t>(), py::arg("image_size"),
|
|
py::arg("n_sigma") = 5.0, py::arg("capacity") = 1'000'000)
|
|
|
|
.def_property(
|
|
"nSigma",
|
|
&ClusterFinder<ClusterType, uint16_t, pd_type>::get_nSigma,
|
|
&ClusterFinder<ClusterType, uint16_t, pd_type>::set_nSigma,
|
|
R"(number of sigma above the pedestal to consider a photon during cluster finding.)")
|
|
|
|
.def("push_pedestal_frame",
|
|
[](ClusterFinder<ClusterType, uint16_t, pd_type> &self,
|
|
py::array_t<uint16_t> frame) {
|
|
auto view = make_view_2d(frame);
|
|
self.push_pedestal_frame(view);
|
|
})
|
|
.def("clear_pedestal",
|
|
&ClusterFinder<ClusterType, uint16_t, pd_type>::clear_pedestal)
|
|
.def_property_readonly(
|
|
"pedestal",
|
|
[](ClusterFinder<ClusterType, uint16_t, pd_type> &self) {
|
|
auto pd = new NDArray<pd_type, 2>{};
|
|
*pd = self.pedestal();
|
|
return return_image_data(pd);
|
|
})
|
|
.def_property_readonly(
|
|
"noise",
|
|
[](ClusterFinder<ClusterType, uint16_t, pd_type> &self) {
|
|
auto arr = new NDArray<pd_type, 2>{};
|
|
*arr = self.noise();
|
|
return return_image_data(arr);
|
|
})
|
|
.def(
|
|
"steal_clusters",
|
|
[](ClusterFinder<ClusterType, uint16_t, pd_type> &self,
|
|
bool realloc_same_capacity) {
|
|
ClusterVector<ClusterType> clusters =
|
|
self.steal_clusters(realloc_same_capacity);
|
|
return clusters;
|
|
},
|
|
py::arg("realloc_same_capacity") = false)
|
|
.def(
|
|
"find_clusters",
|
|
[](ClusterFinder<ClusterType, uint16_t, pd_type> &self,
|
|
py::array_t<uint16_t> frame, uint64_t frame_number,
|
|
bool update_pedestal) {
|
|
auto view = make_view_2d(frame);
|
|
self.find_clusters(view, frame_number, update_pedestal);
|
|
return;
|
|
},
|
|
py::arg(), py::arg("frame_number") = 0,
|
|
py::arg("update_pedestal") = true);
|
|
}
|
|
|
|
#pragma GCC diagnostic pop |