Bounds.h

Go to the documentation of this file.

#ifndef Geom_Bounds_H
#define Geom_Bounds_H

#include "DataFormats/GeometryVector/interface/LocalPoint.h"
#include "DataFormats/GeometrySurface/interface/LocalError.h"

#include "DataFormats/GeometrySurface/interface/BoundSpan.h"

class Bounds {
public:
  virtual ~Bounds() {}

  virtual float length() const = 0;

  virtual float width() const = 0;

  virtual float thickness() const = 0;

  virtual float widthAtHalfLength() const { return width(); }

  virtual bool inside(const Local3DPoint&) const = 0;
  virtual bool inside(const Local2DPoint& p) const { return inside(Local3DPoint(p.x(), p.y(), 0)); }

  virtual bool inside(const Local2DPoint& p, float tollerance) const {
    return inside(p, LocalError(tollerance * tollerance, 0.f, tollerance * tollerance));
  }

  // For propagation with uncertainties - one has to know by how
  // much one missed a surface
  // virtual Local2DVector<float> howFar( const Local2DPoint&) = 0;
  //
  // Or maybe a better way of determining if a surface is "touched"
  // by a trajectory is to have a method

  virtual bool inside(const Local3DPoint&, const LocalError&, float scale = 1.f) const = 0;

  virtual bool inside(const Local2DPoint& p, const LocalError& err, float scale = 1.f) const {
    return inside(Local3DPoint(p.x(), p.y(), 0), err, scale);
  }

  virtual float significanceInside(const Local3DPoint&, const LocalError&) const;

  virtual Bounds* clone() const = 0;

  std::pair<float, float> const& phiSpan() const { return m_span.phiSpan(); }
  std::pair<float, float> const& zSpan() const { return m_span.zSpan(); }
  std::pair<float, float> const& rSpan() const { return m_span.rSpan(); }

  void computeSpan(Surface const& plane) { m_span.compute(plane); }

private:
  BoundSpan m_span;
};

#endif