GeometricSorting.h

Go to the documentation of this file.

#ifndef Surface_GeometricSorting_h
#define Surface_GeometricSorting_h

#include <functional>
#include <DataFormats/GeometryVector/interface/Phi.h>

namespace geomsort{

  template <class T, class Scalar = typename T::Scalar>
  struct ExtractR {
    typedef Scalar result_type;
    Scalar operator()(const T* p) const {return p->position().perp();}
    Scalar operator()(const T& p) const {return p.position().perp();}
  };


  template <class T, class Scalar = typename T::Scalar>
  struct ExtractPhi {
    typedef Geom::Phi<Scalar> result_type;
    Geom::Phi<Scalar> operator()(const T* p) const {return p->position().phi();}
    Geom::Phi<Scalar> operator()(const T& p) const {return p.position().phi();}
  };



  template <class T, class Scalar = typename T::Scalar>
  struct ExtractZ {
    typedef Scalar result_type;
    Scalar operator()(const T* p) const {return p->position().z();}
    Scalar operator()(const T& p) const {return p.position().z();}
  };

  template <class T, class Scalar = typename T::Scalar>
  struct ExtractAbsZ {
    typedef Scalar result_type;
    Scalar operator()(const T* p) const {return fabs(p->position().z());}
    Scalar operator()(const T& p) const {return fabs(p.position().z());}
  };


}
#endif