Line Code
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64
/**
 *   
 * \author G. Franzoni
 *
 */

#include "CaloOnlineTools/EcalTools/plugins/EcalHexDisplay.h"

EcalHexDisplay::EcalHexDisplay(const edm::ParameterSet& ps)
    : verbosity_(ps.getUntrackedParameter<int>("verbosity", 1)),
      beg_fed_id_(ps.getUntrackedParameter<int>("beg_fed_id", 0)),
      end_fed_id_(ps.getUntrackedParameter<int>("end_fed_id", 654)),
      first_event_(ps.getUntrackedParameter<int>("first_event", 1)),
      last_event_(ps.getUntrackedParameter<int>("last_event", 9999999)),
      event_(0),
      writeDcc_(ps.getUntrackedParameter<bool>("writeDCC", false)),
      filename_(ps.getUntrackedParameter<std::string>("filename", "dump.bin")),
      fedRawDataCollectionToken_(
          consumes<FEDRawDataCollection>(ps.getParameter<edm::InputTag>("fedRawDataCollectionTag"))) {}

void EcalHexDisplay::analyze(const edm::Event& e, const edm::EventSetup& c) {
  event_++;
  if (event_ < first_event_ || last_event_ < event_)
    return;

  const edm::Handle<FEDRawDataCollection>& rawdata = e.getHandle(fedRawDataCollectionToken_);

  std::ofstream dumpFile(filename_.c_str(), std::ios::app);

  for (int id = 0; id <= FEDNumbering::MAXFEDID; ++id) {
    if (id < beg_fed_id_ || end_fed_id_ < id)
      continue;

    const FEDRawData& data = rawdata->FEDData(id);

    if (data.size() > 4) {
      edm::LogVerbatim("EcalTools") << "\n\n\n[EcalHexDumperModule] Event: " << std::dec << event_ << " fed_id: " << id
                                    << " size_fed: " << data.size() << "\n";

      if ((data.size() % 16) != 0) {
        edm::LogVerbatim("EcalTools") << "***********************************************";
        edm::LogVerbatim("EcalTools") << "Fed size in bits not multiple of 64, strange.";
        edm::LogVerbatim("EcalTools") << "***********************************************";
      }

      int length = data.size();
      const unsigned long* pData = (reinterpret_cast<unsigned long*>(const_cast<unsigned char*>(data.data())));
      edm::LogVerbatim("EcalTools") << std::setfill('0');
      for (int words = 0; words < length / 4; (words += 2)) {
        edm::LogVerbatim("EcalTools") << std::setw(8) << std::hex << pData[words + 1] << " " << std::setw(8) << std::hex
                                      << pData[words];
      }

      edm::LogVerbatim("EcalTools") << "\n";

      if (beg_fed_id_ <= id && id <= end_fed_id_ && writeDcc_) {
        dumpFile.write(reinterpret_cast<const char*>(pData), length);
      }
    }
  }
  dumpFile.close();
  if (!writeDcc_)
    remove(filename_.c_str());
}