CMHistograms.cc

Go to the documentation of this file.

#include "DataFormats/FEDRawData/interface/FEDNumbering.h"
#include "DataFormats/SiStripCommon/interface/SiStripFedKey.h"
#include "DataFormats/SiStripCommon/interface/ConstantsForHardwareSystems.h"

#include "DQM/SiStripMonitorHardware/interface/CMHistograms.hh"

CMHistograms::CMHistograms() {
  for (unsigned int i(0); i < 500; i++) {
    doFed_[i] = false;
  }
}

CMHistograms::~CMHistograms() {}

void CMHistograms::initialise(const edm::ParameterSet& iConfig, std::ostringstream* pDebugStream) {
  getConfigForHistogram(medianAPV1vsAPV0_, "MedianAPV1vsAPV0", iConfig, pDebugStream);
  getConfigForHistogram(medianAPV0minusAPV1_, "MedianAPV0minusAPV1", iConfig, pDebugStream);

  getConfigForHistogram(meanCMPerFedvsFedId_, "MeanCMPerFedvsFedId", iConfig, pDebugStream);
  getConfigForHistogram(meanCMPerFedvsTime_, "MeanCMPerFedvsTime", iConfig, pDebugStream);
  getConfigForHistogram(variationsPerFedvsFedId_, "VariationsPerFedvsFedId", iConfig, pDebugStream);
  getConfigForHistogram(variationsPerFedvsTime_, "VariationsPerFedvsTime", iConfig, pDebugStream);

  getConfigForHistogram(medianAPV1vsAPV0perFED_, "MedianAPV1vsAPV0perFED", iConfig, pDebugStream);
  getConfigForHistogram(medianAPV0minusAPV1perFED_, "MedianAPV0minusAPV1perFED", iConfig, pDebugStream);

  getConfigForHistogram(medianperChannel_, "MedianperChannel", iConfig, pDebugStream);
  getConfigForHistogram(medianAPV0minusAPV1perChannel_, "MedianAPV0minusAPV1perChannel", iConfig, pDebugStream);

  getConfigForHistogram(tkMapConfig_, "TkHistoMap", iConfig, pDebugStream);

  if (iConfig.exists("FedIdVec")) {
    std::vector<unsigned int> lIdVec = iConfig.getUntrackedParameter<std::vector<unsigned int> >("FedIdVec");
    for (unsigned int i(0); i < lIdVec.size(); i++) {
      if (lIdVec.at(i) < 500)
        doFed_[lIdVec.at(i)] = true;
    }
  }
}

void CMHistograms::fillHistograms(const std::vector<CMvalues>& aVec, float aTime, unsigned int aFedId) {
  float lMean = 0;
  float lPrevMean = 0;

  for (unsigned iEle(0); iEle < aVec.size(); iEle++) {  //loop on elements

    CMvalues lVal = aVec[iEle];

    if (lVal.Medians.first >= 1024 || lVal.Medians.second >= 1024) {
      std::cout << "----- WARNING ! New max found: " << lVal.Medians.first << " " << lVal.Medians.second << " "
                << __FILE__ << " " << __LINE__ << std::endl;
    }

    fillHistogram(medianAPV1vsAPV0_, lVal.Medians.first, lVal.Medians.second);
    fillHistogram(medianAPV0minusAPV1_, lVal.Medians.first - lVal.Medians.second);

    lMean += lVal.Medians.first + lVal.Medians.second;
    lPrevMean += lVal.PreviousMedians.first + lVal.PreviousMedians.second;

    if (doFed_[aFedId]) {
      fillHistogram(medianAPV1vsAPV0perFEDMap_[aFedId], lVal.Medians.first, lVal.Medians.second);
      fillHistogram(medianAPV0minusAPV1perFEDMap_[aFedId], lVal.Medians.first - lVal.Medians.second);
      fillHistogram(medianperChannelMap_[aFedId][lVal.ChannelID], lVal.Medians.first);
      fillHistogram(medianAPV0minusAPV1perChannelMap_[aFedId][lVal.ChannelID],
                    lVal.Medians.first - lVal.Medians.second);
    }

  }  //loop on elements

  if (!aVec.empty()) {
    lMean = lMean / (2 * aVec.size());
    lPrevMean = lPrevMean / (2 * aVec.size());
  }

  fillHistogram(meanCMPerFedvsFedId_, aFedId, lMean);
  fillHistogram(meanCMPerFedvsTime_, aTime, lMean);
  fillHistogram(variationsPerFedvsFedId_, aFedId, lMean - lPrevMean);
  fillHistogram(variationsPerFedvsTime_, aTime, lMean - lPrevMean);
}

void CMHistograms::bookTopLevelHistograms(DQMStore::IBooker& ibooker, const TkDetMap* tkDetMap) {
  //book FED level histograms
  //get FED IDs
  //const unsigned int siStripFedIdMin = FEDNumbering::MINSiStripFEDID;
  const unsigned int siStripFedIdMax = FEDNumbering::MAXSiStripFEDID;
  histosBooked_.resize(siStripFedIdMax + 1, false);

  std::string lDir = ibooker.pwd() + "/";

  ibooker.cd(lDir);

  ibooker.setCurrentFolder(lDir);

  book2DHistogram(ibooker,
                  medianAPV1vsAPV0_,
                  "MedianAPV1vsAPV0",
                  "median APV1 vs APV0",
                  250,
                  0,
                  1024,
                  250,
                  0,
                  1024,
                  "median APV0",
                  "median APV1");

  bookHistogram(ibooker,
                medianAPV0minusAPV1_,
                "MedianAPV0minusAPV1",
                "median APV0 - median APV1",
                500,
                -500,
                500,
                "median APV0 - median APV1");

  bookProfile(ibooker,
              meanCMPerFedvsFedId_,
              "MeanCMPerFedvsFedId",
              "<CM> vs fedID",
              440,
              50,
              490,
              -1000,
              1000,
              "fedID",
              "<CM>^{FED}");

  bookProfile(ibooker, meanCMPerFedvsTime_, "MeanCMPerFedvsTime", "<CM> vs time", 0, 1000, "Time", "<CM>^{FED}");

  bookProfile(ibooker,
              variationsPerFedvsFedId_,
              "VariationsPerFedvsFedId",
              "<CM> vs fedID",
              440,
              50,
              490,
              -1000,
              1000,
              "fedID",
              "<CM>^{FED}_{t}-<CM>^{FED}_{t-1}");

  bookProfile(ibooker,
              variationsPerFedvsTime_,
              "VariationsPerFedvsTime",
              "<CM> vs time",
              0,
              1000,
              "Time",
              "<CM>^{FED}_{t}-<CM>^{FED}_{t-1}");

  ibooker.cd(lDir);

  //book map after, as it creates a new folder...
  if (tkMapConfig_.enabled) {
    tkmapCM_[0] = std::make_unique<TkHistoMap>(tkDetMap, "SiStrip/TkHisto", "TkHMap_MeanCMAPV", 0., true);
    tkmapCM_[1] = std::make_unique<TkHistoMap>(tkDetMap, "SiStrip/TkHisto", "TkHMap_RmsCMAPV", 0., true);
    tkmapCM_[2] = std::make_unique<TkHistoMap>(tkDetMap, "SiStrip/TkHisto", "TkHMap_MeanCMAPV0minusAPV1", -500., true);
    tkmapCM_[3] = std::make_unique<TkHistoMap>(tkDetMap, "SiStrip/TkHisto", "TkHMap_RmsCMAPV0minusAPV1", -500., true);
  } else {
    tkmapCM_[0] = nullptr;
    tkmapCM_[1] = nullptr;
    tkmapCM_[2] = nullptr;
    tkmapCM_[3] = nullptr;
  }

  for (unsigned int i = sistrip::FED_ID_MIN; i <= sistrip::FED_ID_MAX; i++)
    if (doFed_[i])
      bookFEDHistograms(ibooker, i);
}

void CMHistograms::bookFEDHistograms(DQMStore::IBooker& ibooker, unsigned int fedId) {
  if (!histosBooked_[fedId]) {
    //will do that only once

    SiStripFedKey fedKey(fedId, 0, 0, 0);
    std::stringstream fedIdStream;
    fedIdStream << fedId;

    ibooker.setCurrentFolder(fedKey.path());

    book2DHistogram(ibooker,
                    medianAPV1vsAPV0perFED_,
                    medianAPV1vsAPV0perFEDMap_[fedId],
                    "MedianAPV1vsAPV0forFED" + fedIdStream.str(),
                    "median APV1 vs APV0 for FED " + fedIdStream.str(),
                    250,
                    0,
                    500,
                    250,
                    0,
                    500,
                    "APV0",
                    "APV1");

    bookHistogram(ibooker,
                  medianAPV0minusAPV1perFED_,
                  medianAPV0minusAPV1perFEDMap_[fedId],
                  "MedianAPV0minusAPV1forFED" + fedIdStream.str(),
                  "median APV0 - median APV1 for FED " + fedIdStream.str(),
                  500,
                  -500,
                  500,
                  "#Delta(medians)");

    bookChannelsHistograms(ibooker, fedId);

    histosBooked_[fedId] = true;
  }
}

void CMHistograms::bookChannelsHistograms(DQMStore::IBooker& ibooker, unsigned int fedId) {
  SiStripFedKey fedKey(fedId, 0, 0, 0);
  std::stringstream fedIdStream;
  fedIdStream << fedId;

  ibooker.setCurrentFolder(fedKey.path());
  medianperChannelMap_[fedId].resize(sistrip::FEDCH_PER_FED, nullptr);
  medianAPV0minusAPV1perChannelMap_[fedId].resize(sistrip::FEDCH_PER_FED, nullptr);

  for (unsigned int iCh(0); iCh < sistrip::FEDCH_PER_FED; iCh++) {
    std::ostringstream lName0, lTitle0, lName1, lTitle1, lName2, lTitle2;
    lName0 << "MedianForFed" << fedId << "Channel" << iCh;
    lTitle0 << "Median for FED/Ch " << fedId << "/" << iCh;
    lName2 << "MedianAPV0minusAPV1ForFed" << fedId << "Channel" << iCh;
    lTitle2 << "Median APV0-APV1 for FED/Ch " << fedId << "/" << iCh;

    bookHistogram(ibooker,
                  medianperChannel_,
                  medianperChannelMap_[fedId][iCh],
                  lName0.str(),
                  lTitle0.str(),
                  250,
                  0,
                  500,
                  "median APVs");

    bookHistogram(ibooker,
                  medianAPV0minusAPV1perChannel_,
                  medianAPV0minusAPV1perChannelMap_[fedId][iCh],
                  lName2.str(),
                  lTitle2.str(),
                  250,
                  -500,
                  500,
                  "median APV0-APV1");
  }
}

void CMHistograms::bookAllFEDHistograms(DQMStore::IBooker& ibooker) {
  //get FED IDs
  const unsigned int siStripFedIdMin = FEDNumbering::MINSiStripFEDID;
  const unsigned int siStripFedIdMax = FEDNumbering::MAXSiStripFEDID;
  //book them
  for (unsigned int iFed = siStripFedIdMin; iFed <= siStripFedIdMax; iFed++) {
    bookFEDHistograms(ibooker, iFed);
  }
}

bool CMHistograms::tkHistoMapEnabled(unsigned int aIndex) { return tkMapConfig_.enabled; }

TkHistoMap* CMHistograms::tkHistoMapPointer(unsigned int aIndex) {
  assert(aIndex < 4);
  return tkmapCM_[aIndex].get();
}