AMPT Users' Guide
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
12/2018 test version v1.26t9b/v2.26t9b:
* Fixed bugs that can cause segmentation fault (especially for default
AMPT at high energies):
exclude endpoints of 0. and 1. in random values from RANART()
in amptsub.f in order to avoid crash in case of 0 branching ratio,
resample random value from RLU(0) when TEVB in PYSSPA does not change
(this may happen when RLU(0)>~0.9999 and thus
ex==|LOG(RLU(0))*B0/(5.*WTSUM)|<~3e-8 and thus EXP(ex)=1.00000000,
which may cause IS(2) to be undefined and be >10^8 (out of bound)).
The above modification can be found by searching "clin-12/2018".
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
9/2018 test version v1.26t9/v2.26t9:
* Freeze may happen to a run at low energies (e.g. at ~10A-20A GeV)
when an event goes into an infinite loop in HIJHRD/PYTHIA/PYREMN;
fixed now, and a message will be written to nohup.out.
* Also included 12/2016 changes for test version v1.26t8/v2.26t8:
HIJING occasionally produces a particle with negative energy, which
leads to no hadron record for that event in ampt.dat (e.g. pPb at LHC);
made correction in this rare case by recalculating the particle energy;
print message to nohup.out in case of unphysical particle
energy/momentum values;
reset FT(5) to |GZ5(I)| if the former is smaller in order to avoid
Floating-Point Exception in the eta variable.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
5/2018 test version v1.26t7b/v2.26t7b:
* Fixed a bug that causes fatal HIJHRD or HIJSFT error and leads to
an incorrect impact parameter distribution (at sqrt(sNN)>=5.023 TeV)
by resetting the energy in case energy < |pz| (due to finite precision
at very high energies for particles with energy ~= beam energy).
* Updated some explanations in input.ampt.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
5/2016 test version v1.26t7/v2.26t7:
* (w/ GL Ma) Changed the following line in art1f.f to avoid the bug
that may lead to loss of hadrons inside ART (the bug is due to
finite accuracy that results in (ct-dt) + dt != ct exactly):
& FT1(NP1, IRUN) .GT. (CT - DT) .AND.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
8-9/2015 test version v1.26t6/v2.26t6:
* Put in safe conditions to avoid Floating-Point Exception in amptsub.f.
* Fixed a bug that may cause the first few events to be missing in
ampt.dat (mostly for low-multiplicity events such as PP collisions;
this error would also be indicated by file fort.99 that writes
"Unexpected perturbative particles").
* Changed ptwo(2,5) and related variables to double precision
to avoid IEEE_DIVIDE_BY_ZERO or IEEE_INVALID or IEEE_OVERFLOW_FLAG
(especially at high energy such as LHC).
* For String Melting events, added lines to record hadron information
in case it will be used for the constituent partons of the hadron.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4/2015 test version v1.26t5/v2.26t5:
* Random seed for HIJING is modified in main.f, so that a different
random seed will always lead to a different random number sequence
(in earlier versions, an even integer leads to the same random number
sequence as the odd integer that is bigger than it by 1).
* set upper limit for gamma_s=P(s)/P(u) to 0.4 (to limit strangeness
enhancement when string tension is strongly increased due to
a very low value of parameter b in Lund symmetric splitting function
as done in arXiv:1403.6321).
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
8/2014 test version v1.26t4/v2.26t4:
* Crash can come from finite precision of the defined asinh(x),
where the argument of LOG() can be 0 for large negative x;
defined asinh(x) differently in amptsub.f.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2/05/2013 test version v1.26t3/v2.26t3:
* Modified amptsub.f so that spectator target nucleons in the LAB frame
will have the correct rapidity value of 0.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
12/19/2012 test version v1.26t2/v2.26t2:
* Modified format of ana/npart-xy.dat:
added last 2 columns on the present and original flavor code of nucleon.
* Corrected original nucleon flavor in the calculation of NN differential
cross sections in HIJING.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
9/2012 test version v1.26t1/v2.26t1:
* Use double precision in art1f.f when calculating sqrt(s) to avoid crash
(segmentation fault due to s<0) at high energy such as LHC at large NT;
also use double precision in linana.f to calculate beta to avoid NaN.
* Use a more general formula in amptsub.f to calculate rapidity
(valid for hadrons at large rapidities that have |Pz|>=E
due to finite precision)
5/2012
* Removed " $(LIBS)" from Makefile.
* Fixed data-type mismatches (mostly in hipyset1.35.f and linana.f).
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4/2012 test version v1.25t7d/v2.25t7d:
* Added pi0 electromagnetic decay flag (ipi0dcy) in file `input.ampt';
set to 1 to turn on pi0 decays after the hadron cascade.
* Modified format of ana/npart-xy.dat:
write impact parameter at the end of the first line of each event;
events in npart-xy.dat correspond to those in ana/ampt.dat (in that
now events that eventually failed fragmentation and thus get repeated
are not written to ana/npart-xy.dat).
* changed format of zpc.dat: now for both default and String Melting runs,
each line after the 1st line of an event gives the same following info:
PYTHIA particle ID number, three-momentum(Px,Py,Pz), mass, and
space-time coordinates(x,y,z,t) of one final parton at freeze-out.
* corrected Kshort weak-decay width (used when Kshort decay flag is on).
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4/2012 test version v1.25t7c/v2.25t7c:
* fixed a bug (made in test version v1.25t7b/v2.25t7b),
which incorrectly sets the reaction plane at 5Pi/4.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2/2012 test version v1.25t7b/v2.25t7b:
* modified code to enable random orientation of reaction plane by setting
iphirp in input.ampt to 1. The value of the angle (0 to 2pi) is written
out as the last value in the first line of each event in ana/ampt.dat.
* fixed a bug related to event triggering using pttrig and maxmiss.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
9/2011 test version v1.25t7/v2.25t7:
* added Pt tolerance etc to determining spectator nucleons consistently,
this is found to only affect string melting runs when LAB frame is used
(in that 1: npart values in ampt.dat are now correct and
2: the run will not generate incorrect events that have no particles)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
7/2011 test version v1.25t6/v2.25t6:
* Turned off heavy hadron decays in HIJING, so that the net-HeavyQuark
number will be 0 in produced hadrons from HIJING (e.g. #c=#cbar).
* Added a missing /para7/ in hijing1.383_ampt.f for OSCAR output,
this adds 1 output line to ana/parton-after-coalescence.dat for
String Melting at the beginning of each event in the following format:
event #, number of partons in the event.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
5/13/2011 test version v1.25t4/v2.25t5:
* resonances (rho,omega,eta,K*,phi,Delta) produced from quark coalescence
in the String Melting version now have a random mass within a range
(~ +- twice the full Breit-Wigner width);
added for the purpose of resonance reconstruction using invariant mass.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4/15/2010 test version v1.25t4/v2.25t4:
* Darwin added to the list of operating system in `exec';
* crash (especially at low energies) was reported by P. Bhaduri,
modified a statement in subroutine HIJING to solve the problem.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
8/2009 test version v1.25t3/v2.25t3:
* enable embedding of a given number of soft pions along each parton
of the embedded back-to-back jet pair by adding three new
input values in file `input.ampt': nsembd, psembd, tmaxembd.
nsembd: number of soft pions to be embedded with each high-Pt parton
in the embedded jet pair;
psembd: the momentum of each embedded soft pion in GeV/c;
tmaxembd: maximum angle (in radian) of embedded soft pions relative to
the associated high-Pt parton, then the soft pions are generated
uniformly in relative polar and azimutal angle in this cone around the
high-Pt parton. Note that soft pions are not generated uniformly
in solid angle in this cone because that gives a valley at theta=0,
unlike the primitive jet-like correlation with a peak at theta=0.
* enable the embedded jet pair to have random azimuthal angle (still
back-to-back) and to take positions according to a user file
"embed-jet-xy.txt"; enabled by additional values of iembed:
iembed=
1: on with fixed position(xembd,pembd) and Pt(pxqembd,pyqembd);
2: on with fixed position(xembd,pembd) and random azimuthal angle
with Pt-magnitude given by sqrt(pxqembd^2+pyqembd^2);
3: on with random position and fixed Pt(pxqembd,pyqembd);
4: on with random position and random random azimuthal angle
with Pt-magnitude given by sqrt(pxqembd^2+pyqembd^2);
for iembed=3 or 4: need a position file "embed-jet-xy.txt";
Other integers: turn off embedding.
The above modification can be found by searching "clin-8/2009".
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
7/2009 test version v1.25t2/v2.25t2:
* added two new input values in file `input.ampt': ishadow and dshadow.
Setting ishadow=1 (when shadowing flag is on) enables users to set
dshadow (between 0.d0 & 1.d0) to get an intermediate nuclear shadowing
between no-shadowing (when dshadow=0.d0)
and the default HIJING shadowing (when dshadow=1.d0).
The above modification can be found by searching "clin-7/2009".
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
7/2009 test version v1.25/v2.25:
corrected the explanations for deuteron cross section assumptions in
input.ampt and art1f.f
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
6/2009 test version v1.25/v2.25 (==v1.25t1/v2.25t1 meaning "test version 1"):
* added two new input values in file `input.ampt': pttrig and maxmiss.
When enabled, event selection takes place so that each event will have
at least 1 minijet parton with Pt above pttrig in the initial condition
(before parton & hadron cascades).
For this purpose, the HIJING initial condition is repeated until
such an event is generated (but no more than maxmiss times);
* added two new input values in file `input.ampt' so that one can turn off
initial and/or final state radiation as well as the Pt kick
in the HIJING event initial condition; with these effects turned off,
the selected events using pttrig are more likely to have
back-to-back initial minijets;
* added five new input values in file `input.ampt': iembed, pxqembd,
pyqembd, xembd, and yembd. Setting iembed=1 embeds one back-to-back
q/qbar jet pair per event at the specified transverse momentum and
transverse position before parton cascade starts;
* changed the format of zpc.dat, including the following:
event number is given as the first value, and an event-iteration
flag is added as the 2nd value in 1st line of each event in zpc.dat;
an event-iteration flag other than 0 means the event is being repeated
(detailed message is given in the nohup.out file);
* ana/initial_parton_sm.dat is renamed as
ana/parton-initial-afterPropagation.dat,
which outputs the complete information of partons that enter the
parton cascade; i.e., it gives the minijet gluon information
for default runs, or the quark and anti-quark information
after string melting (before parton cascade) for string melting runs.
This option is activated by setting ioscar to 2.
It has the same format as ana/ampt.dat except for the 1st line
for each event, which gives the following:
For isoft=1 (default AMPT):
event number, event-iteration flag, number of partons;
For isoft=4 (String Melting):
event number, event-iteration flag, number of partons,
number of formed baryons, number of formed mesons,
total number of initial particles,
total number of initial particles that cannot enter ZPC.
* added ana/npart-xy.dat, which gives the transverse positions
of all initial nucleons and their status; format is:
For each event, the first line gives:
event #, event-iteration flag, atomic masses of projectile and target,
(impact parameter, added in 4/2012).
The first four values of each line that follows give (for each nucleon):
x,y, sequence number in nucleus (positive for projectile, negative
for target), status (0: spectator, 1 or 2: wounded due to
elastic collisions, 3: wounded due to inelastic collisions);
the 5th value of the line (added in 2/2012) gives the z position
of each nucleon inside a nucleus, with the center of each nucleus
defined as z=0 and the z values of nucleons correspond to those in a
spherical nucleus (i.e. the Lorentz contraction factor of the nucleus
along the z-axis is not included in the z values here).
* ioscar=3 now enables the following output files in addition to
ana/parton-initial-afterPropagation.dat:
all parton collision history in ana/parton-collisionsHistory.dat,
minijet initial condition in ana/minijet-initial-beforePropagation.dat,
ana/parton-after-coalescence.dat for String Melting.
* format of ana/parton-collisionsHistory.dat:
For each collision, the first line gives:
event #, event-iteration flag, sequence numbers of the two partons;
Each of the following four lines gives:
PYTHIA particle ID number, three-momentum(Px,Py,Pz), mass, and
space-time coordinates(x,y,z,t)
for parton1&2 before the collision parton1&2 after the collision.
* format of ana/minijet-initial-beforePropagation.dat:
For each event, the first line gives:
event #, event-iteration flag, atomic masses of projectile and target.
Each of the following lines gives:
PYTHIA particle ID number, three-momentum(Px,Py,Pz), mass, and
space-time coordinates(x,y,z,t) of one minijet parton at production,
ID of origin (from 1: projectile, 2: target; 3: independent string).
* format of ana/parton-after-coalescence.dat:
For each event, the first line gives:
event number, number of partons, number of formed baryons,
number of formed mesons, impact-parameter,
number of participant nucleons in projectile due to elastic
collisions, number of participant nucleons in projectile due to
inelastic collisions, and corresponding numbers in target.
Each of the following lines gives:
PYTHIA parton ID number, three-momentum(Px,Py,Pz), mass,
sequence number of formed hadron, PYTHIA ID number of formed hadron;
* `common /para7/ ioscar' is changed to include string melting info;
* info under "in HJANA1" in file nohup.out is corrected.
The above modifications can be found by searching "clin-6/2009".
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
5/2009 test version v1.24/v2.24:
* Changes made:
* corrected a bug on freezeout time (values in `ampt.dat') of particles
whose formation time is larger than the hadron cascade termination time;
* freezeout time of spectator projectile and target nucleons should be ~0
but it was not correctly updated for string melting runs; now corrected.
The above modifications can be found by searching "clin-5/2009".
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4/2009 update
corrected the explanation on the structure of zpc.dat.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
3/2009 test version v1.23/v2.23:
* Changes made:
* included a subroutine addhad() in amptsub.f to insert
user-defined hadrons before the start of the hadron cascade,
user must set IDPERT to either 1 or 2, user can set NPERTD to 0
if perturbative deuteron productions are not needed;
* add a comment in input.dat saying that setting NPERTD to 0
turns off perturbative deuteron productions;
* moved the location of the subroutine hbtout() in art1f.f,
and as a result there are small changes in hadron freezeout
space-time values, e.g., the freezeout space values may change
by v*DT with v being the particle momentum, and the freezeout time
of spectator projectile or target nucleons change by DT;
* increased the output accuracy of momentum Pz by 1 more digit
in ampt.dat and ampt_pert.dat.
The above modifications can be found by searching "clin-3/2009".
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
10/2008 test version v1.22/v2.22:
* Changes made in the default AMPT model (version 1.22)
and the AMPT model with string melting (version 2.22):
We have included deuteron(d) production and destruction
in the hadronic cascade via d+M <-> B+B, where M represents a meson
including pi, rho, omega and eta; and B represents a baryon
including proton, neutron, Delta, N*(1440) and N*(1535),
Anti-deuteron processes are included similarly.
Elastic collisions of d+M and d+B are also included.
The cross sections for d+pi <-> N+N are based on experimental data;
for other reactions, 4 different assumptions on their cross sections
are available by choosing the value of idxsec in the file `input.ampt'.
Perturbative production/destruction of deuterons can be
turned on by using idpert=1 (or 2) and 1<=NPERTD<=10000 in input.ampt,
and perturbative deuteron results are stored in ana/ampt_pert.dat,
which has the same data format as ana/ampt.dat
(except that ampt_pert.dat does not have the column for deuteron mass
but has the perturbative probability in the last column).
Regular (i.e. non-perturbative) results can be found in `ana/ampt.dat'.
Note that you can use either the regular or perturbative results for
deuterons; you can not use both (that would be double-counting);
regular or perturbative results are close but not exactly the same
due to the biase from triggering perturbative productions.
A deuteron (anti-deuteron) has the particle ID of 42 (-42) in
ana/ampt.dat and ana/ampt_pert.dat.
The above modifications can be found by searching "clin-5/2008",
"clin-6/2008", "clin-8/2008" and "clin-9/2008".
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
10/01/2008 v1.21/v2.21:
* Changes made in the default AMPT model (version 1.21)
and the AMPT model with string melting (version 2.21):
We have added an option to write out the complete parton information
right after string melting (before parton cascade).
Results are stored in ana/initial_parton_sm.dat in the same format
as ana/ampt.dat except for the first few lines that provide
information for each event. This option is only available for
string melting and can be activated by setting ioscar to 2.
We have added an option to turn off phi meson decays at the end
of hadron cascade, i.e., at NT=NTMAX.
This option can be activated by setting iphidcy to 0.
Note that phi decays during hadron cascade are always enabled.
The above two modifications can be found by searching "clin-5b/2008".
Note: the following physics extensions are preliminary and under test,
therefore they have been disabled in v1.21/v2.21:
We have included deuteron(d) production and destruction
in the hadronic cascade via d+M <-> B+B, where M represents a meson
including pi, rho, omega and eta; and B represents a baryon
including proton, neutron, Delta, N*(1440) and N*(1535),
Anti-deuteron processes are included similarly.
Elastic collisions of d+M and d+B are also included.
The above modifications can be found by searching "clin-5/2008",
"clin-6/2008", "clin-8/2008" and "clin-9/2008".
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
04/29/2008 v1.13/v2.13:
* Changes made in the default AMPT model (version 1.13)
and the AMPT model with string melting (version 2.13):
multiple array arguments are checked against out-of-bound in amptsub.f;
commented out the unused "CALL HJAN1A" in zpc.f;
a pause statement in ran1() in zpc.f is modified;
array sizes of ekaon() and sekaon() are increased from 200 to 2000;
array sizes in /HJJET2/ and /xydr/ are increased from 900 to 150001;
multiple compound IF statements are broken up in hipyset1.35.f;
DARWIN added to the list of operating system in `exec';
added a check on the range of IX, IY, IZ in art1f.f and modified
other such checks;
added a check on the range of npion in art1f.f;
RAN() renamed to RANART() to avoid conflict with system functions;
bugs on initializations of xlast() and plast() are fixed;
the variable ISS is modified to avoid out-of-bound error in EKAON();
"IF(K(I,3).EQ.0 .OR. K(K(I,3),2).EQ.IDSTR)" is modified
to avoid out-of-bound error in K();
"DATA NSEED/74769375/" in hijing1.383_ampt.f is commented out;
PYWIDT() subroutine is modified according to pythia-6115.f to avoid
undefined values for variables GGF,GZF,GZPF,ZZF,ZZPF,ZPZPF,API;
"MDCY(KFPR(ISUB,1),1)" is changed to "MDCY(LUCOMP(KFPR(ISUB,1)),1)"
to avoid invalid values for the 1st argument of MDCY();
"if (jscat .ne. 0 .and. next(jscat) .ne. iscat)" is modified
to avoid out-of-bound error in next().
The above modifications can be found by searching "clin-4/2008".
They are not found to change ampt.dat from a few tests on a Linux OS.
We thank A. Vander Molen and G. Westfall for pointing out these issues.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2005 v1.12/v2.12:
* Changes made in the default AMPT model (version 1.12)
and the AMPT model with string melting (version 2.12):
Freezeout time of spectator projectile and target nucleons should be ~0
but it was not correctly updated in `ampt.dat'; now corrected.
The above modifications can be found by searching "clin-12/14/03".
We have corrected a typo in `input.ampt':
"IZT (target A number)" is changed to "IZT (target Z number)".
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
11/30/2004 v1.11/v2.11:
* The default AMPT model (version 1.11) and the AMPT model with
string melting (version 2.11) both use an initialization file
`input.ampt'. The analysis directory `ana/' contains the resulting
data files. The final particle record file is `ana/ampt.dat'.
The version number of AMPT is written to both the file `ana/version'
and `nohup.out'. The AMPT source code has been tested for both the
f77 and the pgf77 compilers on the Unix, Linux, and OSF1 operating systems.
* To run the AMPT program, one needs to:
1. set the initial parameters in `input.ampt'. If one prefers to
use run-time random number seed, set `ihjsed=11', In this way, every
run is different even with the same `input.ampt' file.
2. type `sh exec &' to compile and run the executable `ampt'
with some general information written in `nohup.out'.
* Key initial parameters in `input.ampt' are:
EFRM: sqrt(s_NN) in GeV if FRAME is CMS, e.g. 200 for the maximum RHIC
energy; projectile energy per nucleon if FRAME is LAB.
NEVNT: the total number of events.
BMIN, BMAX: the minimum and maximum impact parameter (in fm)
for all events with BMAX having an upper limit of HIPR1(34)+HIPR1(35)
(=19.87 fm for d+Au collisions and 25.60 fm for Au+Au collisions).
ISOFT: choice of parton-hadron conversion scenario.
=1: default AMPT model (version 1.x);
=4: the AMPT model with string melting (version 2.y).
Note that values of 2, 3, and 5 have never been used for
publications. They are tests of other string melting scenarios:
=2: a string is decomposed into q+dq+minijet partons instead of
using the Lund fragmentation;
=3: a baryon is decomposed into q+qq instead of 3 quarks;
=5: same as 4 but partons freeze out according to
local energy density.
NTMAX: the number of time-steps for hadron cascade, default(D)=150.
Note that NTMAX=3 effectively turns off hadron cascade,
and a larger value than default is usually necessary
for observables at large rapidity or large pseudorapidity.
We use NTMAX=1000 for HBT studies in central Au+Au
collisions due to the need for the last interaction points
and for LHC calculations due to the longer lifetime of the formed matter.
DT: value of the time-step (in fm/c) for hadron cascade, D=0.2.
Note that t_cut=NTMAX*DT is the termination time of hadron cascade.
PARJ(41): parameter a in the Lund symmetric splitting function.
PARJ(42): parameter b in the Lund symmetric splitting function
(in GeV**(-2)). Note that we use default value in HIJING
(a=0.5 and b=0.9) for d+Au collisions,
and a=2.2 and b=0.5 for collisions of heavy nuclei.
flag for popcorn mechanism: D=1(Yes) turns on the popcorn mechanism.
In general, it increases baryon stopping.
PARJ(5): controls BMBbar vs. BBbar in the popcorn mechanism, D=1.0.
shadowing flag: D=1(Yes) turns on nuclear shadowing.
quenching flag: D=0(No) turns off jet quenching
since the parton cascade ZPC simulates final-state effects.
p0 cutoff: D=2.0 (in GeV/c) for p0 in HIJING for minijet production.
parton screening mass: controls the parton cross section,
D=3.2264 (in fm**(-1)). Its square is inversely proportional to
the parton cross section. Use D=3.2264d0 for 3mb, and 2.2814d0 for 6mb.
ihjsed: choice of the random number seed, D=0.
=0: take the `Ran Seed for HIJING' in `input.ampt'
and disregard the random value generated in the file `exec'.
=11: take the HIJING random seed at runtime from the file `exec',
with the seed written in `nohup.out' and `ana/version'.
Ran Seed for HIJING: random number seed for HIJING when ihjsed=0.
Kshort decay flag: depends on the experimental correction procedure,
D=0 turns off Kshort decays after the hadron cascade.
Note that decays of the following resonances and their
antiparticles are always included:
rho, omega, eta, K*, phi, Delta, N*(1440), N*(1535),
Sigma0 (in order to include its feed down to Lambda).
optional OSCAR output: if set to 1, outputs in OSCAR1997A format
are written in `ana/parton.oscar' and `ana/hadron.oscar'.
dpcoal: parton coalescence distance in momentum space (in GeV/c).
drcoal: parton coalescence distance in coordinate space (in fm).
dpcoal, drcoal both have D=10**6 for nearest-neighbor coalescence
in the AMPT model with string melting.
* Key output file are:
ana/ampt.dat: It contains particle records at hadron kinetic freeze-out,
i.e., at the last interaction point.
For each event, the first line gives:
event number, test number(=1), number of particles in the event,
impact parameter, total number of participant nucleons in projectile,
total number of participant nucleons in target, number of participant
nucleons in projectile due to elastic collisions, number of
participant nucleons in projectile due to inelastic collisions,
and corresponding numbers in target.
Note that participant nucleon numbers include nucleons participating
in both elastic and inelastic collisions.
Each of the following lines gives:
PYTHIA particle ID number, three-momentum(Px,Py,Pz), mass, and
space-time coordinates(x,y,z,t) of one final particle at freeze-out.
ana/zpc.dat: similar to `ana/ampt.dat' but for partons at freeze-out.
The first line of each event gives:
event number, (event-iteration flag, added in 6/2009),
number of partons in the event, impact-parameter,
number of participant nucleons in projectile due to elastic
collisions, number of participant nucleons in projectile due to
inelastic collisions, and corresponding numbers in target.
Each of the following lines gives:
For isoft=1 (default AMPT):
PYTHIA particle ID number, three-momentum(Px,Py,Pz), mass, and
space-time coordinates(x,y,z,t) of one final parton at freeze-out.
For isoft=4 (String Melting, but format changed in 4/2012):
PYTHIA particle ID number, three-momentum(Px,Py,Pz), mass,
the parent hadron sequence number that the parton comes from,
the parton sequence number of this parton in the parent hadron
(1-2 for a meson, 1-3 for a baryon), and freeze-out time (t).
Note that momenta are in units of GeV/c, mass in GeV/c**2,
space in fm, and time in fm/c.
If a particle comes from the decay of a resonance which still exists
at the termination time of hadron cascade, then its space-time
corresponds to the decay point of the parent resonance.
Also note that the x-axis in AMPT is defined as the direction along
the impact parameter, and the z-axis is defined as the beam direction.
Please do not hesitate to contact us if needed. Have fun!
Zi-Wei Lin (linz@ecu.edu)
Che-Ming Ko (ko@comp.tamu.edu)
Bao-An Li (Bao-An.Li@tamuc.edu)
Subrata Pal (spal@tifr.res.in)
Bin Zhang (bzhang@astate.edu)
9/16/2018