Skip to content

Commit ed25661

Browse files
njacazionjacazio
committed
Fixes
1 parent 70395f5 commit ed25661

2 files changed

Lines changed: 156 additions & 4 deletions

File tree

‎MC/config/common/ampt/generator/pp_536TeV_AMPT.ampt‎

Lines changed: 155 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,155 @@
1+
5360 ! EFRM (sqrt(S_NN) in GeV if FRAME is CMS)
2+
CMS ! FRAME
3+
A ! PROJ
4+
A ! TARG
5+
1 ! IAP (projectile A number)
6+
1 ! IZP (projectile Z number)
7+
1 ! IAT (target A number)
8+
1 ! IZT (target Z number)
9+
1 ! NEVNT (total number of events) - overwritten by generator_AMPT.C
10+
0.0 ! BMIN (mininum impact parameter in fm)
11+
2.5 ! BMAX (maximum impact parameter in fm, also see below)
12+
4 ! ISOFT (D=4): select Default AMPT or String Melting(see below)
13+
150 ! NTMAX: number of timesteps (D=150), see below
14+
0.2 ! DT: timestep in fm (hadron cascade time= DT*NTMAX) (D=0.2)
15+
0.30 ! PARJ(41): parameter a in Lund symmetric splitting function
16+
0.15 ! PARJ(42): parameter b in Lund symmetric splitting function
17+
1 ! (D=1,yes;0,no) flag for popcorn mechanism(netbaryon stopping)
18+
1.0 ! PARJ(5) to control BMBbar vs BBbar in popcorn (D=1.0)
19+
1 ! shadowing flag (Default=1,yes; 0,no)
20+
0 ! quenching flag (D=0,no; 1,yes)
21+
2.0 ! quenching parameter -dE/dx (GeV/fm) in case quenching flag=1
22+
2.0 ! p0 cutoff in HIJING for minijet productions (D=2.0)
23+
2.265d0 ! parton screening mass in fm^(-1) (D=2.265d0), see below
24+
0 ! IZPC: (D=0 forward-angle parton scatterings; 100,isotropic)
25+
0.33d0 ! alpha in parton cascade (D=0.33d0), see parton screening mass
26+
1d6 ! dpcoal in GeV
27+
1d6 ! drcoal in fm
28+
0 ! ihjsed: take HIJING seed from below (D=0)or at runtime(11) - overwritten by generator_AMPT.C
29+
13150909 ! random seed for HIJING - overwritten by generator_AMPT.C
30+
8 ! random seed for parton cascade - overwritten by generator_AMPT.C
31+
0 ! flag for K0s weak decays (D=0,no; 1,yes)
32+
1 ! flag for phi decays at end of hadron cascade (D=1,yes; 0,no)
33+
0 ! flag for pi0 decays at end of hadron cascade (D=0,no; 1,yes)
34+
0 ! optional OSCAR output (D=0,no; 1,yes; 2&3,more parton info)
35+
0 ! flag for perturbative deuteron calculation (D=0,no; 1or2,yes)
36+
1 ! integer factor for perturbative deuterons(>=1 & <=10000)
37+
1 ! choice of cross section assumptions for deuteron reactions
38+
-7. ! Pt in GeV: generate events with >=1 minijet above this value
39+
1000 ! maxmiss (D=1000): maximum # of tries to repeat a HIJING event
40+
3 ! flag on initial and final state radiation (D=3,both yes; 0,no)
41+
1 ! flag on Kt kick (D=1,yes; 0,no)
42+
0 ! flag to turn on quark pair embedding (D=0,no; 1,yes)
43+
7., 0. ! Initial Px and Py values (GeV) of the embedded quark (u or d)
44+
0., 0. ! Initial x & y values (fm) of the embedded back-to-back q/qbar
45+
1, 5., 0. ! nsembd(D=0), psembd (in GeV),tmaxembd (in radian).
46+
0 ! Flag to enable users to modify shadowing (D=0,no; 1,yes)
47+
1.d0 ! Factor used to modify nuclear shadowing
48+
0 ! Flag for random orientation of reaction plane (D=0,no; 1,yes)
49+
0 ! iexotnp, Flag for exotic nuclear structure of projectile (=0, no; =1,yes)
50+
0 ! iexotnt, Flag for exotic nuclear structure of target (=0, no; =1,yes)
51+
0 ! nstrucp, nucl. struct. of projectile, C12, 1(Triangle), 2(Chain); O16, 1(Tetrahedron), 2(Square), 3(Kite), 4(Chain)
52+
0 ! nstruct, nucl. struct. of target, C12, 1(Triangle), 2(Chain); O16, 1(Tetrahedron), 2(Square), 3(Kite), 4(Chain)
53+
54+
%%%%%%%%%% O2DPG notes:
55+
Pb-Pb at sqrt(s_NN) = 5.36 TeV, minimum bias (b = 0-20 fm), String Melting
56+
with the LHC settings of arXiv:1403.6321 (a=0.30, b=0.15/GeV^2, 1.5 mb parton
57+
cross section: alpha=0.33 and screening mass 3.2264/fm).
58+
Values are read by AMPT line by line: do NOT add or remove lines above.
59+
NEVNT (line 9), ihjsed (line 28) and the two seeds (lines 29-30) are
60+
overwritten by MC/config/common/external/generator/generator_AMPT.C.
61+
%%%%%%%%%% Further explanations:
62+
BMAX: the upper limit HIPR1(34)+HIPR1(35)=19.87fm (dAu), 25.60fm(AuAu).
63+
ISOFT: 1 Default,
64+
4 String Melting.
65+
PARJ(41) & (42): for string melting AMPT, 0.55 & 0.15/GeV^2 are recommended
66+
for top RHIC energies and 0.30 & 0.15/GeV^2 are recommended for
67+
LHC energies (see arXiv:1403.6321 for details).
68+
NTMAX: number of time-steps for hadron cascade.
69+
Use a large value (e.g. 1000) for LHC studies or HBT studies at RHIC.
70+
Using NTMAX=2 or 3 effectively turns off hadronic cascade.
71+
parton screening mass (in 1/fm): its square is inversely proportional to
72+
the parton cross section. Use D=2.265d0 for 3mb cross section
73+
when alpha in parton cascade is set to 0.33;
74+
(note: 3.2264d0 for 3mb cross section when alpha is set to 0.47).
75+
Using 1d4 effectively turns off parton cascade.
76+
ihjsed: if =11, take HIJING random seed at runtime so that
77+
every run may be automatically different (see file 'exec').
78+
iksdcy: flag for K0s weak decays for comparison with data.
79+
iphidcy: flag for phi meson decays at the end of hadron cascade for comparison
80+
with data; default is yes; use 0 to turn off these decays.
81+
Note: phi meson decay during hadron cascade is always enabled.
82+
ipi0dcy: flag for pi0 electromagnetic decays at the end of hadron cascade for
83+
comparison with data; set to 1 to turn on pi0 decays.
84+
ioscar: 0 Dafault,
85+
1 Write output in the OSCAR format,
86+
2 Write out the complete parton information
87+
(ana/parton-initial-afterPropagation.dat)
88+
right after string melting (before parton cascade),
89+
3 Write out several more files on parton information (see readme).
90+
idpert: flag for perturbative deuteron and antideuteron calculations
91+
with results in ana/ampt_pert.dat:
92+
0 No perturbative calculations,
93+
1 Trigger a production of NPERTD perturbative deuterons
94+
in each NN collision,
95+
2 Trigger a production of NPERTD perturbative deuterons only in
96+
an NN collision where a conventional deuteron is produced.
97+
Note: conventional deuteron calculations are always performed
98+
with results in ana/ampt.dat.
99+
NPERTD: number of perturbative deuterons produced in each triggered collision;
100+
setting it to 0 turns off perturbative deuteron productions.
101+
idxsec: choose a cross section model for deuteron inelastic/elastic collisions:
102+
1: same |matrix element|**2/s (after averaging over initial spins
103+
and isospins) for B+B -> deuteron+meson at the same sqrt(s);
104+
2: same |matrix element|**2/s for B+B -> deuteron+meson
105+
at the same sqrt(s)-threshold;
106+
3: same |matrix element|**2/s for deuteron+meson -> B+B
107+
at the same sqrt(s);
108+
4: same |matrix element|**2/s for deuteron+meson -> B+B
109+
at the same sqrt(s)-threshold;
110+
1 or 3 also chooses the same cross section for deuteron+meson or baryon
111+
elastic collision at the same sqrt(s);
112+
2 or 4 also chooses the same cross section for deuteron+meson or baryon
113+
elastic collision at the same sqrt(s)-threshold.
114+
%%%%%%%%%% For jet studies:
115+
pttrig: generate events with at least 1 initial minijet parton above this Pt
116+
value, otherwise repeat HIJING event until reaching maxmiss tries;
117+
use a negative value to disable this requirement and get normal events.
118+
maxmiss: maximum number of tries for the repetition of a HIJING event to obtain
119+
a minijet above the Pt value of pttrig; increase maxmiss if some events
120+
fail to generate at least 1 initial minijet parton above pttrig.
121+
it is safer to set a large value for high pttrig and/or large b value
122+
and/or smaller colliding nuclei.
123+
IHPR2(2): flag to turn off initial and final state radiation:
124+
0 both radiation off, 1 only final off, 2 only initial off, 3 both on.
125+
IHPR2(5): flag to turn off Pt kick due to soft interactions: 0 off, 1 on.
126+
Setting both IHPR2(2) and IHPR2(5) to zero makes it more likely to
127+
have two high-Pt minijet partons that are close to back-to-back.
128+
%%%%%%%%%% To embed a back-to-back light q/qbar jet pair
129+
%%%%%%%%%% and a given number of soft pions along each jet into each event:
130+
iembed: flag to turn on quark pair embedding:
131+
1: on with fixed position(xembd,pembd) and Pt(pxqembd,pyqembd);
132+
2: on with fixed position(xembd,pembd) and random azimuthal angle
133+
with Pt-magnitude given by sqrt(pxqembd^2+pyqembd^2);
134+
3: on with random position and fixed Pt(pxqembd,pyqembd);
135+
4: on with random position and random random azimuthal angle
136+
with Pt-magnitude given by sqrt(pxqembd^2+pyqembd^2);
137+
for iembed=3 or 4: need a position file "embed-jet-xy.txt";
138+
Other integers: off.
139+
pxqembd, pyqembd: sqrt(pxqembd^2+pyqembd^2) > 70MeV/c is required;
140+
the embedded quark and antiquark have pz=0.
141+
xembd, yembd: the embedded quark and antiquark jets have z=0 initially. Note:
142+
the x-axis is defined as the direction along the impact parameter.
143+
nsembd: number of soft pions to be embedded with each high-Pt parton
144+
in the embedded jet pair.
145+
psembd: Momentum of each embedded soft pion in GeV.
146+
tmaxembd: maximum angle(rad) of embedded soft pions relative to high-Pt parton.
147+
%%%%%%%%%% User modification of nuclear shadowing:
148+
ishadow: set to 1 to enable users to adjust nuclear shadowing
149+
provided the shadowing flag IHPR2(6) is turned on; default value is 0.
150+
dshadow: valid when ishadow=1; this parameter modifies the HIJING shadowing
151+
parameterization Ra(x,r)==1+fa(x,r) via Ra(x,r)==1+fa(x,r)*dshadow,
152+
so the value of 0.d0 turns off shadowing
153+
and the value of 1.d0 uses the default HIJING shadowing;
154+
currently limited to 0.d0<=dshadow<=1.d0 to make sure Ra(x,r)>0.
155+
iphirp: set to 1 to turn on random orientation of reaction plane (D=0)

‎MC/config/common/ini/GeneratorAMPTpp536TeV.ini‎

Lines changed: 1 addition & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -1,6 +1,3 @@
1-
#---> GeneratorAMPTPbPb536TeV
2-
#NEV_TEST> 1000
31
[GeneratorExternal]
42
fileName=${O2DPG_MC_CONFIG_ROOT}/MC/config/common/external/generator/generator_AMPT.C
5-
# Grid path to the AMPT input file
6-
funcName=generateAMPT("alien:///alice/cern.ch/user/n/njacazio/ampt_nucleiSim/input_pp5360_sm.ampt")
3+
funcName=generateAMPT("${O2DPG_MC_CONFIG_ROOT}/MC/config/common/ampt/generator/pp_536TeV_AMPT.ampt")

0 commit comments

Comments
 (0)