HcalHardcodeParameters.cc

Go to the documentation of this file.

#include "CalibCalorimetry/HcalAlgos/interface/HcalHardcodeParameters.h"

HcalHardcodeParameters::HcalHardcodeParameters(double pedestal, double pedestalWidth, const std::vector<double>& gain, const std::vector<double>& gainWidth, int zsThreshold,
                                               int qieType, const std::vector<double>& qieOffset, const std::vector<double>& qieSlope, int mcShape, int recoShape,
                                               double photoelectronsToAnalog, const std::vector<double>& darkCurrent)
:   pedestal_(pedestal),
    pedestalWidth_(pedestalWidth),
    gain_(gain),
    gainWidth_(gainWidth),
    zsThreshold_(zsThreshold),
    qieType_(qieType),
    qieOffset_(qieOffset),
    qieSlope_(qieSlope),
    mcShape_(mcShape),
    recoShape_(recoShape),
    photoelectronsToAnalog_(photoelectronsToAnalog),
    darkCurrent_(darkCurrent),
    doSipmRadiationDamage_(false)
{
}

HcalHardcodeParameters::HcalHardcodeParameters(const edm::ParameterSet & p)
:   pedestal_(p.getParameter<double>("pedestal")),
    pedestalWidth_(p.getParameter<double>("pedestalWidth")),
    gain_(p.getParameter<std::vector<double>>("gain")),
    gainWidth_(p.getParameter<std::vector<double>>("gainWidth")),
    zsThreshold_(p.getParameter<int>("zsThreshold")),
    qieType_(p.getParameter<int>("qieType")),
    qieOffset_(p.getParameter<std::vector<double>>("qieOffset")),
    qieSlope_(p.getParameter<std::vector<double>>("qieSlope")),
    mcShape_(p.getParameter<int>("mcShape")),
    recoShape_(p.getParameter<int>("recoShape")),
    photoelectronsToAnalog_(p.getParameter<double>("photoelectronsToAnalog")),
    darkCurrent_(p.getParameter<std::vector<double>>("darkCurrent")),
    doSipmRadiationDamage_(p.getParameter<bool>("doRadiationDamage"))
{
    if(doSipmRadiationDamage_) sipmRadiationDamage_ = HcalSiPMRadiationDamage(darkCurrent_,p.getParameter<edm::ParameterSet>("radiationDamage"));
}

const double HcalHardcodeParameters::darkCurrent(unsigned index, double intlumi) const {
    if(doSipmRadiationDamage_ and intlumi>0) return sipmRadiationDamage_.getDarkCurrent(intlumi,index);
    return darkCurrent_.at(index);
}