hypersurfacecrossingaction.h

Go to the documentation of this file.

/*
 *
 *    Copyright (c) 2019-2020,2022
 *      SMASH Team
 *
 *    GNU General Public License (GPLv3 or later)
 *
 */
 
#ifndef SRC_INCLUDE_SMASH_HYPERSURFACECROSSINGACTION_H_
#define SRC_INCLUDE_SMASH_HYPERSURFACECROSSINGACTION_H_
 
#include <vector>
 
#include "action.h"
#include "actionfinderfactory.h"
 
namespace smash {
 
class HypersurfacecrossingAction : public Action {
 public:
  HypersurfacecrossingAction(const ParticleData &in_part,
                             const ParticleData &out_part,
                             const double time_until)
      : Action(in_part, out_part, time_until,
               ProcessType::HyperSurfaceCrossing) {}
  double get_total_weight() const override { return 0.0; };
  double get_partial_weight() const override { return 0.0; };
  void format_debug_output(std::ostream &out) const override {
    out << "Hypersurface crossing of " << incoming_particles_;
  }
 
  void generate_final_state() override;
 
  double check_conservation(const uint32_t id_process) const override;
};
 
class HyperSurfaceCrossActionsFinder : public ActionFinderInterface {
 public:
  explicit HyperSurfaceCrossActionsFinder(double tau, double y, double pT)
      : prop_time_{tau}, rap_cut_{y}, pT_cut_{pT} {};
 
  ActionList find_actions_in_cell(
      const ParticleList &plist, double dt, const double,
      const std::vector<FourVector> &beam_momentum) const override;
 
  ActionList find_actions_with_neighbors(
      const ParticleList &, const ParticleList &, double,
      const std::vector<FourVector> &) const override {
    return {};
  }
 
  ActionList find_actions_with_surrounding_particles(
      const ParticleList &, const Particles &, double,
      const std::vector<FourVector> &) const override {
    return {};
  }
 
  ActionList find_final_actions(const Particles &, bool) const override {
    return {};
  }
 
 private:
  const double prop_time_;
 
  const double rap_cut_;
 
  const double pT_cut_;
 
  bool crosses_hypersurface(ParticleData &pdata_before_propagation,
                            ParticleData &pdata_after_propagation,
                            const double tau) const;
 
  FourVector coordinates_on_hypersurface(ParticleData &pdata_before_propagation,
                                         ParticleData &pdata_after_propagation,
                                         const double tau) const;
};
 
}  // namespace smash
 
#endif  // SRC_INCLUDE_SMASH_HYPERSURFACECROSSINGACTION_H_