Causal viscous hydro code for central heavy-ion collisions
by Paul Romatschke, 1/2007

Preprints: 
nucl-th/0610108
nucl-th-0701032


Overview

*Information about the actual hydro C code can be found
in /viscous_code/readme.txt


------------------------------------------------------------
* Freeze-out is done in Mathematica (for historical reasons)
It's done a bit from behind: first you have to copy
the contents of the decays directory to the data/ directory
where the freezeout.dat file resides
In the freeze-out subdirectory, you then 
have to change the "header" declaration in 

getspectra_new3.m

such that it points to your directory and in the

dolist.m 

file, change the Tf to the correct freeze-out temperature. The
"oheader" variable in dolist.m is optional, but be sure to
set it to oheader="" if you don't want to use it
Now, a simple command 

nohup math < dolist.m > out.txt &

should read in all the mass resonances specified in the data file
(see below), calculate their spectra for the freeze-out surface
that you specified and put out the whole stuff into the file

phipspectra.dat

in the somewhat funny AZHYDRO format

I'm sorry that this is so complicated!
-----------------------------------------------------------------

* Resonance decays are calculated using the AZHYDRO decay routines
The slightly modified routines can be found in the 
decay_source subdirectory and -- once compiled -- give the executeable
reso that would replace the copy in the decays/ subdirectory

To calculate the decays, make a subdirectory results/ in the directory
where the file phipspectra.dat is located. The file
pdg_weak.dat controls which resonances are to be included,
and pasim.dat,pasin.dat and pasinames.dat are derivatives of this master
file that are required by the resonance decay code. See the
AZHYDRO code or the source in decay_source for details.
Issuing ./reso should start the resonance decay code, which results
in files like spec_211.dat and the like, which correspond to the 
inclusive spectra. To obtain plottable tables, modify the header variable
in "extract_spec.m" file and run 

math < extract.m 

This should produce the files ps.dat, Ks.dat and pis.dat in the results/
subdirectory that you created earlier, which are the proton, kaon and pion
spectra with resonance-decays included


-------------------------------------------------------------------------

*HBT Radii are also calculated in Mathematica, see the 
Mathematica demo.nb notebook in the HBT subdirectory which hopefully should
be self-explanatory

--------------------------------------------------------------------------

*For convenience, I also include results I calculated in the 
results/ subdirectory. There are results for ideal hydro 
(actually very small eta/s), with directory names like "nwE0.0001".
Then there are the viscous hydro results, which are characterized
by the value of the relaxation time used (directory names "C0.75"
correspond to \tau_\Pi=\eta/s*1.5/T, while directory names "C3.0"
correspond to \tau_\Pi=\eta/s*6/T) Following this number, the
directory names include E0.08, E0.16 and so on which specifies the
ratio \eta/s that was used. The next number T0.34, T0.32 specifies
the initial central temperature, T0 and finally Tf0.18, Tf0.17 and so
on specify the freeze-out temperature Tf, both in GeV.
To double-check, see the params.txt included in the respective directories.