hijing/src/hijini.f

0001 C
0002 C
0003 C
0004         SUBROUTINE HIJINI
0005         COMMON/HIPARNT/HIPR1(100),IHPR2(50),HINT1(100),IHNT2(50)
0006         SAVE  /HIPARNT/
0007         COMMON/HISTRNG/NFP(300,15),PP(300,15),NFT(300,15),PT(300,15)
0008         SAVE  /HISTRNG/
0009         COMMON/HIJJET1/NPJ(300),KFPJ(300,500),PJPX(300,500),
0010      &                PJPY(300,500),PJPZ(300,500),PJPE(300,500),
0011      &                PJPM(300,500),NTJ(300),KFTJ(300,500),
0012      &                PJTX(300,500),PJTY(300,500),PJTZ(300,500),
0013      &                PJTE(300,500),PJTM(300,500)
0014         SAVE  /HIJJET1/
0015         COMMON/HIJJET2/NSG,NJSG(900),IASG(900,3),K1SG(900,100),
0016      &          K2SG(900,100),PXSG(900,100),PYSG(900,100),
0017      &          PZSG(900,100),PESG(900,100),PMSG(900,100)
0018         SAVE  /HIJJET2/
0019         COMMON/HIJJET4/NDR,IADR(900,2),KFDR(900),PDR(900,5)
0020         SAVE  /HIJJET4/
0021         COMMON/RANSEED/NSEED
0022         SAVE  /RANSEED/
0023 C****************Reset the momentum of initial particles************
0024 C             and assign flavors to the proj and targ string       *
0025 C*******************************************************************
0026         NSG=0
0027         NDR=0
0028         IPP=2212
0029         IPT=2212
0030         IF(IHNT2(5).NE.0) IPP=IHNT2(5)
0031         IF(IHNT2(6).NE.0) IPT=IHNT2(6)
0032 C               ********in case the proj or targ is a hadron.
0033 C
0034         DO 100 I=1,IHNT2(1)
0035         PP(I,1)=0.0
0036         PP(I,2)=0.0
0037         PP(I,3)=SQRT(HINT1(1)**2/4.0-HINT1(8)**2)
0038         PP(I,4)=HINT1(1)/2
0039         PP(I,5)=HINT1(8)
0040         PP(I,6)=0.0
0041         PP(I,7)=0.0
0042         PP(I,8)=0.0
0043         PP(I,9)=0.0
0044         PP(I,10)=0.0
0045         NFP(I,3)=IPP
0046         NFP(I,4)=IPP
0047         NFP(I,5)=0
0048         NFP(I,6)=0
0049         NFP(I,7)=0
0050         NFP(I,8)=0
0051         NFP(I,9)=0
0052         NFP(I,10)=0
0053         NFP(I,11)=0
0054         NPJ(I)=0
0055         IF(I.GT.ABS(IHNT2(2))) NFP(I,3)=2112
0056         CALL ATTFLV(NFP(I,3),IDQ,IDQQ)
0057         NFP(I,1)=IDQ
0058         NFP(I,2)=IDQQ
0059         NFP(I,15)=-1
0060         IF(ABS(IDQ).GT.1000.OR.(ABS(IDQ*IDQQ).LT.100.AND.
0061      &          ATL_RAN(NSEED).LT.0.5)) NFP(I,15)=1
0062         PP(I,14)=ULMASS(IDQ)
0063         PP(I,15)=ULMASS(IDQQ)
0064 100     CONTINUE
0065 C
0066         DO 200 I=1,IHNT2(3)
0067         PT(I,1)=0.0
0068         PT(I,2)=0.0
0069         PT(I,3)=-SQRT(HINT1(1)**2/4.0-HINT1(9)**2)
0070         PT(I,4)=HINT1(1)/2.0
0071         PT(I,5)=HINT1(9)
0072         PT(I,6)=0.0
0073         PT(I,7)=0.0
0074         PT(I,8)=0.0
0075         PT(I,9)=0.0
0076         PT(I,10)=0.0
0077         NFT(I,3)=IPT
0078         NFT(I,4)=IPT
0079         NFT(I,5)=0
0080         NFT(I,6)=0
0081         NFT(I,7)=0
0082         NFT(I,8)=0
0083         NFT(I,9)=0
0084         NFT(I,10)=0
0085         NFT(I,11)=0
0086         NTJ(I)=0
0087         IF(I.GT.ABS(IHNT2(4))) NFT(I,3)=2112
0088         CALL ATTFLV(NFT(I,3),IDQ,IDQQ)
0089         NFT(I,1)=IDQ
0090         NFT(I,2)=IDQQ
0091         NFT(I,15)=1
0092         IF(ABS(IDQ).GT.1000.OR.(ABS(IDQ*IDQQ).LT.100.AND.
0093      &                  ATL_RAN(NSEED).LT.0.5)) NFT(I,15)=-1
0094         PT(I,14)=ULMASS(IDQ)
0095         PT(I,15)=ULMASS(IDQQ)
0096 200     CONTINUE
0097         RETURN
0098         END