HcalUHTRData.cc

Go to the documentation of this file.

#include "EventFilter/HcalRawToDigi/interface/HcalUHTRData.h"
#include <cstring>

static const int HEADER_LENGTH_16BIT=2*sizeof(uint64_t)/sizeof(uint16_t);

HcalUHTRData::const_iterator::const_iterator(const uint16_t* ptr, const uint16_t* limit) : 
  m_ptr(ptr), m_limit(limit), m_stepclass(0), m_technicalDataType(0)
{
  if (isHeader()) determineMode();
}

HcalUHTRData::const_iterator& HcalUHTRData::const_iterator::operator++() {
  if (m_ptr==m_limit) return *this;
  if (m_stepclass==0) m_ptr++;
  else if (m_stepclass==1) {
    if (m_microstep==0) { m_ptr++; m_microstep++; }
    else { m_microstep--; }
  } 
  else if (m_stepclass==2) {
    if (isHeader()) { m_ptr++; }
    else { m_ptr+=2; }
  }

  if (isHeader()) {
    determineMode();
    m_header_ptr = m_ptr;
    m_0th_data_ptr = m_header_ptr + 1;
  }
  return *this;
}

void HcalUHTRData::const_iterator::determineMode() {
  if (!isHeader()) return;
  m_flavor=flavor();
  m_stepclass=0;
  if (m_flavor==5) { m_stepclass=1; m_microstep=0; }
  else if (m_flavor == 2) { m_stepclass=2; }
  if (m_flavor==7) { m_technicalDataType = technicalDataType(); }
}

int HcalUHTRData::const_iterator::errFlags() const {
  if ((m_flavor==7 && m_technicalDataType==15) && !isHeader()) return ((*m_ptr)>>11)&0x1;
  else return ((*m_ptr)>>10)&0x3;
}

bool HcalUHTRData::const_iterator::dataValid() const {
  if ((m_flavor==7 && m_technicalDataType==15) && !isHeader()) return ((*m_ptr)>>10)&0x1;
  else return !(errFlags()&0x2);
}

int HcalUHTRData::const_iterator::technicalDataType() const {
  if (m_flavor==7) return ((*m_ptr)>>8)&0xF;
  else return 0;
}

uint8_t HcalUHTRData::const_iterator::adc() const {
  if (m_flavor==5 && m_microstep==0) return ((*m_ptr)>>8)&0x7F;
  else if (m_flavor==7 && m_technicalDataType==15) return (*m_ptr)&0x7F;
  else return (*m_ptr)&0xFF;
}

uint8_t HcalUHTRData::const_iterator::le_tdc() const {
  if (m_flavor==5 || (m_flavor==7 && m_technicalDataType==15)) return 0x80;
  else if (m_flavor == 2) return (m_ptr[1]&0x3F);
  else return (((*m_ptr)&0x3F00)>>8);
}

bool HcalUHTRData::const_iterator::soi() const {
  if (m_flavor==5 || (m_flavor==7 && m_technicalDataType==15)) return false;
  else if (m_flavor == 2) return (m_ptr[0]&0x2000);
  else return (((*m_ptr)&0x4000));
}

uint8_t HcalUHTRData::const_iterator::te_tdc() const {
  if (m_flavor==2) return(m_ptr[1]>>6)&0x1F;
  else return 0x80;
}

uint8_t HcalUHTRData::const_iterator::capid() const {
  if (m_flavor==2) return(m_ptr[1]>>12)&0x3;
  else if (m_flavor==7 && m_technicalDataType==15) {
    return ((*m_ptr)>>8)&0x3;
  }
  else if (m_flavor == 1 || m_flavor == 0) {
    // For flavor 0,1 we only get the first capid in the header, and so we need
    // to count the number of data rows and figure out which cap we want,
    // knowing that they go 0->1->2->3->0
    return 0;
  }
  else { return 0; }
}

bool HcalUHTRData::const_iterator::ok() const {
  if (m_flavor == 2) { return (m_ptr[0]>>12)&0x1; }
  else if (m_flavor == 4) { return (m_ptr[0]>>13)&0x1; }
  else { return false; }
}

HcalUHTRData::const_iterator HcalUHTRData::begin() const {
  return HcalUHTRData::const_iterator(m_raw16+HEADER_LENGTH_16BIT,m_raw16+(m_rawLength64-1)*sizeof(uint64_t)/sizeof(uint16_t));
}

HcalUHTRData::const_iterator HcalUHTRData::end() const {
  return HcalUHTRData::const_iterator(m_raw16+(m_rawLength64-1)*sizeof(uint64_t)/sizeof(uint16_t),m_raw16+(m_rawLength64-1)*sizeof(uint64_t)/sizeof(uint16_t));
}

HcalUHTRData::HcalUHTRData() : m_formatVersion(-2), m_rawLength64(0), m_raw64(nullptr), m_raw16(nullptr), m_ownData(nullptr) { }

HcalUHTRData::HcalUHTRData(const uint64_t* data, int length) : m_rawLength64(length),m_raw64(data),m_raw16((const uint16_t*)(data)),m_ownData(nullptr) {
  m_formatVersion=(m_raw16[6]>>12)&0xF;
}

HcalUHTRData::HcalUHTRData(const HcalUHTRData& hd) : m_formatVersion(hd.m_formatVersion), m_rawLength64(hd.m_rawLength64), m_raw64(hd.m_raw64), m_raw16(hd.m_raw16), m_ownData(nullptr) { }

HcalUHTRData::HcalUHTRData(int version_to_create) : m_formatVersion(version_to_create) {

  // the needed space is for the biggest possible event...
  // fibers*maxsamples/fiber
  const int needed=(0x20+FIBERS_PER_UHTR*CHANNELS_PER_FIBER_MAX*(10+1))*sizeof(uint16_t)/sizeof(uint64_t);

  m_ownData=new uint64_t[needed];
  memset(m_ownData,0,sizeof(uint64_t)*needed);
  m_rawLength64=0;
  m_raw64=m_ownData;
  m_raw16=(const uint16_t*)m_raw64;
}

HcalUHTRData& HcalUHTRData::operator=(const HcalUHTRData& hd) {
  if (m_ownData==nullptr) {
    m_formatVersion=hd.m_formatVersion;
    m_rawLength64=hd.m_rawLength64;
    m_raw64=hd.m_raw64;
    m_raw16=hd.m_raw16;
  }
  return (*this);
}