L1CaloHcalScaleConfigOnlineProd.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    L1EmEtScaleOnlineProd
// Class:      L1EmEtScaleOnlineProd
//
//
// Original Author:  Werner Man-Li Sun
//         Created:  Tue Sep 16 22:43:22 CEST 2008
//
//

// system include files

// user include files
#include "CondTools/L1Trigger/interface/L1ConfigOnlineProdBase.h"
#include "CondFormats/L1TObjects/interface/L1CaloHcalScale.h"
#include "CondFormats/DataRecord/interface/L1CaloHcalScaleRcd.h"
#include "Geometry/HcalTowerAlgo/interface/HcalTrigTowerGeometry.h"
#include "Geometry/Records/interface/CaloGeometryRecord.h"
#include "CalibCalorimetry/CaloTPG/src/CaloTPGTranscoderULUT.h"
#include "CondTools/L1Trigger/interface/OMDSReader.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include <cmath>
#include <iostream>
#include <iomanip>

//
// class declaration
//

class L1CaloHcalScaleConfigOnlineProd : public L1ConfigOnlineProdBase<L1CaloHcalScaleRcd, L1CaloHcalScale> {
public:
  L1CaloHcalScaleConfigOnlineProd(const edm::ParameterSet& iConfig);
  ~L1CaloHcalScaleConfigOnlineProd() override;

  std::unique_ptr<L1CaloHcalScale> produce(const L1CaloHcalScaleRcd& iRecord) override;

  std::unique_ptr<L1CaloHcalScale> newObject(const std::string& objectKey) override;

private:
  const HcalTrigTowerGeometry* theTrigTowerGeometry;
  CaloTPGTranscoderULUT* caloTPG;
  typedef std::vector<double> RCTdecompression;
  std::vector<RCTdecompression> hcaluncomp;

  //  HcaluLUTTPGCoder* tpgCoder;// = new    HcaluLUTTPGCoder();

  HcalTrigTowerDetId* ttDetId;

  // ----------member data ---------------------------
};

//
// constants, enums and typedefs
//

//
// static data member definitions
//

//
// constructors and destructor
//
L1CaloHcalScaleConfigOnlineProd::L1CaloHcalScaleConfigOnlineProd(const edm::ParameterSet& iConfig)
    : L1ConfigOnlineProdBase<L1CaloHcalScaleRcd, L1CaloHcalScale>(iConfig), theTrigTowerGeometry(nullptr) {
  caloTPG = new CaloTPGTranscoderULUT();
}

L1CaloHcalScaleConfigOnlineProd::~L1CaloHcalScaleConfigOnlineProd() {
  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)

  if (caloTPG != nullptr)
    delete caloTPG;
}

std::unique_ptr<L1CaloHcalScale> L1CaloHcalScaleConfigOnlineProd::newObject(const std::string& objectKey) {
  assert(theTrigTowerGeometry != nullptr);

  edm::LogInfo("L1CaloHcalScaleConfigOnlineProd") << "object Key " << objectKey;

  if (objectKey == "NULL" || objectKey.empty()) {  // return default blank ecal scale
    return std::make_unique<L1CaloHcalScale>(0);
  }
  if (objectKey == "IDENTITY") {  // return identity ecal scale
    return std::make_unique<L1CaloHcalScale>(1);
  }

  std::vector<unsigned int> analyticalLUT(1024, 0);
  std::vector<unsigned int> identityLUT(1024, 0);

  // Compute compression LUT
  for (unsigned int i = 0; i < 1024; i++) {
    analyticalLUT[i] = (unsigned int)(sqrt(14.94 * log(1. + i / 14.94) * i) + 0.5);
    identityLUT[i] = std::min(i, 0xffu);
  }

  hcaluncomp.clear();
  for (int i = 0; i < 4176; i++) {
    RCTdecompression decompressionTable(256, 0);
    hcaluncomp.push_back(decompressionTable);
  }

  std::vector<std::string> mainStrings;
  mainStrings.push_back("HCAL_LUT_METADATA");
  mainStrings.push_back("HCAL_LUT_CHAN_DATA");

  // ~~~~~~~~~ Cut values ~~~~~~~~~

  std::vector<std::string> metaStrings;
  metaStrings.push_back("RCTLSB");
  metaStrings.push_back("NOMINAL_GAIN");

  l1t::OMDSReader::QueryResults paramResults =
      m_omdsReader.basicQueryView(metaStrings,
                                  "CMS_HCL_HCAL_COND",
                                  "V_HCAL_LUT_METADATA_V1",
                                  "V_HCAL_LUT_METADATA_V1.TAG_NAME",
                                  m_omdsReader.basicQuery("HCAL_LUT_METADATA",
                                                          "CMS_RCT",
                                                          "HCAL_SCALE_KEY",
                                                          "HCAL_SCALE_KEY.HCAL_TAG",
                                                          m_omdsReader.singleAttribute(objectKey)));

  if (paramResults.queryFailed() || (paramResults.numberRows() != 1))  // check query successful
  {
    edm::LogError("L1-O2O") << "Problem with L1CaloHcalScale key.  Unable to find lutparam dat table";
    return std::unique_ptr<L1CaloHcalScale>();
  }

  double hcalLSB, nominal_gain;
  paramResults.fillVariable("RCTLSB", hcalLSB);
  paramResults.fillVariable("NOMINAL_GAIN", nominal_gain);

  float rctlsb = hcalLSB;

  l1t::OMDSReader::QueryResults chanKey = m_omdsReader.basicQuery("HCAL_LUT_CHAN_DATA",
                                                                  "CMS_RCT",
                                                                  "HCAL_SCALE_KEY",
                                                                  "HCAL_SCALE_KEY.HCAL_TAG",
                                                                  m_omdsReader.singleAttribute(objectKey));

  //coral::AttributeList myresult;
  //    myresult.extend(

  std::string schemaName("CMS_HCL_HCAL_COND");
  coral::ISchema& schema = m_omdsReader.dbSession().coralSession().schema(schemaName);
  coral::IQuery* query = schema.newQuery();
  ;

  std::vector<std::string> channelStrings;
  channelStrings.push_back("IPHI");
  channelStrings.push_back("IETA");
  channelStrings.push_back("DEPTH");
  channelStrings.push_back("LUT_GRANULARITY");
  channelStrings.push_back("OUTPUT_LUT_THRESHOLD");
  channelStrings.push_back("OBJECTNAME");

  std::vector<std::string>::const_iterator it = channelStrings.begin();
  std::vector<std::string>::const_iterator end = channelStrings.end();
  for (; it != end; ++it) {
    query->addToOutputList(*it);
  }

  std::string ob = "OBJECTNAME";
  coral::AttributeList myresult;
  myresult.extend("IPHI", typeid(int));
  myresult.extend("IETA", typeid(int));
  myresult.extend("DEPTH", typeid(int));
  myresult.extend("LUT_GRANULARITY", typeid(int));
  myresult.extend("OUTPUT_LUT_THRESHOLD", typeid(int));
  myresult.extend(ob, typeid(std::string));  //, typeid(std::string));

  query->defineOutput(myresult);

  query->addToTableList("V_HCAL_LUT_CHAN_DATA_V1");

  query->setCondition("V_HCAL_LUT_CHAN_DATA_V1.TAG_NAME = :" + chanKey.columnNames().front(),
                      chanKey.attributeLists().front());

  coral::ICursor& cursor = query->execute();

  // when the query goes out of scope.
  std::vector<coral::AttributeList> atts;
  while (cursor.next()) {
    atts.push_back(cursor.currentRow());
  };

  delete query;

  l1t::OMDSReader::QueryResults chanResults(channelStrings, atts);
  if (chanResults.queryFailed() || (chanResults.numberRows() == 0))  // check query successful
  {
    edm::LogError("L1-O2O") << "Problem with L1CaloHcalScale key.  Unable to find lutparam dat table nrows"
                            << chanResults.numberRows();
    return std::unique_ptr<L1CaloHcalScale>();
  }

  chanResults.attributeLists();
  for (int i = 0; i < chanResults.numberRows(); ++i) {
    std::string objectName;
    chanResults.fillVariableFromRow("OBJECTNAME", i, objectName);
    //       int
    if (objectName == "HcalTrigTowerDetId") {  //trig tower
      int ieta, iphi, depth, lutGranularity, threshold;

      chanResults.fillVariableFromRow("LUT_GRANULARITY", i, lutGranularity);
      chanResults.fillVariableFromRow("IPHI", i, iphi);
      chanResults.fillVariableFromRow("IETA", i, ieta);
      chanResults.fillVariableFromRow("DEPTH", i, depth);
      chanResults.fillVariableFromRow("OUTPUT_LUT_THRESHOLD", i, threshold);

      unsigned int outputLut[1024];

      const int tp_version = depth / 10;
      uint32_t lutId = caloTPG->getOutputLUTId(ieta, iphi, tp_version);

      double eta_low = 0., eta_high = 0.;
      theTrigTowerGeometry->towerEtaBounds(ieta, tp_version, eta_low, eta_high);
      double cosh_ieta = fabs(cosh((eta_low + eta_high) / 2.));

      if (!caloTPG->HTvalid(ieta, iphi, tp_version))
        continue;
      double factor = 0.;
      if (abs(ieta) >= theTrigTowerGeometry->firstHFTower(tp_version))
        factor = rctlsb;
      else
        factor = nominal_gain / cosh_ieta * lutGranularity;
      for (int k = 0; k < threshold; ++k)
        outputLut[k] = 0;

      for (unsigned int k = threshold; k < 1024; ++k)
        outputLut[k] = (abs(ieta) < theTrigTowerGeometry->firstHFTower(tp_version)) ? analyticalLUT[k] : identityLUT[k];

      // tpg - compressed value
      unsigned int tpg = outputLut[0];

      int low = 0;

      for (unsigned int k = 0; k < 1024; ++k) {
        if (outputLut[k] != tpg) {
          unsigned int mid = (low + k) / 2;
          hcaluncomp[lutId][tpg] = (tpg == 0 ? low : factor * mid);
          low = k;
          tpg = outputLut[k];
        }
      }
      hcaluncomp[lutId][tpg] = factor * low;
    }
  }

  auto hcalScale = std::make_unique<L1CaloHcalScale>(0);

  // XXX L1CaloHcalScale is only setup for 2x3 TP
  const int tp_version = 0;
  for (unsigned short ieta = 1; ieta <= L1CaloHcalScale::nBinEta; ++ieta) {
    for (int pos = 0; pos <= 1; pos++) {
      for (unsigned short irank = 0; irank < L1CaloHcalScale::nBinRank; ++irank) {
        int zside = (int)pow(-1, pos);
        int nphi = 0;
        double etvalue = 0.;

        for (int iphi = 1; iphi <= 72; iphi++) {
          if (!caloTPG->HTvalid(ieta, iphi, tp_version))
            continue;
          uint32_t lutId = caloTPG->getOutputLUTId(ieta, iphi, tp_version);
          nphi++;
          etvalue += (double)hcaluncomp[lutId][irank];

        }  // phi
        if (nphi > 0)
          etvalue /= nphi;

        hcalScale->setBin(irank, ieta, zside, etvalue);

      }  // rank
    }    // zside
  }      // eta

  std::stringstream s;
  s << std::setprecision(10);
  hcalScale->print(s);
  edm::LogInfo("L1CaloHcalScaleConfigOnlineProd") << s.str();
  // ------------ method called to produce the data  ------------
  return hcalScale;
}

std::unique_ptr<L1CaloHcalScale> L1CaloHcalScaleConfigOnlineProd::produce(const L1CaloHcalScaleRcd& iRecord) {
  edm::ESHandle<HcalTrigTowerGeometry> pG;
  iRecord.getRecord<CaloGeometryRecord>().get(pG);
  theTrigTowerGeometry = pG.product();

  return (L1ConfigOnlineProdBase<L1CaloHcalScaleRcd, L1CaloHcalScale>::produce(iRecord));
}

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