EcalPreshowerTopology.h

Go to the documentation of this file.

#ifndef GEOMETRY_CALOTOPOLOGY_ECALPRESHOWERTOPOLOGY_H
#define GEOMETRY_CALOTOPOLOGY_ECALPRESHOWERTOPOLOGY_H 1

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

class EcalPreshowerTopology final : public CaloSubdetectorTopology {

 public:
  EcalPreshowerTopology() : theGeom_(0) {};

  virtual ~EcalPreshowerTopology() { }  
  
  EcalPreshowerTopology(edm::ESHandle<CaloGeometry> theGeom) : theGeom_(theGeom)
    {
    }

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

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

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

  virtual DetId goWest(const DetId& id) const {
    return decrementIx(ESDetId(id));
  }
  virtual std::vector<DetId> west(const DetId& id) const
  { 
    ESDetId nextId=goWest(id);
    std::vector<DetId> vNeighborsDetId;
    if (! (nextId==ESDetId(0)))
      vNeighborsDetId.push_back(DetId(nextId.rawId()));
    return vNeighborsDetId;
  }
  
  virtual DetId goUp(const DetId& id) const {
    return incrementIz(ESDetId(id));
  }
  virtual std::vector<DetId> up(const DetId& id) const
  {
    ESDetId nextId=goUp(id);
    std::vector<DetId> vNeighborsDetId;
    if (! (nextId==ESDetId(0)))
      vNeighborsDetId.push_back(DetId(nextId.rawId()));
    return  vNeighborsDetId;
  }
  
  virtual DetId goDown(const DetId& id) const {
    return decrementIz(ESDetId(id));
  }
  virtual std::vector<DetId> down(const DetId& id) const
  {
    ESDetId nextId=goDown(id);
    std::vector<DetId> vNeighborsDetId;
    if (! (nextId==ESDetId(0)))
      vNeighborsDetId.push_back(DetId(nextId.rawId()));
    return  vNeighborsDetId;
  }

 private:

  ESDetId incrementIx(const ESDetId& id) const ;

  ESDetId decrementIx(const ESDetId& id) const ;

  ESDetId incrementIy(const ESDetId& id) const ;

  ESDetId decrementIy(const ESDetId& id) const;

  ESDetId incrementIz(const ESDetId& id) const;

  ESDetId decrementIz(const ESDetId& id) const;

  edm::ESHandle<CaloGeometry> theGeom_;
};

#endif