ESDigiToRaw.cc

Go to the documentation of this file.

#include "DataFormats/FEDRawData/interface/FEDRawDataCollection.h"
#include "DataFormats/FEDRawData/interface/FEDRawData.h"
#include "DataFormats/FEDRawData/interface/FEDNumbering.h"

#include "DataFormats/EcalDigi/interface/ESDataFrame.h"

#include "EventFilter/ESDigiToRaw/interface/ESDigiToRaw.h"
#include "EventFilter/ESDigiToRaw/src/ESDataFormatterV1_1.h"
#include "EventFilter/ESDigiToRaw/src/ESDataFormatterV4.h"

using namespace std;
using namespace edm;

ESDigiToRaw::ESDigiToRaw(const edm::ParameterSet& ps)
    : ESDataFormatter_(nullptr),
      label_(ps.getParameter<string>("Label")),
      instanceName_(ps.getParameter<string>("InstanceES")),
      ESDigiToken_(consumes<ESDigiCollection>(edm::InputTag(label_, instanceName_))),
      lookup_(ps.getUntrackedParameter<FileInPath>("LookupTable")),
      debug_(ps.getUntrackedParameter<bool>("debugMode", false)),
      formatMajor_(ps.getUntrackedParameter<int>("formatMajor", 4)),
      formatMinor_(ps.getUntrackedParameter<int>("formatMinor", 1)) {
  if (formatMajor_ == 4)
    ESDataFormatter_ = new ESDataFormatterV4(ps);
  else
    ESDataFormatter_ = new ESDataFormatterV1_1(ps);

  produces<FEDRawDataCollection>();
  // initialize look-up table
  for (int i = 0; i < 2; ++i)
    for (int j = 0; j < 2; ++j)
      for (int k = 0; k < 40; ++k)
        for (int m = 0; m < 40; m++)
          fedId_[i][j][k][m] = -1;

  // read in look-up table
  int nLines, iz, ip, ix, iy, fed, kchip, pace, bundle, fiber, optorx;
  ifstream file;
  file.open(lookup_.fullPath().c_str());
  if (file.is_open()) {
    file >> nLines;
    for (int i = 0; i < nLines; ++i) {
      file >> iz >> ip >> ix >> iy >> fed >> kchip >> pace >> bundle >> fiber >> optorx;
      fedId_[(3 - iz) / 2 - 1][ip - 1][ix - 1][iy - 1] = fed;
    }
  } else {
    cout << "[ESDigiToRaw] Look up table file can not be found in " << lookup_.fullPath().c_str() << endl;
  }

  file.close();
}

ESDigiToRaw::~ESDigiToRaw() {
  if (ESDataFormatter_)
    delete ESDataFormatter_;
}

void ESDigiToRaw::produce(edm::StreamID, edm::Event& ev, const edm::EventSetup& es) const {
  ESDataFormatter::Meta_Data meta_data;
  meta_data.lv1 = ev.id().event();
  meta_data.run_number = ev.id().run();
  meta_data.orbit_number = meta_data.lv1 / LHC_BX_RANGE;
  meta_data.bx = (meta_data.lv1 % LHC_BX_RANGE);

  meta_data.kchip_ec = (meta_data.lv1 % KCHIP_EC_RANGE);
  meta_data.kchip_bc = (meta_data.lv1 % KCHIP_BC_RANGE);

  edm::Handle<ESDigiCollection> digis;
  ev.getByToken(ESDigiToken_, digis);

  int ifed;
  ESDataFormatter::Digis Digis;
  Digis.clear();

  for (ESDigiCollection::const_iterator it = digis->begin(); it != digis->end(); ++it) {
    const ESDataFrame& df = *it;
    const ESDetId& detId = it->id();

    ifed = fedId_[(3 - detId.zside()) / 2 - 1][detId.plane() - 1][detId.six() - 1][detId.siy() - 1];
    if (ifed < 0)
      continue;

    Digis[ifed].push_back(df);
  }

  auto productRawData = std::make_unique<FEDRawDataCollection>();

  ESDataFormatter::Digis::const_iterator itfed;
  for (itfed = Digis.begin(); itfed != Digis.end(); ++itfed) {
    int fId = (*itfed).first;

    FEDRawData& fedRawData = productRawData->FEDData(fId);
    ESDataFormatter_->DigiToRaw(fId, Digis, fedRawData, meta_data);

    if (debug_)
      cout << "FED : " << fId << " Data size : " << fedRawData.size() << " (Bytes)" << endl;
  }

  ev.put(std::move(productRawData));

  return;
}