Line Code
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57
#include "DQM/SiStripCommissioningSources/interface/FastFedCablingTask.h"
#include "DataFormats/SiStripCommon/interface/SiStripConstants.h"
#include "DataFormats/SiStripCommon/interface/SiStripHistoTitle.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <algorithm>
#include <sstream>
#include <iomanip>

using namespace sistrip;

// -----------------------------------------------------------------------------
//
FastFedCablingTask::FastFedCablingTask(DQMStore* dqm, const FedChannelConnection& conn)
    : CommissioningTask(dqm, conn, "FastFedCablingTask"), histo_() {}

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

// -----------------------------------------------------------------------------
//
void FastFedCablingTask::book() {
  std::string title = SiStripHistoTitle(sistrip::EXPERT_HISTO,
                                        sistrip::FAST_CABLING,
                                        sistrip::FED_KEY,
                                        fedKey(),
                                        sistrip::LLD_CHAN,
                                        connection().lldChannel())
                          .title();

  uint16_t nbins = 34;
  histo_.histo(dqm()->bookProfile(title, title, nbins, -0.5, nbins * 1. - 0.5, 1025, 0., 1025.));

  histo_.vNumOfEntries_.resize(nbins, 0);
  histo_.vSumOfContents_.resize(nbins, 0);
  histo_.vSumOfSquares_.resize(nbins, 0);
}

// -----------------------------------------------------------------------------
//
void FastFedCablingTask::fill(const SiStripEventSummary& summary, const edm::DetSet<SiStripRawDigi>& digis) {
  if (digis.data.empty()) {
    edm::LogWarning(mlDqmSource_) << "[FastFedCablingTask::" << __func__ << "]"
                                  << " No digis found!";
    return;
  }

  uint32_t bin = summary.binNumber();
  for (uint16_t ibin = 0; ibin < digis.data.size(); ibin++) {
    updateHistoSet(histo_, bin, digis.data[ibin].adc());
  }
}

// -----------------------------------------------------------------------------
//
void FastFedCablingTask::update() { updateHistoSet(histo_); }