parseHBHEPhase1AlgoDescription.cc File Reference

#include <cassert>
#include "RecoLocalCalo/HcalRecAlgos/interface/parseHBHEPhase1AlgoDescription.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "RecoLocalCalo/HcalRecAlgos/interface/PulseShapeFitOOTPileupCorrection.h"
#include "RecoLocalCalo/HcalRecAlgos/interface/HcalDeterministicFit.h"
#include "RecoLocalCalo/HcalRecAlgos/interface/SimpleHBHEPhase1Algo.h"

Go to the source code of this file.

Functions
static std::unique_ptr< PulseShapeFitOOTPileupCorrection >	parseHBHEMethod2Description (const edm::ParameterSet &conf)
static std::unique_ptr< HcalDeterministicFit >	parseHBHEMethod3Description (const edm::ParameterSet &conf)
std::unique_ptr< AbsHBHEPhase1Algo >	parseHBHEPhase1AlgoDescription (const edm::ParameterSet &ps)

Function Documentation

static std::unique_ptr<PulseShapeFitOOTPileupCorrection> parseHBHEMethod2Description ( const edm::ParameterSet &  conf )

static

Definition at line 14 of file parseHBHEPhase1AlgoDescription.cc.

References corr, edm::ParameterSet::getParameter(), and HcalTimeSlew::Medium.

Referenced by parseHBHEPhase1AlgoDescription().

{
    const bool iPedestalConstraint = conf.getParameter<bool>  ("applyPedConstraint");
    const bool iDCConstraint = conf.getParameter<bool>  ("applyDCConstraint");
    const bool iTimeConstraint =     conf.getParameter<bool>  ("applyTimeConstraint");
    const bool iAddPulseJitter =     conf.getParameter<bool>  ("applyPulseJitter");
    const bool iApplyTimeSlew =      conf.getParameter<bool>  ("applyTimeSlew");
    const double iTS4Min =           conf.getParameter<double>("ts4Min");
    const std::vector<double> iTS4Max =           conf.getParameter<std::vector<double>>("ts4Max");
    const double iPulseJitter =      conf.getParameter<double>("pulseJitter");
    const double iTimeMean =         conf.getParameter<double>("meanTime");
    const double iTimeSigHPD =       conf.getParameter<double>("timeSigmaHPD");
    const double iTimeSigSiPM =      conf.getParameter<double>("timeSigmaSiPM");
    const double iPedMean =          conf.getParameter<double>("meanPed");
    const double iPedSigHPD =        conf.getParameter<double>("pedSigmaHPD");
    const double iPedSigSiPM =       conf.getParameter<double>("pedSigmaSiPM");
    const double iNoiseHPD =         conf.getParameter<double>("noiseHPD");
    const double iNoiseSiPM =        conf.getParameter<double>("noiseSiPM");
    const double iTMin =             conf.getParameter<double>("timeMin");
    const double iTMax =             conf.getParameter<double>("timeMax");
    const std::vector<double> its4Chi2 =           conf.getParameter<std::vector<double>>("ts4chi2");
    const int iFitTimes =            conf.getParameter<int>   ("fitTimes");

    if (iPedestalConstraint) assert(iPedSigHPD);
    if (iPedestalConstraint) assert(iPedSigSiPM);
    if (iTimeConstraint) assert(iTimeSigHPD);
    if (iTimeConstraint) assert(iTimeSigSiPM);

    std::unique_ptr<PulseShapeFitOOTPileupCorrection> corr =
        std::make_unique<PulseShapeFitOOTPileupCorrection>();

    corr->setPUParams(iPedestalConstraint, iTimeConstraint, iAddPulseJitter,
                      iApplyTimeSlew, iTS4Min, iTS4Max,
                      iPulseJitter,
              iTimeMean, iTimeSigHPD, iTimeSigSiPM,
              iPedMean, iPedSigHPD, iPedSigSiPM,
                      iNoiseHPD, iNoiseSiPM,
              iTMin, iTMax, its4Chi2,
                      HcalTimeSlew::Medium, iFitTimes, iDCConstraint);

    return corr;
}

static std::unique_ptr<HcalDeterministicFit> parseHBHEMethod3Description ( const edm::ParameterSet &  conf )

static

Definition at line 59 of file parseHBHEPhase1AlgoDescription.cc.

References trackingPlots::fit, edm::ParameterSet::getParameter(), PedestalSub::init(), and HcalTimeSlew::Medium.

Referenced by parseHBHEPhase1AlgoDescription().

{
    const bool iApplyTimeSlew  =  conf.getParameter<bool>  ("applyTimeSlewM3");
    const float iPedSubThreshold =  conf.getParameter<double>("pedestalUpperLimit");
    const int iTimeSlewParsType  =  conf.getParameter<int>   ("timeSlewParsType");
    const double irespCorrM3 =     conf.getParameter<double>("respCorrM3");
    const std::vector<double>& iTimeSlewPars =
                     conf.getParameter<std::vector<double> >("timeSlewPars");

    PedestalSub pedSubFxn;
    pedSubFxn.init(0, iPedSubThreshold, 0.0);

    std::unique_ptr<HcalDeterministicFit> fit = std::make_unique<HcalDeterministicFit>();
    fit->init( (HcalTimeSlew::ParaSource)iTimeSlewParsType,
           HcalTimeSlew::Medium, iApplyTimeSlew,
           pedSubFxn, iTimeSlewPars, irespCorrM3);
    return fit;
}

std::unique_ptr<AbsHBHEPhase1Algo> parseHBHEPhase1AlgoDescription

(

const edm::ParameterSet &

ps

)

Definition at line 80 of file parseHBHEPhase1AlgoDescription.cc.

References patPFMETCorrections_cff::algo, className(), edm::ParameterSet::getParameter(), eostools::move(), parseHBHEMethod2Description(), parseHBHEMethod3Description(), and AlCaHLTBitMon_QueryRunRegistry::string.

{
    std::unique_ptr<AbsHBHEPhase1Algo> algo;

    const std::string& className = ps.getParameter<std::string>("Class");

    if (className == "SimpleHBHEPhase1Algo")
    {
        std::unique_ptr<PulseShapeFitOOTPileupCorrection> m2;
        if (ps.getParameter<bool>("useM2"))
            m2 = parseHBHEMethod2Description(ps);

        std::unique_ptr<HcalDeterministicFit> detFit;
        if (ps.getParameter<bool>("useM3"))
            detFit = parseHBHEMethod3Description(ps);

        algo = std::unique_ptr<AbsHBHEPhase1Algo>(
            new SimpleHBHEPhase1Algo(ps.getParameter<int>   ("firstSampleShift"),
                                     ps.getParameter<int>   ("samplesToAdd"),
                                     ps.getParameter<double>("correctionPhaseNS"),
                                     ps.getParameter<double>("tdcTimeShift"),
                                     ps.getParameter<bool>  ("correctForPhaseContainment"),
                                     std::move(m2), std::move(detFit))
            );
    }

    return algo;
}