/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_13_patch3/src/RecoCaloTools/Navigation/interface/CaloNavigator.h

Go to the documentation of this file.

#ifndef RECOCALOTOOLS_NAVIGATION_CALONAVIGATOR_H
#define RECOCALOTOOLS_NAVIGATION_CALONAVIGATOR_H 1

#include "Geometry/CaloTopology/interface/CaloSubdetectorTopology.h"

template <class T>
class CaloNavigator {
 public:

  CaloNavigator() : myTopology_(0)
    {
    }

  CaloNavigator(const T& home) : myTopology_(0)
    {
      setHome(home);
    };


  CaloNavigator(const T& home, const CaloSubdetectorTopology* topology) : myTopology_(topology) 
    {
      setHome(home);
    };

  virtual ~CaloNavigator(){};

  void setHome(const T& startingPoint);

  void setTopology(const CaloSubdetectorTopology*);

  const CaloSubdetectorTopology* getTopology() const
    {
      return myTopology_;
    }

  void home() const ;

  T pos() const { return currentPoint_; }

  T operator*() const { return currentPoint_; } 

  virtual T north() const 
    { 
      if ((myTopology_->north(currentPoint_)).size()==1)
        currentPoint_=(myTopology_->north(currentPoint_))[0];
      else
        currentPoint_=T(0);
      return currentPoint_;
    } ;

  virtual T south()  const 
    { 
      if ((myTopology_->south(currentPoint_)).size()==1)
        currentPoint_=(myTopology_->south(currentPoint_))[0];
      else
        currentPoint_=T(0);
      return currentPoint_;
    } ;

  virtual T east() const
    { 
      if ((myTopology_->east(currentPoint_)).size()==1)
        currentPoint_=(myTopology_->east(currentPoint_))[0];
      else
        currentPoint_=T(0);
      return currentPoint_;
    } ;

  virtual T west() const
    { 
      if ((myTopology_->west(currentPoint_)).size()==1)
        currentPoint_=(myTopology_->west(currentPoint_))[0];
      else
        currentPoint_=T(0);
      return currentPoint_;
    } ;

  virtual T up() const
    { 
      if ((myTopology_->up(currentPoint_)).size()==1)
        currentPoint_=(myTopology_->up(currentPoint_))[0];
      else
        currentPoint_=T(0);
      return currentPoint_;
    } ;

  virtual T down() const
    { 
      if ((myTopology_->down(currentPoint_)).size()==1)
        currentPoint_=(myTopology_->down(currentPoint_))[0];
      else
        currentPoint_=T(0);
      return currentPoint_;
    } ;

  virtual T offsetBy(int deltaX, int deltaY) const
    {
      for(int x=0; x < abs(deltaX) && currentPoint_ != T(0); x++)
        {
          if(deltaX > 0) east();
          else           west();
        }

      for(int y=0; y < abs(deltaY) && currentPoint_ != T(0); y++)
        {
          if(deltaY > 0) north();
          else           south();
        }

      return currentPoint_;

    }

 protected:
  
  const CaloSubdetectorTopology* myTopology_;
  mutable T startingPoint_, currentPoint_;
};

template <class T>
void CaloNavigator<T>::setHome(const T& startingPoint)
{
  startingPoint_=startingPoint;
  home();
}

template <class T>
void CaloNavigator<T>::home() const
{
  currentPoint_=startingPoint_;
}

template <class T>
void CaloNavigator<T>::setTopology(const CaloSubdetectorTopology* topology) 
{
  if (myTopology_ == 0)
    myTopology_=topology;
  else
    return;
}

#endif