CTPPSPixelRawToDigi.cc

Go to the documentation of this file.

/*
 * adapting to CTPPS pixel detector March 2017 - F.Ferro
 */

#include "EventFilter/CTPPSRawToDigi/interface/CTPPSPixelRawToDigi.h"

#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "DataFormats/CTPPSDigi/interface/CTPPSPixelDigi.h"

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

#include "EventFilter/CTPPSRawToDigi/interface/CTPPSPixelDataFormatter.h"

#include "FWCore/Framework/interface/ConsumesCollector.h"

using namespace std;

CTPPSPixelRawToDigi::CTPPSPixelRawToDigi(const edm::ParameterSet& conf)
    : config_(conf)

{
  FEDRawDataCollection_ = consumes<FEDRawDataCollection>(config_.getParameter<edm::InputTag>("inputLabel"));

  produces<edm::DetSetVector<CTPPSPixelDigi>>();

  includeErrors_ = config_.getParameter<bool>("includeErrors");
  mappingLabel_ = config_.getParameter<std::string>("mappingLabel");

  if (includeErrors_) {
    produces<edm::DetSetVector<CTPPSPixelDataError>>();
  }
}

CTPPSPixelRawToDigi::~CTPPSPixelRawToDigi() {
  edm::LogInfo("CTPPSPixelRawToDigi") << " CTPPSPixelRawToDigi destructor!";
}

void CTPPSPixelRawToDigi::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
  edm::ParameterSetDescription desc;
  desc.add<bool>("includeErrors", true);
  desc.add<edm::InputTag>("inputLabel", edm::InputTag("rawDataCollector"));
  desc.add<std::string>("mappingLabel", "RPix");
  descriptions.add("ctppsPixelDigis", desc);
}

void CTPPSPixelRawToDigi::produce(edm::Event& ev, const edm::EventSetup& es) {
  edm::Handle<FEDRawDataCollection> buffers;
  ev.getByToken(FEDRawDataCollection_, buffers);

  edm::ESHandle<CTPPSPixelDAQMapping> mapping;

  bool data_exist = false;
  for (int fed = FEDNumbering::MINCTPPSPixelsFEDID; fed <= FEDNumbering::MAXCTPPSPixelsFEDID; fed++) {
    const FEDRawData& tempRawData = buffers->FEDData(fed);
    if (tempRawData.size() != 0) {
      data_exist = true;
      break;
    }
  }
  auto collection = std::make_unique<edm::DetSetVector<CTPPSPixelDigi>>();

  auto errorcollection = std::make_unique<edm::DetSetVector<CTPPSPixelDataError>>();

  if (data_exist) {
    es.get<CTPPSPixelDAQMappingRcd>().get(mapping);

    fedIds_ = mapping->fedIds();

    CTPPSPixelDataFormatter formatter(mapping->ROCMapping);
    formatter.setErrorStatus(includeErrors_);

    bool errorsInEvent = false;
    CTPPSPixelDataFormatter::DetErrors nodeterrors;

    for (auto aFed = fedIds_.begin(); aFed != fedIds_.end(); ++aFed) {
      int fedId = *aFed;

      edm::LogInfo("CTPPSPixelRawToDigi") << " PRODUCE DIGI FOR FED: " << dec << fedId << endl;

      CTPPSPixelDataFormatter::Errors errors;
      const FEDRawData& fedRawData = buffers->FEDData(fedId);

      formatter.interpretRawData(errorsInEvent, fedId, fedRawData, *collection, errors);

      if (includeErrors_) {
        for (auto const& is : errors) {
          uint32_t errordetid = is.first;
          if (errordetid == RPixErrorChecker::dummyDetId) {
            nodeterrors.insert(nodeterrors.end(), errors[errordetid].begin(), errors[errordetid].end());
          } else {
            edm::DetSet<CTPPSPixelDataError>& errorDetSet = errorcollection->find_or_insert(errordetid);
            errorDetSet.data.insert(errorDetSet.data.end(), is.second.begin(), is.second.end());
          }
        }
      }
    }

    if (includeErrors_) {
      errorcollection->find_or_insert(RPixErrorChecker::dummyDetId).data = nodeterrors;
    }
    if (errorsInEvent)
      LogDebug("CTPPSPixelRawToDigi") << "Error words were stored in this event";
  }
  ev.put(std::move(collection));

  if (includeErrors_) {
    ev.put(std::move(errorcollection));
  }
}

DEFINE_FWK_MODULE(CTPPSPixelRawToDigi);