GctDigiToPsbText.cc

Go to the documentation of this file.

#include "GctDigiToPsbText.h"

#include <iomanip>
using std::setfill;
using std::setw;

GctDigiToPsbText::GctDigiToPsbText(const edm::ParameterSet &iConfig)
    : m_gctInputLabel(iConfig.getParameter<edm::InputTag>("GctInputLabel")),
      m_textFileName(iConfig.getParameter<std::string>("TextFileName")),
      m_hexUpperCase(
          iConfig.getUntrackedParameter<bool>("HexUpperCase", false)) {
  for (unsigned i = 0; i < 4; i++) {
    std::stringstream fileStream;
    int ii = (i < 2) ? i : i + 4;
    fileStream << m_textFileName << ii << ".txt";
    std::string fileName(fileStream.str());
    m_file[i].open(fileName.c_str(), std::ios::out);
    if (!m_file[i].good()) {
      throw cms::Exception("GctDigiToPsbTextTextFileOpenError")
          << "GctDigiToPsbText::GctDigiToPsbText : "
          << " couldn't create the file " << fileName << std::endl;
    }
  }
}

GctDigiToPsbText::~GctDigiToPsbText() {
  for (unsigned i = 0; i < 4; i++)
    m_file[i].close();
}

void GctDigiToPsbText::analyze(const edm::Event &iEvent,
                               const edm::EventSetup &iSetup) {

  // static int nevt = -1; nevt++;

  // get digis
  edm::Handle<L1GctEmCandCollection> gctIsolaEm;
  edm::Handle<L1GctEmCandCollection> gctNoIsoEm;
  iEvent.getByLabel(m_gctInputLabel.label(), "isoEm", gctIsolaEm);
  iEvent.getByLabel(m_gctInputLabel.label(), "nonIsoEm", gctNoIsoEm);

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

  std::stringstream sstrm;
  if (m_hexUpperCase)
    sstrm << std::uppercase;
  else
    sstrm.unsetf(std::ios::uppercase);

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

  unsigned int ifile;
  unsigned int cycle;
  unsigned iIsola, iNoIso;

  for (int i = 0; i < 4; i++) {
    cycle = i / 2;
    ifile = i % 2;
    iIsola = ifile + 2;
    iNoIso = ifile;

    // get the data
    data[iIsola][cycle] = gctIsolaEm->at(i).raw();
    data[iNoIso][cycle] = gctNoIsoEm->at(i).raw();

    // print electrons
    sstrm.str("");
    sstrm << setw(4) << setfill('0') << std::hex
          << (data[iIsola][cycle] & 0x7fff) + ((cbs[cycle] & 0x1) << 15);
    m_file[iIsola] << sstrm.str() << std::endl;
    sstrm.str("");
    sstrm << setw(4) << setfill('0') << std::hex
          << (data[iNoIso][cycle] & 0x7fff) + ((cbs[cycle] & 0x1) << 15);
    m_file[iNoIso] << sstrm.str() << std::endl;
  }

  for (unsigned i = 0; i < 4; i++)
    m_file[i] << std::flush;
}