CalibrationHistograms.cc

Go to the documentation of this file.

#include "DQM/SiStripCommissioningClients/interface/CalibrationHistograms.h"
#include "CondFormats/SiStripObjects/interface/CalibrationAnalysis.h"
#include "DQM/SiStripCommissioningAnalysis/interface/CalibrationAlgorithm.h"
#include "DQM/SiStripCommissioningSummary/interface/CalibrationSummaryFactory.h"
#include "DataFormats/SiStripCommon/interface/SiStripConstants.h"
#include "DQM/SiStripCommon/interface/ExtractTObject.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <iostream>
#include <sstream>
#include <iomanip>
#include "TH1F.h"

using namespace std;
using namespace sistrip;

// -----------------------------------------------------------------------------
CalibrationHistograms::CalibrationHistograms( const edm::ParameterSet& pset,
                                              DQMStore* bei,
                                              const sistrip::RunType& task )
  : CommissioningHistograms( pset.getParameter<edm::ParameterSet>("CalibrationParameters"),
                             bei,
                             task ),
    calchan_(0),
    isha_(-1),
    vfs_(-1)
{
  LogTrace(mlDqmClient_) 
       << "[CalibrationHistograms::" << __func__ << "]"
       << " Constructing object...";
  std::string temp = std::string(sistrip::collate_) + "/" + sistrip::root_;
  factory_ = auto_ptr<CalibrationSummaryFactory>( new CalibrationSummaryFactory );
  MonitorElement* calchanElement = bei->get(temp + "/calchan");
  if(calchanElement) {
    calchan_ = calchanElement->getIntValue();
    edm::LogVerbatim(mlDqmClient_)
      << "[CalibrationHistograms::" << __func__ << "]"
      << "CalChan value is " << calchan_;
  } else {
    edm::LogWarning(mlDqmClient_)
      << "[CalibrationHistograms::" << __func__ << "]"
      << "CalChan value not found. Using " << calchan_;
  }
  MonitorElement* ishaElement = bei->get(temp + "/isha");
  if(ishaElement) isha_ = ishaElement->getIntValue() ;
  MonitorElement* vfsElement = bei->get(temp + "/vfs");
  if(vfsElement) vfs_ = vfsElement->getIntValue() ;
}

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

// -----------------------------------------------------------------------------     
void CalibrationHistograms::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_)
       << "[CalibrationHistograms::" << __func__ << "]"
       << " Zero collation histograms found!";
      continue;
    }
    
    // Retrieve pointers to 1D histos for this FED channel 
    vector<TH1*> profs;
    Histos::const_iterator ihis = iter->second.begin();
    for ( ; ihis != iter->second.end(); ihis++ ) {
      TH1F* prof = ExtractTObject<TH1F>().extract( (*ihis)->me_ );
      if ( prof ) { profs.push_back(prof); }
    } 
    // Perform histo analysis 
    CalibrationAnalysis* anal = new CalibrationAnalysis( iter->first, (task()==sistrip::CALIBRATION_DECO), calchan_ );
    CalibrationAlgorithm algo( this->pset(), anal );
    algo.analysis( profs );
    data()[iter->first] = anal; 
    
 }
 
}

// -----------------------------------------------------------------------------
void CalibrationHistograms::printAnalyses() {
  Analyses::iterator ianal = data().begin();
  Analyses::iterator janal = data().end();
  for ( ; ianal != janal; ++ianal ) {
    if ( ianal->second ) {
      std::stringstream ss;
      ianal->second->print( ss,0 );
      ianal->second->print( ss,1 );
      if ( ianal->second->isValid() ) { LogTrace(mlDqmClient_) << ss.str();
      } else { edm::LogWarning(mlDqmClient_) << ss.str(); }
    }
  }
}