Output Classes

Classes
class	smash::BinaryOutputBase
	Base class for SMASH binary output. More...
class	smash::BinaryOutputCollisions
	Saves SMASH collision history to binary file. More...
class	smash::BinaryOutputParticles
	Writes the particle list at specific times to the binary file. More...
class	smash::BinaryOutputInitialConditions
	Writes the particles when crossing the hypersurface to the binary file. More...
class	smash::HepMcInterface
	Base class for output handlers that need the HepMC3 structure. More...
class	smash::HepMcOutput
	SMASH output to HepMC file. More...
class	smash::ICOutput
	SMASH output in ASCII format containing initial conditions for hydrodynamic codes. More...
class	smash::OscarOutput< Format, Contents >
struct	smash::EventInfo
	Structure to contain custom data for output. More...
class	smash::OutputInterface
	Abstraction of generic output. More...
class	smash::RivetOutput
	SMASH output to Rivet analyses. More...
class	smash::RootOutput
class	smash::ThermodynamicLatticeOutput
	Writes the thermodynamic quantities at lattice points versus time. More...
class	smash::ThermodynamicOutput
	Writes the thermodynamic quantities at a specified point versus time. More...
class	smash::VtkOutput
	SMASH output in a paraview format, intended for simple visualization. More...

Enumerations
enum	smash::OscarOutputFormat { smash::OscarFormat2013 , smash::OscarFormat2013Extended , smash::OscarFormat1999 }
	Selector for the output format of OscarOutput. More...
enum	smash::OscarOutputContents {   smash::OscarInteractions = 0x001 , smash::OscarTimesteps = 0x002 , smash::OscarAtEventstart = 0x004 , smash::OscarParticlesAtEventend = 0x008 ,   smash::OscarParticlesAtEventendIfNotEmpty = 0x010 , smash::OscarParticlesIC = 0x020 }
	Flags for the Contents template parameter of OscarOutput. More...

Functions
	smash::OscarOutput< Format, Contents >::OscarOutput (const std::filesystem::path &path, const std::string &name)
	Create oscar output. More...
void	smash::OscarOutput< Format, Contents >::at_eventstart (const Particles &particles, const int event_number, const EventInfo &event) override
	Writes the initial particle information of an event to the oscar output. More...
void	smash::OscarOutput< Format, Contents >::at_eventend (const Particles &particles, const int event_number, const EventInfo &) override
	Writes the final particle information of an event to the oscar output. More...
void	smash::OscarOutput< Format, Contents >::at_interaction (const Action &action, const double density) override
	Writes a interaction prefix line and a line for every incoming and outgoing particle to the oscar output. More...
void	smash::OscarOutput< Format, Contents >::at_intermediate_time (const Particles &particles, const std::unique_ptr< Clock > &clock, const DensityParameters &dens_param, const EventInfo &event) override
	Writes a prefix line then write out all current particles. More...
void	smash::OscarOutput< Format, Contents >::write_particledata (const ParticleData &data)
	Write single particle information line to output. More...
void	smash::OscarOutput< Format, Contents >::write (const Particles &particles)
	Write the particle information of a list of particles to the output. More...
std::unique_ptr< OutputInterface >	smash::create_oscar_output (const std::string &format, const std::string &content, const std::filesystem::path &path, const OutputParameters &out_par)

Variables
int	smash::OscarOutput< Format, Contents >::current_event_ = 0
	Keep track of event number. More...
RenamingFilePtr	smash::OscarOutput< Format, Contents >::file_
	Full filepath of the output file. More...

Enumeration Type Documentation

OscarOutputFormat

enum smash::OscarOutputFormat

Selector for the output format of OscarOutput.

Enumerator
OscarFormat2013
OscarFormat2013Extended
OscarFormat1999

Definition at line 29 of file oscaroutput.h.

                       {
  OscarFormat2013,
  OscarFormat2013Extended,
  OscarFormat1999
};

OscarOutputContents

enum smash::OscarOutputContents

Flags for the Contents template parameter of OscarOutput.

Flags can be combined with binary OR operators to some arbitrary int. That's why the values of the enumerators are written out (in hexadecimal), to ensure every flag occupies a single bit.

Enumerator
OscarInteractions	store interaction information (write_interaction)
OscarTimesteps	store the state after N timesteps (after_Nth_timestep)
OscarAtEventstart	store the state at the start of each event (at_eventstart)
OscarParticlesAtEventend	store the state at the end of each event (at_eventend)
OscarParticlesAtEventendIfNotEmpty	store the state at the end of each event if it is not empty (at_eventend)
OscarParticlesIC	store the particles that are removed on the hypersurface

Definition at line 42 of file oscaroutput.h.

                         {
  OscarInteractions = 0x001,
  OscarTimesteps = 0x002,
  OscarAtEventstart = 0x004,
  OscarParticlesAtEventend = 0x008,
  OscarParticlesAtEventendIfNotEmpty = 0x010,
  OscarParticlesIC = 0x020
};

Function Documentation

OscarOutput()

template<OscarOutputFormat Format, int Contents>

smash::OscarOutput< Format, Contents >::OscarOutput	(	const std::filesystem::path &	path,
		const std::string &	name
	)

Create oscar output.

Parameters

[in]	path	Output path.
[in]	name	Name of the ouput.

at_eventstart()

template<OscarOutputFormat Format, int Contents>

void smash::OscarOutput< Format, Contents >::at_eventstart ( const Particles &  particles, const int  event_number, const EventInfo &  event  )

overridevirtual

Writes the initial particle information of an event to the oscar output.

Parameters

[in]	particles	Current list of all particles.
[in]	event_number	Number of event.
[in]	event	Event info, see event_info

Reimplemented from smash::OutputInterface.

Definition at line 126 of file oscaroutput.cc.

                                                                     {
  current_event_ = event_number;
  if (Contents & OscarAtEventstart) {
    if (Format == OscarFormat2013 || Format == OscarFormat2013Extended) {
      std::fprintf(file_.get(), "# event %i in %zu\n", event_number,
                   particles.size());
    } else {
      /* OSCAR line prefix : initial particles; final particles; event id
       * First block of an event: initial = 0, final = number of particles
       */
      const size_t zero = 0;
      std::fprintf(file_.get(), "%zu %zu %i\n", zero, particles.size(),
                   event_number);
    }
    if (!(Contents & OscarParticlesIC)) {
      // We do not want the inital particle list to be printed in case of IC
      // output
      write(particles);
    }
  }
}

at_eventend()

template<OscarOutputFormat Format, int Contents>

void smash::OscarOutput< Format, Contents >::at_eventend ( const Particles &  particles, const int  event_number, const EventInfo &  event  )

overridevirtual

Writes the final particle information of an event to the oscar output.

Parameters

[in]	particles	Current list of particles.
[in]	event_number	Number of event.
[in]	event	Event info, see event_info

Reimplemented from smash::OutputInterface.

Definition at line 151 of file oscaroutput.cc.

                                                                        {
  if (Format == OscarFormat2013 || Format == OscarFormat2013Extended) {
    if (Contents & OscarParticlesAtEventend ||
        (Contents & OscarParticlesAtEventendIfNotEmpty && !event.empty_event)) {
      std::fprintf(file_.get(), "# event %i out %zu\n", event_number,
                   particles.size());
      write(particles);
    }
    // Comment end of an event
    const char *empty_event_str = event.empty_event ? "no" : "yes";
    std::fprintf(
        file_.get(),
        "# event %i end 0 impact %7.3f scattering_projectile_target %s\n",
        event_number, event.impact_parameter, empty_event_str);
  } else {
    /* OSCAR line prefix : initial particles; final particles; event id
     * Last block of an event: initial = number of particles, final = 0
     * Block ends with null interaction. */
    const size_t zero = 0;
    if (Contents & OscarParticlesAtEventend ||
        (Contents & OscarParticlesAtEventendIfNotEmpty && !event.empty_event)) {
      std::fprintf(file_.get(), "%zu %zu %i\n", particles.size(), zero,
                   event_number);
      write(particles);
    }
    // Null interaction marks the end of an event
    std::fprintf(file_.get(), "%zu %zu %i %7.3f\n", zero, zero, event_number,
                 event.impact_parameter);
  }
  // Flush to disk
  std::fflush(file_.get());
 
  if (Contents & OscarParticlesIC) {
    // If the runtime is too short some particles might not yet have
    // reached the hypersurface. Warning is printed.
    if (particles.size() != 0 && !event.impose_kinematic_cut_for_SMASH_IC) {
      logg[LHyperSurfaceCrossing].warn(
          "End time might be too small for initial conditions output. "
          "Hypersurface has not yet been crossed by ",
          particles.size(), " particle(s).");
    }
  }
}

at_interaction()

template<OscarOutputFormat Format, int Contents>

void smash::OscarOutput< Format, Contents >::at_interaction ( const Action &  action, const double  density  )

overridevirtual

Writes a interaction prefix line and a line for every incoming and outgoing particle to the oscar output.

Parameters

[in]	action	Action that holds the information of the interaction.
[in]	density	Density at the interaction point.

Reimplemented from smash::OutputInterface.

Definition at line 198 of file oscaroutput.cc.

                                                                         {
  if (Contents & OscarInteractions) {
    if (Format == OscarFormat2013 || Format == OscarFormat2013Extended) {
      std::fprintf(file_.get(),
                   "# interaction in %zu out %zu rho %12.7f weight %12.7g"
                   " partial %12.7f type %5i\n",
                   action.incoming_particles().size(),
                   action.outgoing_particles().size(), density,
                   action.get_total_weight(), action.get_partial_weight(),
                   static_cast<int>(action.get_type()));
    } else {
      /* OSCAR line prefix : initial final
       * particle creation: 0 1
       * particle 2<->2 collision: 2 2
       * resonance formation: 2 1
       * resonance decay: 1 2
       * etc.*/
      std::fprintf(file_.get(), "%zu %zu %12.7f %12.7f %12.7f %5i\n",
                   action.incoming_particles().size(),
                   action.outgoing_particles().size(), density,
                   action.get_total_weight(), action.get_partial_weight(),
                   static_cast<int>(action.get_type()));
    }
    for (const auto &p : action.incoming_particles()) {
      write_particledata(p);
    }
    for (const auto &p : action.outgoing_particles()) {
      write_particledata(p);
    }
  } else if (Contents & OscarParticlesIC) {
    for (const auto &p : action.incoming_particles()) {
      write_particledata(p);
    }
  }
}

at_intermediate_time()

template<OscarOutputFormat Format, int Contents>

void smash::OscarOutput< Format, Contents >::at_intermediate_time ( const Particles &  particles, const std::unique_ptr< Clock > &  clock, const DensityParameters &  dens_param, const EventInfo &  event  )

overridevirtual

Writes a prefix line then write out all current particles.

Parameters

[in]	particles	Current list of particles.
[in]	clock	Unused, needed since inherited.
[in]	dens_param	Unused, needed since inherited.
[in]	event	Event info, see event_info

Reimplemented from smash::OutputInterface.

Definition at line 236 of file oscaroutput.cc.

                                                  {
  if (Contents & OscarTimesteps) {
    if (Format == OscarFormat2013 || Format == OscarFormat2013Extended) {
      std::fprintf(file_.get(), "# event %i out %zu\n", current_event_,
                   particles.size());
    } else {
      const size_t zero = 0;
      std::fprintf(file_.get(), "%zu %zu %i\n", particles.size(), zero,
                   current_event_);
    }
    write(particles);
  }
}

write_particledata()

template<OscarOutputFormat Format, int Contents>

void smash::OscarOutput< Format, Contents >::write_particledata ( const ParticleData &  data )

private

Write single particle information line to output.

Parameters

[in]

data

Data of particle.

Definition at line 746 of file oscaroutput.cc.

                              {
  const FourVector pos = data.position();
  const FourVector mom = data.momentum();
  if (Format == OscarFormat2013) {
    std::fprintf(file_.get(), "%g %g %g %g %g %.9g %.9g %.9g %.9g %s %i %i\n",
                 pos.x0(), pos.x1(), pos.x2(), pos.x3(), data.effective_mass(),
                 mom.x0(), mom.x1(), mom.x2(), mom.x3(),
                 data.pdgcode().string().c_str(), data.id(),
                 data.type().charge());
  } else if (Format == OscarFormat2013Extended) {
    const auto h = data.get_history();
    std::fprintf(
        file_.get(),
        "%g %g %g %g %g %.9g %.9g %.9g"
        " %.9g %s %i %i %i %g %g %i %i %g %s %s %i %i\n",
        pos.x0(), pos.x1(), pos.x2(), pos.x3(), data.effective_mass(), mom.x0(),
        mom.x1(), mom.x2(), mom.x3(), data.pdgcode().string().c_str(),
        data.id(), data.type().charge(), h.collisions_per_particle,
        data.formation_time(), data.xsec_scaling_factor(), h.id_process,
        static_cast<int>(h.process_type), h.time_last_collision,
        h.p1.string().c_str(), h.p2.string().c_str(),
        data.type().baryon_number(), data.type().strangeness());
  } else {
    std::fprintf(file_.get(), "%i %s %i %g %g %g %g %g %g %g %g %g\n",
                 data.id(), data.pdgcode().string().c_str(), 0, mom.x1(),
                 mom.x2(), mom.x3(), mom.x0(), data.effective_mass(), pos.x1(),
                 pos.x2(), pos.x3(), pos.x0());
  }
}

write()

template<OscarOutputFormat Format, int Contents>

void smash::OscarOutput< Format, Contents >::write ( const Particles &  particles )

inlineprivate

Write the particle information of a list of particles to the output.

One line per particle.

Parameters

[in]

particles

List of particles to be written

Definition at line 119 of file oscaroutput.cc.

                                                                           {
  for (const ParticleData &data : particles) {
    write_particledata(data);
  }
}

create_oscar_output()

std::unique_ptr< OutputInterface > smash::create_oscar_output	(	const std::string &	format,
		const std::string &	content,
		const std::filesystem::path &	path,
		const OutputParameters &	out_par
	)

Returns

A new OscarOutput object using information from config to select the correct implementation.

Parameters

[in]	format	A string: "Oscar2013" or "Oscar1999"
[in]	content	A string: "Particles", "Collisions", "Photons" or "Dileptons".
[in]	path	The path to the output directory where the file(s) will be placed.
[in]	out_par	A structure containing parameters of the output, in particular if it is extended or not, if printing only final particles in event, etc.

Definition at line 811 of file oscaroutput.cc.

                                                                    {
  if (format != "Oscar2013" && format != "Oscar1999") {
    throw std::invalid_argument("Creating Oscar output: unknown format");
  }
  const bool modern_format = (format == "Oscar2013");
  if (content == "Particles") {
    if (out_par.part_only_final == OutputOnlyFinal::Yes) {
      return create_select_format<OscarParticlesAtEventend>(
          modern_format, path, out_par, "particle_lists");
    } else if (out_par.part_only_final == OutputOnlyFinal::IfNotEmpty) {
      return create_select_format<OscarParticlesAtEventendIfNotEmpty>(
          modern_format, path, out_par, "particle_lists");
 
    } else {  // out_par.part_only_final == OutputOnlyFinal::No
      return create_select_format<OscarTimesteps | OscarAtEventstart |
                                  OscarParticlesAtEventend>(
          modern_format, path, out_par, "particle_lists");
    }
  } else if (content == "Collisions") {
    if (out_par.coll_printstartend) {
      return create_select_format<OscarInteractions | OscarAtEventstart |
                                  OscarParticlesAtEventend>(
          modern_format, path, out_par, "full_event_history");
    } else {
      return create_select_format<OscarInteractions>(
          modern_format, path, out_par, "full_event_history");
    }
  } else if (content == "Dileptons") {
    if (modern_format && out_par.dil_extended) {
      return std::make_unique<
          OscarOutput<OscarFormat2013Extended, OscarInteractions>>(path,
                                                                   "Dileptons");
    } else if (modern_format && !out_par.dil_extended) {
      return std::make_unique<OscarOutput<OscarFormat2013, OscarInteractions>>(
          path, "Dileptons");
    } else if (!modern_format && !out_par.dil_extended) {
      return std::make_unique<OscarOutput<OscarFormat1999, OscarInteractions>>(
          path, "Dileptons");
    } else if (!modern_format && out_par.dil_extended) {
      logg[LOutput].warn()
          << "Creating Oscar output: "
          << "There is no extended Oscar1999 (dileptons) format.";
    }
  } else if (content == "Photons") {
    if (modern_format && !out_par.photons_extended) {
      return std::make_unique<OscarOutput<OscarFormat2013, OscarInteractions>>(
          path, "Photons");
    } else if (modern_format && out_par.photons_extended) {
      return std::make_unique<
          OscarOutput<OscarFormat2013Extended, OscarInteractions>>(path,
                                                                   "Photons");
    } else if (!modern_format && !out_par.photons_extended) {
      return std::make_unique<OscarOutput<OscarFormat1999, OscarInteractions>>(
          path, "Photons");
    } else if (!modern_format && out_par.photons_extended) {
      logg[LOutput].warn()
          << "Creating Oscar output: "
          << "There is no extended Oscar1999 (photons) format.";
    }
  } else if (content == "Initial_Conditions") {
    if (modern_format && !out_par.ic_extended) {
      return std::make_unique<
          OscarOutput<OscarFormat2013, OscarParticlesIC | OscarAtEventstart>>(
          path, "SMASH_IC");
    } else if (modern_format && out_par.ic_extended) {
      return std::make_unique<OscarOutput<
          OscarFormat2013Extended, OscarParticlesIC | OscarAtEventstart>>(
          path, "SMASH_IC");
    } else if (!modern_format && !out_par.ic_extended) {
      return std::make_unique<
          OscarOutput<OscarFormat1999, OscarParticlesIC | OscarAtEventstart>>(
          path, "SMASH_IC");
    } else if (!modern_format && out_par.ic_extended) {
      logg[LOutput].warn()
          << "Creating Oscar output: "
          << "There is no extended Oscar1999 (initial conditions) format.";
    }
  }
 
  throw std::invalid_argument("Create_oscar_output got unknown content.");
}

Variable Documentation

current_event_

template<OscarOutputFormat Format, int Contents>

int smash::OscarOutput< Format, Contents >::current_event_ = 0

private

Keep track of event number.

Definition at line 129 of file oscaroutput.h.

file_

template<OscarOutputFormat Format, int Contents>

RenamingFilePtr smash::OscarOutput< Format, Contents >::file_

private

Full filepath of the output file.

Definition at line 132 of file oscaroutput.h.