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() 
{
   for( TopMap::iterator i ( theTopologies_.begin() ) ;
    i != theTopologies_.end() ; ++i )
   {
      delete i->second ;
   }
}

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

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

const CaloSubdetectorTopology* CaloTopology::getSubdetectorTopology(const DetId& id) const {
  std::map<int, const CaloSubdetectorTopology*>::const_iterator i=theTopologies_.find(makeIndex(id.det(),id.subdetId()));
  return (i==theTopologies_.end())?(0):(i->second);
}

const CaloSubdetectorTopology* CaloTopology::getSubdetectorTopology(DetId::Detector det, int subdet) const {
    std::map<int, const CaloSubdetectorTopology*>::const_iterator i=theTopologies_.find(makeIndex(det,subdet));
    return (i==theTopologies_.end())?(0):(i->second);
}

static const std::vector<DetId> emptyDetIdVector;

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

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

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

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

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

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

std::vector<DetId> CaloTopology::getNeighbours(const DetId& id,const CaloDirection& dir) const {
    const CaloSubdetectorTopology* topology=getSubdetectorTopology(id);
    return (topology==0) ? (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==0) ? (emptyDetIdVector):(topology->getWindow(id,northSouthSize, eastWestSize));
}

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

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