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 edm::DetSet< SiStripRawDigi > &, const edm::DetSet< SiStripRawDigi > &, const std::vector< uint16_t > &)
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	fill (const SiStripEventSummary &, const edm::DetSet< SiStripRawDigi > &, const edm::DetSet< SiStripRawDigi > &, const std::vector< uint16_t > &) override
void	update () override

Private Attributes
std::vector< HistoSet >	cm_
HistoSet	highHeader_
HistoSet	lowHeader_
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 12 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(), highHeader_, CommissioningTask::HistoSet::histo(), CommissioningTask::HistoSet::isProfile_, sistrip::LLD_CHAN, LogTrace, lowHeader_, sistrip::mlDqmSource_, pileupCalc::nbins, nBins_, nBinsSpy_, sistrip::extrainfo::noise_, parameters_, sistrip::extrainfo::pedestals_, peds_, scopeFrame_, sistrip::extrainfo::scopeModeFrame_, sistrip::extrainfo::scopeModeHeaderHigh_, sistrip::extrainfo::scopeModeHeaderLow_, AlCaHLTBitMon_QueryRunRegistry::string, SiStripHistoTitle::title(), runGCPTkAlMap::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);      
    }

    // high and low header histograms
    title = SiStripHistoTitle( sistrip::EXPERT_HISTO,
                sistrip::DAQ_SCOPE_MODE,
                sistrip::FED_KEY,
                fedKey(),
                sistrip::LLD_CHAN,
                connection().lldChannel(),
                sistrip::extrainfo::scopeModeHeaderLow_ ).title();

    lowHeader_.histo( dqm()->book1D( title, title, nbins, -0.5, 1024*1.-0.5 ) );
    lowHeader_.isProfile_ = false;
    lowHeader_.vNumOfEntries_.resize(nbins,0);


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

    highHeader_.histo( dqm()->book1D( title, title, nbins, -0.5, 1024*1.-0.5 ) );
    highHeader_.isProfile_ = false;
    highHeader_.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 167 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_, jetUpdater_cfi::sort, 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 236 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(), highHeader_, mps_fire::i, lowHeader_, sistrip::mlDqmSource_, pileupCalc::nbins, nBins_, nBinsSpy_, parameters_, peds_, scopeFrame_, jetUpdater_cfi::sort, 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() ); 
    }

     // Header low and high for both APVs
    std::vector<uint32_t> adcHeader_high;
    std::vector<uint32_t> adcHeader_low;

    float threshold_high = (digisAlt.data[286].adc() + digisAlt.data[287].adc())/4;
    float threshold_low  = 100;
    int   minNumberForHeader = 4;
    bool  goodHeaderFound   = false;
    int   nConsecutiveHigh  = 0;
    adcHeader_high.clear();
    adcHeader_high.reserve(30);
    adcHeader_low.clear();
    adcHeader_low.reserve(30);

    for ( uint16_t ibin = 6; ibin < 11; ibin++ ) {
        if (digisAlt.data[ibin].adc() > threshold_high) {
            nConsecutiveHigh++;
        }
    }

    if(nConsecutiveHigh > minNumberForHeader) goodHeaderFound = true; // if nConsecutiveHigh > 4 --> good header found
    if(goodHeaderFound == false) return;
    for ( uint16_t ibin = 0; ibin < 30; ibin++ ) {
      if(digisAlt.data[ibin].adc() > threshold_high && goodHeaderFound){ // save of samples above avg(trailing ticks)/4
        adcHeader_high.push_back( digisAlt.data[ibin].adc() );
      }
      if (digisAlt.data[ibin].adc() < threshold_low && goodHeaderFound){
        adcHeader_low.push_back( digisAlt.data[ibin].adc() );
      }
    }
    if(adcHeader_low.empty() || adcHeader_high.empty()) {
        return;
    }
    for ( uint16_t i = 0; i<adcHeader_low.size(); i++){
        updateHistoSet( lowHeader_, adcHeader_low[i]);
    }
    for ( uint16_t i = 0; i<adcHeader_high.size(); i++){
       updateHistoSet( highHeader_, adcHeader_high[i]);
    }



   
  }      
}

void DaqScopeModeTask::fill ( const SiStripEventSummary &  summary, const edm::DetSet< SiStripRawDigi > &  digis, const edm::DetSet< SiStripRawDigi > &  digisAlt, const std::vector< uint16_t > &  stripOnCluster  )

overrideprivatevirtual

Calculate pedestal

Reimplemented from CommissioningTask.

Definition at line 355 of file DaqScopeModeTask.cc.

References ecalMGPA::adc(), create_public_pileup_plots::bins, cm_, edm::DetSet< T >::data, diffTreeTool::diff, edm::ParameterSet::existsAs(), CommissioningTask::fillCntr(), spr::find(), edm::ParameterSet::getParameter(), highHeader_, mps_fire::i, lowHeader_, sistrip::mlDqmSource_, pileupCalc::nbins, nBins_, nBinsSpy_, parameters_, peds_, scopeFrame_, jetUpdater_cfi::sort, 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
    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 ) { 
      if(std::find(stripOnCluster.begin(),stripOnCluster.end(),(iapv*128)+ibin) == stripOnCluster.end()) // if not found, strip is good
        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++ ) {
      if(std::find(stripOnCluster.begin(),stripOnCluster.end(),ibin) != stripOnCluster.end()){
    continue;
      }
      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() ); 
    }    
     // Header low and high for both APVs
    std::vector<uint32_t> adcHeader_high;
    std::vector<uint32_t> adcHeader_low;

    float threshold_high = (digisAlt.data[286].adc() + digisAlt.data[287].adc())/4;
    float threshold_low  = 120;
    int   minNumberForHeader = 4;
    bool  goodHeaderFound   = false;
    int   nConsecutiveHigh  = 0;
    adcHeader_high.clear();
    adcHeader_high.reserve(30);
    adcHeader_low.clear();
    adcHeader_low.reserve(30);

    for ( uint16_t ibin = 6; ibin < 11; ibin++ ) {
        if (digisAlt.data[ibin].adc() > threshold_high) {
            nConsecutiveHigh++;
        }
    }

    if(nConsecutiveHigh > minNumberForHeader) goodHeaderFound = true; // if nConsecutiveHigh > 4 --> good header found
    if(goodHeaderFound == false) return;
    for ( uint16_t ibin = 0; ibin < 30; ibin++ ) {
      if(digisAlt.data[ibin].adc() > threshold_high && goodHeaderFound){ // save of samples above avg(trailing ticks)/4
        adcHeader_high.push_back( digisAlt.data[ibin].adc() );
      }
      if (digisAlt.data[ibin].adc() < threshold_low && goodHeaderFound){
        adcHeader_low.push_back( digisAlt.data[ibin].adc() );
      }
    }
    if(adcHeader_low.empty() || adcHeader_high.empty()) {
        return;
    }
    for ( uint16_t i = 0; i<adcHeader_low.size(); i++){
        updateHistoSet( lowHeader_, adcHeader_low[i]);
    }
    for ( uint16_t i = 0; i<adcHeader_high.size(); i++){
       updateHistoSet( highHeader_, adcHeader_high[i]);
    }


  }      
}

void DaqScopeModeTask::update ( )

overrideprivatevirtual

Reimplemented from CommissioningTask.

Definition at line 471 of file DaqScopeModeTask.cc.

References cm_, relativeConstraints::error, edm::ParameterSet::existsAs(), extract(), edm::ParameterSet::getParameter(), highHeader_, trackerHits::histo, cuy::ii, lowHeader_, 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_ );
    updateHistoSet( lowHeader_ );
    updateHistoSet( highHeader_ );
  }
}

Member Data Documentation

std::vector<HistoSet> DaqScopeModeTask::cm_

private

Definition at line 42 of file DaqScopeModeTask.h.

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

HistoSet DaqScopeModeTask::highHeader_

private

Definition at line 46 of file DaqScopeModeTask.h.

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

HistoSet DaqScopeModeTask::lowHeader_

private

Definition at line 45 of file DaqScopeModeTask.h.

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

uint16_t DaqScopeModeTask::nBins_

private

Definition at line 48 of file DaqScopeModeTask.h.

Referenced by book(), and fill().

uint16_t DaqScopeModeTask::nBinsSpy_

private

Definition at line 49 of file DaqScopeModeTask.h.

Referenced by book(), and fill().

edm::ParameterSet DaqScopeModeTask::parameters_

private

parameters useful for the spy

Definition at line 52 of file DaqScopeModeTask.h.

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

std::vector<HistoSet> DaqScopeModeTask::peds_

private

Definition at line 41 of file DaqScopeModeTask.h.

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

HistoSet DaqScopeModeTask::scopeFrame_

private

Definition at line 38 of file DaqScopeModeTask.h.

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