PtLUTWriter.cc

Go to the documentation of this file.

#include "L1Trigger/L1TMuonEndCap/interface/PtLUTWriter.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include <fstream>
#include <iostream>
#include <stdexcept>
#include <cassert>

#define PTLUT_SIZE (1<<30)

PtLUTWriter::PtLUTWriter() :
    ptlut_(),
    version_(7), // Initial version, but not hard-coded: gets set by "set_version"
    ok_(false)
{
  ptlut_.reserve(PTLUT_SIZE / 64); // Hack! Hard-code / manually set denom_ for now - AWB 24.05.17
}

PtLUTWriter::~PtLUTWriter() {

}

void PtLUTWriter::write(const std::string& lut_full_path, const uint16_t num_, const uint16_t denom_) const {
  //if (ok_)  return;
  if (not(denom_ == 64)) // Check consistency for temporary hack - AWB 24.05.17
    { edm::LogError("L1T") << "denom_ = " << denom_; return; }

  std::cout << "Writing LUT, this might take a while..." << std::endl;

  std::ofstream outfile(lut_full_path, std::ios::binary);
  if (!outfile.good()) {
    char what[256];
    snprintf(what, sizeof(what), "Fail to open %s", lut_full_path.c_str());
    throw std::invalid_argument(what);
  }

  if (ptlut_.size() != (PTLUT_SIZE / denom_)) {
    char what[256];
    snprintf(what, sizeof(what), "ptlut_.size() is %lu != %i", ptlut_.size(), PTLUT_SIZE);
    throw std::invalid_argument(what);
  }

  if (num_ == 1) 
    ptlut_.at(0) = version_;  // address 0 is the pT LUT version number

  typedef uint64_t full_word_t;
  full_word_t full_word;
  full_word_t sub_word[4] = {0, 0, 0, 0};

  table_t::const_iterator ptlut_it  = ptlut_.begin();
  table_t::const_iterator ptlut_end = ptlut_.end();

  while (ptlut_it != ptlut_end) {
    sub_word[0] = *ptlut_it++;
    sub_word[1] = *ptlut_it++;
    sub_word[2] = *ptlut_it++;
    sub_word[3] = *ptlut_it++;

    full_word = 0;
    full_word |= ((sub_word[0] & 0x1FF) << 0);
    full_word |= ((sub_word[1] & 0x1FF) << 9);
    full_word |= ((sub_word[2] & 0x1FF) << 32);
    full_word |= ((sub_word[3] & 0x1FF) << (32+9));

    outfile.write(reinterpret_cast<char*>(&full_word), sizeof(full_word_t));
  }
  outfile.close();

  //ok_ = true;
  return;
}

void PtLUTWriter::push_back(const content_t& pt) {
  ptlut_.push_back(pt);
}