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>
 
namespace HepMC3 {
  class GenParticle;
}
 
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( HepMC3::GenParticle& p );
 
    protected:

      ThreeThreeMatrix matrix_;
  };
 
}