PedestalsTask Class Reference

#include <PedestalsTask.h>

Inheritance diagram for PedestalsTask:

Public Member Functions
	PedestalsTask (DQMStore *, const FedChannelConnection &)
	~PedestalsTask () override
Public Member Functions inherited from CommissioningTask
void	bookHistograms ()
	CommissioningTask (DQMStore *, const FedChannelConnection &, const std::string &my_name)
void	eventSetup (const edm::EventSetup *)
const uint32_t &	fillCntr () const
void	fillHistograms (const SiStripEventSummary &, const edm::DetSet< SiStripRawDigi > &)
void	fillHistograms (const SiStripEventSummary &, const edm::DetSet< SiStripRawDigi > &, const edm::DetSet< SiStripRawDigi > &)
void	fillHistograms (const SiStripEventSummary &, const edm::DetSet< SiStripRawDigi > &, const edm::DetSet< SiStripRawDigi > &, const std::vector< uint16_t > &)
void	fillHistograms (const SiStripEventSummary &, const uint16_t &fed_id, const std::map< uint16_t, float > &fed_ch)
const std::string &	myName () const
const uint32_t &	updateFreq () const
void	updateFreq (const uint32_t &)
void	updateHistograms ()
virtual	~CommissioningTask ()

Private Member Functions
void	book () override
void	fill (const SiStripEventSummary &, const edm::DetSet< SiStripRawDigi > &) override
void	update () override

Private Attributes
std::vector< HistoSet >	cm_
std::vector< HistoSet >	peds_

Additional Inherited Members
Public Types inherited from CommissioningTask
typedef dqm::legacy::DQMStore	DQMStore
typedef dqm::legacy::MonitorElement	MonitorElement
Protected Member Functions inherited from CommissioningTask
const FedChannelConnection &	connection () const
DQMStore *const	dqm () const
const edm::EventSetup *const	eventSetup () const
const uint32_t &	fecKey () const
const uint32_t &	fedKey () const
void	updateHistoSet (HistoSet &, const uint32_t &bin, const float &value)
void	updateHistoSet (CompactHistoSet &, const uint32_t &bin, const short &value)
void	updateHistoSet (HistoSet &, const uint32_t &bin)
void	updateHistoSet (CompactHistoSet &, const uint32_t &bin)
void	updateHistoSet (HistoSet &, const float &value)
void	updateHistoSet (CompactHistoSet &)
void	updateHistoSet (HistoSet &)

Detailed Description

Definition at line 10 of file PedestalsTask.h.

Constructor & Destructor Documentation

PedestalsTask::PedestalsTask	(	DQMStore *	dqm,
		const FedChannelConnection &	conn
	)

Definition at line 15 of file PedestalsTask.cc.

References LogTrace, and sistrip::mlDqmSource_.

    : CommissioningTask(dqm, conn, "PedestalsTask"), peds_(), cm_() {
  LogTrace(mlDqmSource_) << "[PedestalsTask::" << __func__ << "]"
                         << " Constructing object...";
}

PedestalsTask::~PedestalsTask ( )

override

Definition at line 23 of file PedestalsTask.cc.

References LogTrace, and sistrip::mlDqmSource_.

                              {
  LogTrace(mlDqmSource_) << "[PedestalsTask::" << __func__ << "]"
                         << " Destructing object...";
}

Member Function Documentation

void PedestalsTask::book ( )

overrideprivatevirtual

Reimplemented from CommissioningTask.

Definition at line 30 of file PedestalsTask.cc.

References sistrip::APV, cm_, sistrip::extrainfo::commonMode_, CommissioningTask::connection(), CommissioningTask::dqm(), sistrip::EXPERT_HISTO, sistrip::FED_KEY, CommissioningTask::fedKey(), sistrip::LLD_CHAN, LogTrace, sistrip::mlDqmSource_, hlt_dqm_clientPB-live_cfg::nbins, sistrip::extrainfo::noise_, sistrip::PEDESTALS, sistrip::extrainfo::pedestals_, peds_, AlCaHLTBitMon_QueryRunRegistry::string, SiStripHistoTitle::title(), and runGCPTkAlMap::title.

                         {
  LogTrace(mlDqmSource_) << "[PedestalsTask::" << __func__ << "]";

  uint16_t nbins;
  std::string title;
  std::string extra_info;
  peds_.resize(2);
  nbins = 256;

  // Pedestals histogram
  extra_info = sistrip::extrainfo::pedestals_;
  peds_[0].isProfile_ = true;

  title = SiStripHistoTitle(sistrip::EXPERT_HISTO,
                            sistrip::PEDESTALS,
                            sistrip::FED_KEY,
                            fedKey(),
                            sistrip::LLD_CHAN,
                            connection().lldChannel(),
                            extra_info)
              .title();

  peds_[0].histo(dqm()->bookProfile(title, title, nbins, -0.5, nbins * 1. - 0.5, 1025, 0., 1025.));

  peds_[0].vNumOfEntries_.resize(nbins, 0);
  peds_[0].vSumOfContents_.resize(nbins, 0);
  peds_[0].vSumOfSquares_.resize(nbins, 0);

  // Noise histogram
  extra_info = sistrip::extrainfo::noise_;
  peds_[1].isProfile_ = true;

  title = SiStripHistoTitle(sistrip::EXPERT_HISTO,
                            sistrip::PEDESTALS,
                            sistrip::FED_KEY,
                            fedKey(),
                            sistrip::LLD_CHAN,
                            connection().lldChannel(),
                            extra_info)
              .title();

  peds_[1].histo(dqm()->bookProfile(title, title, nbins, -0.5, nbins * 1. - 0.5, 1025, 0., 1025.));

  peds_[1].vNumOfEntries_.resize(nbins, 0);
  peds_[1].vSumOfContents_.resize(nbins, 0);
  peds_[1].vSumOfSquares_.resize(nbins, 0);

  // Common mode histograms
  cm_.resize(2);
  nbins = 1024;
  for (uint16_t iapv = 0; iapv < 2; iapv++) {
    title = SiStripHistoTitle(sistrip::EXPERT_HISTO,
                              sistrip::PEDESTALS,
                              sistrip::FED_KEY,
                              fedKey(),
                              sistrip::APV,
                              connection().i2cAddr(iapv),
                              sistrip::extrainfo::commonMode_)
                .title();

    cm_[iapv].histo(dqm()->book1D(title, title, nbins, -0.5, nbins * 1. - 0.5));
    cm_[iapv].isProfile_ = false;

    cm_[iapv].vNumOfEntries_.resize(nbins, 0);
    cm_[iapv].vNumOfEntries_.resize(nbins, 0);
  }
}

void PedestalsTask::fill ( const SiStripEventSummary &  summary, const edm::DetSet< SiStripRawDigi > &  digis  )

overrideprivatevirtual

Reimplemented from CommissioningTask.

Definition at line 100 of file PedestalsTask.cc.

References gpuClustering::adc, cm_, edm::DetSet< T >::data, change_name::diff, sistrip::mlDqmSource_, hlt_dqm_clientPB-live_cfg::nbins, peds_, and CommissioningTask::updateHistoSet().

                                                                                                   {
  if (digis.data.size() != peds_[0].vNumOfEntries_.size()) {
    edm::LogWarning(mlDqmSource_) << "[PedestalsTask::" << __func__ << "]"
                                  << " Unexpected number of digis: " << digis.data.size();
    return;
  }

  // Check number of digis
  uint16_t nbins = peds_[0].vNumOfEntries_.size();
  if (digis.data.size() < nbins) {
    nbins = digis.data.size();
  }

  //@@ Inefficient!!!
  uint16_t napvs = nbins / 128;
  std::vector<uint32_t> cm;
  cm.resize(napvs, 0);

  // Calc common mode for both APVs
  std::vector<uint16_t> adc;
  for (uint16_t iapv = 0; iapv < napvs; iapv++) {
    adc.clear();
    adc.reserve(128);
    for (uint16_t ibin = 0; ibin < 128; ibin++) {
      if ((iapv * 128) + ibin < nbins) {
        adc.push_back(digis.data[(iapv * 128) + ibin].adc());  //@@ VIRGIN RAW DATA (MUX, APV READOUT)
      }
    }
    sort(adc.begin(), adc.end());
    uint16_t index = adc.size() % 2 ? adc.size() / 2 : adc.size() / 2 - 1;
    if (!adc.empty()) {
      cm[iapv] = static_cast<uint32_t>(adc[index]);
    }
  }

  for (uint16_t ibin = 0; ibin < nbins; ibin++) {
    updateHistoSet(peds_[0], ibin, digis.data[ibin].adc());  // peds and raw noise
    float diff = static_cast<float>(digis.data[ibin].adc()) - static_cast<float>(cm[ibin / 128]);
    updateHistoSet(peds_[1], ibin, diff);  // residuals and real noise
  }

  if (cm.size() < cm_.size()) {
    edm::LogWarning(mlDqmSource_) << "[PedestalsTask::" << __func__ << "]"
                                  << " Fewer CM values than expected: " << cm.size();
  }

  updateHistoSet(cm_[0], cm[0]);
  updateHistoSet(cm_[1], cm[1]);
}

void PedestalsTask::update ( )

overrideprivatevirtual

Reimplemented from CommissioningTask.

Definition at line 152 of file PedestalsTask.cc.

References cm_, relativeConstraints::error, extract(), timingPdfMaker::histo, cuy::ii, SiStripPI::mean, gpuVertexFinder::noise, peds_, UpdateTProfile::setBinContent(), mathSSE::sqrt(), and CommissioningTask::updateHistoSet().

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), dqm-mbProfile.Profile::finish(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

                           {
  // Pedestals
  updateHistoSet(peds_[0]);

  // Noise (cannot use HistoSet directly, as want to plot noise as "contents", not "error")
  TProfile* histo = ExtractTObject<TProfile>().extract(peds_[1].histo());
  for (uint16_t ii = 0; ii < peds_[1].vNumOfEntries_.size(); ++ii) {
    float mean = 0.;
    float spread = 0.;
    float entries = peds_[1].vNumOfEntries_[ii];
    if (entries > 0.) {
      mean = peds_[1].vSumOfContents_[ii] / entries;
      spread = sqrt(peds_[1].vSumOfSquares_[ii] / entries - mean * mean);
    }

    float noise = spread;
    float error = 0;  // sqrt(entries) / entries;

    UpdateTProfile::setBinContent(histo, ii + 1, entries, noise, error);
  }

  // Common mode
  updateHistoSet(cm_[0]);
  updateHistoSet(cm_[1]);
}

Member Data Documentation

std::vector<HistoSet> PedestalsTask::cm_

private

Definition at line 21 of file PedestalsTask.h.

Referenced by book(), fill(), and update().

std::vector<HistoSet> PedestalsTask::peds_

private

Definition at line 20 of file PedestalsTask.h.

Referenced by book(), fill(), and update().