smash::ThermodynamicOutput Class Reference

#include <thermodynamicoutput.h>

Writes the thermodynamic quantities at a specified point versus time.

This class is a temporary solution to write thermodynamic quantities out. Calculations are called directly inside the output functions. In future it should be substituted by some more general output.

Definition at line 39 of file thermodynamicoutput.h.

Inheritance diagram for smash::ThermodynamicOutput:

Collaboration diagram for smash::ThermodynamicOutput:

Public Member Functions
	ThermodynamicOutput (const bf::path &path, const std::string &name, const OutputParameters &out_par)
	Construct Output param[in] path Path to output param[in] name Filename param[in] out_par Parameters of output. More...
	~ThermodynamicOutput ()
	Default destructor. More...
void	at_eventstart (const std::vector< Particles > &ensembles, const int event_number) override
	writes the event header More...
void	at_eventend (const std::vector< Particles > &ensembles, const int event_number) override
	only flushes the output the file More...
void	at_intermediate_time (const std::vector< Particles > &ensembles, const std::unique_ptr< Clock > &clock, const DensityParameters &dens_param) override
	Writes thermodynamics every fixed time interval. More...
void	density_along_line (const char *file_name, const ParticleList &plist, const DensityParameters &param, DensityType dens_type, const ThreeVector &line_start, const ThreeVector &line_end, int n_points)
	Prints density along the specified line. More...
Public Member Functions inherited from smash::OutputInterface
	OutputInterface (std::string name)
	Construct output interface. More...
virtual	~OutputInterface ()=default
virtual void	at_eventstart (const Particles &particles, const int event_number, const EventInfo &info)
	Output launched at event start after initialization, when particles are generated but not yet propagated. More...
virtual void	at_eventstart (const int event_number, const ThermodynamicQuantity tq, const DensityType dens_type, RectangularLattice< DensityOnLattice > lattice)
	Output launched at event start after initialization, when particles are generated but not yet propagated. More...
virtual void	at_eventstart (const int event_number, const ThermodynamicQuantity tq, const DensityType dens_type, RectangularLattice< EnergyMomentumTensor > lattice)
	Output launched atevent start after initialization, when particles are generated but not yet propagated. More...
virtual void	at_eventend (const int event_number, const ThermodynamicQuantity tq, const DensityType dens_type)
	Output launched at event end. More...
virtual void	at_eventend (const ThermodynamicQuantity tq)
	Output launched at event end. More...
virtual void	at_eventend (const Particles &particles, const int event_number, const EventInfo &info)
	Output launched at event end. More...
virtual void	at_interaction (const Action &action, const double density)
	Called whenever an action modified one or more particles. More...
virtual void	at_intermediate_time (const Particles &particles, const std::unique_ptr< Clock > &clock, const DensityParameters &dens_param, const EventInfo &info)
	Output launched after every N'th timestep. More...
virtual void	thermodynamics_output (const ThermodynamicQuantity tq, const DensityType dt, RectangularLattice< DensityOnLattice > &lattice)
	Output to write thermodynamics from the lattice. More...
virtual void	thermodynamics_output (const ThermodynamicQuantity tq, const DensityType dt, RectangularLattice< EnergyMomentumTensor > &lattice)
	Output to write energy-momentum tensor and related quantities from the lattice. More...
virtual void	thermodynamics_lattice_output (RectangularLattice< DensityOnLattice > &lattice, const double current_time)
	Output to write thermodynamics from the lattice. More...
virtual void	thermodynamics_lattice_output (RectangularLattice< DensityOnLattice > &lattice, const double current_time, const std::vector< Particles > &ensembles, const DensityParameters &dens_param)
	Output to write thermodynamics from the lattice. More...
virtual void	thermodynamics_lattice_output (const ThermodynamicQuantity tq, RectangularLattice< EnergyMomentumTensor > &lattice, const double current_time)
	Output to write energy-momentum tensor and related quantities from the lattice. More...
virtual void	thermodynamics_output (const GrandCanThermalizer &gct)
	Output to write energy-momentum tensor and related quantities from the thermalizer class. More...
virtual void	fields_output (const std::string name1, const std::string name2, RectangularLattice< std::pair< ThreeVector, ThreeVector >> &lat)
	Write fields in vtk output. More...
bool	is_dilepton_output () const
	Get, whether this is the dilepton output? More...
bool	is_photon_output () const
	Get, whether this is the photon output? More...
bool	is_IC_output () const
	Get, whether this is the IC output? More...
const char *	to_string (const ThermodynamicQuantity tq)
	Convert thermodynamic quantities to strings. More...
const char *	to_string (const DensityType dens_type)
	Convert density types to strings. More...

Private Attributes
RenamingFilePtr	file_
	Pointer to output file. More...
const OutputParameters	out_par_
	Structure that holds all the information about what to printout. More...

Additional Inherited Members
Protected Attributes inherited from smash::OutputInterface
const bool	is_dilepton_output_
	Is this the dilepton output? More...
const bool	is_photon_output_
	Is this the photon output? More...
const bool	is_IC_output_
	Is this the IC output? More...

Constructor & Destructor Documentation

ThermodynamicOutput()

smash::ThermodynamicOutput::ThermodynamicOutput	(	const bf::path &	path,
		const std::string &	name,
		const OutputParameters &	out_par
	)

Construct Output param[in] path Path to output param[in] name Filename param[in] out_par Parameters of output.

Definition at line 116 of file thermodynamicoutput.cc.

    : OutputInterface(name),
      file_{path / "thermodynamics.dat", "w"},
      out_par_(out_par) {
  std::fprintf(file_.get(), "# %s thermodynamics output\n", SMASH_VERSION);
  const ThreeVector r = out_par.td_position;
  if (out_par_.td_only_participants) {
    std::fprintf(file_.get(), "# only participants are taken into account\n");
  }
  if (out_par_.td_smearing) {
    std::fprintf(file_.get(), "# @ point (%6.2f, %6.2f, %6.2f) [fm]\n", r.x1(),
                 r.x2(), r.x3());
  } else {
    std::fprintf(file_.get(), "# averaged over the entire volume\n");
  }
  std::fprintf(file_.get(), "# %s\n", to_string(out_par.td_dens_type));
  std::fprintf(file_.get(), "# time [fm/c], ");
  if (out_par_.td_rho_eckart) {
    std::fprintf(file_.get(), "%s [fm^-3], ",
                 to_string(ThermodynamicQuantity::EckartDensity));
  }
  if (out_par_.td_tmn) {
    if (out_par_.td_smearing) {
      std::fprintf(file_.get(), "%s [GeV/fm^3] 00 01 02 03 11 12 13 22 23 33, ",
                   to_string(ThermodynamicQuantity::Tmn));
    } else {
      std::fprintf(file_.get(), "%s [GeV] 00 01 02 03 11 12 13 22 23 33, ",
                   to_string(ThermodynamicQuantity::Tmn));
    }
  }
  if (out_par_.td_tmn_landau) {
    if (out_par_.td_smearing) {
      std::fprintf(file_.get(), "%s [GeV/fm^3] 00 01 02 03 11 12 13 22 23 33, ",
                   to_string(ThermodynamicQuantity::TmnLandau));
    } else {
      std::fprintf(file_.get(), "%s [GeV] 00 01 02 03 11 12 13 22 23 33, ",
                   to_string(ThermodynamicQuantity::TmnLandau));
    }
  }
  if (out_par_.td_v_landau) {
    std::fprintf(file_.get(), "%s x y z, ",
                 to_string(ThermodynamicQuantity::LandauVelocity));
  }
  if (out_par_.td_jQBS) {
    if (out_par_.td_smearing) {
      std::fprintf(file_.get(), "j_QBS [(Q,B,S)/fm^3] (0 1 2 3)x3");
    } else {
      std::fprintf(file_.get(), "j_QBS [(Q,B,S)] (0 1 2 3)x3");
    }
  }
  std::fprintf(file_.get(), "\n");
}

Here is the call graph for this function:

~ThermodynamicOutput()

smash::ThermodynamicOutput::~ThermodynamicOutput

(

)

Default destructor.

Definition at line 171 of file thermodynamicoutput.cc.

{}

Member Function Documentation

at_eventstart()

void smash::ThermodynamicOutput::at_eventstart ( const std::vector< Particles > &  ensembles, const int  event_number  )

overridevirtual

writes the event header

Parameters

[in]	ensembles	Dummy, is just here to satisfy inheritance
[in]	event_number	Current event number, that will be written to the header

Reimplemented from smash::OutputInterface.

Definition at line 173 of file thermodynamicoutput.cc.

                                                            {
  std::fprintf(file_.get(), "# event %i\n", event_number);
}

Here is the call graph for this function:

at_eventend()

void smash::ThermodynamicOutput::at_eventend ( const std::vector< Particles > &  ensembles, const int  event_number  )

overridevirtual

only flushes the output the file

Parameters

[in]	ensembles	Dummy, is just here to satisfy inheritance
[in]	event_number	Dummy, is just here to satisfy inheritance

Reimplemented from smash::OutputInterface.

Definition at line 178 of file thermodynamicoutput.cc.

                                                 {
  std::fflush(file_.get());
}

Here is the call graph for this function:

at_intermediate_time()

void smash::ThermodynamicOutput::at_intermediate_time ( const std::vector< Particles > &  ensembles, const std::unique_ptr< Clock > &  clock, const DensityParameters &  dens_param  )

overridevirtual

Writes thermodynamics every fixed time interval.

For configuring the output see ASCII Thermodynamics Output.

Parameters

[in]	ensembles	Particles, from which thermodynamic variables are computed
[in]	clock	Clock, needed to get current time
[in]	dens_param	set of parameters, defining smearing. For more info about smearing see ASCII Thermodynamics Output.

Reimplemented from smash::OutputInterface.

Definition at line 183 of file thermodynamicoutput.cc.

                                                                            {
  std::fprintf(file_.get(), "%6.2f ", clock->current_time());
  constexpr bool compute_gradient = false;
  if (out_par_.td_rho_eckart) {
    FourVector jmu = FourVector();
    for (const Particles &particles : ensembles) {
      jmu += std::get<1>(current_eckart(
          out_par_.td_position, particles, dens_param, out_par_.td_dens_type,
          compute_gradient, out_par_.td_smearing));
    }
    std::fprintf(file_.get(), "%15.12f ", jmu.abs());
  }
  if (out_par_.td_tmn || out_par_.td_tmn_landau || out_par_.td_v_landau) {
    EnergyMomentumTensor Tmn;
    for (const Particles &particles : ensembles) {
      for (const auto &p : particles) {
        if (dens_param.only_participants()) {
          // if this condition holds, the hadron is a spectator and we skip it
          if (p.get_history().collisions_per_particle == 0) {
            continue;
          }
        }
        const double dens_factor =
            density_factor(p.type(), out_par_.td_dens_type);
        if (std::abs(dens_factor) < really_small) {
          continue;
        }
        if (out_par_.td_smearing) {
          const auto sf =
              unnormalized_smearing_factor(
                  p.position().threevec() - out_par_.td_position, p.momentum(),
                  1.0 / p.momentum().abs(), dens_param, compute_gradient)
                  .first;
          if (sf < really_small) {
            continue;
          }
          Tmn.add_particle(p, dens_factor * sf * dens_param.norm_factor_sf());
        } else {
          Tmn.add_particle(p, dens_factor);
        }
      }
    }
    const FourVector u = Tmn.landau_frame_4velocity();
    const EnergyMomentumTensor Tmn_L = Tmn.boosted(u);
    if (out_par_.td_tmn) {
      for (int i = 0; i < 10; i++) {
        std::fprintf(file_.get(), "%15.12f ", Tmn[i]);
      }
    }
    if (out_par_.td_tmn_landau) {
      for (int i = 0; i < 10; i++) {
        std::fprintf(file_.get(), "%15.12f ", Tmn_L[i]);
      }
    }
    if (out_par_.td_v_landau) {
      std::fprintf(file_.get(), "%15.12f %15.12f %15.12f ", -u[1] / u[0],
                   -u[2] / u[0], -u[3] / u[0]);
    }
  }
  if (out_par_.td_jQBS) {
    FourVector jQ = FourVector(), jB = FourVector(), jS = FourVector();
    for (const Particles &particles : ensembles) {
      jQ += std::get<1>(current_eckart(out_par_.td_position, particles,
                                       dens_param, DensityType::Charge,
                                       compute_gradient, out_par_.td_smearing));
      jB += std::get<1>(current_eckart(out_par_.td_position, particles,
                                       dens_param, DensityType::Baryon,
                                       compute_gradient, out_par_.td_smearing));
      jS += std::get<1>(current_eckart(out_par_.td_position, particles,
                                       dens_param, DensityType::Strangeness,
                                       compute_gradient, out_par_.td_smearing));
    }
    std::fprintf(file_.get(), "%15.12f %15.12f %15.12f %15.12f ", jQ[0], jQ[1],
                 jQ[2], jQ[3]);
    std::fprintf(file_.get(), "%15.12f %15.12f %15.12f %15.12f ", jB[0], jB[1],
                 jB[2], jB[3]);
    std::fprintf(file_.get(), "%15.12f %15.12f %15.12f %15.12f ", jS[0], jS[1],
                 jS[2], jS[3]);
  }
  std::fprintf(file_.get(), "\n");
}

Here is the call graph for this function:

density_along_line()

void smash::ThermodynamicOutput::density_along_line	(	const char *	file_name,
		const ParticleList &	plist,
		const DensityParameters &	param,
		DensityType	dens_type,
		const ThreeVector &	line_start,
		const ThreeVector &	line_end,
		int	n_points
	)

Prints density along the specified line.

Useful to make 1D plots of density profiles.

Parameters

[in]	file_name	name of the file to print out
[in]	param	Parameters for density calculation
[in]	plist	particles, from which density is computed
[in]	dens_type	type of density
[in]	line_start	starting point of the line
[in]	line_end	ending point of the line
[in]	n_points	number of points along the line, where density is printed out

Definition at line 267 of file thermodynamicoutput.cc.

                                                                              {
  ThreeVector r;
  std::ofstream a_file;
  a_file.open(file_name, std::ios::out);
  const bool compute_gradient = false;
  const bool smearing = true;
 
  for (int i = 0; i <= n_points; i++) {
    r = line_start + (line_end - line_start) * (1.0 * i / n_points);
    double rho_eck = std::get<0>(
        current_eckart(r, plist, param, dens_type, compute_gradient, smearing));
    a_file << r.x1() << " " << r.x2() << " " << r.x3() << " " << rho_eck
           << "\n";
  }
}

Here is the call graph for this function:

Member Data Documentation

file_

RenamingFilePtr smash::ThermodynamicOutput::file_

private

Pointer to output file.

Definition at line 103 of file thermodynamicoutput.h.

out_par_

const OutputParameters smash::ThermodynamicOutput::out_par_

private

Structure that holds all the information about what to printout.

Definition at line 105 of file thermodynamicoutput.h.

---

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

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