-
Notifications
You must be signed in to change notification settings - Fork 7
/
Copy pathiSS_parameters.dat
193 lines (165 loc) · 11 KB
/
iSS_parameters.dat
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
hydro_mode = 2 # mode for reading in freeze out information
# 0: reads outputs from VISH2+1
# 1: reads outputs from MUSIC assuming
# boost-invariant
# 2: reads outputs from MUSIC full (3+1)-d
afterburner_type = 1 # 0: PDG_Decay
# 1: UrQMD
# 2: SMASH
turn_on_bulk = 0 # read in bulk viscous pressure
turn_on_rhob = 0 # read in net baryon chemical potential
turn_on_diff = 0 # read in baryon diffusion current
include_deltaf_shear = 1 # include delta f contribution from shear
include_deltaf_bulk = 0 # include delta f contribution from bulk
include_deltaf_diffusion = 0 # include delta f contribution from diffusion
bulk_deltaf_kind = 1 # 0 : 14-moment approximation (s95p-v0-PCE)
# 1 : relaxation time approximation [default]
# 11: OSU 14-moment
# 20: WSU 22-moment (NEoS-BQS) shear and bulk
# 21: WSU Chapman-Enskog (NEoS-BQS) shear and bulk
restrict_deltaf = 1 # flag to apply restriction on the size of delta f
deltaf_max_ratio = 1.0 # the maximum allowed size of delta f w.r.t f0
quantum_statistics = 1 # include quantum statistics or not (1: yes, 0: no)
output_samples_into_files = 0 # output particle samples into individual files
# for individual particle species
store_samples_in_memory = 1 # flag to store particle samples in memory
use_OSCAR_format = 0 # output results in OSCAR format
use_gzip_format = 0 # output results in gzip format (only works with
# store_samples_in_memory = 1)
use_binary_format = 1 # output results in binary format
perform_decays = 0 # flag to perform resonance decay
perform_checks = 0 # perform tests for particle samples
include_spectators = 0 # include spectators (filename: spectators.dat)
local_charge_conservation = 0 # flag to impose local charge conservation
calculate_vn = 0 # 1/0: whether to calculate the
# dN/(pt dpt dphi dy) and v_n flows
# (they not required for MC-sampling)
MC_sampling = 4 # 0/1/2/3: whether to perform Monte-Carlo sampling
# (not required for spectra calculation).
# 0: No sampling.
# 1: use dN_dxtdetady to sample.
# 2: use dN_dxtdy to sample.
# 3: use dN_pTdpTdphidy to sample
# (overwrites calculate_vn to be 1).
# 4: use FSSW (fast)
# Since this parameter controls whether to
# calculate certain arrays, it controls whether to
# perform other related calculations (see below).
dN_dy_sampling_model = 30 # Controls how an non-integer dN_dy is sampled to
# produce an integer that can be used in actual
# sampling
# -- all numbers below 100 are reserved for
# "conventional" sampling where actual particle
# yields are used to constrain the sampling.
# 101-200 are reserved for sampling using total
# energy flux from freeze-out surface.
# -- 1: The fractional part of dN_dy is used as a
# probability to determine whether there is 1
# more particle
# -- 10: Use NBD to sample the fractional part
# of dN_dy. The probability is calculated so
# that p=nb/(nb+k) and k is assumed to be
# proportional to n, and the proportionality
# constant is given by the parameter
# dN_dy_sampling_para1.
# -- 20: Use NBD to sample the whole dN_dy.
# k is assumed to be proportional to n,
# with coefficient dN_dy_sampling_para1.
# -- 30: Use Poisson distribution to sample the
# whole dN_dy. The mean value is set to be
# dN_dy
# -- 110: Total energy flux from freeze-out
# surface is used to constrain dN_dy.
# Whenever the total energy carried away from
# all the sampled particles exceed the total
# energy flux, the sampling procedure stops.
dN_dy_sampling_para1 = 0.16 # Additional parameters for dN/dy sampling.
# -- For dN_dy_sampling_model==10 or 20,
# this parameter sets the ratio k/n for NBD,
# see explanation for dN_dy_sampling_model.
y_LB = -5.0 # lower bound for y-sampling;
# used in "conventional" sampling
y_RB = 5.0 # upper bound for y-sampling; used in
# "conventional" sampling
eta_s_LB = -0.5 # lower bound for eta_s sampling; used only when
# sampling using total energy flux
eta_s_RB = 0.5 # upper bound for eta_s sampling.
use_dynamic_maximum = 0 # 0/1: Whether to automatically reduce the
# guessed maximum after some calculations.
# Work only when MC_sampling is set to 2.
adjust_maximum_after = 100000 # Used only when use_dynamic_maximum=1.
# After the number of sampling given by
# this parameter the guessed maximum is
# adjusted.
adjust_maximum_to = 1.2 # [1,inf]: When guessed maximum is adjusted,
# it is adjusted to the "observed maximum"
# multiplied by this value. Note that the
# "observed maximum" is measured relative to
# the guessed maximum. See code for details.
grouping_particles = 1 # 0/1: Particles will be re-order according to
# their mass. This parameter combined with
# grouping_tolerance parameter can make particles
# with similar mass and chemical potentials to be
# sampled together.
grouping_tolerance = 0.01 # If two particles adjacent in the table have
# mass and chemical potentials close within this
# relative tolerance, they are considered to be
# identical and will be sampled successively
# without regenerating the dN / (dxt deta dy)
# matrix for efficiency.
use_historic_flow_output_format = 0 # 1/0: The "historical flow format"
# means to output flows for all
# particles in a single v2** file and
# to add lines with particles names to
# the v2** file. Turn this option off
# to enbrace the new way of outputting
# which allows the using of parameter
# grouping_particles to speed up the
# calculations.
calculate_vn_to_order = 9 # v_n's are calculated up to this order
sample_upto_desired_particle_number = 0 # flag to run sampling until desired
# particle numbers is reached
number_of_particles_needed = 100000 # number of hadrons to sample
number_of_repeated_sampling = 1000 # How many times should the sampling be
# repeated.
maximum_sampling_events = 10000
sample_pT_up_to = -1 # Up to this value will pT be sampled;
# if<0 then use the largest value in
# the pT table.
sample_y_minus_eta_s_range = 3 # y_minus_eta_s will be sampled between
# +- this value. It is used only when
# sampling using
# sample_using_dN_dxtdy_4all_particles
# function.
use_pos_dN_only = 0 # 1/0: When set to 1, all negative
# emission functions will be skipped.
# Effects the both dN_ptdptdphidy and
# dN_dxtdetady calculations.
minimum_emission_function_val = 1e-30 # If dN/(dx_t deta dy) is evaluated to
# be smaller than this value, then it
# is replaced by this value.
calculate_dN_dtau = 0 # Output dN_dtau table. Only applicable
# if MC_sampling parameter is set to 1.
bin_tau0 = 0.6 # used to generate bins for
# calculate_dN_dtau_using_dN_dxtdeta function
bin_dtau = 0.2 # used to generate bins for
# calculate_dN_dtau_using_dN_dxtdeta function
bin_tau_max = 17.0 # used to generate bins for
# calculate_dN_dtau_using_dN_dxtdeta function
calculate_dN_dx = 0 # Output dN_dx table. Only applicable
# if MC_sampling parameter is set to 1.
bin_x_min = -10.0 # used to generate bins for
# calculate_dN_dx_using_dN_dxtdeta function
bin_dx = 0.5 # used to generate bins
# for calculate_dN_dx_using_dN_dxtdeta function
bin_x_max = 10.0 # used to generate bins for
# calculate_dN_dx_using_dN_dxtdeta function
calculate_dN_dphi = 0 # Output dN_dphi table. Only applicable
# if calculate_vn parameter is set to 1.
calculate_dN_deta = 1 # Output dN_deta table. Only applicable
# if MC_sampling parameter is set to 1.
calculate_dN_dxt = 1 # Output dN_dxt table. Only applicable
# if MC_sampling parameter is set to 1.
output_dN_dxtdy_4all = 0 # Output dN_dxtdy table. Only applicable
# if MC_sampling parameter is set to 2.
randomSeed = -1 # If <0, use system clock.