CastorMergerData.cc

Go to the documentation of this file.

/*  
 *  \author A. Campbell - DESY
 */
#ifndef HTBDAQ_DATA_STANDALONE
#include "EventFilter/CastorRawToDigi/interface/CastorMergerData.h"
#else
#include "CastorMergerData.h"
#endif
#include <cstring>
#include <cstdio>
#include <algorithm>

CastorMergerData::CastorMergerData() : m_formatVersion(-2), m_rawLength(0), m_rawConst(nullptr), m_ownData(nullptr) { }
CastorMergerData::CastorMergerData(const unsigned short* data, int length) {
  adoptData(data,length);
  m_ownData=nullptr;
}
CastorMergerData::CastorMergerData(const CastorMergerData& hd) : m_formatVersion(hd.m_formatVersion), m_rawLength(hd.m_rawLength), m_rawConst(hd.m_rawConst), m_ownData(nullptr) { }

CastorMergerData::CastorMergerData(int version_to_create) : m_formatVersion(version_to_create) {
  allocate(version_to_create);
}

void CastorMergerData::allocate(int version_to_create) {
  m_formatVersion=version_to_create;
  // the needed space is for the biggest possible event...
  const int needed=0x200;
  // create a buffer big enough...
  m_ownData=new unsigned short[needed];
  m_rawLength=0;
  m_rawConst=m_ownData;
}

CastorMergerData& CastorMergerData::operator=(const CastorMergerData& hd) {
  if (m_ownData==nullptr) {
    m_formatVersion=hd.m_formatVersion;
    m_rawLength=hd.m_rawLength;
    m_rawConst=hd.m_rawConst;
  }
  return (*this);
}

void CastorMergerData::adoptData(const unsigned short* data, int length) {
  m_rawLength=length;
  m_rawConst=data;
  if (m_rawLength<5) {
    m_formatVersion=-2; // invalid!
  } else {
    m_formatVersion=(m_rawConst[4]>>12)&0xF;
  }
}

// check :: not EE, length is reasonable, length matches wordcount
//          length required for tp+daq is correct

bool CastorMergerData::check() const {
     // length checks
    //  minimum length
    if (m_rawLength<6+12) return false;
    //  matches wordcount
    if (m_rawLength!=m_rawConst[m_rawLength-3]) return false;
 
  return true;
}


void CastorMergerData::unpack(
             unsigned char* tp_lengths, unsigned short* tp_samples) const {

  if (tp_lengths!=nullptr) memset(tp_lengths,0,1);

  int tp_words_total,headerLen,trailerLen;
  determineSectionLengths(tp_words_total,headerLen,trailerLen);

  int wordPtr;
  const unsigned short* tpBase=m_rawConst+headerLen;
  // process the trigger primitive words
  if (tp_lengths!=nullptr) {
    for (wordPtr=0; wordPtr<tp_words_total; wordPtr++) {
       tp_samples[tp_lengths[0]]=tpBase[wordPtr];
       tp_lengths[0]++;
    }
  }
 
}
void CastorMergerData::determineSectionLengths(int& tpWords, int& headerWords, int& trailerWords) const {
 
    tpWords=m_rawConst[5]>>8; // should be 8 but could be up to 12
    headerWords=8;
    trailerWords=0; // minimum, may be more...
}

void CastorMergerData::determineStaticLengths(int& headerWords, int& trailerWords) const {
   headerWords=8;
   trailerWords=0; // minimum, may be more...

}
void CastorMergerData::pack(
               unsigned char* tp_lengths, unsigned short* tp_samples) {
  
  int tp_words_total=0, headerLen, trailerLen;
  determineStaticLengths(headerLen,trailerLen);

  tp_words_total=0;
  int isample;

  // trigger words
  unsigned short* ptr=m_ownData+headerLen;
  if (tp_samples!=nullptr && tp_lengths!=nullptr) {
      for (isample=0; isample<tp_lengths[0] && isample<12; isample++) {
    ptr[tp_words_total]=tp_samples[isample];
    tp_words_total++;
      }
  }

   m_ownData[5]=(tp_words_total<<8)|0x1;
   unsigned short totalLen=headerLen+tp_words_total+trailerLen;

   m_rawLength=totalLen;
   m_ownData[totalLen-2]=totalLen/2; // 32-bit words
   m_ownData[totalLen-3]=totalLen;
   m_ownData[totalLen-4]=tp_words_total;
 
}

void CastorMergerData::packHeaderTrailer(int L1Anumber, int bcn, int submodule, int orbitn, int pipeline, int ndd, int nps, int firmwareRev) {
  m_ownData[0]=L1Anumber&0xFF;
  m_ownData[1]=(L1Anumber&0xFFFF00)>>8;

  m_ownData[2]=0x8000; // Version is valid, no error bits - status bits need definition
  m_ownData[3]=((orbitn&0x1F)<<11)|(submodule&0x7FF);
  m_ownData[4]=((m_formatVersion&0xF)<<12)|(bcn&0xFFF);
  m_ownData[5]|=((nps&0xF)<<4)|0x1;
  m_ownData[6]=((firmwareRev&0x70000)>>3)|(firmwareRev&0x1FFF);
  m_ownData[7]=(pipeline&0xFF) | ((ndd&0x1F)<<8);
  m_ownData[m_rawLength-4]&=0x7FF;
  m_ownData[m_rawLength-4]|=(ndd&0x1F)<<11;
  
  m_ownData[m_rawLength-2]=m_rawLength/2; // 32-bit words
  m_ownData[m_rawLength-1]=(L1Anumber&0xFF)<<8;
}

unsigned int CastorMergerData::getOrbitNumber() const { 
  return (m_rawConst[3]>>11);
}

unsigned int CastorMergerData::getFirmwareRevision() const {
  return (m_rawConst[6]);
}