cluster information

[u/mrichter/AliRoot.git] / HIJING / hipyset1_35 / luchge_hijing.F

* $Id$
    
C*********************************************************************  
    
      FUNCTION LUCHGE_HIJING(KF)   
    
C...Purpose: to give three times the charge for a particle/parton.  
#include "ludat2_hijing.inc"
    
C...Initial values. Simple case of direct readout.  
      LUCHGE_HIJING=0  
      KFA=IABS(KF)  
      KC=LUCOMP_HIJING(KFA)    
      IF(KC.EQ.0) THEN  
      ELSEIF(KFA.LE.100.OR.KC.LE.80.OR.KC.GT.100) THEN  
        LUCHGE_HIJING=KCHG(KC,1)   
    
C...Construction from quark content for heavy meson, diquark, baryon.   
      ELSEIF(MOD(KFA/1000,10).EQ.0) THEN    
        LUCHGE_HIJING=(KCHG(MOD(KFA/100,10),1)-KCHG(MOD(KFA/10,10),1))*    
     &  (-1)**MOD(KFA/100,10)   
      ELSEIF(MOD(KFA/10,10).EQ.0) THEN  
        LUCHGE_HIJING=KCHG(MOD(KFA/1000,10),1)+KCHG(MOD(KFA/100,10),1) 
      ELSE  
        LUCHGE_HIJING=KCHG(MOD(KFA/1000,10),1)+KCHG(MOD(KFA/100,10),1)+    
     &  KCHG(MOD(KFA/10,10),1)  
      ENDIF 
    
C...Add on correct sign.    
      LUCHGE_HIJING=LUCHGE_HIJING*ISIGN(1,KF) 
    
      RETURN    
      END