RPCLBLinkMapHandler.cc

Go to the documentation of this file.

#include "CondTools/RPC/interface/RPCLBLinkMapHandler.h"

#include <cstdint>
#include <fstream>
#include <memory>
#include <sstream>

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/Exception.h"

#include "CondCore/CondDB/interface/Types.h"

#include "RelationalAccess/ICursor.h"
#include "RelationalAccess/IQuery.h"
#include "RelationalAccess/IQueryDefinition.h"
#include "RelationalAccess/ISchema.h"
#include "RelationalAccess/ISessionProxy.h"
#include "RelationalAccess/ITransaction.h"

#include "CoralBase/Attribute.h"
#include "CoralBase/AttributeList.h"

#include "CondTools/RPC/interface/RPCLBLinkNameParser.h"

RPCDetId RPCLBLinkMapHandler::getRPCDetId(
    int region, int disk_or_wheel, int layer, int sector, std::string subsector_string, std::string partition) {
  int station(0), ring(0), subsector(0), roll(0);
  // region well-defined
  if (!region) {  // barrel
    switch (layer) {
      case 1:
      case 2:
        station = 1;
        break;
      case 3:
      case 4:
        station = 2;
        layer -= 2;
        break;
      case 5:
        station = 3;
        layer = 1;
        break;
      case 6:
        station = 4;
        layer = 1;
        break;
    }
    ring = disk_or_wheel;
    // sector well-defined
    subsector = 1;
    if (subsector_string == "+" && (station == 3 || (station == 4 && (sector != 4 && sector != 9 && sector != 11))))
      subsector = 2;
    if (station == 4 && sector == 4) {
      if (subsector_string == "-")
        subsector = 2;
      else if (subsector_string == "+")
        subsector = 3;
      else if (subsector_string == "++")
        subsector = 4;
    }
  } else {  // endcap
    station = std::abs(disk_or_wheel);
    ring = layer;
    layer = 1;
    if (ring > 1 || station == 1) {
      subsector = (sector - 1) % 6 + 1;
      sector = (sector - 1) / 6 + 1;
    } else {
      subsector = (sector - 1) % 3 + 1;
      sector = (sector - 1) / 3 + 1;
    }
  }
  // roll
  if (partition == "Backward" || partition == "A")
    roll = 1;
  else if (partition == "Central" || partition == "B")
    roll = 2;
  else if (partition == "Forward" || partition == "C")
    roll = 3;
  else if (partition == "D")
    roll = 4;
  else
    throw cms::Exception("RPCLBLinkMapHandler") << "Unexpected partition name " << partition;

  return RPCDetId(region, ring, station, sector, layer, subsector, roll);
}

RPCLBLinkMapHandler::RPCLBLinkMapHandler(edm::ParameterSet const& config)
    : id_(config.getParameter<std::string>("identifier")),
      data_tag_(config.getParameter<std::string>("dataTag")),
      since_run_(config.getParameter<unsigned long long>("sinceRun")),
      txt_file_(config.getUntrackedParameter<std::string>("txtFile", "")),
      connect_(config.getParameter<std::string>("connect")) {
  edm::LogInfo("RPCDCCLinkMapHandler") << "Configuring Input Connection";
  connection_.setParameters(config.getParameter<edm::ParameterSet>("DBParameters"));
  connection_.configure();
}

RPCLBLinkMapHandler::~RPCLBLinkMapHandler() {}

void RPCLBLinkMapHandler::getNewObjects() {
  edm::LogInfo("RPCLBLinkMapHandler") << "getNewObjects";
  cond::TagInfo const& tag_info = tagInfo();
  if (since_run_ < tag_info.lastInterval.first)
    throw cms::Exception("RPCLBLinkMapHandler") << "Refuse to create RPCLBLinkMap for run " << since_run_
                                                << ", older than most recent tag" << tag_info.lastInterval.first;

  edm::LogInfo("RPCDCCLinkMapHandler") << "Opening read-only Input Session";
  auto input_session = connection_.createCoralSession(connect_, false);  // writeCapable
  edm::LogInfo("RPCDCCLinkMapHandler") << "Started Input Transaction";
  input_session->transaction().start(true);  // readOnly

  std::unique_ptr<coral::IQuery> query(input_session->schema("CMS_RPC_CONF").newQuery());
  query->addToTableList("BOARD");
  query->addToTableList("CHAMBERSTRIP");
  query->addToTableList("CHAMBERLOCATION");
  query->addToTableList("FEBLOCATION");
  query->addToTableList("FEBCONNECTOR");
  coral::IQueryDefinition& subquery_min_channel(query->defineSubQuery("MIN_CHANNEL"));
  query->addToTableList("MIN_CHANNEL");
  coral::IQueryDefinition& subquery_max_strip(query->defineSubQuery("MAX_STRIP"));
  query->addToTableList("MAX_STRIP");

  query->addToOutputList("BOARD.NAME", "LB_NAME");
  query->addToOutputList("FEBCONNECTOR.LINKBOARDINPUTNUM", "CONNECTOR");
  query->addToOutputList("CHAMBERSTRIP.CHAMBERSTRIPNUMBER", "FIRST_STRIP");
  query->addToOutputList("CAST(DECODE(SIGN(MAX_STRIP.STRIP - CHAMBERSTRIP.CHAMBERSTRIPNUMBER), 1, 1, -1) AS INTEGER)",
                         "SLOPE");
  query->addToOutputList("MIN_CHANNEL.CHANNELS", "CHANNELS");
  query->addToOutputList(
      "CAST(DECODE(CHAMBERLOCATION.BARRELORENDCAP, 'Barrel', 0, DECODE(SIGN(CHAMBERLOCATION.DISKORWHEEL), 1, 1, -1)) "
      "AS INTEGER)",
      "REGION");
  query->addToOutputList("CHAMBERLOCATION.DISKORWHEEL", "DISKORWHEEL");
  query->addToOutputList("CHAMBERLOCATION.LAYER", "LAYER");
  query->addToOutputList("CHAMBERLOCATION.SECTOR", "SECTOR");
  query->addToOutputList("CHAMBERLOCATION.SUBSECTOR", "SUBSECTOR");
  query->addToOutputList("FEBLOCATION.FEBLOCALETAPARTITION", "ETAPARTITION");

  subquery_min_channel.addToTableList("CHAMBERSTRIP");
  subquery_min_channel.addToOutputList("FC_FEBCONNECTORID");
  subquery_min_channel.addToOutputList("MIN(CABLECHANNELNUM)", "CHANNEL");
  subquery_min_channel.addToOutputList("CAST(SUM(POWER(2, CABLECHANNELNUM-1)) AS INTEGER)", "CHANNELS");
  subquery_min_channel.groupBy("FC_FEBCONNECTORID");
  coral::AttributeList subquery_min_channel_condition_data;
  subquery_min_channel.setCondition("CABLECHANNELNUM IS NOT NULL", subquery_min_channel_condition_data);

  subquery_max_strip.addToTableList("CHAMBERSTRIP");
  coral::IQueryDefinition& subquery_max_channel(subquery_max_strip.defineSubQuery("MAX_CHANNEL"));
  subquery_max_strip.addToTableList("MAX_CHANNEL");
  subquery_max_strip.addToOutputList("CHAMBERSTRIP.FC_FEBCONNECTORID", "FC_FEBCONNECTORID");
  subquery_max_strip.addToOutputList("CHAMBERSTRIP.CHAMBERSTRIPNUMBER", "STRIP");
  coral::AttributeList subquery_max_strip_condition_data;
  subquery_max_strip.setCondition(
      "CHAMBERSTRIP.FC_FEBCONNECTORID=MAX_CHANNEL.FC_FEBCONNECTORID"
      " AND CHAMBERSTRIP.CABLECHANNELNUM=MAX_CHANNEL.CHANNEL",
      subquery_max_strip_condition_data);

  subquery_max_channel.addToTableList("CHAMBERSTRIP");
  subquery_max_channel.addToOutputList("FC_FEBCONNECTORID");
  subquery_max_channel.addToOutputList("MAX(CABLECHANNELNUM)", "CHANNEL");
  subquery_max_channel.groupBy("FC_FEBCONNECTORID");
  coral::AttributeList subquery_max_channel_condition_data;
  subquery_max_channel.setCondition("CABLECHANNELNUM IS NOT NULL", subquery_max_channel_condition_data);

  coral::AttributeList query_condition_data;
  query->setCondition(
      "CHAMBERSTRIP.FC_FEBCONNECTORID=MIN_CHANNEL.FC_FEBCONNECTORID"
      " AND CHAMBERSTRIP.CABLECHANNELNUM=MIN_CHANNEL.CHANNEL"
      " AND CHAMBERSTRIP.FC_FEBCONNECTORID=MAX_STRIP.FC_FEBCONNECTORID"
      " AND CHAMBERSTRIP.FC_FEBCONNECTORID=FEBCONNECTOR.FEBCONNECTORID"
      " AND FEBCONNECTOR.FL_FEBLOCATIONID=FEBLOCATION.FEBLOCATIONID"
      " AND BOARD.BOARDID=FEBLOCATION.LB_LINKBOARDID"
      " AND CHAMBERLOCATION.CHAMBERLOCATIONID=FEBLOCATION.CL_CHAMBERLOCATIONID",
      query_condition_data);

  std::string lb_name("");
  int first_strip(0), slope(1);
  std::uint16_t channels(0x0);

  std::unique_ptr<RPCLBLinkMap> lb_link_map_object(new RPCLBLinkMap());
  RPCLBLinkMap::map_type& lb_link_map = lb_link_map_object->getMap();
  RPCLBLink lb_link;
  RPCDetId det_id;
  std::string subsector;

  edm::LogInfo("RPCLBLinkMapHandler") << "Execute query";
  coral::ICursor& cursor(query->execute());
  while (cursor.next()) {
    coral::AttributeList const& row(cursor.currentRow());

    // RPCLBLink
    lb_name = row["LB_NAME"].data<std::string>();
    RPCLBLinkNameParser::parse(lb_name, lb_link);
    if (lb_name != lb_link.getName())
      edm::LogWarning("RPCLBLinkMapHandler") << "Mismatch LinkBoard Name: " << lb_name << " vs " << lb_link;
    lb_link.setConnector(row["CONNECTOR"].data<short>() - 1);  // 1-6 to 0-5

    // RPCDetId
    if (row["SUBSECTOR"].isNull())
      subsector = "";
    else
      subsector = row["SUBSECTOR"].data<std::string>();
    det_id = getRPCDetId(row["REGION"].data<long long>(),
                         row["DISKORWHEEL"].data<short>(),
                         row["LAYER"].data<short>(),
                         row["SECTOR"].data<short>(),
                         subsector,
                         row["ETAPARTITION"].data<std::string>());

    // RPCFebConnector
    first_strip = row["FIRST_STRIP"].data<int>();
    slope = row["SLOPE"].data<long long>();
    channels = (std::uint16_t)(row["CHANNELS"].data<long long>());

    lb_link_map.insert(
        std::pair<RPCLBLink, RPCFebConnector>(lb_link, RPCFebConnector(det_id, first_strip, slope, channels)));
  }
  cursor.close();

  input_session->transaction().commit();

  if (!txt_file_.empty()) {
    edm::LogInfo("RPCLBLinkMapHandler") << "Fill txtFile";
    std::ofstream ofstream(txt_file_);
    for (auto const& link_connector : lb_link_map) {
      ofstream << link_connector.first << ": " << link_connector.second << std::endl;
    }
  }

  edm::LogInfo("RPCLBLinkMapHandler") << "Add to transfer list";
  m_to_transfer.push_back(std::make_pair(lb_link_map_object.release(), since_run_));
}

std::string RPCLBLinkMapHandler::id() const { return id_; }