CastorElectronicsMap.cc

Go to the documentation of this file.

#include <iostream>
#include <set>

#include "CondFormats/CastorObjects/interface/CastorElectronicsMap.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

CastorElectronicsMap::CastorElectronicsMap()
    : mPItems(CastorElectronicsId::maxLinearIndex + 1),
      mTItems(CastorElectronicsId::maxLinearIndex + 1),
      mPItemsById(nullptr),
      mTItemsByTrigId(nullptr) {}

namespace castor_impl {
  class LessById {
  public:
    bool operator()(const CastorElectronicsMap::PrecisionItem* a, const CastorElectronicsMap::PrecisionItem* b) {
      return a->mId < b->mId;
    }
  };
  class LessByTrigId {
  public:
    bool operator()(const CastorElectronicsMap::TriggerItem* a, const CastorElectronicsMap::TriggerItem* b) {
      return a->mTrigId < b->mTrigId;
    }
  };
}  // namespace castor_impl

CastorElectronicsMap::~CastorElectronicsMap() {
  delete mPItemsById.load();
  delete mTItemsByTrigId.load();
}
// copy-ctor
CastorElectronicsMap::CastorElectronicsMap(const CastorElectronicsMap& src)
    : mPItems(src.mPItems), mTItems(src.mTItems), mPItemsById(nullptr), mTItemsByTrigId(nullptr) {}
// copy assignment operator
CastorElectronicsMap& CastorElectronicsMap::operator=(const CastorElectronicsMap& rhs) {
  CastorElectronicsMap temp(rhs);
  temp.swap(*this);
  return *this;
}
// public swap function
void CastorElectronicsMap::swap(CastorElectronicsMap& other) {
  std::swap(mPItems, other.mPItems);
  std::swap(mTItems, other.mTItems);
  other.mTItemsByTrigId.exchange(mTItemsByTrigId.exchange(other.mTItemsByTrigId));
  other.mPItemsById.exchange(mPItemsById.exchange(other.mPItemsById));
}
// move constructor
CastorElectronicsMap::CastorElectronicsMap(CastorElectronicsMap&& other) : CastorElectronicsMap() { other.swap(*this); }

const CastorElectronicsMap::PrecisionItem* CastorElectronicsMap::findById(unsigned long fId) const {
  PrecisionItem target(fId, 0);
  std::vector<const CastorElectronicsMap::PrecisionItem*>::const_iterator item;

  sortById();

  item = std::lower_bound((*mPItemsById).begin(), (*mPItemsById).end(), &target, castor_impl::LessById());
  if (item == (*mPItemsById).end() || (*item)->mId != fId)
    //    throw cms::Exception ("Conditions not found") << "Unavailable Electronics map for cell " << fId;
    return nullptr;
  return *item;
}

const CastorElectronicsMap::PrecisionItem* CastorElectronicsMap::findPByElId(unsigned long fElId) const {
  CastorElectronicsId eid(fElId);
  const PrecisionItem* i = &(mPItems[eid.linearIndex()]);

  if (i != nullptr && i->mElId != fElId)
    i = nullptr;
  return i;
}

const CastorElectronicsMap::TriggerItem* CastorElectronicsMap::findTByElId(unsigned long fElId) const {
  CastorElectronicsId eid(fElId);
  const TriggerItem* i = &(mTItems[eid.linearIndex()]);

  if (i != nullptr && i->mElId != fElId)
    i = nullptr;
  return i;
}

const CastorElectronicsMap::TriggerItem* CastorElectronicsMap::findByTrigId(unsigned long fTrigId) const {
  TriggerItem target(fTrigId, 0);
  std::vector<const CastorElectronicsMap::TriggerItem*>::const_iterator item;

  sortByTriggerId();

  item = std::lower_bound((*mTItemsByTrigId).begin(), (*mTItemsByTrigId).end(), &target, castor_impl::LessByTrigId());
  if (item == (*mTItemsByTrigId).end() || (*item)->mTrigId != fTrigId)
    //    throw cms::Exception ("Conditions not found") << "Unavailable Electronics map for cell " << fId;
    return nullptr;
  return *item;
}

const DetId CastorElectronicsMap::lookup(CastorElectronicsId fId) const {
  const PrecisionItem* item = findPByElId(fId.rawId());
  return DetId(item ? item->mId : 0);
}

const CastorElectronicsId CastorElectronicsMap::lookup(DetId fId) const {
  const PrecisionItem* item = findById(fId.rawId());
  return CastorElectronicsId(item ? item->mElId : 0);
}

const DetId CastorElectronicsMap::lookupTrigger(CastorElectronicsId fId) const {
  const TriggerItem* item = findTByElId(fId.rawId());
  return DetId(item ? item->mTrigId : 0);
}

const CastorElectronicsId CastorElectronicsMap::lookupTrigger(DetId fId) const {
  const TriggerItem* item = findByTrigId(fId.rawId());
  return CastorElectronicsId(item ? item->mElId : 0);
}

bool CastorElectronicsMap::lookup(const CastorElectronicsId pid,
                                  CastorElectronicsId& eid,
                                  HcalGenericDetId& did) const {
  const PrecisionItem* i = &(mPItems[pid.linearIndex()]);
  if (i != nullptr && i->mId != 0) {
    eid = CastorElectronicsId(i->mElId);
    did = HcalGenericDetId(i->mId);
    return true;
  } else
    return false;
}

bool CastorElectronicsMap::lookup(const CastorElectronicsId pid,
                                  CastorElectronicsId& eid,
                                  HcalTrigTowerDetId& did) const {
  const TriggerItem* i = &(mTItems[pid.linearIndex()]);
  if (i != nullptr && i->mTrigId != 0) {
    eid = CastorElectronicsId(i->mElId);
    did = HcalGenericDetId(i->mTrigId);
    return true;
  } else
    return false;
}

std::vector<CastorElectronicsId> CastorElectronicsMap::allElectronicsId() const {
  std::vector<CastorElectronicsId> result;
  for (std::vector<PrecisionItem>::const_iterator item = mPItems.begin(); item != mPItems.end(); item++)
    if (item->mElId)
      result.push_back(CastorElectronicsId(item->mElId));
  for (std::vector<TriggerItem>::const_iterator item = mTItems.begin(); item != mTItems.end(); item++)
    if (item->mElId)
      result.push_back(CastorElectronicsId(item->mElId));

  return result;
}

std::vector<CastorElectronicsId> CastorElectronicsMap::allElectronicsIdPrecision() const {
  std::vector<CastorElectronicsId> result;
  for (std::vector<PrecisionItem>::const_iterator item = mPItems.begin(); item != mPItems.end(); item++)
    if (item->mElId)
      result.push_back(CastorElectronicsId(item->mElId));
  return result;
}

std::vector<CastorElectronicsId> CastorElectronicsMap::allElectronicsIdTrigger() const {
  std::vector<CastorElectronicsId> result;
  for (std::vector<TriggerItem>::const_iterator item = mTItems.begin(); item != mTItems.end(); item++)
    if (item->mElId)
      result.push_back(CastorElectronicsId(item->mElId));

  return result;
}

std::vector<HcalGenericDetId> CastorElectronicsMap::allPrecisionId() const {
  std::vector<HcalGenericDetId> result;
  std::set<unsigned long> allIds;
  for (std::vector<PrecisionItem>::const_iterator item = mPItems.begin(); item != mPItems.end(); item++)
    if (item->mId)
      allIds.insert(item->mId);
  for (std::set<unsigned long>::const_iterator channel = allIds.begin(); channel != allIds.end(); channel++) {
    result.push_back(HcalGenericDetId(*channel));
  }
  return result;
}

std::vector<HcalTrigTowerDetId> CastorElectronicsMap::allTriggerId() const {
  std::vector<HcalTrigTowerDetId> result;
  std::set<unsigned long> allIds;
  for (std::vector<TriggerItem>::const_iterator item = mTItems.begin(); item != mTItems.end(); item++)
    if (item->mTrigId)
      allIds.insert(item->mTrigId);
  for (std::set<unsigned long>::const_iterator channel = allIds.begin(); channel != allIds.end(); channel++)
    result.push_back(HcalTrigTowerDetId(*channel));
  return result;
}

bool CastorElectronicsMap::mapEId2tId(CastorElectronicsId fElectronicsId, HcalTrigTowerDetId fTriggerId) {
  TriggerItem& item = mTItems[fElectronicsId.linearIndex()];
  if (item.mElId == 0)
    item.mElId = fElectronicsId.rawId();
  if (item.mTrigId == 0) {
    item.mTrigId = fTriggerId.rawId();  // just cast avoiding long machinery
  } else if (item.mTrigId != fTriggerId.rawId()) {
    edm::LogWarning("CASTOR") << "CastorElectronicsMap::mapEId2tId-> Electronics channel " << fElectronicsId
                              << " already mapped to trigger channel " << (HcalTrigTowerDetId(item.mTrigId))
                              << ". New value " << fTriggerId << " is ignored";
    return false;
  }
  return true;
}

bool CastorElectronicsMap::mapEId2chId(CastorElectronicsId fElectronicsId, DetId fId) {
  PrecisionItem& item = mPItems[fElectronicsId.linearIndex()];

  if (item.mElId == 0)
    item.mElId = fElectronicsId.rawId();
  if (item.mId == 0) {
    item.mId = fId.rawId();
  } else if (item.mId != fId.rawId()) {
    edm::LogWarning("CASTOR") << "CastorElectronicsMap::mapEId2tId-> Electronics channel " << fElectronicsId
                              << " already mapped to channel " << HcalGenericDetId(item.mId) << ". New value "
                              << HcalGenericDetId(fId) << " is ignored";
    return false;
  }
  return true;
}

void CastorElectronicsMap::sortById() const {
  if (!mPItemsById) {
    auto ptr = new std::vector<const PrecisionItem*>;
    for (auto i = mPItems.begin(); i != mPItems.end(); ++i) {
      if (i->mElId)
        (*ptr).push_back(&(*i));
    }
    std::sort((*ptr).begin(), (*ptr).end(), castor_impl::LessById());
    //atomically try to swap this to become mPItemsById
    std::vector<const PrecisionItem*>* expect = nullptr;
    bool exchanged = mPItemsById.compare_exchange_strong(expect, ptr);
    if (!exchanged) {
      delete ptr;
    }
  }
}

void CastorElectronicsMap::sortByTriggerId() const {
  if (!mTItemsByTrigId) {
    auto ptr = new std::vector<const TriggerItem*>;
    for (auto i = mTItems.begin(); i != mTItems.end(); ++i) {
      if (i->mElId)
        (*ptr).push_back(&(*i));
    }

    std::sort((*ptr).begin(), (*ptr).end(), castor_impl::LessByTrigId());
    //atomically try to swap this to become mTItemsByTrigId
    std::vector<const TriggerItem*>* expect = nullptr;
    bool exchanged = mTItemsByTrigId.compare_exchange_strong(expect, ptr);
    if (!exchanged) {
      delete ptr;
    }
  }
}