freeforallaction.h

Go to the documentation of this file.

/*
 *
 *    Copyright (c) 2022-2023
 *      SMASH Team
 *
 *    GNU General Public License (GPLv3 or later)
 *
 */
 
#ifndef SRC_INCLUDE_SMASH_FREEFORALLACTION_H_
#define SRC_INCLUDE_SMASH_FREEFORALLACTION_H_
 
#include "action.h"
 
namespace smash {
 
class FreeforallAction : public Action {
 public:
  FreeforallAction(const ParticleList &in_part, const ParticleList &out_part,
                   double absolute_labframe_time)
      : Action(in_part, out_part, absolute_labframe_time,
               ProcessType::Freeforall) {
    if (!in_part.empty() && !out_part.empty()) {
      throw std::invalid_argument(
          "The FreeforallAction was used with two non-empty lists. This is not "
          "allowed. Please use the action twice, if particles should be added "
          "and removed.");
    }
  }
 
  double get_total_weight() const override { return 0.0; }
  double get_partial_weight() const override { return 0.0; }
 
  void generate_final_state() override {
    // Set time for arbitrary outgoing particles to time of action
    for (auto &particle : outgoing_particles_) {
      const double t = particle.position().x0();
      const FourVector u(1.0, particle.velocity());
      particle.set_formation_time(t);
      particle.set_4position(particle.position() +
                             u * (time_of_execution() - t));
    }
  }
 
  void format_debug_output(std::ostream &out) const override {
    out << "Freeforall action of " << incoming_particles_.size() << " to "
        << outgoing_particles_.size() << " particles.";
  }
};
 
}  // namespace smash
 
#endif  // SRC_INCLUDE_SMASH_FREEFORALLACTION_H_