SamplingHistograms.cc

Go to the documentation of this file.

#include "DQM/SiStripCommissioningClients/interface/SamplingHistograms.h"
#include "CondFormats/SiStripObjects/interface/SamplingAnalysis.h"
#include "DataFormats/SiStripCommon/interface/SiStripConstants.h"
#include "DataFormats/SiStripCommon/interface/SiStripEnumsAndStrings.h"
#include "DQM/SiStripCommissioningAnalysis/interface/SamplingAlgorithm.h"
#include "DQM/SiStripCommissioningSummary/interface/SummaryGenerator.h"
#include "DQM/SiStripCommon/interface/ExtractTObject.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <iostream>
#include <sstream>
#include <iomanip>
#include "TProfile.h"

using namespace std;
using namespace sistrip;

// -----------------------------------------------------------------------------
SamplingHistograms::SamplingHistograms(const edm::ParameterSet& pset, DQMStore* bei, const sistrip::RunType& task)
    : CommissioningHistograms(pset.getParameter<edm::ParameterSet>("SamplingParameters"), bei, task), sOnCut_(3) {
  LogTrace(mlDqmClient_) << "[SamplingHistograms::" << __func__ << "]"
                         << " Constructing object...";
  factory_ = unique_ptr<SamplingSummaryFactory>(new SamplingSummaryFactory);
  // retreive the latency code from the root file
  std::string dataPath = std::string(sistrip::collate_) + "/" + sistrip::root_ + "/latencyCode";
  MonitorElement* codeElement = bei->get(dataPath);
  if (codeElement)
    latencyCode_ = codeElement->getIntValue();
  else
    latencyCode_ = 0;
}

// -----------------------------------------------------------------------------
SamplingHistograms::~SamplingHistograms() {
  LogTrace(mlDqmClient_) << "[SamplingHistograms::" << __func__ << "]"
                         << " Deleting object...";
}

// -----------------------------------------------------------------------------
void SamplingHistograms::histoAnalysis(bool debug) {
  // Clear map holding analysis objects
  Analyses::iterator ianal;
  for (ianal = data().begin(); ianal != data().end(); ianal++) {
    if (ianal->second) {
      delete ianal->second;
    }
  }
  data().clear();

  // Iterate through map containing vectors of profile histograms
  HistosMap::const_iterator iter = histos().begin();
  for (; iter != histos().end(); iter++) {
    // Check vector of histos is not empty (should be 1 histo)
    if (iter->second.empty()) {
      edm::LogWarning(mlDqmClient_) << "[SamplingHistograms::" << __func__ << "]"
                                    << " Zero collation histograms found!";
      continue;
    }

    // Retrieve pointers to profile histos for this FED channel
    vector<TH1*> profs;
    Histos::const_iterator ihis = iter->second.begin();
    for (; ihis != iter->second.end(); ihis++) {
      TProfile* prof = ExtractTObject<TProfile>().extract((*ihis)->me_);
      if (prof) {
        profs.push_back(prof);
      }
    }

    // Perform histo analysis
    SamplingAnalysis* anal = new SamplingAnalysis(iter->first);
    anal->setSoNcut(sOnCut_);
    SamplingAlgorithm algo(this->pset(), anal, latencyCode_);
    algo.analysis(profs);
    data()[iter->first] = anal;
  }
}

void SamplingHistograms::configure(const edm::ParameterSet& pset, const edm::EventSetup&) {
  //TODO: should use the parameter set. Why is this crashing ???
  //  sOnCut_ = pset.getParameter<double>("SignalToNoiseCut");
  sOnCut_ = 3.;
}