vhlle/src/icGlauber.h

0001
0002 class EoS;
0003 class Fluid;
0004 class TF1;
0005
0006 // this class takes care of the initial conditions for hydrodynamic evolution
0007 class ICGlauber {
0008   TF1 *iff;  // instance of TF1 class used for numerical integration relevant to
0009              // the initial conditions from optical Glauber approach
0010   double rho0, prms[4];  // relevant to the initial conditions from optical
0011                          // Glauber approach
0012   double *_rphi;
0013   void findRPhi(void);
0014   double rPhi(double phi);
0015
0016   double WoodSaxon(double *x, double *p);  // Wood-Saxon profile
0017   double Thickness(double *x, double *p);  // nuclear thickness profile function
0018                                            // for optical Glauber approach
0019   // epsilon: normalization of the initial energy density
0020   // alpha: parameter relevant to the initial transverse flow
0021   // b: impact parameter (for optical Glauber)
0022   double epsilon, b, tau0;
0023
0024  public:
0025   ICGlauber(double e, double impactPar, double _tau0);
0026   ~ICGlauber(void);
0027   // energy density profile at given point in transverse plane
0028   double eProfile(double x, double y);
0029   // setIC: initializes entire hydro grid at a given initial proper time tau
0030   void setIC(Fluid *f, EoS *eos);
0031 };