Digit.cxx

// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
// This software is distributed under the terms of the GNU General Public
// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
//
// In applying this license CERN does not waive the privileges and immunities
// granted to it by virtue of its status as an Intergovernmental Organization
// or submit itself to any jurisdiction.

#include "DataFormatsTRD/Digit.h"
#include <iostream>
#include <algorithm>
#include "fairlogger/Logger.h"

namespace o2::trd
{

using namespace constants;

Digit::Digit(int det, int row, int pad, ArrayADC adc, int pretrigphase)
{
  setDetector(det);
  setROB(row, pad);
  setMCM(row, pad);
  setADC(adc);
  setChannel(NADCMCM - 2 - (pad % NCOLMCM));
  setPreTrigPhase(pretrigphase);
}

Digit::Digit(int det, int row, int pad) // add adc data in a seperate step
{
  setDetector(det);
  setROB(row, pad);
  setMCM(row, pad);
  setChannel(NADCMCM - 2 - (pad % NCOLMCM));
}

Digit::Digit(int det, int rob, int mcm, int channel, ArrayADC adc, int pretrigphase)
{
  setDetector(det);
  setROB(rob);
  setMCM(mcm);
  setChannel(channel);
  setADC(adc);
  setPreTrigPhase(pretrigphase);
}

Digit::Digit(int det, int rob, int mcm, int channel, int pretrigphase) // add adc data in a seperate step
{
  setDetector(det);
  setROB(rob);
  setMCM(mcm);
  setChannel(channel);
  setPreTrigPhase(pretrigphase);
}

void Digit::setPreTrigPhase(int phase)
{
  mDetector = ((((phase) & 0x3) << 12) | (mDetector & 0xfff));
}

bool Digit::isSharedDigit() const
{
  if (mChannel == 0 || mChannel == 1 || mChannel == NADCMCM - 1) {
    return 1;
  } else {
    return 0;
  }
}

bool Digit::isNeighbour(const Digit& other) const
{
  return (getDetector() == other.getDetector() && getROB() == other.getROB() && getMCM() == other.getMCM() && std::abs(getChannel() - other.getChannel()) == 1);
}

ADC_t Digit::getADCmax(int& idx) const
{
  auto itMax = std::max_element(mADC.begin(), mADC.end());
  idx = std::distance(mADC.begin(), itMax);
  return *itMax;
}

std::ostream& operator<<(std::ostream& stream, const Digit& d)
{
  stream << "Digit Det: " << HelperMethods::getSector(d.getDetector()) << "_" << HelperMethods::getStack(d.getDetector()) << "_" << HelperMethods::getLayer(d.getDetector()) << " pad row: " << d.getPadRow() << " pad column: " << d.getPadCol() << " Channel: " << d.getChannel() << " ADCs:";
  for (int i = 0; i < constants::TIMEBINS; i++) {
    stream << "[" << d.getADC()[i] << "]";
  }
  return stream;
}

} // namespace o2::trd