CMS 3D CMS Logo
Classes | Public Types | Public Member Functions | Private Member Functions | Private Attributes | Friends

sistrip::RawToDigiUnpacker Class Reference

#include <SiStripRawToDigiUnpacker.h>

List of all members.

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
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.
void quiet (bool)
 verbosity
 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
void triggerFed (const FEDRawDataCollection &, SiStripEventSummary &, const uint32_t &event)
 trigger info
void useDaqRegister (bool)
 EventSummary update request -> not yet implemented for FEDBuffer class.
 ~RawToDigiUnpacker ()
 default constructor

Private Member Functions

void cleanupWorkVectors ()
 method to clear registries and digi collections
void dumpRawData (uint16_t fed_id, const FEDRawData &, std::stringstream &)
 dumps raw data to stdout (NB: payload is byte-swapped,headers/trailer are not).
sistrip::FedBufferFormat fedBufferFormat (const uint16_t &register_value)
 returns buffer format
sistrip::FedReadoutMode fedReadoutMode (const uint16_t &register_value)
 returns buffer readout mode
void handleException (std::string method_name, std::string extra_info="")
 catches all possible exceptions and rethrows them as cms::Exceptions
void physicalOrder (uint16_t &readout_order, uint16_t &physical_order)
 order of strips
 RawToDigiUnpacker ()
 private default constructor
void readoutOrder (uint16_t &physical_order, uint16_t &readout_order)
 order of strips
void update (RawDigis &scope_mode, RawDigis &virgin_raw, RawDigis &proc_raw, Digis &zero_suppr, RawDigis &common_mode)
 fill DetSetVectors using registries
void updateEventSummary (const sistrip::FEDBuffer &, SiStripEventSummary &)
 sets the SiStripEventSummary -> not yet implemented for FEDBuffer class

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
bool extractCm_
int16_t fedBufferDumpFreq_
int16_t fedEventDumpFreq_
bool first_
int16_t headerBytes_
 configurables
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
std::vector< Registry > zs_work_registry_
 registries

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_, SiStripFedCabling::connections(), FEDRawData::data(), doAPVEmulatorCheck_, doFullCorruptBufferChecks_, dumpRawData(), alignCSCRings::e, errorThreshold_, event_, Exception, extractCm_, sistrip::FED_CABLING, fedBufferDumpFreq_, sistrip::FEDCH_PER_FEUNIT, FEDRawDataCollection::FEDData(), fedEventDumpFreq_, SiStripFedCabling::feds(), sistrip::FINE_DELAY, first_, sistrip::FEDRawChannelUnpacker::hasData(), sistrip::FEDZSChannelUnpacker::hasData(), i, LaserDQM_cfg::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(), FEDRawData::size(), edm::EDCollection< T >::size(), 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.feds().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.feds().begin();
    for ( ; ifed != cabling.feds().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
      const std::vector<FedChannelConnection>& conns = cabling.connections(*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 
) [private]

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(), and sistrip::RawToClustersLazyUnpacker::fill().

                                                                                                      {

    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) [inline, private]

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) [inline, private]

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, exception, edm::isDebugEnabled(), and sistrip::mlRawToDigi_.

                                                                                       { 

    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(), LogTrace, sistrip::mlRawToDigi_, evf::evtn::offset(), FEDRawData::resize(), and FEDRawData::size().

                                                                                                                      {
  
    // 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 
) [inline, private]

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 
) [inline, private]

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(), summarizeEdmComparisonLogfiles::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, 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 createDigis().

                                                                                                                                           {
  
    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_, summarizeEdmComparisonLogfiles::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().

Generated for CMSSW Reference Manual by doxygen 1.7.3