PreequilibriumDynamics.h

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file PreequilibriumDynamics.h

/*******************************************************************************
 * Copyright (c) The JETSCAPE Collaboration, 2018
 *
 * Modular, task-based framework for simulating all aspects of heavy-ion collisions
 * 
 * For the list of contributors see AUTHORS.
 *
 * Report issues at https://github.com/JETSCAPE/JETSCAPE/issues
 *
 * or via email to bugs.jetscape@gmail.com
 *
 * Distributed under the GNU General Public License 3.0 (GPLv3 or later).
 * See COPYING for details.
 ******************************************************************************/

#ifndef PREEQUILDYNAMICS_H
#define PREEQUILDYNAMICS_H

#include <vector>
#include "InitialState.h"
#include "JetScapeModuleBase.h"
#include "RealType.h"

namespace Jetscape {
// Flags for preequilibrium dynamics status.
enum PreequilibriumStatus { NOT_STARTED, INIT, DONE, ERR };

class PreEquilibriumParameterFile {
public:
  // preequilibrium dynamics parameters file name.
  char *preequilibrium_input_filename;
};

// Interface for the Preequilibrium Dynamics of the medium
class PreequilibriumDynamics : public JetScapeModuleBase {
private:
  PreEquilibriumParameterFile parameter_list_;
  // record preequilibrium start and end proper time [fm/c]
  real preequilibrium_tau_0_, preequilibrium_tau_max_;

public:
  PreequilibriumDynamics();

  virtual ~PreequilibriumDynamics();

  void Init();

  void Exec();

  virtual void Clear();

  virtual void
  InitializePreequilibrium(PreEquilibriumParameterFile parameter_list) {}
  virtual void EvolvePreequilibrium() {}

  // add initial state shared pointer
  std::shared_ptr<InitialState> ini;

  PreEquilibriumParameterFile &GetParameterList() { return parameter_list_; }

  int GetPreequilibriumStatus() { return (preequilibrium_status_); }

  // @return Start time (or tau) for hydrodynamic evolution
  real GetPreequilibriumStartTime() { return (preequilibrium_tau_0_); }

  // @return End time (or tau) for hydrodynamic evolution.
  real GetPreequilibriumEndTime() { return (preequilibrium_tau_max_); }

  // record preequilibrium running status
  PreequilibriumStatus preequilibrium_status_;

  std::vector<double> e_;
  std::vector<double> P_;
  std::vector<double> utau_;
  std::vector<double> ux_;
  std::vector<double> uy_;
  std::vector<double> ueta_;
  std::vector<double> pi00_;
  std::vector<double> pi01_;
  std::vector<double> pi02_;
  std::vector<double> pi03_;
  std::vector<double> pi11_;
  std::vector<double> pi12_;
  std::vector<double> pi13_;
  std::vector<double> pi22_;
  std::vector<double> pi23_;
  std::vector<double> pi33_;
  std::vector<double> bulk_Pi_;
};

} // end namespace Jetscape

#endif