smash::Action Class Reference

#include <action.h>

Action is the base class for a generic process that takes a number of incoming particles and transforms them into any number of outgoing particles.

Currently such an action can be either a decay, a two-body collision, a wallcrossing or a thermalization. (see derived classes).

Definition at line 34 of file action.h.

Inheritance diagram for smash::Action:

Collaboration diagram for smash::Action:

Classes
class	InvalidResonanceFormation
	Thrown for example when ScatterAction is called to perform with a wrong number of final-state particles or when the energy is too low to produce the resonance. More...

Public Member Functions
	Action (const ParticleList &in_part, double time)
	Construct an action object with incoming particles and relative time. More...
	Action (const ParticleData &in_part, const ParticleData &out_part, double time, ProcessType type)
	Construct an action object with the incoming particles, relative time, and the already known outgoing particles and type of the process. More...
	Action (const ParticleList &in_part, const ParticleList &out_part, double absolute_execution_time, ProcessType type)
	Construct an action object with the incoming particles, absolute time, and the already known outgoing particles and type of the process. More...
	Action (const Action &)=delete
	Copying is disabled. Use pointers or create a new Action. More...
virtual	~Action ()
	Virtual Destructor. More...
bool	operator< (const Action &rhs) const
	Determine whether one action takes place before another in time. More...
virtual double	get_total_weight () const =0
	Return the total weight value, which is mainly used for the weight output entry. More...
virtual double	get_partial_weight () const =0
	Return the specific weight for the chosen outgoing channel, which is mainly used for the partial weight output entry. More...
virtual ProcessType	get_type () const
	Get the process type. More...
template<typename Branch >
void	add_process (ProcessBranchPtr< Branch > &p, ProcessBranchList< Branch > &subprocesses, double &total_weight)
	Add a new subprocess. More...
template<typename Branch >
void	add_processes (ProcessBranchList< Branch > pv, ProcessBranchList< Branch > &subprocesses, double &total_weight)
	Add several new subprocesses at once. More...
virtual void	generate_final_state ()=0
	Generate the final state for this action. More...
virtual void	perform (Particles *particles, uint32_t id_process)
	Actually perform the action, e.g. More...
bool	is_valid (const Particles &particles) const
	Check whether the action still applies. More...
bool	is_pauli_blocked (const Particles &particles, const PauliBlocker &p_bl) const
	Check if the action is Pauli-blocked. More...
const ParticleList &	incoming_particles () const
	Get the list of particles that go into the action. More...
void	update_incoming (const Particles &particles)
	Update the incoming particles that are stored in this action to the state they have in the global particle list. More...
const ParticleList &	outgoing_particles () const
	Get the list of particles that resulted from the action. More...
double	time_of_execution () const
	Get the time at which the action is supposed to be performed. More...
void	check_conservation (const uint32_t id_process) const
	Check various conservation laws. More...
double	sqrt_s () const
	Determine the total energy in the center-of-mass frame [GeV]. More...
FourVector	total_momentum_of_outgoing_particles () const
	Calculate the total kinetic momentum of the outgoing particles. More...
FourVector	get_interaction_point () const
	Get the interaction point. More...
std::pair< FourVector, FourVector >	get_potential_at_interaction_point () const
	Get the skyrme and asymmetry potential at the interaction point. More...

Protected Member Functions
FourVector	total_momentum () const
	Sum of 4-momenta of incoming particles. More...
template<typename Branch >
const Branch *	choose_channel (const ProcessBranchList< Branch > &subprocesses, double total_weight)
	Decide for a particular final-state channel via Monte-Carlo and return it as a ProcessBranch. More...
virtual std::pair< double, double >	sample_masses (double kinetic_energy_cm) const
	Sample final-state masses in general X->2 processes (thus also fixing the absolute c.o.m. More...
virtual void	sample_angles (std::pair< double, double > masses, double kinetic_energy_cm)
	Sample final-state momenta in general X->2 processes (here: using an isotropical angular distribution). More...
void	sample_2body_phasespace ()
	Sample the full 2-body phase-space (masses, momenta, angles) in the center-of-mass frame for the final state particles. More...
virtual void	sample_3body_phasespace ()
	Sample the full 3-body phase-space (masses, momenta, angles) in the center-of-mass frame for the final state particles. More...
virtual void	format_debug_output (std::ostream &out) const =0
	Writes information about this action to the out stream. More...

Protected Attributes
ParticleList	incoming_particles_
	List with data of incoming particles. More...
ParticleList	outgoing_particles_
	Initially this stores only the PDG codes of final-state particles. More...
const double	time_of_execution_
	Time at which the action is supposed to be performed (absolute time in the lab frame in fm/c). More...
ProcessType	process_type_
	type of process More...

Private Member Functions
const ParticleType &	type_of_pout (const ParticleData &p_out) const
	Get the type of a given particle. More...
const ParticleType &	type_of_pout (const ParticleTypePtr &p_out) const
	Get the particle type for given pointer to a particle type. More...

Friends
std::ostream &	operator<< (std::ostream &out, const Action &action)
	Dispatches formatting to the virtual Action::format_debug_output function. More...

Constructor & Destructor Documentation

smash::Action::Action ( const ParticleList &  in_part, double  time  )

inline

Construct an action object with incoming particles and relative time.

Parameters

[in]	in_part	list of incoming particles
[in]	time	time at which the action is supposed to take place (relative to the current time of the incoming particles)

Definition at line 43 of file action.h.

      : incoming_particles_(in_part),
        time_of_execution_(time + in_part[0].position().x0()) {}

Here is the caller graph for this function:

smash::Action::Action ( const ParticleData &  in_part, const ParticleData &  out_part, double  time, ProcessType  type  )

inline

Construct an action object with the incoming particles, relative time, and the already known outgoing particles and type of the process.

Parameters

[in]	in_part	list of incoming particles
[in]	out_part	list of outgoing particles
[in]	time	time at which the action is supposed to take place (relative to the current time of the incoming particles)
[in]	type	type of the interaction

Definition at line 57 of file action.h.

      : incoming_particles_({in_part}),
        outgoing_particles_({out_part}),
        time_of_execution_(time + in_part.position().x0()),
        process_type_(type) {}

smash::Action::Action ( const ParticleList &  in_part, const ParticleList &  out_part, double  absolute_execution_time, ProcessType  type  )

inline

Construct an action object with the incoming particles, absolute time, and the already known outgoing particles and type of the process.

Parameters

[in]	in_part	list of incoming particles
[in]	out_part	list of outgoing particles
[in]	absolute_execution_time	absolute time at which the action is supposed to take place
[in]	type	type of the interaction

Definition at line 74 of file action.h.

      : incoming_particles_(std::move(in_part)),
        outgoing_particles_(std::move(out_part)),
        time_of_execution_(absolute_execution_time),
        process_type_(type) {}

Here is the call graph for this function:

smash::Action::Action ( const Action &  )

delete

Copying is disabled. Use pointers or create a new Action.

smash::Action::~Action ( )

virtualdefault

Virtual Destructor.

Destructor.

The declaration of the destructor is necessary to make it virtual.

Here is the caller graph for this function:

Member Function Documentation

bool smash::Action::operator< ( const Action &  rhs ) const

inline

Determine whether one action takes place before another in time.

Returns

if the first argument action takes place before the other

Definition at line 95 of file action.h.

                                          {
    return time_of_execution_ < rhs.time_of_execution_;
  }

Here is the call graph for this function:

virtual double smash::Action::get_total_weight ( ) const

pure virtual

Return the total weight value, which is mainly used for the weight output entry.

It has different meanings depending of the type of action. It is the total cross section in case of a ScatterAction, the total decay width in case of a DecayAction and the shining weight in case of a DecayActionDilepton.

Prefer to use a more specific function. If there is no weight for the action type, 0 should be returned.

Returns

total cross section, decay width or shining weight

Implemented in smash::ScatterAction, smash::DecayAction, smash::ScatterActionPhoton, smash::DecayActionDilepton, smash::ThermalizationAction, and smash::WallcrossingAction.

Here is the caller graph for this function:

virtual double smash::Action::get_partial_weight ( ) const

pure virtual

Return the specific weight for the chosen outgoing channel, which is mainly used for the partial weight output entry.

For scatterings it will be the partial cross section, for decays (including dilepton decays) the partial decay width.

If there is no weight for the action type, 0 should be returned.

Returns

specific weight for the chosen output channel.

Implemented in smash::ScatterAction, smash::DecayAction, smash::ThermalizationAction, and smash::WallcrossingAction.

Here is the caller graph for this function:

virtual ProcessType smash::Action::get_type ( ) const

inlinevirtual

Get the process type.

Returns

type of the process

Definition at line 130 of file action.h.

{ return process_type_; }

Here is the caller graph for this function:

template<typename Branch >

void smash::Action::add_process ( ProcessBranchPtr< Branch > &  p, ProcessBranchList< Branch > &  subprocesses, double &  total_weight  )

inline

Add a new subprocess.

Parameters

[in]	p	process to be added
[out]	subprocesses	processes, where p is added to
[out]	total_weight	summed weights of all the subprocesses

Definition at line 140 of file action.h.

                                         {
    if (p->weight() > 0) {
      total_weight += p->weight();
      subprocesses.emplace_back(std::move(p));
    }
  }

template<typename Branch >

void smash::Action::add_processes ( ProcessBranchList< Branch >  pv, ProcessBranchList< Branch > &  subprocesses, double &  total_weight  )

inline

Add several new subprocesses at once.

Parameters

[in]	pv	processes list to be added
[out]	subprocesses	processes, where pv are added to
[out]	total_weight	summed weights of all the subprocesses

Definition at line 157 of file action.h.

                                           {
    subprocesses.reserve(subprocesses.size() + pv.size());
    for (auto &proc : pv) {
      if (proc->weight() > 0) {
        total_weight += proc->weight();
        subprocesses.emplace_back(std::move(proc));
      }
    }
  }

Here is the call graph for this function:

virtual void smash::Action::generate_final_state ( )

pure virtual

Generate the final state for this action.

This function selects a subprocess by Monte-Carlo decision and sets up the final-state particles in phase space.

Implemented in smash::ScatterAction, smash::ScatterActionPhoton, smash::DecayAction, smash::WallcrossingAction, and smash::ThermalizationAction.

Here is the caller graph for this function:

void smash::Action::perform ( Particles *  particles, uint32_t  id_process  )

virtual

Actually perform the action, e.g.

carry out a decay or scattering by updating the particle list.

This function removes the initial-state particles from the particle list and then inserts the final-state particles. It does not do any sanity checks, but assumes that is_valid has been called to determine if the action is still valid.

Parameters

[in]	id_process	unique id of the performed process
[out]	particles	particle list that is updated

Note that you are required to increase id_process before the next call, such that you get unique numbers.

Definition at line 94 of file action.cc.

                                                              {
  assert(id_process != 0);
  const auto &log = logger<LogArea::Action>();

  for (ParticleData &p : outgoing_particles_) {
    // store the history info
    if (process_type_ != ProcessType::Wall) {
      p.set_history(p.get_history().collisions_per_particle + 1, id_process,
                    process_type_, time_of_execution_, incoming_particles_);
    }
  }

  /* For elastic collisions and box wall crossings it is not necessary to remove
   * particles from the list and insert new ones, it is enough to update their
   * properties. */
  particles->update(incoming_particles_, outgoing_particles_,
                    (process_type_ != ProcessType::Elastic) &&
                        (process_type_ != ProcessType::Wall));

  log.debug("Particle map now has ", particles->size(), " elements.");

  /* Check the conservation laws if the modifications of the total kinetic
   * energy of the outgoing particles by the mean field potentials are not
   * taken into account. */
  if (UB_lat_pointer == nullptr && UI3_lat_pointer == nullptr) {
    check_conservation(id_process);
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

bool smash::Action::is_valid

(

const Particles &

particles

)

const

Check whether the action still applies.

It can happen that a different action removed the incoming_particles from the set of existing particles in the experiment, or that the particle has scattered elastically in the meantime. In this case the Action doesn't apply anymore and should be discarded.

Parameters

[in]

particles

current particle list

Returns

true, if action still applies; false otherwise

Definition at line 29 of file action.cc.

                                                      {
  return std::all_of(
      incoming_particles_.begin(), incoming_particles_.end(),
      [&particles](const ParticleData &p) { return particles.is_valid(p); });
}

Here is the call graph for this function:

Here is the caller graph for this function:

bool smash::Action::is_pauli_blocked	(	const Particles &	particles,
		const PauliBlocker &	p_bl
	)		const

Check if the action is Pauli-blocked.

If there are baryons in the final state then blocking probability is \( 1 - \Pi (1-f_i) \), where the product is taken by all fermions in the final state and \( f_i \) denotes the phase-space density at the position of i-th final-state fermion.

Parameters

[in]	particles	current particle list
[in]	p_bl	PauliBlocker that stores the configurations concerning Pauli-blocking.

Returns

true, if the action is Pauli-blocked, false otherwise

Definition at line 35 of file action.cc.

                                                              {
  // Wall-crossing actions should never be blocked: currently
  // if the action is blocked, a particle continues to propagate in a straight
  // line. This would simply bring it out of the box.
  if (process_type_ == ProcessType::Wall) {
    return false;
  }
  const auto &log = logger<LogArea::PauliBlocking>();
  for (const auto &p : outgoing_particles_) {
    if (p.is_baryon()) {
      const auto f =
          p_bl.phasespace_dens(p.position().threevec(), p.momentum().threevec(),
                               particles, p.pdgcode(), incoming_particles_);
      if (f > random::uniform(0., 1.)) {
        log.debug("Action ", *this, " is pauli-blocked with f = ", f);
        return true;
      }
    }
  }
  return false;
}

Here is the call graph for this function:

Here is the caller graph for this function:

const ParticleList & smash::Action::incoming_particles

(

)

const

Get the list of particles that go into the action.

Returns

a list of incoming particles

Definition at line 58 of file action.cc.

                                                     {
  return incoming_particles_;
}

Here is the caller graph for this function:

void smash::Action::update_incoming

(

const Particles &

particles

)

Update the incoming particles that are stored in this action to the state they have in the global particle list.

Parameters

[in]

particles

current particle list

Definition at line 62 of file action.cc.

                                                       {
  for (auto &p : incoming_particles_) {
    p = particles.lookup(p);
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

const ParticleList& smash::Action::outgoing_particles ( ) const

inline

Get the list of particles that resulted from the action.

Returns

list of outgoing particles

Definition at line 244 of file action.h.

{ return outgoing_particles_; }

Here is the caller graph for this function:

double smash::Action::time_of_execution ( ) const

inline

Get the time at which the action is supposed to be performed.

Returns

absolute time in the calculation frame in fm/c

Definition at line 251 of file action.h.

{ return time_of_execution_; }

Here is the call graph for this function:

void smash::Action::check_conservation

(

const uint32_t

id_process

)

const

Check various conservation laws.

Parameters

[in]

id_process

process id only used for debugging output

Definition at line 247 of file action.cc.

                                                               {
  QuantumNumbers before(incoming_particles_);
  QuantumNumbers after(outgoing_particles_);
  if (before != after) {
    std::stringstream particle_names;
    for (const auto &p : incoming_particles_) {
      particle_names << p.type().name();
    }
    particle_names << " vs. ";
    for (const auto &p : outgoing_particles_) {
      particle_names << p.type().name();
    }
    particle_names << "\n";
    const auto &log = logger<LogArea::Action>();
    std::string err_msg = before.report_deviations(after);
    log.error() << particle_names.str() << err_msg;
    /* Pythia does not conserve energy and momentum at high energy, so we just
     * print the error and continue. */
    if ((is_string_soft_process(process_type_)) ||
        (process_type_ == ProcessType::StringHard)) {
      return;
    }
    if (id_process == ID_PROCESS_PHOTON) {
      throw std::runtime_error("Conservation laws violated in photon process");
    } else {
      throw std::runtime_error("Conservation laws violated in process " +
                               std::to_string(id_process));
    }
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

double smash::Action::sqrt_s ( ) const

inline

Determine the total energy in the center-of-mass frame [GeV].

Returns

\( \sqrt{s}\) of incoming particles

Definition at line 265 of file action.h.

{ return total_momentum().abs(); }

Here is the call graph for this function:

Here is the caller graph for this function:

FourVector smash::Action::total_momentum_of_outgoing_particles

(

)

const

Calculate the total kinetic momentum of the outgoing particles.

This function is used when the species of the outgoing particles are already determined.

Returns

total kinetic momentum of the outgoing particles [GeV]

Definition at line 123 of file action.cc.

                                                              {
  const auto potentials = get_potential_at_interaction_point();
  /* scale_B returns the difference of the total force scales of the skyrme
   * potential between the initial and final states. */
  double scale_B = 0.0;
  /* scale_I3 returns the difference of the total force scales of the symmetry
   * potential between the initial and final states. */
  double scale_I3 = 0.0;
  for (const auto &p_in : incoming_particles_) {
    // Get the force scale of the incoming particle.
    const auto scale =
        ((pot_pointer != nullptr) ? pot_pointer->force_scale(p_in.type())
                                  : std::make_pair(0.0, 0));
    scale_B += scale.first;
    scale_I3 += scale.second * p_in.type().isospin3_rel();
  }
  for (const auto &p_out : outgoing_particles_) {
    // Get the force scale of the outgoing particle.
    const auto scale = ((pot_pointer != nullptr)
                            ? pot_pointer->force_scale(type_of_pout(p_out))
                            : std::make_pair(0.0, 0));
    scale_B -= scale.first;
    scale_I3 -= scale.second * type_of_pout(p_out).isospin3_rel();
  }
  /* Rescale to get the potential difference between the
   * initial and final state, and thus get the total momentum
   * of the outgoing particles*/
  return total_momentum() + potentials.first * scale_B +
         potentials.second * scale_I3;
}

Here is the call graph for this function:

Here is the caller graph for this function:

FourVector smash::Action::get_interaction_point

(

)

const

Get the interaction point.

Returns

four vector of interaction point

Definition at line 68 of file action.cc.

                                               {
  // Estimate for the interaction point in the calculational frame
  FourVector interaction_point = FourVector(0., 0., 0., 0.);
  for (const auto &part : incoming_particles_) {
    interaction_point += part.position();
  }
  interaction_point /= incoming_particles_.size();
  return interaction_point;
}

Here is the caller graph for this function:

std::pair< FourVector, FourVector > smash::Action::get_potential_at_interaction_point

(

)

const

Get the skyrme and asymmetry potential at the interaction point.

Returns

skyrme and asymmetry potential [GeV]

Definition at line 78 of file action.cc.

          {
  const ThreeVector r = get_interaction_point().threevec();
  FourVector UB = FourVector();
  FourVector UI3 = FourVector();
  /* Check:
   * Lattice is turned on. */
  if (UB_lat_pointer != nullptr) {
    UB_lat_pointer->value_at(r, UB);
  }
  if (UI3_lat_pointer != nullptr) {
    UI3_lat_pointer->value_at(r, UI3);
  }
  return std::make_pair(UB, UI3);
}

Here is the call graph for this function:

Here is the caller graph for this function:

FourVector smash::Action::total_momentum ( ) const

inlineprotected

Sum of 4-momenta of incoming particles.

Definition at line 323 of file action.h.

                                    {
    FourVector mom(0.0, 0.0, 0.0, 0.0);
    for (const auto &p : incoming_particles_) {
      mom += p.momentum();
    }
    return mom;
  }

Here is the caller graph for this function:

template<typename Branch >

const Branch* smash::Action::choose_channel ( const ProcessBranchList< Branch > &  subprocesses, double  total_weight  )

inlineprotected

Decide for a particular final-state channel via Monte-Carlo and return it as a ProcessBranch.

Template Parameters

Branch

Type of processbranch

Parameters

[in]	subprocesses	list of possible processes
[in]	total_weight	summed weight of all processes

Returns

ProcessBranch that is sampled

Definition at line 341 of file action.h.

                                                    {
    const auto &log = logger<LogArea::Action>();
    double random_weight = random::uniform(0., total_weight);
    double weight_sum = 0.;
    /* Loop through all subprocesses and select one by Monte Carlo, based on
     * their weights.  */
    for (const auto &proc : subprocesses) {
      weight_sum += proc->weight();
      if (random_weight <= weight_sum) {
        /* Return the full process information. */
        return proc.get();
      }
    }
    /* Should never get here. */
    log.fatal(source_location,
              "Problem in choose_channel: ", subprocesses.size(), " ",
              weight_sum, " ", total_weight, " ",
              //          random_weight, "\n", *this);
              random_weight, "\n");
    std::abort();
  }

Here is the call graph for this function:

std::pair< double, double > smash::Action::sample_masses ( double  kinetic_energy_cm ) const

protectedvirtual

Sample final-state masses in general X->2 processes (thus also fixing the absolute c.o.m.

momentum).

Parameters

[in]

kinetic_energy_cm

total kinetic energy of the outgoing particles in their center of mass frame [GeV]

Exceptions

InvalidResonanceFormation

Returns

masses of final state particles

Reimplemented in smash::DecayAction.

Definition at line 154 of file action.cc.

                                          {
  const ParticleType &t_a = outgoing_particles_[0].type();
  const ParticleType &t_b = outgoing_particles_[1].type();
  // start with pole masses
  std::pair<double, double> masses = {t_a.mass(), t_b.mass()};

  if (kinetic_energy_cm < t_a.min_mass_kinematic() + t_b.min_mass_kinematic()) {
    const std::string reaction = incoming_particles_[0].type().name() +
                                 incoming_particles_[1].type().name() + "→" +
                                 t_a.name() + t_b.name();
    throw InvalidResonanceFormation(
        reaction + ": not enough energy, " + std::to_string(kinetic_energy_cm) +
        " < " + std::to_string(t_a.min_mass_kinematic()) + " + " +
        std::to_string(t_b.min_mass_kinematic()));
  }

  /* If one of the particles is a resonance, sample its mass. */
  if (!t_a.is_stable() && t_b.is_stable()) {
    masses.first = t_a.sample_resonance_mass(t_b.mass(), kinetic_energy_cm);
  } else if (!t_b.is_stable() && t_a.is_stable()) {
    masses.second = t_b.sample_resonance_mass(t_a.mass(), kinetic_energy_cm);
  } else if (!t_a.is_stable() && !t_b.is_stable()) {
    // two resonances in final state
    masses = t_a.sample_resonance_masses(t_b, kinetic_energy_cm);
  }
  return masses;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void smash::Action::sample_angles ( std::pair< double, double >  masses, double  kinetic_energy_cm  )

protectedvirtual

Sample final-state momenta in general X->2 processes (here: using an isotropical angular distribution).

Parameters

[in]	kinetic_energy_cm	total kinetic energy of the outgoing particles in their center of mass frame [GeV]
[in]	masses	masses of each of the final state particles

Reimplemented in smash::ScatterAction.

Definition at line 183 of file action.cc.

                                                           {
  const auto &log = logger<LogArea::Action>();

  ParticleData *p_a = &outgoing_particles_[0];
  ParticleData *p_b = &outgoing_particles_[1];

  const double pcm = pCM(kinetic_energy_cm, masses.first, masses.second);
  if (!(pcm > 0.0)) {
    log.warn("Particle: ", p_a->pdgcode(), " radial momentum: ", pcm);
    log.warn("Ektot: ", kinetic_energy_cm, " m_a: ", masses.first,
             " m_b: ", masses.second);
  }
  /* Here we assume an isotropic angular distribution. */
  Angles phitheta;
  phitheta.distribute_isotropically();

  p_a->set_4momentum(masses.first, phitheta.threevec() * pcm);
  p_b->set_4momentum(masses.second, -phitheta.threevec() * pcm);
  /* Debug message is printed before boost, so that p_a and p_b are
   * the momenta in the center of mass frame and thus opposite to
   * each other.*/
  log.debug("p_a: ", *p_a, "\np_b: ", *p_b);
}

Here is the call graph for this function:

Here is the caller graph for this function:

void smash::Action::sample_2body_phasespace ( )

protected

Sample the full 2-body phase-space (masses, momenta, angles) in the center-of-mass frame for the final state particles.

Definition at line 208 of file action.cc.

                                     {
  /* This function only operates on 2-particle final states. */
  assert(outgoing_particles_.size() == 2);
  const FourVector p_tot = total_momentum_of_outgoing_particles();
  const double cm_kin_energy = p_tot.abs();
  // first sample the masses
  const std::pair<double, double> masses = sample_masses(cm_kin_energy);
  // after the masses are fixed (and thus also pcm), sample the angles
  sample_angles(masses, cm_kin_energy);
}

Here is the call graph for this function:

Here is the caller graph for this function:

void smash::Action::sample_3body_phasespace ( )

protectedvirtual

Sample the full 3-body phase-space (masses, momenta, angles) in the center-of-mass frame for the final state particles.

Reimplemented in smash::DecayActionDilepton.

Definition at line 219 of file action.cc.

                                     {
  assert(outgoing_particles_.size() == 3);
  const double m_a = outgoing_particles_[0].type().mass(),
               m_b = outgoing_particles_[1].type().mass(),
               m_c = outgoing_particles_[2].type().mass();
  const double sqrts = sqrt_s();
  // sample mab from pCM(sqrt, mab, mc) pCM (mab, ma, mb) <= sqrts^2/4
  double mab, r, probability, pcm_ab, pcm;
  do {
    mab = random::uniform(m_a + m_b, sqrts - m_c);
    r = random::canonical();
    pcm = pCM(sqrts, mab, m_c);
    pcm_ab = pCM(mab, m_a, m_b);
    probability = pcm * pcm_ab * 4 / (sqrts * sqrts);
  } while (r > probability);
  Angles phitheta;
  phitheta.distribute_isotropically();
  outgoing_particles_[2].set_4momentum(m_c, pcm * phitheta.threevec());
  const ThreeVector beta_cm =
      pcm * phitheta.threevec() / std::sqrt(pcm * pcm + mab * mab);

  phitheta.distribute_isotropically();
  outgoing_particles_[0].set_4momentum(m_a, pcm_ab * phitheta.threevec());
  outgoing_particles_[1].set_4momentum(m_b, -pcm_ab * phitheta.threevec());
  outgoing_particles_[0].boost_momentum(beta_cm);
  outgoing_particles_[1].boost_momentum(beta_cm);
}

Here is the call graph for this function:

Here is the caller graph for this function:

const ParticleType& smash::Action::type_of_pout ( const ParticleData &  p_out ) const

inlineprivate

Get the type of a given particle.

Parameters

[in]

p_out

particle of which the type will be returned

Returns

type of given particle

Definition at line 425 of file action.h.

                                                                    {
    return p_out.type();
  }

Here is the call graph for this function:

Here is the caller graph for this function:

const ParticleType& smash::Action::type_of_pout ( const ParticleTypePtr &  p_out ) const

inlineprivate

Get the particle type for given pointer to a particle type.

Helper function for total_momentum_of_outgoing_particles

Parameters

[in]

p_out

pointer to a particle type

Returns

particle type

Definition at line 436 of file action.h.

                                                                       {
    return *p_out;
  }

Member Data Documentation

ParticleList smash::Action::incoming_particles_

protected

List with data of incoming particles.

Definition at line 303 of file action.h.

ParticleList smash::Action::outgoing_particles_

protected

Initially this stores only the PDG codes of final-state particles.

After perform was called it contains the complete particle data of the outgoing particles.

Definition at line 311 of file action.h.

const double smash::Action::time_of_execution_

protected

Time at which the action is supposed to be performed (absolute time in the lab frame in fm/c).

Definition at line 317 of file action.h.

ProcessType smash::Action::process_type_

protected

type of process

Definition at line 320 of file action.h.

---

The documentation for this class was generated from the following files:

• src/include/smash/action.h
• src/action.cc