FourPointPlaneBounds.h

Go to the documentation of this file.

#ifndef FourPointPlaneBounds_H
#define FourPointPlaneBounds_H

#include "DataFormats/GeometrySurface/interface/Bounds.h"
#include "DataFormats/GeometryVector/interface/LocalVector.h"
#include "DataFormats/GeometryVector/interface/Vector2DBase.h"

typedef Vector2DBase< float, LocalTag>    Local2DVector;

class FourPointPlaneBounds : public Bounds {
public:

    typedef LocalPoint::ScalarType    Scalar;

    FourPointPlaneBounds( const LocalPoint& a, const LocalPoint& b, 
              const LocalPoint& c, const LocalPoint& d); 
    ~FourPointPlaneBounds() override {}

  float length() const override;
  float width() const override;
  float thickness() const override;

    // basic bounds function
    bool inside( const Local3DPoint& lp) const override;
    
    bool inside( const Local3DPoint& lp , const LocalError& e, float scale) const override {
      return inside( lp);
    }

    Bounds* clone() const override {return new FourPointPlaneBounds(*this);}

private:

    Local2DPoint corners_[4];

    const Local2DPoint& corner(int i) const {return corners_[i%4];}
    double checkSide( int i, const Local2DPoint& lp) const;

    double checkSide( int i, Scalar x, Scalar y) const {
    const Local2DPoint& cor( corner(i));
    Local2DVector v( corner(i+1) - cor);
    // Local2DVector normal( -v.y(), v.x());  //  90 deg rotated
    return -v.y() * (x-cor.x()) + v.x() * (y-cor.y()); // == normal.dot(LP(x,y)-cor))
    }

};

#endif