CaloTopology.cc

Go to the documentation of this file.

#include "Geometry/CaloTopology/interface/CaloSubdetectorTopology.h"
#include "Geometry/CaloTopology/interface/CaloTopology.h"


CaloTopology::CaloTopology() {
}

CaloTopology::~CaloTopology() = default;


int CaloTopology::makeIndex(DetId::Detector det, int subdet) const {
  return (int(det)<<4) | (subdet&0xF);
}

void CaloTopology::setSubdetTopology(DetId::Detector det, int subdet, std::unique_ptr<const CaloSubdetectorTopology> geom) {
  int index=makeIndex(det,subdet);
  theTopologies_[index]=std::move(geom);
}

const CaloSubdetectorTopology* CaloTopology::getSubdetectorTopology(const DetId& id) const {
  auto i=theTopologies_.find(makeIndex(id.det(),id.subdetId()));
  return (i==theTopologies_.end())?(nullptr):(i->second.get());
}

const CaloSubdetectorTopology* CaloTopology::getSubdetectorTopology(DetId::Detector det, int subdet) const {
    auto i=theTopologies_.find(makeIndex(det,subdet));
    return (i==theTopologies_.end())?(nullptr):(i->second.get());
}

static const std::vector<DetId> emptyDetIdVector;

std::vector<DetId> CaloTopology::east(const DetId& id) const {
  const CaloSubdetectorTopology* topology=getSubdetectorTopology(id);
  return (topology==nullptr) ? (emptyDetIdVector):(topology->east(id));
}

std::vector<DetId> CaloTopology::west(const DetId& id) const {
  const CaloSubdetectorTopology* topology=getSubdetectorTopology(id);
  return (topology==nullptr) ? (emptyDetIdVector):(topology->west(id));
}

std::vector<DetId> CaloTopology::north(const DetId& id) const {
  const CaloSubdetectorTopology* topology=getSubdetectorTopology(id);
  return (topology==nullptr) ? (emptyDetIdVector):(topology->north(id));
}

std::vector<DetId> CaloTopology::south(const DetId& id) const {
  const CaloSubdetectorTopology* topology=getSubdetectorTopology(id);
  return (topology==nullptr) ? (emptyDetIdVector):(topology->south(id));
}

std::vector<DetId> CaloTopology::up(const DetId& id) const {
  const CaloSubdetectorTopology* topology=getSubdetectorTopology(id);
  return (topology==nullptr) ? (emptyDetIdVector):(topology->up(id));
}

std::vector<DetId> CaloTopology::down(const DetId& id) const {
    const CaloSubdetectorTopology* topology=getSubdetectorTopology(id);
    return (topology==nullptr) ? (emptyDetIdVector):(topology->down(id));
}

std::vector<DetId> CaloTopology::getNeighbours(const DetId& id,const CaloDirection& dir) const {
    const CaloSubdetectorTopology* topology=getSubdetectorTopology(id);
    return (topology==nullptr) ? (emptyDetIdVector):(topology->getNeighbours(id,dir));
}

std::vector<DetId> CaloTopology::getWindow(const DetId& id, const int& northSouthSize, const int& eastWestSize) const {
    const CaloSubdetectorTopology* topology=getSubdetectorTopology(id);
    return (topology==nullptr) ? (emptyDetIdVector):(topology->getWindow(id,northSouthSize, eastWestSize));
}

std::vector<DetId> CaloTopology::getAllNeighbours(const DetId& id) const {
    const CaloSubdetectorTopology* topology=getSubdetectorTopology(id);
    return (topology==nullptr) ? (emptyDetIdVector):(topology->getAllNeighbours(id));
}

bool CaloTopology::valid(const DetId& id) const {
  const CaloSubdetectorTopology* geom=getSubdetectorTopology(id);
  return (geom==nullptr)?(false):(geom->valid(id));
}