Activity for MUSIC
-
Anonymous posted a comment on discussion Hydrodynamics - numerics
ohh !!! Thanks for the clarification.
-
Chun Shen modified a comment on discussion Hydrodynamics - numerics
Hi, Thank you for your reporting! Because MUSIC internally uses the mostly plus-sign metric g = diag(-1, 1, 1, 1), there will be an minus sign when computing Grad^\mu(\mu_B/T). The implementation in the current code should be correct. Chun
-
Hi, Thank you for your reporting! Because MUSIC internally use the mostly plus-sign metric g = diag(-1, 1, 1, 1), there will be an minus sign when computing Grad^\mu(\mu_B/T). The implementation in the current code should be correct. Chun
-
Anonymous posted a comment on discussion Hydrodynamics - numerics
In the current version of MUSIC, in the line 823 of dissipative.cpp, there should be a positive sign before NS. But possibly there is a mistake. Kindly correct me if I am wrong. q[mu] and the Navier stokes term q[mu]_NS should be of opposite sign in the RHS of diffusion current evolution equation.
-
Ivan Gaspar posted a comment on discussion General Discussion
Hello dear Chun, I am trying to use BEST EoS in MUSIC. For this purpose I used the tables from the EPOS/src/MSt directory and use the eos_best.cpp code for reading the tables. The thing is that I obtined the following: | __ __ _ _ _ | | U|' \/ '|uU |"|u| | / __"| u ___ U /"___| | | | |\/| |/ | || |<___ \/ |"| | | u | | | | | | | || | u) | | | | |/ | | || |_| <<_/ |_/>> U/| |\u _| | | <<,-,,-. () )( )( (__).-,||,-.// \ | | (./ .) () () _)-' '-(_/()() | ====================================================...
-
Thanks a lot for the reply Chun Cheers, Ivan Gaspar
-
Yes, the events saved in particle_list.dat are from the same hydro surface, referred as oversampled events.
-
Thanks for the info Chun, I was worried about he number of events that I could generate from one single freeze-ou surface, because we would like to compute HBT correlation functions at lower energies. So, I was a little confused. Additionally, when you refer to the oversampled events you mean to the events saved on the particle_list.dat that are generated after iSS , or to the particle_list.dat oversampled files?
-
Thanks for the info Chun, I was worried about he number of events that I could generate from one single freeze-ou surface, because we would like to compute HBT correlation functions at lower energies. So, I was a little confused. Additionally, when you refer to the oversampled events you mean to the events saved on the particle_list.dat that are generated after iSS , or to the particle_list.dat oversampled files?
-
Thanks for the info Chun, I was worried about he number of events that I could generate from one single freeze-ou surface, because we would like to compute HBT correlation functions at lower energies. So, I was a little confused. Additionally, when you refer to the oversampled events you mean to the events saved on the particle_list.dat that are generated after iSS , or to the particle_list.dat oversampled files?
-
Thanks for the info Chun, I was worried about he number of events that I could generate from one single freeze-ou surface, because we would like to compute two-pion correlation function at lower energies. So, I was a little confused. Additionally, when you refer to the oversampled events you mean to the events saved on the particle_list.dat that are generated after iSS , or to the particle_list.dat oversampled files?
-
Dear Ivan, I am not sure about what kind of bias you have in mind. If you want to treat individual hadronic events as independent events, you can set oversampled events to 200-500 typically for the top RHIC and LHC energies. For lower collision energies or small collision systems, you should increase the oversampled events to 1,000 or more. Chun
-
Hello dear developers, I am interested on generate about 5 millions of events implementing iEBE-MUSIC framework. I noticed that this framework have the feature of generating several particle_lists.dat files coming from one freeze-out surface or one hydro dynamical event with MUSIC. My question is directed to how many oversampling files one could obtain from one freeze-out surface, and not get biased. How many oversampling files do you recommend to get from one freeze-out surface?
-
Ankit kumar panda posted a comment on ticket #2
i have set the mentioned set of parameters:
-
You can look at the source code, https://github.com/MUSIC-fluid/MUSIC/blob/15e428efc926abfcace641cf6d2a67ea3b13c85b/src/evolve.cpp#L129 to see which function was used to output the evolution data. The output functions are defined the grid_info.cpp.
-
Which value do you set for the parameter "output_evolution_data"?
-
Grid information on the evolution file
-
Than you Dr. Chen for the support, I am interested in exploring lower energies using iEBE-MUSIC framework. I noticed that many parameters for each stage of the analysis are tuned for RHIC and LHC energies, and some adjustmets are available for 19.6 GeV in the Glauber model. But, I would like to explore even lower energies while using IPGLasma as initial condition model. Could you provide guidance on which parameters should be adjusted to describe lower energies? Additionally, I also noticed that...
-
UrQMD was used as an afterburner instead of a full-fledged event generator. For afterburner mode, UrQMD will directly read in particle lists from particle sampler and run with them. The impact parameter information is not necessary for the UrQMD simualtions. We did not propagate for all the collision information from one module in the chain to the other. In the simualtion chain, the initial state module determines the impact parameter for the collision.
-
Thank you very much for the information. It has helped me to generate initial conditions for my research purposes. I have another question, and I hope you can clarify it for me. After spectra sampling of the freeze out hyper surface I noticed that the impact parameter is missing from the headers in OSCAR.DAT and all subsequent files. However, I do see the impact parameter info in the run.log file generated by IPGlasma analysis. is there a parameter I did not consider or is this the expected behavior?...
-
If you want to generate initial condition for nuclear collisions at certain energy, you should assign the center of mass energy to the roots parameter in the input file. The files under the folder nucleusConfigurations provide event by event nucleus configurations from ab-initio approach from the low-energy nuclear structure community. Certain ab-initio approaches are limited to generate nucleus configurations for light nuclei. If the configurations are not available for some species of nuclei, the...
-
Hello dear colleagues, I am trying to generate initial conditions with IPGlasma, and two doubts arouse. First, I am concerned by previous posts that one can set the collisions energy by changing the parameter "roots" on the IPGlasma input file, is that true? On the other hand, when I try to change the projectile species and targets, such as XeXe, this output appears: [Info] 1238 MB Sampling nucleon positions ... [Info] 1238 MB configuration file for A = 129 is not available, generate the nucleus...
-
iSS performs Monte-Carlo sampling of the particle spectra from the Cooper-Frye formula. It generate particle samples. You can set use_OSCAR_format to 1 in the iSS_parameters.dat file to output the results in the OSCAR format. The explanation of the OSCAR format can be found here: https://theory.gsi.de/~smash/userguide/1.5/format_oscar_particlelist.html
-
iSS perform Monte-Carlo sampling of the particle spectra from the Cooper-Frye formula. It generate particle samples. You can set use_OSCAR_format to 1 in the iSS_parameters.dat file to output the results in the OSCAR format. The explanation of the OSCAR format can be found here: https://theory.gsi.de/~smash/userguide/1.5/format_oscar_particlelist.html
-
Anonymous posted a comment on discussion Cooper-Frye
what are the output files for those means where will the v's and spectra would be stored.
-
Hi, I would suggest to use iSS (https://github.com/chunshen1987/iSS) with the generated hyper-surface to perform the Cooper-Frye particlization. It is numerical more efficient and maintained up-to-date. Chun
-
Anonymous posted a comment on discussion Cooper-Frye
i attach the attachement here .
-
Anonymous posted a comment on discussion Cooper-Frye
Even after using the following input file with assigning mode-1 in the file attached below, my spectra and flows are not calculated only the evolution part is done.
-
I think you are using an obsolete version of MUSIC on sourceforge. The latest version MUSIC can be downloaded from https://github.com/MUSIC-fluid/MUSIC Because of limited manpower, we no longer support optical Glauber model in the repository. You can generate initial condition from any external Glauber model and set MUSIC to read in initial density from file. You can find an example for IP-Glasma initial condition in the current version. If you want to run full-fledge end-to-end simulation for heavy-ion...
-
Anonymous posted a comment on discussion Technical issues with running the code
I have been running the 2+1 D music code with the initial Glauber models, My first issue is: (1) How to change the center of mass energy of the collisions, is it through \sigma_NN (inelastic cross-section). Next if it is so, then i want to see the flow harmonics at different \s_{NN}, so i tried to get the results of flow harmonics and spectra but i have been supplied with the error (2) "[1m[33m[Warning] Can not create folder ./outputs ...[0m" and "[1m[33m[Warning] compute_thermal_particle_spectra_and_vn:...
-
Anonymous posted a comment on discussion General Discussion
Hi, I would like to ask if there is a simple way to use the TRENTo initial conditions instead of the IPglasma? The code is running with the ipglasma conditions, but we prefer trento, and i do not know how it can be done easily. Regards
-
Anonymous posted a comment on discussion Technical issues with running the code
Hi Chun, My apologies, I sent you the wrong paper (it had almost the same name). You helped me out a lot, thank you very much! Lieuwe
-
arXiv:1511.06289 is not my paper. The parameters in arXiv:2005.14682 are the latest ones for IPGlasma + MUSIC + UrQMD simualtions.
-
Anonymous posted a comment on discussion Technical issues with running the code
Hi Chun, The code repositery you sent does not seem to include the article with the arxiv reference "arXiv:1511.06289". Is that right? Lieuwe
-
The model parameters can be found in our code repository, https://github.com/chunshen1987/iEBE-MUSIC/blob/main/config/arXiv_2005.14682/parameters_dict_user_IPGlasma_gamut.py Chun
-
Anonymous posted a comment on discussion Technical issues with running the code
Hi Chun, I did include the resonance decays in iSS. I used "perform_decays = 1" as you can see in "202.par" in my previously send attachments. I have fixed the issue of the MUSIC plots. It is now similar to the iSS results. The difference from ALICE data still persists, but I probably would have to tweak MUSIC or IPGlasma parameters. Since you are one of the authors of the previously linked article, I am curious if the parameter files used for that are available.
-
Hi Lieuwe, Have you included contributions from resonance decays in iSS? Chun
-
Anonymous posted a comment on discussion Technical issues with running the code
Hi Chun, I have now corrected some mistakes on my end such that it agrees with the check_211_spectra.dat. However, when I compare the results with the output of MUSIC, I get completely different yields. The MUSIC data was gained using the output of "music_input_4", where I - as you previously mentioned - set pseudofreeze to 0. The iSS data seems to be a certain factor too low compared to ALICE data, but the shape is similar to the curve mentioned in https://journals.aps.org/prc/abstract/10.1103/PhysRevC.95.064913....
-
I ran your surface file with iSS. The result looks fine. The pion spectra (in check_211_spectra.dat) has a maximum value > 1 and it is reasonable to dN/dy ~ 35 for pion yield.
-
Anonymous posted a comment on discussion Technical issues with running the code
I couldn't add them in one go, but I used the attached music input files to use the music builtin particlisation code, to gain particle yields without using iSS.
-
Anonymous posted a comment on discussion Technical issues with running the code
Hi Chun, thank you for your quick answer. I added the MUSIC input file and the surface file. I ran IPGlasma to create the initial state. I also added my iSS parameters.
-
Can you post your input files and the MUSIC surface file (or a truncated one if it is too large)?
-
Anonymous posted a comment on discussion Technical issues with running the code
Hi Chun, Thanks for you suggestions. I tried using iSS with MC_sampling = 4, but the file "check_211_spectra.dat" does not seem to give a proper yield (maximum is below 1). When I calculate the 211 yield from "particle_samples.bin" it gives a similar result. The shape of the curve however seems to be as expected (if comparing to the Cooper-Frye result directly from MUSIC). When i use MC_sampling = 3, i get a lot of "samples_???.dat" files. But i cannot deduce the format from the code. Could you maybe...
-
Hi Lieuwe, First of all, the Cooper-Fyre part of the code (mode 3) is no longer under active maintaince. It is a lot more numerical efficient if you run particle sampler like iSS with the hypersurface file output from MUSIC mode 2. Your problem is caused by fact the grid range in rapidity or pseudo-rapidity was not big enough for the gsl function to do integeration after translating rapidity to pseudo-rapidity or the other way around. If you set pseudofreeze to 0, then the code will compute particle...
-
Anonymous posted a comment on discussion Technical issues with running the code
Hi Chun, Currently I am running with IP-glasma initial conditions and I get the same error as I got a while ago. 🚫 Error: called out of range rapidity in rap_integrated_flow, -2 < minimum -1.19629 -5 0.01 0.49368 I get a similar error when I had set dNdy_y(eta)(min/max) and dNdyptdpt_y(eta)(min/max) to +/- 2.5 and 5. I've looked in the code at the function Freeze::rapidity_integrated_flow, but I don't understand what is wrong with my settings. Could you maybe explain this error? Before I ran the...
-
Anonymous posted a comment on discussion Technical issues with running the code
Hi Chun, Thank you very much. That seemed to be indeed a problem. Lieuwe
-
Hi Lieuwe, Which version of MUSIC are you using? In your input file, 'Initial_Profile' was set to 3. This option was no longer supported in the latest version. Chun
-
Anonymous posted a comment on discussion Technical issues with running the code
I am currently trying to test run MUSIC, but got the following error: 🚫 Error: called out of range rapidity in rap_integrated_flow, -0.5 < minimum -0.299021 -0.5 0.1 0.13957 Can I get some help with this error? I attached my input file. Kind regards, Lieuwe
-
Anonymous posted a comment on discussion General Discussion
Dear Chun, That worked. Thanks a lot. Anurag,
-
Hi Anurag, You need to set batch_density_output to 0 in the input file. Chun
-
Anonymous posted a comment on discussion General Discussion
Dear Chun, I am trying to run MUSIC, to generate the external hydro profiles needed for JETSCAPE. For that, I need to provide it with some initial energy distribution, for which I am using the 3dMCGlauber repository (v1). I ran it for the UnitTests first and they are running fine. However, when I try to run the main program ./3dMCGlb.e with the following command "./3dMCGlb.e 1" it produces #.dat files (such as ed_etas_distribution_N72.dat, events_summary.dat) which does not have any data in it. I...
-
Hi Monideepa, For the latest version of MUSIC, you can simply run ./compile.sh in the MUSIC root directory to compile the code. The makefile_guillimin is out-of-date and no longer maintained. Chun
-
Anonymous posted a comment on discussion Technical issues with running the code
Hi, I am trying to implement the first example given here https://webhome.phy.duke.edu/~jp401/old_music_manual/music_software.html#running-music-a-simple-example. I am running the code in an Ubuntu system. I have installed GSL as well as MPI libraries, in the music-hydro-code directory. However I am getting the following error. Could you please help? monideepa@monideepa-hp-z240-tower-workstation:~/Desktop$ cd music-hydro-code monideepa@monideepa-hp-z240-tower-workstation:~/Desktop/music-hydro-code$...
-
Thanks for reporting the error! Can you try to use GNU compiler? I will try to improve the cmake options to ensure the code can compile with Apple Clang compiler later. Chun
-
Anonymous posted a comment on discussion Technical issues with running the code
I have downloaded the code and proceeded with the instructions mentioned in https://webhome.phy.duke.edu/~jp401/music_manual/ . I get the following error: command1: cmake .. -- The CXX compiler identification is AppleClang 13.1.6.13160021 -- Detecting CXX compiler ABI info -- Detecting CXX compiler ABI info - done -- Check for working CXX compiler: /Library/Developer/CommandLineTools/usr/bin/c++ - skipped -- Detecting CXX compile features -- Detecting CXX compile features - done -- Found PkgConfig:...
-
Fabian Zhou posted a comment on ticket #1
Is there a difference to the evolve.cpp file from MUSIC-public_stable? I cant see any
-
Hi Egor, We are updating the 3D Glauber model. We will have a refresh of the code packages and make sure all of them are compatible with each other once our paper is released. It should be in one month scale. Regarding to the LHAPDF error, it was from the LHAPDF side. We now use a different approach than the one implemented in the old version. Chun
-
The hyper-surface is output with the following code, https://github.com/MUSIC-fluid/MUSIC/blob/a91e08e933e61169ed6fde1ae7a20233bec8ef35/src/evolve.cpp#L887-L952
-
Egor Nedorezov posted a comment on discussion General Discussion
Dear Chun, I'm trying to use initial conditions from 3D MCGlauber. I installed v1.0 from the 3DMCGlauber git repository as you said early for compatibility with the public version MUSIC and the latest version of MUSIC from the github repository. I launched 3DMCGlauber and got files named "strings_event_XX.dat". I try to use one of them as initial conditions for MUSIC but I got an error: "read_in_QCD_strings_and_partons: the format of filestrings_event_100.datis wrong~". Can you tell me what I'm doing...
-
surface file
-
Hi Fabian, There is no standard formula to compute the average transverse velocity because different weighting function can be chosen in the averaging process. I think the most useful one would be to use the local energy density or entropy density as the weighting function to compute the averaged transverse flow. In the zero net baryon density case, the freeze-out hypersurface is an isothermal surface where local energy density or entropy density is constant. In this situations, many definitions...
-
Hello, I have a question considering the hypersurface. With the surface.dat file we are able to compute for example the final total Energy or the final total Entropy via the formulas: E_f = \int(T^{0\mu} d\Sigma_{\mu}) or S_f = \int(s u^{\mu} d\Sigma_{\mu}), respectively. Since there are also columns in the surface.dat file containing u_x and u_y, is there a formula similar to the above ones, to compute the average transverse velocity on the freeze-out hypersurface? Thanks in advance. Best regards,...
-
Yes, you can use it. There is an alternative parameterization of the lattice EoS from the Duke group (EOS==9) at zero mu_B. In our recent review paper https://arxiv.org/abs/2101.11591, we compared a few of the latest EoSs in the same setting. Chun
-
Anonymous posted a comment on discussion General Discussion
Thanks Chun, Things are fine now. I have one more doubt. Suppose I want to give a hydro run with zero net-baryon density and I want to use a latest EoS instead of s95pv1. Can I use the zero baryon density slice of the EoS developed in the below paper ? A. Monnai, B. Schenke, C. Shen, Phys. Rev. C 100, 024907 (2019).
-
You should set freeze_eps_flag 0 Chun
-
Anonymous posted a comment on discussion General Discussion
Dear Chun, As you instructed, I made the changes in the input file. But following problem is arising. Attaching the input file also. 🎵 2.637 MB For a (2+1)-d simulation, reset deta = 0.1 and eta_size = 0.0 🎵 2.637 MB Finished checking input parameter list. Everything looks reasonable so far 😁👍🍺🍺🍻🍻 🎵 2.637 MB Using lattice EOS at finite muB from A. Monnai 🎵 2.637 MB reading EOS neos_b ... 🎵 2.637 MB from path ./EOS/neos_b/ 🎵 9.527 MB Done reading EOS. 🎵 10.2 MB OpenMP: using 80 threads. 🎵 10.2 MB...
-
Dear Tribhuban, I think the error is caused becuase you tried to freeze-out at a constant temperature. At finite baryon density, the local temperature is a function of energy density and net baryon density. The current program supports to freeze-out at a constant energy density in the finite net baryon density case. Can you try to use a constant energy density to freeze-out? For T_frez = 150 MeV at n_B = 0, the local energy density is around e = 0.25 GeV/fm^3. Please set the following parameters...
-
Anonymous posted a comment on discussion General Discussion
Dear Authors, I tried to give a hydro run with the latest EoS , neos_b . As I know this EoS will be activated in MUSIC when eos_id==12. The EoS/neos_b directory was missing. I downloaded it from the site https://sites.google.com/view/qcdneos/ . I intend to use the zero baryon baryon density part of this EoS so I set the baryon density to be zero from the beginning. The hydro run started from tau=0.4 and the maximum temperature was 429MeV initially. But the problem is, all the cells got frozen out...
-
Great!
-
Anonymous posted a comment on discussion Cooper-Frye
Hi chun, I found out where the problem was. Actually the generated hypersurface file name was "surface_0.24.dat" where "0.24" was the energy density at the freezeout hypersurface. I renamed the file as "surface.dat" and things are fine now. Thanks Tribhuban
-
Hi Tribhuban, I haven't run the Cooper-Frye part of calculations for many years. Can you post the hydrodynamic hypersurface so that I can ran a quick check on mode 3, which generate the file "yptphiSpectra.dat"? Thanks! Chun
-
Anonymous posted a comment on discussion Cooper-Frye
sorry, one mistake in above information. I had set "s_factor" to be 1.0 in the input file.
-
Anonymous posted a comment on discussion Cooper-Frye
Hi authors, I have taken an external input file as an initial condition. The file is attached below. The initial condition I have taken belongs to 20-30% centrality of Au+Au collision. It is a two component optical glauber initial entropy distribution with hardness factor 0.14. In the input file the entropy density is scaled, which will approximately reproduce the required final multiplicity . So I have taken "s_factor" to be zero in the MUSIC input file. The input file can be found in the attachment...
-
Anonymous posted a comment on discussion General Discussion
Thanks Chun for your reply.
-
Hi, The WB eos parameterization was taken from the hydrodynamic code GPU-VH developed by Dennis Bazow at Ohio State University. You can refer to their paper for some details, https://arxiv.org/pdf/1608.06577.pdf I haven't used this equation of state for any production simulations myself.
-
Anonymous posted a comment on discussion General Discussion
Dear authors, 'eos_wb.cpp' file contains a parameterised form of lattice EoS. In this case, does the Wuppertal-Budapest collaboration lattice result is matched with HRG EoS in low temperature region ? Can this EoS be used for very low T (say below 100 MeV ) ? From where I can read more about this parametrisation ? Thanks
-
Hi Martin, The public_stable branch is more current than v3.0 as there are some updates since the v3.0 release. We will gradually add v3.x release to incorporate recent updates in the public_stable branch. The 3dMCGlauber model is under active developement. I change the inferface very frequently. The connection to the IPGlasma model should be stable. Can you provide more information about how you run simulations? Since IP-Glasma, 3dMCGlauber, and MUSIC are individual modules for simulating heavy-ion...
-
Anonymous posted a comment on discussion General Discussion
Hello everyone, thanks for making the wonderful program available. I am still a little confused about which version to use. There are the versions with tags v2.5 and v3.0 available, which would be reasonable. However, the public_stable branch seems more current. The version with the current updates seems to belong to v2.5 however: v2.5-9-g7be8d74. I have noticed for example that in v3 if I activate freeze-out by temperature I get a message that this is not supported yet. On the other hand it seems...
-
Dario Chaires posted a comment on discussion Technical issues with running the code
Yes, I'm using the public_stable branch. Here goes the surface file, also I have another note: even though I've set the thime evolution to Tau = 30 fm, the outpust message reads that all cells are frozen out, leaving the final time evolution to 6 fm, I'm not sure what parameters to use in order to get the full-evolution time, assuming that I'm using the temprerature T as the criteria to stop hydrodynamic evolution. Thanks, Chun! Dario
-
Hi Dario, Can you post the surface file from mode 2? I would like to try to reproduce your error. And you are running the latest version of MUSIC on the Github repository (the public_stable branch), right? Thanks! Chun
-
Hello, I have several days working with the latest version of MUSIC, so far just mode 2 and 3. I have two questions about the outputs: 1) From mode 2, I manage to run an example with ideal hydeodynamica using the IP -Glasma initial conditions of Au-Au at 200 GeV. All the outputs look fine, although I'm not sure how the surface data should be readed, I made a graph by rougly plotting all the elements of every vector coming from the surface file to give me an idea of the output. Looking at the data...
-
Hello, I have several days working with the latest version of MUSIC, so far just mode 2 and 3. I have two questions about the outputs: 1) From mode 2, I manage to run an example with ideal hydeodynamica using the IP -Glasma initial conditions of Au-Au at 200 GeV. All the outputs look fine, although I'm not sure how the surface data should be readed, I made a graph by rougly plotting all the elements of every vector coming from the surface file to give me an idea of the output. Looking at the data...
-
Hello, I have several days working with the latest version of MUSIC, so far just mode 2 and 3. I have two questions about the outputs: 1) From mode 2, I manage to run an example with ideal hydeodynamica using the IP -Glasma initial conditions of Au-Au at 200 GeV. All the outputs look fine, although I'm not sure how the surface data should be readed, I made a graph by rougly plotting all the elements of every vector coming from the surface file to give me an idea of the output. Looking at the data...
-
Hello, I have several days working with the latest version of MUSIC, so far just mode 2 and 3. I have two questions about the outputs: 1) From mode 2, I manage to run an example with ideal hydeodynamica using the IP -Glasma initial conditions of Au-Au at 200 GeV. All the outputs look fine, although I'm not sure how the surface data should be readed, I made a graph by rougly plotting all the elements of every vector coming from the surface file to give me an idea of the output. Looking at the data...
-
Hi, You can use "freeze_out_method 4" (Cornelius) to handle complicate freeze-out geometry. Chun
-
Anonymous posted a comment on discussion General Discussion
I gave a run in MUSIC for 3+1D evolution of a optical glauber initial condition. The input file that I am using is attached below. In runtime I am getting following kind of warning message in terminal. [.....] [Warning] NUMBER OF SINGLE EDGES IS NOT A MULTIPLE OF 3, number=5 percent error=8.27606 [Warning] NUMBER OF SINGLE EDGES IS NOT A MULTIPLE OF 3, number=4 percent error=7.32281 [Warning] NSurfaces=12 but number of tetrahedra=18 tauf=1.07005 xf=6.94539 yf=-0.64898 etaf=4.09063 single edges=4...
-
The NEoS family can be downloaded from https://sites.google.com/view/qcdneos/home The NEoS is a family of Equation of States at finite baryon density. If you just want to use EoS at zero mu_B, you can set EOS_to_use to 9 in the parameter file.
-
The NEoS family can be downloaded from https://sites.google.com/view/qcdneos/home The NEoS is a family of Equation of States at finite baryon density. If you want to just use EoS at zero mu_B, you can set EOS_to_use to 9 in the parameter file.
-
Yes, you should use the github repository for the latest version of MUSIC.
-
Indranil Das posted a comment on discussion General Discussion
Hi MUSIC team, I guess I found the answer by following the discussion at https://sourceforge.net/p/music-hydro/discussion/issues/thread/efde9a81df/, that it is now encouraged to use github as official version. Am I correct ? As mentioned in the post above the github still has problem with input files. It will be very helpful if working example input files are uploaded in the git. regards Indra
-
Dear Ana, The required combination (IP-Glasma + MUSICS + URQMD) with other possible combinations is now available as iEBE-MUSIC, check at, https://github.com/chunshen1987 regards Indra
-
Dear MUSIC team, First a big thank you for building this 3+1D heavy-ion hydro-MC and keeping this in public domain. One can download the source codes either from github as mentioned at https://www.physics.mcgill.ca/music/#Download or the package available in this site (sourceforge) as mentioned in the manual https://webhome.phy.duke.edu/~jp401/music_manual/music_software.html Which one is the Official and Stable version ? The reason for this question is, If I judge by commit date then github is most...
-
Wei-Liang Qian posted a comment on discussion Technical issues with running the code
Hi, Chun, Thanks a lot for all the helps. W.L.
-
Hi Wei-Liang, (1) Not all the Initial conditions provide the grid information in the file in the same format. It is hard to write a generic function to deal with this. Right now, MUISC support to override the grid information for the IP-Glasma initial conditions. Please see examples in the init.cpp (Initial_profile == 8 & 9). (2) You can write a script to loop over events outside MUSIC. (3) Cooper Frye and resonance decays part of code were developed before hydrodynamics + hadronic transport approach....
-
Hi, Chun, After some trial and error, also after getting a lot of helps from Matt, I managed to run the version of MUSIC 3.0 recently released. I summarize below some minor tuning I did in order to run the code, and also post some further doubts and questions. As a first test, I am trying to run a couple of events by feeding the code a NEXUS/EPOS IC. I still have a few questions, hopefully it does not too much of your time, really appreciated! Minor issues I encountered. (1) It seems that one has...
-
Hi Chun, According the above instructions, as well as the reading of some other threads in the fourm. I tried to run one event of IPGlasma IC for Au+Au 200GeV, downloaded from Bjoern's page. Also by using the following options mode 1 Initial_profile 92 Initial_Distribution_input_filename initial/epsilon-u-Hydro.dat boost_invariant 1 EOS_to_use 9 where I guess "boost_invariant 1" for used for 2D+1 caluclations. The hydrodynamical evolution seems work, since the "eccentricities_evo_eta-0.50.5.dat"...
-
Hi Chun, According the above instructions, as well as the reading of some other threads in the fourm. I tried to run one event of IPGlasma IC for Au+Au 200GeV, downloaded from Bjoern's page. Also by using the following options mode 1 Initial_profile 92 Initial_Distribution_input_filename initial/epsilon-u-Hydro.dat boost_invariant 1 EOS_to_use 9 where I guess "boost_invariant 1" for used for 2D+1 caluclations. The hydrodynamical evolution seems work, since the "eccentricities_evo_eta-0.50.5.dat"...
-
Hi Chun, In both tests, I used the music_input_2, generated by the script generate_music_inputfile.py Its content is attached. Then, I realized it seems that the folder "Initial_TA_Distribution_Filename initial" does not exist, so I need to generate/download them myself? Many thanks! Best Regards Wei-Liang echo_level 1 mode 2 eta_rhob_plateau_height 0.3 Eta_plateau_size 2.7 initial_eta_rhob_profile 2 eta_rhob_width 1.0 Eta_fall_off 0.3 Initial_rhob_TB_Distribution_Filename initial/rhob_TB_event_1_block.dat...
-
Hi Chun, More question. I am trying to follow you suggestion of using enviroment variable $HYDROPROGRAMPATH by setting it to the MUSIC root folder. Then try to run the code where a third folder using $HYDROPROGRAMPATH/MUSIChydro music_input_mode_2 where I have copied music_input_mode_2 to the spot. I got the following error message. (The same command seems to work when executed in the root folder of MUSIC) Is there a way to work around? Many thanks! [wlqian@admin b20200629.music.tester.tube]$$$HYDROPROGRAMPATH/MUSIChydro...
-
Hi Chun, Just to confirm that after copying the file containing the information on the initial conditions, it does work. I should have look more carefully at the option file. Thanks! -- Best Regards Wei-Liang
1 >
