Merging the VirtualMC branch to the main development branch (HEAD)
[u/mrichter/AliRoot.git] / 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