HcalTimeSlewEP.cc

Go to the documentation of this file.

#include "HcalTimeSlewEP.h"
#include "CalibCalorimetry/HcalAlgos/interface/HcalTimeSlew.h"
#include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
#include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
#include "FWCore/ParameterSet/interface/FileInPath.h"
#include "FWCore/Framework/interface/SourceFactory.h"

#include <string>
#include <vector>
#include <map>
#include <iostream>

HcalTimeSlewEP::HcalTimeSlewEP(const edm::ParameterSet& pset) : pset_(pset)
{
  setWhatProduced(this);
  findingRecord<HcalTimeSlewRecord>();
}

HcalTimeSlewEP::~HcalTimeSlewEP()
{}

void HcalTimeSlewEP::setIntervalFor(const edm::eventsetup::EventSetupRecordKey& iKey, const edm::IOVSyncValue& iTime, edm::ValidityInterval& oInterval) {
  oInterval = edm::ValidityInterval(edm::IOVSyncValue::beginOfTime(), edm::IOVSyncValue::endOfTime());
}

void HcalTimeSlewEP::fillDescriptions( edm::ConfigurationDescriptions & descriptions ) {
  edm::ParameterSetDescription desc;

  edm::ParameterSetDescription desc_M2;
  desc_M2.add<double>("tzero");
  desc_M2.add<double>("slope");
  desc_M2.add<double>("tmax");
  std::vector<edm::ParameterSet> default_M2(1);
  desc.addVPSet("timeSlewParametersM2",desc_M2,default_M2);  

  edm::ParameterSetDescription desc_M3;
  desc_M3.add<double>("cap"); 
  desc_M3.add<double>("tspar0");
  desc_M3.add<double>("tspar1");
  desc_M3.add<double>("tspar2");
  desc_M3.add<double>("tspar0_siPM");
  desc_M3.add<double>("tspar1_siPM");
  desc_M3.add<double>("tspar2_siPM");
  std::vector<edm::ParameterSet> default_M3(1);
  desc.addVPSet("timeSlewParametersM3",desc_M3,default_M3);  

  descriptions.addDefault(desc);

}

// ------------ method called to produce the data  ------------
HcalTimeSlewEP::ReturnType
HcalTimeSlewEP::produce(const HcalTimeSlewRecord& iRecord){

  //Two pset sets for M2/Simulation and M3
  std::vector<edm::ParameterSet> p_TimeSlewM2 = pset_.getParameter<std::vector<edm::ParameterSet>>("timeSlewParametersM2");
  std::vector<edm::ParameterSet> p_TimeSlewM3 = pset_.getParameter<std::vector<edm::ParameterSet>>("timeSlewParametersM3");

  ReturnType hcalTimeSlew = std::make_unique<HcalTimeSlew>();

  //loop over the VPSets  
  for(const auto& p_timeslew : p_TimeSlewM2){
    double t0       = p_timeslew.getParameter<double>("tzero");
    double m        = p_timeslew.getParameter<double>("slope");
    double tmaximum = p_timeslew.getParameter<double>("tmax");
    hcalTimeSlew->addM2ParameterSet(t0, m, tmaximum);                     
  }

  for(const auto& p_timeslew : p_TimeSlewM3){
    double cap_         = p_timeslew.getParameter<double>("cap");
    double tspar0_      = p_timeslew.getParameter<double>("tspar0");
    double tspar1_      = p_timeslew.getParameter<double>("tspar1");
    double tspar2_      = p_timeslew.getParameter<double>("tspar2");
    double tspar0_siPM_ = p_timeslew.getParameter<double>("tspar0_siPM");
    double tspar1_siPM_ = p_timeslew.getParameter<double>("tspar1_siPM");
    double tspar2_siPM_ = p_timeslew.getParameter<double>("tspar2_siPM");
    hcalTimeSlew->addM3ParameterSet(cap_,tspar0_,tspar1_,tspar2_,tspar0_siPM_,tspar1_siPM_,tspar2_siPM_);
  }
  
  return hcalTimeSlew;
}

DEFINE_FWK_EVENTSETUP_SOURCE(HcalTimeSlewEP);