RotationMatrix.hh

//
// This file is part of MARLEY (Model of Argon Reaction Low Energy Yields)
//
// MARLEY is free software: you can redistribute it and/or modify it under the
// terms of version 3 of the GNU General Public License as published by the
// Free Software Foundation.
//
// For the full text of the license please see COPYING or
// visit http://opensource.org/licenses/GPL-3.0
//
// Please respect the MCnet academic usage guidelines. See GUIDELINES
// or visit https://www.montecarlonet.org/GUIDELINES for details.
 
#pragma once
 
// Standard library includes
#include <array>
 
// MARLEY includes
#include "marley/Particle.hh"
 
namespace marley {
 
  class RotationMatrix {
 
    using ThreeVector = std::array<double, 3>;
    using ThreeThreeMatrix = std::array<ThreeVector, 3>;
 
    public:
 
      RotationMatrix();
 
      RotationMatrix(const ThreeVector& from_vec, const ThreeVector& to_vec);
 
      ThreeVector rotate_copy(const ThreeVector& v);
 
      // Returns a copy of the 3-vector v normalized to have unit magnitude
      static ThreeVector normalize(const ThreeVector& v);
 
      void rotate_inplace(ThreeVector& v);
 
      void rotate_particle_inplace(marley::Particle& p);
 
    protected:
 
      ThreeThreeMatrix matrix_;
  };
 
}