6 C****************************************************************
7 C conduct soft radiation according to dipole approxiamtion
8 C****************************************************************
9 SUBROUTINE ATTRAD(IERROR)
11 #include "hiparnt.inc"
13 #include "lujets_hijing.inc"
17 C.....S INVARIANT MASS-SQUARED BETWEEN PARTONS I AND I+1......
18 C.....SM IS THE LARGEST MASS-SQUARED....
22 S=2.*(P(I,4)*P(I+1,4)-P(I,1)*P(I+1,1)-P(I,2)*P(I+1,2)
23 & -P(I,3)*P(I+1,3))+P(I,5)**2+P(I+1,5)**2
25 WP=SQRT(S)-1.5*(P(I,5)+P(I+1,5))
31 BTT=(PBT1**2+PBT2**2+PBT3**2)/PBT4**2
32 IF(BTT.GE.1.0-1.0E-10) GO TO 30
33 IF((I.NE.1.OR.I.NE.N-1).AND.
34 & (K(I,2).NE.21.AND.K(I+1,2).NE.21)) GO TO 30
39 S=(SM+1.5*(P(JL,5)+P(JL+1,5)))**2
40 IF(SM.LT.HIPR1(5)) GOTO 2
42 C.....MAKE PLACE FOR ONE GLUON.....
43 IF(JL+1.EQ.N) GOTO 190
52 C.....BOOST TO REST SYSTEM FOR PARTICLES JL AND JL+1.....
62 CALL ATROBO(0.,0.,BEX,BEY,BEZ,IMIN,IMAX,IERROR)
63 IF(IERROR.NE.0) RETURN
64 C.....ROTATE TO Z-AXIS....
65 CTH=P(JL,3)/SQRT(P(JL,4)**2-P(JL,5)**2)
66 IF(ABS(CTH).GT.1.0) CTH=MAX(-1.,MIN(1.,CTH))
68 PHI=ULANGL_HIJING(P(JL,1),P(JL,2))
69 CALL ATROBO(0.,-PHI,0.,0.,0.,IMIN,IMAX,IERROR)
70 CALL ATROBO(-THETA,0.,0.,0.,0.,IMIN,IMAX,IERROR)
72 C.....CREATE ONE GLUON AND ORIENTATE.....
74 1 CALL AR3JET(S,X1,X3,JL)
75 CALL ARORIE(S,X1,X3,JL)
76 IF(HIDAT(2).GT.0.0) THEN
77 PTG1=SQRT(P(JL,1)**2+P(JL,2)**2)
78 PTG2=SQRT(P(JL+1,1)**2+P(JL+1,2)**2)
79 PTG3=SQRT(P(JL+2,1)**2+P(JL+2,2)**2)
80 PTG=MAX(PTG1,PTG2,PTG3)
81 IF(PTG.GT.HIDAT(2)) THEN
82 FMFACT=EXP(-(PTG**2-HIDAT(2)**2)/HIPR1(2)**2)
83 IF(RLU_HIJING(0).GT.FMFACT) GO TO 1
86 C.....ROTATE AND BOOST BACK.....
89 CALL ATROBO(THETA,PHI,-BEX,-BEY,-BEZ,IMIN,IMAX,IERROR)
90 IF(IERROR.NE.0) RETURN
91 C.....ENUMERATE THE GLUONS.....
104 C----THETA FUNCTION DAMPING OF THE EMITTED GLUONS. FOR HADRON-HADRON.
106 C IF(VFR(2).GT.0.) THEN
107 C PTG=SQRT(P(JL+1,1)**2+P(JL+1,2)**2)
109 C DOPT=SQRT((4.*PAR(71)*VFR(2))/WSTRI)
110 C PTOPT=(DOPT*WSTRI)/(2.*VFR(2))
111 C IF(PTG.GT.PTOPT) IORDER=IORDER-1
112 C IF(PTG.GT.PTOPT) GOTO 1
115 IF(SM.GE.HIPR1(5)) GOTO 40