]>
Commit | Line | Data |
---|---|---|
e74335a4 | 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 |