CalibrationScanSummaryFactory.cc

Go to the documentation of this file.

#include "DQM/SiStripCommissioningSummary/interface/CalibrationScanSummaryFactory.h"
#include "CondFormats/SiStripObjects/interface/CalibrationScanAnalysis.h"
#include "DataFormats/SiStripCommon/interface/SiStripEnumsAndStrings.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <iostream>
#include <sstream>

using namespace sistrip;

// -----------------------------------------------------------------------------
//
void CalibrationScanSummaryFactory::extract( Iterator iter ) {

  CalibrationScanAnalysis* anal = dynamic_cast<CalibrationScanAnalysis*>( iter->second );
  if ( !anal ) { return; }

  SiStripFecKey lldKey = SiStripFecKey(iter->first);

  uint32_t key1 = SiStripFecKey(lldKey.fecCrate(),
                                lldKey.fecSlot(),
                                lldKey.fecRing(),
                                lldKey.ccuAddr(),
                                lldKey.ccuChan(),
                                lldKey.lldChan(),
                                lldKey.i2cAddr(lldKey.lldChan(),true)).key();

  uint32_t key2 = SiStripFecKey(lldKey.fecCrate(),
                                lldKey.fecSlot(),
                                lldKey.fecRing(),
                                lldKey.ccuAddr(),
                                lldKey.ccuChan(),
                                lldKey.lldChan(),
                                lldKey.i2cAddr(lldKey.lldChan(),false)).key();

  std::vector<float> value (2,1*sistrip::invalid_);
  
  // search for the best isha and VFS values cose to the optimal ones
  if ( mon_ == sistrip::CALIBRATION_AMPLITUDE_TUNED) {
    value[0] = anal->tunedAmplitude()[0];
    value[1] = anal->tunedAmplitude()[1];
  } 
  else if ( mon_ == sistrip::CALIBRATION_BASELINE_TUNED) {
    value[0] = anal->tunedBaseline()[0];
    value[1] = anal->tunedBaseline()[1];
  } 
  else if ( mon_ == sistrip::CALIBRATION_TURNON_TUNED) {
    value[0] = anal->tunedTurnOn()[0];
    value[1] = anal->tunedTurnOn()[1];
  } 
  else if ( mon_ == sistrip::CALIBRATION_RISETIME_TUNED) {
    value[0] = anal->tunedRiseTime()[0];
    value[1] = anal->tunedRiseTime()[1];
  } 
  else if ( mon_ == sistrip::CALIBRATION_DECAYTIME_TUNED) {
    value[0] = anal->tunedDecayTime()[0];
    value[1] = anal->tunedDecayTime()[1];
  } 
  else if ( mon_ == sistrip::CALIBRATION_PEAKTIME_TUNED) {
    value[0] = anal->tunedPeakTime()[0];
    value[1] = anal->tunedPeakTime()[1];
  } 
  else if ( mon_ == sistrip::CALIBRATION_UNDERSHOOT_TUNED) {
    value[0] = anal->tunedUndershoot()[0];
    value[1] = anal->tunedUndershoot()[1];
  } 
  else if ( mon_ == sistrip::CALIBRATION_TAIL_TUNED) {
    value[0] = anal->tunedTail()[0];
    value[1] = anal->tunedTail()[1];
  } 
  else if ( mon_ == sistrip::CALIBRATION_SMEARING_TUNED) {
    value[0] = anal->tunedSmearing()[0];
    value[1] = anal->tunedSmearing()[1];
  } 
  else if ( mon_ == sistrip::CALIBRATION_CHI2_TUNED) {
    value[0] = anal->tunedChi2()[0];
    value[1] = anal->tunedChi2()[1];
  } 
  else if ( mon_ == sistrip::CALIBRATION_ISHA_TUNED) {
    value[0] = anal->tunedISHA()[0];
    value[1] = anal->tunedISHA()[1];
  } 
  else if ( mon_ == sistrip::CALIBRATION_VFS_TUNED) {
    value[0] = anal->tunedVFS()[0];
    value[1] = anal->tunedVFS()[1];
  } 
  else if ( mon_ == sistrip::CALIBRATION_ISHA_BEST) {
    value[0] = anal->bestISHA()[0];
    value[1] = anal->bestISHA()[1];
  } 
  else if ( mon_ == sistrip::CALIBRATION_VFS_BEST) {
    value[0] = anal->bestVFS()[0];
    value[1] = anal->bestVFS()[1];
  } 
  else { 
    edm::LogWarning(mlSummaryPlots_)
      << "[SummaryPlotFactory::" << __func__ << "]"
      << " Unexpected monitorable: "
      << SiStripEnumsAndStrings::monitorable( SummaryPlotFactoryBase::mon_ );
    return; 
  }
  
  SummaryPlotFactoryBase::generator_->fillMap( SummaryPlotFactoryBase::level_,
                           SummaryPlotFactoryBase::gran_,
                           key1,
                           value[0] );
    
  SummaryPlotFactoryBase::generator_->fillMap( SummaryPlotFactoryBase::level_,
                           SummaryPlotFactoryBase::gran_,
                           key2,
                           value[1] );

  // set the x-axis
  format();
  
}

//------------------------------------------------------------------------------
//
void CalibrationScanSummaryFactory::format() {

  // Histogram formatting
  if ( mon_ == sistrip::CALIBRATION_AMPLITUDE_TUNED) 
    generator_->axisLabel( "Amplitude (ADC)" );

  else if ( mon_ == sistrip::CALIBRATION_BASELINE_TUNED)
    generator_->axisLabel( "Baseline (ADC)" );

  else if ( mon_ == sistrip::CALIBRATION_TURNON_TUNED)
    generator_->axisLabel( "Turn-On (ns)" );

  else if ( mon_ == sistrip::CALIBRATION_TAIL_TUNED)
    generator_->axisLabel( "Tail (%)" );

  else if ( mon_ == sistrip::CALIBRATION_RISETIME_TUNED)
    generator_->axisLabel( "Rise Time (ns)" );

  else if ( mon_ == sistrip::CALIBRATION_PEAKTIME_TUNED)
    generator_->axisLabel( "Peak Time (ns)" );
   
  else if ( mon_ == sistrip::CALIBRATION_DECAYTIME_TUNED)
    generator_->axisLabel( "Decay Time (ns)" );

  else if ( mon_ == sistrip::CALIBRATION_SMEARING_TUNED)
    generator_->axisLabel( "Smearing (ns)" );
  
  else if ( mon_ == sistrip::CALIBRATION_CHI2_TUNED)
    generator_->axisLabel( "Chi2/ndf" );

  else if ( mon_ == sistrip::CALIBRATION_VFS_TUNED or
        mon_ == sistrip::CALIBRATION_VFS_BEST)
    generator_->axisLabel( "VFS" );

  else if ( mon_ == sistrip::CALIBRATION_ISHA_TUNED or
        mon_ == sistrip::CALIBRATION_ISHA_BEST)
    generator_->axisLabel( "ISHA" );
  else { 
    edm::LogWarning(mlSummaryPlots_) 
      << "[SummaryPlotFactory::" << __func__ << "]"
      <<  " Unexpected SummaryHisto value:"
      << SiStripEnumsAndStrings::monitorable( SummaryPlotFactoryBase::mon_ ) ;
  } 
}