L1TMuonGlobalParamsESProducer.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    L1Trigger/L1TMuonGlobalParamsESProducer
// Class:      L1TMuonGlobalParamsESProducer
//
//
// Original Author:  Thomas Reis
//         Created:  Mon, 21 Sep 2015 13:28:49 GMT
//
//

// system include files
#include <memory>
#include <strstream>

// user include files
#include "FWCore/Framework/interface/ModuleFactory.h"
#include "FWCore/Framework/interface/ESProducer.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/ESProducts.h"

#include "CondFormats/L1TObjects/interface/L1TMuonGlobalParams.h"
#include "CondFormats/DataRecord/interface/L1TMuonGlobalParamsRcd.h"
#include "L1Trigger/L1TMuon/interface/L1TMuonGlobalParamsHelper.h"
#include "L1Trigger/L1TMuon/interface/MicroGMTLUTFactories.h"
#include "L1Trigger/L1TCommon/interface/TriggerSystem.h"
#include "L1Trigger/L1TCommon/interface/Parameter.h"
#include "L1Trigger/L1TCommon/interface/Mask.h"

//
// class declaration
//

class L1TMuonGlobalParamsESProducer : public edm::ESProducer {
public:
  L1TMuonGlobalParamsESProducer(const edm::ParameterSet&);
  ~L1TMuonGlobalParamsESProducer() override;

  using ReturnType = std::unique_ptr<L1TMuonGlobalParams>;

  ReturnType produce(const L1TMuonGlobalParamsRcd&);

private:
  L1TMuonGlobalParams m_params;
};

//
// constants, enums and typedefs
//

//
// static data member definitions
//

//
// constructors and destructor
//
L1TMuonGlobalParamsESProducer::L1TMuonGlobalParamsESProducer(const edm::ParameterSet& iConfig) {
  //the following line is needed to tell the framework what
  // data is being produced
  setWhatProduced(this);

  L1TMuonGlobalParamsHelper m_params_helper;

  // get configuration from DB
  if (iConfig.getParameter<bool>("configFromXml")) {
    l1t::TriggerSystem trgSys;
    edm::FileInPath hwXmlFile(iConfig.getParameter<std::string>("hwXmlFile"));
    edm::FileInPath topCfgXmlFile(iConfig.getParameter<std::string>("topCfgXmlFile"));
    // These xml files are for testing the configuration from the online DB
    trgSys.configureSystemFromFiles(hwXmlFile.fullPath().c_str(),
                                    topCfgXmlFile.fullPath().c_str(),
                                    iConfig.getParameter<std::string>("xmlCfgKey").c_str());

    m_params_helper.loadFromOnline(trgSys, iConfig.getParameter<std::string>("uGmtProcessorId"));
  } else {
    // Firmware version
    unsigned fwVersion = iConfig.getParameter<unsigned>("fwVersion");
    m_params_helper.setFwVersion(fwVersion);

    // LUTs
    m_params_helper.setFwdPosSingleMatchQualLUTMaxDR(iConfig.getParameter<double>("FwdPosSingleMatchQualLUTMaxDR"),
                                                     iConfig.getParameter<double>("FwdPosSingleMatchQualLUTfEta"),
                                                     iConfig.getParameter<double>("FwdPosSingleMatchQualLUTfPhi"));
    m_params_helper.setFwdNegSingleMatchQualLUTMaxDR(iConfig.getParameter<double>("FwdNegSingleMatchQualLUTMaxDR"),
                                                     iConfig.getParameter<double>("FwdNegSingleMatchQualLUTfEta"),
                                                     iConfig.getParameter<double>("FwdNegSingleMatchQualLUTfPhi"));
    m_params_helper.setOvlPosSingleMatchQualLUTMaxDR(iConfig.getParameter<double>("OvlPosSingleMatchQualLUTMaxDR"),
                                                     iConfig.getParameter<double>("OvlPosSingleMatchQualLUTfEta"),
                                                     iConfig.getParameter<double>("OvlPosSingleMatchQualLUTfEtaCoarse"),
                                                     iConfig.getParameter<double>("OvlPosSingleMatchQualLUTfPhi"));
    m_params_helper.setOvlNegSingleMatchQualLUTMaxDR(iConfig.getParameter<double>("OvlNegSingleMatchQualLUTMaxDR"),
                                                     iConfig.getParameter<double>("OvlNegSingleMatchQualLUTfEta"),
                                                     iConfig.getParameter<double>("OvlNegSingleMatchQualLUTfEtaCoarse"),
                                                     iConfig.getParameter<double>("OvlNegSingleMatchQualLUTfPhi"));
    m_params_helper.setBOPosMatchQualLUTMaxDR(iConfig.getParameter<double>("BOPosMatchQualLUTMaxDR"),
                                              iConfig.getParameter<double>("BOPosMatchQualLUTfEta"),
                                              iConfig.getParameter<double>("BOPosMatchQualLUTfEtaCoarse"),
                                              iConfig.getParameter<double>("BOPosMatchQualLUTfPhi"));
    m_params_helper.setBONegMatchQualLUTMaxDR(iConfig.getParameter<double>("BONegMatchQualLUTMaxDR"),
                                              iConfig.getParameter<double>("BONegMatchQualLUTfEta"),
                                              iConfig.getParameter<double>("BONegMatchQualLUTfEtaCoarse"),
                                              iConfig.getParameter<double>("BONegMatchQualLUTfPhi"));
    m_params_helper.setFOPosMatchQualLUTMaxDR(iConfig.getParameter<double>("FOPosMatchQualLUTMaxDR"),
                                              iConfig.getParameter<double>("FOPosMatchQualLUTfEta"),
                                              iConfig.getParameter<double>("FOPosMatchQualLUTfEtaCoarse"),
                                              iConfig.getParameter<double>("FOPosMatchQualLUTfPhi"));
    m_params_helper.setFONegMatchQualLUTMaxDR(iConfig.getParameter<double>("FONegMatchQualLUTMaxDR"),
                                              iConfig.getParameter<double>("FONegMatchQualLUTfEta"),
                                              iConfig.getParameter<double>("FONegMatchQualLUTfEtaCoarse"),
                                              iConfig.getParameter<double>("FONegMatchQualLUTfPhi"));

    unsigned sortRankLUTPtFactor = iConfig.getParameter<unsigned>("SortRankLUTPtFactor");
    unsigned sortRankLUTQualFactor = iConfig.getParameter<unsigned>("SortRankLUTQualFactor");
    m_params_helper.setSortRankLUTFactors(sortRankLUTPtFactor, sortRankLUTQualFactor);

    auto absIsoCheckMemLUT = l1t::MicroGMTAbsoluteIsolationCheckLUTFactory::create(
        iConfig.getParameter<std::string>("AbsIsoCheckMemLUTPath"), fwVersion);
    auto relIsoCheckMemLUT = l1t::MicroGMTRelativeIsolationCheckLUTFactory::create(
        iConfig.getParameter<std::string>("RelIsoCheckMemLUTPath"), fwVersion);
    auto idxSelMemPhiLUT = l1t::MicroGMTCaloIndexSelectionLUTFactory::create(
        iConfig.getParameter<std::string>("IdxSelMemPhiLUTPath"), l1t::MicroGMTConfiguration::PHI, fwVersion);
    auto idxSelMemEtaLUT = l1t::MicroGMTCaloIndexSelectionLUTFactory::create(
        iConfig.getParameter<std::string>("IdxSelMemEtaLUTPath"), l1t::MicroGMTConfiguration::ETA, fwVersion);
    auto fwdPosSingleMatchQualLUT = l1t::MicroGMTMatchQualLUTFactory::create(
        iConfig.getParameter<std::string>("FwdPosSingleMatchQualLUTPath"),
        iConfig.getParameter<double>("FwdPosSingleMatchQualLUTMaxDR"),
        iConfig.getParameter<double>("FwdPosSingleMatchQualLUTfEta"),
        iConfig.getParameter<double>("FwdPosSingleMatchQualLUTfEta"),  // set the coarse eta factor = fine eta factor
        iConfig.getParameter<double>("FwdPosSingleMatchQualLUTfPhi"),
        l1t::cancel_t::emtf_emtf_pos,
        fwVersion);
    auto fwdNegSingleMatchQualLUT = l1t::MicroGMTMatchQualLUTFactory::create(
        iConfig.getParameter<std::string>("FwdNegSingleMatchQualLUTPath"),
        iConfig.getParameter<double>("FwdNegSingleMatchQualLUTMaxDR"),
        iConfig.getParameter<double>("FwdNegSingleMatchQualLUTfEta"),
        iConfig.getParameter<double>("FwdNegSingleMatchQualLUTfEta"),  // set the coarse eta factor = fine eta factor
        iConfig.getParameter<double>("FwdNegSingleMatchQualLUTfPhi"),
        l1t::cancel_t::emtf_emtf_neg,
        fwVersion);
    auto ovlPosSingleMatchQualLUT =
        l1t::MicroGMTMatchQualLUTFactory::create(iConfig.getParameter<std::string>("OvlPosSingleMatchQualLUTPath"),
                                                 iConfig.getParameter<double>("OvlPosSingleMatchQualLUTMaxDR"),
                                                 iConfig.getParameter<double>("OvlPosSingleMatchQualLUTfEta"),
                                                 iConfig.getParameter<double>("OvlPosSingleMatchQualLUTfEtaCoarse"),
                                                 iConfig.getParameter<double>("OvlPosSingleMatchQualLUTfPhi"),
                                                 l1t::cancel_t::omtf_omtf_pos,
                                                 fwVersion);
    auto ovlNegSingleMatchQualLUT =
        l1t::MicroGMTMatchQualLUTFactory::create(iConfig.getParameter<std::string>("OvlNegSingleMatchQualLUTPath"),
                                                 iConfig.getParameter<double>("OvlNegSingleMatchQualLUTMaxDR"),
                                                 iConfig.getParameter<double>("OvlNegSingleMatchQualLUTfEta"),
                                                 iConfig.getParameter<double>("OvlNegSingleMatchQualLUTfEtaCoarse"),
                                                 iConfig.getParameter<double>("OvlNegSingleMatchQualLUTfPhi"),
                                                 l1t::cancel_t::omtf_omtf_neg,
                                                 fwVersion);
    auto bOPosMatchQualLUT =
        l1t::MicroGMTMatchQualLUTFactory::create(iConfig.getParameter<std::string>("BOPosMatchQualLUTPath"),
                                                 iConfig.getParameter<double>("BOPosMatchQualLUTMaxDR"),
                                                 iConfig.getParameter<double>("BOPosMatchQualLUTfEta"),
                                                 iConfig.getParameter<double>("BOPosMatchQualLUTfEtaCoarse"),
                                                 iConfig.getParameter<double>("BOPosMatchQualLUTfPhi"),
                                                 l1t::cancel_t::omtf_bmtf_pos,
                                                 fwVersion);
    auto bONegMatchQualLUT =
        l1t::MicroGMTMatchQualLUTFactory::create(iConfig.getParameter<std::string>("BONegMatchQualLUTPath"),
                                                 iConfig.getParameter<double>("BONegMatchQualLUTMaxDR"),
                                                 iConfig.getParameter<double>("BONegMatchQualLUTfEta"),
                                                 iConfig.getParameter<double>("BONegMatchQualLUTfEtaCoarse"),
                                                 iConfig.getParameter<double>("BONegMatchQualLUTfPhi"),
                                                 l1t::cancel_t::omtf_bmtf_neg,
                                                 fwVersion);
    auto fOPosMatchQualLUT =
        l1t::MicroGMTMatchQualLUTFactory::create(iConfig.getParameter<std::string>("FOPosMatchQualLUTPath"),
                                                 iConfig.getParameter<double>("FOPosMatchQualLUTMaxDR"),
                                                 iConfig.getParameter<double>("FOPosMatchQualLUTfEta"),
                                                 iConfig.getParameter<double>("FOPosMatchQualLUTfEtaCoarse"),
                                                 iConfig.getParameter<double>("FOPosMatchQualLUTfPhi"),
                                                 l1t::cancel_t::omtf_emtf_pos,
                                                 fwVersion);
    auto fONegMatchQualLUT =
        l1t::MicroGMTMatchQualLUTFactory::create(iConfig.getParameter<std::string>("FONegMatchQualLUTPath"),
                                                 iConfig.getParameter<double>("FONegMatchQualLUTMaxDR"),
                                                 iConfig.getParameter<double>("FONegMatchQualLUTfEta"),
                                                 iConfig.getParameter<double>("FONegMatchQualLUTfEtaCoarse"),
                                                 iConfig.getParameter<double>("FONegMatchQualLUTfPhi"),
                                                 l1t::cancel_t::omtf_emtf_neg,
                                                 fwVersion);
    auto bPhiExtrapolationLUT = l1t::MicroGMTExtrapolationLUTFactory::create(
        iConfig.getParameter<std::string>("BPhiExtrapolationLUTPath"), l1t::MicroGMTConfiguration::PHI_OUT, fwVersion);
    auto oPhiExtrapolationLUT = l1t::MicroGMTExtrapolationLUTFactory::create(
        iConfig.getParameter<std::string>("OPhiExtrapolationLUTPath"), l1t::MicroGMTConfiguration::PHI_OUT, fwVersion);
    auto fPhiExtrapolationLUT = l1t::MicroGMTExtrapolationLUTFactory::create(
        iConfig.getParameter<std::string>("FPhiExtrapolationLUTPath"), l1t::MicroGMTConfiguration::PHI_OUT, fwVersion);
    auto bEtaExtrapolationLUT = l1t::MicroGMTExtrapolationLUTFactory::create(
        iConfig.getParameter<std::string>("BEtaExtrapolationLUTPath"), l1t::MicroGMTConfiguration::ETA_OUT, fwVersion);
    auto oEtaExtrapolationLUT = l1t::MicroGMTExtrapolationLUTFactory::create(
        iConfig.getParameter<std::string>("OEtaExtrapolationLUTPath"), l1t::MicroGMTConfiguration::ETA_OUT, fwVersion);
    auto fEtaExtrapolationLUT = l1t::MicroGMTExtrapolationLUTFactory::create(
        iConfig.getParameter<std::string>("FEtaExtrapolationLUTPath"), l1t::MicroGMTConfiguration::ETA_OUT, fwVersion);
    auto rankPtQualityLUT =
        l1t::MicroGMTRankPtQualLUTFactory::create(iConfig.getParameter<std::string>("SortRankLUTPath"),
                                                  fwVersion,
                                                  sortRankLUTPtFactor,
                                                  sortRankLUTQualFactor);  // LUTs defined from config file
    m_params_helper.setAbsIsoCheckMemLUT(*absIsoCheckMemLUT);
    m_params_helper.setRelIsoCheckMemLUT(*relIsoCheckMemLUT);
    m_params_helper.setIdxSelMemPhiLUT(*idxSelMemPhiLUT);
    m_params_helper.setIdxSelMemEtaLUT(*idxSelMemEtaLUT);
    m_params_helper.setFwdPosSingleMatchQualLUT(*fwdPosSingleMatchQualLUT);
    m_params_helper.setFwdNegSingleMatchQualLUT(*fwdNegSingleMatchQualLUT);
    m_params_helper.setOvlPosSingleMatchQualLUT(*ovlPosSingleMatchQualLUT);
    m_params_helper.setOvlNegSingleMatchQualLUT(*ovlNegSingleMatchQualLUT);
    m_params_helper.setBOPosMatchQualLUT(*bOPosMatchQualLUT);
    m_params_helper.setBONegMatchQualLUT(*bONegMatchQualLUT);
    m_params_helper.setFOPosMatchQualLUT(*fOPosMatchQualLUT);
    m_params_helper.setFONegMatchQualLUT(*fONegMatchQualLUT);
    m_params_helper.setBPhiExtrapolationLUT(*bPhiExtrapolationLUT);
    m_params_helper.setOPhiExtrapolationLUT(*oPhiExtrapolationLUT);
    m_params_helper.setFPhiExtrapolationLUT(*fPhiExtrapolationLUT);
    m_params_helper.setBEtaExtrapolationLUT(*bEtaExtrapolationLUT);
    m_params_helper.setOEtaExtrapolationLUT(*oEtaExtrapolationLUT);
    m_params_helper.setFEtaExtrapolationLUT(*fEtaExtrapolationLUT);
    m_params_helper.setSortRankLUT(*rankPtQualityLUT);

    // LUT paths
    m_params_helper.setAbsIsoCheckMemLUTPath(iConfig.getParameter<std::string>("AbsIsoCheckMemLUTPath"));
    m_params_helper.setRelIsoCheckMemLUTPath(iConfig.getParameter<std::string>("RelIsoCheckMemLUTPath"));
    m_params_helper.setIdxSelMemPhiLUTPath(iConfig.getParameter<std::string>("IdxSelMemPhiLUTPath"));
    m_params_helper.setIdxSelMemEtaLUTPath(iConfig.getParameter<std::string>("IdxSelMemEtaLUTPath"));
    m_params_helper.setFwdPosSingleMatchQualLUTPath(iConfig.getParameter<std::string>("FwdPosSingleMatchQualLUTPath"));
    m_params_helper.setFwdNegSingleMatchQualLUTPath(iConfig.getParameter<std::string>("FwdNegSingleMatchQualLUTPath"));
    m_params_helper.setOvlPosSingleMatchQualLUTPath(iConfig.getParameter<std::string>("OvlPosSingleMatchQualLUTPath"));
    m_params_helper.setOvlNegSingleMatchQualLUTPath(iConfig.getParameter<std::string>("OvlNegSingleMatchQualLUTPath"));
    m_params_helper.setBOPosMatchQualLUTPath(iConfig.getParameter<std::string>("BOPosMatchQualLUTPath"));
    m_params_helper.setBONegMatchQualLUTPath(iConfig.getParameter<std::string>("BONegMatchQualLUTPath"));
    m_params_helper.setFOPosMatchQualLUTPath(iConfig.getParameter<std::string>("FOPosMatchQualLUTPath"));
    m_params_helper.setFONegMatchQualLUTPath(iConfig.getParameter<std::string>("FONegMatchQualLUTPath"));
    m_params_helper.setBPhiExtrapolationLUTPath(iConfig.getParameter<std::string>("BPhiExtrapolationLUTPath"));
    m_params_helper.setOPhiExtrapolationLUTPath(iConfig.getParameter<std::string>("OPhiExtrapolationLUTPath"));
    m_params_helper.setFPhiExtrapolationLUTPath(iConfig.getParameter<std::string>("FPhiExtrapolationLUTPath"));
    m_params_helper.setBEtaExtrapolationLUTPath(iConfig.getParameter<std::string>("BEtaExtrapolationLUTPath"));
    m_params_helper.setOEtaExtrapolationLUTPath(iConfig.getParameter<std::string>("OEtaExtrapolationLUTPath"));
    m_params_helper.setFEtaExtrapolationLUTPath(iConfig.getParameter<std::string>("FEtaExtrapolationLUTPath"));
    m_params_helper.setSortRankLUTPath(iConfig.getParameter<std::string>("SortRankLUTPath"));

    // uGMT disabled inputs
    bool disableCaloInputs = iConfig.getParameter<bool>("caloInputsDisable");
    std::vector<unsigned> bmtfInputsToDisable = iConfig.getParameter<std::vector<unsigned> >("bmtfInputsToDisable");
    std::vector<unsigned> omtfInputsToDisable = iConfig.getParameter<std::vector<unsigned> >("omtfInputsToDisable");
    std::vector<unsigned> emtfInputsToDisable = iConfig.getParameter<std::vector<unsigned> >("emtfInputsToDisable");

    if (disableCaloInputs) {
      m_params_helper.setCaloInputsToDisable(std::bitset<28>(0xFFFFFFF));
    } else {
      m_params_helper.setCaloInputsToDisable(std::bitset<28>());
    }

    std::bitset<12> bmtfDisables;
    for (size_t i = 0; i < bmtfInputsToDisable.size(); ++i) {
      bmtfDisables.set(i, bmtfInputsToDisable[i] > 0);
    }
    m_params_helper.setBmtfInputsToDisable(bmtfDisables);

    std::bitset<6> omtfpDisables;
    std::bitset<6> omtfnDisables;
    for (size_t i = 0; i < omtfInputsToDisable.size(); ++i) {
      if (i < 6) {
        omtfpDisables.set(i, omtfInputsToDisable[i] > 0);
      } else {
        omtfnDisables.set(i - 6, omtfInputsToDisable[i] > 0);
      }
    }
    m_params_helper.setOmtfpInputsToDisable(omtfpDisables);
    m_params_helper.setOmtfnInputsToDisable(omtfnDisables);

    std::bitset<6> emtfpDisables;
    std::bitset<6> emtfnDisables;
    for (size_t i = 0; i < emtfInputsToDisable.size(); ++i) {
      if (i < 6) {
        emtfpDisables.set(i, emtfInputsToDisable[i] > 0);
      } else {
        emtfnDisables.set(i - 6, emtfInputsToDisable[i] > 0);
      }
    }
    m_params_helper.setEmtfpInputsToDisable(emtfpDisables);
    m_params_helper.setEmtfnInputsToDisable(emtfnDisables);

    // masked inputs
    bool caloInputsMasked = iConfig.getParameter<bool>("caloInputsMasked");
    std::vector<unsigned> maskedBmtfInputs = iConfig.getParameter<std::vector<unsigned> >("maskedBmtfInputs");
    std::vector<unsigned> maskedOmtfInputs = iConfig.getParameter<std::vector<unsigned> >("maskedOmtfInputs");
    std::vector<unsigned> maskedEmtfInputs = iConfig.getParameter<std::vector<unsigned> >("maskedEmtfInputs");

    if (caloInputsMasked) {
      m_params_helper.setMaskedCaloInputs(std::bitset<28>(0xFFFFFFF));
    } else {
      m_params_helper.setMaskedCaloInputs(std::bitset<28>());
    }

    std::bitset<12> bmtfMasked;
    for (size_t i = 0; i < maskedBmtfInputs.size(); ++i) {
      bmtfMasked.set(i, maskedBmtfInputs[i] > 0);
    }
    m_params_helper.setMaskedBmtfInputs(bmtfMasked);

    std::bitset<6> omtfpMasked;
    std::bitset<6> omtfnMasked;
    for (size_t i = 0; i < maskedOmtfInputs.size(); ++i) {
      if (i < 6) {
        omtfpMasked.set(i, maskedOmtfInputs[i] > 0);
      } else {
        omtfnMasked.set(i - 6, maskedOmtfInputs[i] > 0);
      }
    }
    m_params_helper.setMaskedOmtfpInputs(omtfpMasked);
    m_params_helper.setMaskedOmtfnInputs(omtfnMasked);

    std::bitset<6> emtfpMasked;
    std::bitset<6> emtfnMasked;
    for (size_t i = 0; i < maskedEmtfInputs.size(); ++i) {
      if (i < 6) {
        emtfpMasked.set(i, maskedEmtfInputs[i] > 0);
      } else {
        emtfnMasked.set(i - 6, maskedEmtfInputs[i] > 0);
      }
    }
    m_params_helper.setMaskedEmtfpInputs(emtfpMasked);
    m_params_helper.setMaskedEmtfnInputs(emtfnMasked);
  }

  // temp hack to avoid ALCA/DB signoff:
  m_params = cast_to_L1TMuonGlobalParams((L1TMuonGlobalParams_PUBLIC)m_params_helper);
}

L1TMuonGlobalParamsESProducer::~L1TMuonGlobalParamsESProducer() {}

//
// member functions
//

// ------------ method called to produce the data  ------------
L1TMuonGlobalParamsESProducer::ReturnType L1TMuonGlobalParamsESProducer::produce(const L1TMuonGlobalParamsRcd& iRecord) {
  return std::make_unique<L1TMuonGlobalParams>(m_params);
}

//define this as a plug-in
DEFINE_FWK_EVENTSETUP_MODULE(L1TMuonGlobalParamsESProducer);