3 C*********************************************************************
5 FUNCTION LUCOMP_HIJING(KF)
7 C...Purpose: to compress the standard KF codes for use in mass and decay
8 C...arrays; also to check whether a given code actually is defined.
9 #include "ludat2_hijing.inc"
11 C...Subdivide KF code into constituent pieces.
18 KFLR=MOD(KFA/10000,10)
20 C...Simple cases: direct translation or special codes.
21 IF(KFA.EQ.0.OR.KFA.GE.100000) THEN
22 ELSEIF(KFA.LE.100) THEN
24 IF(KF.LT.0.AND.KCHG(KFA,3).EQ.0) LUCOMP_HIJING=0
25 ELSEIF(KFLS.EQ.0) THEN
26 IF(KF.EQ.130) LUCOMP_HIJING=221
27 IF(KF.EQ.310) LUCOMP_HIJING=222
28 IF(KFA.EQ.210) LUCOMP_HIJING=281
29 IF(KFA.EQ.2110) LUCOMP_HIJING=282
30 IF(KFA.EQ.2210) LUCOMP_HIJING=283
33 ELSEIF(KFA-10000*KFLR.LT.1000) THEN
34 IF(KFLB.EQ.0.OR.KFLB.EQ.9.OR.KFLC.EQ.0.OR.KFLC.EQ.9) THEN
35 ELSEIF(KFLB.LT.KFLC) THEN
36 ELSEIF(KF.LT.0.AND.KFLB.EQ.KFLC) THEN
37 ELSEIF(KFLB.EQ.KFLC) THEN
38 IF(KFLR.EQ.0.AND.KFLS.EQ.1) THEN
39 LUCOMP_HIJING=110+KFLB
40 ELSEIF(KFLR.EQ.0.AND.KFLS.EQ.3) THEN
41 LUCOMP_HIJING=130+KFLB
42 ELSEIF(KFLR.EQ.1.AND.KFLS.EQ.3) THEN
43 LUCOMP_HIJING=150+KFLB
44 ELSEIF(KFLR.EQ.1.AND.KFLS.EQ.1) THEN
45 LUCOMP_HIJING=170+KFLB
46 ELSEIF(KFLR.EQ.2.AND.KFLS.EQ.3) THEN
47 LUCOMP_HIJING=190+KFLB
48 ELSEIF(KFLR.EQ.0.AND.KFLS.EQ.5) THEN
49 LUCOMP_HIJING=210+KFLB
51 ELSEIF(KFLB.LE.5.AND.KFLC.LE.3) THEN
52 IF(KFLR.EQ.0.AND.KFLS.EQ.1) THEN
53 LUCOMP_HIJING=100+((KFLB-1)*(KFLB-2))/2+KFLC
54 ELSEIF(KFLR.EQ.0.AND.KFLS.EQ.3) THEN
55 LUCOMP_HIJING=120+((KFLB-1)*(KFLB-2))/2+KFLC
56 ELSEIF(KFLR.EQ.1.AND.KFLS.EQ.3) THEN
57 LUCOMP_HIJING=140+((KFLB-1)*(KFLB-2))/2+KFLC
58 ELSEIF(KFLR.EQ.1.AND.KFLS.EQ.1) THEN
59 LUCOMP_HIJING=160+((KFLB-1)*(KFLB-2))/2+KFLC
60 ELSEIF(KFLR.EQ.2.AND.KFLS.EQ.3) THEN
61 LUCOMP_HIJING=180+((KFLB-1)*(KFLB-2))/2+KFLC
62 ELSEIF(KFLR.EQ.0.AND.KFLS.EQ.5) THEN
63 LUCOMP_HIJING=200+((KFLB-1)*(KFLB-2))/2+KFLC
65 ELSEIF((KFLS.EQ.1.AND.KFLR.LE.1).OR.(KFLS.EQ.3.AND.KFLR.LE.2).
66 & OR.(KFLS.EQ.5.AND.KFLR.EQ.0)) THEN
71 ELSEIF((KFLR.EQ.0.OR.KFLR.EQ.1).AND.KFLC.EQ.0) THEN
72 IF(KFLS.NE.1.AND.KFLS.NE.3) THEN
73 ELSEIF(KFLA.EQ.9.OR.KFLB.EQ.0.OR.KFLB.EQ.9) THEN
74 ELSEIF(KFLA.LT.KFLB) THEN
75 ELSEIF(KFLS.EQ.1.AND.KFLA.EQ.KFLB) THEN
81 ELSEIF(KFLR.EQ.0.AND.KFLS.EQ.2) THEN
82 IF(KFLA.EQ.9.OR.KFLB.EQ.0.OR.KFLB.EQ.9.OR.KFLC.EQ.9) THEN
83 ELSEIF(KFLA.LE.KFLC.OR.KFLA.LT.KFLB) THEN
84 ELSEIF(KFLA.GE.6.OR.KFLB.GE.4.OR.KFLC.GE.4) THEN
86 ELSEIF(KFLB.LT.KFLC) THEN
87 LUCOMP_HIJING=300+((KFLA+1)*KFLA*(KFLA-1))/6+(KFLC*(KFLC-1))
90 LUCOMP_HIJING=330+((KFLA+1)*KFLA*(KFLA-1))/6+(KFLB*(KFLB-1))
95 ELSEIF(KFLR.EQ.0.AND.KFLS.EQ.4) THEN
96 IF(KFLA.EQ.9.OR.KFLB.EQ.0.OR.KFLB.EQ.9.OR.KFLC.EQ.9) THEN
97 ELSEIF(KFLA.LT.KFLB.OR.KFLB.LT.KFLC) THEN
98 ELSEIF(KFLA.GE.6.OR.KFLB.GE.4) THEN
101 LUCOMP_HIJING=360+((KFLA+1)*KFLA*(KFLA-1))/6+(KFLB*(KFLB-1))