DUNE_FelixAvxReorder_t.cc

Go to the documentation of this file.

#include <fstream>
#include <iostream>
#include <string>
#include <vector>

#include "dune-raw-data/Overlays/FelixFragment.hh"
#include "dune-raw-data/Overlays/FelixReordererFacility.hh"

#pragma GCC diagnostic push
#pragma GCC diagnostic ignored "-Wsign-compare"

#define BOOST_TEST_MODULE(MilliSlice_t)
#include "cetlib/quiet_unit_test.hpp"

BOOST_AUTO_TEST_SUITE(FelixAVXReorder_test)

BOOST_AUTO_TEST_CASE(BaselineTest) {
  // // Get all files.
  // std::vector<int> event_nums = {7258, 7263, 7264, 7269, 7276, 7283, 7284, 7287, 7294, 7296};
  // MC event numbers
  //std::vector<int> event_nums = {3059515, 3059537, 3059542, 3059574, 3059575,
  //                               3059577, 3059599, 3059603, 3059620, 3059622};
  std::vector<std::string> filenames;
  // for (auto event : event_nums) {
  //   for (unsigned p = 0; p < 3; ++p) {
  //     for (unsigned f = 1; f < 10; ++f) {
  //       filenames.push_back(
  //           "/dune/app/users/milov/kevlar/run/uBdat/Run_8700-SubRun_145-Event_" +
  //           std::to_string(event) + "-Plane_" + std::to_string(p) + "-Frame_" +
  //           std::to_string(f) + ".dat");
  //     }
  //     if (p == 2) {
  //       for (unsigned f = 10; f < 14; ++f) {
  //         filenames.push_back(
  //             "/dune/app/users/milov/kevlar/run/uBdat/"
  //             "Run_8700-SubRun_145-Event_" +
  //             std::to_string(event) + "-Plane_" + std::to_string(p) +
  //             "-Frame_" + std::to_string(f) + ".dat");
  //       }
  //     }
  //   }
  // }
  filenames.push_back("/nfs/home/np04daq/milo/frametests/channelid/felix-data-milo.dat");

  for (auto filename : filenames) {
    // Create a regular fragment from file.
    std::ifstream in(filename, std::ios::binary);
    if (!in.is_open()) {
      std::cout << "Could not open file " << filename << ".\n";
      continue;
    }
    std::string contents((std::istreambuf_iterator<char>(in)),
                         (std::istreambuf_iterator<char>()));
    // WARNING: inserting buffer of 32 bytes to correspond with raw data files.
    std::string bufferstring = "abcdefghijklmnopqrstuvwxyzabcdef";
    contents.insert(0, bufferstring);

    dune::FelixFragmentBase::Metadata meta;
    std::unique_ptr<artdaq::Fragment> frag_ptr(artdaq::Fragment::FragmentBytes(
        contents.size(), 1, 1, dune::toFragmentType("FELIX"), meta));
    frag_ptr->resizeBytes(contents.size()-bufferstring.size());
    memcpy(frag_ptr->dataBeginBytes(), contents.c_str(), contents.size()-bufferstring.size());
    in.close();

    dune::FelixFragment flxfrg(*frag_ptr);
    const size_t frames = flxfrg.total_frames();

    // Test whether FelixFragment and FelixFragmentReordered correspond.
    std::cout << "### MEOW -> Reordered FELIX Fragment test.\n";

    dune::ReorderFacility my_reorderer(false/* dune::ReordererType::TypeMILO */);
    artdaq::Fragment reordfrg;
    reordfrg.resizeBytes(frames * dune::FelixReorder::m_num_bytes_per_reord_frame);
    my_reorderer.do_reorder(reordfrg.dataBeginBytes(), frag_ptr->dataBeginBytes() + bufferstring.size(), frames);
    dune::FelixFragment reordflxfrg(reordfrg);

    std::cout << "  -> Total words: " << reordflxfrg.total_words() << '\n';
    std::cout << "  -> Total frames: " << reordflxfrg.total_frames() << '\n';
    std::cout << "  -> Total adc values: " << reordflxfrg.total_adc_values()
              << "\n\n";

    std::cout << "### WOOF -> Reordered WIB frame test.\n";
    std::cout << " -> SOF: " << unsigned(reordflxfrg.sof(0)) << "\n";
    std::cout << " -> Version: " << unsigned(reordflxfrg.version(0)) << "\n";
    std::cout << " -> FiberNo: " << unsigned(reordflxfrg.fiber_no(0)) << "\n";
    std::cout << " -> SlotNo: " << unsigned(reordflxfrg.slot_no(0)) << "\n";
    std::cout << " -> CrateNo: " << unsigned(reordflxfrg.crate_no(0)) << "\n";
    std::cout << " -> Timestamp: " << std::hex << reordflxfrg.timestamp(0)
              << std::dec;
    std::cout << "\n\n";

    std::cout << "### MEOW -> Comparing " << frames << " frames.\n";
    auto compare_begin = std::chrono::high_resolution_clock::now();
    for (unsigned i = 0; i < frames; ++i) {
      BOOST_REQUIRE_EQUAL(flxfrg.sof(i), reordflxfrg.sof(i));
      BOOST_REQUIRE_EQUAL(flxfrg.version(i), reordflxfrg.version(i));
      BOOST_REQUIRE_EQUAL(flxfrg.fiber_no(i), reordflxfrg.fiber_no(i));
      BOOST_REQUIRE_EQUAL(flxfrg.slot_no(i), reordflxfrg.slot_no(i));
      BOOST_REQUIRE_EQUAL(flxfrg.crate_no(i), reordflxfrg.crate_no(i));
      BOOST_REQUIRE_EQUAL(flxfrg.timestamp(i), reordflxfrg.timestamp(i));
      // BOOST_REQUIRE_EQUAL(flxfrg.CRC(i), reordflxfrg.CRC(i));
      for (unsigned j = 0; j < 4; ++j) {
        BOOST_REQUIRE_EQUAL(flxfrg.s1_error(i, j), reordflxfrg.s1_error(i, j));
        BOOST_REQUIRE_EQUAL(flxfrg.s2_error(i, j), reordflxfrg.s2_error(i, j));
        BOOST_REQUIRE_EQUAL(flxfrg.checksum_a(i, j),
                            reordflxfrg.checksum_a(i, j));
        BOOST_REQUIRE_EQUAL(flxfrg.checksum_b(i, j),
                            reordflxfrg.checksum_b(i, j));
        BOOST_REQUIRE_EQUAL(flxfrg.coldata_convert_count(i, j),
                            reordflxfrg.coldata_convert_count(i, j));
        BOOST_REQUIRE_EQUAL(flxfrg.error_register(i, j),
                            reordflxfrg.error_register(i, j));
        for (unsigned h = 0; h < 8; ++h) {
          BOOST_REQUIRE_EQUAL(flxfrg.hdr(i, j, h), reordflxfrg.hdr(i, j, h));
        }
      }
      for (unsigned ch = 0; ch < 256; ++ch) {
        BOOST_REQUIRE_EQUAL(flxfrg.get_ADC(i, ch), reordflxfrg.get_ADC(i, ch));
      }
    }
    auto compare_end = std::chrono::high_resolution_clock::now();
    std::cout << "### MEOW -> Tests successful.\n";
    std::cout << "Took "
              << std::chrono::duration_cast<std::chrono::milliseconds>(
                     compare_end - compare_begin)
                     .count()
              << " ms.\n";

    // // Write reordered fragment to file for compression testing.
    // std::string outname = filename;
    // outname.insert(outname.size() - 4, "_reordered_prev-subtracted");
    // std::ofstream ofile(outname);
    // ofile.write(reinterpret_cast<char const*>(reordflxfrg.dataBeginBytes()),
    //             reordflxfrg.dataSizeBytes());
    // ofile.close();
  } // Loop over files.
}

BOOST_AUTO_TEST_SUITE_END()

#pragma GCC diagnostic pop