/data/refman/pasoursint/CMSSW_4_2_9_HLT1_bphpatch4/src/Geometry/CaloTopology/interface/EcalEndcapHardcodedTopology.h

Go to the documentation of this file.

#ifndef GEOMETRY_CALOTOPOLOGY_ECALENDCAPHARDCODEDTOPOLOGY_H
#define GEOMETRY_CALOTOPOLOGY_ECALENDCAPHARDCODEDTOPOLOGY_H 1

#include <vector>
#include <iostream>
#include "DataFormats/EcalDetId/interface/EEDetId.h"
#include "Geometry/CaloTopology/interface/CaloSubdetectorTopology.h"

class EcalEndcapHardcodedTopology : public CaloSubdetectorTopology
{

 public:
  EcalEndcapHardcodedTopology() {};

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

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

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

  virtual std::vector<DetId> west(const DetId& id) const
    { 
      EEDetId nextId=decrementIx(EEDetId(id));
      std::vector<DetId> vNeighborsDetId;
      if (! (nextId==EEDetId(0)))
        vNeighborsDetId.push_back(DetId(nextId.rawId()));
      return vNeighborsDetId;
    }
  
  virtual std::vector<DetId> up(const DetId& id) const
    {
      std::cout << "EcalEndcapHardcodedTopology::up() not yet implemented" << std::endl; 
      std::vector<DetId> vNeighborsDetId;
      return  vNeighborsDetId;
    }
  
  virtual std::vector<DetId> down(const DetId& id) const
    {
      std::cout << "EcalEndcapHardcodedTopology::down() not yet implemented" << std::endl; 
      std::vector<DetId> vNeighborsDetId;
      return  vNeighborsDetId;
    }

 private:

  EEDetId incrementIx(const EEDetId&) const ;

  EEDetId decrementIx(const EEDetId&) const ;

  EEDetId incrementIy(const EEDetId&) const ;

  EEDetId decrementIy(const EEDetId&) const;

};

#endif