VpspScanAlgorithm.cc

Go to the documentation of this file.

#include "DQM/SiStripCommissioningAnalysis/interface/VpspScanAlgorithm.h"
#include "CondFormats/SiStripObjects/interface/VpspScanAnalysis.h"
#include "DataFormats/SiStripCommon/interface/SiStripHistoTitle.h"
#include "DataFormats/SiStripCommon/interface/SiStripEnumsAndStrings.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "TProfile.h"
#include "TH1.h"
#include <iostream>
#include <sstream>
#include <iomanip>
#include <cmath>

using namespace sistrip;

// -----------------------------------------------------------------------------
//
VpspScanAlgorithm::VpspScanAlgorithm(const edm::ParameterSet& pset, VpspScanAnalysis* const anal)
    : CommissioningAlgorithm(anal), histos_(2, Histo(nullptr, "")) {
  ;
}

// ----------------------------------------------------------------------------
//
void VpspScanAlgorithm::extract(const std::vector<TH1*>& histos) {
  if (!anal()) {
    edm::LogWarning(mlCommissioning_) << "[VpspScanAlgorithm::" << __func__ << "]"
                                      << " NULL pointer to Analysis object!";
    return;
  }

  // Check number of histograms
  if (histos.size() != 2) {
    anal()->addErrorCode(sistrip::numberOfHistos_);
  }

  // Extract FED key from histo title
  if (!histos.empty()) {
    anal()->fedKey(extractFedKey(histos.front()));
  }

  // Extract histograms
  std::vector<TH1*>::const_iterator ihis = histos.begin();
  for (; ihis != histos.end(); ihis++) {
    // Check pointer
    if (!(*ihis)) {
      continue;
    }

    // Check name
    SiStripHistoTitle title((*ihis)->GetName());
    if (title.runType() != sistrip::VPSP_SCAN) {
      anal()->addErrorCode(sistrip::unexpectedTask_);
      continue;
    }

    // Extract APV number
    uint16_t apv = sistrip::invalid_;
    if (title.extraInfo().find(sistrip::apv_) != std::string::npos) {
      std::stringstream ss;
      ss << title.extraInfo().substr(title.extraInfo().find(sistrip::apv_) + (sizeof(sistrip::apv_) - 1), 1);
      ss >> std::dec >> apv;
    }

    if (apv <= 1) {
      histos_[apv].first = *ihis;
      histos_[apv].second = (*ihis)->GetName();
    } else {
      anal()->addErrorCode(sistrip::unexpectedExtraInfo_);
    }
  }
}

// -----------------------------------------------------------------------------
//
void VpspScanAlgorithm::analyse() {
  if (!anal()) {
    edm::LogWarning(mlCommissioning_) << "[VpspScanAlgorithm::" << __func__ << "]"
                                      << " NULL pointer to base Analysis object!";
    return;
  }

  CommissioningAnalysis* tmp = const_cast<CommissioningAnalysis*>(anal());
  VpspScanAnalysis* anal = dynamic_cast<VpspScanAnalysis*>(tmp);
  if (!anal) {
    edm::LogWarning(mlCommissioning_) << "[VpspScanAlgorithm::" << __func__ << "]"
                                      << " NULL pointer to derived Analysis object!";
    return;
  }

  // from deprecated()

  std::vector<const TProfile*> histos;
  std::vector<uint16_t> monitorables;
  for (uint16_t iapv = 0; iapv < 2; iapv++) {
    monitorables.clear();
    monitorables.resize(7, sistrip::invalid_);

    histos.clear();
    histos.push_back(const_cast<const TProfile*>(dynamic_cast<TProfile*>(histos_[iapv].first)));

    if (!histos[0]) {
      anal->addErrorCode(sistrip::nullPtr_);
      continue;
    }

    // Find "top" plateau
    int first = sistrip::invalid_;
    float top = -1. * sistrip::invalid_;
    ;
    for (int k = 5; k < 55; k++) {
      if (!histos[0]->GetBinEntries(k)) {
        continue;
      }
      if (histos[0]->GetBinContent(k) >= top) {
        first = k;
        top = histos[0]->GetBinContent(k);
      }
    }
    if (top < -1. * sistrip::valid_) {
      top = sistrip::invalid_;
    }  //@@ just want +ve invalid number here
    if (top > 1. * sistrip::valid_) {
      anal->addErrorCode(sistrip::noTopPlateau_);
      continue;
    }
    monitorables[5] = static_cast<uint16_t>(top);
    monitorables[3] = first;

    // Find "bottom" plateau
    int last = sistrip::invalid_;
    float bottom = 1. * sistrip::invalid_;
    for (int k = 55; k > 5; k--) {
      if (!histos[0]->GetBinEntries(k)) {
        continue;
      }
      if (histos[0]->GetBinContent(k) <= bottom) {
        last = k;
        bottom = histos[0]->GetBinContent(k);
      }
    }
    if (bottom > 1. * sistrip::valid_) {
      anal->addErrorCode(sistrip::noBottomPlateau_);
      continue;
    }
    monitorables[6] = static_cast<uint16_t>(bottom);
    monitorables[4] = last;

    // Set optimum baseline level
    float opt = bottom + (top - bottom) * 1. / 3.;
    monitorables[1] = static_cast<uint16_t>(opt);

    // Find optimum VPSP setting
    uint16_t vpsp = sistrip::invalid_;
    if (opt < 1. * sistrip::valid_) {
      uint16_t ivpsp = 5;
      for (; ivpsp < 55; ivpsp++) {
        if (histos[0]->GetBinContent(ivpsp) < opt) {
          break;
        }
      }
      if (ivpsp != 54) {
        vpsp = ivpsp;
        monitorables[0] = vpsp;
      } else {
        anal->addErrorCode(sistrip::noVpspSetting_);
        continue;
      }

    } else {
      anal->addErrorCode(sistrip::noBaselineLevel_);
      continue;
    }

    // Set analysis results for both APVs
    anal->vpsp_[iapv] = monitorables[0];
    anal->adcLevel_[iapv] = monitorables[1];
    anal->fraction_[iapv] = monitorables[2];
    anal->topEdge_[iapv] = monitorables[3];
    anal->bottomEdge_[iapv] = monitorables[4];
    anal->topLevel_[iapv] = monitorables[5];
    anal->bottomLevel_[iapv] = monitorables[6];
  }
}