BoundingBox.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2015 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Oliver Frost                                             *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
#pragma once
 
#include <algorithm>
#include <cmath>
 
namespace Belle2 {
  namespace TrackFindingCDC {
 
    class BoundingBox {
 
    public:
 
      BoundingBox() :
        m_left(NAN),
        m_bottom(NAN),
        m_right(NAN),
        m_top(NAN)
      {}
 
      BoundingBox(float x1, float y1, float x2, float y2) :
        m_left(std::min(x1, x2)),
        m_bottom(std::min(y1, y2)),
        m_right(std::max(x1, x2)),
        m_top(std::max(y1, y2))
      {}
 
      void operator&=(const BoundingBox& other)
      {
 
        m_left = std::isnan(getLeft()) ? other.getLeft() : std::min(getLeft(), other.getLeft());
        m_bottom = std::isnan(getBottom()) ? other.getBottom() : std::min(getBottom(), other.getBottom());
 
        m_right = std::isnan(getRight()) ? other.getRight() : std::max(getRight(), other.getRight());
        m_top = std::isnan(getTop()) ? other.getTop() : std::max(getTop(), other.getTop());
 
      }
 
      void expand(float delta)
      {
        m_left -= delta;
        m_bottom -= delta;
        m_right += delta;
        m_top += delta;
      }
 
      float getWidth() const
      { return getRight() - getLeft(); }
 
      float getHeight() const
      { return getTop() - getBottom(); }
 
      float getLeft() const
      { return m_left; }
 
      float getBottom() const
      { return m_bottom; }
 
      float getRight() const
      { return m_right; }
 
      float getTop() const
      { return m_top; }
 
      void clear()
      {
        m_left = NAN;
        m_bottom = NAN;
        m_right = NAN;
        m_top = NAN;
      }
 
    private:
      float m_left;
 
      float m_bottom;
 
      float m_right;
 
      float m_top;
 
    };
 
  }
}