]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 1 | *CMZ : 17/07/98 15.44.33 by Federico Carminati |
2 | *-- Author : | |
3 | C********************************************************************* | |
4 | ||
5 | FUNCTION LUCHGE(KF) | |
6 | ||
7 | C...Purpose: to give three times the charge for a particle/parton. | |
8 | *KEEP,LUDAT2. | |
9 | COMMON /LUDAT2/ KCHG(500,3),PMAS(500,4),PARF(2000),VCKM(4,4) | |
10 | SAVE /LUDAT2/ | |
11 | *KEND. | |
12 | ||
13 | C...Initial values. Simple case of direct readout. | |
14 | LUCHGE=0 | |
15 | KFA=IABS(KF) | |
16 | KC=LUCOMP(KFA) | |
17 | IF(KC.EQ.0) THEN | |
18 | ELSEIF(KFA.LE.100.OR.KC.LE.80.OR.KC.GT.100) THEN | |
19 | LUCHGE=KCHG(KC,1) | |
20 | ||
21 | C...Construction from quark content for heavy meson, diquark, baryon. | |
22 | ELSEIF(MOD(KFA/1000,10).EQ.0) THEN | |
23 | LUCHGE=(KCHG(MOD(KFA/100,10),1)-KCHG(MOD(KFA/10,10),1))* | |
24 | & (-1)**MOD(KFA/100,10) | |
25 | ELSEIF(MOD(KFA/10,10).EQ.0) THEN | |
26 | LUCHGE=KCHG(MOD(KFA/1000,10),1)+KCHG(MOD(KFA/100,10),1) | |
27 | ELSE | |
28 | LUCHGE=KCHG(MOD(KFA/1000,10),1)+KCHG(MOD(KFA/100,10),1)+ | |
29 | & KCHG(MOD(KFA/10,10),1) | |
30 | ENDIF | |
31 | ||
32 | C...Add on correct sign. | |
33 | LUCHGE=LUCHGE*ISIGN(1,KF) | |
34 | ||
35 | RETURN | |
36 | END |