DAQEventUnpack.cxx

Go to the documentation of this file.

///////////////////////////////////////////////////////////////////////////////
/// \file     DAQEventUnpack.cxx
/// \brief    Producer module to extract DAQ header variables relevant to
///           data quality
/// \author   Justin Vasel <jvasel@indiana.edu>
///           Mark Messier <messier@indiana.edu>
///////////////////////////////////////////////////////////////////////////////

#include <list>

// NOvA Online
#include "DAQDataFormats/RawEvent.h"
#include "DAQDataFormats/RawTrigger.h"
#include "DAQDataFormats/RawTriggerMask.h"
#include "DAQDataFormats/RawDataBlock.h"
#include "DAQDataFormats/RawMicroBlock.h"
#include "DAQDataFormats/RawMicroSlice.h"
#include "DAQDataFormats/RawNanoSlice.h"
#include "DAQChannelMap/DAQChannelMap.h"
#include "RawData/RawSumDropMB.h"

// NOvA Offline
#include "DataQuality/DAQHeader/DAQEventUnpack.h"
#include "DataQuality/DAQHeader/DAQEventSummary.h"


namespace dq {
  //...........................................................................
  DAQEventUnpack::DAQEventUnpack() { }


  //...........................................................................
  DAQEventUnpack::~DAQEventUnpack() { }


  //...........................................................................
  /// Extract raw information at event level and unpack conatined datablocks
  /// \param raw : Raw event
  /// \param daqsum : DAQEventSummary object
  void DAQEventUnpack::UnpackRawEvent(daqdataformats::RawEvent& raw,
              dq::DAQEventSummary*  daqsum)
  {
    if (raw.getHeader()->isEventIncomplete()) daqsum->fEventIncomplete = true;
    unsigned int nDataBlock = raw.getDataBlockNumber();
    unsigned int eventNumber = raw.getEventNumber();

    daqsum->fEventNumber = eventNumber;

    // Now unpack the associated data blocks
    unsigned int idata;
    for (idata = 0; idata < nDataBlock; ++idata) {
      this->UnpackDataBlock(raw, idata, daqsum);
    }
  }


  //...........................................................................
  /// Extract raw information at datablock level and unpack contained microblocks
  /// \param raw : Raw event
  /// \param idata : Data block number
  /// \param daqsum : DAQEventSummary object
  void DAQEventUnpack::UnpackDataBlock(daqdataformats::RawEvent& raw,
               unsigned int          idata,
               dq::DAQEventSummary*  daqsum)
  {
    raw.setFloatingDataBlock(idata);
    daqdataformats::RawDataBlock& datablock = (*raw.getFloatingDataBlock());
    
    // Bail if this is the DAQ status block
    if (datablock.getHeader()->getMarker() == daqdataformats::datablockheader::SummaryBlock_Marker){
      if (datablock.getHeader()->checkMarker()) {
        return;
      }
    }

    if (datablock.isMissingData()) ++daqsum->fNdataBlockMissingData;

    unsigned int nMicroBlock = datablock.getNumMicroBlocks();
    daqsum->fNmicroBlocks = nMicroBlock;

    if (nMicroBlock == 0) ++daqsum->fNemptyDataBlock;

    // Now unpack the associated micro slices
    unsigned int imicro;
    for (imicro = 0; imicro < nMicroBlock; ++imicro) {
      this->UnpackMicroSlice(datablock, imicro, daqsum);
    }
  }


  //...........................................................................
  /// Extract raw information at microblock level and unpack contained nanoslices
  /// \param datablock : Datablock
  /// \param imicro : Microblock number
  /// \param daqsum : DAQEventSummary object
  void DAQEventUnpack::UnpackMicroSlice(daqdataformats::RawDataBlock& datablock,
                unsigned int          imicro,
                dq::DAQEventSummary*  daqsum)
  {
    datablock.setFloatingMicroBlock(imicro);

    daqdataformats::RawMicroBlock&
      microblock = (*datablock.getFloatingMicroBlock());

    daqdataformats::RawMicroBlockHeader&
      microblock_header = (*microblock.getHeader());

    // Get diblock number from microblock header
    daqsum->fDiblocks.push_back(microblock_header.getDiblockId());

    daqdataformats::RawMicroSlice& microslice = (*microblock.getMicroSlice());
    if (!microslice.getDataPresent()) ++daqsum->fNmicroSliceDataNotPresent;

    unsigned int nNanoSlice = microslice.getNumNanoSlices();
    if (nNanoSlice == 0) ++daqsum->fNemptyMicroSlice;
    daqsum->fNtotalNanoSlices += nNanoSlice;

    // Now unpack all the associated nano slices
    unsigned int inano;
    for (inano = 0; inano < nNanoSlice; ++inano) {
      this->UnpackNanoSlice(microslice, inano, daqsum);
    }
  }


  //...........................................................................
  /// Extract raw information at nanoslice level
  /// \param microslice : Microslice
  /// \param inano : Nanoslice number
  /// \param daqsum : DAQEventSummary object
  void DAQEventUnpack::UnpackNanoSlice(daqdataformats::RawMicroSlice& microslice,
               unsigned int          inano,
               dq::DAQEventSummary*  daqsum)
  {
    microslice.setFloatingNanoSlice(inano);

    daqdataformats::RawNanoSlice&       nanoslice    = (*microslice.getFloatingNanoSlice());
    daqdataformats::RawNanoSliceHeader& nanoslicehdr = (*nanoslice.getHeader());

    if (!nanoslicehdr.getDataPresent())   ++daqsum->fNanoSliceDataNotPresent;
    if (!nanoslicehdr.getLinkStatus())    ++daqsum->fNanoSliceNoLinkStatus;
    if ( nanoslicehdr.getBufferEmpty())   ++daqsum->fNanoSliceBufferEmpty;
    if ( nanoslicehdr.getBufferFull())    ++daqsum->fNanoSliceBufferFull;
    if ( nanoslicehdr.getCommError())     ++daqsum->fNanoSliceCommError;
    if ( nanoslicehdr.getPacketError())   ++daqsum->fNanoSlicePacketError;
    if ( nanoslicehdr.getOverflowError()) ++daqsum->fNanoSliceOverflowError;
    if ( nanoslicehdr.getADCError())      ++daqsum->fNanoSliceADCError;
  }
}

////////////////////////////////////////////////////////////////////////////////