hypersurfacecrossingfinder.h

Go to the documentation of this file.

/*
 *
 *    Copyright (c) 2019-2020,2022-2025
 *      SMASH Team
 *
 *    GNU General Public License (GPLv3 or later)
 *
 */
 
#ifndef SRC_INCLUDE_SMASH_HYPERSURFACECROSSINGFINDER_H_
#define SRC_INCLUDE_SMASH_HYPERSURFACECROSSINGFINDER_H_
 
#include <vector>
 
#include "action.h"
#include "actionfinderfactory.h"
 
namespace smash {
 
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 &) const override { return {}; }
 
  static void warn_if_some_particles_did_not_cross(
      const size_t number_of_particles, bool impose_kinematic_cut);
 
 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_HYPERSURFACECROSSINGFINDER_H_