MatacqRawEvent.cc

Go to the documentation of this file.

// -*- Mode: C++; c-basic-offset: 2; indent-tabs-mode: t; tab-width: 8; -*-
/*
 * Original author: Ph. Gras CEA/Saclay 
 */

#include <unistd.h>
#include <cstdlib>
#include <fstream>
#include <vector>
#include <iostream>
#include <iomanip>
#include <ctime>
#include <limits>

#define CMSSW

#ifdef CMSSW //compilation within CMSSW framework

#include "EventFilter/EcalRawToDigi/interface/MatacqRawEvent.h"
#include "FWCore/Utilities/interface/Exception.h"

static inline void throwExcept(const std::string& s){
  throw cms::Exception("Matacq") << s;
}

#else //compilation outside CMSSW framework (e.g. online)

#include "MatacqRawEvent.h"
#include <stdexcept>
static inline void throwExcept(const std::string& s){
  throw std::runtime_error(s.c_str());
}

#endif //CMSSW not defined

using namespace std;

//DAQ header fields:
const MatacqRawEvent::field32spec_t MatacqRawEvent::fov32               = {0, 0x000000F0};
const MatacqRawEvent::field32spec_t MatacqRawEvent::fedId32             = {0, 0x000FFF00};
const MatacqRawEvent::field32spec_t MatacqRawEvent::bxId32              = {0, 0xFFF00000};
const MatacqRawEvent::field32spec_t MatacqRawEvent::lv132               = {1, 0x00FFFFFF};
const MatacqRawEvent::field32spec_t MatacqRawEvent::triggerType32       = {1, 0x0F000000};
const MatacqRawEvent::field32spec_t MatacqRawEvent::boeType32           = {1, 0xF0000000};
const MatacqRawEvent::field32spec_t MatacqRawEvent::dccLen32            = {2, 0x00FFFFFF};
const MatacqRawEvent::field32spec_t MatacqRawEvent::dccErrors32         = {2, 0xFF000000};
const MatacqRawEvent::field32spec_t MatacqRawEvent::runNum32            = {3, 0x00FFFFFF};
const MatacqRawEvent::field32spec_t MatacqRawEvent::h1Marker32          = {3, 0x3F000000};

//Matacq header fields:
const MatacqRawEvent::field32spec_t MatacqRawEvent::formatVersion32     = {4, 0x0000FFFF};
const MatacqRawEvent::field32spec_t MatacqRawEvent::freqGHz32           = {4, 0x00FF0000};
const MatacqRawEvent::field32spec_t MatacqRawEvent::channelCount32      = {4, 0xFF000000};
const MatacqRawEvent::field32spec_t MatacqRawEvent::timeStamp32         = {5, 0xFFFFFFFF};
//  for data format version >=2:
const MatacqRawEvent::field32spec_t MatacqRawEvent::tTrigPs32           = {6, 0xFFFFFFFF};
//  for data format version >=3:
const MatacqRawEvent::field32spec_t MatacqRawEvent::orbitId32           = {7, 0xFFFFFFFF};
const MatacqRawEvent::field32spec_t MatacqRawEvent::vernier0_32         = {8, 0x0000FFFF};
const MatacqRawEvent::field32spec_t MatacqRawEvent::vernier1_32         = {8, 0xFFFF0000};
const MatacqRawEvent::field32spec_t MatacqRawEvent::vernier2_32         = {9, 0x0000FFFF};
const MatacqRawEvent::field32spec_t MatacqRawEvent::vernier3_32         = {9, 0xFFFF0000};
const MatacqRawEvent::field32spec_t MatacqRawEvent::timeStampMicroSec32 = {10,0xFFFFFFFF};
const MatacqRawEvent::field32spec_t MatacqRawEvent::trigRec32           = {11,0xFF000000};
const MatacqRawEvent::field32spec_t MatacqRawEvent::postTrig32          = {11,0x0000FFFF};
const MatacqRawEvent::field32spec_t MatacqRawEvent::laserPower32        = {12,0x000000FF};
const MatacqRawEvent::field32spec_t MatacqRawEvent::attenuation_dB32    = {12,0x00000F00};
const MatacqRawEvent::field32spec_t MatacqRawEvent::emtcPhase32         = {12,0x0000F000};
const MatacqRawEvent::field32spec_t MatacqRawEvent::emtcDelay32         = {12,0xFFFF0000};
const MatacqRawEvent::field32spec_t MatacqRawEvent::delayA32            = {13,0x0000FFFF};
const MatacqRawEvent::field32spec_t MatacqRawEvent::dccId32             = {13,0x003F0000}; 
const MatacqRawEvent::field32spec_t MatacqRawEvent::color32             = {13,0x00600000};
const MatacqRawEvent::field32spec_t MatacqRawEvent::trigType32          = {13,0x07000000};
const MatacqRawEvent::field32spec_t MatacqRawEvent::side32              = {13,0x08000000};

void MatacqRawEvent::setRawData(const unsigned char* pData, size_t maxSize){
  error = 0;
  const int16le_t *begin16 = (const int16le_t *) pData;
  const uint32le_t *begin32 = (const uint32le_t *) pData;
  if(maxSize < 6*4){
    error = errorLength;
    return;
  }
  
  daqHeader = begin32;
  matacqDataFormatVersion = read32(begin32, formatVersion32);
  freqGHz             = read32(begin32, freqGHz32);
  channelCount        = read32(begin32, channelCount32);
  timeStamp.tv_sec        = read32(begin32, timeStamp32);
  int headerLen = 24; //in bytes
  if(matacqDataFormatVersion>=2){
    tTrigPs       = read32(begin32, tTrigPs32);
    headerLen += 4;
  } else{
    tTrigPs = numeric_limits<int>::max();
  }

  if(matacqDataFormatVersion>=3){
    orbitId           = read32(begin32, orbitId32);
    vernier[0]        = read32(begin32, vernier0_32);
    vernier[1]        = read32(begin32, vernier1_32);
    vernier[2]        = read32(begin32, vernier2_32);
    vernier[3]        = read32(begin32, vernier3_32);
    timeStamp.tv_usec = read32(begin32, timeStampMicroSec32);
    trigRec           = read32(begin32, trigRec32, true);
    postTrig          = read32(begin32, postTrig32);
    delayA            = read32(begin32, delayA32, true);
    emtcDelay         = read32(begin32, emtcDelay32, true);
    emtcPhase         = read32(begin32, emtcPhase32, true);
    attenuation_dB    = read32(begin32, attenuation_dB32, true);
    laserPower        = read32(begin32, laserPower32, true);
    headerLen = 64;
  } else{
    orbitId = 0;
    vernier[0] = -1;
    vernier[1] = -1;
    vernier[2] = -1;
    vernier[3] = -1;
    trigRec = -1;
    postTrig = -1;
    delayA = -1;
    emtcDelay = -1;
    emtcPhase = -1;
    attenuation_dB = -1;
    laserPower = -1;
  }
    
  const int nCh = getChannelCount();
  channelData.resize(nCh);

  const int16le_t *pData16 = (const int16le_t *) (pData+headerLen); 

  for(int iCh=0; iCh<nCh; ++iCh){
    if((size_t)(pData16-begin16)>maxSize){
      throwExcept(string("Corrupted or truncated data"));
    }
    //channel id:
    channelData[iCh].chId = *(pData16++);
    //number of time samples for this channel:
    channelData[iCh].nSamples = *(pData16++);
    //pointer to time sample data of this channel:
    channelData[iCh].samples = pData16;
    //moves to next channel data block:
    if(channelData[iCh].nSamples<0){
      throwExcept(string("Corrupted or truncated data"));
    }
    pData16 += channelData[iCh].nSamples;
  }
  
  //data trailer chekes:
  //FED header is aligned on 64-bit=>padding to skip
  int padding = (4-(pData16-begin16))%4;
  if(padding<0) padding+=4;
  pData16 += padding;
  if((size_t)(pData16-begin16)>maxSize){
    throwExcept(string("Corrupted or truncated data"));
  }
  const uint32le_t* trailer32 = (const uint32le_t*)(pData16);
  fragLen = trailer32[1]&0xFFFFFF;
  
  //cout << "Event fragment length including headers: " << fragLen
  //     << " 64-bit words\n";
  
  //FIXME: I am expecting the event length specifies in the header to
  //include the header, while it is not the case in current TB 2006 data
  const int nHeaders = 3;
  if(fragLen!=read32(begin32,dccLen32)+nHeaders
     && fragLen != read32(begin32,dccLen32)){
    //cout << "Error: fragment length is not consistent with DCC "
    //  "length\n";
    error |= errorLengthConsistency;
  }
  
  //skip trailers
  const int trailerLen = 4;
  pData16 += trailerLen;
  
  parsedLen = (pData16-begin16) / 4;

  if((pData16-begin16)!=(4*fragLen)){
    error |= errorLength;
  }

  if((size_t)(pData16-begin16)>maxSize){
    throwExcept(string("Corrupted or truncated data"));
  }

  //some checks
  if(getBoe()!=0x5){
    error |= errorWrongBoe;
  }
}

int MatacqRawEvent::read32(const uint32le_t* pData, field32spec_t spec32,
               bool ovfTrans){
  uint32_t result =  pData[spec32.offset] & spec32.mask;
  uint32_t mask = spec32.mask;
  while((mask&0x1) == 0){
      mask >>= 1;
      result >>= 1;
  }
  if(ovfTrans){
    //overflow bit (MSB) mask:
    mask = ((mask >>1) + 1);
    if(result & mask)  result = (uint32_t)-1;
  }
  return result;
}