FluidCellInfo.h

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file FluidCellInfo.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.
 ******************************************************************************/
// This is a general basic class for hydrodynamics

#ifndef FLUIDCELLINFO_H
#define FLUIDCELLINFO_H

#include "RealType.h"

namespace Jetscape {

class FluidCellInfo {
public:
  // data structure for outputing fluid cell information
  Jetscape::real energy_density;  
  Jetscape::real entropy_density; 
  Jetscape::real temperature;     
  Jetscape::real pressure;        
  Jetscape::real
      qgp_fraction; 
  Jetscape::real mu_B;       
  Jetscape::real mu_C;       
  Jetscape::real mu_S;       
  Jetscape::real vx, vy, vz; 
  Jetscape::real pi[4][4];   
  Jetscape::real bulk_Pi;    

  FluidCellInfo();

  FluidCellInfo inline operator*=(Jetscape::real b);

  //void Print();
};

//overload +-*/ for easier linear interpolation
inline FluidCellInfo operator+(FluidCellInfo a, const FluidCellInfo &b) {
  a.energy_density += b.energy_density;
  a.entropy_density += b.entropy_density;
  a.temperature += b.temperature;
  a.pressure += b.pressure;
  a.qgp_fraction += b.qgp_fraction;
  a.mu_B += b.mu_B;
  a.mu_C += b.mu_C;
  a.mu_S += b.mu_S;
  a.vx += b.vx;
  a.vy += b.vy;
  a.vz += b.vz;
  for (int i = 0; i < 4; i++) {
    for (int j = 0; j < 4; j++) {
      a.pi[i][j] += b.pi[i][j];
    }
  }
  a.bulk_Pi += b.bulk_Pi;
  return a;
}

// Multiply the fluid cell with a scalar factor
FluidCellInfo inline FluidCellInfo::operator*=(Jetscape::real b) {
  this->energy_density *= b;
  this->entropy_density *= b;
  this->temperature *= b;
  this->pressure *= b;
  this->qgp_fraction *= b;
  this->mu_B *= b;
  this->mu_C *= b;
  this->mu_S *= b;
  this->vx *= b;
  this->vy *= b;
  this->vz *= b;
  for (int i = 0; i < 4; i++) {
    for (int j = 0; j < 4; j++) {
      this->pi[i][j] *= b;
    }
  }
  this->bulk_Pi *= b;
  return *this;
}

inline FluidCellInfo operator*(Jetscape::real a, FluidCellInfo b) {
  b *= a;
  return b;
}

inline FluidCellInfo operator*(FluidCellInfo a, Jetscape::real b) {
  a *= b;
  return a;
}

inline FluidCellInfo operator/(FluidCellInfo a, Jetscape::real b) {
  a *= 1.0 / b;
  return a;
}

// print the fluid cell information for debuging
// this function has bugs
//std::ostream &operator<<(std::ostream &os, const FluidCellInfo &cell) {
//    os << "energy_density=" << cell.energy_density << std::endl;
//    os << "entropy_density=" << cell.entropy_density << std::endl;
//    os << "temperature=" << cell.temperature << std::endl;
//    os << "pressure=" << cell.pressure << std::endl;
//    os << "qgp_fraction=" << cell.qgp_fraction << std::endl;
//    os << "mu_B=" << cell.mu_B << std::endl;
//    os << "mu_C=" << cell.mu_C << std::endl;
//    os << "mu_S=" << cell.mu_S << std::endl;
//    os << "vx=" << cell.vx << std::endl;
//    os << "vy=" << cell.vy << std::endl;
//    os << "vz=" << cell.vz << std::endl;
//    os << "pi[mu][nu]=" << std::endl;
//    for (int i = 0; i < 4; i++) {
//        for (int j = 0; j < 4; j++) {
//            os << cell.pi[i][j] << ' ';
//        }
//        os << std::endl;
//    }
//    os << "bulk_Pi=" << cell.bulk_Pi;
//    return os << std::endl;
//}

} // end namespace Jetscape

#endif // FluidCellInfo