COMBoost.hpp

/*
 * (c) Copyright 2020 CORSIKA Project, corsika-project@lists.kit.edu
 *
 * This software is distributed under the terms of the GNU General Public
 * Licence version 3 (GPL Version 3). See file LICENSE for a full version of
 * the license.
 */

#pragma once

#include <corsika/framework/geometry/CoordinateSystem.hpp>
#include <corsika/framework/geometry/FourVector.hpp>
#include <corsika/framework/geometry/PhysicalGeometry.hpp>
#include <corsika/framework/core/PhysicalUnits.hpp>
#include <corsika/framework/core/Logging.hpp>

#include <Eigen/Dense>

namespace corsika {

  class COMBoost {

  public:
    COMBoost(FourMomentum const& P4projectile, HEPEnergyType const massTarget);

    COMBoost(MomentumVector const& momentum, HEPEnergyType const mass);

    COMBoost(FourMomentum const& P4projectile, FourMomentum const& P4target);

    template <typename FourVector>
    FourVector toCoM(FourVector const& p4) const;

    template <typename FourVector>
    FourVector fromCoM(FourVector const& p4) const;

    CoordinateSystemPtr const& getRotatedCS() const;

    CoordinateSystemPtr const& getOriginalCS() const;

  protected:
    void setBoost(double const coshEta, double const sinhEta);

  public:
    Eigen::Matrix2d boost_;
    Eigen::Matrix2d inverseBoost_;
    CoordinateSystemPtr const originalCS_;
    CoordinateSystemPtr rotatedCS_;
  };
} // namespace corsika

#include <corsika/detail/framework/utility/COMBoost.inl>