BDSTwoVector.hh

/* 
Beam Delivery Simulation (BDSIM) Copyright (C) Royal Holloway, 
University of London 2001 - 2023.
 
This file is part of BDSIM.
 
BDSIM is free software: you can redistribute it and/or modify 
it under the terms of the GNU General Public License as published 
by the Free Software Foundation version 3 of the License.
 
BDSIM is distributed in the hope that it will be useful, but 
WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
GNU General Public License for more details.
 
You should have received a copy of the GNU General Public License
along with BDSIM.  If not, see <http://www.gnu.org/licenses/>.
*/
#ifndef BDSTWOVECTOR_H
#define BDSTWOVECTOR_H
#include "BDSException.hh"
 
#include "globals.hh" // geant4 types / globals
 
#include <ostream>
 
template <typename T>
class BDSTwoVector
{
public:
  BDSTwoVector()
  {values[0] = 0; values[1] = 0;}
 
  BDSTwoVector(T x, T y)
  {values[0] = x; values[1] = y;}
 
  BDSTwoVector(const BDSTwoVector& other)
  {
    for (G4int i : {0, 1})
      {values[i] = other.values[i];}
  }
 
  BDSTwoVector& operator=(const BDSTwoVector& rhs)
  {
    if (this == &rhs)
      {return *this;}
    values[0] = rhs.values[0];
    values[1] = rhs.values[1];
    return *this;
  }
 
  T& operator[](const G4int index)
  {
    if (index > 1)
      {throw BDSException("BDSTwoVector::operator[]> index outside array");}
    return values[index];
  }
 
  // Access a single element.
  const T& operator[](const G4int index) const
  {
    if (index > 1)
      {throw BDSException("BDSTwoVector::operator[]> index outside array");}
    return values[index];
  }
 
  inline const T& x() const {return values[0];}
  inline const T& y() const {return values[1];}
  
  friend std::ostream& operator<< (std::ostream& out, BDSTwoVector const &v)
  {out << "(" << v.values[0] << ", " << v.values[1] << ")"; return out;}
 
private:
  T values[2];
};
 
 
#endif