HIJING/hijing1_36/hijini.F

   1 * $Id$
   2 C
   3 C
   4 C
   5         SUBROUTINE HIJINI
   6 #include "hiparnt.inc"
   7 #include "histrng.inc"
   8 #include "hijjet1.inc"
   9 #include "hijjet2.inc"
  10 #include "hijjet4.inc"
  11         SAVE
  12 C****************Reset the momentum of initial particles************
  13 C             and assign flavors to the proj and targ string       *
  14 C*******************************************************************
  15         NSG=0
  16         NDR=0
  17         IPP=2212
  18         IPT=2212
  19         IF(IHNT2(5).NE.0) IPP=IHNT2(5)
  20         IF(IHNT2(6).NE.0) IPT=IHNT2(6)
  21 C               ********in case the proj or targ is a hadron.
  22 C
  23         DO 100 I=1,IHNT2(1)
  24         PP(I,1)=0.0
  25         PP(I,2)=0.0
  26         PP(I,3)=SQRT(HINT1(1)**2/4.0-HINT1(8)**2)
  27         PP(I,4)=HINT1(1)/2
  28         PP(I,5)=HINT1(8)
  29         PP(I,6)=0.0
  30         PP(I,7)=0.0
  31         PP(I,8)=0.0
  32         PP(I,9)=0.0
  33         PP(I,10)=0.0
  34         NFP(I,3)=IPP
  35         NFP(I,4)=IPP
  36         NFP(I,5)=0
  37         NFP(I,6)=0
  38         NFP(I,7)=0
  39         NFP(I,8)=0
  40         NFP(I,9)=0
  41         NFP(I,10)=0
  42         NFP(I,11)=0
  43         NPJ(I)=0
  44         IF(I.GT.ABS(IHNT2(2))) NFP(I,3)=2112
  45         CALL ATTFLV(NFP(I,3),IDQ,IDQQ)
  46         NFP(I,1)=IDQ
  47         NFP(I,2)=IDQQ
  48         NFP(I,15)=-1
  49         IF(ABS(IDQ).GT.1000.OR.(ABS(IDQ*IDQQ).LT.100.AND.
  50      &          RLU_HIJING(0).LT.0.5)) NFP(I,15)=1
  51         PP(I,14)=ULMASS_HIJING(IDQ)
  52         PP(I,15)=ULMASS_HIJING(IDQQ)
  53 100     CONTINUE
  54 C
  55         DO 200 I=1,IHNT2(3)
  56         PT(I,1)=0.0
  57         PT(I,2)=0.0
  58         PT(I,3)=-SQRT(HINT1(1)**2/4.0-HINT1(9)**2)
  59         PT(I,4)=HINT1(1)/2.0
  60         PT(I,5)=HINT1(9)
  61         PT(I,6)=0.0
  62         PT(I,7)=0.0
  63         PT(I,8)=0.0
  64         PT(I,9)=0.0
  65         PT(I,10)=0.0
  66         NFT(I,3)=IPT
  67         NFT(I,4)=IPT
  68         NFT(I,5)=0
  69         NFT(I,6)=0
  70         NFT(I,7)=0
  71         NFT(I,8)=0
  72         NFT(I,9)=0
  73         NFT(I,10)=0
  74         NFT(I,11)=0
  75         NTJ(I)=0
  76         IF(I.GT.ABS(IHNT2(4))) NFT(I,3)=2112
  77         CALL ATTFLV(NFT(I,3),IDQ,IDQQ)
  78         NFT(I,1)=IDQ
  79         NFT(I,2)=IDQQ
  80         NFT(I,15)=1
  81         IF(ABS(IDQ).GT.1000.OR.(ABS(IDQ*IDQQ).LT.100.AND.
  82      &                  RLU_HIJING(0).LT.0.5)) NFT(I,15)=-1
  83         PT(I,14)=ULMASS_HIJING(IDQ)
  84         PT(I,15)=ULMASS_HIJING(IDQQ)
  85 200     CONTINUE
  86         RETURN
  87         END