VpspScanAlgorithm Class Reference

Histogram-based analysis for VPSP scan. More...

#include <VpspScanAlgorithm.h>

Inheritance diagram for VpspScanAlgorithm:

Public Member Functions
const Histo &	histo (const uint16_t &apv) const
	VpspScanAlgorithm (const edm::ParameterSet &pset, VpspScanAnalysis *const)
	~VpspScanAlgorithm () override
Public Member Functions inherited from CommissioningAlgorithm
void	analysis (const std::vector< TH1 * > &)
	CommissioningAlgorithm ()
	CommissioningAlgorithm (CommissioningAnalysis *const)
virtual	~CommissioningAlgorithm ()

Private Member Functions
void	analyse () override
void	extract (const std::vector< TH1 * > &) override
	VpspScanAlgorithm ()

Private Attributes
std::vector< Histo >	histos_

Additional Inherited Members
Public Types inherited from CommissioningAlgorithm
typedef std::pair< TH1 *, std::string >	Histo
Protected Member Functions inherited from CommissioningAlgorithm
CommissioningAnalysis *const	anal () const
uint32_t	extractFedKey (const TH1 *const)

Detailed Description

Histogram-based analysis for VPSP scan.

Author

M. Wingham, R.Bainbridge

Definition at line 16 of file VpspScanAlgorithm.h.

Constructor & Destructor Documentation

VpspScanAlgorithm() [1/2]

VpspScanAlgorithm::VpspScanAlgorithm	(	const edm::ParameterSet &	pset,
		VpspScanAnalysis * const	anal
	)

Definition at line 17 of file VpspScanAlgorithm.cc.

    : CommissioningAlgorithm(anal), histos_(2, Histo(nullptr, "")) {
  ;
}

~VpspScanAlgorithm()

VpspScanAlgorithm::~VpspScanAlgorithm ( )

inlineoverride

Definition at line 20 of file VpspScanAlgorithm.h.

{ ; }

VpspScanAlgorithm() [2/2]

VpspScanAlgorithm::VpspScanAlgorithm ( )

inlineprivate

Definition at line 26 of file VpspScanAlgorithm.h.

{ ; }

Member Function Documentation

analyse()

void VpspScanAlgorithm::analyse ( )

overrideprivatevirtual

Performs histogram anaysis.

Implements CommissioningAlgorithm.

Definition at line 75 of file VpspScanAlgorithm.cc.

                                {
  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];
  }
}

References CommissioningAnalysis::addErrorCode(), CommissioningAlgorithm::anal(), dqmdumpme::first, combine::histos, histos_, sistrip::invalid_, dqmdumpme::k, dqmdumpme::last, sistrip::mlCommissioning_, sistrip::noBaselineLevel_, sistrip::noBottomPlateau_, sistrip::noTopPlateau_, sistrip::noVpspSetting_, sistrip::nullPtr_, runTheMatrix::opt, createJobs::tmp, and sistrip::valid_.

extract()

void VpspScanAlgorithm::extract ( const std::vector< TH1 * > &  histos )

overrideprivatevirtual

Extracts and organises histograms.

Implements CommissioningAlgorithm.

Definition at line 24 of file VpspScanAlgorithm.cc.

                                                             {
  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_);
    }
  }
}

References CommissioningAnalysis::addErrorCode(), CommissioningAlgorithm::anal(), sistrip::apv_, TauDecayModes::dec, CommissioningAlgorithm::extractFedKey(), CommissioningAnalysis::fedKey(), combine::histos, histos_, sistrip::invalid_, sistrip::mlCommissioning_, sistrip::numberOfHistos_, contentValuesCheck::ss, runGCPTkAlMap::title, sistrip::unexpectedExtraInfo_, sistrip::unexpectedTask_, and sistrip::VPSP_SCAN.

histo()

const Histo& VpspScanAlgorithm::histo

(

const uint16_t &

apv

)

const

Histogram pointer and title.

Member Data Documentation

histos_

std::vector<Histo> VpspScanAlgorithm::histos_

private

Pointers and titles for histograms.

Definition at line 36 of file VpspScanAlgorithm.h.

Referenced by analyse(), and extract().