scatteractionmulti.h

Go to the documentation of this file.

/*
 *
 *    Copyright (c) 2020-2022
 *      SMASH Team
 *
 *    GNU General Public License (GPLv3 or later)
 *
 */
 
#ifndef SRC_INCLUDE_SMASH_SCATTERACTIONMULTI_H_
#define SRC_INCLUDE_SMASH_SCATTERACTIONMULTI_H_
 
#include "action.h"
 
namespace smash {
 
class ScatterActionMulti : public Action {
 public:
  ScatterActionMulti(const ParticleList& in_plist, double time);
 
  void generate_final_state() override;
 
  double get_total_weight() const override;
 
  double get_partial_weight() const override;
 
  void add_possible_reactions(double dt, const double gcell_vol,
                              const MultiParticleReactionsBitSet incl_multi);
 
  const CollisionBranchList& reaction_channels() { return reaction_channels_; }
 
  class InvalidScatterActionMulti : public std::invalid_argument {
    using std::invalid_argument::invalid_argument;
  };
 
  double calculate_I3(const double sqrts) const;
 
  double parametrizaton_phi4(const double man_s) const;
 
 protected:
  /*
   * \ingroup logging
   * Writes information about this action to the \p out stream.
   */
  void format_debug_output(std::ostream& out) const override;
 
 private:
  void add_reaction(CollisionBranchPtr p);
 
  void add_reactions(CollisionBranchList pv);
 
  void annihilation();
 
  void n_to_two();
 
  double probability_three_to_one(const ParticleType& type_out, double dt,
                                  const double gcell_vol,
                                  const int degen_sym_factor = 1) const;
 
  double probability_three_to_two(const ParticleType& type_out1,
                                  const ParticleType& type_out2, double dt,
                                  const double gcell_vol,
                                  const double degen_sym_factor = 1.0) const;
  double probability_four_to_two(const ParticleType& type_out1,
                                 const ParticleType& type_out2, double dt,
                                 const double gcell_vol,
                                 const double degen_sym_factor = 1.0) const;
  double probability_five_to_two(const double m_out, double dt,
                                 const double gcell_vol,
                                 const double degen_sym_factor = 1.0) const;
 
  double parametrizaton_phi5_pions(const double man_s) const;
 
  double react_degen_factor(const int spin_factor_inc,
                            const int spin_degen_out1,
                            const int spin_degen_out2) const {
    return static_cast<double>(spin_degen_out1 * spin_degen_out2) /
           static_cast<double>(spin_factor_inc);
  }
 
  bool three_different_pions(const ParticleData& data_a,
                             const ParticleData& data_b,
                             const ParticleData& data_c) const;
 
  bool two_pions_eta(const ParticleData& data_a, const ParticleData& data_b,
                     const ParticleData& data_c) const;
 
  bool all_incoming_particles_are_pions_have_zero_charge_only_one_piz() const;
 
  double total_probability_;
 
  double partial_probability_;
 
  CollisionBranchList reaction_channels_;
};
 
}  // namespace smash
 
#endif  // SRC_INCLUDE_SMASH_SCATTERACTIONMULTI_H_