EcalBarrelHardcodedTopology.h

Go to the documentation of this file.

#ifndef GEOMETRY_CALOTOPOLOGY_ECALBARRELHARDCODEDTOPOLOGY_H
#define GEOMETRY_CALOTOPOLOGY_ECALBARRELHARDCODEDTOPOLOGY_H 1
 
#include <vector>
#include <iostream>
#include "DataFormats/EcalDetId/interface/EBDetId.h"
#include "Geometry/CaloTopology/interface/CaloSubdetectorTopology.h"
 
class EcalBarrelHardcodedTopology final : public CaloSubdetectorTopology {
public:
  EcalBarrelHardcodedTopology(){};
 
  ~EcalBarrelHardcodedTopology() override{};
 
  DetId goNorth(const DetId& id) const override { return incrementIphi(EBDetId(id)); }
  std::vector<DetId> north(const DetId& id) const override {
    EBDetId nextId = goNorth(id);
    std::vector<DetId> vNeighborsDetId;
    if (!(nextId == EBDetId(0)))
      vNeighborsDetId.emplace_back(DetId(nextId.rawId()));
    return vNeighborsDetId;
  }
 
  DetId goSouth(const DetId& id) const override { return decrementIphi(EBDetId(id)); }
  std::vector<DetId> south(const DetId& id) const override {
    EBDetId nextId = goSouth(id);
    std::vector<DetId> vNeighborsDetId;
    if (!(nextId == EBDetId(0)))
      vNeighborsDetId.emplace_back(DetId(nextId.rawId()));
    return vNeighborsDetId;
  }
 
  DetId goEast(const DetId& id) const override { return decrementIeta(EBDetId(id)); }
  std::vector<DetId> east(const DetId& id) const override {
    EBDetId nextId = goEast(id);
    std::vector<DetId> vNeighborsDetId;
    if (!(nextId == EBDetId(0)))
      vNeighborsDetId.emplace_back(DetId(nextId.rawId()));
    return vNeighborsDetId;
  }
 
  DetId goWest(const DetId& id) const override { return incrementIeta(EBDetId(id)); }
  std::vector<DetId> west(const DetId& id) const override {
    EBDetId nextId = goWest(id);
    std::vector<DetId> vNeighborsDetId;
    if (!(nextId == EBDetId(0)))
      vNeighborsDetId.emplace_back(DetId(nextId.rawId()));
    return vNeighborsDetId;
  }
 
  std::vector<DetId> up(const DetId& /*id*/) const override {
    std::cout << "EcalBarrelHardcodedTopology::up() not yet implemented" << std::endl;
    std::vector<DetId> vNeighborsDetId;
    return vNeighborsDetId;
  }
 
  std::vector<DetId> down(const DetId& /*id*/) const override {
    std::cout << "EcalBarrelHardcodedTopology::down() not yet implemented" << std::endl;
    std::vector<DetId> vNeighborsDetId;
    return vNeighborsDetId;
  }
 
private:
  EBDetId incrementIeta(const EBDetId&) const;
 
  EBDetId decrementIeta(const EBDetId&) const;
 
  EBDetId incrementIphi(const EBDetId&) const;
 
  EBDetId decrementIphi(const EBDetId&) const;
};
 
#endif