WIBtoWIB2.hpp

Go to the documentation of this file.

 
#ifndef RAWDATAUTILS_INCLUDE_WIBTOWIB2_HPP_
#define RAWDATAUTILS_INCLUDE_WIBTOWIB2_HPP_
 
#include <cstdint>
#include <iostream>
#include <fstream>
#include <filesystem>
#include "fddetdataformats/WIBFrame.hpp"
#include "fddetdataformats/WIB2Frame.hpp"
 
namespace dunedaq {
namespace rawdatautils {
 
fddetdataformats::WIB2Frame
wibtowib2(fddetdataformats::WIBFrame* fr, uint64_t timestamp=0) {
  fddetdataformats::WIB2Frame res;
  for (int i = 0; i < 256; ++i) {
    res.set_adc(i, fr->get_channel(i));
  }
  auto header = fr->get_wib_header();
  res.header.version = header->version;
  res.header.crate = header->crate_no;
  res.header.slot = header->slot_no;
  res.header.link = header->fiber_no;
  res.set_timestamp(timestamp);
  return res;
}
 
void
wib_binary_to_wib2_binary(std::string& filename, std::string& output) {
  std::ifstream file(filename.c_str(), std::ios::binary);
  std::ofstream out(output.c_str(), std::ios::binary);
  std::cout << "Transforming " << filename << " to " << output << '\n';
  auto size = std::filesystem::file_size(filename);
  std::vector<char> v(size);
  file.read(v.data(), size);
  file.close();
  int num_frames = size / sizeof(fddetdataformats::WIBFrame);
  std::cout << "Number of frames found: "<< num_frames << '\n';
  auto ptr = reinterpret_cast<fddetdataformats::WIBFrame*>(v.data());
  uint64_t timestamp = ptr->get_timestamp();
  while(num_frames--){
    auto wib2fr = wibtowib2(ptr, timestamp);
    timestamp += 32;
    ptr++;
    out.write(reinterpret_cast<char*>(&wib2fr), sizeof(wib2fr));
  }
  out.close();
}
 
void
wib_hdf5_to_wib2_binary(std::string& /*filename*/, std::string& /*output*/) {
}
 
 
} // namespace dunedaq::rawdatautils
}
 
#endif // RAWDATAUTILS_INCLUDE_WIBTOWIB2_HPP_