GammaLib

#ifndef _TOOLS_POLYGON_H_

#define _TOOLS_POLYGON_H_

#include <vector>

#include <TGeoPolygon.h>

#include <stash/Coordinate.hh>

#include <stash/Lambda.hh>

#include <TArrayD.h>

#include "Region.hh"

namespace Tools {

  /**

   * A polygonal region, defined by a set of vertex coordinates (clockwise)

   */

  class Polygon : public Region

  {

  public:

    Polygon( std::vector<Stash::Coordinate> vertices, 

             std::string name="poly");

    Polygon(const Tools::Polygon &cpy);

    Polygon();

    virtual double GetArea() const;

    virtual bool IsInside( const Stash::Coordinate& dir ) const;

    virtual std::string AsFITSRegionString() const;

    virtual VertexList AsVertexList() const;

    virtual double GetBound1() const;

    virtual double GetBound2() const;

    virtual double GetRadialBound()const;

    virtual ~Polygon();

    virtual Tools::Polygon* GrowBoundary(const double degrees,

                                        const std::string name="_grown") const;

    virtual void MoveTo( Stash::Coordinate newpos );

    /**

     * Move the internal "center" position to the Center-of-gravity

     * (doesn't move the region itself, just the point from which

     * bounding boxes are calculated.) Use MoveTo() instead to move the region

     */

    void RecenterOnCOG() { fCentre = GetCentreOfGravity(); }

    int GetNumVertices() { return fPoly->GetNvert(); }

    TGeoPolygon* GetPoly() {return fPoly; }

  private:

    void UpdatePoly();

    TGeoPolygon* fPoly;  

    TArrayD *fX; // needed for TGeoPolygon

    TArrayD *fY; // needed for TGeoPolygon

    ClassDef(Tools::Polygon, 4);

  };

}

#endif /* _TOOLS_POLYGON_H_ */