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 37 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 Particles &particles, const int event_number) override
	writes the event header More...
void	at_eventend (const Particles &particles, const int event_number, double impact_parameter, bool empty_event) override
	only flushes the output the file More...
void	at_intermediate_time (const Particles &particles, 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_interaction (const Action &action, const double density)
	Called whenever an action modified one or more particles. 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_output (const GrandCanThermalizer &gct)
	Output to write energy-momentum tensor and related quantities from the thermalizer class. 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

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 109 of file thermodynamicoutput.cc.

    : OutputInterface(name),
      file_{path / "thermodynamics.dat", "w"},
      out_par_(out_par) {
  std::fprintf(file_.get(), "# %s thermodynamics output\n", VERSION_MAJOR);
  const ThreeVector r = out_par.td_position;
  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:

smash::ThermodynamicOutput::~ThermodynamicOutput

(

)

Default destructor.

Definition at line 161 of file thermodynamicoutput.cc.

{}

Member Function Documentation

void smash::ThermodynamicOutput::at_eventstart ( const Particles &  particles, const int  event_number  )

overridevirtual

writes the event header

Parameters

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

Implements smash::OutputInterface.

Definition at line 163 of file thermodynamicoutput.cc.

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

Here is the call graph for this function:

void smash::ThermodynamicOutput::at_eventend ( const Particles &  particles, const int  event_number, double  impact_parameter, bool  empty_event  )

overridevirtual

only flushes the output the file

Parameters

[in]	particles	Dummy, is just here to satisfy inheritance
[in]	event_number	Dummy, is just here to satisfy inheritance
[in]	impact_parameter	Dummy, is just here to satisfy inheritance
[in]	empty_event	Dummy, is just here to satisfy inheritance

Implements smash::OutputInterface.

Definition at line 168 of file thermodynamicoutput.cc.

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

Here is the call graph for this function:

void smash::ThermodynamicOutput::at_intermediate_time ( const Particles &  particles, 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]	particles	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 175 of file thermodynamicoutput.cc.

                                         {
  std::fprintf(file_.get(), "%6.2f ", clock->current_time());
  constexpr bool compute_gradient = false;
  if (out_par_.td_rho_eckart) {
    const double rho = std::get<0>(current_eckart(
        out_par_.td_position, particles, dens_param, out_par_.td_dens_type,
        compute_gradient, out_par_.td_smearing));
    std::fprintf(file_.get(), "%7.4f ", rho);
  }
  if (out_par_.td_tmn || out_par_.td_tmn_landau || out_par_.td_v_landau) {
    EnergyMomentumTensor Tmn;
    for (const auto &p : particles) {
      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(), "%7.4f ", Tmn_L[i]);
      }
    }
    if (out_par_.td_v_landau) {
      std::fprintf(file_.get(), "%7.4f %7.4f %7.4f ", -u[1] / u[0],
                   -u[2] / u[0], -u[3] / u[0]);
    }
  }
  if (out_par_.td_jQBS) {
    FourVector jQ = std::get<1>(current_eckart(
        out_par_.td_position, particles, dens_param, DensityType::Charge,
        compute_gradient, out_par_.td_smearing));
    FourVector jB = std::get<1>(current_eckart(
        out_par_.td_position, particles, dens_param, DensityType::Baryon,
        compute_gradient, out_par_.td_smearing));
    FourVector 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:

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 245 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

RenamingFilePtr smash::ThermodynamicOutput::file_

private

Pointer to output file.

Definition at line 100 of file thermodynamicoutput.h.

const OutputParameters smash::ThermodynamicOutput::out_par_

private

Structure that holds all the information about what to printout.

Definition at line 102 of file thermodynamicoutput.h.

---

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

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