/usr/share/doc/liblhapdf-dev/examples/Example1.f is in liblhapdf-dev 5.9.1-3+b1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
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 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 | program example1
implicit double precision (a-h,o-z)
character name*64
double precision f(-6:6)
character*20 lparm
logical has_photon
Dimension Z(10), XX(10)
Data (Z(I), I=1,10) /.05, .1, .2, .3, .4, .5, .6, .7, .8, .9/
Do I = 1, 10
XX(I) = Z(I) **3
EndDo
name='cteq61.LHgrid'
! name='cteq65.LHgrid'
! name='cteq66.LHgrid'
! name='MRST2006nnlo.LHgrid'
! name='MRST2001E.LHgrid'
! name='H12000ms.LHgrid'
call InitPDFsetByName(name)
QMZ=91.18d0
write(*,*)
call numberPDF(N)
print *,'There are ',N,' PDF sets'
do i=0,N
write(*,*) '---------------------------------------------'
call InitPDF(i)
write(*,*) 'PDF set ',i
call GetXmin(i,xmin)
call GetXmax(i,xmax)
call GetQ2min(i,q2min)
call GetQ2max(i,q2max)
print *,'xmin=',xmin,' xmax=',xmax,' Q2min=',q2min,' Q2max=',q2max
call GetMinMax(i,xmin,xmax,q2min,q2max)
print *,'xmin=',xmin,' xmax=',xmax,' Q2min=',q2min,' Q2max=',q2max
call setlhaparm('EXTRAPOLATE')
call getlhaparm(18,lparm)
print *,'lhaparm(18)=',lparm
write(*,*)
a=alphasPDF(QMZ)
write(*,*) 'alpha_S(M_Z) = ',a
call getLam4M(1,i,xlam4)
call getLam5M(1,i,xlam5)
print *,' lambda5: ',xlam5, ' lambda4: ',xlam4
write(*,*)
write(*,*) 'x*up'
write(*,*) ' x Q=10 GeV Q=100 GeV Q=1000 GeV'
! q2 = 10.0d0
! q = dsqrt(q2)
q = 50.0d0
print *,q
do ix=1,10
! x = (ix-0.5d0)/10.0d0
x = xx(ix)
! x = z(ix)
if(has_photon()) then
print *,"This set has a photon"
call evolvePDFphoton(x,Q,f,photon)
else
call evolvePDF(x,Q,f)
endif
g = f(0)
u = f(2)
d = f(1)
s = f(3)
c = f(4)
b = f(5)
ubar = f(-2)
dbar = f(-1)
sbar = f(-3)
cbar = f(-4)
bbar = f(-5)
write(*,'(F7.4,13(1pE10.3))') x,u,d,ubar,Dbar,s,sbar,c,cbar,b,bbar,g,photon
enddo
enddo
end program example1
|