7 #ifndef SRC_INCLUDE_GRANDCAN_THERMALIZER_H_ 8 #define SRC_INCLUDE_GRANDCAN_THERMALIZER_H_ 94 double nb()
const {
return nb_; }
96 double ns()
const {
return ns_; }
98 double e()
const {
return e_; }
100 double p()
const {
return p_; }
104 double T()
const {
return T_; }
222 const std::array<int, 3> n_cells,
223 const std::array<double, 3> origin,
bool periodicity,
224 double e_critical,
double t_start,
double delta_t,
228 const std::array<double, 3> lat_sizes,
229 const std::array<double, 3> origin,
bool periodicity)
231 lat_sizes, conf.take({
"Cell_Number"}), origin, periodicity,
232 conf.
take({
"Critical_Edens"}), conf.take({
"Start_Time"}),
233 conf.take({
"Timestep"}),
255 bool ignore_cells_under_threshold =
true);
313 template <
typename F>
319 const double gamma = 1.0 / std::sqrt(1.0 - cell.
v().
sqr());
322 if (condition(i->strangeness(), i->baryon_number(), i->charge())) {
344 template <
typename F>
349 double partial_sum = 0.0;
350 int index_only_thermalized = -1;
351 while (partial_sum < r) {
352 index_only_thermalized++;
353 partial_sum +=
N_in_cells_[index_only_thermalized];
358 const double gamma = 1.0 / std::sqrt(1.0 - cell.
v().
sqr());
359 const double N_in_cell =
N_in_cells_[index_only_thermalized];
365 if (!condition(i->strangeness(), i->baryon_number(), i->charge())) {
378 const double m = type_to_sample->
mass();
385 particle.
set_4momentum(m, phitheta.threevec() * momentum_radial);
431 ParticleTypePtrList res;
434 res.push_back(&ptype);
444 const int B =
eos_typelist_[typelist_index]->baryon_number();
468 std::unique_ptr<RectangularLattice<ThermLatticeNode>>
lat_;
512 #endif // SRC_INCLUDE_GRANDCAN_THERMALIZER_H_ double p() const
Get pressure in the rest frame.
const double e_crit_
Critical energy density above which cells are thermalized.
void add_particle(const ParticleData &p, double factor)
Add particle contribution to Tmu0, nb and ns May look like unused at first glance, but it is actually used by update_lattice, where the node type of the lattice is templated.
A class to pre-calculate and store parameters relevant for density calculation.
ThermalizationAlgorithm
Defines the algorithm used for the forced thermalization.
GrandCanThermalizer(const std::array< double, 3 > lat_sizes, const std::array< int, 3 > n_cells, const std::array< double, 3 > origin, bool periodicity, double e_critical, double t_start, double delta_t, ThermalizationAlgorithm algo)
Default constructor for the GranCanThermalizer to allocate the lattice.
const size_t N_sorts_
Number of different species to be sampled.
ParticleList particles_to_remove() const
List of particles to be removed from the simulation.
The ThreeVector class represents a physical three-vector with the components .
void sample_multinomial(HadronClass particle_class, int N)
The sample_multinomial function samples integer numbers n_i distributed according to the multinomial ...
ThermLatticeNode()
Default constructor of thermal quantities on the lattice returning thermodynamic quantities in comput...
double mub() const
Get the net baryon chemical potential.
Class to handle the equation of state (EoS) of the hadron gas, consisting of all hadrons included int...
FourVector Tmu0_
Four-momentum flow of the cell.
void thermalize(const Particles &particles, double time, int ntest)
Main thermalize function, that chooses the algorithm to follow (BF or mode sampling).
double timestep_duration() const
bool is_time_to_thermalize(const Clock &clock) const
Check that the clock is close to n * period of thermalization, since the thermalization only happens ...
void add_particle_for_derivatives(const ParticleData &, double, ThreeVector)
dummy function for update_lattice
void thermalize_BF_algo(QuantumNumbers &conserved_initial, double time, int ntest)
Samples particles according to the BF algorithm by making use of the.
ParticleTypePtrList list_eos_particles() const
Extracts the particles in the hadron gas equation of state from the complete list of particle types i...
const double period_
Defines periodicity of the lattice in fm.
double ns() const
Get net strangeness density of the cell in the computational frame.
Interface to the SMASH configuration files.
static bool is_eos_particle(const ParticleType &ptype)
Check if a particle belongs to the EoS.
const ThermalizationAlgorithm algorithm_
Algorithm to choose for sampling of particles obeying conservation laws.
void set_formation_time(const double &form_time)
Set the absolute formation time.
double nb() const
Get net baryon density of the cell in the computational frame.
void print_statistics(const Clock &clock) const
Generates standard output with information about the thermodynamic properties of the lattice...
void set_rest_frame_quantities(double T0, double mub0, double mus0, const ThreeVector v0)
Set all the rest frame quantities to some values, this is useful for testing.
A container class to hold all the arrays on the lattice and access them.
void update_thermalizer_lattice(const Particles &particles, const DensityParameters &par, bool ignore_cells_under_threshold=true)
Compute all the thermodynamical quantities on the lattice from particles.
ThreeVector v_
Velocity of the rest frame.
double p_
Pressure in the rest frame.
static const ParticleTypeList & list_all()
HadronClass get_class(size_t typelist_index) const
Defines the class of hadrons by quantum numbers.
double ns_
Net strangeness density of the cell in the computational frame.
double mub_
Net baryon chemical potential.
GrandCanThermalizer(Configuration &conf, const std::array< double, 3 > lat_sizes, const std::array< double, 3 > origin, bool periodicity)
ParticleList sampled_list_
Newly generated particles by thermalizer.
Particle type contains the static properties of a particle species.
void set_4momentum(const FourVector &momentum_vector)
Set the particle's 4-momentum directly.
Non-strange mesons (S = 0) with electric cherge Q < 0.
A container for storing conserved values.
void compute_rest_frame_quantities(HadronGasEos &eos)
Temperature, chemical potentials and rest frame velocity are calculated given the hadron gas equation...
double e() const
Get energy density in the rest frame.
double sample_momenta_from_thermal(const double temperature, const double mass)
Samples a momentum from the Maxwell-Boltzmann (thermal) distribution in a faster way, given by Scott Pratt.
void compute_N_in_cells_mode_algo(F &&condition)
Computes average number of particles in each cell for the mode algorithm.
std::vector< size_t > cells_to_sample_
Cells above critical energy density.
Mesons with strangeness S < 0.
double N_total_in_cells_
Total number of particles over all cells in thermalization region.
void thermalize_mode_algo(QuantumNumbers &conserved_initial, double time)
Samples particles to the according to the mode algorithm.
The GrandCanThermalizer class implements the following functionality:
void renormalize_momenta(ParticleList &plist, const FourVector required_total_momentum)
Changes energy and momenta of the particles in plist to match the required_total_momentum.
Clock tracks the time in the simulation.
Neutral non-strange mesons.
Non-strange mesons (S = 0) with electric cherge Q > 0.
ParticleData sample_in_random_cell_mode_algo(const double time, F &&condition)
Samples one particle and the species, cell, momentum and coordinate are chosen from the corresponding...
std::vector< double > mult_sort_
Real number multiplicity for each particle type.
ParticleList to_remove_
Particles to be removed after this thermalization step.
double e_crit() const
Get the critical energy density.
Mesons with strangeness S > 0.
std::vector< double > N_in_cells_
Number of particles to be sampled in one cell.
double nb_
Net baryon density of the cell in the computational frame.
ThreeVector uniform_in_cell() const
RectangularLattice< ThermLatticeNode > & lattice() const
Getter function for the lattice.
static double partial_density(const ParticleType &ptype, double T, double mub, double mus)
Compute partial density of one hadron sort.
FourVector Tmu0() const
Get Four-momentum flow of the cell.
void set_4position(const FourVector &pos)
Set the particle's 4-position directly.
double e_
Energy density in the rest frame.
The ThermLatticeNode class is intended to compute thermodynamical quantities in a cell given a set of...
HadronClass
Specifier to classify the different hadron species according to their quantum numbers.
std::vector< int > mult_int_
Integer multiplicity for each particle type.
double mus_
Net strangeness chemical potential.
Value take(std::initializer_list< const char *> keys)
The default interface for SMASH to read configuration values.
void sample_in_random_cell_BF_algo(ParticleList &plist, const double time, size_t type_index)
The total number of particles of species type_index is defined by mult_int_ array that is returned by...
double T() const
Get the temperature.
double mus() const
Get the net strangeness chemical potential.
std::array< double, 7 > mult_classes_
The different hadron species according to the enum defined in.
A pointer-like interface to global references to ParticleType objects.
Angles provides a common interface for generating directions: i.e., two angles that should be interpr...
std::unique_ptr< RectangularLattice< ThermLatticeNode > > lat_
The lattice on which the thermodynamic quantities are calculated.
void distribute_isotropically()
Populate the object with a new direction.
The Particles class abstracts the storage and manipulation of particles.
double current_time() const
std::ostream & operator<<(std::ostream &out, const ActionPtr &action)
Convenience: dereferences the ActionPtr to Action.
ParticleList particles_to_insert() const
List of newly created particles to be inserted in the simulation.
The FourVector class holds relevant values in Minkowski spacetime with (+, −, −, −) metric signature.
double cell_volume_
Volume of a single cell, necessary to convert thermal densities to actual particle numbers...
HadronGasEos eos_
Hadron gas equation of state.
ParticleData contains the dynamic information of a certain particle.
void boost_momentum(const ThreeVector &v)
Apply a Lorentz-boost to only the momentum.
const ParticleTypePtrList eos_typelist_
List of particle types from which equation of state is computed.
double mult_class(const HadronClass cl) const
const double t_start_
Starting time of the simulation.
ThreeVector v() const
Get 3-velocity of the rest frame.