Output Classes

Classes
class	smash::BinaryOutputBase
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::ICOutput
class	smash::OscarOutput< Format, Contents >
class	smash::OutputInterface
	Abstraction of generic output. More...
class	smash::RootOutput
class	smash::ThermodynamicOutput
	Writes the thermodynamic quantities at a specified point versus time. More...
class	smash::VtkOutput

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 bf::path &path, const std::string &name)
	Create oscar output. More...
void	smash::OscarOutput< Format, Contents >::at_eventstart (const Particles &particles, const int event_number) 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, double impact_parameter, bool empty_event) 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) 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 bf::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 bf::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  )

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.

Implements smash::OutputInterface.

Definition at line 128 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 + 1,
                   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 + 1);
    }
    if (!(Contents & OscarParticlesIC)) {
      // We do not want the inital particle list to be printed in case of IC
      // output
      write(particles);
    }
  }
}

Here is the call graph for this function:

at_eventend()

template<OscarOutputFormat Format, int Contents>

void smash::OscarOutput< Format, Contents >::at_eventend ( const Particles &  particles, const int  event_number, double  impact_parameter, bool  empty_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]	impact_parameter	Impact parameter of this event.
[in]	empty_event	Whether there was no interaction between the target and the projectile.

Implements smash::OutputInterface.

Definition at line 152 of file oscaroutput.cc.

                                                                  {
  if (Format == OscarFormat2013 || Format == OscarFormat2013Extended) {
    if (Contents & OscarParticlesAtEventend ||
        (Contents & OscarParticlesAtEventendIfNotEmpty && !empty_event)) {
      std::fprintf(file_.get(), "# event %i out %zu\n", event_number + 1,
                   particles.size());
      write(particles);
    }
    // Comment end of an event
    const char *empty_event_str = empty_event ? "yes" : "no";
    std::fprintf(file_.get(), "# event %i end 0 impact %7.3f empty %s\n",
                 event_number + 1, 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 && !empty_event)) {
      std::fprintf(file_.get(), "%zu %zu %i\n", particles.size(), zero,
                   event_number + 1);
      write(particles);
    }
    // Null interaction marks the end of an event
    std::fprintf(file_.get(), "%zu %zu %i %7.3f\n", zero, zero,
                 event_number + 1, 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) {
      logg[LHyperSurfaceCrossing].warn(
          "End time might be too small for initial conditions output. "
          "Hypersurface has not yet been crossed by ",
          particles.size(), " particle(s).");
    }
  }
}

Here is the call graph for this function:

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);
    }
  }
}

Here is the call graph for this function:

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  )

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.

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_ + 1,
                   particles.size());
    } else {
      const size_t zero = 0;
      std::fprintf(file_.get(), "%zu %zu %i\n", particles.size(), zero,
                   current_event_ + 1);
    }
    write(particles);
  }
}

Here is the call graph for this function:

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 704 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\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());
  } 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());
  }
}

Here is the call graph for this function:

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 121 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 bf::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 768 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 make_unique<
          OscarOutput<OscarFormat2013Extended, OscarInteractions>>(path,
                                                                   "Dileptons");
    } else if (modern_format && !out_par.dil_extended) {
      return make_unique<OscarOutput<OscarFormat2013, OscarInteractions>>(
          path, "Dileptons");
    } else if (!modern_format && !out_par.dil_extended) {
      return 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 make_unique<OscarOutput<OscarFormat2013, OscarInteractions>>(
          path, "Photons");
    } else if (modern_format && out_par.photons_extended) {
      return make_unique<
          OscarOutput<OscarFormat2013Extended, OscarInteractions>>(path,
                                                                   "Photons");
    } else if (!modern_format && !out_par.photons_extended) {
      return 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 make_unique<
          OscarOutput<OscarFormat2013, OscarParticlesIC | OscarAtEventstart>>(
          path, "SMASH_IC");
    } else if (modern_format && out_par.ic_extended) {
      return make_unique<OscarOutput<OscarFormat2013Extended,
                                     OscarParticlesIC | OscarAtEventstart>>(
          path, "SMASH_IC");
    } else if (!modern_format && !out_par.ic_extended) {
      return 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.");
}

Here is the call graph for this function:

Here is the caller graph for this function:

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 128 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 131 of file oscaroutput.h.