VpspScanTask.cc

Go to the documentation of this file.

#include "DQM/SiStripCommissioningSources/interface/VpspScanTask.h"
#include "DataFormats/SiStripCommon/interface/SiStripConstants.h"
#include "DataFormats/SiStripCommon/interface/SiStripFecKey.h"
#include "DataFormats/SiStripCommon/interface/SiStripHistoTitle.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <algorithm>

using namespace sistrip;

// -----------------------------------------------------------------------------
//
VpspScanTask::VpspScanTask( DQMStore* dqm,
                const FedChannelConnection& conn ) :
  CommissioningTask( dqm, conn, "VpspScanTask" ),
  vpsp_()
{}

// -----------------------------------------------------------------------------
//
VpspScanTask::~VpspScanTask() {
}

// -----------------------------------------------------------------------------
//
void VpspScanTask::book() {
  
  uint16_t nbins = 60;
 
  std::string title;

  vpsp_.resize(2);
  for ( uint16_t iapv = 0; iapv < 2; iapv++ ) {
    if ( connection().i2cAddr(iapv) ) { 

      std::stringstream extra_info; 
      extra_info << sistrip::apv_ << iapv;
      
      title = SiStripHistoTitle( sistrip::EXPERT_HISTO, 
                 sistrip::VPSP_SCAN, 
                 sistrip::FED_KEY, 
                 fedKey(),
                 sistrip::LLD_CHAN, 
                 connection().lldChannel(),
                 extra_info.str() ).title();
      
      vpsp_[iapv].histo( dqm()->bookProfile( title, title, 
                         nbins, -0.5, nbins*1.-0.5,
                         1025, 0., 1025. ) );
      
      vpsp_[iapv].vNumOfEntries_.resize(nbins,0);
      vpsp_[iapv].vSumOfContents_.resize(nbins,0);
      vpsp_[iapv].vSumOfSquares_.resize(nbins,0);
      
    }
  }
  
}

// -----------------------------------------------------------------------------
//
void VpspScanTask::fill( const SiStripEventSummary& summary,
             const edm::DetSet<SiStripRawDigi>& digis ) {

  // Retrieve VPSP setting and CCU channel
  uint32_t vpsp = const_cast<SiStripEventSummary&>(summary).vpsp();
  uint32_t ccu_chan = const_cast<SiStripEventSummary&>(summary).vpspCcuChan();

  // Check CCU channel from EventSummary is consistent with this module
  if ( SiStripFecKey( fecKey() ).ccuChan() != ccu_chan ) { return; }

  if ( digis.data.size() != 256 ) {
    edm::LogWarning(mlDqmSource_)
      << "[VpspScanTask::" << __func__ << "]"
      << " Unexpected number of digis! " 
      << digis.data.size(); 
    return;
  }

  // Fill histo with baseline(calc'ed from median value of data)
  for ( uint16_t iapv = 0; iapv < 2; iapv++ ) {
    
    if ( vpsp >= vpsp_[iapv].vNumOfEntries_.size() ) { 
      edm::LogWarning(mlDqmSource_)
    << "[VpspScanTask::" << __func__ << "]"
    << " Unexpected VPSP value! " << vpsp;
      return;
    }
    
    std::vector<uint16_t> baseline;
    baseline.reserve(128); 
    for ( uint16_t idigi = 128*iapv; idigi < 128*(iapv+1); idigi++ ) {
      baseline.push_back( digis.data[idigi].adc() ); 
    }
    sort( baseline.begin(), baseline.end() ); 
    uint16_t index = baseline.size()%2 ? baseline.size()/2 : baseline.size()/2-1;
    
    if ( !baseline.empty() ) { 
      updateHistoSet( vpsp_[iapv], vpsp, baseline[index] );
    }
    
  }

}

// -----------------------------------------------------------------------------
//
void VpspScanTask::update() {
  for ( uint32_t iapv = 0; iapv < vpsp_.size(); iapv++ ) {
    updateHistoSet( vpsp_[iapv] );
  }
}