pythiaeRHIC/pythia/pymrun.f

0001
0002 C*********************************************************************
0003
0004 C...PYMRUN
0005 C...Gives the running, current-algebra mass of a d, u, s, c or b quark,
0006 C...for Higgs couplings. Everything else sent on to PYMASS.
0007
0008       FUNCTION PYMRUN(KF,Q2)
0009
0010 C...Double precision and integer declarations.
0011       IMPLICIT DOUBLE PRECISION(A-H, O-Z)
0012       IMPLICIT INTEGER(I-N)
0013       INTEGER PYK,PYCHGE,PYCOMP
0014 C...Commonblocks.
0015       COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
0016       COMMON/PYDAT2/KCHG(500,4),PMAS(500,4),PARF(2000),VCKM(4,4)
0017       COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
0018       SAVE /PYDAT1/,/PYDAT2/,/PYPARS/
0019
0020 C...Most masses not handled here.
0021       KFA=IABS(KF)
0022       IF(KFA.EQ.0.OR.KFA.GT.6) THEN
0023         PYMRUN=PYMASS(KF)
0024
0025 C...Current-algebra masses, but no Q2 dependence.
0026       ELSEIF(MSTP(37).NE.1.OR.MSTP(2).LE.0) THEN
0027         PYMRUN=PARF(90+KFA)
0028
0029 C...Running current-algebra masses.
0030       ELSE
0031         AS=PYALPS(Q2)
0032         PYMRUN=PARF(90+KFA)*
0033      &  (LOG(MAX(4D0,PARP(37)**2*PARF(90+KFA)**2/PARU(117)**2))/
0034      &  LOG(MAX(4D0,Q2/PARU(117)**2)))**(12D0/(33D0-2D0*MSTU(118)))
0035       ENDIF
0036
0037       RETURN
0038       END