EnergyMomentumOperations.hpp

Go to the documentation of this file.

/*
 * (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/core/PhysicalUnits.hpp>
#include <corsika/framework/core/PhysicalGeometry.hpp>
#include <corsika/framework/core/PhysicalConstants.hpp>

namespace corsika {

  namespace detail {
    using HEPEnergyTypeSqr = decltype(1_GeV * 1_GeV);
  }

  auto constexpr calculate_mass_sqr(HEPEnergyType const E, HEPMomentumType const p) {
    return (E - p) * (E + p);
  }

  HEPEnergyType constexpr calculate_mass(HEPEnergyType const E, HEPMomentumType const p) {
    return sqrt(calculate_mass_sqr(E, p));
  }

  auto constexpr calculate_momentum_sqr(HEPEnergyType const E, HEPMassType const m) {
    return (E - m) * (E + m);
  }

  HEPEnergyType constexpr calculate_momentum(HEPEnergyType const E, HEPMassType const m) {
    return sqrt(calculate_momentum_sqr(E, m));
  }

  auto constexpr calculate_total_energy_sqr(HEPMomentumType const p,
                                            HEPMassType const m) {
    return p * p + m * m;
  }

  HEPEnergyType constexpr calculate_total_energy(HEPMomentumType const p,
                                                 HEPMassType const m) {
    return sqrt(calculate_total_energy_sqr(p, m));
  }

  HEPEnergyType constexpr calculate_kinetic_energy(HEPMomentumType const p,
                                                   HEPMassType const m) {
    return calculate_total_energy(p, m) - m;
  }

  HEPEnergyType constexpr calculate_lab_energy(detail::HEPEnergyTypeSqr sqrtS_sqr,
                                               HEPMassType const m_proj,
                                               HEPMassType const m_targ) {
    return (sqrtS_sqr - static_pow<2>(m_proj) - static_pow<2>(m_targ)) / (2 * m_targ);
  }

  HEPEnergyType constexpr calculate_com_energy(HEPEnergyType Elab,
                                               HEPMassType const m_proj,
                                               HEPMassType const m_targ) {
    return sqrt(2 * Elab * m_targ + static_pow<2>(m_proj) + static_pow<2>(m_targ));
  }

} // namespace corsika