HBHEDarkeningEP.cc

Go to the documentation of this file.

#include "HBHEDarkeningEP.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>

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

HBHEDarkeningEP::~HBHEDarkeningEP()
{}

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

void HBHEDarkeningEP::fillDescriptions( edm::ConfigurationDescriptions & descriptions ) {
    edm::ParameterSetDescription desc;
    desc.add<int>("ieta_shift");
    desc.add<double>("drdA");
    desc.add<double>("drdB");
    edm::ParameterSetDescription desc_dosemaps;
    desc_dosemaps.add<int>("energy");
    desc_dosemaps.add<edm::FileInPath>("file");
    std::vector<edm::ParameterSet> default_dosemap(1);
    desc.addVPSet("dosemaps",desc_dosemaps,default_dosemap);
    edm::ParameterSetDescription desc_years;
    desc_years.add<std::string>("year");
    desc_years.add<double>("intlumi");
    desc_years.add<double>("lumirate");
    desc_years.add<int>("energy");
    std::vector<edm::ParameterSet> default_year(1);
    desc.addVPSet("years",desc_years,default_year);

    descriptions.addDefault(desc);
}

// ------------ method called to produce the data  ------------
HBHEDarkeningEP::ReturnType
HBHEDarkeningEP::produce(const HBHEDarkeningRecord& iRecord) {
    //initialize dose maps
    std::vector<edm::ParameterSet> p_dosemaps = pset_.getParameter<std::vector<edm::ParameterSet>>("dosemaps");
    std::map<int,std::vector<std::vector<float>>> dosemaps;
    for(const auto& p_dosemap : p_dosemaps){
        edm::FileInPath fp = p_dosemap.getParameter<edm::FileInPath>("file");
        int file_energy = p_dosemap.getParameter<int>("energy");
        dosemaps.emplace(file_energy,HBHEDarkening::readDoseMap(fp.fullPath()));
    }
    
    //initialize years
    std::vector<edm::ParameterSet> p_years = pset_.getParameter<std::vector<edm::ParameterSet>>("years");
    std::vector<HBHEDarkening::LumiYear> years;
    years.reserve(p_years.size());
    for(const auto& p_year : p_years){
        years.emplace_back(
            p_year.getParameter<std::string>("year"),
            p_year.getParameter<double>("intlumi"),
            p_year.getParameter<double>("lumirate"),
            p_year.getParameter<int>("energy")
        );
    }
  
    ReturnType myResult(
        new HBHEDarkening(
            pset_.getParameter<int>("ieta_shift"),
            pset_.getParameter<double>("drdA"),
            pset_.getParameter<double>("drdB"),
            dosemaps,
            years
        )
    );

    return myResult;
}

DEFINE_FWK_EVENTSETUP_SOURCE(HBHEDarkeningEP);