FourPointPlaneBounds.cc

Go to the documentation of this file.

// #include "Utilities/Configuration/interface/Architecture.h"

#include "MagneticField/VolumeGeometry/interface/FourPointPlaneBounds.h"
#include "DataFormats/GeometrySurface/interface/GeomExceptions.h"
#include <algorithm>
#include <iostream>

FourPointPlaneBounds::FourPointPlaneBounds( const LocalPoint& a, const LocalPoint& b, 
                        const LocalPoint& c, const LocalPoint& d) {
    corners_[0] = Local2DPoint( a.x(), a.y());
    corners_[1] = Local2DPoint( b.x(), b.y());
    corners_[2] = Local2DPoint( c.x(), c.y());
    corners_[3] = Local2DPoint( d.x(), d.y());

// check for convexity
    for (int i=0; i<4; ++i) {
    if (checkSide( i, corner(i+2)) *  checkSide( i, corner(i+3)) < 0) { // not on same side
        throw GeometryError("FourPointPlaneBounds: coners not in order or not convex");
    }
    }

    double side = checkSide( 0, corners_[2]); // - for clockwise corners, + for counterclockwise
    if (side < 0) {
    std::cout << "FourPointPlaneBounds: Changing order of corners to counterclockwise" << std::endl;
    std::swap( corners_[1], corners_[3]);
    }
}

double FourPointPlaneBounds::checkSide( int i, const Local2DPoint& lp) const {
    return checkSide(i, lp.x(), lp.y());
}

bool FourPointPlaneBounds::inside( const Local3DPoint& lp) const
{
    for (int i=0; i<4; ++i) {
    if (checkSide(i,lp.x(),lp.y()) < 0) return false;
    }
    return true;
}

float FourPointPlaneBounds::length() const {return 0;}
float FourPointPlaneBounds::width() const {return 0;}
float FourPointPlaneBounds::thickness() const {return 0;}