/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_5/src/EventFilter/HcalRawToDigi/plugins/HcalDigiToRaw.cc

Go to the documentation of this file.

#include "DataFormats/FEDRawData/interface/FEDRawDataCollection.h"
#include "DataFormats/FEDRawData/interface/FEDNumbering.h"
#include "DataFormats/HcalDigi/interface/HcalDigiCollections.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()))
{
  produces<FEDRawDataCollection>();
}

// Virtual destructor needed.
HcalDigiToRaw::~HcalDigiToRaw() { }  

// Functions that gets called by framework every event
void HcalDigiToRaw::produce(edm::Event& e, const edm::EventSetup& es)
{
  HcalPacker::Collections colls;

  
  // Step A: Get Inputs 
  edm::Handle<HBHEDigiCollection> hbhe;
  if (!hbheTag_.label().empty()) {
    e.getByLabel(hbheTag_,hbhe);
    colls.hbhe=hbhe.product();
  }
  edm::Handle<HODigiCollection> ho;
  if (!hoTag_.label().empty()) {
    e.getByLabel(hoTag_,ho);
    colls.hoCont=ho.product();
  }
  edm::Handle<HFDigiCollection> hf;
  if (!hfTag_.label().empty()) {
    e.getByLabel(hfTag_,hf);
    colls.hfCont=hf.product();
  }
  edm::Handle<HcalCalibDigiCollection> Calib;
  if (!calibTag_.label().empty()) {
    e.getByLabel(calibTag_,Calib);
    colls.calibCont=Calib.product();
  }
  edm::Handle<ZDCDigiCollection> zdc;
  if (!zdcTag_.label().empty()) {
    e.getByLabel(zdcTag_,zdc);
    colls.zdcCont=zdc.product();
  }
  edm::Handle<HcalTrigPrimDigiCollection> htp;
  if (!trigTag_.label().empty()) {
    e.getByLabel(trigTag_,htp);
    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
  std::auto_ptr<FEDRawDataCollection> raw=std::auto_ptr<FEDRawDataCollection>(new 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(raw);
}