PYTHIA/jetset/luchge.F

   1
   2 C*********************************************************************
   3
   4       FUNCTION LUCHGE(KF)
   5
   6 C...Purpose: to give three times the charge for a particle/parton.
   7       COMMON/LUDAT2/KCHG(500,3),PMAS(500,4),PARF(2000),VCKM(4,4)
   8       SAVE /LUDAT2/
   9
  10 C...Initial values. Simple case of direct readout.
  11       LUCHGE=0
  12       KFA=IABS(KF)
  13       KC=LUCOMP(KFA)
  14       IF(KC.EQ.0) THEN
  15       ELSEIF(KFA.LE.100.OR.KC.LE.80.OR.KC.GT.100) THEN
  16         LUCHGE=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=(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=KCHG(MOD(KFA/1000,10),1)+KCHG(MOD(KFA/100,10),1)
  24       ELSE
  25         LUCHGE=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=LUCHGE*ISIGN(1,KF)
  31
  32       RETURN
  33       END