sistrip::RawToDigiUnpacker Class Reference

#include <SiStripRawToDigiUnpacker.h>

Classes
class	Registry
	private class to register start and end index of digis in a collection More...

Public Types
typedef edm::DetSetVector < SiStripDigi >	Digis
typedef edm::DetSetVector < SiStripRawDigi >	RawDigis

Public Member Functions
void	createDigis (const SiStripFedCabling &, const FEDRawDataCollection &, SiStripEventSummary &, RawDigis &scope_mode, RawDigis &virgin_raw, RawDigis &proc_raw, Digis &zero_suppr, DetIdCollection &, RawDigis &common_mode)
	creates digis More...
void	doAPVEmulatorCheck (bool)
void	doFullCorruptBufferChecks (bool)
void	extractCm (bool)
void	locateStartOfFedBuffer (const uint16_t &fed_id, const FEDRawData &input, FEDRawData &output)
	Removes any data appended prior to FED buffer and reorders 32-bit words if swapped. More...
void	quiet (bool)
	verbosity More...
	RawToDigiUnpacker (int16_t appended_bytes, int16_t fed_buffer_dump_freq, int16_t fed_event_dump_freq, int16_t trigger_fed_id, bool using_fed_key, bool unpack_bad_channels, bool mark_missing_feds, const uint32_t errorThreshold)
	constructor More...
void	triggerFed (const FEDRawDataCollection &, SiStripEventSummary &, const uint32_t &event)
	trigger info More...
void	useDaqRegister (bool)
	EventSummary update request -> not yet implemented for FEDBuffer class. More...
	~RawToDigiUnpacker ()
	default constructor More...

Private Member Functions
void	cleanupWorkVectors ()
	method to clear registries and digi collections More...
sistrip::FedBufferFormat	fedBufferFormat (const uint16_t &register_value)
	returns buffer format More...
sistrip::FedReadoutMode	fedReadoutMode (const uint16_t &register_value)
	returns buffer readout mode More...
void	handleException (std::string method_name, std::string extra_info="")
	catches all possible exceptions and rethrows them as cms::Exceptions More...
void	physicalOrder (uint16_t &readout_order, uint16_t &physical_order)
	order of strips More...
	RawToDigiUnpacker ()
	private default constructor More...
void	readoutOrder (uint16_t &physical_order, uint16_t &readout_order)
	order of strips More...
void	update (RawDigis &scope_mode, RawDigis &virgin_raw, RawDigis &proc_raw, Digis &zero_suppr, RawDigis &common_mode)
	fill DetSetVectors using registries More...
void	updateEventSummary (const sistrip::FEDBuffer &, SiStripEventSummary &)
	sets the SiStripEventSummary -> not yet implemented for FEDBuffer class More...

Static Private Member Functions
static void	dumpRawData (uint16_t fed_id, const FEDRawData &, std::stringstream &)
	dumps raw data to stdout (NB: payload is byte-swapped,headers/trailer are not). More...

Private Attributes
std::vector< SiStripRawDigi >	cm_work_digis_
std::vector< Registry >	cm_work_registry_
bool	doAPVEmulatorCheck_
bool	doFullCorruptBufferChecks_
uint32_t	errorThreshold_
uint32_t	event_
	other values More...
bool	extractCm_
int16_t	fedBufferDumpFreq_
int16_t	fedEventDumpFreq_
bool	first_
int16_t	headerBytes_
	configurables More...
bool	markMissingFeds_
bool	once_
std::vector< SiStripRawDigi >	proc_work_digis_
std::vector< Registry >	proc_work_registry_
bool	quiet_
std::vector< SiStripRawDigi >	scope_work_digis_
std::vector< Registry >	scope_work_registry_
int16_t	triggerFedId_
bool	unpackBadChannels_
bool	useDaqRegister_
bool	useFedKey_
std::vector< SiStripRawDigi >	virgin_work_digis_
std::vector< Registry >	virgin_work_registry_
std::vector< SiStripDigi >	zs_work_digis_
	digi collections More...
std::vector< Registry >	zs_work_registry_
	registries More...

Friends
class	RawToClustersLazyUnpacker

Detailed Description

Definition at line 29 of file SiStripRawToDigiUnpacker.h.

Member Typedef Documentation

typedef edm::DetSetVector<SiStripDigi> sistrip::RawToDigiUnpacker::Digis

Definition at line 35 of file SiStripRawToDigiUnpacker.h.

typedef edm::DetSetVector<SiStripRawDigi> sistrip::RawToDigiUnpacker::RawDigis

Definition at line 36 of file SiStripRawToDigiUnpacker.h.

Constructor & Destructor Documentation

sistrip::RawToDigiUnpacker::RawToDigiUnpacker	(	int16_t	appended_bytes,
		int16_t	fed_buffer_dump_freq,
		int16_t	fed_event_dump_freq,
		int16_t	trigger_fed_id,
		bool	using_fed_key,
		bool	unpack_bad_channels,
		bool	mark_missing_feds,
		const uint32_t	errorThreshold
	)

constructor

Definition at line 22 of file SiStripRawToDigiUnpacker.cc.

References edm::isDebugEnabled(), LogTrace, and unpackBadChannels_.

                                                                                                                                              :
    headerBytes_( appended_bytes ),
    fedBufferDumpFreq_( fed_buffer_dump_freq ),
    fedEventDumpFreq_( fed_event_dump_freq ),
    triggerFedId_( trigger_fed_id ),
    useFedKey_( using_fed_key ),
    unpackBadChannels_( unpack_bad_channels ),
    markMissingFeds_( mark_missing_feds ),
    event_(0),
    once_(true),
    first_(true),
    useDaqRegister_(false),
    quiet_(true),
    extractCm_(false),
    doFullCorruptBufferChecks_(false),
    doAPVEmulatorCheck_(true),
    errorThreshold_(errorThreshold)
  {
    if ( edm::isDebugEnabled() ) {
      LogTrace("SiStripRawToDigi")
    << "[sistrip::RawToDigiUnpacker::"<<__func__<<"]"
    <<" Constructing object...";
    }
    if (unpackBadChannels_) {
      edm::LogWarning("SiStripRawToDigi") << "Warning: Unpacking of bad channels enabled. Only enable this if you know what you are doing. " << std::endl;
    }
  }

sistrip::RawToDigiUnpacker::~RawToDigiUnpacker

(

)

default constructor

Definition at line 51 of file SiStripRawToDigiUnpacker.cc.

References edm::isDebugEnabled(), and LogTrace.

                                        {
    if ( edm::isDebugEnabled() ) {
      LogTrace("SiStripRawToDigi")
    << "[sistrip::RawToDigiUnpacker::"<<__func__<<"]"
    << " Destructing object...";
    }
  }

sistrip::RawToDigiUnpacker::RawToDigiUnpacker ( )

private

private default constructor

Member Function Documentation

void sistrip::RawToDigiUnpacker::cleanupWorkVectors ( )

private

method to clear registries and digi collections

Definition at line 746 of file SiStripRawToDigiUnpacker.cc.

References cm_work_digis_, cm_work_registry_, proc_work_digis_, proc_work_registry_, scope_work_digis_, scope_work_registry_, virgin_work_digis_, virgin_work_registry_, zs_work_digis_, and zs_work_registry_.

Referenced by createDigis().

                                             {
    // Clear working areas and registries
    zs_work_registry_.clear();      zs_work_digis_.clear();
    virgin_work_registry_.clear();  virgin_work_digis_.clear();
    proc_work_registry_.clear();    proc_work_digis_.clear();
    scope_work_registry_.clear();   scope_work_digis_.clear();
    cm_work_registry_.clear();      cm_work_digis_.clear();
  }

void sistrip::RawToDigiUnpacker::createDigis	(	const SiStripFedCabling &	cabling,
		const FEDRawDataCollection &	buffers,
		SiStripEventSummary &	summary,
		RawDigis &	scope_mode,
		RawDigis &	virgin_raw,
		RawDigis &	proc_raw,
		Digis &	zero_suppr,
		DetIdCollection &	detids,
		RawDigis &	common_mode
	)

creates digis

extract readout mode

FED channel

create unpacker

unpack -> add check to make sure strip < nstrips && strip > last strip......

create unpacker

unpack -> add check to make sure strip < nstrips && strip > last strip......

create unpacker

unpack -> add check to make sure strip < nstrips && strip > last strip......

create unpacker

unpack -> add check to make sure strip < nstrips && strip > last strip......

create unpacker

unpack -> add check to make sure strip < nstrips && strip > last strip......

create unpacker

unpack -> add check to make sure strip < nstrips && strip > last strip......

Definition at line 59 of file SiStripRawToDigiUnpacker.cc.

References sistrip::FEDZSChannelUnpacker::adc(), sistrip::FEDRawChannelUnpacker::adc(), sistrip::APV_LATENCY, cleanupWorkVectors(), cm_work_digis_, cm_work_registry_, FEDRawData::data(), doAPVEmulatorCheck_, doFullCorruptBufferChecks_, dumpRawData(), alignCSCRings::e, errorThreshold_, event_, edm::hlt::Exception, extractCm_, sistrip::FED_CABLING, fedBufferDumpFreq_, sistrip::FEDCH_PER_FEUNIT, SiStripFedCabling::fedConnections(), FEDRawDataCollection::FEDData(), fedEventDumpFreq_, SiStripFedCabling::fedIds(), sistrip::FINE_DELAY, first_, sistrip::FEDZSChannelUnpacker::hasData(), sistrip::FEDRawChannelUnpacker::hasData(), i, input, sistrip::invalid32_, sistrip::invalid_, edm::isDebugEnabled(), SiStripEventSummary::isSet(), combine::key, LogTrace, FEDNumbering::MAXSiStripFEDID, FEDNumbering::MINSiStripFEDID, sistrip::mlRawToDigi_, alignBH_cfg::mode, n, proc_work_digis_, proc_work_registry_, sistrip::FEDRawChannelUnpacker::procRawModeUnpacker(), edm::EDCollection< T >::push_back(), quiet_, sistrip::READOUT_MODE_PROC_RAW, sistrip::READOUT_MODE_SCOPE, sistrip::READOUT_MODE_VIRGIN_RAW, sistrip::READOUT_MODE_ZERO_SUPPRESSED, sistrip::READOUT_MODE_ZERO_SUPPRESSED_LITE, readoutOrder(), edm::EDCollection< T >::reserve(), SiStripEventSummary::runType(), sistrip::FEDZSChannelUnpacker::sampleNumber(), scope_work_digis_, scope_work_registry_, sistrip::FEDRawChannelUnpacker::scopeModeUnpacker(), edm::EDCollection< T >::size(), FEDRawData::size(), contentValuesCheck::ss, triggerFedId_, unpackBadChannels_, update(), updateEventSummary(), useDaqRegister_, useFedKey_, SiStripEventSummary::valid(), virgin_work_digis_, virgin_work_registry_, sistrip::FEDRawChannelUnpacker::virginRawModeUnpacker(), cms::Exception::what(), sistrip::FEDZSChannelUnpacker::zeroSuppressedLiteModeUnpacker(), sistrip::FEDZSChannelUnpacker::zeroSuppressedModeUnpacker(), zs_work_digis_, and zs_work_registry_.

Referenced by sistrip::RawToDigiModule::produce().

                                                                                                                                                                                                                                                                            {

    // Clear working areas and registries
    cleanupWorkVectors();
    // Reserve space in bad module list
    detids.reserve(100);
  
    // Check if FEDs found in cabling map and event data
    if ( edm::isDebugEnabled() ) {
      if ( cabling.fedIds().empty() ) {
    edm::LogWarning(sistrip::mlRawToDigi_)
      << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
      << " No FEDs found in cabling map!";
    // Check which FED ids have non-zero size buffers
    std::vector<uint16_t> feds;
    for ( uint16_t ifed = FEDNumbering::MINSiStripFEDID; ifed < FEDNumbering::MAXSiStripFEDID; ifed++ ) {
      if ( ifed != triggerFedId_ && 
           buffers.FEDData( static_cast<int>(ifed) ).size() ) {
        feds.push_back(ifed);
      }
    }
    LogTrace("SiStripRawToDigi")
      << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
      << " Found " 
      << feds.size() 
      << " FED buffers with non-zero size!";
      }
    }

    // Flag for EventSummary update using DAQ register  
    bool first_fed = true;
  
    // Retrieve FED ids from cabling map and iterate through 
    std::vector<uint16_t>::const_iterator ifed = cabling.fedIds().begin();
    for ( ; ifed != cabling.fedIds().end(); ifed++ ) {

      // ignore trigger FED
      if ( *ifed == triggerFedId_ ) { continue;  }
    
      // Retrieve FED raw data for given FED 
      const FEDRawData& input = buffers.FEDData( static_cast<int>(*ifed) );
    
      // Some debug on FED buffer size
      if ( edm::isDebugEnabled() ) {
    if ( first_ && input.data() ) {
      std::stringstream ss;
      ss << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
         << " Found FED id " 
         << std::setw(4) << std::setfill(' ') << *ifed 
         << " in FEDRawDataCollection"
         << " with non-zero pointer 0x" 
         << std::hex
         << std::setw(8) << std::setfill('0') 
         << reinterpret_cast<uint32_t*>( const_cast<uint8_t*>(input.data()))
         << std::dec
         << " and size " 
         << std::setw(5) << std::setfill(' ') << input.size()
         << " chars";
      LogTrace("SiStripRawToDigi") << ss.str();
    }   
      }
    
      // Dump of FEDRawData to stdout
      if ( edm::isDebugEnabled() ) {
    if ( fedBufferDumpFreq_ && !(event_%fedBufferDumpFreq_) ) {
      std::stringstream ss;
      dumpRawData( *ifed, input, ss );
      edm::LogVerbatim(sistrip::mlRawToDigi_) << ss.str();
    }
      }
      
      // get the cabling connections for this FED
      auto conns = cabling.fedConnections(*ifed);
    
      // Check on FEDRawData pointer
      if ( !input.data() ) {
    if ( edm::isDebugEnabled() ) {
      edm::LogWarning(sistrip::mlRawToDigi_)
        << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
        << " NULL pointer to FEDRawData for FED id " 
        << *ifed;
    }
        // Mark FED modules as bad
        detids.reserve(detids.size()+conns.size());
        std::vector<FedChannelConnection>::const_iterator iconn = conns.begin();
        for ( ; iconn != conns.end(); iconn++ ) {
          if ( !iconn->detId() || iconn->detId() == sistrip::invalid32_ ) continue;
          detids.push_back(iconn->detId()); //@@ Possible multiple entries (ok for Giovanni)
        }
    continue;
      } 
    
      // Check on FEDRawData size
      if ( !input.size() ) {
    if ( edm::isDebugEnabled() ) {
      edm::LogWarning(sistrip::mlRawToDigi_)
        << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
        << " FEDRawData has zero size for FED id " 
        << *ifed;
    }
        // Mark FED modules as bad
        detids.reserve(detids.size()+conns.size());
        std::vector<FedChannelConnection>::const_iterator iconn = conns.begin();
        for ( ; iconn != conns.end(); iconn++ ) {
          if ( !iconn->detId() || iconn->detId() == sistrip::invalid32_ ) continue;
          detids.push_back(iconn->detId()); //@@ Possible multiple entries (ok for Giovanni)
        }
        continue;
      }
      
      // construct FEDBuffer
      std::auto_ptr<sistrip::FEDBuffer> buffer;
      try {
        buffer.reset(new sistrip::FEDBuffer(input.data(),input.size()));
        if (!buffer->doChecks()) {
          if (!unpackBadChannels_ || !buffer->checkNoFEOverflows() )
            throw cms::Exception("FEDBuffer") << "FED Buffer check fails for FED ID " << *ifed << ".";
        }
        if (doFullCorruptBufferChecks_ && !buffer->doCorruptBufferChecks()) {
          throw cms::Exception("FEDBuffer") << "FED corrupt buffer check fails for FED ID " << *ifed << ".";
        }
      }
      catch (const cms::Exception& e) { 
    if ( edm::isDebugEnabled() ) {
      edm::LogWarning("sistrip::RawToDigiUnpacker") << "Exception caught when creating FEDBuffer object for FED " << *ifed << ": " << e.what();
    }
        // FED buffer is bad and should not be unpacked. Skip this FED and mark all modules as bad. 
        std::vector<FedChannelConnection>::const_iterator iconn = conns.begin();
        for ( ; iconn != conns.end(); iconn++ ) {
          if ( !iconn->detId() || iconn->detId() == sistrip::invalid32_ ) continue;
          detids.push_back(iconn->detId()); //@@ Possible multiple entries (ok for Giovanni)
        }
        continue;
      }

      // Check if EventSummary ("trigger FED info") needs updating
      if ( first_fed && useDaqRegister_ ) { updateEventSummary( *buffer, summary ); first_fed = false; }
    
      // Check to see if EventSummary info is set
      if ( edm::isDebugEnabled() ) {
    if ( !quiet_ && !summary.isSet() ) {
      std::stringstream ss;
      ss << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
         << " EventSummary is not set correctly!"
         << " Missing information from both \"trigger FED\" and \"DAQ registers\"!";
      edm::LogWarning(sistrip::mlRawToDigi_) << ss.str();
    }
      }
    
      // Check to see if event is to be analyzed according to EventSummary
      if ( !summary.valid() ) { 
    if ( edm::isDebugEnabled() ) {
      LogTrace("SiStripRawToDigi")
        << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
        << " EventSummary is not valid: skipping...";
    }
    continue; 
      }
    
      sistrip::FEDReadoutMode mode = buffer->readoutMode(); 

      // Retrive run type
      sistrip::RunType runType_ = summary.runType();
      if( runType_ == sistrip::APV_LATENCY || runType_ == sistrip::FINE_DELAY ) { useFedKey_ = false; } 
     
      // Dump of FED buffer
      if ( edm::isDebugEnabled() ) {
    if ( fedEventDumpFreq_ && !(event_%fedEventDumpFreq_) ) {
      std::stringstream ss;
      buffer->dump( ss );
      edm::LogVerbatim(sistrip::mlRawToDigi_) << ss.str();
    }
      }
    
      // Iterate through FED channels, extract payload and create Digis
      std::vector<FedChannelConnection>::const_iterator iconn = conns.begin();
      for ( ; iconn != conns.end(); iconn++ ) {

    uint16_t chan = iconn->fedCh();

    // Check if fed connection is valid
    if ( !iconn->isConnected() ) { continue; }
        
        // Check DetId is valid (if to be used as key)
    if ( !useFedKey_ && ( !iconn->detId() || iconn->detId() == sistrip::invalid32_ ) ) { continue; }
      
    // Check FED channel
    if (!buffer->channelGood(iconn->fedCh(),doAPVEmulatorCheck_)) {
          if (!unpackBadChannels_ || !(buffer->fePresent(iconn->fedCh()/FEDCH_PER_FEUNIT) && buffer->feEnabled(iconn->fedCh()/FEDCH_PER_FEUNIT)) ) {
            detids.push_back(iconn->detId()); //@@ Possible multiple entries (ok for Giovanni)
            continue;
          }
    }

    // Determine whether FED key is inferred from cabling or channel loop
    uint32_t fed_key = ( summary.runType() == sistrip::FED_CABLING ) ? ( ( *ifed & sistrip::invalid_ ) << 16 ) | ( chan & sistrip::invalid_ ) : ( ( iconn->fedId() & sistrip::invalid_ ) << 16 ) | ( iconn->fedCh() & sistrip::invalid_ );

    // Determine whether DetId or FED key should be used to index digi containers
    uint32_t key = ( useFedKey_ || mode == sistrip::READOUT_MODE_SCOPE ) ? fed_key : iconn->detId();
      
    // Determine APV std::pair number (needed only when using DetId)
    uint16_t ipair = ( useFedKey_ || mode == sistrip::READOUT_MODE_SCOPE ) ? 0 : iconn->apvPairNumber();
      

    if (mode == sistrip::READOUT_MODE_ZERO_SUPPRESSED ) { 
    
      Registry regItem(key, 0, zs_work_digis_.size(), 0);
    
          try {
        sistrip::FEDZSChannelUnpacker unpacker = sistrip::FEDZSChannelUnpacker::zeroSuppressedModeUnpacker(buffer->channel(iconn->fedCh()));
        
            
        while (unpacker.hasData()) {zs_work_digis_.push_back(SiStripDigi(unpacker.sampleNumber()+ipair*256,unpacker.adc()));unpacker++;}
          } catch (const cms::Exception& e) {
            if ( edm::isDebugEnabled() ) {
              edm::LogWarning(sistrip::mlRawToDigi_)
                << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
                << " Clusters are not ordered for FED "
                << *ifed << " channel " << iconn->fedCh()
                << ": " << e.what();
            }
            detids.push_back(iconn->detId()); //@@ Possible multiple entries (ok for Giovanni)
            continue;
          }
          
      regItem.length = zs_work_digis_.size() - regItem.index;
      if (regItem.length > 0) {
        regItem.first = zs_work_digis_[regItem.index].strip();
        zs_work_registry_.push_back(regItem);
      }

        
      // Common mode values
      if ( extractCm_ ) {
        try {
          Registry regItem2( key, 2*ipair, cm_work_digis_.size(), 2 );
          cm_work_digis_.push_back( SiStripRawDigi( buffer->channel(iconn->fedCh()).cmMedian(0) ) );
          cm_work_digis_.push_back( SiStripRawDigi( buffer->channel(iconn->fedCh()).cmMedian(1) ) );
          cm_work_registry_.push_back( regItem2 );
        } catch (const cms::Exception& e) {
          if ( edm::isDebugEnabled() ) {
        edm::LogWarning(sistrip::mlRawToDigi_)
          << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
          << " Problem extracting common modes for FED id "
          << *ifed << " and channel " << iconn->fedCh()
          << ": " << std::endl << e.what();
          }
        }
      }
      
    }

    else if (mode == sistrip::READOUT_MODE_ZERO_SUPPRESSED_LITE ) { 

      Registry regItem(key, 0, zs_work_digis_.size(), 0);
    
      try {
        sistrip::FEDZSChannelUnpacker unpacker = sistrip::FEDZSChannelUnpacker::zeroSuppressedLiteModeUnpacker(buffer->channel(iconn->fedCh()));
        
        while (unpacker.hasData()) {zs_work_digis_.push_back(SiStripDigi(unpacker.sampleNumber()+ipair*256,unpacker.adc()));unpacker++;}
      } catch (const cms::Exception& e) {
            if ( edm::isDebugEnabled() ) {
              edm::LogWarning(sistrip::mlRawToDigi_)
                << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
                << " Clusters are not ordered for FED "
                << *ifed << " channel " << iconn->fedCh()
                << ": " << e.what();
            }
            detids.push_back(iconn->detId()); //@@ Possible multiple entries (ok for Giovanni)
            continue;
          }  

      regItem.length = zs_work_digis_.size() - regItem.index;
      if (regItem.length > 0) {
        regItem.first = zs_work_digis_[regItem.index].strip();
        zs_work_registry_.push_back(regItem);
      }
          

    } 
     
    else if ( mode == sistrip::READOUT_MODE_VIRGIN_RAW ) {

      std::vector<uint16_t> samples; 

      sistrip::FEDRawChannelUnpacker unpacker = sistrip::FEDRawChannelUnpacker::virginRawModeUnpacker(buffer->channel(iconn->fedCh()));

      while (unpacker.hasData()) {samples.push_back(unpacker.adc());unpacker++;}

      if ( !samples.empty() ) { 
        Registry regItem(key, 256*ipair, virgin_work_digis_.size(), samples.size());
        uint16_t physical;
        uint16_t readout; 
        for ( uint16_t i = 0, n = samples.size(); i < n; i++ ) {
          physical = i%128;
          readoutOrder( physical, readout );                 // convert index from physical to readout order
          (i/128) ? readout=readout*2+1 : readout=readout*2; // un-multiplex data
          virgin_work_digis_.push_back(  SiStripRawDigi( samples[readout] ) );
        }
        virgin_work_registry_.push_back( regItem );
      }
    } 
    
    else if ( mode == sistrip::READOUT_MODE_PROC_RAW ) {
    
      std::vector<uint16_t> samples; 
    
      sistrip::FEDRawChannelUnpacker unpacker = sistrip::FEDRawChannelUnpacker::procRawModeUnpacker(buffer->channel(iconn->fedCh()));
    
      while (unpacker.hasData()) {samples.push_back(unpacker.adc());unpacker++;}
    
      if ( !samples.empty() ) { 
        Registry regItem(key, 256*ipair, proc_work_digis_.size(), samples.size());
        for ( uint16_t i = 0, n = samples.size(); i < n; i++ ) {
          proc_work_digis_.push_back(  SiStripRawDigi( samples[i] ) );
        }
        proc_work_registry_.push_back( regItem );
      }
    } 

    else if ( mode == sistrip::READOUT_MODE_SCOPE ) {
    
      std::vector<uint16_t> samples; 
    
      sistrip::FEDRawChannelUnpacker unpacker = sistrip::FEDRawChannelUnpacker::scopeModeUnpacker(buffer->channel(iconn->fedCh()));

      while (unpacker.hasData()) {samples.push_back(unpacker.adc());unpacker++;}
    
      if ( !samples.empty() ) { 
        Registry regItem(key, 0, scope_work_digis_.size(), samples.size());
        for ( uint16_t i = 0, n = samples.size(); i < n; i++ ) {
          scope_work_digis_.push_back(  SiStripRawDigi( samples[i] ) );
        }
        scope_work_registry_.push_back( regItem );
      }
    } 
    
    else { // Unknown readout mode! => assume scope mode

      if ( edm::isDebugEnabled() ) {
        std::stringstream ss;
        ss << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
           << " Unknown FED readout mode (" << mode
           << ")! Assuming SCOPE MODE..."; 
        edm::LogWarning(sistrip::mlRawToDigi_) << ss.str();
      }
    
      std::vector<uint16_t> samples; 
    
      sistrip::FEDRawChannelUnpacker unpacker = sistrip::FEDRawChannelUnpacker::scopeModeUnpacker(buffer->channel(iconn->fedCh()));
    
      while (unpacker.hasData()) {samples.push_back(unpacker.adc());unpacker++;}
    
      if ( !samples.empty() ) { 
        Registry regItem(key, 0, scope_work_digis_.size(), samples.size());
        for ( uint16_t i = 0, n = samples.size(); i < n; i++ ) {
          scope_work_digis_.push_back(  SiStripRawDigi( samples[i] ) );
        }
        scope_work_registry_.push_back( regItem );
      
        if ( edm::isDebugEnabled() ) {
          std::stringstream ss;
          ss << "Extracted " << samples.size() 
         << " SCOPE MODE digis (samples[0] = " 
         << samples[0] 
         << ") from FED id/ch " 
         << iconn->fedId() 
         << "/" 
         << iconn->fedCh();
          LogTrace("SiStripRawToDigi") << ss.str();
        }
      }
      else if ( edm::isDebugEnabled() ) {
        edm::LogWarning(sistrip::mlRawToDigi_)
          << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
          << " No SM digis found!"; 
      }
    } 
      } // channel loop
    } // fed loop

    // bad channels warning
    unsigned int detIdsSize = detids.size();
    if ( edm::isDebugEnabled() && detIdsSize ) {
      std::ostringstream ss;
      ss << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
         << " Problems were found in data and " << detIdsSize << " channels could not be unpacked. "
         << "See output of FED Hardware monitoring for more information. ";
      edm::LogWarning(sistrip::mlRawToDigi_) << ss.str();
    }
    if( (errorThreshold_ != 0) && (detIdsSize > errorThreshold_) ) {
      edm::LogError("TooManyErrors") << "Total number of errors = " << detIdsSize;
    }

    // update DetSetVectors
    update(scope_mode, virgin_raw, proc_raw, zero_suppr, cm_values);

    // increment event counter
    event_++;
  
    // no longer first event!
    if ( first_ ) { first_ = false; }
  
    // final cleanup, just in case
    cleanupWorkVectors();
  }

void sistrip::RawToDigiUnpacker::doAPVEmulatorCheck ( bool  do_APVEmulator_check )

inline

Definition at line 186 of file SiStripRawToDigiUnpacker.h.

{ doAPVEmulatorCheck_ = do_APVEmulator_check; }

void sistrip::RawToDigiUnpacker::doFullCorruptBufferChecks ( bool  do_full_corrupt_buffer_checks )

inline

Definition at line 184 of file SiStripRawToDigiUnpacker.h.

{ doFullCorruptBufferChecks_ = do_full_corrupt_buffer_checks; }

void sistrip::RawToDigiUnpacker::dumpRawData ( uint16_t  fed_id, const FEDRawData &  buffer, std::stringstream &  ss  )

staticprivate

dumps raw data to stdout (NB: payload is byte-swapped,headers/trailer are not).

Definition at line 1051 of file SiStripRawToDigiUnpacker.cc.

References FEDRawData::data(), relativeConstraints::empty, i, sistrip::invalid32_, and FEDRawData::size().

Referenced by createDigis().

                                                                                                      {

    ss << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
       << " Dump of buffer for FED id " <<  fed_id << std::endl
       << " Buffer contains " << buffer.size()
       << " bytes (NB: payload is byte-swapped)" << std::endl;
    uint32_t* buffer_u32 = reinterpret_cast<uint32_t*>( const_cast<unsigned char*>( buffer.data() ) );
    unsigned int empty = 0;

    if ( 0 ) { 

      ss << "Byte->   4 5 6 7 0 1 2 3\n";
      for ( uint32_t i = 0; i < buffer.size()/8; i++ ) {
    unsigned int temp0 = buffer_u32[i*2] & sistrip::invalid32_;
    unsigned int temp1 = buffer_u32[i*2+1] & sistrip::invalid32_;
    if ( !temp0 && !temp1 ) { empty++; }
    else { 
      if ( empty ) { 
        ss << "        [ empty  words ]" << std::endl; 
        empty = 0; 
      }
      ss << std::dec
         << std::setfill(' ')  << std::setw(6) << i*8 << ": " 
         << std::hex 
         << std::setfill('0') << std::setw(8) << temp0 
         << std::setfill('0') << std::setw(8) << temp1 
         << std::dec
         << std::endl;
    }
      }

    } else {
    
      ss << "  Byte |  <---- Byte order ----<  | Byte" << std::endl;
      ss << "  cntr |  7  6  5  4  3  2  1  0  | cntr" << std::endl;
      for ( uint32_t i = 0; i < buffer.size()/8; i++ ) {
    //if ( i>=20 && ((i+4)<(buffer.size()/8)) ) { continue; }
    uint16_t tmp0 = buffer.data()[i*8+0] & 0xFF;
    uint16_t tmp1 = buffer.data()[i*8+1] & 0xFF;
    uint16_t tmp2 = buffer.data()[i*8+2] & 0xFF;
    uint16_t tmp3 = buffer.data()[i*8+3] & 0xFF;
    uint16_t tmp4 = buffer.data()[i*8+4] & 0xFF;
    uint16_t tmp5 = buffer.data()[i*8+5] & 0xFF;
    uint16_t tmp6 = buffer.data()[i*8+6] & 0xFF;
    uint16_t tmp7 = buffer.data()[i*8+7] & 0xFF;
//  if ( !tmp0 && !tmp1 && !tmp2 && !tmp3 &&
//       !tmp4 && !tmp5 && !tmp6 && !tmp7 ) { empty++; }
//  else { 
//    if ( empty ) { 
//      ss << "         [.." 
//         << std::dec << std::setfill('.') << std::setw(4) << empty 
//         << " null words....]" << std::endl; 
//      empty = 0; 
//    }
      ss << std::dec
         << std::setfill(' ')  << std::setw(6) << i*8+7 << " : " 
         << std::hex 
         << std::setfill('0') << std::setw(2) << tmp7 << " " 
         << std::setfill('0') << std::setw(2) << tmp6 << " " 
         << std::setfill('0') << std::setw(2) << tmp5 << " " 
         << std::setfill('0') << std::setw(2) << tmp4 << " " 
         << std::setfill('0') << std::setw(2) << tmp3 << " " 
         << std::setfill('0') << std::setw(2) << tmp2 << " " 
         << std::setfill('0') << std::setw(2) << tmp1 << " " 
         << std::setfill('0') << std::setw(2) << tmp0 
         << std::dec
         << " :" << std::setfill(' ')  << std::setw(6) << i*8 
         << std::endl;
//  }
      }

    }
    ss << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
       << " End of FED buffer";
  }

void sistrip::RawToDigiUnpacker::extractCm ( bool  extract_cm )

inline

Definition at line 182 of file SiStripRawToDigiUnpacker.h.

{ extractCm_ = extract_cm; }

sistrip::FedBufferFormat sistrip::RawToDigiUnpacker::fedBufferFormat ( const uint16_t &  register_value )

inlineprivate

returns buffer format

Definition at line 160 of file SiStripRawToDigiUnpacker.h.

References sistrip::APV_ERROR_FORMAT, sistrip::FULL_DEBUG_FORMAT, sistrip::UNDEFINED_FED_BUFFER_FORMAT, and sistrip::UNKNOWN_FED_BUFFER_FORMAT.

{
  if ( (register_value&0xF) == 0x1 ) { return sistrip::FULL_DEBUG_FORMAT; }
  else if ( (register_value&0xF) == 0x2 ) { return sistrip::APV_ERROR_FORMAT; }
  else if ( (register_value&0xF) == 0x0 ) { return sistrip::UNDEFINED_FED_BUFFER_FORMAT; }
  else { return sistrip::UNKNOWN_FED_BUFFER_FORMAT; }
}

sistrip::FedReadoutMode sistrip::RawToDigiUnpacker::fedReadoutMode ( const uint16_t &  register_value )

inlineprivate

returns buffer readout mode

Definition at line 168 of file SiStripRawToDigiUnpacker.h.

References sistrip::FED_PROC_RAW, sistrip::FED_SCOPE_MODE, sistrip::FED_VIRGIN_RAW, sistrip::FED_ZERO_SUPPR, sistrip::FED_ZERO_SUPPR_LITE, and sistrip::UNKNOWN_FED_READOUT_MODE.

{
  if ( ((register_value>>1)&0x7) == 0x0 ) { return sistrip::FED_SCOPE_MODE; }
  else if ( ((register_value>>1)&0x7) == 0x1 ) { return sistrip::FED_VIRGIN_RAW; }
  else if ( ((register_value>>1)&0x7) == 0x3 ) { return sistrip::FED_PROC_RAW; }
  else if ( ((register_value>>1)&0x7) == 0x5 ) { return sistrip::FED_ZERO_SUPPR; }
  else if ( ((register_value>>1)&0x7) == 0x6 ) { return sistrip::FED_ZERO_SUPPR_LITE; }
  else { return sistrip::UNKNOWN_FED_READOUT_MODE; }
}

void sistrip::RawToDigiUnpacker::handleException ( std::string  method_name, std::string  extra_info = ""  )

private

catches all possible exceptions and rethrows them as cms::Exceptions

Definition at line 1127 of file SiStripRawToDigiUnpacker.cc.

References alignCSCRings::e, cppFunctionSkipper::exception, edm::isDebugEnabled(), sistrip::mlRawToDigi_, and contentValuesCheck::ss.

                                                                                       { 

    method_name = "sistrip::RawToDigiUnpacker::" + method_name;
    try {
      throw; // rethrow caught exception to be dealt with below
    } 
    catch ( const cms::Exception& e ) { 
      //throw e; // rethrow cms::Exception to be caught by framework
    }
    catch ( const std::exception& e ) {
      if ( edm::isDebugEnabled() ) {
    std::stringstream ss;
    ss << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
       << " Caught std::exception!" << std::endl;
    if ( extra_info != "" ) { 
      ss << " Information: " << extra_info << std::endl;
    }
    ss << " Caught std::exception in ["
       << method_name << "] with message:" << std::endl 
       << e.what();
    edm::LogWarning(sistrip::mlRawToDigi_) << ss.str();
      }
      //throw cms::Exception(sistrip::mlRawToDigi_) << ss.str();
    }
    catch (...) {
      if ( edm::isDebugEnabled() ) {
    std::stringstream ss;
    ss << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
       << " Caught unknown exception!" << std::endl;
    if ( extra_info != "" ) { 
      ss << " Information: " << extra_info << std::endl;
    }
    ss << "Caught unknown exception in ["
       << method_name << "]" << std::endl;
    edm::LogWarning(sistrip::mlRawToDigi_) << ss.str();
      }
      //throw cms::Exception(sistrip::mlRawToDigi_) << ss.str();
    }
  }

void sistrip::RawToDigiUnpacker::locateStartOfFedBuffer	(	const uint16_t &	fed_id,
		const FEDRawData &	input,
		FEDRawData &	output
	)

Removes any data appended prior to FED buffer and reorders 32-bit words if swapped.

Definition at line 873 of file SiStripRawToDigiUnpacker.cc.

References FEDRawData::data(), newFWLiteAna::found, headerBytes_, edm::isDebugEnabled(), getDQMSummary::iter, LogTrace, sistrip::mlRawToDigi_, evf::evtn::offset(), FEDRawData::resize(), FEDRawData::size(), and contentValuesCheck::ss.

                                                                                                                      {
  
    // Check size of input buffer
    if ( input.size() < 24 ) { 
      output.resize( input.size() ); // Return UNadjusted buffer start position and size
      memcpy( output.data(), input.data(), input.size() );
      if ( edm::isDebugEnabled() ) {
    std::stringstream ss; 
    ss << "[sistrip::RawToDigiUnpacker::" << __func__ << "] "
       << "Input FEDRawData with FED id " << fed_id 
       << " has size " << input.size();
    edm::LogWarning(sistrip::mlRawToDigi_) << ss.str();
      }
      return;
    } 
  
    // Iterator through buffer to find DAQ header 
    bool found = false;
    uint16_t ichar = 0;
    while ( ichar < input.size()-16 && !found ) {
      uint16_t offset = headerBytes_ < 0 ? ichar : headerBytes_; // Negative value means use "search mode" to find DAQ header
      uint32_t* input_u32   = reinterpret_cast<uint32_t*>( const_cast<unsigned char*>( input.data() ) + offset );
      uint32_t* fed_trailer = reinterpret_cast<uint32_t*>( const_cast<unsigned char*>( input.data() ) + input.size() - 8 );
    
      // see info on FED 32-bit swapping at end-of-file

      bool old_vme_header = ( input_u32[0] & 0xF0000000 ) == 0x50000000 && ( fed_trailer[0]  & 0xF0000000 ) == 0xA0000000 && ( (fed_trailer[0] & 0x00FFFFFF)*0x8 ) == (input.size() - offset);
    
      bool old_slink_header = ( input_u32[1] & 0xF0000000 ) == 0x50000000 && ( fed_trailer[1]  & 0xF0000000 ) == 0xA0000000 && ( (fed_trailer[1] & 0x00FFFFFF)*0x8 ) == (input.size() - offset);
    
      bool old_slink_payload = ( input_u32[3] & 0xFF000000 ) == 0xED000000;
    
      bool new_buffer_format = ( input_u32[2] & 0xFF000000 ) == 0xC5000000;
    
      if ( old_vme_header )  {
      
    // Found DAQ header at byte position 'offset'
    found = true;
    output.resize( input.size()-offset );
    memcpy( output.data(),         // target
        input.data()+offset,   // source
        input.size()-offset ); // nbytes
    if ( headerBytes_ < 0 ) {
      if ( edm::isDebugEnabled() ) {
        std::stringstream ss;
        ss << "[sistrip::RawToDigiUnpacker::" << __func__ << "]" 
           << " Buffer for FED id " << fed_id 
           << " has been found at byte position " << offset
           << " with a size of " << input.size()-offset << " bytes."
           << " Adjust the configurable 'AppendedBytes' to " << offset;
        LogTrace("SiStripRawToDigi") << ss.str();
      }
    }
      
      } else if ( old_slink_header ) {
      
    if ( old_slink_payload ) {
      
      // Found DAQ header (with MSB and LSB 32-bit words swapped) at byte position 'offset' 
      found = true;
      output.resize( input.size()-offset );
      uint32_t* output_u32 = reinterpret_cast<uint32_t*>( const_cast<unsigned char*>( output.data() ) );
      uint16_t iter = offset; 
      while ( iter < output.size() / sizeof(uint32_t) ) {
        output_u32[iter] = input_u32[iter+1];
        output_u32[iter+1] = input_u32[iter];
        iter+=2;
      }
      if ( headerBytes_ < 0 ) {
        if ( edm::isDebugEnabled() ) {
          std::stringstream ss;
          ss << "[sistrip::RawToDigiUnpacker::" << __func__ << "]" 
         << " Buffer (with MSB and LSB 32-bit words swapped) for FED id " << fed_id 
         << " has been found at byte position " << offset
         << " with a size of " << output.size() << " bytes."
         << " Adjust the configurable 'AppendedBytes' to " << offset;
          LogTrace("SiStripRawToDigi") << ss.str();
        }
      }

    } else if ( new_buffer_format ) {
    
      // Found DAQ header at byte position 'offset'
      found = true;
      output.resize( input.size()-offset );
      memcpy( output.data(),         // target
          input.data()+offset,   // source
          input.size()-offset ); // nbytes
      if ( headerBytes_ < 0 ) {
        if ( edm::isDebugEnabled() ) {
          std::stringstream ss;
          ss << "[sistrip::RawToDigiUnpacker::" << __func__ << "]" 
         << " Buffer for FED id " << fed_id 
         << " has been found at byte position " << offset
         << " with a size of " << input.size()-offset << " bytes."
         << " Adjust the configurable 'AppendedBytes' to " << offset;
          LogTrace("SiStripRawToDigi") << ss.str();
        }
      }
    
    } else { headerBytes_ < 0 ? found = false : found = true; }
      } else { headerBytes_ < 0 ? found = false : found = true; }
      ichar++;
    }      
  
    // Check size of output buffer
    if ( output.size() == 0 ) { 
    
      // Did not find DAQ header after search => return buffer with null size
      output.resize( 0 ); //@@ NULL SIZE
      memcpy( output.data(), input.data(), 0 ); //@@ NULL SIZE
      if ( edm::isDebugEnabled() ) {
    std::stringstream ss;
    if ( headerBytes_ < 0 ) {
      ss << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
         << " DAQ header not found within buffer for FED id: " << fed_id;
    } else {
      uint32_t* input_u32 = reinterpret_cast<uint32_t*>( const_cast<unsigned char*>( input.data() ) );
      ss << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
         << " DAQ header not found at expected location for FED id: " << fed_id << std::endl
         << " First 64-bit word of buffer is 0x"
         << std::hex 
         << std::setfill('0') << std::setw(8) << input_u32[0] 
         << std::setfill('0') << std::setw(8) << input_u32[1] 
         << std::dec << std::endl
         << " Adjust 'AppendedBytes' configurable to '-1' to activate 'search mode'";
    }
    edm::LogWarning(sistrip::mlRawToDigi_) << ss.str();
      }
    
    } else if ( output.size() < 24 ) { // Found DAQ header after search, but too few words
    
      if ( edm::isDebugEnabled() ) {
    std::stringstream ss; 
    ss << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
       << " Unexpected buffer size! FEDRawData with FED id " << fed_id 
       << " has size " << output.size();
    edm::LogWarning(sistrip::mlRawToDigi_) << ss.str();
      }
    }   
  }

void sistrip::RawToDigiUnpacker::physicalOrder ( uint16_t &  readout_order, uint16_t &  physical_order  )

inlineprivate

order of strips

Definition at line 155 of file SiStripRawToDigiUnpacker.h.

{
  physical_order = ( (32 * (readout_order%4)) + (8 * static_cast<uint16_t>(static_cast<float>(readout_order)/4.0)) - (31 * static_cast<uint16_t>(static_cast<float>(readout_order)/16.0)) );
}

void sistrip::RawToDigiUnpacker::quiet ( bool  quiet )

inline

verbosity

Definition at line 178 of file SiStripRawToDigiUnpacker.h.

{ quiet_ = quiet; }

void sistrip::RawToDigiUnpacker::readoutOrder ( uint16_t &  physical_order, uint16_t &  readout_order  )

inlineprivate

order of strips

Definition at line 150 of file SiStripRawToDigiUnpacker.h.

Referenced by createDigis().

{
  readout_order = ( 4*((static_cast<uint16_t>((static_cast<float>(physical_order)/8.0)))%4) + static_cast<uint16_t>(static_cast<float>(physical_order)/32.0) + 16*(physical_order%8) );
}

void sistrip::RawToDigiUnpacker::triggerFed	(	const FEDRawDataCollection &	buffers,
		SiStripEventSummary &	summary,
		const uint32_t &	event
	)

trigger info

Definition at line 755 of file SiStripRawToDigiUnpacker.cc.

References SiStripEventSummary::bx(), SiStripEventSummary::commissioningInfo(), fedt_struct::conscheck, FEDRawData::data(), SiStripEventSummary::event(), FEDRawDataCollection::FEDData(), TFHeaderDescription::getBunchCrossing(), TFHeaderDescription::getFedEventNumber(), edm::isDebugEnabled(), SiStripEventSummary::isSet(), FEDNumbering::lastFEDId(), LogTrace, sistrip::mlRawToDigi_, once_, FEDRawData::size(), contentValuesCheck::ss, edmLumisInFiles::summary, groupFilesInBlocks::temp, SiStripEventSummary::triggerFed(), and triggerFedId_.

Referenced by sistrip::RawToDigiModule::produce().

                                                                                                                               {
  
    // Pointer to data (recast as 32-bit words) and number of 32-bit words
    uint32_t* data_u32 = 0;
    uint32_t  size_u32 = 0;
  
    // Search mode
    if ( triggerFedId_ < 0 ) { 
    
      uint16_t ifed = 0;
      while ( triggerFedId_ < 0 && 
          ifed < 1 + FEDNumbering::lastFEDId() ) {
    const FEDRawData& trigger_fed = buffers.FEDData( ifed );
    if ( trigger_fed.data() && trigger_fed.size() ) {
      uint8_t*  temp = const_cast<uint8_t*>( trigger_fed.data() );
      data_u32 = reinterpret_cast<uint32_t*>( temp ) + sizeof(fedh_t)/sizeof(uint32_t) + 1;
      size_u32 = trigger_fed.size()/sizeof(uint32_t) - sizeof(fedh_t)/sizeof(uint32_t) - 1;
      fedt_t* fed_trailer = reinterpret_cast<fedt_t*>( temp + trigger_fed.size() - sizeof(fedt_t) );
      if ( fed_trailer->conscheck == 0xDEADFACE ) { 
        triggerFedId_ = ifed; 
        if ( edm::isDebugEnabled() ) {
          std::stringstream ss;
          ss << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
         << " Search mode for 'trigger FED' activated!"
         << " Found 'trigger FED' info with id " << triggerFedId_;
          LogTrace("SiStripRawToDigi") << ss.str();
        }
      }
    }
    ifed++;
      }
      if ( triggerFedId_ < 0 ) {
    triggerFedId_ = 0;
    if ( edm::isDebugEnabled() ) {
      std::stringstream ss;
      ss << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
         << " Search mode for 'trigger FED' activated!"
         << " 'Trigger FED' info not found!";
      edm::LogWarning(sistrip::mlRawToDigi_) << ss.str();
    }
      }  
    } 

    // "Trigger FED" id given in .cfg file
    else if ( triggerFedId_ > 0 ) { 
    
      const FEDRawData& trigger_fed = buffers.FEDData( triggerFedId_ );
      if ( trigger_fed.data() && trigger_fed.size() ) {
    uint8_t*  temp = const_cast<uint8_t*>( trigger_fed.data() );
    data_u32 = reinterpret_cast<uint32_t*>( temp ) + sizeof(fedh_t)/sizeof(uint32_t) + 1;
    size_u32 = trigger_fed.size()/sizeof(uint32_t) - sizeof(fedh_t)/sizeof(uint32_t) - 1;
    fedt_t* fed_trailer = reinterpret_cast<fedt_t*>( temp + trigger_fed.size() - sizeof(fedt_t) );
    if ( fed_trailer->conscheck != 0xDEADFACE ) { 
      if ( edm::isDebugEnabled() ) {
        edm::LogWarning(sistrip::mlRawToDigi_) 
          << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
          << " Unexpected stamp found in DAQ trailer (ie, not 0xDEADFACE)!"
          << " Buffer appears not to contain 'trigger FED' data!";
      }
      triggerFedId_ = 0; 
    }
      }
      
    } else { 
      triggerFedId_ = 0; 
      data_u32 = 0;
      size_u32 = 0;
    }
  
    // Populate summary object with commissioning information
    if ( triggerFedId_ > 0 ) { 

      // Some checks
      if ( !data_u32 ) {
    if ( edm::isDebugEnabled() ) {
      std::stringstream ss;
      ss << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
         << " NULL pointer to 'trigger FED' data";
      edm::LogWarning(sistrip::mlRawToDigi_) << ss.str();
    }
    return;
      } 
      if ( size_u32 < sizeof(TFHeaderDescription)/sizeof(uint32_t) ) {
    if ( edm::isDebugEnabled() ) {
      std::stringstream ss;
      ss << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
         << " Unexpected 'Trigger FED' data size [32-bit words]: " << size_u32;
      edm::LogWarning(sistrip::mlRawToDigi_) << ss.str();
    }
    return;
      }
    
      // Write event-specific data to event
      TFHeaderDescription* header = (TFHeaderDescription*) data_u32;
      summary.event( static_cast<uint32_t>( header->getFedEventNumber()) );
      summary.bx( static_cast<uint32_t>( header->getBunchCrossing()) );
    
      // Write commissioning information to event 
      uint32_t hsize = sizeof(TFHeaderDescription)/sizeof(uint32_t);
      uint32_t* head = &data_u32[hsize];
      summary.commissioningInfo( head, event );
      summary.triggerFed( triggerFedId_ );
    
    }

    // Some debug
    if ( summary.isSet() && once_ ) {
      if ( edm::isDebugEnabled() ) {
    std::stringstream ss;
    ss << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
       << " EventSummary built from \"trigger FED\":" 
       << std::endl << summary;
    LogTrace("SiStripRawToDigi") << ss.str();
      }
      once_ = false;
    }
  }

void sistrip::RawToDigiUnpacker::update ( RawDigis &  scope_mode, RawDigis &  virgin_raw, RawDigis &  proc_raw, Digis &  zero_suppr, RawDigis &  common_mode  )

private

fill DetSetVectors using registries

Definition at line 480 of file SiStripRawToDigiUnpacker.cc.

References cm_work_digis_, cm_work_registry_, filterCSVwithJSON::copy, cond::rpcobgas::detid, end, edm::hlt::Exception, extractCm_, min(), proc_work_digis_, proc_work_registry_, scope_work_digis_, scope_work_registry_, python.multivaluedict::sort(), edm::DetSetVector< T >::swap(), virgin_work_digis_, virgin_work_registry_, zs_work_digis_, and zs_work_registry_.

Referenced by progressbar.ProgressBar::__next__(), relval_steps.Matrix::__setitem__(), relval_steps.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), createDigis(), 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(), relval_steps.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 ( ! zs_work_registry_.empty() ) {
      std::sort( zs_work_registry_.begin(), zs_work_registry_.end() );
      std::vector< edm::DetSet<SiStripDigi> > sorted_and_merged;
      sorted_and_merged.reserve(  std::min(zs_work_registry_.size(), size_t(17000)) );
    
      bool errorInData = false;
      std::vector<Registry>::iterator it = zs_work_registry_.begin(), it2 = it+1, end = zs_work_registry_.end();
      while (it < end) {
    sorted_and_merged.push_back( edm::DetSet<SiStripDigi>(it->detid) );
    std::vector<SiStripDigi> & digis = sorted_and_merged.back().data;
    // first count how many digis we have
    size_t len = it->length;
    for (it2 = it+1; (it2 != end) && (it2->detid == it->detid); ++it2) { len += it2->length; }
    // reserve memory 
    digis.reserve(len);
    // push them in
    for (it2 = it+0; (it2 != end) && (it2->detid == it->detid); ++it2) {
      digis.insert( digis.end(), & zs_work_digis_[it2->index], & zs_work_digis_[it2->index + it2->length] );
    }
    it = it2;
      }
    
      // check sorting
      if (!__gnu_cxx::is_sorted( sorted_and_merged.begin(), sorted_and_merged.end() )) {
    // this is an error in the code: i DID sort it already!
    throw cms::Exception("Bug Found") 
      << "Container must be already sorted!\nat " 
      << __FILE__ 
      << ", line " 
      << __LINE__ 
      <<"\n";
      }
    
      std::vector< edm::DetSet<SiStripDigi> >::iterator iii = sorted_and_merged.begin();
      std::vector< edm::DetSet<SiStripDigi> >::iterator jjj = sorted_and_merged.end();
      for ( ; iii != jjj; ++iii ) { 
    if ( ! __gnu_cxx::is_sorted( iii->begin(), iii->end() ) ) {
      // this might be an error in the data, if the raws from one FED are not sorted
      iii->clear(); 
      errorInData = true;
    }
      } 
    
      // output error
      if (errorInData) edm::LogWarning("CorruptData") << "Some modules contained corrupted ZS raw data, and have been skipped in unpacking\n";
    
      // make output DetSetVector
      edm::DetSetVector<SiStripDigi> zero_suppr_dsv( sorted_and_merged, true ); 
      zero_suppr.swap( zero_suppr_dsv );
    } 
  
    // Populate final DetSetVector container with VR data 
    if ( !virgin_work_registry_.empty() ) {

      std::sort( virgin_work_registry_.begin(), virgin_work_registry_.end() );
    
      std::vector< edm::DetSet<SiStripRawDigi> > sorted_and_merged;
      sorted_and_merged.reserve( std::min(virgin_work_registry_.size(), size_t(17000)) );
    
      bool errorInData = false;
      std::vector<Registry>::iterator it = virgin_work_registry_.begin(), it2, end = virgin_work_registry_.end();
      while (it < end) {
    sorted_and_merged.push_back( edm::DetSet<SiStripRawDigi>(it->detid) );
    std::vector<SiStripRawDigi> & digis = sorted_and_merged.back().data;
      
    bool isDetOk = true; 
    // first count how many digis we have
    int maxFirstStrip = it->first;
    for (it2 = it+1; (it2 != end) && (it2->detid == it->detid); ++it2) { 
      // duplicated APV or data corruption. DO NOT 'break' here!
      if (it2->first <= maxFirstStrip) { isDetOk = false; continue; } 
      maxFirstStrip = it2->first;                           
    }
    if (!isDetOk) { errorInData = true; it = it2; continue; } // skip whole det
      
    // make room for 256 * (max_apv_pair + 1) Raw Digis
    digis.resize(maxFirstStrip + 256);
    // push them in
    for (it2 = it+0; (it2 != end) && (it2->detid == it->detid); ++it2) {
      // data corruption. DO NOT 'break' here
      if (it->length != 256)  { isDetOk = false; continue; } 
      std::copy( & virgin_work_digis_[it2->index], & virgin_work_digis_[it2->index + it2->length], & digis[it2->first] );
    }
    if (!isDetOk) { errorInData = true; digis.clear(); it = it2; continue; } // skip whole det
    it = it2;
      }
    
      // output error
      if (errorInData) edm::LogWarning("CorruptData") << "Some modules contained corrupted virgin raw data, and have been skipped in unpacking\n"; 
    
      // check sorting
      if ( !__gnu_cxx::is_sorted( sorted_and_merged.begin(), sorted_and_merged.end()  ) ) {
    // this is an error in the code: i DID sort it already!
    throw cms::Exception("Bug Found") 
      << "Container must be already sorted!\nat " 
      << __FILE__ 
      << ", line " 
      << __LINE__ 
      <<"\n";
      }
    
      // make output DetSetVector
      edm::DetSetVector<SiStripRawDigi> virgin_raw_dsv( sorted_and_merged, true ); 
      virgin_raw.swap( virgin_raw_dsv );
    }
  
    // Populate final DetSetVector container with VR data 
    if ( !proc_work_registry_.empty() ) {
      std::sort( proc_work_registry_.begin(), proc_work_registry_.end() );
    
      std::vector< edm::DetSet<SiStripRawDigi> > sorted_and_merged;
      sorted_and_merged.reserve( std::min(proc_work_registry_.size(), size_t(17000)) );
    
      bool errorInData = false;
      std::vector<Registry>::iterator it = proc_work_registry_.begin(), it2, end = proc_work_registry_.end();
      while (it < end) {
    sorted_and_merged.push_back( edm::DetSet<SiStripRawDigi>(it->detid) );
    std::vector<SiStripRawDigi> & digis = sorted_and_merged.back().data;
      
    bool isDetOk = true; 
    // first count how many digis we have
    int maxFirstStrip = it->first;
    for (it2 = it+1; (it2 != end) && (it2->detid == it->detid); ++it2) { 
      // duplicated APV or data corruption. DO NOT 'break' here!
      if (it2->first <= maxFirstStrip) { isDetOk = false; continue; } 
      maxFirstStrip = it2->first;                           
    }
    // skip whole det
    if (!isDetOk) { errorInData = true; it = it2; continue; } 
      
    // make room for 256 * (max_apv_pair + 1) Raw Digis
    digis.resize(maxFirstStrip + 256);
    // push them in
    for (it2 = it+0; (it2 != end) && (it2->detid == it->detid); ++it2) {
      // data corruption. DO NOT 'break' here
      if (it->length != 256)  { isDetOk = false; continue; } 
      std::copy( & proc_work_digis_[it2->index], & proc_work_digis_[it2->index + it2->length], & digis[it2->first] );
    }
    // skip whole det
    if (!isDetOk) { errorInData = true; digis.clear(); it = it2; continue; } 
    it = it2;
      }
    
      // output error
      if (errorInData) edm::LogWarning("CorruptData") << "Some modules contained corrupted proc raw data, and have been skipped in unpacking\n";
    
      // check sorting
      if ( !__gnu_cxx::is_sorted( sorted_and_merged.begin(), sorted_and_merged.end()  ) ) {
    // this is an error in the code: i DID sort it already!
    throw cms::Exception("Bug Found") 
      << "Container must be already sorted!\nat " 
      << __FILE__ 
      << ", line " 
      << __LINE__ 
      <<"\n";
      }
    
      // make output DetSetVector
      edm::DetSetVector<SiStripRawDigi> proc_raw_dsv( sorted_and_merged, true ); 
      proc_raw.swap( proc_raw_dsv );
    }
  
    // Populate final DetSetVector container with SM data 
    if ( !scope_work_registry_.empty() ) {
      std::sort( scope_work_registry_.begin(), scope_work_registry_.end() );
    
      std::vector< edm::DetSet<SiStripRawDigi> > sorted_and_merged;
      sorted_and_merged.reserve( scope_work_registry_.size() );
    
      bool errorInData = false;
      std::vector<Registry>::iterator it, end;
      for (it = scope_work_registry_.begin(), end = scope_work_registry_.end() ; it != end; ++it) {
    sorted_and_merged.push_back( edm::DetSet<SiStripRawDigi>(it->detid) );
    std::vector<SiStripRawDigi> & digis = sorted_and_merged.back().data;
    digis.insert( digis.end(), & scope_work_digis_[it->index], & scope_work_digis_[it->index + it->length] );
      
    if ( (it +1 != end) && (it->detid == (it+1)->detid) ) {
      errorInData = true; 
      // let's skip *all* the detsets for that key, as we don't know which is the correct one!
      do { ++it; } while ( ( it+1 != end) && (it->detid == (it+1)->detid) );
    }
      }

      // output error
      if (errorInData) edm::LogWarning("CorruptData") << "Some fed keys contained corrupted scope mode data, and have been skipped in unpacking\n"; 

      // check sorting
      if ( !__gnu_cxx::is_sorted( sorted_and_merged.begin(), sorted_and_merged.end()  ) ) {
    // this is an error in the code: i DID sort it already!
    throw cms::Exception("Bug Found") 
      << "Container must be already sorted!\nat " 
      << __FILE__ 
      << ", line " 
      << __LINE__ 
      <<"\n";
      }

      // make output DetSetVector
      edm::DetSetVector<SiStripRawDigi> scope_mode_dsv( sorted_and_merged, true ); 
      scope_mode.swap( scope_mode_dsv );
    }

    // Populate DetSetVector with Common Mode values 
    if ( extractCm_ ) {

      // Populate final DetSetVector container with VR data 
      if ( !cm_work_registry_.empty() ) {

    std::sort( cm_work_registry_.begin(), cm_work_registry_.end() );
    
    std::vector< edm::DetSet<SiStripRawDigi> > sorted_and_merged;
    sorted_and_merged.reserve( std::min(cm_work_registry_.size(), size_t(17000)) );
    
    bool errorInData = false;
    std::vector<Registry>::iterator it = cm_work_registry_.begin(), it2, end = cm_work_registry_.end();
    while (it < end) {
      sorted_and_merged.push_back( edm::DetSet<SiStripRawDigi>(it->detid) );
      std::vector<SiStripRawDigi> & digis = sorted_and_merged.back().data;
      
      bool isDetOk = true; 
      // first count how many digis we have
      int maxFirstStrip = it->first;
      for (it2 = it+1; (it2 != end) && (it2->detid == it->detid); ++it2) { 
        // duplicated APV or data corruption. DO NOT 'break' here!
        if (it2->first <= maxFirstStrip) { isDetOk = false; continue; } 
        maxFirstStrip = it2->first;                           
      }
      if (!isDetOk) { errorInData = true; it = it2; continue; } // skip whole det
      
      // make room for 2 * (max_apv_pair + 1) Common mode values
      digis.resize(maxFirstStrip + 2);
      // push them in
      for (it2 = it+0; (it2 != end) && (it2->detid == it->detid); ++it2) {
        // data corruption. DO NOT 'break' here
        if (it->length != 2)  { isDetOk = false; continue; } 
        std::copy( & cm_work_digis_[it2->index], & cm_work_digis_[it2->index + it2->length], & digis[it2->first] );
      }
      if (!isDetOk) { errorInData = true; digis.clear(); it = it2; continue; } // skip whole det
      it = it2;
    }
    
    // output error
    if (errorInData) edm::LogWarning("CorruptData") << "Some modules contained corrupted common mode data, and have been skipped in unpacking\n"; 
    
    // check sorting
    if ( !__gnu_cxx::is_sorted( sorted_and_merged.begin(), sorted_and_merged.end()  ) ) {
      // this is an error in the code: i DID sort it already!
      throw cms::Exception("Bug Found") 
        << "Container must be already sorted!\nat " 
        << __FILE__ 
        << ", line " 
        << __LINE__ 
        <<"\n";
    }
    
    // make output DetSetVector
    edm::DetSetVector<SiStripRawDigi> common_mode_dsv( sorted_and_merged, true ); 
    common_mode.swap( common_mode_dsv );
      }
      
    }
    
  }

void sistrip::RawToDigiUnpacker::updateEventSummary ( const sistrip::FEDBuffer &  fed, SiStripEventSummary &  summary  )

private

sets the SiStripEventSummary -> not yet implemented for FEDBuffer class

Definition at line 1015 of file SiStripRawToDigiUnpacker.cc.

References SiStripEventSummary::commissioningInfo(), sistrip::FEDBufferBase::daqEventType(), sistrip::FEDFullDebugHeader::daqRegister(), sistrip::FEDFullDebugHeader::daqRegister2(), SiStripEventSummary::fedReadoutMode(), sistrip::FEDBuffer::feHeader(), sistrip::HEADER_TYPE_FULL_DEBUG, sistrip::FEDBufferBase::headerType(), sistrip::invalid32_, edm::isDebugEnabled(), SiStripEventSummary::isSet(), LogTrace, once_, contentValuesCheck::ss, edmLumisInFiles::summary, SiStripEventSummary::triggerFed(), and triggerFedId_.

Referenced by createDigis().

                                                                                                        {
  
    // Retrieve contents of DAQ registers

    sistrip::FEDDAQEventType readout_mode = fed.daqEventType(); 
    uint32_t daq1 = sistrip::invalid32_;
    uint32_t daq2 = sistrip::invalid32_;

    if (fed.headerType() == sistrip::HEADER_TYPE_FULL_DEBUG) {
      const sistrip::FEDFullDebugHeader* header = 0;
      header = dynamic_cast<const sistrip::FEDFullDebugHeader*>(fed.feHeader());
      daq1 = static_cast<uint32_t>( header->daqRegister() ); 
      daq2 = static_cast<uint32_t>( header->daqRegister2() ); 
    }

  
    // If FED DAQ registers contain info, update (and possibly overwrite) EventSummary 
    if ( daq1 != 0 && daq1 != sistrip::invalid32_ ) {
    
      summary.triggerFed( triggerFedId_ );
      summary.fedReadoutMode( readout_mode );
      summary.commissioningInfo( daq1, daq2 );
    
      if ( summary.isSet() && once_ ) {
    if ( edm::isDebugEnabled() ) {
      std::stringstream ss;
      ss << "[sistrip::RawToDigiUnpacker::" << __func__ << "]"
         << " EventSummary built from FED DAQ registers:"
         << std::endl << summary;
      LogTrace("SiStripRawToDigi") << ss.str();
    }
    once_ = false;
      }
    }
  }

void sistrip::RawToDigiUnpacker::useDaqRegister ( bool  use )

inline

EventSummary update request -> not yet implemented for FEDBuffer class.

Definition at line 180 of file SiStripRawToDigiUnpacker.h.

{ useDaqRegister_ =  use; }

Friends And Related Function Documentation

friend class RawToClustersLazyUnpacker

friend

Definition at line 31 of file SiStripRawToDigiUnpacker.h.

Member Data Documentation

std::vector<SiStripRawDigi> sistrip::RawToDigiUnpacker::cm_work_digis_

private

Definition at line 145 of file SiStripRawToDigiUnpacker.h.

Referenced by cleanupWorkVectors(), createDigis(), and update().

std::vector<Registry> sistrip::RawToDigiUnpacker::cm_work_registry_

private

Definition at line 138 of file SiStripRawToDigiUnpacker.h.

Referenced by cleanupWorkVectors(), createDigis(), and update().

bool sistrip::RawToDigiUnpacker::doAPVEmulatorCheck_

private

Definition at line 130 of file SiStripRawToDigiUnpacker.h.

Referenced by createDigis().

bool sistrip::RawToDigiUnpacker::doFullCorruptBufferChecks_

private

Definition at line 129 of file SiStripRawToDigiUnpacker.h.

Referenced by createDigis().

uint32_t sistrip::RawToDigiUnpacker::errorThreshold_

private

Definition at line 131 of file SiStripRawToDigiUnpacker.h.

Referenced by createDigis().

uint32_t sistrip::RawToDigiUnpacker::event_

private

other values

Definition at line 123 of file SiStripRawToDigiUnpacker.h.

Referenced by createDigis().

bool sistrip::RawToDigiUnpacker::extractCm_

private

Definition at line 128 of file SiStripRawToDigiUnpacker.h.

Referenced by createDigis(), and update().

int16_t sistrip::RawToDigiUnpacker::fedBufferDumpFreq_

private

Definition at line 115 of file SiStripRawToDigiUnpacker.h.

Referenced by createDigis().

int16_t sistrip::RawToDigiUnpacker::fedEventDumpFreq_

private

Definition at line 116 of file SiStripRawToDigiUnpacker.h.

Referenced by createDigis().

bool sistrip::RawToDigiUnpacker::first_

private

Definition at line 125 of file SiStripRawToDigiUnpacker.h.

Referenced by createDigis().

int16_t sistrip::RawToDigiUnpacker::headerBytes_

private

configurables

Definition at line 114 of file SiStripRawToDigiUnpacker.h.

Referenced by locateStartOfFedBuffer().

bool sistrip::RawToDigiUnpacker::markMissingFeds_

private

Definition at line 120 of file SiStripRawToDigiUnpacker.h.

bool sistrip::RawToDigiUnpacker::once_

private

Definition at line 124 of file SiStripRawToDigiUnpacker.h.

Referenced by triggerFed(), and updateEventSummary().

std::vector<SiStripRawDigi> sistrip::RawToDigiUnpacker::proc_work_digis_

private

Definition at line 144 of file SiStripRawToDigiUnpacker.h.

Referenced by cleanupWorkVectors(), createDigis(), and update().

std::vector<Registry> sistrip::RawToDigiUnpacker::proc_work_registry_

private

Definition at line 137 of file SiStripRawToDigiUnpacker.h.

Referenced by cleanupWorkVectors(), createDigis(), and update().

bool sistrip::RawToDigiUnpacker::quiet_

private

Definition at line 127 of file SiStripRawToDigiUnpacker.h.

Referenced by createDigis().

std::vector<SiStripRawDigi> sistrip::RawToDigiUnpacker::scope_work_digis_

private

Definition at line 143 of file SiStripRawToDigiUnpacker.h.

Referenced by cleanupWorkVectors(), createDigis(), and update().

std::vector<Registry> sistrip::RawToDigiUnpacker::scope_work_registry_

private

Definition at line 136 of file SiStripRawToDigiUnpacker.h.

Referenced by cleanupWorkVectors(), createDigis(), and update().

int16_t sistrip::RawToDigiUnpacker::triggerFedId_

private

Definition at line 117 of file SiStripRawToDigiUnpacker.h.

Referenced by createDigis(), triggerFed(), and updateEventSummary().

bool sistrip::RawToDigiUnpacker::unpackBadChannels_

private

Definition at line 119 of file SiStripRawToDigiUnpacker.h.

Referenced by createDigis(), and RawToDigiUnpacker().

bool sistrip::RawToDigiUnpacker::useDaqRegister_

private

Definition at line 126 of file SiStripRawToDigiUnpacker.h.

Referenced by createDigis().

bool sistrip::RawToDigiUnpacker::useFedKey_

private

Definition at line 118 of file SiStripRawToDigiUnpacker.h.

Referenced by createDigis().

std::vector<SiStripRawDigi> sistrip::RawToDigiUnpacker::virgin_work_digis_

private

Definition at line 142 of file SiStripRawToDigiUnpacker.h.

Referenced by cleanupWorkVectors(), createDigis(), and update().

std::vector<Registry> sistrip::RawToDigiUnpacker::virgin_work_registry_

private

Definition at line 135 of file SiStripRawToDigiUnpacker.h.

Referenced by cleanupWorkVectors(), createDigis(), and update().

std::vector<SiStripDigi> sistrip::RawToDigiUnpacker::zs_work_digis_

private

digi collections

Definition at line 141 of file SiStripRawToDigiUnpacker.h.

Referenced by cleanupWorkVectors(), createDigis(), and update().

std::vector<Registry> sistrip::RawToDigiUnpacker::zs_work_registry_

private

registries

Definition at line 134 of file SiStripRawToDigiUnpacker.h.

Referenced by cleanupWorkVectors(), createDigis(), and update().