DataFormats/CSCDigi/src/CSCComparatorDigi.cc - master - githubmirror/cmssw

DataFormats/CSCDigi/src/CSCComparatorDigi.cc
122 строки4 KB
Matti Kortelainen
Move stored types in io_v1 namespace, reset class versions to 3
22 май 2026, 18:18
22 май 2026, 18:18909f324
/** \file
 *
 *
 * \author M.Schmitt, Northwestern
 */
#include "DataFormats/CSCDigi/interface/CSCComparatorDigi.h"
#include "DataFormats/CSCDigi/interface/CSCConstants.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <iostream>
#include <algorithm>
#include <iterator>

using namespace std;

namespace io_v1 {

  // Constructors
  CSCComparatorDigi::CSCComparatorDigi(int strip, int comparator, int timeBinWord)
      : strip_(strip), comparator_(comparator), timeBinWord_(timeBinWord) {}

  CSCComparatorDigi::CSCComparatorDigi() : strip_(0), comparator_(0), timeBinWord_(0) {}

  // Comparison

  bool CSCComparatorDigi::operator==(const CSCComparatorDigi& digi) const {
    if (getStrip() != digi.getStrip())
      return false;
    if (getComparator() != digi.getComparator())
      return false;
    if (getTimeBinWord() != digi.getTimeBinWord())
      return false;
    return true;
  }

  //@@ If one wanted to order comparator digis how would one want op< to behave?
  // I don't know...
  // I think LHS < RHS only makes sense if
  // i) time(LHS) .eq. time(RHS)
  // AND
  // ii) strip(LHS) .lt. strip(RHS)
  // But I don't see how this can be useful.

  bool CSCComparatorDigi::operator<(const CSCComparatorDigi& digi) const {
    bool result = false;
    if (getTimeBin() == digi.getTimeBin()) {
      result = (getStrip() < digi.getStrip());
    }
    return result;
  }

  // Getters

  int CSCComparatorDigi::getTimeBin() const {
    // Find first bin which fired, counting from 0
    uint16_t tbit = 1;
    int tbin = -1;
    for (int i = 0; i < 16; ++i) {
      if (tbit & timeBinWord_) {
        tbin = i;
        break;
      }
      tbit = tbit << 1;
    }
    return tbin;
  }

  /// Get the distrip number. Counts from 0.
  // originally defined in EventFilter/CSCRawToDigi/src/CSCComparatorData.cc
  int CSCComparatorDigi::getDiStrip() const { return ((strip_ - 1) % CSCConstants::NUM_STRIPS_PER_CFEB) / 2; }

  /// Get the CFEB number. Counts from 0.
  // originally defined in EventFilter/CSCRawToDigi/src/CSCComparatorData.cc
  int CSCComparatorDigi::getCFEB() const { return (strip_ - 1) / CSCConstants::NUM_STRIPS_PER_CFEB; }

  // This definition is consistent with the one used in
  // the function CSCComparatorData::add() in EventFilter/CSCRawToDigi
  // The halfstrip counts from 0!
  int CSCComparatorDigi::getHalfStrip() const { return (getStrip() - 1) * 2 + getComparator(); }

  // Return the fractional half-strip
  float CSCComparatorDigi::getFractionalStrip() const { return getStrip() + getComparator() * 0.5f - 0.75f; }

  std::vector<int> CSCComparatorDigi::getTimeBinsOn() const {
    std::vector<int> tbins;
    uint16_t tbit = timeBinWord_;
    const uint16_t one = 1;
    for (int i = 0; i < 16; ++i) {
      if (tbit & one)
        tbins.push_back(i);
      tbit = tbit >> 1;
      if (tbit == 0)
        break;  // end already if no more bits set
    }
    return tbins;
  }

  // Setters
  //@@ No way to set time word?

  void CSCComparatorDigi::setStrip(int strip) { strip_ = strip; }
  void CSCComparatorDigi::setComparator(int comparator) { comparator_ = comparator; }

  // Output

  void CSCComparatorDigi::print() const {
    std::ostringstream ost;
    ost << "CSCComparatorDigi | strip " << getStrip() << " | comparator " << getComparator() << " | first time bin "
        << getTimeBin() << " | time bins on ";
    std::vector<int> tbins = getTimeBinsOn();
    std::copy(tbins.begin(), tbins.end(), std::ostream_iterator<int>(ost, " "));
    edm::LogVerbatim("CSCDigi") << ost.str();
  }

  std::ostream& operator<<(std::ostream& o, const CSCComparatorDigi& digi) {
    o << "CSCComparatorDigi Strip:" << digi.getStrip() << ", Comparator: " << digi.getComparator()
      << ", First Time Bin On: " << digi.getTimeBin() << ", Time Bins On: ";
    std::vector<int> tbins = digi.getTimeBinsOn();
    std::copy(tbins.begin(), tbins.end(), std::ostream_iterator<int>(o, " "));
    return o;
  }

}  // namespace io_v1