parametrizations.h

Go to the documentation of this file.

/*
 *
 *    Copyright (c) 2013-2018,2020,2023-2024
 *      SMASH Team
 *
 *    GNU General Public License (GPLv3 or later)
 *
 */
#ifndef SRC_INCLUDE_SMASH_PARAMETRIZATIONS_H_
#define SRC_INCLUDE_SMASH_PARAMETRIZATIONS_H_
 
#include <unordered_map>
#include <utility>
 
#include "particletype.h"
 
/* All quantities in this file use they same units as the rest of SMASH.
 * That is: GeV for energies and momenta, fm for distances and time, and mb for
 * cross-sections. */
 
namespace smash {
 
bool parametrization_exists(const PdgCode& pdg_a, const PdgCode& pdg_b);
double xs_high_energy(double mandelstam_s, bool is_opposite_charge, double ma,
                      double mb, double P, double R1, double R2);
 
double pp_high_energy(double mandelstam_s);
 
double ppbar_high_energy(double mandelstam_s);
 
double np_high_energy(double mandelstam_s);
 
double npbar_high_energy(double mandelstam_s);
 
double piplusp_high_energy(double mandelstam_s);
 
double piminusp_high_energy(double mandelstam_s);
 
double xs_ppbar_annihilation(double mandelstam_s);
 
double xs_string_hard(double mandelstam_s, double xs_0, double e_0,
                      double lambda_pow);
 
double NN_string_hard(double mandelstam_s);
 
double Npi_string_hard(double mandelstam_s);
 
double pipi_string_hard(double mandelstam_s);
 
double pipluspiminus_total(double sqrts);
 
double pizeropizero_total(double sqrts);
 
double piplusp_total(double sqrts);
 
double piplusp_elastic_high_energy(double mandelstam_s, double m1, double m2);
 
double piplusp_elastic_AQM(double mandelstam_s, double m1, double m2);
 
double piplusp_elastic(double mandelstam_s);
 
double piplusp_sigmapluskplus_pdg(double mandelstam_s);
 
double piminusp_total(double sqrts);
 
double piminusp_elastic(double mandelstam_s);
 
double piminusp_lambdak0_pdg(double mandelstam_s);
 
double piminusp_sigmaminuskplus_pdg(double mandelstam_s);
 
double piminusp_sigma0k0_res(double mandelstam_s);
 
double pp_elastic(double mandelstam_s);
 
double pp_elastic_high_energy(double mandelstam_s, double m1, double m2);
 
double pp_total(double mandelstam_s);
 
double np_elastic(double mandelstam_s);
 
double np_total(double mandelstam_s);
 
double ppbar_elastic(double mandelstam_s);
 
double ppbar_total(double mandelstam_s);
 
double deuteron_pion_elastic(double mandelstam_s);
 
double deuteron_nucleon_elastic(double mandelstam_s);
 
double kplusp_total(double mandelstam_s);
 
double kplusn_total(double mandelstam_s);
 
double kminusn_total(double mandelstam_s);
 
double kminusp_total(double mandelstam_s);
 
double kplusp_elastic_background(double mandelstam_s);
 
double kplusn_elastic_background(double mandelstam_s);
 
double kplusn_k0p(double mandelstam_s);
 
double kminusp_elastic_background(double mandelstam_s);
 
double kminusn_elastic_background(double mandelstam_s);
 
double k0p_elastic_background(double mandelstam_s);
 
double k0n_elastic_background(double mandelstam_s);
 
double kbar0p_elastic_background(double mandelstam_s);
 
double kbar0n_elastic_background(double mandelstam_s);
 
double kplusp_inelastic_background(double mandelstam_s);
 
double kplusn_inelastic_background(double mandelstam_s);
 
struct pair_hash {
  std::size_t operator()(const std::pair<uint64_t, uint64_t>& p) const {
    auto h1 = std::hash<uint64_t>{}(p.first);
    auto h2 = std::hash<uint64_t>{}(p.second);
 
    /* In our case the integers are PDG codes. We know they are different
     * and their order is defined, so we can simply combine the hashes
     * using XOR. Note that this yields 0 for h1 == h2. Also,
     * std::swap(h1, h2) does not not change the final hash. */
    assert(h1 != h2);
    return h1 ^ h2;
  }
};
 
class KaonNucleonRatios {
 private:
  mutable std::unordered_map<std::pair<uint64_t, uint64_t>, double, pair_hash>
      ratios_;
 
 public:
  KaonNucleonRatios() : ratios_({}) {}
 
  double get_ratio(const ParticleType& a, const ParticleType& b,
                   const ParticleType& c, const ParticleType& d) const;
};
 
extern /*thread_local (see #3075)*/ KaonNucleonRatios kaon_nucleon_ratios;
 
double kminusp_kbar0n(double mandelstam_s);
 
double kminusp_piminussigmaplus(double sqrts);
 
double kminusp_piplussigmaminus(double sqrts);
 
double kminusp_pi0sigma0(double sqrts);
 
double kminusp_pi0lambda(double sqrts);
 
double kminusn_piminussigma0(double sqrts);
 
double kminusn_pi0sigmaminus(double sqrts);
 
double kminusn_piminuslambda(double sqrts);
 
double lambdalambda_ximinusp(double sqrts_sqrts0, double p_N, double p_lambda);
 
double lambdalambda_xi0n(double sqrts_sqrts0, double p_N, double p_lambda);
 
double lambdasigmaplus_xi0p(double sqrts_sqrts0);
 
double lambdasigmaminus_ximinusn(double sqrts_sqrts0);
 
double lambdasigma0_ximinusp(double sqrts_sqrts0);
 
double lambdasigma0_xi0n(double sqrts_sqrts0);
 
double sigma0sigma0_ximinusp(double sqrts_sqrts0);
 
double sigma0sigma0_xi0n(double sqrts_sqrts0);
 
double sigmaplussigmaminus_xi0p(double sqrts_sqrts0);
 
double sigma0sigmaminus_ximinusn(double sqrts_sqrts0);
 
double sigmaplussigmaminus_ximinusp(double sqrts_sqrts0);
 
double sigmaplussigmaminus_xi0n(double sqrts_sqrts0);
 
}  // namespace smash
 
#endif  // SRC_INCLUDE_SMASH_PARAMETRIZATIONS_H_