bremsstrahlungaction.h

Go to the documentation of this file.

/*
 *
 *    Copyright (c) 2020 -
 *      SMASH Team
 *
 *    GNU General Public License (GPLv3 or later)
 *
 */
 
#ifndef SRC_INCLUDE_BREMSSTRAHLUNG_H_
#define SRC_INCLUDE_BREMSSTRAHLUNG_H_
 
#include "scatteraction.h"
 
namespace smash {
class BremsstrahlungAction : public ScatterAction {
 public:
  BremsstrahlungAction(const ParticleList &in, const double time,
                       const int n_frac_photons,
                       const double hadronic_cross_section_input);
  void perform_bremsstrahlung(const OutputsList &outputs);
 
  void generate_final_state() override;
 
  void sample_3body_phasespace() override;
 
  double get_total_weight() const override { return weight_; }
 
  double hadronic_cross_section() const { return hadronic_cross_section_; }
 
  void add_dummy_hadronic_process(double reaction_cross_section);
 
  void add_single_process() {
    add_processes<CollisionBranch>(brems_cross_sections(),
                                   collision_processes_bremsstrahlung_,
                                   cross_section_bremsstrahlung_);
  }
 
  enum class ReactionType {
    no_reaction,
    pi_z_pi_m,
    pi_z_pi_p,
    pi_p_pi_m,
    pi_m_pi_m,
    pi_p_pi_p,
    pi_z_pi_z
  };
 
  static ReactionType bremsstrahlung_reaction_type(const ParticleList &in);
 
  static bool is_bremsstrahlung_reaction(const ParticleList &in) {
    return bremsstrahlung_reaction_type(in) != ReactionType::no_reaction;
  }
 
 private:
  CollisionBranchList collision_processes_bremsstrahlung_;
 
  const ReactionType reac_;
 
  const int number_of_fractional_photons_;
 
  double weight_ = 0.0;
 
  double cross_section_bremsstrahlung_ = 0.0;
 
  const double hadronic_cross_section_;
 
  double k_;
 
  double theta_;
 
  void create_interpolations();
 
  CollisionBranchList brems_cross_sections();
 
  std::pair<double, double> brems_diff_cross_sections();
};
 
}  // namespace smash
 
#endif  // SRC_INCLUDE_BREMSSTRAHLUNG_H_