// SPDX-License-Identifier: MPL-2.0
#include "aare/RawFile.hpp"
#include "aare/File.hpp"
#include "aare/RawMasterFile.hpp" //needed for ROI

#include <catch2/catch_test_macros.hpp>
#include <catch2/generators/catch_generators.hpp>
#include <filesystem>

#include "test_config.hpp"
#include "test_macros.hpp"

using aare::File;
using aare::RawFile;
using namespace aare;

TEST_CASE("Read number of frames from a jungfrau raw file", "[.with-data]") {

    auto fpath =
        test_data_path() / "raw/jungfrau/jungfrau_single_master_0.json";
    REQUIRE(std::filesystem::exists(fpath));

    File f(fpath, "r");
    REQUIRE(f.total_frames() == 10);
}

TEST_CASE("Read frame numbers from a jungfrau raw file", "[.with-data]") {
    auto fpath =
        test_data_path() / "raw/jungfrau/jungfrau_single_master_0.json";
    REQUIRE(std::filesystem::exists(fpath));

    File f(fpath, "r");

    // we know this file has 10 frames with frame numbers 1 to 10
    // f0 1,2,3
    // f1 4,5,6
    // f2 7,8,9
    // f3 10
    for (size_t i = 0; i < 10; i++) {
        CHECK(f.frame_number(i) == i + 1);
    }
}

TEST_CASE("Read a frame number too high throws", "[.with-data]") {
    auto fpath =
        test_data_path() / "raw/jungfrau/jungfrau_single_master_0.json";
    REQUIRE(std::filesystem::exists(fpath));

    File f(fpath, "r");

    // we know this file has 10 frames with frame numbers 1 to 10
    // f0 1,2,3
    // f1 4,5,6
    // f2 7,8,9
    // f3 10
    REQUIRE_THROWS(f.frame_number(10));
}

TEST_CASE("Read a frame numbers where the subfile is missing throws",
          "[.with-data]") {
    auto fpath = test_data_path() / "raw/jungfrau" /
                 "jungfrau_missing_subfile_master_0.json";

    REQUIRE(std::filesystem::exists(fpath));

    File f(fpath, "r");

    // we know this file has 10 frames with frame numbers 1 to 10
    // f0 1,2,3
    // f1 4,5,6 - but files f1-f3 are missing
    // f2 7,8,9 - gone
    // f3 10    - gone
    REQUIRE(f.frame_number(0) == 1);
    REQUIRE(f.frame_number(1) == 2);
    REQUIRE(f.frame_number(2) == 3);
    REQUIRE_THROWS(f.frame_number(4));
    REQUIRE_THROWS(f.frame_number(7));
    REQUIRE_THROWS(f.frame_number(937));
    REQUIRE_THROWS(f.frame_number(10));
}

TEST_CASE("Read data from a jungfrau 500k single port raw file",
          "[.with-data]") {
    auto fpath =
        test_data_path() / "raw/jungfrau/jungfrau_single_master_0.json";
    REQUIRE(std::filesystem::exists(fpath));

    File f(fpath, "r");

    // we know this file has 10 frames with pixel 0,0 being: 2123, 2051, 2109,
    // 2117, 2089, 2095, 2072, 2126, 2097, 2102
    std::vector<uint16_t> pixel_0_0 = {2123, 2051, 2109, 2117, 2089,
                                       2095, 2072, 2126, 2097, 2102};
    for (size_t i = 0; i < 10; i++) {
        auto frame = f.read_frame();
        CHECK(frame.rows() == 512);
        CHECK(frame.cols() == 1024);
        CHECK(frame.view<uint16_t>()(0, 0) == pixel_0_0[i]);
    }
}

TEST_CASE("Read frame numbers from a raw file", "[.with-data]") {
    auto fpath =
        test_data_path() / "raw/eiger" / "eiger_500k_16bit_master_0.json";
    REQUIRE(std::filesystem::exists(fpath));

    // we know this file has 3 frames with frame numbers 14, 15, 16
    std::vector<size_t> frame_numbers = {14, 15, 16};

    File f(fpath, "r");
    for (size_t i = 0; i < 3; i++) {
        CHECK(f.frame_number(i) == frame_numbers[i]);
    }
}

TEST_CASE("Compare reading from a numpy file with a raw file", "[.with-data]") {

    SECTION("jungfrau data") {
        auto fpath_raw =
            test_data_path() / "raw/jungfrau" / "jungfrau_single_master_0.json";
        REQUIRE(std::filesystem::exists(fpath_raw));

        auto fpath_npy =
            test_data_path() / "raw/jungfrau" / "jungfrau_single_0.npy";
        REQUIRE(std::filesystem::exists(fpath_npy));

        File raw(fpath_raw, "r");
        File npy(fpath_npy, "r");

        CHECK(raw.total_frames() == 10);
        CHECK(npy.total_frames() == 10);

        for (size_t i = 0; i < 10; ++i) {
            CHECK(raw.tell() == i);
            auto raw_frame = raw.read_frame();
            auto npy_frame = npy.read_frame();
            CHECK((raw_frame.view<uint16_t>() == npy_frame.view<uint16_t>()));
        }
    }

    SECTION("eiger quad data") {
        auto fpath_raw =
            test_data_path() / "raw/eiger_quad_data" /
            "W13_vrpreampscan_m21C_300V_800eV_vthre2000_master_0.json";
        REQUIRE(std::filesystem::exists(fpath_raw));

        auto fpath_npy = test_data_path() / "raw/eiger_quad_data" /
                         "W13_vrpreampscan_m21C_300V_800eV_vthre2000.npy";
        REQUIRE(std::filesystem::exists(fpath_npy));

        File raw(fpath_raw, "r");
        File npy(fpath_npy, "r");

        raw.seek(20);
        auto raw_frame = raw.read_frame();

        auto npy_frame = npy.read_frame();
        CHECK((raw_frame.view<uint32_t>() == npy_frame.view<uint32_t>()));
    }
    SECTION("eiger data") {
        auto fpath_raw = test_data_path() / "raw/eiger" /
                         "Lab6_20500eV_2deg_20240629_master_7.json";
        REQUIRE(std::filesystem::exists(fpath_raw));

        auto fpath_npy =
            test_data_path() / "raw/eiger" / "Lab6_20500eV_2deg_20240629_7.npy";
        REQUIRE(std::filesystem::exists(fpath_npy));

        File raw(fpath_raw, "r");
        File npy(fpath_npy, "r");

        auto raw_frame = raw.read_frame();

        auto npy_frame = npy.read_frame();
        CHECK((raw_frame.view<uint32_t>() == npy_frame.view<uint32_t>()));
    }
}

TEST_CASE("Read multipart files", "[.with-data]") {
    auto fpath =
        test_data_path() / "raw/jungfrau" / "jungfrau_double_master_0.json";
    REQUIRE(std::filesystem::exists(fpath));

    File f(fpath, "r");

    // we know this file has 10 frames check read_multiport.py for the values
    std::vector<uint16_t> pixel_0_0 = {2099, 2121, 2108, 2084, 2084,
                                       2118, 2066, 2108, 2112, 2116};
    std::vector<uint16_t> pixel_0_1 = {2842, 2796, 2865, 2798, 2805,
                                       2817, 2852, 2789, 2792, 2833};
    std::vector<uint16_t> pixel_255_1023 = {2149, 2037, 2115, 2102, 2118,
                                            2090, 2036, 2071, 2073, 2142};
    std::vector<uint16_t> pixel_511_1023 = {3231, 3169, 3167, 3162, 3168,
                                            3160, 3171, 3171, 3169, 3171};
    std::vector<uint16_t> pixel_1_0 = {2748, 2614, 2665, 2629, 2618,
                                       2630, 2631, 2634, 2577, 2598};

    for (size_t i = 0; i < 10; i++) {
        auto frame = f.read_frame();
        CHECK(frame.rows() == 512);
        CHECK(frame.cols() == 1024);
        CHECK(frame.view<uint16_t>()(0, 0) == pixel_0_0[i]);
        CHECK(frame.view<uint16_t>()(0, 1) == pixel_0_1[i]);
        CHECK(frame.view<uint16_t>()(1, 0) == pixel_1_0[i]);
        CHECK(frame.view<uint16_t>()(255, 1023) == pixel_255_1023[i]);
        CHECK(frame.view<uint16_t>()(511, 1023) == pixel_511_1023[i]);
    }
}

struct TestParameters {
    const std::string master_filename{};
    const uint8_t num_ports{};
    const size_t modules_x{};
    const size_t modules_y{};
    const size_t pixels_x{};
    const size_t pixels_y{};
    std::vector<ModuleGeometry> module_geometries{};
};

TEST_CASE("check find_geometry", "[.with-data]") {

    auto test_parameters = GENERATE(
        TestParameters{"raw/jungfrau_2modules_version6.1.2/run_master_0.raw", 2,
                       1, 2, 1024, 1024,
                       std::vector<ModuleGeometry>{
                           ModuleGeometry{0, 0, 512, 1024, 0, 0},
                           ModuleGeometry{0, 512, 512, 1024, 0, 1}}},
        TestParameters{
            "raw/eiger_1_module_version7.0.0/eiger_1mod_master_7.json", 4, 2, 2,
            1024, 512,
            std::vector<ModuleGeometry>{
                ModuleGeometry{0, 0, 256, 512, 0, 0},
                ModuleGeometry{512, 0, 256, 512, 0, 1},
                ModuleGeometry{0, 256, 256, 512, 1, 0},
                ModuleGeometry{512, 256, 256, 512, 1, 1}}},

        TestParameters{"raw/jungfrau_2modules_2interfaces/run_master_0.json", 4,
                       1, 4, 1024, 1024,
                       std::vector<ModuleGeometry>{
                           ModuleGeometry{0, 0, 256, 1024, 0, 0},
                           ModuleGeometry{0, 256, 256, 1024, 1, 0},
                           ModuleGeometry{0, 512, 256, 1024, 2, 0},
                           ModuleGeometry{0, 768, 256, 1024, 3, 0}}},
        TestParameters{
            "raw/eiger_quad_data/"
            "W13_vthreshscan_m21C_300V_800eV_vrpre3400_master_0.json",
            2, 1, 2, 512, 512,
            std::vector<ModuleGeometry>{ModuleGeometry{0, 256, 256, 512, 1, 0},
                                        ModuleGeometry{0, 0, 256, 512, 0, 0}}});

    auto fpath = test_data_path() / test_parameters.master_filename;

    REQUIRE(std::filesystem::exists(fpath));

    RawMasterFile master_file(fpath);

    auto geometry = DetectorGeometry(
        master_file.geometry(), master_file.pixels_x(), master_file.pixels_y(),
        master_file.udp_interfaces_per_module(), master_file.quad());

    CHECK(geometry.modules_x() == test_parameters.modules_x);
    CHECK(geometry.modules_y() == test_parameters.modules_y);
    CHECK(geometry.pixels_x() == test_parameters.pixels_x);
    CHECK(geometry.pixels_y() == test_parameters.pixels_y);

    REQUIRE(geometry.get_module_geometries().size() ==
            test_parameters.num_ports);

    // compare to data stored in header
    RawFile f(fpath, "r");
    for (size_t i = 0; i < test_parameters.num_ports; ++i) {

        auto subfile1_path = f.master().data_fname(i, 0);
        REQUIRE(std::filesystem::exists(subfile1_path));

        auto header = RawFile::read_header(subfile1_path);

        CHECK(header.column == geometry.get_module_geometries(i).col_index);
        CHECK(header.row == geometry.get_module_geometries(i).row_index);

        CHECK(geometry.get_module_geometries(i).height ==
              test_parameters.module_geometries[i].height);
        CHECK(geometry.get_module_geometries(i).width ==
              test_parameters.module_geometries[i].width);

        CHECK(geometry.get_module_geometries(i).origin_x ==
              test_parameters.module_geometries[i].origin_x);
        CHECK(geometry.get_module_geometries(i).origin_y ==
              test_parameters.module_geometries[i].origin_y);
    }
}

TEST_CASE("Open multi module file with ROI", "[.with-data]") {

    auto fpath =
        test_data_path() / "raw/ROITestData/SingleChipROI/Data_master_0.json";
    REQUIRE(std::filesystem::exists(fpath));

    RawFile f(fpath, "r");

    SECTION("read 2 frames") {
        REQUIRE(f.master().roi().value().width() == 256);
        REQUIRE(f.master().roi().value().height() == 256);

        CHECK(f.n_modules() == 2);

        CHECK(f.n_modules_in_roi() == 1);

        auto frames = f.read_n(2);

        CHECK(frames.size() == 2);

        CHECK(frames[0].rows() == 256);
        CHECK(frames[1].cols() == 256);
    }
}

TEST_CASE("Read file with unordered frames", "[.with-data]") {
    // TODO! Better explanation and error message
    auto fpath = test_data_path() / "raw/mythen/scan242_master_3.raw";
    REQUIRE(std::filesystem::exists(fpath));
    File f(fpath);
    REQUIRE_THROWS((f.read_frame()));
}

TEST_CASE("Read Mythenframe", "[.with-data]") {
    auto fpath = test_data_path() / "raw/newmythen03/run_2_master_1.json";
    REQUIRE(std::filesystem::exists(fpath));
    RawFile f(fpath);
    REQUIRE(f.master().roi().value().width() == 2560);
    REQUIRE(f.master().roi().value().height() == 1);
    auto frame = f.read_frame();
    REQUIRE(frame.cols() == 2560);
}