pymrun.f

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

 
C*********************************************************************
 
C...PYMRUN
C...Gives the running, current-algebra mass of a d, u, s, c or b quark,
C...for Higgs couplings. Everything else sent on to PYMASS.
 
      FUNCTION pymrun(KF,Q2)
 
C...Double precision and integer declarations.
      IMPLICIT DOUBLE PRECISION(a-h, o-z)
      IMPLICIT INTEGER(i-n)
      INTEGER pyk,pychge,pycomp
C...Commonblocks.
      common/pydat1/mstu(200),paru(200),mstj(200),parj(200)
      common/pydat2/kchg(500,4),pmas(500,4),parf(2000),vckm(4,4)
      common/pypars/mstp(200),parp(200),msti(200),pari(200)
      SAVE /pydat1/,/pydat2/,/pypars/
 
C...Most masses not handled here.
      kfa=iabs(kf)
      IF(kfa.EQ.0.OR.kfa.GT.6) THEN
        pymrun=pymass(kf)
 
C...Current-algebra masses, but no Q2 dependence.
      ELSEIF(mstp(37).NE.1.OR.mstp(2).LE.0) THEN
        pymrun=parf(90+kfa)
 
C...Running current-algebra masses.
      ELSE
        as=pyalps(q2)
        pymrun=parf(90+kfa)*
     &  (log(max(4d0,parp(37)**2*parf(90+kfa)**2/paru(117)**2))/
     &  log(max(4d0,q2/paru(117)**2)))**(12d0/(33d0-2d0*mstu(118)))
      ENDIF
 
      RETURN
      END