#include <thermalizationaction.h>
ThermalizationAction implements forced thermalization as an Action class.
Particles before thermalization are treated as incoming, after thermalization as outgoing. This is an N->M action.
Definition at line 24 of file thermalizationaction.h.
Public Member Functions | |
ThermalizationAction (const GrandCanThermalizer &gct, double absolute_labframe_time) | |
The inherited class. More... | |
void | generate_final_state () override |
No need to do anything, because outgoing particles are set in constructor. More... | |
double | get_total_weight () const override |
Return the total weight value, which is mainly used for the weight output entry. More... | |
double | get_partial_weight () const override |
Return the specific weight for the chosen outgoing channel, which is mainly used for the partial weight output entry. More... | |
bool | any_particles_thermalized () const |
This method checks, if there are particles in the region to be thermalized. More... | |
void | format_debug_output (std::ostream &out) const override |
Function for debug output of incoming and outgoing particles from thermalization action. More... | |
Public Member Functions inherited from smash::Action | |
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 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 double | 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 std::vector< Particles > &ensembles, 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... | |
virtual double | 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... | |
void | set_stochastic_pos_idx () |
Setter function that stores a random incoming particle index latter used to determine the interaction point. More... | |
Additional Inherited Members | |
Static Public Member Functions inherited from smash::Action | |
static double | lambda_tilde (double a, double b, double c) |
Little helper function that calculates the lambda function (sometimes written with a tilde to better distinguish it) that appears e.g. More... | |
static void | sample_manybody_phasespace_impl (double sqrts, const std::vector< double > &m, std::vector< FourVector > &sampled_momenta) |
Implementation of the full n-body phase-space sampling (masses, momenta, angles) in the center-of-mass frame for the final state particles. More... | |
Protected Member Functions inherited from smash::Action | |
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_manybody_phasespace () |
Sample the full n-body phase-space (masses, momenta, angles) in the center-of-mass frame for the final state particles. More... | |
void | assign_formation_time_to_outgoing_particles () |
Assign the formation time to the outgoing particles. More... | |
Protected Attributes inherited from smash::Action | |
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). More... | |
ProcessType | process_type_ |
type of process More... | |
double | box_length_ = -1.0 |
Box length: needed to determine coordinates of collision correctly in case of collision through the wall. More... | |
int | stochastic_position_idx_ = -1 |
This stores a randomly-chosen index to an incoming particle. More... | |
smash::ThermalizationAction::ThermalizationAction | ( | const GrandCanThermalizer & | gct, |
double | absolute_labframe_time | ||
) |
The inherited class.
[in] | gct | The thermalization object taking care of removing and sampling new particles |
[in] | absolute_labframe_time | Current time in the computational frame. |
Definition at line 12 of file thermalizationaction.cc.
|
inlineoverridevirtual |
No need to do anything, because outgoing particles are set in constructor.
Implements smash::Action.
Definition at line 35 of file thermalizationaction.h.
|
inlineoverridevirtual |
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.
Implements smash::Action.
Definition at line 36 of file thermalizationaction.h.
|
inlineoverridevirtual |
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.
Implements smash::Action.
Definition at line 37 of file thermalizationaction.h.
|
inline |
This method checks, if there are particles in the region to be thermalized.
Definition at line 39 of file thermalizationaction.h.
|
inlineoverridevirtual |
Function for debug output of incoming and outgoing particles from thermalization action.
[in] | out | Location of the output stream |
Implements smash::Action.
Definition at line 47 of file thermalizationaction.h.