Reading QA thresholds from external file II
[u/mrichter/AliRoot.git] / LHAPDF / lhapdf5.3.1 / evolution.f
CommitLineData
4e9e3152 1 subroutine GetOrderPDF(order)
2 implicit none
3 integer nset,order
4 real*8 Q,Q2
5 nset = 1
6 call GetOrderPDFM(nset,order)
7 return
8c
9 entry GetRenFac(Q)
10 nset = 1
11 call GetRenFacM(nset,Q)
12 return
13c
14 entry GetQ2fit(Q2)
15 nset = 1
16 call GetQ2fitM(nset,Q2)
17 return
18 end
19c
20 subroutine InitEvolve(nset)
21 implicit none
22 character*16 s2
23 real*8 x,Q,Q2,Q2fit,muR
24 integer nset,order,EvlOrd
25c
26 save EvlOrd,Q2fit,muR
27*
28 EvlOrd=-1
29 read(1,*) s2,Q2fit,muR
30 if (index(s2,'lo').eq.1) EvlOrd=0
31 if (index(s2,'nlo').eq.1) EvlOrd=1
32 if (index(s2,'nnlo').eq.1) EvlOrd=2
33 if (EvlOrd.lt.0) then
34 write(*,*) 'File description error:'
35 write(*,*) 'Unknown PDF evolution order ',s2
36 stop
37 endif
38 if (muR.ne.1.0d0) then
39 write(*,*) '***********************************************'
40 write(*,*) '* Note than the renormalization scale is *'
41 write(*,*) '* unequal to the factorization scale for this *'
42 write(*,*) '* particular PDF set. *'
43 write(*,*) '* See manual for proper use. *'
44 write(*,*) '***********************************************'
45 endif
46c print *,'calling readevolve', nset
47 call readevolve(nset)
48 return
49*
50 entry InitEvolveCode(nset)
51c print *,'calling initevolution', nset,Evlord,Q2fit
52 call initevolution(nset,EvlOrd,Q2fit)
53 return
54*
55 entry GetOrderPDFM(nset,order)
56 order=EvlOrd
57 return
58*
59 entry GetRenFacM(nset,Q)
60 Q=muR
61 return
62*
63 entry GetQ2fitM(nset,Q2)
64 Q2=Q2fit
65 return
66*
67 end