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
attrad.f
Go to the documentation of this file.
Or view
the newest version in sPHENIX GitHub for file attrad.f
1
2
C
3
C
4
C
5
C****************************************************************
6
C conduct soft radiation according to dipole approxiamtion
7
C****************************************************************
8
SUBROUTINE
attrad
(IERROR)
9
C
10
common
/hiparnt/
hipr1
(100),
ihpr2
(50),
hint1
(100),
ihnt2
(50)
11
SAVE
/hiparnt/
12
common
/hijdat/hidat0(10,10),hidat(10)
13
SAVE
/hijdat/
14
common
/lujets/
n
,
k
(9000,5),
p
(9000,5),
v
(9000,5)
15
SAVE
/lujets/
16
ierror=0
17
18
C.....S INVARIANT MASS-SQUARED BETWEEN PARTONS I AND I+1......
19
C.....SM IS THE LARGEST MASS-SQUARED....
20
21
40 sm=0.
22
jl=1
23
DO
30 i=1,
n
-1
24
s
=2.*(
p
(i,4)*
p
(i+1,4)-
p
(i,1)*
p
(i+1,1)-
p
(i,2)*
p
(i+1,2)
25
& -
p
(i,3)*
p
(i+1,3))+
p
(i,5)**2+
p
(i+1,5)**2
26
IF
(
s
.LT.0.)
s
=0.
27
wp
=
sqrt
(
s
)-1.5*(
p
(i,5)+
p
(i+1,5))
28
IF
(
wp
.GT.sm)
THEN
29
pbt1=
p
(i,1)+
p
(i+1,1)
30
pbt2=
p
(i,2)+
p
(i+1,2)
31
pbt3=
p
(i,3)+
p
(i+1,3)
32
pbt4=
p
(i,4)+
p
(i+1,4)
33
btt=(pbt1**2+pbt2**2+pbt3**2)/pbt4**2
34
IF
(btt.GE.1.0-1.0
e
-10) go to 30
35
IF
((i.NE.1.OR.i.NE.
n
-1).AND.
36
& (
k
(i,2).NE.21.AND.
k
(i+1,2).NE.21)) go to 30
37
jl=i
38
sm=
wp
39
ENDIF
40
30
CONTINUE
41
s
=(sm+1.5*(
p
(jl,5)+
p
(jl+1,5)))**2
42
IF
(sm.LT.
hipr1
(5)) goto 2
43
44
C.....MAKE PLACE FOR ONE GLUON.....
45
IF
(jl+1.EQ.
n
) goto 190
46
DO
160 j=
n
,jl+2,-1
47
k
(j+1,1)=
k
(j,1)
48
k
(j+1,2)=
k
(j,2)
49
DO
150
m
=1,5
50
C+++BAC
51
v
(j+1,
m
) =
v
(j,
m
)
52
C---BAC
53
150
p
(j+1,
m
)=
p
(j,
m
)
54
160
CONTINUE
55
190
n
=
n
+1
56
57
C.....BOOST TO REST SYSTEM FOR PARTICLES JL AND JL+1.....
58
p1
=
p
(jl,1)+
p
(jl+1,1)
59
p2
=
p
(jl,2)+
p
(jl+1,2)
60
p3
=
p
(jl,3)+
p
(jl+1,3)
61
p4=
p
(jl,4)+
p
(jl+1,4)
62
bex=-
p1
/p4
63
bey=-
p2
/p4
64
bez=-
p3
/p4
65
imin=jl
66
imax
=jl+1
67
CALL
atrobo
(0.,0.,bex,bey,bez,imin,
imax
,ierror)
68
IF
(ierror.NE.0)
RETURN
69
C.....ROTATE TO Z-AXIS....
70
cth=
p
(jl,3)/
sqrt
(
p
(jl,4)**2-
p
(jl,5)**2)
71
IF
(abs(cth).GT.1.0) cth=
max
(-1.,
min
(1.,cth))
72
theta
=acos(cth)
73
phi
=
ulangl
(
p
(jl,1),
p
(jl,2))
74
CALL
atrobo
(0.,-
phi
,0.,0.,0.,imin,
imax
,ierror)
75
CALL
atrobo
(-
theta
,0.,0.,0.,0.,imin,
imax
,ierror)
76
77
C.....CREATE ONE GLUON AND ORIENTATE.....
78
79
1 CALL
ar3jet
(
s
,
x1
,
x3
,jl)
80
CALL
arorie
(
s
,
x1
,
x3
,jl)
81
IF
(hidat(2).GT.0.0)
THEN
82
ptg1=
sqrt
(
p
(jl,1)**2+
p
(jl,2)**2)
83
ptg2=
sqrt
(
p
(jl+1,1)**2+
p
(jl+1,2)**2)
84
ptg3=
sqrt
(
p
(jl+2,1)**2+
p
(jl+2,2)**2)
85
ptg=
max
(ptg1,ptg2,ptg3)
86
IF
(ptg.GT.hidat(2))
THEN
87
fmfact=
exp
(-(ptg**2-hidat(2)**2)/
hipr1
(2)**2)
88
IF
(
atl_ran
(
nseed
).GT.fmfact) go to 1
89
ENDIF
90
ENDIF
91
C.....ROTATE AND BOOST BACK.....
92
imin=jl
93
imax
=jl+2
94
CALL
atrobo
(
theta
,
phi
,-bex,-bey,-bez,imin,
imax
,ierror)
95
IF
(ierror.NE.0)
RETURN
96
C.....ENUMERATE THE GLUONS.....
97
k
(jl+2,1)=
k
(jl+1,1)
98
k
(jl+2,2)=
k
(jl+1,2)
99
k
(jl+2,3)=
k
(jl+1,3)
100
k
(jl+2,4)=
k
(jl+1,4)
101
k
(jl+2,5)=
k
(jl+1,5)
102
p
(jl+2,5)=
p
(jl+1,5)
103
k
(jl+1,1)=2
104
k
(jl+1,2)=21
105
k
(jl+1,3)=0
106
k
(jl+1,4)=0
107
k
(jl+1,5)=0
108
p
(jl+1,5)=0.
109
110
v
(jl+1,1) = 0.
111
v
(jl+1,2) = 0.
112
v
(jl+1,3) = 0.
113
v
(jl+1,4) = 0.
114
v
(jl+1,5) = 0.
115
116
C----THETA FUNCTION DAMPING OF THE EMITTED GLUONS. FOR HADRON-HADRON.
117
C----R0=VFR(2)
118
C IF(VFR(2).GT.0.) THEN
119
C PTG=SQRT(P(JL+1,1)**2+P(JL+1,2)**2)
120
C PTGMAX=WSTRI/2.
121
C DOPT=SQRT((4.*PAR(71)*VFR(2))/WSTRI)
122
C PTOPT=(DOPT*WSTRI)/(2.*VFR(2))
123
C IF(PTG.GT.PTOPT) IORDER=IORDER-1
124
C IF(PTG.GT.PTOPT) GOTO 1
125
C ENDIF
126
C-----
127
IF
(sm.GE.
hipr1
(5)) goto 40
128
129
2
k
(1,1)=2
130
k
(1,3)=0
131
k
(1,4)=0
132
k
(1,5)=0
133
k
(
n
,1)=1
134
k
(
n
,3)=0
135
k
(
n
,4)=0
136
k
(
n
,5)=0
137
138
RETURN
139
END
coresoftware
blob
master
generators
hijing
src
attrad.f
Built by
Jin Huang
. updated:
Wed Jun 29 2022 17:24:29
using
1.8.2 with
ECCE GitHub integration