PPSTimingCalibration.cc

Go to the documentation of this file.

/****************************************************************************
 *
 * This is a part of CTPPS offline software.
 * Authors:
 *   Filip Dej
 *   Laurent Forthomme
 *
 ****************************************************************************/

#include "CondFormats/CTPPSReadoutObjects/interface/PPSTimingCalibration.h"

//--------------------------------------------------------------------------

bool PPSTimingCalibration::Key::operator<(const PPSTimingCalibration::Key& rhs) const {
  if (db == rhs.db) {
    if (sampic == rhs.sampic) {
      if (channel == rhs.channel)
        return cell < rhs.cell;
      return channel < rhs.channel;
    }
    return sampic < rhs.sampic;
  }
  return db < rhs.db;
}

std::ostream& operator<<(std::ostream& os, const PPSTimingCalibration::Key& key) {
  return os << key.db << " " << key.sampic << " " << key.channel << " " << key.cell;
}

//--------------------------------------------------------------------------

std::vector<double> PPSTimingCalibration::parameters(int key1, int key2, int key3, int key4) const {
  Key key{key1, key2, key3, key4};
  auto out = parameters_.find(key);
  if (out == parameters_.end())
    return {};
  return out->second;
}

double PPSTimingCalibration::timeOffset(int key1, int key2, int key3, int key4) const {
  Key key{key1, key2, key3, key4};
  auto out = timeInfo_.find(key);
  if (out == timeInfo_.end())
    return 0.;
  return out->second.first;
}

double PPSTimingCalibration::timePrecision(int key1, int key2, int key3, int key4) const {
  Key key{key1, key2, key3, key4};
  auto out = timeInfo_.find(key);
  if (out == timeInfo_.end())
    return 0.;
  return out->second.second;
}

std::ostream& operator<<(std::ostream& os, const PPSTimingCalibration& data) {
  os << "FORMULA: " << data.formula_ << "\nDB SAMPIC CHANNEL CELL PARAMETERS TIME_OFFSET\n";
  for (const auto& kv : data.parameters_) {
    os << kv.first << " [";
    for (size_t i = 0; i < kv.second.size(); ++i)
      os << (i > 0 ? ", " : "") << kv.second.at(i);
    const auto& time = data.timeInfo_.at(kv.first);
    os << "] " << time.first << " " << time.second << "\n";
  }
  return os;
}