2 C*********************************************************************
4 SUBROUTINE LUNAME(KF,CHAU)
6 C...Purpose: to give the particle/parton name as a character string.
7 COMMON/LUDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
8 COMMON/LUDAT2/KCHG(500,3),PMAS(500,4),PARF(2000),VCKM(4,4)
9 COMMON/LUDAT4/CHAF(500)
11 SAVE /LUDAT1/,/LUDAT2/,/LUDAT4/
14 C...Initial values. Charge. Subdivide code.
24 KFLR=MOD(KFA/10000,10)
26 C...Read out root name and spin for simple particle.
27 IF(KFA.LE.100.OR.(KFA.GT.100.AND.KC.GT.100)) THEN
31 IF(CHAU(LEM:LEM).NE.' ') LEN=LEM
34 C...Construct root name for diquark. Add on spin.
35 ELSEIF(KFLC.EQ.0) THEN
36 CHAU(1:2)=CHAF(KFLA)(1:1)//CHAF(KFLB)(1:1)
37 IF(KFLS.EQ.1) CHAU(3:4)='_0'
38 IF(KFLS.EQ.3) CHAU(3:4)='_1'
41 C...Construct root name for heavy meson. Add on spin and heavy flavour.
42 ELSEIF(KFLA.EQ.0) THEN
43 IF(KFLB.EQ.5) CHAU(1:1)='B'
44 IF(KFLB.EQ.6) CHAU(1:1)='T'
45 IF(KFLB.EQ.7) CHAU(1:1)='L'
46 IF(KFLB.EQ.8) CHAU(1:1)='H'
48 IF(KFLR.EQ.0.AND.KFLS.EQ.1) THEN
49 ELSEIF(KFLR.EQ.0.AND.KFLS.EQ.3) THEN
52 ELSEIF(KFLR.EQ.1.AND.KFLS.EQ.3) THEN
55 ELSEIF(KFLR.EQ.1.AND.KFLS.EQ.1) THEN
58 ELSEIF(KFLR.EQ.2) THEN
61 ELSEIF(KFLS.EQ.5) THEN
65 IF(KFLC.GE.3.AND.KFLR.EQ.0.AND.KFLS.LE.3) THEN
66 CHAU(LEN+1:LEN+2)='_'//CHAF(KFLC)(1:1)
68 ELSEIF(KFLC.GE.3) THEN
69 CHAU(LEN+1:LEN+1)=CHAF(KFLC)(1:1)
73 C...Construct root name and spin for heavy baryon.
75 IF(KFLB.LE.2.AND.KFLC.LE.2) THEN
77 IF(KFLC.GT.KFLB) CHAU='Lambda'
78 IF(KFLS.EQ.4) CHAU='Sigma*'
80 IF(CHAU(6:6).NE.' ') LEN=6
81 ELSEIF(KFLB.LE.2.OR.KFLC.LE.2) THEN
83 IF(KFLA.GT.KFLB.AND.KFLB.GT.KFLC) CHAU='Xi'''
84 IF(KFLS.EQ.4) CHAU='Xi*'
86 IF(CHAU(3:3).NE.' ') LEN=3
89 IF(KFLA.GT.KFLB.AND.KFLB.GT.KFLC) CHAU='Omega'''
90 IF(KFLS.EQ.4) CHAU='Omega*'
92 IF(CHAU(6:6).NE.' ') LEN=6
95 C...Add on heavy flavour content for heavy baryon.
96 CHAU(LEN+1:LEN+2)='_'//CHAF(KFLA)(1:1)
98 IF(KFLB.GE.KFLC.AND.KFLC.GE.4) THEN
99 CHAU(LEN+1:LEN+2)=CHAF(KFLB)(1:1)//CHAF(KFLC)(1:1)
101 ELSEIF(KFLB.GE.KFLC.AND.KFLB.GE.4) THEN
102 CHAU(LEN+1:LEN+1)=CHAF(KFLB)(1:1)
104 ELSEIF(KFLC.GT.KFLB.AND.KFLB.GE.4) THEN
105 CHAU(LEN+1:LEN+2)=CHAF(KFLC)(1:1)//CHAF(KFLB)(1:1)
107 ELSEIF(KFLC.GT.KFLB.AND.KFLC.GE.4) THEN
108 CHAU(LEN+1:LEN+1)=CHAF(KFLC)(1:1)
113 C...Add on bar sign for antiparticle (where necessary).
114 IF(KF.GT.0.OR.LEN.EQ.0) THEN
115 ELSEIF(KFA.GT.10.AND.KFA.LE.40.AND.KQ.NE.0.AND.MOD(KQ,3).EQ.0)
117 ELSEIF(KFA.EQ.89.OR.(KFA.GE.91.AND.KFA.LE.99)) THEN
118 ELSEIF(KFA.GT.100.AND.KFLA.EQ.0.AND.KQ.NE.0) THEN
119 ELSEIF(MSTU(15).LE.1) THEN
120 CHAU(LEN+1:LEN+1)='~'
123 CHAU(LEN+1:LEN+3)='bar'
127 C...Add on charge where applicable (conventional cases skipped).
128 IF(KQ.EQ.6) CHAU(LEN+1:LEN+2)='++'
129 IF(KQ.EQ.-6) CHAU(LEN+1:LEN+2)='--'
130 IF(KQ.EQ.3) CHAU(LEN+1:LEN+1)='+'
131 IF(KQ.EQ.-3) CHAU(LEN+1:LEN+1)='-'
132 IF(KQ.EQ.0.AND.(KFA.LE.22.OR.LEN.EQ.0)) THEN
133 ELSEIF(KQ.EQ.0.AND.(KFA.GE.81.AND.KFA.LE.100)) THEN
134 ELSEIF(KFA.EQ.28.OR.KFA.EQ.29) THEN
135 ELSEIF(KFA.GT.100.AND.KFLA.EQ.0.AND.KFLB.EQ.KFLC.AND.
138 CHAU(LEN+1:LEN+1)='0'