luchge.f

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file luchge.f

    
C*********************************************************************  
    
      FUNCTION luchge(KF)   
    
C...Purpose: to give three times the charge for a particle/parton.  
      common/ludat2/kchg(500,3),pmas(500,4),parf(2000),vckm(4,4)    
      SAVE /ludat2/ 
    
C...Initial values. Simple case of direct readout.  
      luchge=0  
      kfa=iabs(kf)  
      kc=lucomp(kfa)    
      IF(kc.EQ.0) THEN  
      ELSEIF(kfa.LE.100.OR.kc.LE.80.OR.kc.GT.100) THEN  
        luchge=kchg(kc,1)   
    
C...Construction from quark content for heavy meson, diquark, baryon.   
      ELSEIF(mod(kfa/1000,10).EQ.0) THEN    
        luchge=(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=kchg(mod(kfa/1000,10),1)+kchg(mod(kfa/100,10),1) 
      ELSE  
        luchge=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=luchge*isign(1,kf) 
    
      RETURN    
      END