DaqScopeModeTask Class Reference

#include <DaqScopeModeTask.h>

Inheritance diagram for DaqScopeModeTask:

Public Member Functions
	DaqScopeModeTask (DQMStore *, const FedChannelConnection &, const edm::ParameterSet &)
	~DaqScopeModeTask () override
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 edm::DetSet< SiStripRawDigi > &, 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
void	book () override
void	fill (const SiStripEventSummary &, const edm::DetSet< SiStripRawDigi > &) override
void	fill (const SiStripEventSummary &, const edm::DetSet< SiStripRawDigi > &, const edm::DetSet< SiStripRawDigi > &) override
void	update () override

Private Attributes
std::vector< HistoSet >	cm_
uint16_t	nBins_
uint16_t	nBinsSpy_
edm::ParameterSet	parameters_
	parameters useful for the spy More...
std::vector< HistoSet >	peds_
HistoSet	scopeFrame_

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 DaqScopeModeTask.h.

Constructor & Destructor Documentation

DaqScopeModeTask::DaqScopeModeTask	(	DQMStore *	dqm,
		const FedChannelConnection &	conn,
		const edm::ParameterSet &	pset
	)

Definition at line 13 of file DaqScopeModeTask.cc.

                                                  :
  CommissioningTask( dqm, conn, "DaqScopeModeTask" ),
  scopeFrame_(),
  nBins_(256), //@@ number of strips per FED channel
  nBinsSpy_(298), //@@ in case of spy events includes header and trailing tick
  parameters_(pset)
{}

DaqScopeModeTask::~DaqScopeModeTask ( )

override

Definition at line 25 of file DaqScopeModeTask.cc.

                                    {
}

Member Function Documentation

void DaqScopeModeTask::book ( )

overrideprivatevirtual

Reimplemented from CommissioningTask.

Definition at line 30 of file DaqScopeModeTask.cc.

References sistrip::APV, cm_, sistrip::extrainfo::commonMode_, CommissioningTask::connection(), sistrip::DAQ_SCOPE_MODE, CommissioningTask::dqm(), edm::ParameterSet::existsAs(), sistrip::EXPERT_HISTO, sistrip::FED_KEY, CommissioningTask::fedKey(), edm::ParameterSet::getParameter(), CommissioningTask::HistoSet::histo(), CommissioningTask::HistoSet::isProfile_, sistrip::LLD_CHAN, LogTrace, sistrip::mlDqmSource_, pileupCalc::nbins, nBins_, nBinsSpy_, sistrip::extrainfo::noise_, parameters_, sistrip::extrainfo::pedestals_, peds_, scopeFrame_, sistrip::extrainfo::scopeModeFrame_, AlCaHLTBitMon_QueryRunRegistry::string, fftjetcommon_cfi::title, SiStripHistoTitle::title(), CommissioningTask::HistoSet::vNumOfEntries_, CommissioningTask::HistoSet::vSumOfContents_, and CommissioningTask::HistoSet::vSumOfSquares_.

                            {
  LogTrace(mlDqmSource_) << "[CommissioningTask::" << __func__ << "]";

  std::string title;
  if(not (parameters_.existsAs<bool>("isSpy") and parameters_.getParameter<bool>("isSpy"))){

    title = SiStripHistoTitle( sistrip::EXPERT_HISTO, 
                   sistrip::DAQ_SCOPE_MODE, 
                   sistrip::FED_KEY, 
                   fedKey(),
                   sistrip::LLD_CHAN, 
                   connection().lldChannel(),
                   sistrip::extrainfo::scopeModeFrame_).title();
    
    scopeFrame_.histo( dqm()->book1D( title, title, 
                      nBins_, -0.5, nBins_-0.5 ) );
    
    scopeFrame_.vNumOfEntries_.resize(nBins_,0);
    scopeFrame_.vSumOfContents_.resize(nBins_,0);
    scopeFrame_.vSumOfSquares_.resize(nBins_,0);
    scopeFrame_.isProfile_ = false; 
  }
  else{ // use spy run to measure pedestal and tick height
    
    std::string extra_info;
    peds_.resize(2);
    extra_info = sistrip::extrainfo::pedestals_;

    peds_[0].isProfile_ = true;
    title = SiStripHistoTitle( sistrip::EXPERT_HISTO,
                   sistrip::DAQ_SCOPE_MODE,
                   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::DAQ_SCOPE_MODE,
                   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);
    int nbins = 1024;

    for ( uint16_t iapv = 0; iapv < 2; iapv++ ) {
      title = SiStripHistoTitle( sistrip::EXPERT_HISTO,
                 sistrip::DAQ_SCOPE_MODE,
                 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);      
    }


    title = SiStripHistoTitle( sistrip::EXPERT_HISTO,
                               sistrip::DAQ_SCOPE_MODE,
                               sistrip::FED_KEY,
                               fedKey(),
                               sistrip::LLD_CHAN,
                               connection().lldChannel(),
                   sistrip::extrainfo::scopeModeFrame_).title();

    scopeFrame_.histo( dqm()->bookProfile( title, title,
                       nBinsSpy_, -0.5, nBinsSpy_*1.-0.5,
                       1025, 0., 1025. ) );

    scopeFrame_.vNumOfEntries_.resize(nBinsSpy_,0);
    scopeFrame_.vSumOfContents_.resize(nBinsSpy_,0);
    scopeFrame_.vSumOfSquares_.resize(nBinsSpy_,0);
    scopeFrame_.isProfile_ = true;
    

  }
}

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

overrideprivatevirtual

Reimplemented from CommissioningTask.

Definition at line 141 of file DaqScopeModeTask.cc.

References ecalMGPA::adc(), create_public_pileup_plots::bins, cm_, edm::DetSet< T >::data, diffTreeTool::diff, edm::ParameterSet::existsAs(), CommissioningTask::fillCntr(), edm::ParameterSet::getParameter(), sistrip::mlDqmSource_, pileupCalc::nbins, nBins_, parameters_, peds_, scopeFrame_, CommissioningTask::updateFreq(), and CommissioningTask::updateHistoSet().

                                                          {
  
  // Only fill every 'N' events 
  if(not (parameters_.existsAs<bool>("isSpy") and parameters_.getParameter<bool>("isSpy"))){
    if ( !updateFreq() || fillCntr()%updateFreq() ) { return; }
  }

  if ( digis.data.size() != nBins_ ) { //@@ check scope mode length?  
    edm::LogWarning(mlDqmSource_)
      << "[DaqScopeModeTask::" << __func__ << "]"
      << " Unexpected number of digis (" 
      << digis.data.size()
      << ") wrt number of histogram bins ("
      << nBins_ << ")!";
  }

 
  if(not (parameters_.existsAs<bool>("isSpy") and parameters_.getParameter<bool>("isSpy"))){
    uint16_t bins = digis.data.size() < nBins_ ? digis.data.size() : nBins_;
    for ( uint16_t ibin = 0; ibin < bins; ibin++ ) {
      updateHistoSet( scopeFrame_, ibin, digis.data[ibin].adc() ); 
    }
  }
  else{
    // fill the pedestal histograms as done in the pedestal task                                                                                                                                     
    // Check number of digis                                                                                                                                                                          
    uint16_t nbins = peds_[0].vNumOfEntries_.size();
    if ( digis.data.size() < nbins ) { nbins = digis.data.size(); }

    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() );
        }
      }

      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++ ) {
      float digiVal = digis.data[ibin].adc();
      updateHistoSet( peds_[0], ibin, digiVal ); // peds and raw noise                                                                                                                                
      float diff = digiVal - 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 DaqScopeModeTask::fill ( const SiStripEventSummary &  summary, const edm::DetSet< SiStripRawDigi > &  digis, const edm::DetSet< SiStripRawDigi > &  digisAlt  )

overrideprivatevirtual

Calculate pedestal

Reimplemented from CommissioningTask.

Definition at line 210 of file DaqScopeModeTask.cc.

References ecalMGPA::adc(), create_public_pileup_plots::bins, cm_, edm::DetSet< T >::data, diffTreeTool::diff, edm::ParameterSet::existsAs(), CommissioningTask::fillCntr(), edm::ParameterSet::getParameter(), sistrip::mlDqmSource_, pileupCalc::nbins, nBins_, nBinsSpy_, parameters_, peds_, scopeFrame_, CommissioningTask::updateFreq(), and CommissioningTask::updateHistoSet().

                                                            {
  
  // Only fill every 'N' events 
  if(not (parameters_.existsAs<bool>("isSpy") and parameters_.getParameter<bool>("isSpy"))){
    if ( !updateFreq() || fillCntr()%updateFreq() ) { return; }
  }

  if ( digis.data.size() != nBins_ ) { //@@ check scope mode length?  
    edm::LogWarning(mlDqmSource_)
      << "[DaqScopeModeTask::" << __func__ << "]"
      << " Unexpected number of digis (" 
      << digis.data.size()
      << ") wrt number of histogram bins ("
      << nBins_ << ")!";
  }

  if(not (parameters_.existsAs<bool>("isSpy") and parameters_.getParameter<bool>("isSpy"))){
    uint16_t bins = digis.data.size() < nBins_ ? digis.data.size() : nBins_;
    for ( uint16_t ibin = 0; ibin < bins; ibin++ ) {
      updateHistoSet( scopeFrame_, ibin, digis.data[ibin].adc() ); 
    }
  }
  else{
    // fill the pedestal histograms as done in the pedestal task                                                                                                                                      
    // Check number of digis                                                                                                                                                                         
    uint16_t nbins = peds_[0].vNumOfEntries_.size();
    if ( digis.data.size() < nbins ) { nbins = digis.data.size(); }

    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() );
        }
      }
      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++ ) {
      float digiVal = digis.data[ibin].adc();
      updateHistoSet( peds_[0], ibin, digiVal ); // peds and raw noise                                                                                                                                
      float diff = digiVal - 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] );

    uint16_t bins = digisAlt.data.size() < nBinsSpy_ ? digisAlt.data.size() : nBinsSpy_;
    for ( uint16_t ibin = 0; ibin < bins; ibin++ ) {
      updateHistoSet( scopeFrame_, ibin, digisAlt.data[ibin].adc() );
    }
  }
}

void DaqScopeModeTask::update ( )

overrideprivatevirtual

Reimplemented from CommissioningTask.

Definition at line 283 of file DaqScopeModeTask.cc.

References cm_, relativeConstraints::error, edm::ParameterSet::existsAs(), extract(), edm::ParameterSet::getParameter(), trackerHits::histo, cuy::ii, SiStripPI::mean, parameters_, peds_, scopeFrame_, 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(), dqm-mbProfile.Profile::finish(), 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().

                              {

  if(not (parameters_.existsAs<bool>("isSpy") and parameters_.getParameter<bool>("isSpy")))
    updateHistoSet( scopeFrame_ );
  else{
    
    updateHistoSet( peds_[0] );   
    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 );      
    }
    updateHistoSet( cm_[0] );
    updateHistoSet( cm_[1] );
    updateHistoSet( scopeFrame_ );

  }
}

Member Data Documentation

std::vector<HistoSet> DaqScopeModeTask::cm_

private

Definition at line 34 of file DaqScopeModeTask.h.

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

uint16_t DaqScopeModeTask::nBins_

private

Definition at line 36 of file DaqScopeModeTask.h.

Referenced by book(), and fill().

uint16_t DaqScopeModeTask::nBinsSpy_

private

Definition at line 37 of file DaqScopeModeTask.h.

Referenced by book(), and fill().

edm::ParameterSet DaqScopeModeTask::parameters_

private

parameters useful for the spy

Definition at line 40 of file DaqScopeModeTask.h.

Referenced by book(), Mixins._TypedParameterizable::clone(), Types.PSet::clone(), Mixins._TypedParameterizable::copy(), fill(), and update().

std::vector<HistoSet> DaqScopeModeTask::peds_

private

Definition at line 33 of file DaqScopeModeTask.h.

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

HistoSet DaqScopeModeTask::scopeFrame_

private

Definition at line 30 of file DaqScopeModeTask.h.

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