HIJING/hipyset1_35/luname_hijing.F

   1 * $Id$
   2
   3 C*********************************************************************
   4
   5       SUBROUTINE LUNAME_HIJING(KF,CHAU)
   6
   7 C...Purpose: to give the particle/parton name as a character string.
   8 #include "ludat1_hijing.inc"
   9 #include "ludat2_hijing.inc"
  10 #include "ludat4_hijing.inc"
  11       CHARACTER CHAU*16
  12
  13 C...Initial values. Charge. Subdivide code.
  14       CHAU=' '
  15       KFA=IABS(KF)
  16       KC=LUCOMP_HIJING(KF)
  17       IF(KC.EQ.0) RETURN
  18       KQ=LUCHGE_HIJING(KF)
  19       KFLA=MOD(KFA/1000,10)
  20       KFLB=MOD(KFA/100,10)
  21       KFLC=MOD(KFA/10,10)
  22       KFLS=MOD(KFA,10)
  23       KFLR=MOD(KFA/10000,10)
  24
  25 C...Read out root name and spin for simple particle.
  26       IF(KFA.LE.100.OR.(KFA.GT.100.AND.KC.GT.100)) THEN
  27         CHAU=CHAF(KC)
  28         LEN=0
  29         DO 100 LEM=1,8
  30   100   IF(CHAU(LEM:LEM).NE.' ') LEN=LEM
  31
  32 C...Construct root name for diquark. Add on spin.
  33       ELSEIF(KFLC.EQ.0) THEN
  34         CHAU(1:2)=CHAF(KFLA)(1:1)//CHAF(KFLB)(1:1)
  35         IF(KFLS.EQ.1) CHAU(3:4)='_0'
  36         IF(KFLS.EQ.3) CHAU(3:4)='_1'
  37         LEN=4
  38
  39 C...Construct root name for heavy meson. Add on spin and heavy flavour.
  40       ELSEIF(KFLA.EQ.0) THEN
  41         IF(KFLB.EQ.5) CHAU(1:1)='B'
  42         IF(KFLB.EQ.6) CHAU(1:1)='T'
  43         IF(KFLB.EQ.7) CHAU(1:1)='L'
  44         IF(KFLB.EQ.8) CHAU(1:1)='H'
  45         LEN=1
  46         IF(KFLR.EQ.0.AND.KFLS.EQ.1) THEN
  47         ELSEIF(KFLR.EQ.0.AND.KFLS.EQ.3) THEN
  48           CHAU(2:2)='*'
  49           LEN=2
  50         ELSEIF(KFLR.EQ.1.AND.KFLS.EQ.3) THEN
  51           CHAU(2:3)='_1'
  52           LEN=3
  53         ELSEIF(KFLR.EQ.1.AND.KFLS.EQ.1) THEN
  54           CHAU(2:4)='*_0'
  55           LEN=4
  56         ELSEIF(KFLR.EQ.2) THEN
  57           CHAU(2:4)='*_1'
  58           LEN=4
  59         ELSEIF(KFLS.EQ.5) THEN
  60           CHAU(2:4)='*_2'
  61           LEN=4
  62         ENDIF
  63         IF(KFLC.GE.3.AND.KFLR.EQ.0.AND.KFLS.LE.3) THEN
  64           CHAU(LEN+1:LEN+2)='_'//CHAF(KFLC)(1:1)
  65           LEN=LEN+2
  66         ELSEIF(KFLC.GE.3) THEN
  67           CHAU(LEN+1:LEN+1)=CHAF(KFLC)(1:1)
  68           LEN=LEN+1
  69         ENDIF
  70
  71 C...Construct root name and spin for heavy baryon.
  72       ELSE
  73         IF(KFLB.LE.2.AND.KFLC.LE.2) THEN
  74           CHAU='Sigma '
  75           IF(KFLC.GT.KFLB) CHAU='Lambda'
  76           IF(KFLS.EQ.4) CHAU='Sigma*'
  77           LEN=5
  78           IF(CHAU(6:6).NE.' ') LEN=6
  79         ELSEIF(KFLB.LE.2.OR.KFLC.LE.2) THEN
  80           CHAU='Xi '
  81           IF(KFLA.GT.KFLB.AND.KFLB.GT.KFLC) CHAU='Xi'''
  82           IF(KFLS.EQ.4) CHAU='Xi*'
  83           LEN=2
  84           IF(CHAU(3:3).NE.' ') LEN=3
  85         ELSE
  86           CHAU='Omega '
  87           IF(KFLA.GT.KFLB.AND.KFLB.GT.KFLC) CHAU='Omega'''
  88           IF(KFLS.EQ.4) CHAU='Omega*'
  89           LEN=5
  90           IF(CHAU(6:6).NE.' ') LEN=6
  91         ENDIF
  92
  93 C...Add on heavy flavour content for heavy baryon.
  94         CHAU(LEN+1:LEN+2)='_'//CHAF(KFLA)(1:1)
  95         LEN=LEN+2
  96         IF(KFLB.GE.KFLC.AND.KFLC.GE.4) THEN
  97           CHAU(LEN+1:LEN+2)=CHAF(KFLB)(1:1)//CHAF(KFLC)(1:1)
  98           LEN=LEN+2
  99         ELSEIF(KFLB.GE.KFLC.AND.KFLB.GE.4) THEN
 100           CHAU(LEN+1:LEN+1)=CHAF(KFLB)(1:1)
 101           LEN=LEN+1
 102         ELSEIF(KFLC.GT.KFLB.AND.KFLB.GE.4) THEN
 103           CHAU(LEN+1:LEN+2)=CHAF(KFLC)(1:1)//CHAF(KFLB)(1:1)
 104           LEN=LEN+2
 105         ELSEIF(KFLC.GT.KFLB.AND.KFLC.GE.4) THEN
 106           CHAU(LEN+1:LEN+1)=CHAF(KFLC)(1:1)
 107           LEN=LEN+1
 108         ENDIF
 109       ENDIF
 110
 111 C...Add on bar sign for antiparticle (where necessary).
 112       IF(KF.GT.0.OR.LEN.EQ.0) THEN
 113       ELSEIF(KFA.GT.10.AND.KFA.LE.40.AND.KQ.NE.0) THEN
 114       ELSEIF(KFA.EQ.89.OR.(KFA.GE.91.AND.KFA.LE.99)) THEN
 115       ELSEIF(KFA.GT.100.AND.KFLA.EQ.0.AND.KQ.NE.0) THEN
 116       ELSEIF(MSTU(15).LE.1) THEN
 117         CHAU(LEN+1:LEN+1)='~'
 118         LEN=LEN+1
 119       ELSE
 120         CHAU(LEN+1:LEN+3)='bar'
 121         LEN=LEN+3
 122       ENDIF
 123
 124 C...Add on charge where applicable (conventional cases skipped).
 125       IF(KQ.EQ.6) CHAU(LEN+1:LEN+2)='++'
 126       IF(KQ.EQ.-6) CHAU(LEN+1:LEN+2)='--'
 127       IF(KQ.EQ.3) CHAU(LEN+1:LEN+1)='+'
 128       IF(KQ.EQ.-3) CHAU(LEN+1:LEN+1)='-'
 129       IF(KQ.EQ.0.AND.(KFA.LE.22.OR.LEN.EQ.0)) THEN
 130       ELSEIF(KQ.EQ.0.AND.(KFA.GE.81.AND.KFA.LE.100)) THEN
 131       ELSEIF(KFA.GT.100.AND.KFLA.EQ.0.AND.KFLB.EQ.KFLC.AND.
 132      &KFLB.NE.1) THEN
 133       ELSEIF(KQ.EQ.0) THEN
 134         CHAU(LEN+1:LEN+1)='0'
 135       ENDIF
 136
 137       RETURN
 138       END