CalorimeterFragmentReader.hh

#ifndef mu2e_artdaq_core_Overlays_CalorimeterFragmentReader_hh
#define mu2e_artdaq_core_Overlays_CalorimeterFragmentReader_hh

#include "artdaq-core/Data/Fragment.hh"
#include "cetlib/exception.h"

#include "mu2e-artdaq-core/Overlays/DetectorFragment.hh"

#include <iostream>
#include <ostream>
#include <vector>
#include <bitset>

namespace mu2e {
  class CalorimeterFragmentReader;
}

class mu2e::CalorimeterFragmentReader: public mu2e::DetectorFragment {

  public:

  CalorimeterFragmentReader(artdaq::Fragment const & f); 
  ~CalorimeterFragmentReader(); 

  // Calorimeter-specific methods:
  mu2e::DetectorFragment::adc_t crystalID();
  mu2e::DetectorFragment::adc_t apdID();
  mu2e::DetectorFragment::adc_t time();
  mu2e::DetectorFragment::adc_t numSamples();
  std::vector<mu2e::DetectorFragment::adc_t> calorimeterADC();
  void printAll();

};

mu2e::CalorimeterFragmentReader::CalorimeterFragmentReader(artdaq::Fragment const & f ) :
  DetectorFragment(f) {
}

mu2e::CalorimeterFragmentReader::~CalorimeterFragmentReader() {
}

mu2e::DetectorFragment::adc_t mu2e::CalorimeterFragmentReader::crystalID() {
  // The calorimeter data packets begin 8*16=128 bits after the beginning of the fragment data
  // Bits 0 to 11
  return mu2e::DetectorFragment::convertFromBinary(bitArray(dataBlockBegin()+8),127-12,127-0);
}

mu2e::DetectorFragment::adc_t mu2e::CalorimeterFragmentReader::apdID() {
  // Bits 12 to 15
  return mu2e::DetectorFragment::convertFromBinary(bitArray(dataBlockBegin()+8),127-16,127-12);
}

mu2e::DetectorFragment::adc_t mu2e::CalorimeterFragmentReader::time() {
  return mu2e::DetectorFragment::convertFromBinary(bitArray(dataBlockBegin()+8),127-32,127-16);
}

mu2e::DetectorFragment::adc_t mu2e::CalorimeterFragmentReader::numSamples() {
  return mu2e::DetectorFragment::convertFromBinary(bitArray(dataBlockBegin()+8),127-48,127-32);
}

std::vector<mu2e::DetectorFragment::adc_t> mu2e::CalorimeterFragmentReader::calorimeterADC() {
  std::vector<adc_t> theVector;
  // Start after the 3rd 16-bit block of the first calorimeter DTC data packet:
  for(size_t adcIdx=0; adcIdx<numSamples(); adcIdx++) {
    theVector.push_back( *((adc_t const *)( (dataBlockBegin()+8+3) + adcIdx )) );
  }
  return theVector;
}

void mu2e::CalorimeterFragmentReader::printAll() {
  std::cout << "\t\t" << "Binary Representation: ";
  printBitArray(bitArray(dataBlockBegin()+8));
  std::cout << "\t\t" << "Crystal IDX: " << (int)crystalID()  << std::endl;
  std::cout << "\t\t" << "APD ID:      " << (int)apdID()      << std::endl;
  std::cout << "\t\t" << "Time:        " << (int)time()       << std::endl;
  std::cout << "\t\t" << "Num Samples: " << (int)numSamples() << std::endl;
  std::cout << "\t\t" << "ADC:         {[";
  std::vector<mu2e::DetectorFragment::adc_t> ADCarray = calorimeterADC();
  for(size_t i=0; i<ADCarray.size(); i++) {
    std::cout << (int)ADCarray[i];
    if(i<ADCarray.size()-1) {
      std::cout << ",";
    }
  }
  std::cout << "]}" << std::endl;
}

#endif /* mu2e_artdaq_core_Overlays_CalorimeterFragmentReader_hh */