SiStripNoisesDQM.cc

Go to the documentation of this file.

#include "DQM/SiStripMonitorSummary/interface/SiStripNoisesDQM.h"

#include "DQMServices/Core/interface/DQMStore.h"

#include "TCanvas.h"

// -----
SiStripNoisesDQM::SiStripNoisesDQM(const edm::EventSetup &eSetup,
                                   edm::RunNumber_t iRun,
                                   edm::ParameterSet const &hPSet,
                                   edm::ParameterSet const &fPSet)
    : SiStripBaseCondObjDQM(eSetup, iRun, hPSet, fPSet) {
  gainRenormalisation_ = hPSet_.getParameter<bool>("GainRenormalisation");
  simGainRenormalisation_ = hPSet_.getParameter<bool>("SimGainRenormalisation");
  if (gainRenormalisation_ && !simGainRenormalisation_) {
    eSetup.get<SiStripApvGainRcd>().get(gainHandle_);
  }
  if (simGainRenormalisation_) {
    eSetup.get<SiStripApvGainSimRcd>().get(gainHandle_);
  }

  // Build the Histo_TkMap:
  if (HistoMaps_On_) {
    edm::ESHandle<TkDetMap> tkDetMapHandle;
    eSetup.get<TrackerTopologyRcd>().get(tkDetMapHandle);
    Tk_HM_ = std::make_unique<TkHistoMap>(tkDetMapHandle.product(), "SiStrip/Histo_Map", "MeanNoise_TkMap", 0.);
  }
}
// -----

// -----
SiStripNoisesDQM::~SiStripNoisesDQM() {}
// -----

// -----
void SiStripNoisesDQM::getActiveDetIds(const edm::EventSetup &eSetup) {
  getConditionObject(eSetup);
  noiseHandle_->getDetIds(activeDetIds);
}

// -----
void SiStripNoisesDQM::fillMEsForDet(const ModMEs &_selModME_, uint32_t selDetId_, const TrackerTopology *tTopo) {
  ModMEs selModME_ = _selModME_;
  std::vector<uint32_t> DetIds;
  noiseHandle_->getDetIds(DetIds);

  SiStripNoises::Range noiseRange = noiseHandle_->getRange(selDetId_);

  int nStrip = reader->getNumberOfApvsAndStripLength(selDetId_).first * 128;

  getModMEs(selModME_, selDetId_, tTopo);

  float gainFactor;
  float stripnoise;

  SiStripApvGain::Range gainRange;
  if (gainRenormalisation_ || simGainRenormalisation_) {
    gainRange = gainHandle_->getRange(selDetId_);
  }

  for (int istrip = 0; istrip < nStrip; ++istrip) {
    if (gainRenormalisation_ || simGainRenormalisation_)
      gainFactor = gainHandle_->getStripGain(istrip, gainRange) ? gainHandle_->getStripGain(istrip, gainRange) : 1.;
    else
      gainFactor = 1;

    stripnoise = noiseHandle_->getNoise(istrip, noiseRange) / gainFactor;
    if (CondObj_fillId_ == "onlyProfile" || CondObj_fillId_ == "ProfileAndCumul") {
      selModME_.ProfileDistr->Fill(istrip + 1, stripnoise);
    }
    if (CondObj_fillId_ == "onlyCumul" || CondObj_fillId_ == "ProfileAndCumul") {
      selModME_.CumulDistr->Fill(stripnoise);
    }
  }  // istrip
}

// -----
// FIXME too long. factorize this method.
// FIXME the number of lines of code in the derived classes should be reduced
// ONLY at what cannot be done in the base class because of the specific
// implementation
// FIXME of the derived class. Moreover, several loops on the same quantities
// should be avoided...

void SiStripNoisesDQM::fillMEsForLayer(
    /*std::map<uint32_t, ModMEs> selMEsMap_,*/ uint32_t selDetId_, const TrackerTopology *tTopo) {
  // ----
  int subdetectorId_ = ((selDetId_ >> 25) & 0x7);

  if (subdetectorId_ < 3 || subdetectorId_ > 6) {
    edm::LogError("SiStripNoisesDQM") << "[SiStripNoisesDQM::fillMEsForLayer] WRONG INPUT : no such "
                                         "subdetector type : "
                                      << subdetectorId_ << " no folder set!" << std::endl;
    return;
  }
  // ----

  std::map<uint32_t, ModMEs>::iterator selMEsMapIter_ = SummaryMEsMap_.find(getLayerNameAndId(selDetId_, tTopo).second);
  ModMEs selME_;
  if (selMEsMapIter_ != SummaryMEsMap_.end())
    selME_ = selMEsMapIter_->second;
  getSummaryMEs(selME_, selDetId_, tTopo);

  SiStripNoises::Range noiseRange = noiseHandle_->getRange(selDetId_);
  int nStrip = reader->getNumberOfApvsAndStripLength(selDetId_).first * 128;

  float stripnoise = -1.;
  float meanNoise = 0;
  int Nbadstrips = 0;

  SiStripApvGain::Range gainRange;
  if (gainRenormalisation_ || simGainRenormalisation_) {
    gainRange = gainHandle_->getRange(selDetId_);
  }
  float gainFactor = 1;

  SiStripHistoId hidmanager;

  if (hPSet_.getParameter<bool>("FillSummaryProfileAtLayerLevel")) {
    // --> profile summary
    std::string hSummaryOfProfile_description;
    hSummaryOfProfile_description = hPSet_.getParameter<std::string>("SummaryOfProfile_description");

    std::string hSummaryOfProfile_name;
    hSummaryOfProfile_name = hidmanager.createHistoLayer(
        hSummaryOfProfile_description, "layer", getLayerNameAndId(selDetId_, tTopo).first, "");
  }
  if (hPSet_.getParameter<bool>("FillCumulativeSummaryAtLayerLevel")) {
    std::string hSummaryOfCumul_description;
    hSummaryOfCumul_description = hPSet_.getParameter<std::string>("Cumul_description");

    std::string hSummaryOfCumul_name;
    hSummaryOfCumul_name = hidmanager.createHistoLayer(
        hSummaryOfCumul_description, "layer", getStringNameAndId(selDetId_, tTopo).first, "");
  }
  if (hPSet_.getParameter<bool>("FillSummaryAtLayerLevel")) {
    // --> cumul summary
    std::string hSummary_description;
    hSummary_description = hPSet_.getParameter<std::string>("Summary_description");

    std::string hSummary_name;
    hSummary_name =
        hidmanager.createHistoLayer(hSummary_description, "layer", getLayerNameAndId(selDetId_, tTopo).first, "");
  }

  for (int istrip = 0; istrip < nStrip; ++istrip) {
    if (gainRenormalisation_ || simGainRenormalisation_) {
      gainFactor = gainHandle_->getStripGain(istrip, gainRange) ? gainHandle_->getStripGain(istrip, gainRange) : 1.;
    } else {
      gainFactor = 1.;
    }

    stripnoise = noiseHandle_->getNoise(istrip, noiseRange) / gainFactor;
    meanNoise += stripnoise;
    if (hPSet_.getParameter<bool>("FillSummaryProfileAtLayerLevel")) {
      if (CondObj_fillId_ == "onlyProfile" || CondObj_fillId_ == "ProfileAndCumul") {
        selME_.SummaryOfProfileDistr->Fill(istrip + 1, stripnoise);
      }
    }

    if (hPSet_.getParameter<bool>("FillCumulativeSummaryAtLayerLevel")) {
      if (CondObj_fillId_ == "onlyCumul" || CondObj_fillId_ == "ProfileAndCumul") {
        selME_.SummaryOfCumulDistr->Fill(stripnoise);
      }
    }

    // Fill the TkMap
    if (fPSet_.getParameter<bool>("TkMap_On") || hPSet_.getParameter<bool>("TkMap_On")) {
      fillTkMap(selDetId_, stripnoise);
    }

  }  // istrip

  if (hPSet_.getParameter<bool>("FillSummaryAtLayerLevel")) {
    meanNoise = meanNoise / (nStrip - Nbadstrips);
    // get detIds belonging to same layer to fill X-axis with detId-number

    std::vector<uint32_t> sameLayerDetIds_;
    sameLayerDetIds_ = GetSameLayerDetId(activeDetIds, selDetId_, tTopo);

    std::vector<uint32_t>::const_iterator ibound =
        lower_bound(sameLayerDetIds_.begin(), sameLayerDetIds_.end(), selDetId_);
    if (ibound != sameLayerDetIds_.end() && *ibound == selDetId_)
      selME_.SummaryDistr->Fill(ibound - sameLayerDetIds_.begin() + 1, meanNoise);

    // Fill the Histo_TkMap with the mean Noise:
    if (HistoMaps_On_) {
      Tk_HM_->fill(selDetId_, meanNoise);
    }

    // Check the axis range for tkmap, and in case redefine;
    int intNoise = int(meanNoise);
    if (intNoise + 1 > (int)tkMapScaler.size()) {
      tkMapScaler.resize(intNoise + 1, 0);
    }
    tkMapScaler[intNoise]++;
  }
}