L1MuTriggerScalesOnlineProducer.cc

Go to the documentation of this file.

//-------------------------------------------------
//
//   \class L1MuTriggerScalesOnlineProducer
//
//   Description:  A class to produce the L1 mu emulator scales record in the event setup
//                 from the OMDS database.
//
//
//   Author :
//   Thomas Themel
//
//--------------------------------------------------
#include "L1TriggerConfig/L1ScalesProducers/interface/L1MuTriggerScalesOnlineProducer.h"
#include "L1TriggerConfig/L1ScalesProducers/interface/ScaleRecordHelper.h"
#include <cmath>

using namespace std;

L1MuTriggerScalesOnlineProducer::L1MuTriggerScalesOnlineProducer(const edm::ParameterSet& ps)
    : L1ConfigOnlineProdBase<L1MuTriggerScalesRcd, L1MuTriggerScales>(ps),
      m_scales(ps.getParameter<int>("nbitPackingDTEta"),
               ps.getParameter<bool>("signedPackingDTEta"),
               ps.getParameter<int>("nbinsDTEta"),
               ps.getParameter<double>("minDTEta"),
               ps.getParameter<double>("maxDTEta"),
               ps.getParameter<int>("offsetDTEta"),

               ps.getParameter<int>("nbitPackingCSCEta"),
               ps.getParameter<int>("nbinsCSCEta"),
               ps.getParameter<double>("minCSCEta"),
               ps.getParameter<double>("maxCSCEta"),

               ps.getParameter<std::vector<double> >("scaleRPCEta"),
               ps.getParameter<int>("nbitPackingBrlRPCEta"),
               ps.getParameter<bool>("signedPackingBrlRPCEta"),
               ps.getParameter<int>("nbinsBrlRPCEta"),
               ps.getParameter<int>("offsetBrlRPCEta"),
               ps.getParameter<int>("nbitPackingFwdRPCEta"),
               ps.getParameter<bool>("signedPackingFwdRPCEta"),
               ps.getParameter<int>("nbinsFwdRPCEta"),
               ps.getParameter<int>("offsetFwdRPCEta"),
               // Fields that should now be generated from OMDS:
               // TODO: Adjust m_scales's definition to be a bit
               //       more accessible for the partial initialization.
               //ps.getParameter<int>("nbitPackingGMTEta"),
               0,
               //ps.getParameter<int>("nbinsGMTEta"),
               0,
               //ps.getParameter<std::vector<double> >("scaleGMTEta"),
               std::vector<double>(1),
               //ps.getParameter<int>("nbitPackingPhi"),
               0,
               //ps.getParameter<bool>("signedPackingPhi"),
               false,
               //ps.getParameter<int>("nbinsPhi"),
               0,
               //ps.getParameter<double>("minPhi"),
               0,
               //ps.getParameter<double>("maxPhi")
               0),
      /* Metadata that's not yet in the database. */
      m_nbitPackingPhi(ps.getParameter<int>("nbitPackingPhi")),
      m_nbitPackingEta(ps.getParameter<int>("nbitPackingGMTEta")),
      m_nbinsEta(ps.getParameter<int>("nbinsGMTEta")),
      m_signedPackingPhi(ps.getParameter<bool>("signedPackingPhi")) {}

L1MuTriggerScalesOnlineProducer::~L1MuTriggerScalesOnlineProducer() {}

//
// member functions
//

class PhiScaleHelper {
public:
  static L1MuBinnedScale* makeBinnedScale(l1t::OMDSReader::QueryResults& record, int nBits, bool signedPacking) {
    short nbins = 0;
    record.fillVariable(BinsColumn, nbins);
    float lowMark = 0.;
    record.fillVariable(LowMarkColumn, lowMark);
    float step = 0.;
    record.fillVariable(StepColumn, step);

    return new L1MuBinnedScale(nBits, signedPacking, nbins, deg2rad(lowMark), deg2rad(lowMark + nbins * step));
  }

  static void pushColumnNames(vector<string>& columns) {
    columns.push_back(BinsColumn);
    columns.push_back(LowMarkColumn);
    columns.push_back(StepColumn);
  }

  static double deg2rad(double deg) { return deg * M_PI / 180.0; }
  static double rad2deg(double rad) { return rad / M_PI * 180.0; }

  static const string BinsColumn;
  static const string LowMarkColumn;
  static const string StepColumn;
};

const string PhiScaleHelper::BinsColumn = "PHI_BINS";
const string PhiScaleHelper::LowMarkColumn = "PHI_DEG_BIN_LOW_0";
const string PhiScaleHelper::StepColumn = "PHI_DEG_BIN_STEP";

// ------------ method called to produce the data  ------------
std::unique_ptr<L1MuTriggerScales> L1MuTriggerScalesOnlineProducer::newObject(const std::string& objectKey) {
  // The key we get from the O2O subsystem is the CMS_GMT.L1T_SCALES key,
  // but the eta/phi scales have their own subtables, so let's find
  // out.
  vector<string> foreignKeys;

  const std::string etaKeyColumn("SC_MUON_ETA_FK");
  const std::string phiKeyColumn("SC_MUON_PHI_FK");

  foreignKeys.push_back(etaKeyColumn);
  foreignKeys.push_back(phiKeyColumn);

  l1t::OMDSReader::QueryResults keysRecord = m_omdsReader.basicQuery(
      // SELECTed columns
      foreignKeys,
      // schema name
      "CMS_GT",
      // table name
      "L1T_SCALES",
      // WHERE lhs
      "L1T_SCALES.ID",
      // WHERE rhs
      m_omdsReader.singleAttribute(objectKey));

  if (keysRecord.numberRows() != 1)  // check if query was successful
  {
    throw cond::Exception(
        "Problem finding L1MuTriggerScales associated "
        "with scales key `" +
        objectKey + "'");
  }

  std::string etaKeyValue;
  std::string phiKeyValue;
  keysRecord.fillVariable(etaKeyColumn, etaKeyValue);
  keysRecord.fillVariable(phiKeyColumn, phiKeyValue);

  vector<string> columns;

  // get the eta scales from the database
  ScaleRecordHelper etaHelper("ETA_BIN_LOW", m_nbinsEta);
  etaHelper.pushColumnNames(columns);

  l1t::OMDSReader::QueryResults etaRecord = m_omdsReader.basicQuery(
      // SELECTed columns
      columns,
      // schema name
      "CMS_GT",
      // table name
      "L1T_SCALE_MUON_ETA",
      // WHERE lhs
      "L1T_SCALE_MUON_ETA.ID",
      // WHERE rhs
      m_omdsReader.singleAttribute(etaKeyValue));

  vector<double> etaScales;
  etaHelper.extractScales(etaRecord, etaScales);

  unique_ptr<L1MuSymmetricBinnedScale> ptrEtaScale(
      new L1MuSymmetricBinnedScale(m_nbitPackingEta, m_nbinsEta, etaScales));
  m_scales.setGMTEtaScale(*ptrEtaScale);

  columns.clear();

  // get the phi scales from the database
  PhiScaleHelper phiHelper;

  l1t::OMDSReader::QueryResults phiRecord = m_omdsReader.basicQuery(
      // SELECTed columns
      columns,
      // schema name
      "CMS_GT",
      // table name
      "L1T_SCALE_MUON_PHI",
      // WHERE lhs
      "L1T_SCALE_MUON_PHI.ID",
      // WHERE rhs
      m_omdsReader.singleAttribute(phiKeyValue));

  unique_ptr<L1MuBinnedScale> ptrPhiScale(phiHelper.makeBinnedScale(phiRecord, m_nbitPackingPhi, m_signedPackingPhi));

  m_scales.setPhiScale(*ptrPhiScale);

  return std::make_unique<L1MuTriggerScales>(m_scales);
}