PedestalsTask Class Reference

#include <PedestalsTask.h>

Inheritance diagram for PedestalsTask:

Public Member Functions
	PedestalsTask (DQMStore *, const FedChannelConnection &)
virtual	~PedestalsTask ()
Public Member Functions inherited from CommissioningTask
void	bookHistograms ()
	CommissioningTask (DQMStore *, const FedChannelConnection &, const std::string &my_name)
void	eventSetup (const edm::EventSetup *)
const uint32_t &	fillCntr () const
void	fillHistograms (const SiStripEventSummary &, const edm::DetSet< SiStripRawDigi > &)
void	fillHistograms (const SiStripEventSummary &, const uint16_t &fed_id, const std::map< uint16_t, float > &fed_ch)
const std::string &	myName () const
const uint32_t &	updateFreq () const
void	updateFreq (const uint32_t &)
void	updateHistograms ()
virtual	~CommissioningTask ()

Private Member Functions
virtual void	book ()
virtual void	fill (const SiStripEventSummary &, const edm::DetSet< SiStripRawDigi > &)
virtual void	update ()

Private Attributes
std::vector< HistoSet >	cm_
std::vector< HistoSet >	peds_

Additional Inherited Members
Protected Member Functions inherited from CommissioningTask
const FedChannelConnection &	connection () const
DQMStore *const	dqm () const
const edm::EventSetup *const	eventSetup () const
const uint32_t &	fecKey () const
const uint32_t &	fedKey () const
void	updateHistoSet (HistoSet &, const uint32_t &bin, const float &value)
void	updateHistoSet (CompactHistoSet &, const uint32_t &bin, const short &value)
void	updateHistoSet (HistoSet &, const uint32_t &bin)
void	updateHistoSet (CompactHistoSet &, const uint32_t &bin)
void	updateHistoSet (HistoSet &, const float &value)
void	updateHistoSet (CompactHistoSet &)
void	updateHistoSet (HistoSet &)

Detailed Description

Definition at line 10 of file PedestalsTask.h.

Constructor & Destructor Documentation

PedestalsTask::PedestalsTask	(	DQMStore *	dqm,
		const FedChannelConnection &	conn
	)

Definition at line 15 of file PedestalsTask.cc.

References LogTrace, and sistrip::mlDqmSource_.

                                                       :
  CommissioningTask( dqm, conn, "PedestalsTask" ),
  peds_(),
  cm_()
{
  LogTrace(mlDqmSource_)
    << "[PedestalsTask::" << __func__ << "]"
    << " Constructing object...";
}

PedestalsTask::~PedestalsTask ( )

virtual

Definition at line 28 of file PedestalsTask.cc.

References LogTrace, and sistrip::mlDqmSource_.

                              {
  LogTrace(mlDqmSource_)
    << "[PedestalsTask::" << __func__ << "]"
    << " Destructing object...";
}

Member Function Documentation

void PedestalsTask::book ( )

privatevirtual

Reimplemented from CommissioningTask.

Definition at line 36 of file PedestalsTask.cc.

References sistrip::APV, cm_, sistrip::extrainfo::commonMode_, CommissioningTask::connection(), CommissioningTask::dqm(), sistrip::EXPERT_HISTO, sistrip::FED_KEY, CommissioningTask::fedKey(), sistrip::LLD_CHAN, LogTrace, sistrip::mlDqmSource_, pileupCalc::nbins, sistrip::extrainfo::noise_, sistrip::PEDESTALS, sistrip::extrainfo::pedestals_, peds_, AlCaHLTBitMon_QueryRunRegistry::string, SiStripHistoTitle::title(), and indexGen::title.

                         {
  LogTrace(mlDqmSource_) << "[PedestalsTask::" << __func__ << "]";
  
  uint16_t nbins;
  std::string title;
  std::string extra_info;
  peds_.resize(2);
  nbins = 256;
  
  // Pedestals histogram
  extra_info = sistrip::extrainfo::pedestals_; 
  peds_[0].isProfile_ = true;
  
  title = SiStripHistoTitle( sistrip::EXPERT_HISTO, 
                 sistrip::PEDESTALS, 
                 sistrip::FED_KEY, 
                 fedKey(),
                 sistrip::LLD_CHAN, 
                 connection().lldChannel(),
                 extra_info ).title();
  
  peds_[0].histo( dqm()->bookProfile( title, title, 
                      nbins, -0.5, nbins*1.-0.5,
                      1025, 0., 1025. ) );
  
  peds_[0].vNumOfEntries_.resize(nbins,0);
  peds_[0].vSumOfContents_.resize(nbins,0);
  peds_[0].vSumOfSquares_.resize(nbins,0);

  // Noise histogram
  extra_info = sistrip::extrainfo::noise_; 
  peds_[1].isProfile_ = true;
  
  title = SiStripHistoTitle( sistrip::EXPERT_HISTO, 
                 sistrip::PEDESTALS, 
                 sistrip::FED_KEY, 
                 fedKey(),
                 sistrip::LLD_CHAN, 
                 connection().lldChannel(),
                 extra_info ).title();
  
  peds_[1].histo( dqm()->bookProfile( title, title, 
                      nbins, -0.5, nbins*1.-0.5,
                      1025, 0., 1025. ) );
  
  peds_[1].vNumOfEntries_.resize(nbins,0);
  peds_[1].vSumOfContents_.resize(nbins,0);
  peds_[1].vSumOfSquares_.resize(nbins,0);

  // Common mode histograms
  cm_.resize(2);
  nbins = 1024;
  for ( uint16_t iapv = 0; iapv < 2; iapv++ ) { 
    
    title = SiStripHistoTitle( sistrip::EXPERT_HISTO, 
                   sistrip::PEDESTALS, 
                   sistrip::FED_KEY, 
                   fedKey(),
                   sistrip::APV, 
                   connection().i2cAddr(iapv),
                   sistrip::extrainfo::commonMode_ ).title();
    
    cm_[iapv].histo( dqm()->book1D( title, title, nbins, -0.5, nbins*1.-0.5 ) );
    cm_[iapv].isProfile_ = false;
    
    cm_[iapv].vNumOfEntries_.resize(nbins,0);
    cm_[iapv].vNumOfEntries_.resize(nbins,0);
    
  }
  
}

void PedestalsTask::fill ( const SiStripEventSummary &  summary, const edm::DetSet< SiStripRawDigi > &  digis  )

privatevirtual

Reimplemented from CommissioningTask.

Definition at line 110 of file PedestalsTask.cc.

References ecalMGPA::adc(), cm_, edm::DetSet< T >::data, diffTreeTool::diff, cmsHarvester::index, sistrip::mlDqmSource_, pileupCalc::nbins, peds_, python.multivaluedict::sort(), and CommissioningTask::updateHistoSet().

                                                       {
  
  if ( digis.data.size() != peds_[0].vNumOfEntries_.size() ) {
    edm::LogWarning(mlDqmSource_)
      << "[PedestalsTask::" << __func__ << "]"
      << " Unexpected number of digis: " 
      << digis.data.size(); 
    return;
  }
  
  // Check number of digis
  uint16_t nbins = peds_[0].vNumOfEntries_.size();
  if ( digis.data.size() < nbins ) { nbins = digis.data.size(); }

  //@@ Inefficient!!!
  uint16_t napvs = nbins / 128;
  std::vector<uint32_t> cm; cm.resize(napvs,0);
  
  // Calc common mode for both APVs
  std::vector<uint16_t> adc;
  for ( uint16_t iapv = 0; iapv < napvs; iapv++ ) { 
    adc.clear(); adc.reserve(128);
    for ( uint16_t ibin = 0; ibin < 128; ibin++ ) { 
      if ( (iapv*128)+ibin < nbins ) { 
    adc.push_back( digis.data[(iapv*128)+ibin].adc() );  //@@ VIRGIN RAW DATA (MUX, APV READOUT)
      }
    }
    sort( adc.begin(), adc.end() ); 
    uint16_t index = adc.size()%2 ? adc.size()/2 : adc.size()/2-1;
    if ( !adc.empty() ) { cm[iapv] = static_cast<uint32_t>( adc[index] ); }
  }
  
  for ( uint16_t ibin = 0; ibin < nbins; ibin++ ) {
    updateHistoSet( peds_[0], ibin, digis.data[ibin].adc() ); // peds and raw noise
    float diff = static_cast<float>( digis.data[ibin].adc() ) - static_cast<float>( cm[ibin/128] );
    updateHistoSet( peds_[1], ibin, diff ); // residuals and real noise
  }
  
  if ( cm.size() < cm_.size() ) {
    edm::LogWarning(mlDqmSource_)
      << "[PedestalsTask::" << __func__ << "]"
      << " Fewer CM values than expected: " << cm.size();
  }
  
  updateHistoSet( cm_[0], cm[0] );
  updateHistoSet( cm_[1], cm[1] );
  
}

void PedestalsTask::update ( )

privatevirtual

Reimplemented from CommissioningTask.

Definition at line 162 of file PedestalsTask.cc.

References cm_, python.tagInventory::entries, relativeConstraints::error, extract(), timingPdfMaker::histo, cuy::ii, timingPdfMaker::mean, HLT_25ns14e33_v1_cff::noise, peds_, UpdateTProfile::setBinContent(), mathSSE::sqrt(), and CommissioningTask::updateHistoSet().

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

                           {
  
  // Pedestals 
  updateHistoSet( peds_[0] );
  
  // Noise (cannot use HistoSet directly, as want to plot noise as "contents", not "error")
  TProfile* histo = ExtractTObject<TProfile>().extract( peds_[1].histo() );
  for ( uint16_t ii = 0; ii < peds_[1].vNumOfEntries_.size(); ++ii ) {
    
    float mean    = 0.;
    float spread  = 0.;
    float entries = peds_[1].vNumOfEntries_[ii];
    if ( entries > 0. ) {
      mean = peds_[1].vSumOfContents_[ii] / entries;
      spread = sqrt( peds_[1].vSumOfSquares_[ii] / entries - mean * mean ); 
    }
    
    float noise = spread;
    float error = 0; // sqrt(entries) / entries;
    
    UpdateTProfile::setBinContent( histo, ii+1, entries, noise, error );
    
  }
  
  // Common mode
  updateHistoSet( cm_[0] );
  updateHistoSet( cm_[1] );
  
}

Member Data Documentation

std::vector<HistoSet> PedestalsTask::cm_

private

Definition at line 25 of file PedestalsTask.h.

Referenced by book(), fill(), and update().

std::vector<HistoSet> PedestalsTask::peds_

private

Definition at line 24 of file PedestalsTask.h.

Referenced by book(), fill(), and update().