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(FourVector<HEPEnergyType, MomentumVector> const& Pprojectile,
             HEPEnergyType const massTarget);

    COMBoost(MomentumVector const& momentum, HEPEnergyType mass);

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

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

    CoordinateSystemPtr getRotatedCS() const;

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

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

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