1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
|
program example3
implicit double precision (a-h,o-z)
double precision f(-6:6)
QMZ=91.71d0
Nmem=31
write(*,*)
do i=1,Nmem
write(*,*) '---------------------------------------------'
write(*,*) 'PDF set ',i
write(*,*)
write(*,*) 'x*Gluon'
write(*,*) ' x Q=10 GeV Q=100 GeV Q=1000 GeV'
do ix=10,95,10
x = dfloat(ix)*1.0d-03
Q=10d0
call parden(x,Q,f,i)
g1=f(0)
Q=100d0
call parden(x,Q,f,i)
g2=f(0)
Q=1000d0
call parden(x,Q,f,i)
g3=f(0)
write(*,*) x,g1,g2,g3
enddo
a=alphasPDF(QMZ)
write(*,*)
write(*,*) 'alpha_S(M_Z) = ',a
write(*,*)
enddo
end program example3
subroutine parden(x,Q,f,imem)
implicit none
character*64 name
data name/'Alekhin_100.LHpdf'/
double precision x,Q,f(-6:6)
integer imem,init,lmem
data init/0/
data lmem/-1/
save init,lmem
if (init.eq.0) then
init=1
call InitPDFsetByName(name)
endif
if (imem.ne.lmem) then
lmem=imem
call InitPDF(lmem)
endif
call evolvePDF(x,Q,f)
return
end subroutine parden
|
