lucomp.f

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

    
C*********************************************************************  
    
      FUNCTION lucomp(KF)   
    
C...Purpose: to compress the standard KF codes for use in mass and decay    
C...arrays; also to check whether a given code actually is defined. 
      common/ludat2/kchg(500,3),pmas(500,4),parf(2000),vckm(4,4)    
      SAVE /ludat2/ 
    
C...Subdivide KF code into constituent pieces.  
      lucomp=0  
      kfa=iabs(kf)  
      kfla=mod(kfa/1000,10) 
      kflb=mod(kfa/100,10)  
      kflc=mod(kfa/10,10)   
      kfls=mod(kfa,10)  
      kflr=mod(kfa/10000,10)    
    
C...Simple cases: direct translation or special codes.  
      IF(kfa.EQ.0.OR.kfa.GE.100000) THEN    
      ELSEIF(kfa.LE.100) THEN   
        lucomp=kfa  
        IF(kf.LT.0.AND.kchg(kfa,3).EQ.0) lucomp=0   
      ELSEIF(kfls.EQ.0) THEN    
        IF(kf.EQ.130) lucomp=221    
        IF(kf.EQ.310) lucomp=222    
        IF(kfa.EQ.210) lucomp=281   
        IF(kfa.EQ.2110) lucomp=282  
        IF(kfa.EQ.2210) lucomp=283  
    
C...Mesons. 
      ELSEIF(kfa-10000*kflr.LT.1000) THEN   
        IF(kflb.EQ.0.OR.kflb.EQ.9.OR.kflc.EQ.0.OR.kflc.EQ.9) THEN   
        ELSEIF(kflb.LT.kflc) THEN   
        ELSEIF(kf.LT.0.AND.kflb.EQ.kflc) THEN   
        ELSEIF(kflb.EQ.kflc) THEN   
          IF(kflr.EQ.0.AND.kfls.EQ.1) THEN  
            lucomp=110+kflb 
          ELSEIF(kflr.EQ.0.AND.kfls.EQ.3) THEN  
            lucomp=130+kflb 
          ELSEIF(kflr.EQ.1.AND.kfls.EQ.3) THEN  
            lucomp=150+kflb 
          ELSEIF(kflr.EQ.1.AND.kfls.EQ.1) THEN  
            lucomp=170+kflb 
          ELSEIF(kflr.EQ.2.AND.kfls.EQ.3) THEN  
            lucomp=190+kflb 
          ELSEIF(kflr.EQ.0.AND.kfls.EQ.5) THEN  
            lucomp=210+kflb 
          ENDIF 
        ELSEIF(kflb.LE.5.AND.kflc.LE.3) THEN    
          IF(kflr.EQ.0.AND.kfls.EQ.1) THEN  
            lucomp=100+((kflb-1)*(kflb-2))/2+kflc   
          ELSEIF(kflr.EQ.0.AND.kfls.EQ.3) THEN  
            lucomp=120+((kflb-1)*(kflb-2))/2+kflc   
          ELSEIF(kflr.EQ.1.AND.kfls.EQ.3) THEN  
            lucomp=140+((kflb-1)*(kflb-2))/2+kflc   
          ELSEIF(kflr.EQ.1.AND.kfls.EQ.1) THEN  
            lucomp=160+((kflb-1)*(kflb-2))/2+kflc   
          ELSEIF(kflr.EQ.2.AND.kfls.EQ.3) THEN  
            lucomp=180+((kflb-1)*(kflb-2))/2+kflc   
          ELSEIF(kflr.EQ.0.AND.kfls.EQ.5) THEN  
            lucomp=200+((kflb-1)*(kflb-2))/2+kflc   
          ENDIF 
        ELSEIF((kfls.EQ.1.AND.kflr.LE.1).OR.(kfls.EQ.3.AND.kflr.LE.2).  
     &  or.(kfls.EQ.5.AND.kflr.EQ.0)) THEN  
          lucomp=80+kflb    
        ENDIF   
    
C...Diquarks.   
      ELSEIF((kflr.EQ.0.OR.kflr.EQ.1).AND.kflc.EQ.0) THEN   
        IF(kfls.NE.1.AND.kfls.NE.3) THEN    
        ELSEIF(kfla.EQ.9.OR.kflb.EQ.0.OR.kflb.EQ.9) THEN    
        ELSEIF(kfla.LT.kflb) THEN   
        ELSEIF(kfls.EQ.1.AND.kfla.EQ.kflb) THEN 
        ELSE    
          lucomp=90 
        ENDIF   
    
C...Spin 1/2 baryons.   
      ELSEIF(kflr.EQ.0.AND.kfls.EQ.2) THEN  
        IF(kfla.EQ.9.OR.kflb.EQ.0.OR.kflb.EQ.9.OR.kflc.EQ.9) THEN   
        ELSEIF(kfla.LE.kflc.OR.kfla.LT.kflb) THEN   
        ELSEIF(kfla.GE.6.OR.kflb.GE.4.OR.kflc.GE.4) THEN    
          lucomp=80+kfla    
        ELSEIF(kflb.LT.kflc) THEN   
          lucomp=300+((kfla+1)*kfla*(kfla-1))/6+(kflc*(kflc-1))/2+kflb  
        ELSE    
          lucomp=330+((kfla+1)*kfla*(kfla-1))/6+(kflb*(kflb-1))/2+kflc  
        ENDIF   
    
C...Spin 3/2 baryons.   
      ELSEIF(kflr.EQ.0.AND.kfls.EQ.4) THEN  
        IF(kfla.EQ.9.OR.kflb.EQ.0.OR.kflb.EQ.9.OR.kflc.EQ.9) THEN   
        ELSEIF(kfla.LT.kflb.OR.kflb.LT.kflc) THEN   
        ELSEIF(kfla.GE.6.OR.kflb.GE.4) THEN 
          lucomp=80+kfla    
        ELSE    
          lucomp=360+((kfla+1)*kfla*(kfla-1))/6+(kflb*(kflb-1))/2+kflc  
        ENDIF   
      ENDIF 
    
      RETURN    
      END