Skip to content

Latest commit

 

History

History
39 lines (34 loc) · 2.12 KB

README.md

File metadata and controls

39 lines (34 loc) · 2.12 KB
Raw

This is a collection of scripts I am using to analyze my data and plot my results

vHLLE Setup

In order to run a quasi-2D hydro simulation, the eta mapping in the vHLLE config file must be adjusted. For now, I am using the following three parameters besides the default

    nz:          7  (needs to be odd)
etamin:     -0.075
etamax:      0.075

2D Sampler Code

The sampler code for the 2d setup can be found on the branch nsass/sampler_2D of the repository
https://github.com/smash-transport/smash-hadron-sampler/tree/master

The sampler will take the central cell of vHLLE and expand it continuously in eta. To define the sampling range, set eta_min and eta_max in lines 268-269 of smash-hadron-sampler/src/gen.cpp

Creating Rapidity and Spacetime Rapidity Spectra

From the .oscar output of the sampler (find a test output file in ./data/particle_lists.oscar) the dN/dy and dN/dEta spectra can be created with the bash script make_histograms.sh which automatically calls the y_eta_spectra.py script with 5 arguments as

python3 [path_to_script]/y_eta_spectra.py hist_min hist_max num_bins PATH_OSCAR PATH_OUTPUT

where

   hist_min:    lower limit for histogram (make sure that all particles are contained, I use -5.0)
   hist_max:    upper limit for histogram (                                                   5.0)
   num_bins:    number of bins
 PATH_OSCAR:    path to .oscar output from sampler
PATH_OUTPUT:    path to output directory

These parameters can be set in the bash script. Executing the bash script by bash make_histograms.sh will create to files dNdy.txt and dNdEta.txt containing the corresponding histograms.

Plotting the Histograms

Save the python script plot_histograms.py in the same directory as the histogram .txt files and run

python3 plot_histograms.py

This will create the plots for the corresponding histograms. In order to change the particles plotted, one can change the pdg codes inthe plotting script in lines 122 and 125. The available particles can be read off the dictionary in line 42.