6 C****************************************************************
7 C conduct soft radiation according to dipole approxiamtion
8 C****************************************************************
9 SUBROUTINE ATTRAD(IERROR)
12 #include "hiparnt.inc"
14 #include "lujets_hijing.inc"
18 C.....S INVARIANT MASS-SQUARED BETWEEN PARTONS I AND I+1......
19 C.....SM IS THE LARGEST MASS-SQUARED....
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
27 WP=SQRT(S)-1.5*(P(I,5)+P(I+1,5))
33 BTT=(PBT1**2+PBT2**2+PBT3**2)/PBT4**2
34 IF(BTT.GE.1.0-1.0E-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
41 S=(SM+1.5*(P(JL,5)+P(JL+1,5)))**2
42 IF(SM.LT.HIPR1(5)) GOTO 2
44 C.....MAKE PLACE FOR ONE GLUON.....
45 IF(JL+1.EQ.N) GOTO 190
54 C.....BOOST TO REST SYSTEM FOR PARTICLES JL AND JL+1.....
64 CALL ATROBO(0.,0.,BEX,BEY,BEZ,IMIN,IMAX,IERROR)
65 IF(IERROR.NE.0) RETURN
66 C.....ROTATE TO Z-AXIS....
67 PABS=(P(JL,4)-P(JL,5))*(P(JL,4)+P(JL,5))
69 CTH=P(JL,3)/SQRT(PABS)
73 IF(ABS(CTH).GT.1.0) CTH=MAX(-1.,MIN(1.,CTH))
75 PHI=ULANGL_HIJING(P(JL,1),P(JL,2))
76 CALL ATROBO(0.,-PHI,0.,0.,0.,IMIN,IMAX,IERROR)
77 CALL ATROBO(-THETA,0.,0.,0.,0.,IMIN,IMAX,IERROR)
79 C.....CREATE ONE GLUON AND ORIENTATE.....
81 1 CALL AR3JET(S,X1,X3,JL)
82 CALL ARORIE(S,X1,X3,JL)
83 IF(HIDAT(2).GT.0.0) THEN
84 PTG1=SQRT(P(JL,1)**2+P(JL,2)**2)
85 PTG2=SQRT(P(JL+1,1)**2+P(JL+1,2)**2)
86 PTG3=SQRT(P(JL+2,1)**2+P(JL+2,2)**2)
87 PTG=MAX(PTG1,PTG2,PTG3)
88 IF(PTG.GT.HIDAT(2)) THEN
89 FMFACT=EXP(-(PTG**2-HIDAT(2)**2)/HIPR1(2)**2)
90 IF(RLU_HIJING(0).GT.FMFACT) GO TO 1
93 C.....ROTATE AND BOOST BACK.....
96 CALL ATROBO(THETA,PHI,-BEX,-BEY,-BEZ,IMIN,IMAX,IERROR)
97 IF(IERROR.NE.0) RETURN
98 C.....ENUMERATE THE GLUONS.....
111 C----THETA FUNCTION DAMPING OF THE EMITTED GLUONS. FOR HADRON-HADRON.
113 C IF(VFR(2).GT.0.) THEN
114 C PTG=SQRT(P(JL+1,1)**2+P(JL+1,2)**2)
116 C DOPT=SQRT((4.*PAR(71)*VFR(2))/WSTRI)
117 C PTOPT=(DOPT*WSTRI)/(2.*VFR(2))
118 C IF(PTG.GT.PTOPT) IORDER=IORDER-1
119 C IF(PTG.GT.PTOPT) GOTO 1
122 IF(SM.GE.HIPR1(5)) GOTO 40