PYTHIA/jetset/ulalem.F

   1
   2 C*********************************************************************
   3
   4       FUNCTION ULALEM(Q2)
   5
   6 C...Purpose: to calculate the running alpha_electromagnetic.
   7       COMMON/LUDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
   8       SAVE /LUDAT1/
   9
  10 C...Calculate real part of photon vacuum polarization.
  11 C...For leptons simplify by using asymptotic (Q^2 >> m^2) expressions.
  12 C...For hadrons use parametrization of H. Burkhardt et al.
  13 C...See R. Kleiss et al, CERN 89-08, vol. 3, pp. 129-131.
  14       AEMPI=PARU(101)/(3.*PARU(1))
  15       IF(MSTU(101).LE.0.OR.Q2.LT.2E-6) THEN
  16         RPIGG=0.
  17       ELSEIF(MSTU(101).EQ.2.AND.Q2.LT.PARU(104)) THEN
  18         RPIGG=0.
  19       ELSEIF(MSTU(101).EQ.2) THEN
  20         RPIGG=1.-PARU(101)/PARU(103)
  21       ELSEIF(Q2.LT.0.09) THEN
  22         RPIGG=AEMPI*(13.4916+LOG(Q2))+0.00835*LOG(1.+Q2)
  23       ELSEIF(Q2.LT.9.) THEN
  24         RPIGG=AEMPI*(16.3200+2.*LOG(Q2))+0.00238*LOG(1.+3.927*Q2)
  25       ELSEIF(Q2.LT.1E4) THEN
  26         RPIGG=AEMPI*(13.4955+3.*LOG(Q2))+0.00165+0.00299*LOG(1.+Q2)
  27       ELSE
  28         RPIGG=AEMPI*(13.4955+3.*LOG(Q2))+0.00221+0.00293*LOG(1.+Q2)
  29       ENDIF
  30
  31 C...Calculate running alpha_em.
  32       ULALEM=PARU(101)/(1.-RPIGG)
  33       PARU(108)=ULALEM
  34
  35       RETURN
  36       END