FedTimingHistograms.cc

Go to the documentation of this file.

#include "DQM/SiStripCommissioningClients/interface/FedTimingHistograms.h"
#include "CondFormats/SiStripObjects/interface/FedTimingAnalysis.h"
#include "DataFormats/SiStripCommon/interface/SiStripEnumsAndStrings.h"
#include "DQM/SiStripCommissioningAnalysis/interface/FedTimingAlgorithm.h"
#include "DQM/SiStripCommissioningSummary/interface/SummaryGenerator.h"
#include <iostream>
#include <sstream>
#include <iomanip>

using namespace std;

// -----------------------------------------------------------------------------
FedTimingHistograms::FedTimingHistograms(const edm::ParameterSet& pset, DQMStore* bei)
    : CommissioningHistograms(pset.getParameter<edm::ParameterSet>("FedTimingParameters"), bei, sistrip::FED_TIMING),
      factory_(new Factory),
      optimumSamplingPoint_(15.),
      minDelay_(sistrip::invalid_),
      maxDelay_(-1. * sistrip::invalid_),
      deviceWithMinDelay_(sistrip::invalid_),
      deviceWithMaxDelay_(sistrip::invalid_) {
  cout << "[" << __PRETTY_FUNCTION__ << "]"
       << " Created object for APV TIMING histograms" << endl;
}

// -----------------------------------------------------------------------------
FedTimingHistograms::~FedTimingHistograms() { cout << "[" << __PRETTY_FUNCTION__ << "]" << endl; }

// -----------------------------------------------------------------------------
void FedTimingHistograms::histoAnalysis(bool debug) {
  // Clear std::map holding analysis objects
  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 std::map containing std::vectors of profile histograms
  //   CollationsMap::const_iterator iter = collations().begin();
  //   for ( ; iter != collations().end(); iter++ ) {

  //     // Check std::vector of histos is not empty (should be 2 histos)
  //     if ( iter->second.empty() ) {
  //       cerr << "[" << __PRETTY_FUNCTION__ << "]"
  //       << " Zero collation histograms found!" << endl;
  //       continue;
  //     }

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

  //     // Perform histo analysis
  //     FedTimingAnalysis anal( iter->first );
  //     FedTimingAlgorithm algo( &anal );
  //     algo.analysis( profs );
  //     data_[iter->first] = anal;

  //     // Check tick height is valid
  //     if ( anal.height() < 100. ) {
  //       cerr << "[" << __PRETTY_FUNCTION__ << "]"
  //       << " Tick mark height too small: " << anal.height() << endl;
  //       continue;
  //     }

  //     // Check time of rising edge
  //     if ( anal.time() > sistrip::maximum_ ) { 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;
  //     }

  //   }

  //   cout << "[" << __PRETTY_FUNCTION__ << "]"
  //        << " Analyzed histograms for "
  //        << collations().size()
  //        << " FED channels" << endl;

  // Adjust maximum (and minimum) delay(s) to find optimum sampling point(s)
  if (time_max > sistrip::maximum_ || time_max < -1. * sistrip::maximum_) {
    cerr << "[" << __PRETTY_FUNCTION__ << "]"
         << " Unable to set maximum time! Found unexpected value: " << time_max << endl;
    return;
  }

  SiStripFecKey max(device_max);
  cout << " Device (FEC/slot/ring/CCU/module/channel) " << max.fecCrate() << "/" << max.fecSlot() << "/"
       << max.fecRing() << "/" << max.ccuAddr() << "/" << max.ccuChan() << "/"
       << " has maximum delay (rising edge) [ns]:" << time_max << endl;

  SiStripFecKey min(device_min);
  cout << " Device (FEC/slot/ring/CCU/module/channel): " << min.fecCrate() << "/" << min.fecSlot() << "/"
       << min.fecRing() << "/" << min.ccuAddr() << "/" << min.ccuChan() << "/"
       << " has minimum delay (rising edge) [ns]:" << time_min << endl;

  // Set maximum time for all analysis objects
  std::map<uint32_t, FedTimingAnalysis>::iterator ianal = data_.begin();
  for (; ianal != data_.end(); ianal++) {
    ianal->second.max(time_max);
    static uint16_t cntr = 0;
    if (debug) {
      std::stringstream ss;
      ianal->second.print(ss);
      cout << ss.str() << endl;
      cntr++;
    }
  }
}

// -----------------------------------------------------------------------------
void FedTimingHistograms::createSummaryHisto(const sistrip::Monitorable& histo,
                                             const sistrip::Presentation& type,
                                             const std::string& directory,
                                             const sistrip::Granularity& gran) {
  cout << "[" << __PRETTY_FUNCTION__ << "]" << endl;

  // Check view
  sistrip::View view = SiStripEnumsAndStrings::view(directory);
  if (view == sistrip::UNKNOWN_VIEW) {
    return;
  }

  // Analyze histograms
  histoAnalysis(false);

  // Extract data to be histogrammed
  factory_->init(histo, type, view, directory, gran);
  uint32_t xbins = factory_->extract(data_);

  // Create summary histogram (if it doesn't already exist)
  TH1* summary = histogram(histo, type, view, directory, xbins);

  // Fill histogram with data
  factory_->fill(*summary);
}