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::shared_ptr<L1MuTriggerScales> L1MuTriggerScalesOnlineProducer::newObject(const std::string& objectKey ) 
{
   using namespace edm::es;   

   // 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_shared<L1MuTriggerScales>(m_scales);
}