luname.f

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

    
C*********************************************************************  
    
      SUBROUTINE luname(KF,CHAU)    
    
C...Purpose: to give the particle/parton name as a character string.    
      common/ludat1/mstu(200),paru(200),mstj(200),parj(200) 
      SAVE /ludat1/ 
      common/ludat2/kchg(500,3),pmas(500,4),parf(2000),vckm(4,4)    
      SAVE /ludat2/ 
      common/ludat4/chaf(500)   
      CHARACTER chaf*8  
      SAVE /ludat4/ 
      CHARACTER chau*16 
    
C...Initial values. Charge. Subdivide code. 
      chau=' '  
      kfa=iabs(kf)  
      kc=lucomp(kf) 
      IF(kc.EQ.0) RETURN    
      kq=luchge(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...Read out root name and spin for simple particle.    
      IF(kfa.LE.100.OR.(kfa.GT.100.AND.kc.GT.100)) THEN 
        chau=chaf(kc)   
        len=0   
        DO 100 lem=1,8  
  100   IF(chau(lem:lem).NE.' ') len=lem    
    
C...Construct root name for diquark. Add on spin.   
      ELSEIF(kflc.EQ.0) THEN    
        chau(1:2)=chaf(kfla)(1:1)//chaf(kflb)(1:1)  
        IF(kfls.EQ.1) chau(3:4)='_0'    
        IF(kfls.EQ.3) chau(3:4)='_1'    
        len=4   
    
C...Construct root name for heavy meson. Add on spin and heavy flavour. 
      ELSEIF(kfla.EQ.0) THEN    
        IF(kflb.EQ.5) chau(1:1)='B' 
        IF(kflb.EQ.6) chau(1:1)='T' 
        IF(kflb.EQ.7) chau(1:1)='L' 
        IF(kflb.EQ.8) chau(1:1)='H' 
        len=1   
        IF(kflr.EQ.0.AND.kfls.EQ.1) THEN    
        ELSEIF(kflr.EQ.0.AND.kfls.EQ.3) THEN    
          chau(2:2)='*' 
          len=2 
        ELSEIF(kflr.EQ.1.AND.kfls.EQ.3) THEN    
          chau(2:3)='_1'    
          len=3 
        ELSEIF(kflr.EQ.1.AND.kfls.EQ.1) THEN    
          chau(2:4)='*_0'   
          len=4 
        ELSEIF(kflr.EQ.2) THEN  
          chau(2:4)='*_1'   
          len=4 
        ELSEIF(kfls.EQ.5) THEN  
          chau(2:4)='*_2'   
          len=4 
        ENDIF   
        IF(kflc.GE.3.AND.kflr.EQ.0.AND.kfls.LE.3) THEN  
          chau(len+1:len+2)='_'//chaf(kflc)(1:1)    
          len=len+2 
        ELSEIF(kflc.GE.3) THEN  
          chau(len+1:len+1)=chaf(kflc)(1:1) 
          len=len+1 
        ENDIF   
    
C...Construct root name and spin for heavy baryon.  
      ELSE  
        IF(kflb.LE.2.AND.kflc.LE.2) THEN    
          chau='Sigma ' 
          IF(kflc.GT.kflb) chau='Lambda'    
          IF(kfls.EQ.4) chau='Sigma*'   
          len=5 
          IF(chau(6:6).NE.' ') len=6    
        ELSEIF(kflb.LE.2.OR.kflc.LE.2) THEN 
          chau='Xi '    
          IF(kfla.GT.kflb.AND.kflb.GT.kflc) chau='Xi''' 
          IF(kfls.EQ.4) chau='Xi*'  
          len=2 
          IF(chau(3:3).NE.' ') len=3    
        ELSE    
          chau='Omega ' 
          IF(kfla.GT.kflb.AND.kflb.GT.kflc) chau='Omega'''  
          IF(kfls.EQ.4) chau='Omega*'   
          len=5 
          IF(chau(6:6).NE.' ') len=6    
        ENDIF   
    
C...Add on heavy flavour content for heavy baryon.  
        chau(len+1:len+2)='_'//chaf(kfla)(1:1)  
        len=len+2   
        IF(kflb.GE.kflc.AND.kflc.GE.4) THEN 
          chau(len+1:len+2)=chaf(kflb)(1:1)//chaf(kflc)(1:1)    
          len=len+2 
        ELSEIF(kflb.GE.kflc.AND.kflb.GE.4) THEN 
          chau(len+1:len+1)=chaf(kflb)(1:1) 
          len=len+1 
        ELSEIF(kflc.GT.kflb.AND.kflb.GE.4) THEN 
          chau(len+1:len+2)=chaf(kflc)(1:1)//chaf(kflb)(1:1)    
          len=len+2 
        ELSEIF(kflc.GT.kflb.AND.kflc.GE.4) THEN 
          chau(len+1:len+1)=chaf(kflc)(1:1) 
          len=len+1 
        ENDIF   
      ENDIF 
    
C...Add on bar sign for antiparticle (where necessary). 
      IF(kf.GT.0.OR.len.EQ.0) THEN  
      ELSEIF(kfa.GT.10.AND.kfa.LE.40.AND.kq.NE.0) THEN  
      ELSEIF(kfa.EQ.89.OR.(kfa.GE.91.AND.kfa.LE.99)) THEN   
      ELSEIF(kfa.GT.100.AND.kfla.EQ.0.AND.kq.NE.0) THEN 
      ELSEIF(mstu(15).LE.1) THEN    
        chau(len+1:len+1)='~'   
        len=len+1   
      ELSE  
        chau(len+1:len+3)='bar' 
        len=len+3   
      ENDIF 
    
C...Add on charge where applicable (conventional cases skipped).    
      IF(kq.EQ.6) chau(len+1:len+2)='++'    
      IF(kq.EQ.-6) chau(len+1:len+2)='--'   
      IF(kq.EQ.3) chau(len+1:len+1)='+' 
      IF(kq.EQ.-3) chau(len+1:len+1)='-'    
      IF(kq.EQ.0.AND.(kfa.LE.22.OR.len.EQ.0)) THEN  
      ELSEIF(kq.EQ.0.AND.(kfa.GE.81.AND.kfa.LE.100)) THEN   
      ELSEIF(kfa.GT.100.AND.kfla.EQ.0.AND.kflb.EQ.kflc.AND. 
     &kflb.NE.1) THEN   
      ELSEIF(kq.EQ.0) THEN  
        chau(len+1:len+1)='0'   
      ENDIF 
    
      RETURN    
      END