HIJING/hipyset1_35/lu1ent_hijing.F

   1 * $Id$
   2 C
   3 C
   4 C
   5 C     Modified for HIJING program
   6 c
   7 c    modification July 22, 1997  In pyremnn put an upper limit
   8 c     on the total pt kick the parton can accumulate via multiple
   9 C     scattering. Set the upper limit to be the sqrt(s)/2,
  10 c     this is fix cronin bug for Pb+Pb events at SPS energy.
  11 c
  12 C
  13 C Last modification Oct. 1993 to comply with non-vax
  14 C machines' compiler
  15 C
  16 C
  17       SUBROUTINE LU1ENT_HIJING(IP,KF,PE,THE,PHI)
  18
  19 C...Purpose: to store one parton/particle in commonblock LUJETS.
  20 #include "lujets_hijing.inc"
  21 #include "ludat1_hijing.inc"
  22 #include "ludat2_hijing.inc"
  23
  24 C...Standard checks.
  25       MSTU(28)=0
  26       IF(MSTU(12).GE.1) CALL LULIST_HIJING(0)
  27       IPA=MAX(1,IABS(IP))
  28       IF(IPA.GT.MSTU(4)) CALL LUERRM_HIJING(21,
  29      &'(LU1ENT_HIJING:) writing outside LUJETS_HIJING memory')
  30       KC=LUCOMP_HIJING(KF)
  31       IF(KC.EQ.0) CALL LUERRM_HIJING(12
  32      $     ,'(LU1ENT_HIJING:) unknown flavour code')
  33
  34 C...Find mass. Reset K, P and V vectors.
  35       PM=0.
  36       IF(MSTU(10).EQ.1) PM=P(IPA,5)
  37       IF(MSTU(10).GE.2) PM=ULMASS_HIJING(KF)
  38       DO 100 J=1,5
  39       K(IPA,J)=0
  40       P(IPA,J)=0.
  41   100 V(IPA,J)=0.
  42
  43 C...Store parton/particle in K and P vectors.
  44       K(IPA,1)=1
  45       IF(IP.LT.0) K(IPA,1)=2
  46       K(IPA,2)=KF
  47       P(IPA,5)=PM
  48       P(IPA,4)=MAX(PE,PM)
  49       PA=SQRT(P(IPA,4)**2-P(IPA,5)**2)
  50       P(IPA,1)=PA*SIN(THE)*COS(PHI)
  51       P(IPA,2)=PA*SIN(THE)*SIN(PHI)
  52       P(IPA,3)=PA*COS(THE)
  53
  54 C...Set N. Optionally fragment/decay.
  55       N=IPA
  56       IF(IP.EQ.0) CALL LUEXEC_HIJING
  57
  58       RETURN
  59       END