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"

#include "FWCore/Utilities/interface/GCC11Compatibility.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 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