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);