EcalEndcapTopology.h

Go to the documentation of this file.

#ifndef GEOMETRY_CALOTOPOLOGY_ECALENDCAPTOPOLOGY_H
#define GEOMETRY_CALOTOPOLOGY_ECALENDCAPTOPOLOGY_H 1

#include "DataFormats/EcalDetId/interface/EEDetId.h"
#include "Geometry/CaloTopology/interface/CaloSubdetectorTopology.h"
#include "Geometry/CaloGeometry/interface/CaloGeometry.h"
#include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include <utility>
#include <vector>
#include <iostream>

class EcalEndcapTopology final : public CaloSubdetectorTopology {

public:
  EcalEndcapTopology() : theGeom_(nullptr) {};

  ~EcalEndcapTopology() override { }  
  
  EcalEndcapTopology(CaloGeometry const& theGeom) {
    theGeom_ = theGeom.getSubdetectorGeometry(DetId::Ecal,EcalEndcap);
  }

  DetId  goNorth(const DetId& id) const override {
    return incrementIy(EEDetId(id));
  }
  std::vector<DetId> north(const DetId& id) const override
    { 
      EEDetId nextId= goNorth(id);
      std::vector<DetId> vNeighborsDetId;
      if (! (nextId==EEDetId(0)))
    vNeighborsDetId.emplace_back(DetId(nextId.rawId()));
      return vNeighborsDetId;
    }

  DetId goSouth(const DetId& id) const override {
    return decrementIy(EEDetId(id));
  }
  std::vector<DetId> south(const DetId& id) const override
    { 
      EEDetId nextId= goSouth(id);
      std::vector<DetId> vNeighborsDetId;
      if (! (nextId==EEDetId(0)))
    vNeighborsDetId.emplace_back(DetId(nextId.rawId()));
      return vNeighborsDetId;
    }

  DetId  goEast(const DetId& id) const override {
    return incrementIx(EEDetId(id));
  }
  std::vector<DetId> east(const DetId& id) const override
    { 
      EEDetId nextId=goEast(id);
      std::vector<DetId> vNeighborsDetId;
      if (! (nextId==EEDetId(0)))
    vNeighborsDetId.emplace_back(DetId(nextId.rawId()));
      return vNeighborsDetId;
    }

  DetId goWest(const DetId& id) const override {
    return decrementIx(EEDetId(id));
  }
  std::vector<DetId> west(const DetId& id) const override
    { 
      EEDetId nextId=goWest(id);
      std::vector<DetId> vNeighborsDetId;
      if (! (nextId==EEDetId(0)))
    vNeighborsDetId.emplace_back(DetId(nextId.rawId()));
      return vNeighborsDetId;
    }
  
  std::vector<DetId> up(const DetId& /*id*/) const override
    {
      std::cout << "EcalBarrelTopology::up() not yet implemented" << std::endl; 
      std::vector<DetId> vNeighborsDetId;
      return  vNeighborsDetId;
    }
  
  std::vector<DetId> down(const DetId& /*id*/) const override
    {
      std::cout << "EcalBarrelTopology::down() not yet implemented" << std::endl; 
      std::vector<DetId> vNeighborsDetId;
      return  vNeighborsDetId;
    }

 private:

  EEDetId incrementIx(const EEDetId& id) const ;

  EEDetId decrementIx(const EEDetId& id) const ;

  EEDetId incrementIy(const EEDetId& id) const ;

  EEDetId decrementIy(const EEDetId& id) const;

  const CaloSubdetectorGeometry* theGeom_;
};

#endif