hypersurfacecrossingaction.h

Go to the documentation of this file.

/*
 *
 *    Copyright (c) 2019-2019
 *      SMASH Team
 *
 *    GNU General Public License (GPLv3 or later)
 *
 */

#ifndef SRC_INCLUDE_HYPERSURFACECROSSINGACTION_H_
#define SRC_INCLUDE_HYPERSURFACECROSSINGACTION_H_

#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;

  void check_conservation(const uint32_t id_process) const override;
};

class HyperSurfaceCrossActionsFinder : public ActionFinderInterface {
 public:
  explicit HyperSurfaceCrossActionsFinder(double tau) : prop_time_{tau} {};

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

  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_HYPERSURFACECROSSINGACTION_H_