]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 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 |