SiStripMiscalibrateHelper.cc

Go to the documentation of this file.

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "CondFormats/SiStripObjects/interface/SiStripSummary.h"
#include "CommonTools/TrackerMap/interface/TmModule.h"
#include "CondTools/SiStrip/interface/SiStripMiscalibrateHelper.h"
#include "CondCore/SiStripPlugins/interface/SiStripPayloadInspectorHelper.h"

/*--------------------------------------------------------------------*/
sistripsummary::TrackerRegion SiStripMiscalibrate::getRegionFromString(std::string region)
/*--------------------------------------------------------------------*/
{
  std::map<std::string, sistripsummary::TrackerRegion> mapping = {
      {"Tracker", sistripsummary::TRACKER}, {"TIB", sistripsummary::TIB},       {"TIB_1", sistripsummary::TIB_1},
      {"TIB_2", sistripsummary::TIB_2},     {"TIB_3", sistripsummary::TIB_3},   {"TIB_4", sistripsummary::TIB_4},
      {"TID", sistripsummary::TID},         {"TIDP", sistripsummary::TIDP},     {"TIDP_1", sistripsummary::TIDP_1},
      {"TIDP_2", sistripsummary::TIDP_2},   {"TIDP_3", sistripsummary::TIDP_3}, {"TIDM", sistripsummary::TIDM},
      {"TIDM_1", sistripsummary::TIDM_1},   {"TIDM_2", sistripsummary::TIDM_2}, {"TIDM_3", sistripsummary::TIDM_3},
      {"TOB", sistripsummary::TOB},         {"TOB_1", sistripsummary::TOB_1},   {"TOB_2", sistripsummary::TOB_2},
      {"TOB_3", sistripsummary::TOB_3},     {"TOB_4", sistripsummary::TOB_4},   {"TOB_5", sistripsummary::TOB_5},
      {"TOB_6", sistripsummary::TOB_6},     {"TEC", sistripsummary::TEC},       {"TECP", sistripsummary::TECP},
      {"TECP_1", sistripsummary::TECP_1},   {"TECP_2", sistripsummary::TECP_2}, {"TECP_3", sistripsummary::TECP_3},
      {"TECP_4", sistripsummary::TECP_4},   {"TECP_5", sistripsummary::TECP_5}, {"TECP_6", sistripsummary::TECP_6},
      {"TECP_7", sistripsummary::TECP_7},   {"TECP_8", sistripsummary::TECP_8}, {"TECP_9", sistripsummary::TECP_9},
      {"TECM", sistripsummary::TECM},       {"TECM_1", sistripsummary::TECM_1}, {"TECM_2", sistripsummary::TECM_2},
      {"TECM_3", sistripsummary::TECM_3},   {"TECM_4", sistripsummary::TECM_4}, {"TECM_5", sistripsummary::TECM_5},
      {"TECM_6", sistripsummary::TECM_6},   {"TECM_7", sistripsummary::TECM_7}, {"TECM_8", sistripsummary::TECM_8},
      {"TECM_9", sistripsummary::TECM_9}};

  if (mapping.find(region) == mapping.end()) {
    edm::LogError("SiStripMiscalibrate") << "@SUB=analyze"
                                         << "Unknown partition: " << region;
    throw cms::Exception("Invalid Partition passed");
  } else {
    return mapping[region];
  }
}

/*--------------------------------------------------------------------*/
std::vector<sistripsummary::TrackerRegion> SiStripMiscalibrate::getRegionsFromDetId(
    const TrackerTopology* m_trackerTopo, DetId detid)
/*--------------------------------------------------------------------*/
{
  int layer = 0;
  int side = 0;
  int subdet = 0;
  int detCode = 0;

  std::vector<sistripsummary::TrackerRegion> ret;

  switch (detid.subdetId()) {
    case StripSubdetector::TIB:
      layer = m_trackerTopo->tibLayer(detid);
      subdet = 1;
      break;
    case StripSubdetector::TOB:
      layer = m_trackerTopo->tobLayer(detid);
      subdet = 2;
      break;
    case StripSubdetector::TID:
      // is this module in TID+ or TID-?
      layer = m_trackerTopo->tidWheel(detid);
      side = m_trackerTopo->tidSide(detid);
      subdet = 3 * 10 + side;
      break;
    case StripSubdetector::TEC:
      // is this module in TEC+ or TEC-?
      layer = m_trackerTopo->tecWheel(detid);
      side = m_trackerTopo->tecSide(detid);
      subdet = 4 * 10 + side;
      break;
  }

  detCode = (subdet * 10) + layer;

  ret.push_back(static_cast<sistripsummary::TrackerRegion>(detCode));

  if (subdet / 10 > 0) {
    ret.push_back(static_cast<sistripsummary::TrackerRegion>(subdet / 10));
  }

  ret.push_back(static_cast<sistripsummary::TrackerRegion>(subdet));
  ret.push_back(sistripsummary::TRACKER);

  return ret;
}

/*--------------------------------------------------------------------*/
std::pair<float, float> SiStripMiscalibrate::getTruncatedRange(const TrackerMap* theMap) {
  /*--------------------------------------------------------------------*/

  // ------------ trim the tracker map  ------------

  auto map = theMap->smoduleMap;
  std::map<unsigned int, float> info_per_detid;
  for (int layer = 1; layer < 44; layer++) {
    for (int ring = theMap->firstRing[layer - 1]; ring < theMap->ntotRing[layer - 1] + theMap->firstRing[layer - 1];
         ring++) {
      for (int module = 1; module < 200; module++) {
        int key = layer * 100000 + ring * 1000 + module;
        TmModule* mod = map[key];
        if (mod != nullptr && !mod->notInUse() && mod->count > 0) {
          info_per_detid[key] = mod->value;
        }
      }  // loop on modules
    }  // loop on ring
  }  // loop on layers

  auto range = SiStripPI::getTheRange(info_per_detid, 2);
  return range;
}