Running with IP-Glasma initial conditions

[Anonymous]

Hi!

I am trying to run MUSIC with publically available IP-Glasma input files.
https://quark.phy.bnl.gov/~bschenke/ipglasma-doc-rhic.txt

I am getting
[Warning] u0 varies more than 100 times compared to its value at previous time step
[Warning] e = 5.06842, u[0] = 1.00129, prev_u[0] = 0

Which input parameters should I use?
Here are some of essential changes I did to example input file in music_input_mode_2

< Initial_profile 8 # type of initial condition

Initial_profile 3 # type of initial condition
< Initial_Distribution_input_filename ./inputs/gs/project/cqn-654-ab/jfpaquet/init_cond/rhic_jul2014/RHIC00-05-t04/1/epsilon-u-Hydro.dat
Initial_Distribution_input_filename example_inputs/filename
< s_factor 1.0 # normalization factor read in
s_factor 30.0 # normalization factor read in
< Total_evolution_time_tau 20. # the maximum allowed running
Total_evolution_time_tau 30. # the maximum allowed running
< Delta_Tau 0.01 # time step to use in the evolution [fm/c]
Delta_Tau 0.02 # time step to use in the evolution [fm/c]
< X_grid_size_in_fm 34.0 # spatial range along x direction in the
X_grid_size_in_fm 20.0 # spatial range along x direction in the
< Y_grid_size_in_fm 34.0 # spatial range along y direction in the
Y_grid_size_in_fm 20.0 # spatial range along y direction in the

Thanks!
Aleksas

[Chun Shen]

Hi Aleksas,

It seems that there is a bug in the warning message. prev_u[0] should be always larger than 1. Can you attach the initial condition file and music_input_mode_2 in the reply? I will take a look at why prev_u[0] = 0. Thanks a lot!

Chun

[Anonymous]

Hi!

Maybe the file is not read in properly? I didn't check that, just assumed that it is compatible with public IP-Glasma imput files.

Aleksas

[J-F Paquet]

Hi Aleksas,

"Initial_profile 8" is good, this reads in an external file

"Initial_Distribution_input_filename ./inputs/gs/project/cqn-654-ab/jfpaquet/init_cond/rhic_jul2014/RHIC00-05-t04/1/epsilon-u-Hydro.dat" may not work. Can you put the full directory of the initial condition file, instead of the relative directory?

The grid size doesn't need to be specified, it's read from the initial condition files (whatever is in the input file is overwritten).

The public version of MUSIC is currently set-up so that it read the initial condition files as an 11 columns file, with only four columns used:
" profile >> not_used >> not_used >> not_used
>> energy_density >> u^\tau >> u^x >> u^y
>> not_used >> not_used >> not_used >> not_used;"

I suspect that the RHIC IP-Glasma initial conditions are only 7 columns wide, not 11. Could you check? If that's the case, edit line 1142-1143 of init.cpp to make MUSIC read the file properly.

I think the LHC IP-Glasma might have the expect 11 columns width (there was a small change of format between RHIC and LHC)

I put additional documentation on the wiki: https://sourceforge.net/p/music-hydro/wiki/Initial%20conditions/ . We'll try to add it to the code internal documentation.

Let us know if it works.

Cheers,

J-F

[Anonymous]

Thanks for the reply. Yes, RHIC events are 7 collumn. It seems to work now.
I suggest putting in the read in format in an example imput file.

Aleksas

[Chun Shen]

Thanks! I will add some codes to check the format of the IP-Glasma initial condition in the code to avoid similar problem in the future.