EcalTBH2TDCRecInfoProducer.cc

Go to the documentation of this file.

#include "RecoTBCalo/EcalTBTDCReconstructor/interface/EcalTBH2TDCRecInfoProducer.h"
#include "TBDataFormats/HcalTBObjects/interface/HcalTBTiming.h"
#include "TBDataFormats/HcalTBObjects/interface/HcalTBTriggerData.h"
#include "TBDataFormats/EcalTBObjects/interface/EcalTBTDCRecInfo.h"
#include "DataFormats/Common/interface/EDCollection.h"
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"

EcalTBH2TDCRecInfoProducer::EcalTBH2TDCRecInfoProducer(edm::ParameterSet const& ps) {
  rawInfoCollection_ = ps.getParameter<std::string>("rawInfoCollection");
  rawInfoProducer_ = ps.getParameter<std::string>("rawInfoProducer");
  triggerDataCollection_ = ps.getParameter<std::string>("triggerDataCollection");
  triggerDataProducer_ = ps.getParameter<std::string>("triggerDataProducer");
  recInfoCollection_ = ps.getParameter<std::string>("recInfoCollection");

  std::vector<EcalTBH2TDCRecInfoAlgo::EcalTBH2TDCRanges> tdcRanges;

  typedef std::vector<edm::ParameterSet> Parameters;
  Parameters ranges = ps.getParameter<Parameters>("tdcZeros");
  for (Parameters::iterator itRanges = ranges.begin(); itRanges != ranges.end(); ++itRanges) {
    EcalTBH2TDCRecInfoAlgo::EcalTBH2TDCRanges aRange;
    aRange.runRanges.first = itRanges->getParameter<int>("startRun");
    aRange.runRanges.second = itRanges->getParameter<int>("endRun");
    aRange.tdcZero = itRanges->getParameter<double>("tdcZero");
    tdcRanges.push_back(aRange);
  }

  produces<EcalTBTDCRecInfo>(recInfoCollection_);

  algo_ = new EcalTBH2TDCRecInfoAlgo(tdcRanges);
}

EcalTBH2TDCRecInfoProducer::~EcalTBH2TDCRecInfoProducer() {
  if (algo_)
    delete algo_;
}

void EcalTBH2TDCRecInfoProducer::produce(edm::Event& e, const edm::EventSetup& es) {
  int runNumber = e.id().run();
  // Get input
  edm::Handle<HcalTBTiming> ecalRawTDC;
  const HcalTBTiming* ecalTDCRawInfo = nullptr;

  //evt.getByLabel( digiProducer_, digiCollection_, pDigis);
  e.getByLabel(rawInfoProducer_, ecalRawTDC);
  if (ecalRawTDC.isValid()) {
    ecalTDCRawInfo = ecalRawTDC.product();
  }

  if (!ecalTDCRawInfo) {
    edm::LogError("EcalTBTDCRecInfoError") << "Error! can't get the product " << rawInfoCollection_.c_str();
    return;
  }

  // Get input
  edm::Handle<HcalTBTriggerData> triggerData;
  const HcalTBTriggerData* h2TriggerData = nullptr;
  //evt.getByLabel( digiProducer_, digiCollection_, pDigis);
  e.getByLabel(triggerDataProducer_, triggerData);
  if (triggerData.isValid()) {
    h2TriggerData = triggerData.product();
  }

  if (!h2TriggerData) {
    edm::LogError("EcalTBTDCRecInfoError") << "Error! can't get the product " << triggerDataCollection_.c_str();
    return;
  }

  if (!h2TriggerData->wasBeamTrigger()) {
    e.put(std::make_unique<EcalTBTDCRecInfo>(0.5), recInfoCollection_);
  } else {
    e.put(std::make_unique<EcalTBTDCRecInfo>(algo_->reconstruct(runNumber, *ecalRawTDC)), recInfoCollection_);
  }
}