parseHFPhase1AlgoDescription.cc

Go to the documentation of this file.

#include <cfloat>

#include "RecoLocalCalo/HcalRecAlgos/interface/parseHFPhase1AlgoDescription.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

// Phase 1 HF reco algorithm headers
#include "RecoLocalCalo/HcalRecAlgos/interface/HFSimpleTimeCheck.h"
#include "RecoLocalCalo/HcalRecAlgos/interface/HFFlexibleTimeCheck.h"

std::unique_ptr<AbsHFPhase1Algo>
parseHFPhase1AlgoDescription(const edm::ParameterSet& ps)
{
    std::unique_ptr<AbsHFPhase1Algo> algo;

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

    const bool isHFSimpleTimeCheck = className == "HFSimpleTimeCheck";
    if (isHFSimpleTimeCheck || className == "HFFlexibleTimeCheck")
    {
        const std::vector<double>& energyWeightsVec =
            ps.getParameter<std::vector<double> >("energyWeights");
        const unsigned soiPhase =
            ps.getParameter<unsigned>("soiPhase");
        const float timeShift =
            ps.getParameter<double>("timeShift");
        const float triseIfNoTDC =
            ps.getParameter<double>("triseIfNoTDC");
        const float tfallIfNoTDC =
            ps.getParameter<double>("tfallIfNoTDC");
        const bool rejectAllFailures =
            ps.getParameter<bool>("rejectAllFailures");
        const float minChargeForUndershoot =
            ps.getParameter<double>("minChargeForUndershoot");
        const float minChargeForOvershoot =
            ps.getParameter<double>("minChargeForOvershoot");
        const bool alwaysCalculateQAsymmetry =
            ps.getParameter<bool>("alwaysCalculateQAsymmetry");

        float energyWeights[2*HFAnodeStatus::N_POSSIBLE_STATES-1][2];
        const unsigned sz = sizeof(energyWeights)/sizeof(energyWeights[0][0]);

        if (energyWeightsVec.size() == sz)
        {
            std::pair<float,float> tlimits[2];
            if (isHFSimpleTimeCheck)
            {
                // Must specify the time limits explicitly for this algorithm
                const std::vector<double>& tlimitsVec =
                    ps.getParameter<std::vector<double> >("tlimits");
                if (tlimitsVec.size() == 4)
                {
                    tlimits[0] = std::pair<float,float>(tlimitsVec[0], tlimitsVec[1]);
                    tlimits[1] = std::pair<float,float>(tlimitsVec[2], tlimitsVec[3]);
                }
                else
                    return algo;
            }
            else
            {
                // Use "all pass" time limits values, just in case
                tlimits[0] = std::pair<float,float>(-FLT_MAX, FLT_MAX);
                tlimits[1] = tlimits[0];
            }

            // Same order of elements as in the natural C array mapping
            float* to = &energyWeights[0][0];
            for (unsigned i=0; i<sz; ++i)
                to[i] = energyWeightsVec[i];

            // Create the algorithm object
            if (isHFSimpleTimeCheck)
                algo = std::unique_ptr<AbsHFPhase1Algo>(
                    new HFSimpleTimeCheck(tlimits, energyWeights, soiPhase,
                                          timeShift, triseIfNoTDC, tfallIfNoTDC,
                                          minChargeForUndershoot, minChargeForOvershoot,
                                          rejectAllFailures, alwaysCalculateQAsymmetry));
            else
                algo = std::unique_ptr<AbsHFPhase1Algo>(
                    new HFFlexibleTimeCheck(tlimits, energyWeights, soiPhase,
                                            timeShift, triseIfNoTDC, tfallIfNoTDC,
                                            minChargeForUndershoot, minChargeForOvershoot,
                                            rejectAllFailures, alwaysCalculateQAsymmetry));
        }
    }

    return algo;
}

edm::ParameterSetDescription fillDescriptionForParseHFPhase1AlgoDescription()
{
    edm::ParameterSetDescription desc;

    std::vector<double> allPass{-10000.0, 10000.0, -10000.0, 10000.0};
    desc.add<std::vector<double> >("tlimits", allPass);
    desc.add<std::vector<double> >("energyWeights");
    desc.add<unsigned>("soiPhase", 1U);
    desc.add<double>("timeShift", 0.0);
    desc.add<double>("triseIfNoTDC", -100.0);
    desc.add<double>("tfallIfNoTDC", -101.0);
    desc.add<double>("minChargeForUndershoot", 1.0e10);
    desc.add<double>("minChargeForOvershoot", 1.0e10);
    desc.add<bool>("alwaysCalculateQAsymmetry", true);

    desc.ifValue(edm::ParameterDescription<std::string>("Class", "HFSimpleTimeCheck", true),
                 "HFSimpleTimeCheck" >> edm::ParameterDescription<bool>("rejectAllFailures", false, true) or
                 "HFFlexibleTimeCheck" >> edm::ParameterDescription<bool>("rejectAllFailures", true, true));

    return desc;
}