smash::decaytree Namespace Reference

Classes
struct	Node
	Node of a decay tree, representing a possible action (2-to-2 or 1-to-2). More...

Functions
static std::string	make_decay_name (const std::string &res_name, const DecayBranchPtr &decay, ParticleTypePtrList &final_state)
	Generate name for decay and update final state. More...
static void	add_decays (Node &node)
	Add nodes for all decays possible from the given node and all of its children. More...

Function Documentation

make_decay_name()

static std::string smash::decaytree::make_decay_name ( const std::string &  res_name, const DecayBranchPtr &  decay, ParticleTypePtrList &  final_state  )

static

Generate name for decay and update final state.

Parameters

[in]	res_name	Name of resonance.
[in]	decay	Decay branch.
[out]	final_state	Final state of decay.

Returns

Name of decay.

Definition at line 660 of file scatteractionsfinder.cc.

                                                                     {
  std::stringstream name;
  name << "[" << res_name << "->";
  for (const auto& p : decay->particle_types()) {
    name << p->name();
    final_state.push_back(p);
  }
  name << "]";
  return name.str();
}

Here is the caller graph for this function:

add_decays()

static void smash::decaytree::add_decays ( Node &  node )

static

Add nodes for all decays possible from the given node and all of its children.

Parameters

node	Starting node.

Definition at line 679 of file scatteractionsfinder.cc.

                                   {
  // If there is more than one unstable particle in the current state, then
  // there will be redundant paths in the decay tree, corresponding to
  // reorderings of the decays. To avoid double counting, we normalize by the
  // number of possible decay orderings. Normalizing by the number of unstable
  // particles recursively corresponds to normalizing by the factorial that
  // gives the number of reorderings.
  //
  // Ideally, the redundant paths should never be added to the decay tree, but
  // we never have more than two redundant paths, so it probably does not matter
  // much.
  uint32_t n_unstable = 0;
  for (const ParticleTypePtr ptype : node.state_) {
    if (!ptype->is_stable()) {
      n_unstable += 1;
    }
  }
  const double norm =
      n_unstable != 0 ? 1. / static_cast<double>(n_unstable) : 1.;

  for (const ParticleTypePtr ptype : node.state_) {
    if (!ptype->is_stable()) {
      for (const auto& decay : ptype->decay_modes().decay_mode_list()) {
        ParticleTypePtrList parts;
        const auto name = make_decay_name(ptype->name(), decay, parts);
        auto& new_node = node.add_action(name, norm * decay->weight(), {ptype},
                                         std::move(parts));
        add_decays(new_node);
      }
    }
  }
}

Here is the call graph for this function:

Here is the caller graph for this function: