parametrizations.h

Go to the documentation of this file.

/*
 *
 *    Copyright (c) 2013-2018
 *      SMASH Team
 *
 *    GNU General Public License (GPLv3 or later)
 *
 */
#ifndef SRC_INCLUDE_PARAMETRIZATIONS_H_
#define SRC_INCLUDE_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 {

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 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_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_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_PARAMETRIZATIONS_H_