3 C*********************************************************************
5 SUBROUTINE PYSPLI_HIJING(KF,KFLIN,KFLCH,KFLSP)
7 C...In case of a hadron remnant which is more complicated than just a
8 C...quark or a diquark, split it into two (partons or hadron + parton).
11 C...Preliminaries. Parton composition.
14 KFL(1)=MOD(KFA/1000,10)
15 KFL(2)=MOD(KFA/100,10)
22 KFL(2)=KFL(2)*(-1)**KFL(2)
23 KFL(3)=-KFL(3)*(-1)**IABS(KFL(2))
24 IF(KFLR.EQ.KFL(2)) THEN
26 ELSEIF(KFLR.EQ.KFL(3)) THEN
28 ELSEIF(IABS(KFLR).EQ.21.AND.RLU_HIJING(0).GT.0.5) THEN
31 ELSEIF(IABS(KFLR).EQ.21) THEN
34 ELSEIF(KFLR*KFL(2).GT.0) THEN
35 CALL LUKFDI_HIJING(-KFLR,KFL(2),KFDUMP,KFLCH)
38 CALL LUKFDI_HIJING(-KFLR,KFL(3),KFDUMP,KFLCH)
46 100 IF(KFLR.EQ.KFL(J)) NAGR=NAGR+1
48 RAGR=0.00001+(NAGR-0.00002)*RLU_HIJING(0)
51 IF(KFLR.EQ.KFL(J)) RAGR=RAGR-1.
52 110 IF(IAGR.EQ.0.AND.RAGR.LE.0.) IAGR=J
54 IAGR=1.00001+2.99998*RLU_HIJING(0)
58 IF(IAGR.EQ.1.AND.KFL(3).GT.KFL(2)) ID1=3
61 IF(MOD(KFA,10).EQ.2.AND.KFL(1).EQ.KFL(2)) THEN
62 IF(IAGR.NE.3.AND.RLU_HIJING(0).GT.0.25) KSP=1
63 ELSEIF(MOD(KFA,10).EQ.2.AND.KFL(2).GE.KFL(3)) THEN
64 IF(IAGR.NE.1.AND.RLU_HIJING(0).GT.0.25) KSP=1
65 ELSEIF(MOD(KFA,10).EQ.2) THEN
67 IF(IAGR.NE.1.AND.RLU_HIJING(0).GT.0.75) KSP=1
69 KFLSP=1000*KFL(ID1)+100*KFL(ID2)+KSP
72 ELSEIF(NAGR.EQ.0.AND.KFLR.GT.0) THEN
73 CALL LUKFDI_HIJING(-KFLR,KFL(IAGR),KFDUMP,KFLCH)
74 ELSEIF(NAGR.EQ.0) THEN
75 CALL LUKFDI_HIJING(10000+KFLSP,-KFLR,KFDUMP,KFLCH)
80 C...Add on correct sign for result.