Jungfraujoch/writer/HDF5DataFilePluginROI.cpp

// SPDX-FileCopyrightText: 2024 Filip Leonarski, Paul Scherrer Institute <filip.leonarski@psi.ch>
// SPDX-License-Identifier: GPL-3.0-only

#include "HDF5DataFilePluginROI.h"

void HDF5DataFilePluginROI::OpenFile(HDF5File &data_file, const DataMessage &msg, size_t in_images_per_file) {
    images_per_file = in_images_per_file;

    for (const auto &r: msg.roi) {
        auto &entry = roi_data[r.first];
        entry.max.reserve(images_per_file);
        entry.sum.reserve(images_per_file);
        entry.sum_sq.reserve(images_per_file);
        entry.npixel.reserve(images_per_file);
        entry.x.reserve(images_per_file);
        entry.y.reserve(images_per_file);
    }
}

void HDF5DataFilePluginROI::Write(const DataMessage &msg, uint64_t image_number) {
    for (const auto &r: msg.roi) {
        if (!roi_data.contains(r.first)) {
            roi_data[r.first].max.reserve(images_per_file);
            roi_data[r.first].sum.reserve(images_per_file);
            roi_data[r.first].sum_sq.reserve(images_per_file);
            roi_data[r.first].npixel.reserve(images_per_file);
            roi_data[r.first].x.reserve(images_per_file);
            roi_data[r.first].y.reserve(images_per_file);
        }

        roi_data[r.first].max[image_number] = r.second.max_count;
        roi_data[r.first].sum[image_number] = r.second.sum;
        roi_data[r.first].npixel[image_number] = r.second.pixels;
        roi_data[r.first].sum_sq[image_number] = r.second.sum_square;

        roi_data[r.first].x[image_number] = static_cast<double>(r.second.x_weighted) / static_cast<double>(r.second.sum);
        roi_data[r.first].y[image_number] = static_cast<double>(r.second.y_weighted) / static_cast<double>(r.second.sum);
    }
}

void HDF5DataFilePluginROI::WriteFinal(HDF5File &data_file) {
    if (roi_data.empty())
        return;

    HDF5Group(data_file, "/entry/roi").NXClass("NXcollection");

    for (const auto &r: roi_data) {
        HDF5Group group_roi(data_file, "/entry/roi/" + r.first);
        group_roi.SaveVector("max", r.second.max.vec());
        group_roi.SaveVector("sum", r.second.sum.vec());
        group_roi.SaveVector("sum_sq", r.second.sum_sq.vec());
        group_roi.SaveVector("npixel", r.second.npixel.vec());
        group_roi.SaveVector("x", r.second.x.vec());
        group_roi.SaveVector("y", r.second.y.vec());
    }
}