PhiBorderFinder.h

Go to the documentation of this file.

#ifndef TrackingTools_DetLayers_PhiBorderFinder_H
#define TrackingTools_DetLayers_PhiBorderFinder_H
 
#include <Geometry/CommonDetUnit/interface/GeomDet.h>
 
#include <DataFormats/GeometrySurface/interface/BoundingBox.h>
#include <DataFormats/GeometrySurface/interface/GeometricSorting.h>
 
#include <DataFormats/GeometryVector/interface/Pi.h>
#include <Utilities/General/interface/precomputed_value_sort.h>
#include <TrackingTools/DetLayers/interface/simple_stat.h>
#include <TrackingTools/DetLayers/interface/DetRod.h>
#include <FWCore/Utilities/interface/Exception.h>
 
// FIXME: remove this include
#include "FWCore/MessageLogger/interface/MessageLogger.h"
 
#include <vector>
 
class PhiBorderFinder {
public:
  typedef DetRod Det;  //FIXME!!!
  typedef geomsort::ExtractPhi<Det, float> DetPhi;
 
  PhiBorderFinder(const std::vector<const Det*>& utheDets);
 
  virtual ~PhiBorderFinder(){};
 
  inline unsigned int nBins() { return theNbins; }
 
  inline bool isPhiPeriodic() const { return isPhiPeriodic_; }
 
  inline bool isPhiOverlapping() const { return isPhiOverlapping_; }
 
  inline const std::vector<double>& phiBorders() const { return thePhiBorders; }
 
  inline const std::vector<double>& phiBins() const { return thePhiBins; }
 
  //  inline std::vector<double> etaBorders() {}
  //  inline std::vector<double> zBorders() {}
 
private:
  unsigned int theNbins;
  bool isPhiPeriodic_;
  bool isPhiOverlapping_;
  std::vector<double> thePhiBorders;
  std::vector<double> thePhiBins;
 
  inline double positiveRange(double phi) const { return (phi > 0) ? phi : phi + Geom::twoPi(); }
 
  int binIndex(int i) const {
    int ind = i % (int)theNbins;
    return (ind < 0) ? ind + theNbins : ind;
  }
};
#endif