15 const ParticleList& plist,
double t_max)
const {
16 std::vector<ActionPtr> actions;
20 double time_until_crossing = t_max;
22 for (
int i = 0; i < 3; i++) {
23 double t = t_max + 1.;
25 t = (
l_[i] - r[i]) / v[i];
29 if (t < time_until_crossing) {
30 time_until_crossing = t;
38 FourVector crossing_point(
p.position().x0() + time_until_crossing,
39 r + v * time_until_crossing);
40 crossing_point[i_cross + 1] = ((v[i_cross] > 0.0) ? 0.0 :
l_[i_cross]);
44 ActionPtr action = make_unique<WallcrossingAction>(
p, outgoing_particle,
46 actions.emplace_back(std::move(action));
const std::array< double, 3 > l_
Periods in x,y,z directions in fm.
The ThreeVector class represents a physical three-vector with the components .
constexpr double really_small
Numerical error tolerance.
ActionList find_actions_in_cell(const ParticleList &plist, double t_max) const override
Find the next wall crossings for every particle before time t_max.
void set_4position(const FourVector &pos)
Set the particle's 4-position directly.
The FourVector class holds relevant values in Minkowski spacetime with (+, −, −, −) metric signature.
ParticleData contains the dynamic information of a certain particle.