L1TOMDSHelper.cc

Go to the documentation of this file.

#include "DQM/L1TMonitor/interface/L1TOMDSHelper.h"

#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/Run.h"

using namespace std;

//_____________________________________________________________________
L1TOMDSHelper::L1TOMDSHelper() { m_omdsReader = nullptr; }

//_____________________________________________________________________
L1TOMDSHelper::~L1TOMDSHelper() { delete m_omdsReader; }

//_____________________________________________________________________
bool L1TOMDSHelper::connect(string iOracleDB, string iPathCondDB, int &error) {
  // Handeling inputs
  m_oracleDB = iOracleDB;
  m_pathCondDB = iPathCondDB;
  error = NO_ERROR;

  // Initializing variables
  bool SessionOpen = false;
  bool out = false;

  m_omdsReader = new l1t::OMDSReader();
  m_omdsReader->connect(m_oracleDB, m_pathCondDB);

  // Testing session
  if (!m_omdsReader->dbSession().connectionString().empty()) {
    SessionOpen = true;
  }

  // Defining output and error message if needed
  if (SessionOpen) {
    out = true;
  } else {
    error = WARNING_DB_CONN_FAILED;
  }

  return out;
}

//_____________________________________________________________________
map<string, WbMTriggerXSecFit> L1TOMDSHelper::getWbMTriggerXsecFits(string iTable, int &error) {
  map<string, WbMTriggerXSecFit> out;
  const l1t::OMDSReader::QueryResults qresults;
  error = NO_ERROR;

  // Parameters
  string qSchema = "CMS_WBM";

  // Fields to retrive in the query
  vector<string> qStrings;
  qStrings.push_back("BIT");
  qStrings.push_back("NAME");
  qStrings.push_back("PM1");  //Inverse
  qStrings.push_back("P0");   //Constant
  qStrings.push_back("P1");   //Linear
  qStrings.push_back("P2");   //Quadratic

  l1t::OMDSReader::QueryResults paramResults = m_omdsReader->basicQuery(qStrings, qSchema, iTable, "", qresults);

  if (paramResults.queryFailed()) {
    error = WARNING_DB_QUERY_FAILED;
  } else {
    for (int i = 0; i < paramResults.numberRows(); ++i) {
      WbMTriggerXSecFit tFit;
      tFit.fitFunction = "[0]/x+[1]+[2]*x+[3]*x*x";  // Fitting function hardcoded for now

      string tBitName;

      paramResults.fillVariableFromRow("BIT", i, tFit.bitNumber);
      paramResults.fillVariableFromRow("NAME", i, tBitName);
      paramResults.fillVariableFromRow("PM1", i, tFit.pm1);
      paramResults.fillVariableFromRow("P0", i, tFit.p0);
      paramResults.fillVariableFromRow("P1", i, tFit.p1);
      paramResults.fillVariableFromRow("P2", i, tFit.p2);

      tFit.bitName = tBitName;

      out[tBitName] = tFit;
    }
  }

  return out;
}

//_____________________________________________________________________
map<string, WbMTriggerXSecFit> L1TOMDSHelper::getWbMAlgoXsecFits(int &error) {
  return getWbMTriggerXsecFits("LEVEL1_ALGO_CROSS_SECTION", error);
}

//_____________________________________________________________________
map<string, WbMTriggerXSecFit> L1TOMDSHelper::getWbMTechXsecFits(int &error) {
  return getWbMTriggerXsecFits("LEVEL1_TECH_CROSS_SECTION", error);
}

//_____________________________________________________________________
int L1TOMDSHelper::getNumberCollidingBunches(int lhcFillNumber, int &error) {
  int nCollidingBunches = 0;
  error = NO_ERROR;

  // Parameters
  string rtlSchema = "CMS_RUNTIME_LOGGER";
  string table = "FILL_INITLUMIPERBUNCH";
  string atribute1 = "LHCFILL";

  // Fields to retrive in the query
  vector<std::string> qStrings;
  qStrings.push_back("BUNCH");
  qStrings.push_back("BEAM1CONFIG");
  qStrings.push_back("BEAM2CONFIG");

  l1t::OMDSReader::QueryResults qResults =
      m_omdsReader->basicQuery(qStrings, rtlSchema, table, atribute1, m_omdsReader->singleAttribute(lhcFillNumber));

  // Check query successful
  if (qResults.queryFailed()) {
    error = WARNING_DB_QUERY_FAILED;
  } else {
    if (qResults.numberRows() != 3564) {
      error = WARNING_DB_INCORRECT_NBUNCHES;
    } else {
      // Now we count the number of bunches with both beam 1 and 2 configured
      for (int i = 0; i < qResults.numberRows(); ++i) {
        int bunch;
        bool beam1config, beam2config;
        qResults.fillVariableFromRow("BUNCH", i, bunch);
        qResults.fillVariableFromRow("BEAM1CONFIG", i, beam1config);
        qResults.fillVariableFromRow("BEAM2CONFIG", i, beam2config);

        if (beam1config && beam2config) {
          nCollidingBunches++;
        }
      }
    }
  }

  return nCollidingBunches;
}

//_____________________________________________________________________
BeamConfiguration L1TOMDSHelper::getBeamConfiguration(int lhcFillNumber, int &error) {
  BeamConfiguration bConfig;
  error = NO_ERROR;

  // Fields to retrive in the query
  string rtlSchema = "CMS_RUNTIME_LOGGER";
  string table = "FILL_INITLUMIPERBUNCH";
  string atribute1 = "LHCFILL";

  // Setting the strings we want to recover from OMDS
  vector<std::string> qStrings;
  qStrings.push_back("BUNCH");
  qStrings.push_back("BEAM1CONFIG");
  qStrings.push_back("BEAM2CONFIG");

  l1t::OMDSReader::QueryResults qResults =
      m_omdsReader->basicQuery(qStrings, rtlSchema, table, atribute1, m_omdsReader->singleAttribute(lhcFillNumber));

  if (qResults.queryFailed()) {
    error = WARNING_DB_QUERY_FAILED;
  } else {
    if (qResults.numberRows() != 3564) {
      error = WARNING_DB_INCORRECT_NBUNCHES;
    } else {
      bConfig.m_valid = true;

      int nCollidingBunches = 0;

      for (int i = 0; i < qResults.numberRows(); ++i) {
        int bunch;
        bool beam1config, beam2config;
        qResults.fillVariableFromRow("BUNCH", i, bunch);
        qResults.fillVariableFromRow("BEAM1CONFIG", i, beam1config);
        qResults.fillVariableFromRow("BEAM2CONFIG", i, beam2config);

        if (beam1config) {
          bConfig.beam1.push_back(true);
        } else {
          bConfig.beam1.push_back(false);
        }

        if (beam2config) {
          bConfig.beam2.push_back(true);
        } else {
          bConfig.beam2.push_back(false);
        }

        if (beam1config && beam2config) {
          nCollidingBunches++;
        }
      }

      bConfig.nCollidingBunches = nCollidingBunches;
    }
  }

  return bConfig;
}

//_____________________________________________________________________
vector<bool> L1TOMDSHelper::getBunchStructure(int lhcFillNumber, int &error) {
  vector<bool> BunchStructure;
  error = NO_ERROR;

  // Fields to retrive in the query
  string rtlSchema = "CMS_RUNTIME_LOGGER";
  string table = "FILL_INITLUMIPERBUNCH";
  string atribute1 = "LHCFILL";

  // Setting the strings we want to recover from OMDS
  vector<std::string> qStrings;
  qStrings.push_back("BUNCH");
  qStrings.push_back("BEAM1CONFIG");
  qStrings.push_back("BEAM2CONFIG");

  l1t::OMDSReader::QueryResults qResults =
      m_omdsReader->basicQuery(qStrings, rtlSchema, table, atribute1, m_omdsReader->singleAttribute(lhcFillNumber));

  if (qResults.queryFailed()) {
    error = WARNING_DB_QUERY_FAILED;
  } else {
    if (qResults.numberRows() != 3564) {
      error = WARNING_DB_INCORRECT_NBUNCHES;
    } else {
      for (int i = 0; i < qResults.numberRows(); ++i) {
        int bunch;
        bool beam1config, beam2config;
        qResults.fillVariableFromRow("BUNCH", i, bunch);
        qResults.fillVariableFromRow("BEAM1CONFIG", i, beam1config);
        qResults.fillVariableFromRow("BEAM2CONFIG", i, beam2config);

        if (beam1config && beam2config) {
          BunchStructure.push_back(true);
        } else {
          BunchStructure.push_back(false);
        }
      }
    }
  }

  return BunchStructure;
}

//_____________________________________________________________________
vector<float> L1TOMDSHelper::getInitBunchLumi(int lhcFillNumber, int &error) {
  vector<float> InitBunchLumi;
  error = NO_ERROR;

  // Fields to retrive in the query
  string rtlSchema = "CMS_RUNTIME_LOGGER";
  string table = "FILL_INITLUMIPERBUNCH";
  string atribute1 = "LHCFILL";

  // Setting the strings we want to recover from OMDS
  vector<std::string> qStrings;
  qStrings.push_back("BUNCH");
  qStrings.push_back("INITBUNCHLUMI");

  l1t::OMDSReader::QueryResults qResults =
      m_omdsReader->basicQuery(qStrings, rtlSchema, table, atribute1, m_omdsReader->singleAttribute(lhcFillNumber));

  if (qResults.queryFailed()) {
    error = WARNING_DB_QUERY_FAILED;
  } else {
    if (qResults.numberRows() != 3564) {
      error = WARNING_DB_INCORRECT_NBUNCHES;
    } else {
      for (int i = 0; i < qResults.numberRows(); ++i) {
        int bunch;
        float initbunchlumi;
        qResults.fillVariableFromRow("BUNCH", i, bunch);
        qResults.fillVariableFromRow("INITBUNCHLUMI", i, initbunchlumi);

        InitBunchLumi.push_back(initbunchlumi);
      }
    }
  }

  return InitBunchLumi;
}

//_____________________________________________________________________
vector<double> L1TOMDSHelper::getRelativeBunchLumi(int lhcFillNumber, int &error) {
  vector<double> RelativeBunchLumi;
  error = NO_ERROR;

  // Fields to retrive in the query
  string rtlSchema = "CMS_RUNTIME_LOGGER";
  string table = "FILL_INITLUMIPERBUNCH";
  string atribute1 = "LHCFILL";

  // Setting the strings we want to recover from OMDS
  vector<std::string> qStrings;
  qStrings.push_back("BUNCH");
  qStrings.push_back("INITBUNCHLUMI");

  l1t::OMDSReader::QueryResults qResults =
      m_omdsReader->basicQuery(qStrings, rtlSchema, table, atribute1, m_omdsReader->singleAttribute(lhcFillNumber));

  if (qResults.queryFailed()) {
    error = WARNING_DB_QUERY_FAILED;
  } else {
    if (qResults.numberRows() != 3564) {
      error = WARNING_DB_INCORRECT_NBUNCHES;
    } else {
      //-> Get the inicial bunch luminosity add calculate the total luminosity of the fill
      double InitTotalLumi = 0;
      vector<float> InitBunchLumi;

      for (int i = 0; i < qResults.numberRows(); ++i) {
        int bunch;
        float initbunchlumi;
        qResults.fillVariableFromRow("BUNCH", i, bunch);
        qResults.fillVariableFromRow("INITBUNCHLUMI", i, initbunchlumi);

        InitTotalLumi += initbunchlumi;
        InitBunchLumi.push_back(initbunchlumi);
      }

      //-> We calculate the relative luminosity for each bunch
      for (unsigned int i = 0; i < InitBunchLumi.size(); i++) {
        RelativeBunchLumi.push_back(InitBunchLumi[i] / InitTotalLumi);
      }
    }
  }

  return RelativeBunchLumi;
}

//_____________________________________________________________________
string L1TOMDSHelper::enumToStringError(int iObject) {
  string out;

  switch (iObject) {
    case NO_ERROR:
      out = "NO_ERROR";
      break;
    case WARNING_DB_CONN_FAILED:
      out = "WARNING_DB_CONN_FAILED";
      break;
    case WARNING_DB_QUERY_FAILED:
      out = "WARNING_DB_QUERY_FAILED";
      break;
    case WARNING_DB_INCORRECT_NBUNCHES:
      out = "WARNING_DB_INCORRECT_NBUNCHES";
      break;
    default:
      out = "UNKNOWN";
      break;
  };

  return out;
}