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_1.f
Go to the documentation of this file.
Or view
the newest version in sPHENIX GitHub for file example_1.f
1
C****************************************************************************
2
C
3
C
4
C
5
C The following is an example program for calling HIJING. one should
6
C include all the common blocks and the data values which are listed
7
C below in his own program.
8
C***************************************************************************
9
10
CHARACTER
frame*8,
proj
*8,
targ
*8
11
common
/himain1/
natt
,
eatt
,
jatt
,
nt
,
np
,
n0
,
n01
,
n10
,
n11
12
SAVE
/himain1/
13
14
common
/himain2/
katt
(130000,4),
patt
(130000,4)
15
SAVE
/himain2/
16
17
C ********information of produced particles
18
C
19
common
/hijjet1/
npj
(300),
kfpj
(300,500),
pjpx
(300,500),
pjpy
(300,500)
20
& ,
pjpz
(300,500),
pjpe
(300,500),
pjpm
(300,500)
21
& ,
ntj
(300),
kftj
(300,500),
pjtx
(300,500),
pjty
(300,500)
22
& ,
pjtz
(300,500),
pjte
(300,500),
pjtm
(300,500)
23
SAVE
/hijjet1/
24
25
common
/hijjet2/
nsg
,
njsg
(900),
iasg
(900,3),
k1sg
(900,100)
26
& ,
k2sg
(900,100),
pxsg
(900,100),
pysg
(900,100),
pzsg
(900,100)
27
& ,
pesg
(900,100),
pmsg
(900,100)
28
SAVE
/hijjet2/
29
30
C ********information of produced partons
31
32
common
/ranseed/
nseed
! Uzhi
33
SAVE
/ranseed/
! Uzhi
34
nseed
=0
! Uzhi
35
36
efrm
=200.0
37
frame=
'CMS'
38
proj
=
'P'
! 'A' ! Uzhi
39
targ
=
'P'
! 'A' ! Uzhi
40
iap
=1
! 197 ! Uzhi
41
izp
=1
! 79 ! Uzhi
42
iat
=1
! 197 ! Uzhi
43
izt
=1
! 79 ! Uzhi
44
n_events=100
! Uzhi
45
46
C Simulation of PP-interactions at \sqrt{s}=200 GeV ! Uzhi
47
48
49
CALL
hijset
(
efrm
,frame,
proj
,
targ
,
iap
,
izp
,
iat
,
izt
)
50
C ********Initialize HIJING
51
52
WRITE
(6,*)
' Simulation of interactions with'
! Uzhi
53
WRITE
(6,*)
! Uzhi
54
WRITE
(6,*)
' Proj = '
,
proj
,
' and Targ = '
,
targ
! Uzhi
55
WRITE
(6,*)
' IAP ='
,
iap
,
' IAT ='
,
iat
! Uzhi
56
WRITE
(6,*)
' IZP ='
,
izp
,
' IZT ='
,
izt
! Uzhi
57
WRITE
(6,*)
! Uzhi
58
WRITE
(6,*)
' Reference frame - '
,frame
! Uzhi
59
WRITE
(6,*)
' ENERGY '
,
efrm
,
' GeV'
! Uzhi
60
WRITE
(6,*)
' Number of generated events -'
,n_events
! Uzhi
61
WRITE
(6,*)
! Uzhi
62
63
bmin=0.0
64
bmax=0.0
65
DO
2000 i_event=1,n_events
66
67
WRITE
(6,*)
' Event # '
,i_event,
' ------------------'
! Uzhi
68
69
CALL
hijing
(frame,bmin,bmax)
70
C
71
WRITE
(6,*)
' Multiplicity of produced particles - '
,
natt
! Uzhi
72
write
(6,*)
! Uzhi
73
WRITE
(6,*)
' ID Charge Mass (GeV) Px Py'
,
! Uzhi
74
&
' Pz (GeV/c)'
! Uzhi
75
WRITE
(6,*)
' ----------------------------------------'
! Uzhi
76
77
78
DO
1000 i=1,
natt
79
! Uzhi
80
C IF(LUCHGE(KATT(I,1)).NE.0) THEN ! Uzhi
81
C this select charged particles only
82
C !information of produced particles
83
C !is stored in common blocks HIMAIN1 and
84
C HIMAIN2
85
C ENDIF ! Uzhi
86
87
ich=
luchge
(
katt
(i,1))/3
! Uzhi
88
89
amass=
ulmass
(
katt
(i,1))
! Uzhi
90
91
WRITE
(6,900)
katt
(i,1),ich,amass,
! Uzhi
92
&
patt
(i,1),
patt
(i,2),
patt
(i,3)
! Uzhi
93
900
FORMAT
(1
x
,i6,i4,2
x
,f7.3,3(2
x
,e11.4))
! Uzhi
94
95
1000
CONTINUE
96
97
pause
! Uzhi
98
99
2000
CONTINUE
100
stop
101
END
102
103
FUNCTION
ran
(NSEED) ! Uzhi
104
ran
=
rlu
(
nseed
)
! Uzhi
105
RETURN
! Uzhi
106
END
! Uzhi
coresoftware
blob
master
generators
hijing
examples
example_1.f
Built by
Jin Huang
. updated:
Wed Jun 29 2022 17:24:29
using
1.8.2 with
ECCE GitHub integration