HcalDigiToRaw.cc

Go to the documentation of this file.

#include "DataFormats/FEDRawData/interface/FEDRawDataCollection.h"
#include "DataFormats/FEDRawData/interface/FEDNumbering.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "CalibFormats/HcalObjects/interface/HcalDbService.h"
#include "CalibFormats/HcalObjects/interface/HcalDbRecord.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <iostream>
 
#include "EventFilter/HcalRawToDigi/plugins/HcalDigiToRaw.h"
 
using namespace std;
 
HcalDigiToRaw::HcalDigiToRaw(edm::ParameterSet const& conf)
    : hbheTag_(conf.getUntrackedParameter("HBHE", edm::InputTag())),
      hoTag_(conf.getUntrackedParameter("HO", edm::InputTag())),
      hfTag_(conf.getUntrackedParameter("HF", edm::InputTag())),
      zdcTag_(conf.getUntrackedParameter("ZDC", edm::InputTag())),
      calibTag_(conf.getUntrackedParameter("CALIB", edm::InputTag())),
      trigTag_(conf.getUntrackedParameter("TRIG", edm::InputTag())),
      // register for data access
      tok_hbhe_(consumes<HBHEDigiCollection>(hbheTag_)),
      tok_ho_(consumes<HODigiCollection>(hoTag_)),
      tok_hf_(consumes<HFDigiCollection>(hfTag_)),
      tok_calib_(consumes<HcalCalibDigiCollection>(calibTag_)),
      tok_zdc_(consumes<ZDCDigiCollection>(zdcTag_)),
      tok_htp_(consumes<HcalTrigPrimDigiCollection>(trigTag_)) {
  produces<FEDRawDataCollection>();
}
 
// Virtual destructor needed.
HcalDigiToRaw::~HcalDigiToRaw() {}
 
// Functions that gets called by framework every event
void HcalDigiToRaw::produce(edm::StreamID id, edm::Event& e, const edm::EventSetup& es) const {
  HcalPacker::Collections colls;
 
  // Step A: Get Inputs
  edm::Handle<HBHEDigiCollection> hbhe;
  if (!hbheTag_.label().empty()) {
    e.getByToken(tok_hbhe_, hbhe);
    colls.hbhe = hbhe.product();
  }
  edm::Handle<HODigiCollection> ho;
  if (!hoTag_.label().empty()) {
    e.getByToken(tok_ho_, ho);
    colls.hoCont = ho.product();
  }
  edm::Handle<HFDigiCollection> hf;
  if (!hfTag_.label().empty()) {
    e.getByToken(tok_hf_, hf);
    colls.hfCont = hf.product();
  }
  edm::Handle<HcalCalibDigiCollection> Calib;
  if (!calibTag_.label().empty()) {
    e.getByToken(tok_calib_, Calib);
    colls.calibCont = Calib.product();
  }
  edm::Handle<ZDCDigiCollection> zdc;
  if (!zdcTag_.label().empty()) {
    e.getByToken(tok_zdc_, zdc);
    colls.zdcCont = zdc.product();
  }
  edm::Handle<HcalTrigPrimDigiCollection> htp;
  if (!trigTag_.label().empty()) {
    e.getByToken(tok_htp_, htp);
    if (htp.isValid())
      colls.tpCont = htp.product();
  }
  // get the mapping
  edm::ESHandle<HcalDbService> pSetup;
  es.get<HcalDbRecord>().get(pSetup);
  const HcalElectronicsMap* readoutMap = pSetup->getHcalMapping();
  // Step B: Create empty output
  auto raw = std::make_unique<FEDRawDataCollection>();
 
  const int ifed_first = FEDNumbering::MINHCALFEDID;
  const int ifed_last = FEDNumbering::MAXHCALFEDID;
 
  int orbitN = e.id().event();
  int bcnN = 2000;
 
  // Step C: pack all requested FEDs
  for (int ifed = ifed_first; ifed <= ifed_last; ++ifed) {
    FEDRawData& fed = raw->FEDData(ifed);
    try {
      packer_.pack(ifed, ifed - ifed_first, e.id().event(), orbitN, bcnN, colls, *readoutMap, fed);
    } catch (cms::Exception& e) {
      edm::LogWarning("Unpacking error") << e.what();
    } catch (...) {
      edm::LogWarning("Unpacking exception");
    }
  }
 
  e.put(std::move(raw));
}