hijing/src/attflv.f

0001 C
0002 C
0003 C
0004         SUBROUTINE ATTFLV(ID,IDQ,IDQQ)
0005         COMMON/RANSEED/NSEED
0006         SAVE  /RANSEED/
0007 C
0008         IF(ABS(ID).LT.100) THEN
0009                 NSIGN=1
0010                 IDQ=ID/100
0011                 IDQQ=-ID/10+IDQ*10
0012                 IF(ABS(IDQ).EQ.3) NSIGN=-1
0013                 IDQ=NSIGN*IDQ
0014                 IDQQ=NSIGN*IDQQ
0015                 IF(IDQ.LT.0) THEN
0016                         ID0=IDQ
0017                         IDQ=IDQQ
0018                         IDQQ=ID0
0019                 ENDIF
0020                 RETURN
0021         ENDIF
0022 C               ********return ID of quark(IDQ) and anti-quark(IDQQ)
0023 C                       for pions and kaons
0024 c
0025 C       Return LU ID for quarks and diquarks for proton(ID=2212)
0026 C       anti-proton(ID=-2212) and nuetron(ID=2112)
0027 C       LU ID for d=1,u=2, (ud)0=2101, (ud)1=2103,
0028 C       (dd)1=1103,(uu)1=2203.
0029 C       Use SU(6)  weight  proton=1/3d(uu)1 + 1/6u(ud)1 + 1/2u(ud)0
0030 C                         nurtron=1/3u(dd)1 + 1/6d(ud)1 + 1/2d(ud)0
0031 C
0032         IDQ=2
0033         IF(ABS(ID).EQ.2112) IDQ=1
0034         IDQQ=2101
0035         X=ATL_RAN(NSEED)
0036         IF(X.LE.0.5) GO TO 30
0037         IF(X.GT.0.666667) GO TO 10
0038         IDQQ=2103
0039         GO TO 30
0040 10      IDQ=1
0041         IDQQ=2203
0042         IF(ABS(ID).EQ.2112) THEN
0043                 IDQ=2
0044                 IDQQ=1103
0045         ENDIF
0046 30      IF(ID.LT.0) THEN
0047                 ID00=IDQQ
0048                 IDQQ=-IDQ
0049                 IDQ=-ID00
0050         ENDIF
0051         RETURN
0052         END