HgcrocDigiCollection.cxx

 
#include "Recon/Event/HgcrocDigiCollection.h"
 
ClassImp(ldmx::HgcrocDigiCollection);
 
namespace ldmx {
HgcrocDigiCollection::Sample::Sample(bool tot_progress, bool tot_complete,
                                     int firstMeas, int seconMeas, int toa,
                                     int version) {
  version_ = version;
  if (version_ == 2) {
    // version 2 HGC ROC had a much simpler word structure
    // 12-bit TOT | 10-bit TOA | 10-bit ADCt
    // we assume that firstMeas is adct and secondMeas is tot
    word_ =
        (((seconMeas > 0xfff ? 0xfff : seconMeas) & 0xfff) << 20) +
        (((firstMeas > TEN_BIT_MASK ? TEN_BIT_MASK : firstMeas) & TEN_BIT_MASK)
         << 10) +
        ((toa > TEN_BIT_MASK ? TEN_BIT_MASK : toa) & TEN_BIT_MASK);
 
  } else {
    if (not tot_progress and tot_complete) {
      // the 12 bit internal tot measurement needs to
      //  be packed into a 10 bit int
      if (seconMeas > 512)
        seconMeas =
            512 + seconMeas / 8;  // lost some precision but can go higher
    }
 
    // check if over largest number possible ==> set to largest if over (don't
    // want wrapping) and then do bit shifting nonsense to code the
    // measurements into the 32-bit word set last measurement to TOA
    word_ =
        (tot_progress << FIRSTFLAG_POS) + (tot_complete << SECONFLAG_POS) +
        (((firstMeas > TEN_BIT_MASK ? TEN_BIT_MASK : firstMeas) & TEN_BIT_MASK)
         << FIRSTMEAS_POS) +
        (((seconMeas > TEN_BIT_MASK ? TEN_BIT_MASK : seconMeas) & TEN_BIT_MASK)
         << SECONMEAS_POS) +
        (((toa > TEN_BIT_MASK ? TEN_BIT_MASK : toa) & TEN_BIT_MASK));
  }
}
 
void HgcrocDigiCollection::clear() {
  channel_ids_.clear();
  samples_.clear();
 
  return;
}
 
std::ostream &operator<<(std::ostream &o, const HgcrocDigiCollection &c) {
  return o << "HgcrocDigiCollection { Num Channel IDs: "
           << c.channel_ids_.size() << ", Num Samples: " << c.samples_.size()
           << ", Samples Per Digi: " << c.num_samples_per_digi_
           << ", Index for SOI: " << c.sample_of_interest_ << "}";
}
 
const HgcrocDigiCollection::HgcrocDigi HgcrocDigiCollection::getDigi(
    unsigned int digiIndex) const {
  return HgcrocDigiCollection::HgcrocDigi(
      channel_ids_.at(digiIndex),
      samples_.begin() + digiIndex * getNumSamplesPerDigi(), *this);
}
 
void HgcrocDigiCollection::addDigi(
    unsigned int id, const std::vector<HgcrocDigiCollection::Sample> &digi) {
  if (digi.size() != this->getNumSamplesPerDigi()) {
    std::cerr << "[ WARN ] [ HgcrocDigiCollection ] Input list of samples "
                 "has size '"
              << digi.size()
              << "' that does not match the number of samples per digi '"
              << this->getNumSamplesPerDigi() << "'!." << std::endl;
    return;
  }
 
  channel_ids_.push_back(id);
  for (auto const &s : digi) samples_.push_back(s.raw());
 
  return;
}
 
void HgcrocDigiCollection::addDigi(unsigned int id,
                                   const std::vector<uint32_t> &digi) {
  if (digi.size() != this->getNumSamplesPerDigi()) {
    std::cerr << "[ WARN ] [ HgcrocDigiCollection ] Input list of samples "
                 "has size '"
              << digi.size()
              << "' that does not match the number of samples per digi '"
              << this->getNumSamplesPerDigi() << "'!." << std::endl;
    return;
  }
 
  channel_ids_.push_back(id);
  for (auto const &s : digi) samples_.push_back(s);
 
  return;
}
}  // namespace ldmx
 
std::ostream &operator<<(std::ostream &s,
                         const ldmx::HgcrocDigiCollection::Sample &sample) {
  s << "Sample { " << "tot prog: " << sample.isTOTinProgress() << ", "
    << "tot comp: " << sample.isTOTComplete() << ", ";
  if (sample.isTOTComplete() and sample.isTOTinProgress())
    s << "adc t: " << sample.adcT() << ", " << "tot: " << sample.tot() << ", ";
  else if (sample.isTOTComplete())
    s << "adc t-1: " << sample.adcTm1() << ", " << "tot: " << sample.tot()
      << ", ";
  else
    s << "adc t-1: " << sample.adcTm1() << ", " << "adc t: " << sample.adcT()
      << ", ";
 
  s << "toa: " << sample.toa() << " }";
  return s;
}
 
std::ostream &operator<<(std::ostream &s,
                         const ldmx::HgcrocDigiCollection::HgcrocDigi &digi) {
  s << "HgcrocDigi { ";
 
  s << " Id: 0x" << std::hex << digi.id() << std::dec << " ";
 
  if (digi.isADC())
    s << "ADC Mode -> SOI: " << digi.soi() << " }";
  else
    s << "TOT Mode -> " << digi.tot() << " }";
 
  return s;
}
 
std::ostream &operator<<(std::ostream &s,
                         const ldmx::HgcrocDigiCollection &col) {
  s << "HgcrocDigiCollection { " << std::endl;
  for (unsigned int i_digi = 0; i_digi < col.getNumDigis(); i_digi++)
    s << "  " << col.getDigi(i_digi) << std::endl;
  s << "}";
  return s;
}