ApvTimingHistograms Class Reference

#include <ApvTimingHistograms.h>

Inheritance diagram for ApvTimingHistograms:

Public Member Functions
	ApvTimingHistograms (const edm::ParameterSet &pset, DQMStore *)
void	histoAnalysis (bool debug) override
	~ApvTimingHistograms () override
Public Member Functions inherited from CommissioningHistograms
	CommissioningHistograms (const edm::ParameterSet &pset, DQMStore *const, const sistrip::RunType &)
	CommissioningHistograms ()
virtual void	configure (const edm::ParameterSet &, const edm::EventSetup &)
void	createCollations (const std::vector< std::string > &)
virtual void	createSummaryHisto (const sistrip::Monitorable &, const sistrip::Presentation &, const std::string &top_level_dir, const sistrip::Granularity &)
void	extractHistograms (const std::vector< std::string > &)
virtual void	printAnalyses ()
virtual void	printSummary ()
void	remove (std::string pattern="")
void	save (std::string &filename, uint32_t run_number=0, std::string partitionName="")
virtual	~CommissioningHistograms ()

Additional Inherited Members
Public Types inherited from CommissioningHistograms
typedef std::map< uint32_t, CommissioningAnalysis * >	Analyses
typedef Analyses::iterator	Analysis
typedef dqm::harvesting::DQMStore	DQMStore
typedef SummaryPlotFactory < CommissioningAnalysis * >	Factory
typedef std::map< uint32_t, uint32_t >	FedToFecMap
typedef std::vector< Histo * >	Histos
typedef std::map< uint32_t, Histos >	HistosMap
typedef dqm::harvesting::MonitorElement	MonitorElement
Static Public Member Functions inherited from CommissioningHistograms
static void	copyCustomInformation (DQMStore *const, const std::vector< std::string > &)
static uint32_t	runNumber (DQMStore *const, const std::vector< std::string > &)
static sistrip::RunType	runType (DQMStore *const, const std::vector< std::string > &)
Protected Member Functions inherited from CommissioningHistograms
DQMStore *const	bei () const
void	clearHistosMap ()
Analyses &	data (bool getMaskedData=false)
Factory *const	factory ()
TH1 *	histogram (const sistrip::Monitorable &, const sistrip::Presentation &, const sistrip::View &, const std::string &directory, const uint32_t &xbins, const float &xlow=1.*sistrip::invalid_, const float &xhigh=1.*sistrip::invalid_)
const HistosMap &	histos () const
const FedToFecMap &	mapping () const
void	printHistosMap ()
const edm::ParameterSet &	pset () const
const sistrip::RunType &	task () const
Protected Attributes inherited from CommissioningHistograms
std::unique_ptr< Factory >	factory_

Detailed Description

Definition at line 7 of file ApvTimingHistograms.h.

Constructor & Destructor Documentation

ApvTimingHistograms::ApvTimingHistograms	(	const edm::ParameterSet &	pset,
		DQMStore *	bei
	)

Definition at line 20 of file ApvTimingHistograms.cc.

References CommissioningHistograms::factory_, LogTrace, and sistrip::mlDqmClient_.

    : CommissioningHistograms(pset.getParameter<edm::ParameterSet>("ApvTimingParameters"), bei, sistrip::APV_TIMING) {
  factory_ = std::make_unique<ApvTimingSummaryFactory>();
  LogTrace(mlDqmClient_) << "[ApvTimingHistograms::" << __func__ << "]"
                         << " Constructing object...";
}

ApvTimingHistograms::~ApvTimingHistograms ( )

override

Definition at line 29 of file ApvTimingHistograms.cc.

References LogTrace, and sistrip::mlDqmClient_.

                                          {
  LogTrace(mlDqmClient_) << "[ApvTimingHistograms::" << __func__ << "]"
                         << " Destructing object...";
}

Member Function Documentation

void ApvTimingHistograms::histoAnalysis ( bool  debug )

overridevirtual

Reimplemented from CommissioningHistograms.

Definition at line 36 of file ApvTimingHistograms.cc.

References CommissioningAlgorithm::analysis(), SplitLinear::begin, SiStripFecKey::ccuAddr(), SiStripFecKey::ccuChan(), submitPVResolutionJobs::count, CommissioningHistograms::data(), relativeConstraints::empty, extract(), SiStripFecKey::fecCrate(), SiStripFecKey::fecRing(), SiStripFecKey::fecSlot(), ApvTimingAnalysis::foundTickMark(), CommissioningAnalysis::getErrorCodes(), CommissioningHistograms::histos(), cuy::ii, sistrip::invalid_, ApvTimingAnalysis::isValid(), SiStripFecKey::lldChan(), LogTrace, SiStripPI::max, SiStripPI::min, sistrip::mlDqmClient_, CommissioningHistograms::pset(), ApvTimingAnalysis::refTime(), contentValuesCheck::ss, ApvTimingAnalysis::time(), and sistrip::valid_.

                                                  {
  LogTrace(mlDqmClient_) << "[ApvTimingHistograms::" << __func__ << "]";

  // Some initialisation
  uint16_t valid = 0;
  HistosMap::const_iterator iter;
  Analyses::iterator ianal;
  std::map<std::string, uint16_t> errors;

  // Clear map holding analysis objects
  for (ianal = data().begin(); ianal != data().end(); ianal++) {
    if (ianal->second) {
      delete ianal->second;
    }
  }
  data().clear();

  // Reset minimum / maximum delays
  float time_min = 1. * sistrip::invalid_;
  float time_max = -1. * sistrip::invalid_;
  uint32_t device_min = sistrip::invalid_;
  uint32_t device_max = sistrip::invalid_;

  // Iterate through map containing histograms
  for (iter = histos().begin(); iter != histos().end(); iter++) {
    // Check vector of histos is not empty
    if (iter->second.empty()) {
      edm::LogWarning(mlDqmClient_) << "[ApvTimingHistograms::" << __func__ << "]"
                                    << " Zero histograms found!";
      continue;
    }

    // Retrieve pointers to histos
    std::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
    ApvTimingAnalysis* anal = new ApvTimingAnalysis(iter->first);
    ApvTimingAlgorithm algo(this->pset(), anal);
    algo.analysis(profs);
    data()[iter->first] = anal;

    // Check if tick mark found
    if (!anal->foundTickMark()) {
      continue;
    }

    // Find maximum time
    if (anal->time() > time_max) {
      time_max = anal->time();
      device_max = iter->first;
    }

    // Find minimum time
    if (anal->time() < time_min) {
      time_min = anal->time();
      device_min = iter->first;
    }
  }

  // Adjust maximum (and minimum) delay(s) to find optimum sampling point(s)
  if (time_max > sistrip::valid_ || time_max < -1. * sistrip::valid_) {
    edm::LogWarning(mlDqmClient_) << "[ApvTimingHistograms::" << __func__ << "]"
                                  << " Unable to set maximum time! Found unexpected value: " << time_max;

  } else {
    SiStripFecKey min(device_min);
    edm::LogVerbatim(mlDqmClient_) << "[ApvTimingHistograms::" << __func__ << "]"
                                   << " Crate/FEC/Ring/CCU/module/channel: " << min.fecCrate() << "/" << min.fecSlot()
                                   << "/" << min.fecRing() << "/" << min.ccuAddr() << "/" << min.ccuChan() << "/"
                                   << min.lldChan() << " has minimum time for tick mark rising edge [ns]: " << time_min;

    SiStripFecKey max(device_max);
    edm::LogVerbatim(mlDqmClient_) << "[ApvTimingHistograms::" << __func__ << "]"
                                   << " Crate/FEC/Ring/CCU/module/channel: " << max.fecCrate() << "/" << max.fecSlot()
                                   << "/" << max.fecRing() << "/" << max.ccuAddr() << "/" << max.ccuChan() << "/"
                                   << max.lldChan() << " has maximum time for tick mark rising edge [ns]: " << time_max;

    edm::LogVerbatim(mlDqmClient_) << "[ApvTimingHistograms::" << __func__ << "]"
                                   << " Difference b/w minimum and maximum times"
                                   << " for tick mark rising edges [ns] is: " << (time_max - time_min);
  }

  // Set reference time for all analysis objects
  for (ianal = data().begin(); ianal != data().end(); ianal++) {
    ApvTimingAnalysis* anal = dynamic_cast<ApvTimingAnalysis*>(ianal->second);
    if (!anal) {
      continue;
    }
    anal->refTime(time_max, this->pset().getParameter<int>("TargetDelay"));
    if (anal->isValid()) {
      valid++;
    }
    if (!anal->getErrorCodes().empty()) {
      errors[anal->getErrorCodes()[0]]++;
    }
  }

  if (!histos().empty()) {
    edm::LogVerbatim(mlDqmClient_) << "[ApvTimingHistograms::" << __func__ << "]"
                                   << " Analyzed histograms for " << histos().size() << " FED channels, of which "
                                   << valid << " (" << 100 * valid / histos().size() << "%) are valid.";
  } else {
    edm::LogWarning(mlDqmClient_) << "[ApvTimingHistograms::" << __func__ << "]"
                                  << " No histograms to analyze!";
  }

  if (!histos().empty()) {
    edm::LogVerbatim(mlDqmClient_) << "[ApvTimingHistograms::" << __func__ << "]"
                                   << " Analyzed histograms for " << histos().size() << " FED channels, of which "
                                   << valid << " (" << 100 * valid / histos().size() << "%) are valid.";
    if (!errors.empty()) {
      uint16_t count = 0;
      std::stringstream ss;
      ss << std::endl;
      std::map<std::string, uint16_t>::const_iterator ii;
      for (ii = errors.begin(); ii != errors.end(); ++ii) {
        ss << " " << ii->first << ": " << ii->second << std::endl;
        count += ii->second;
      }
      edm::LogWarning(mlDqmClient_) << "[ApvTimingHistograms::" << __func__ << "]"
                                    << " Found " << count << " errors (" << 100 * count / histos().size()
                                    << "%): " << ss.str();
    }
  } else {
    edm::LogWarning(mlDqmClient_) << "[ApvTimingHistograms::" << __func__ << "]"
                                  << " No histograms to analyze!";
  }
}