GtPsbTextToDigi.cc

Go to the documentation of this file.

#include "GtPsbTextToDigi.h"
// general
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include <iomanip>
#include <iostream>
// gct
#include "DataFormats/L1GlobalCaloTrigger/interface/L1GctCollections.h"
#include "DataFormats/L1GlobalCaloTrigger/interface/L1GctEmCand.h"
#include "DataFormats/L1GlobalCaloTrigger/interface/L1GctEtSums.h"
#include "DataFormats/L1GlobalCaloTrigger/interface/L1GctJetCand.h"
//#include "DataFormats/L1GlobalCaloTrigger/interface/L1GctJetCounts.h"

GtPsbTextToDigi::GtPsbTextToDigi(const edm::ParameterSet &iConfig)
    : m_fileEventOffset(
          iConfig.getUntrackedParameter<int>("FileEventOffset", 0)),
      m_textFileName(iConfig.getParameter<std::string>("TextFileName")),
      m_nevt(0) {
  // Produces collections
  produces<L1GctEmCandCollection>("isoEm");
  produces<L1GctEmCandCollection>("nonIsoEm");
  produces<L1GctJetCandCollection>("cenJets");
  produces<L1GctJetCandCollection>("forJets");
  produces<L1GctJetCandCollection>("tauJets");
  // produces<L1GctEtTotal>();

  // Open the input files
  for (int ifile = 0; ifile < 4; ifile++) {
    // gct em cand coll: (noiso) 0<->0, 1<->1, (iso) 6<->2, 7<->3
    int ii = (ifile < 2) ? ifile : ifile + 4;
    std::stringstream fileStream;
    fileStream << m_textFileName << ii << ".txt";
    std::string fileName(fileStream.str());
    m_file[ifile].open(fileName.c_str(), std::ios::in);
    if (!m_file[ifile].good()) {
      throw cms::Exception("GtPsbTextToDigiTextFileOpenError")
          // LogDebug("GtPsbTextToDigi")
          << "GtPsbTextToDigi::GtPsbTextToDigi : "
          << " couldn't open the file "
          << fileName
          //<< "...skipping!"
          << std::endl;
    }
    std::hex(m_file[ifile]);
  }

  // Initialize bc0 position holder
  for (int i = 0; i < 4; i++)
    m_bc0[i] = -1;
}

GtPsbTextToDigi::~GtPsbTextToDigi() {
  // Close the input files
  for (unsigned i = 0; i < 4; i++) {
    m_file[i].close();
  }
}

void GtPsbTextToDigi::putEmptyDigi(edm::Event &iEvent) {
  std::unique_ptr<L1GctEmCandCollection> gctIsolaEm(
      new L1GctEmCandCollection());
  std::unique_ptr<L1GctEmCandCollection> gctNoIsoEm(
      new L1GctEmCandCollection());
  std::unique_ptr<L1GctJetCandCollection> gctCenJets(
      new L1GctJetCandCollection());
  std::unique_ptr<L1GctJetCandCollection> gctForJets(
      new L1GctJetCandCollection());
  std::unique_ptr<L1GctJetCandCollection> gctTauJets(
      new L1GctJetCandCollection());
  // std::unique_ptr<L1GctEtTotal>           gctEtTotal( new L1GctEtTotal () );
  for (int i = 0; i < 4; i++) {
    gctIsolaEm->push_back(L1GctEmCand(0, true));
    gctNoIsoEm->push_back(L1GctEmCand(0, false));
    gctCenJets->push_back(L1GctJetCand(0, false, false));
    gctForJets->push_back(L1GctJetCand(0, false, true));
    gctTauJets->push_back(L1GctJetCand(0, true, false));
    // gctEtTotal->push_back(());
  }
  iEvent.put(std::move(gctIsolaEm), "isoEm");
  iEvent.put(std::move(gctNoIsoEm), "nonIsoEm");
  iEvent.put(std::move(gctCenJets), "cenJets");
  iEvent.put(std::move(gctForJets), "forJets");
  iEvent.put(std::move(gctTauJets), "tauJets");
  // iEvent.put(std::move(gctEtTotal));

  LogDebug("GtPsbTextToDigi") << "putting empty digi (evt:" << m_nevt << ")\n";
}

void GtPsbTextToDigi::produce(edm::Event &iEvent,
                              const edm::EventSetup &iSetup) {

  // specify clock cycle bit sequence 1 0 1 0... or 0 1 0 1...
  unsigned short cbs[2] = {1, 0};

  // Skip event if required
  if (m_nevt < m_fileEventOffset) {
    putEmptyDigi(iEvent);
    LogDebug("GtPsbTextToDigi")
        << "[GtPsbTextToDigi::produce()] skipping event " << m_nevt
        << std::endl;
    m_nevt++;
    return;
  } else if (m_nevt == 0 && m_fileEventOffset < 0) {
    // skip first fileEventOffset input events
    unsigned long int buff;
    for (int ievt = 0; ievt < abs(m_fileEventOffset); ievt++) {
      for (int ifile = 0; ifile < 4; ifile++) {
        for (int cycle = 0; cycle < 2; cycle++) {
          std::hex(m_file[ifile]);
          m_file[ifile] >> buff;
          unsigned tmp = (buff >> 15) & 0x1;
          if (tmp != cbs[cycle]) {
            if (m_bc0[ifile] == -1 && cycle == 1 && tmp == 1)
              m_bc0[ifile] = ievt;
            else
              throw cms::Exception("GtPsbTextToDigiTextFileFormatError")
                  // std::cout << "GtPsbTextToDigiTextFileFormatError"
                  << "GtPsbTextToDigi::produce : "
                  << " found format inconsistency in file #" << ifile
                  << "\n in skipped line:" << ievt * 2 + 1 << " cycle:" << tmp
                  << " is different from " << cbs[cycle] << std::endl;
          }
        }
      }
      LogDebug("GtPsbTextToDigi")
          << "[GtPsbTextToDigi::produce()] skipping input " << ievt
          << std::endl;
    }
  }
  m_nevt++;

  // New collections
  std::unique_ptr<L1GctEmCandCollection> gctIsolaEm(
      new L1GctEmCandCollection());
  std::unique_ptr<L1GctEmCandCollection> gctNoIsoEm(
      new L1GctEmCandCollection());
  std::unique_ptr<L1GctJetCandCollection> gctCenJets(
      new L1GctJetCandCollection());
  std::unique_ptr<L1GctJetCandCollection> gctForJets(
      new L1GctJetCandCollection());
  std::unique_ptr<L1GctJetCandCollection> gctTauJets(
      new L1GctJetCandCollection());
  // std::unique_ptr<L1GctEtTotal>           gctEtTotal( new L1GctEtTotal () );

  uint16_t data[4][2] = {{0}};
  for (int i = 0; i < 4; i++)
    for (int j = 0; j < 2; j++)
      data[i][j] = 0;

  // Loop over files
  for (int ifile = 0; ifile < 4; ifile++) {
    int ii = (ifile < 2) ? ifile : ifile + 4;

    // Check we're not at the end of the file
    if (m_file[ifile].eof()) {
      LogDebug("GtPsbTextToDigi")
          << "GtPsbTextToDigi::produce : "
          << " unexpected end of file " << m_textFileName << ii << ".txt"
          << std::endl;
      putEmptyDigi(iEvent);
      continue;
    }

    if (!m_file[ifile].good()) {
      LogDebug("GtPsbTextToDigi") << "GtPsbTextToDigi::produce : "
                                  << " problem reading file " << m_textFileName
                                  << ii << ".txt" << std::endl;
      putEmptyDigi(iEvent);
      continue;
    }

    unsigned long int uLongBuffer;

    for (unsigned cycle = 0; cycle < 2; cycle++) {
      m_file[ifile] >> uLongBuffer;
      unsigned tmp = (uLongBuffer >> 15) & 0x1;

      if (false && tmp != cbs[cycle])
        std::cout << "[GtPsbTextToDigi::produce()] asserting "
                  << " evt:" << m_nevt << " ifile:" << ifile
                  << " cycle:" << cbs[cycle] << std::hex
                  << " buffer:" << uLongBuffer << " tmp: " << tmp << std::dec
                  << "\n\n"
                  << std::flush;

      if (tmp != cbs[cycle]) {
        if (m_bc0[ifile] == -1 && cycle == 1 && tmp == 1) {
          m_bc0[ifile] = (m_nevt - m_fileEventOffset);
        } else {
          throw cms::Exception("GtPsbTextToDigiTextFileFormatError")
              // std::cout << "GtPsbTextToDigiTextFileFormatError "
              << "GtPsbTextToDigi::produce : "
              << " found format inconsistency in file #" << ifile
              << "\n in line:" << (m_nevt - m_fileEventOffset) * 2 - 1
              << " cycle:" << tmp << " is different from " << cbs[cycle]
              << std::endl;
        }
      }
      data[ifile][cycle] = (uLongBuffer & 0x7fff);
    } // cycle loop
  }   // ifile loop

  unsigned iIsola, iNoIso;
  for (unsigned cycle = 0; cycle < 2; cycle++) {
    for (unsigned i = 0; i < 2; i++) {
      iIsola = i + 2;
      iNoIso = i;
      gctIsolaEm->push_back(L1GctEmCand(data[iIsola][cycle] & 0x7fff, true));
      gctNoIsoEm->push_back(L1GctEmCand(data[iNoIso][cycle] & 0x7fff, false));
      L1GctEmCand candI(data[iIsola][cycle], true);
      L1GctEmCand candN(data[iNoIso][cycle], false);
    }
  }

  iEvent.put(std::move(gctIsolaEm), "isoEm");
  iEvent.put(std::move(gctNoIsoEm), "nonIsoEm");
  iEvent.put(std::move(gctCenJets), "cenJets");
  iEvent.put(std::move(gctForJets), "forJets");
  iEvent.put(std::move(gctTauJets), "tauJets");
  // iEvent.put(std::move(gctEtTotal));
}

void GtPsbTextToDigi::endJob() {
  int nmem = 4;
  bool match = true;
  for (int i = 0; i < nmem - 1; i++)
    match &= (m_bc0[i] == m_bc0[i + 1]);
  LogDebug("GtPsbTextToDigi") << "[GtPsbTextToDigi::endJob()] ";
  if (!match)
    LogDebug("GtPsbTextToDigi")
        << "did not find matching BC0 in all input files: ";
  else
    LogDebug("GtPsbTextToDigi") << "detected common BC0 in all input files: ";
  for (int i = 0; i < nmem; i++)
    LogDebug("GtPsbTextToDigi") << " " << m_bc0[i];
  LogDebug("GtPsbTextToDigi") << std::flush << std::endl;
}