SiPixelFEDChannelContainer_PayloadInspector.cc

Go to the documentation of this file.

 
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "CondCore/CondDB/interface/ConnectionPool.h"
 
#include "CondCore/Utilities/interface/PayloadInspectorModule.h"
#include "CondCore/Utilities/interface/PayloadInspector.h"
#include "CondCore/CondDB/interface/Time.h"
#include "CondCore/SiPixelPlugins/interface/SiPixelPayloadInspectorHelper.h"
#include "FWCore/ParameterSet/interface/FileInPath.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "CalibTracker/StandaloneTrackerTopology/interface/StandaloneTrackerTopology.h"
 
#include "CondFormats/SiPixelObjects/interface/CablingPathToDetUnit.h"
#include "CondFormats/SiPixelObjects/interface/LocalPixel.h"
#include "CondFormats/SiPixelObjects/interface/PixelIndices.h"
#include "CondFormats/SiPixelObjects/interface/PixelROC.h"
#include "CondFormats/SiPixelObjects/interface/SiPixelFedCabling.h"
#include "CondFormats/SiPixelObjects/interface/SiPixelFedCablingMap.h"
#include "CondFormats/SiPixelObjects/interface/SiPixelFedCablingTree.h"
#include "CondFormats/SiPixelObjects/interface/SiPixelFrameConverter.h"
 
// the data format of the condition to be inspected
#include "CondFormats/SiPixelObjects/interface/SiPixelFEDChannelContainer.h"
#include "DataFormats/SiPixelDetId/interface/PixelSubdetector.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "DQM/TrackerRemapper/interface/Phase1PixelROCMaps.h"
 
#include <memory>
#include <sstream>
#include <iostream>
#include <fmt/printf.h>
 
// include ROOT
#include "TH2F.h"
#include "TLegend.h"
#include "TCanvas.h"
#include "TLine.h"
#include "TGraph.h"
#include "TGaxis.h"
#include "TStyle.h"
#include "TLatex.h"
#include "TPave.h"
#include "TPaveStats.h"
 
namespace {
 
  using namespace cond::payloadInspector;
 
  /************************************************
  1d histogram of SiPixelFEDChannelContainer of 1 IOV 
  *************************************************/
 
  class SiPixelFEDChannelContainerTest : public PlotImage<SiPixelFEDChannelContainer, SINGLE_IOV> {
  public:
    SiPixelFEDChannelContainerTest()
        : PlotImage<SiPixelFEDChannelContainer, SINGLE_IOV>("SiPixelFEDChannelContainer scenarios count"),
          m_trackerTopo{StandaloneTrackerTopology::fromTrackerParametersXMLFile(
              edm::FileInPath("Geometry/TrackerCommonData/data/PhaseI/trackerParameters.xml").fullPath())} {
      // for inputs
      PlotBase::addInputParam("Scenarios");
 
      // hardcoded connection to the MC cabling tag, though luck
      m_condDbCabling = "frontier://FrontierProd/CMS_CONDITIONS";
      m_CablingTagName = "SiPixelFedCablingMap_phase1_v7";
 
      m_connectionPool.setParameters(m_connectionPset);
      m_connectionPool.configure();
    }
 
    bool fill() override {
      std::vector<std::string> the_scenarios = {};
 
      auto paramValues = PlotBase::inputParamValues();
      auto ip = paramValues.find("Scenarios");
      if (ip != paramValues.end()) {
        auto input = boost::lexical_cast<std::string>(ip->second);
        typedef boost::tokenizer<boost::char_separator<char>> tokenizer;
        boost::char_separator<char> sep{","};
        tokenizer tok{input, sep};
        for (const auto& t : tok) {
          the_scenarios.push_back(t);
        }
      } else {
        edm::LogWarning("SiPixelFEDChannelContainerTest")
            << "\n WARNING!!!! \n The needed parameter Scenarios has not been passed. Will use all the scenarios in "
               "the file!"
            << "\n Buckle your seatbelts... this might take a while... \n\n";
        the_scenarios.push_back("all");
      }
 
      Phase1PixelROCMaps theROCMap("");
 
      auto tag = PlotBase::getTag<0>();
      auto iov = tag.iovs.front();
 
      // open db session for the cabling map
      edm::LogPrint("SiPixelFEDChannelContainerTest") << "[SiPixelFEDChannelContainerTest::" << __func__ << "] "
                                                      << "Query the condition database " << m_condDbCabling;
 
      cond::persistency::Session condDbSession = m_connectionPool.createSession(m_condDbCabling);
      condDbSession.transaction().start(true);
 
      // query the database
      edm::LogPrint("SiPixelFEDChannelContainerTest") << "[SiPixelFEDChannelContainerTest::" << __func__ << "] "
                                                      << "Reading IOVs from tag " << m_CablingTagName;
 
      const auto MIN_VAL = cond::timeTypeSpecs[cond::runnumber].beginValue;
      const auto MAX_VAL = cond::timeTypeSpecs[cond::runnumber].endValue;
 
      // get the list of payloads for the Cabling Map
      std::vector<std::tuple<cond::Time_t, cond::Hash>> m_cabling_iovs;
      condDbSession.readIov(m_CablingTagName).selectRange(MIN_VAL, MAX_VAL, m_cabling_iovs);
 
      std::vector<unsigned int> listOfCablingIOVs;
      std::transform(m_cabling_iovs.begin(),
                     m_cabling_iovs.end(),
                     std::back_inserter(listOfCablingIOVs),
                     [](std::tuple<cond::Time_t, cond::Hash> myIOV2) -> unsigned int { return std::get<0>(myIOV2); });
 
      edm::LogPrint("SiPixelFEDChannelContainerTest")
          << " Number of SiPixelFedCablngMap payloads: " << listOfCablingIOVs.size() << std::endl;
 
      auto it = std::find(
          listOfCablingIOVs.begin(), listOfCablingIOVs.end(), closest_from_below(listOfCablingIOVs, std::get<0>(iov)));
      int index = std::distance(listOfCablingIOVs.begin(), it);
 
      edm::LogPrint("SiPixelFEDChannelContainerTest")
          << " using the SiPixelFedCablingMap with hash: " << std::get<1>(m_cabling_iovs.at(index)) << std::endl;
 
      auto theCablingMap = condDbSession.fetchPayload<SiPixelFedCablingMap>(std::get<1>(m_cabling_iovs.at(index)));
      theCablingMap->initializeRocs();
      // auto theCablingTree = (*theCablingMap).cablingTree();
 
      //auto map = theCablingMap->det2fedMap();
      //for (const auto &element : map){
      //    std::cout << element.first << " " << element.second << std::endl;
      //}
 
      std::shared_ptr<SiPixelFEDChannelContainer> payload = fetchPayload(std::get<1>(iov));
      const auto& scenarioMap = payload->getScenarioMap();
 
      auto pIndexConverter = PixelIndices(numColumns, numRows);
 
      for (const auto& scenario : scenarioMap) {
        std::string scenName = scenario.first;
 
        if (std::find_if(the_scenarios.begin(), the_scenarios.end(), compareKeys(scenName)) != the_scenarios.end() ||
            the_scenarios[0] == "all") {
          edm::LogPrint("SiPixelFEDChannelContainerTest") << "\t Found Scenario: " << scenName << " ==> dumping it";
        } else {
          continue;
        }
 
        //if (strcmp(scenName.c_str(),"320824_103") != 0) continue;
 
        const auto& theDetSetBadPixelFedChannels = payload->getDetSetBadPixelFedChannels(scenName);
        for (const auto& disabledChannels : *theDetSetBadPixelFedChannels) {
          const auto t_detid = disabledChannels.detId();
          int subid = DetId(t_detid).subdetId();
          LogDebug("SiPixelFEDChannelContainerTest") << fmt::sprintf("DetId : %i \n", t_detid) << std::endl;
 
          std::bitset<16> badRocsFromFEDChannels;
 
          for (const auto& ch : disabledChannels) {
            std::string toOut_ = fmt::sprintf("fed : %i | link : %2i | roc_first : %2i | roc_last: %2i \n",
                                              ch.fed,
                                              ch.link,
                                              ch.roc_first,
                                              ch.roc_last);
 
            LogDebug("SiPixelFEDChannelContainerTest") << toOut_ << std::endl;
            const std::vector<sipixelobjects::CablingPathToDetUnit>& path =
                theCablingMap->pathToDetUnit(disabledChannels.detId());
            for (unsigned int i_roc = ch.roc_first; i_roc <= ch.roc_last; ++i_roc) {
              for (const auto p : path) {
                const sipixelobjects::PixelROC* myroc = theCablingMap->findItem(p);
                if (myroc->idInDetUnit() == static_cast<unsigned int>(i_roc)) {
                  sipixelobjects::LocalPixel::RocRowCol local = {39, 25};  //corresponding to center of ROC row,col
                  sipixelobjects::GlobalPixel global = myroc->toGlobal(sipixelobjects::LocalPixel(local));
                  int chipIndex(0), colROC(0), rowROC(0);
 
                  pIndexConverter.transformToROC(global.col, global.row, chipIndex, colROC, rowROC);
 
                  LogDebug("SiPixelFEDChannelContainerTest")
                      << " => i_roc:" << i_roc << "  " << global.col << "-" << global.row << " | => " << chipIndex
                      << " : (" << colROC << "," << rowROC << ")" << std::endl;
 
                  badRocsFromFEDChannels[chipIndex] = true;
                }
              }
            }
          }
 
          LogDebug("SiPixelFEDChannelContainerTest") << badRocsFromFEDChannels << std::endl;
 
          auto myDetId = DetId(t_detid);
 
          if (subid == PixelSubdetector::PixelBarrel) {
            theROCMap.fillSelectedRocs(myDetId, badRocsFromFEDChannels, 1.);
          }  // if it's barrel
          else if (subid == PixelSubdetector::PixelEndcap) {
            theROCMap.fillSelectedRocs(myDetId, badRocsFromFEDChannels, 1.);
          }  // if it's endcap
          else {
            throw cms::Exception("LogicError") << "Unknown Pixel SubDet ID " << std::endl;
          }  // else nonsense
        }    // loop on the channels
      }      // loop on the scenarios
 
      gStyle->SetOptStat(0);
      //=========================
      TCanvas canvas("Summary", "Summary", 1200, 1600);
      theROCMap.drawMaps(canvas);
 
      auto unpacked = SiPixelPI::unpack(std::get<0>(iov));
 
      for (unsigned int lay = 1; lay <= n_layers; lay++) {
        canvas.cd(lay);
        auto ltx = TLatex();
        ltx.SetTextFont(62);
        ltx.SetTextColor(kBlue);
        ltx.SetTextSize(0.055);
        ltx.SetTextAlign(11);
        ltx.DrawLatexNDC(gPad->GetLeftMargin(),
                         1 - gPad->GetTopMargin() + 0.01,
                         unpacked.first == 0
                             ? ("IOV:" + std::to_string(unpacked.second)).c_str()
                             : (std::to_string(unpacked.first) + "," + std::to_string(unpacked.second)).c_str());
      }
 
      for (unsigned int ring = 1; ring <= n_rings; ring++) {
        canvas.cd(n_layers + ring);
        auto ltx = TLatex();
        ltx.SetTextFont(62);
        ltx.SetTextColor(kBlue);
        ltx.SetTextSize(0.050);
        ltx.SetTextAlign(11);
        ltx.DrawLatexNDC(gPad->GetLeftMargin(),
                         1 - gPad->GetTopMargin() + 0.01,
                         unpacked.first == 0
                             ? ("IOV:" + std::to_string(unpacked.second)).c_str()
                             : (std::to_string(unpacked.first) + "," + std::to_string(unpacked.second)).c_str());
      }
 
      std::string fileName(m_imageFileName);
      canvas.SaveAs(fileName.c_str());
 
      // close the DB session
      condDbSession.transaction().commit();
 
      return true;
    }
 
  public:
    inline unsigned int closest_from_above(std::vector<unsigned int> const& vec, unsigned int value) {
      auto const it = std::lower_bound(vec.begin(), vec.end(), value);
      return vec.at(it - vec.begin() - 1);
    }
 
    inline unsigned int closest_from_below(std::vector<unsigned int> const& vec, unsigned int value) {
      auto const it = std::upper_bound(vec.begin(), vec.end(), value);
      return vec.at(it - vec.begin() - 1);
    }
 
    // auxilliary check
    struct compareKeys {
      std::string key;
      compareKeys(std::string const& i) : key(i) {}
 
      bool operator()(std::string const& i) { return (key == i); }
    };
 
  private:
    // tough luck, we can only do phase-1...
    static constexpr int numColumns = 416;
    static constexpr int numRows = 160;
    static constexpr int n_rings = 2;
    static constexpr int n_layers = 4;
 
    TrackerTopology m_trackerTopo;
    edm::ParameterSet m_connectionPset;
    cond::persistency::ConnectionPool m_connectionPool;
    std::string m_CablingTagName;
    std::string m_condDbCabling;
  };
 
  /************************************************
  1d histogram of SiPixelFEDChannelContainer of 1 IOV
  *************************************************/
 
  class SiPixelFEDChannelContainerScenarios : public PlotImage<SiPixelFEDChannelContainer, SINGLE_IOV> {
  public:
    SiPixelFEDChannelContainerScenarios()
        : PlotImage<SiPixelFEDChannelContainer, SINGLE_IOV>("SiPixelFEDChannelContainer scenarios count") {}
    bool fill() override {
      auto tag = PlotBase::getTag<0>();
      auto tagname = tag.name;
      auto iov = tag.iovs.front();
 
      TGaxis::SetMaxDigits(3);
 
      std::shared_ptr<SiPixelFEDChannelContainer> payload = fetchPayload(std::get<1>(iov));
      std::vector<std::string> scenarios = payload->getScenarioList();
      sort(scenarios.begin(), scenarios.end());
 
      TCanvas canvas("Canv", "Canv", 1200, 1000);
      canvas.cd();
      canvas.SetGrid();
      auto h1 = std::make_unique<TH1F>("Count",
                                       "SiPixelFEDChannelContainer Bad Roc count;Scenario index;n. of bad ROCs",
                                       scenarios.size(),
                                       1,
                                       scenarios.size());
      h1->SetStats(false);
 
      canvas.SetTopMargin(0.06);
      canvas.SetBottomMargin(0.12);
      canvas.SetLeftMargin(0.12);
      canvas.SetRightMargin(0.05);
      canvas.Modified();
 
      int scenarioIndex = 0;
      for (const auto& scenario : scenarios) {
        scenarioIndex++;
        int badRocCount = 0;
        LogDebug("SiPixelFEDChannelContainerScenarios") << scenario << std::endl;
        auto badChannelCollection = payload->getDetSetBadPixelFedChannels(scenario);
        for (const auto& disabledChannels : *badChannelCollection) {
          for (const auto& ch : disabledChannels) {
            int local_bad_rocs = ch.roc_last - ch.roc_first;
            badRocCount += local_bad_rocs;
          }  // loop on the channels
        }    // loop on the DetSetVector
 
        h1->SetBinContent(scenarioIndex, badRocCount);
      }  // loop on scenarios
 
      TGaxis::SetExponentOffset(-0.1, 0.01, "y");    // Y offset
      TGaxis::SetExponentOffset(-0.03, -0.10, "x");  // Y and Y offset for X axis
 
      h1->SetTitle("");
      h1->GetYaxis()->SetRangeUser(0., h1->GetMaximum() * 1.30);
      h1->SetFillColor(kRed);
      h1->SetMarkerStyle(20);
      h1->SetMarkerSize(1);
      h1->Draw("bar2");
 
      SiPixelPI::makeNicePlotStyle(h1.get());
 
      canvas.Update();
 
      TLegend legend = TLegend(0.30, 0.88, 0.95, 0.94);
      //legend.SetHeader(("#splitline{Payload hash: #bf{" + (std::get<1>(iov)) + "}}{Total Scenarios:"+std::to_string(scenarioIndex)+"}").c_str(),"C");  // option "C" allows to center the header
 
      legend.SetHeader(fmt::sprintf("Payload hash: #bf{%s}", std::get<1>(iov)).c_str(), "C");
      legend.AddEntry(h1.get(), fmt::sprintf("total scenarios: #bf{%s}", std::to_string(scenarioIndex)).c_str(), "F");
      legend.SetTextSize(0.025);
      legend.Draw("same");
 
      auto ltx = TLatex();
      ltx.SetTextFont(62);
      //ltx.SetTextColor(kBlue);
      //ltx.SetTextAlign(11);
      ltx.SetTextSize(0.040);
      ltx.DrawLatexNDC(
          gPad->GetLeftMargin(),
          1 - gPad->GetTopMargin() + 0.01,
          fmt::sprintf("#color[4]{%s} IOV: #color[4]{%s}", tagname, std::to_string(std::get<0>(iov))).c_str());
 
      std::string fileName(m_imageFileName);
      canvas.SaveAs(fileName.c_str());
 
      return true;
 
    }  // fill
  };
 
}  // namespace
 
// Register the classes as boost python plugin
PAYLOAD_INSPECTOR_MODULE(SiPixelFEDChannelContainer) {
  PAYLOAD_INSPECTOR_CLASS(SiPixelFEDChannelContainerTest);
  PAYLOAD_INSPECTOR_CLASS(SiPixelFEDChannelContainerScenarios);
}