1 #include "isajet/pilot.h"
2 DOUBLE PRECISION FUNCTION SSMQCD(DM,DQ)
3 C-----------------------------------------------------------------------
4 C Calculate leading-log running mass for quark with mass DM at
5 C scale Q, using alpha_s which is continuous across thresholds.
6 C See Drees and Hikasa, Phys. Lett. B240: 455-464, Eq. 4.5.
8 C Note the threshold is at Q = 2 m, not at Q = m as in MSbar.
10 C Bisset's QCDRAD, WDHFFC
11 C-----------------------------------------------------------------------
12 #if defined(CERNLIB_IMPNONE)
15 #include "isajet/sssm.inc"
17 DOUBLE PRECISION DM,DQ,DLAM4,DLAM5,DLAM6,DNEFF,POW,RENORM
20 C Do nothing for light quarks
27 C Calculate running mass
35 POW=12.D0/(33.D0-2.*DNEFF)
37 RENORM=(LOG(2*DM/DLAM4)/LOG(DQ/DLAM4))**POW
41 RENORM=(LOG(2*DM/DLAM4)/LOG(DQBT/DLAM4))**POW
45 POW=12.D0/(33.D0-2.*DNEFF)
46 DLAM5=DEXP((25.D0*LOG(DLAM4)-LOG(DQBT**2))/23.D0)
47 IF(DQ.GE.DQBT.AND.DQ.LT.DQTP) THEN
49 $ *(LOG(DQBT/DLAM5)/LOG(DQ/DLAM5))**POW
54 $ *(LOG(DQBT/DLAM5)/LOG(DQTP/DLAM5))**POW
58 POW=12.D0/(33.D0-2.*DNEFF)
59 DLAM6=DEXP((25.D0*LOG(DLAM4)-LOG(DQBT**2)
60 $-LOG(4*AMTP**2))/21.D0)
62 $*(LOG(DQTP/DLAM6)/LOG(DQ/DLAM6))**POW