HIJING/hipyset1_35/luchge_hijing.F

   1 * $Id$
   2
   3 C*********************************************************************
   4
   5       FUNCTION LUCHGE_HIJING(KF)
   6
   7 C...Purpose: to give three times the charge for a particle/parton.
   8 #include "ludat2_hijing.inc"
   9
  10 C...Initial values. Simple case of direct readout.
  11       LUCHGE_HIJING=0
  12       KFA=IABS(KF)
  13       KC=LUCOMP_HIJING(KFA)
  14       IF(KC.EQ.0) THEN
  15       ELSEIF(KFA.LE.100.OR.KC.LE.80.OR.KC.GT.100) THEN
  16         LUCHGE_HIJING=KCHG(KC,1)
  17
  18 C...Construction from quark content for heavy meson, diquark, baryon.
  19       ELSEIF(MOD(KFA/1000,10).EQ.0) THEN
  20         LUCHGE_HIJING=(KCHG(MOD(KFA/100,10),1)-KCHG(MOD(KFA/10,10),1))*
  21      &  (-1)**MOD(KFA/100,10)
  22       ELSEIF(MOD(KFA/10,10).EQ.0) THEN
  23         LUCHGE_HIJING=KCHG(MOD(KFA/1000,10),1)+KCHG(MOD(KFA/100,10),1)
  24       ELSE
  25         LUCHGE_HIJING=KCHG(MOD(KFA/1000,10),1)+KCHG(MOD(KFA/100,10),1)+
  26      &  KCHG(MOD(KFA/10,10),1)
  27       ENDIF
  28
  29 C...Add on correct sign.
  30       LUCHGE_HIJING=LUCHGE_HIJING*ISIGN(1,KF)
  31
  32       RETURN
  33       END