crosssectionsphoton.h

Go to the documentation of this file.

/*
 *
 *    Copyright (c) 2017-2023
 *      SMASH Team
 *
 *    GNU General Public License (GPLv3 or later)
 *
 */
 
#ifndef SRC_INCLUDE_SMASH_CROSSSECTIONSPHOTON_H_
#define SRC_INCLUDE_SMASH_CROSSSECTIONSPHOTON_H_
 
#include "kinematics.h"
 
namespace smash {
double cut_off(const double sigma_mb);
 
enum class ComputationMethod { Analytic };
 
template <ComputationMethod method>
class CrosssectionsPhoton {};
 
template <>
class CrosssectionsPhoton<ComputationMethod::Analytic> {
 public:
 
  static double xs_pi_pi_rho0(const double s, const double m_rho);
  static double xs_pi_pi0_rho(const double s, const double m_rho);
  static double xs_pi0_rho0_pi0(const double s, const double m_rho);
  static double xs_pi_rho0_pi(const double s, const double m_rho);
 
  static double xs_pi_rho_pi0(const double s, const double m_rho);
  static double xs_pi_rho_pi0_rho_mediated(const double s, const double m_rho);
  static double xs_pi_rho_pi0_omega_mediated(const double s,
                                             const double m_rho);
 
  static double xs_pi0_rho_pi(const double s, const double m_rho);
  static double xs_pi0_rho_pi_rho_mediated(const double s, const double m_rho);
  static double xs_pi0_rho_pi_omega_mediated(const double s,
                                             const double m_rho);
 
 
  static double xs_diff_pi_pi_rho0(const double s, const double t,
                                   const double m_rho);
  static double xs_diff_pi_pi0_rho(const double s, const double t,
                                   const double m_rho);
  static double xs_diff_pi0_rho0_pi0(const double s, const double t,
                                     const double m_rho);
  static double xs_diff_pi_rho0_pi(const double s, const double t,
                                   const double m_rho);
 
  static double xs_diff_pi_rho_pi0_rho_mediated(const double s, const double t,
                                                const double m_rho);
  static double xs_diff_pi_rho_pi0_omega_mediated(const double s,
                                                  const double t,
                                                  const double m_rho);
 
  static double xs_diff_pi0_rho_pi_rho_mediated(const double s, const double t,
                                                const double m_rho);
  static double xs_diff_pi0_rho_pi_omega_mediated(const double s,
                                                  const double t,
                                                  const double m_rho);
 
  // Todo{JR}: Needed for tabulation.
  static double s_min, s_max, t_min, t_max;
 
 private:
 
  constexpr static double Const = 0.059;
  constexpr static double g_POR = 22.6;
  constexpr static double ghat = 6.4483;
  constexpr static double eta1 = 2.3920;
  constexpr static double eta2 = 1.9430;
  constexpr static double delta = -0.6426;
  constexpr static double C4 = -0.14095;
  constexpr static double Gammaa1 = 0.4;
  constexpr static double Pi = M_PI;
  constexpr static double a1_mass = 1.26;
};
 
}  // namespace smash
 
#endif  // SRC_INCLUDE_SMASH_CROSSSECTIONSPHOTON_H_