ECCE @ EIC Software
Reference for
ECCE @ EIC
simulation and reconstruction software on GitHub
Home page
Related Pages
Modules
Namespaces
Classes
Files
External Links
File List
File Members
All
Classes
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Groups
Pages
example_4.f
Go to the documentation of this file.
Or view
the newest version in sPHENIX GitHub for file example_4.f
1
C****************************************************************************
2
C Program # 3 from Comp. Phys. Commun. 83 (1994) 307
3
C by M. Gyulassy and X-.N. Wang
4
C Modified by V.Uzhinsky, CERN, Oct. 2003
5
C***************************************************************************
6
7
CHARACTER
frame*8,
proj
*8,
targ
*8
8
9
common
/hiparnt/
hipr1
(100),
ihpr2
(50),
hint1
(100),
ihnt2
(50)
10
SAVE
/hiparnt/
11
12
common
/ranseed/
nseed
13
SAVE
/ranseed/
14
15
CHARACTER *1
key
16
17
nseed
=0
18
19
C....switch on triggered jet production:
20
ihpr2
(3)=1
21
22
frame=
'CMS'
23
24
write
(6,*)
'===================================================='
25
write
(6,*)
' Simulation of events with triggered hard processes '
26
write
(6,*)
' in hh-, hA- and AA-collisions '
27
write
(6,*)
' Calculation will be performed in CM System '
28
write
(6,*)
'===================================================='
29
30
write
(6,*)
31
write
(6,*)
'Enter the energy per NN-collision (GeV)'
32
read
(5,*)
efrm
33
34
write
(6,*)
35
write
(6,*)
'Enter a type of the "projectile" particle'
36
write
(6,*)
37
write
(6,*)
' P proton, PBAR anti-proton,'
38
write
(6,*)
' N neutron, NBAR anti-neutron,'
39
write
(6,*)
' PI+ - positive pion, PI- negative pion,'
40
write
(6,*)
' K+ positive kaon, K- negative kaon'
41
write
(6,*)
42
write
(6,*)
' A - nucleus --------------------------'
43
44
read
(5,1)
proj
45
1
format
(a8)
46
47
if
(
proj
.ne.
'A'
)
then
48
iap
=1
49
if
(
proj
.eq.
'P'
)
izp
= 1
50
if
(
proj
.eq.
'PBAR'
)
izp
=-1
51
if
(
proj
.eq.
'N'
)
izp
= 0
52
if
(
proj
.eq.
'NBAR'
)
izp
= 0
53
if
(
proj
.eq.
'PI+'
)
izp
= 1
54
if
(
proj
.eq.
'PI-'
)
izp
=-1
55
if
(
proj
.eq.
'K+'
)
izp
= 1
56
if
(
proj
.eq.
'K-'
)
izp
=-1
57
else
58
write
(6,*)
59
write
(6,*)
'Enter mass number and charge of the proj. nucleus'
60
read
(5,*)
iap
,
izp
61
endif
62
63
write
(6,*)
64
write
(6,*)
'Enter a type of the "target" particle (same notations)'
65
read
(5,1)
targ
66
67
if
(
targ
.ne.
'A'
)
then
68
iat
=1
69
if
(
targ
.eq.
'P'
)
izt
= 1
70
if
(
targ
.eq.
'PBAR'
)
izt
=-1
71
if
(
targ
.eq.
'N'
)
izt
= 0
72
if
(
targ
.eq.
'NBAR'
)
izt
= 0
73
if
(
targ
.eq.
'PI+'
)
izt
= 1
74
if
(
targ
.eq.
'PI-'
)
izt
=-1
75
if
(
targ
.eq.
'K+'
)
izt
= 1
76
if
(
targ
.eq.
'K-'
)
izt
=-1
77
else
78
write
(6,*)
79
write
(6,*)
'Enter mass number and charge of the target nucleus'
80
read
(5,*)
iat
,
izt
81
endif
82
83
if
(
proj
.eq.
'A'
.or.
targ
.eq.
'A'
)
then
84
write
(6,*)
'Enter Min. and Max. values of impact parameter (fm)'
85
read
(5,*) bmin, bmax
86
else
87
bmin=0.0
88
bmax=0.0
89
endif
90
91
if
(
proj
.eq.
'A'
.or.
targ
.eq.
'A'
)
then
92
write
(6,*)
93
write
(6,*)
' Would you like to take jet quenching into account?'
94
write
(6,*)
' Y - Yes, N - No? '
95
96
read
(5,2)
key
97
2
format
(a1)
98
99
if
(
key
.eq.
'Y'
.or.
key
.eq.
'y'
)
then
100
ihpr2
(4)=1
101
else
102
ihpr2
(4)=0
103
endif
104
endif
105
106
C....set the pt range of the triggered jets:
107
write
(6,*)
108
write
(6,*)
' Set the Pt of triggered jets (in GeV/c)'
109
read
(5,*) pt_trigger
110
111
write
(6,*)
112
write
(6,*)
'Enter number of events'
113
read
(5,*) n_events
114
115
C....initialize HIJING for requested interactions
116
CALL
hijset
(
efrm
,frame,
proj
,
targ
,
iap
,
izp
,
iat
,
izt
)
117
118
write
(6,*)
' Simulation of interactions with'
119
write
(6,*)
120
write
(6,*)
' Proj = '
,
proj
,
' and Targ = '
,
targ
121
write
(6,*)
' IAP ='
,
iap
,
' IAT ='
,
iat
122
write
(6,*)
' IZP ='
,
izp
,
' IZT ='
,
izt
123
write
(6,*)
124
write
(6,*)
' Reference frame - '
,frame
125
write
(6,*)
' ENERGY '
,
efrm
,
' GeV'
126
write
(6,*)
' Number of generated events'
,n_events
127
write
(6,*)
' Triggered jets Pt =>'
,pt_trigger,
' (GeV/c)'
128
write
(6,*)
129
130
DO
100 i_event=1,n_events
131
132
write
(6,*)
'Event # '
,i_event,
' -------------------------'
133
write
(6,*)
134
135
hipr1
(10)=-pt_trigger
136
137
CALL
hijing
(frame,bmin,bmax)
138
139
C....print out flavor code of the first jet:
140
write
(6,*)
'Flavor code of the first jet: '
,
ihnt2
(9)
141
142
C....and its four momentum:
143
write
(6,*)
' Px Py Pz'
,
144
,
' E (GeV/c, GeV)'
145
write
(6,3)
hint1
(21),
hint1
(22),
hint1
(23),
hint1
(24)
146
3
format
(4(e11.4,2
x
))
147
write
(6,*)
148
149
C....print out flavor code of the second jet:
150
write
(6,*)
'Flavor code of the second jet:'
,
ihnt2
(10)
151
152
C....and its four momentum:
153
write
(6,*)
' Px Py Pz'
,
154
,
' E (GeV/c, GeV)'
155
write
(6,3)
hint1
(31),
hint1
(32),
hint1
(33),
hint1
(34)
156
write
(6,*)
157
158
pause
159
100
CONTINUE
160
161
stop
162
END
163
164
FUNCTION
ran
(NSEED)
165
ran
=
rlu
(
nseed
)
166
RETURN
167
END
coresoftware
blob
master
generators
hijing
examples
example_4.f
Built by
Jin Huang
. updated:
Wed Jun 29 2022 17:24:29
using
1.8.2 with
ECCE GitHub integration