Reading QA thresholds from external file II
[u/mrichter/AliRoot.git] / LHAPDF / lhapdf5.3.1 / QCDparams.f
CommitLineData
4e9e3152 1 subroutine GetLam4(mem,lam4)
2 implicit none
3 integer mem,nset
4 real*8 lam4,lam5
5 nset = 1
6 call GetLam4M(nset,mem,lam4)
7 return
8
9 entry GetLam5(mem,lam5)
10 nset = 1
11 call GetLam5M(nset,mem,lam5)
12 return
13
14 end
15
16 subroutine GetXmin(mem,xmin)
17 implicit none
18 integer mem,nset
19 real*8 xmin,xmax,q2min,q2max
20 nset = 1
21 call GetXminM(nset,mem,xmin)
22 return
23
24 entry GetXmax(mem,xmax)
25 nset = 1
26 call GetXmaxM(nset,mem,xmax)
27 return
28
29 entry GetQ2min(mem,q2min)
30 nset = 1
31 call GetQ2minM(nset,mem,q2min)
32 return
33
34 entry GetQ2max(mem,q2max)
35 nset = 1
36 call GetQ2maxM(nset,mem,q2max)
37 return
38
39 entry GetMinMax(mem,xmin,xmax,q2min,q2max)
40 nset = 1
41 call GetMinMaxM(nset,mem,xmin,xmax,q2min,q2max)
42 return
43
44 end
45
46 subroutine initQCDparams(nset)
47 implicit real*8(a-h,o-z)
48 include 'parmsetup.inc'
49 real*8 parmQCD(nmxset,0:noemax,2),lam4,lam5
50 integer nset
51c integer iset,imem
52c common/SET/iset,imem
53 save
54 read(1,*)nmem,nrep
55 if(nrep.eq.0) then
56 do i=0,nmem
57 read(1,*) (parmQCD(nset,i,j),j=1,2)
58 enddo
59 else
60 read(1,*) (parmQCD(nset,0,j),j=1,2)
61 do i=1,nmem
62 do j=1,2
63 parmQCD(nset,i,j)=parmQCD(nset,0,j)
64 enddo
65 enddo
66 endif
67 return
68
69 entry GetLam4M(nset,mem,lam4)
70 lam4 = parmQCD(nset,mem,1)
71 return
72c
73 entry GetLam5M(nset,mem,lam5)
74 lam5 = parmQCD(nset,mem,2)
75 return
76c
77 end
78
79 subroutine initMinMax(nset)
80 implicit real*8(a-h,o-z)
81 include 'parmsetup.inc'
82 real*8 parmXmin(nmxset,0:noemax),xmin
83 real*8 parmXmax(nmxset,0:noemax),xmax
84 real*8 parmQ2min(nmxset,0:noemax),q2min
85 real*8 parmQ2max(nmxset,0:noemax),q2max
86 integer nset
87 save
88 read(1,*)nmem,nrep
89 if(nrep.eq.0) then
90 do i=0,nmem
91 read(1,*) parmXmin(nset,i),
92 + parmXmax(nset,i),
93 + parmQ2min(nset,i),
94 + parmQ2max(nset,i)
95 enddo
96 else
97 read(1,*) parmXmin(nset,0),
98 + parmXmax(nset,0),
99 + parmQ2min(nset,0),
100 + parmQ2max(nset,0)
101 do i=1,nmem
102 parmXmin(nset,i) = parmXmin(nset,0)
103 parmXmax(nset,i) = parmXmax(nset,0)
104 parmQ2min(nset,i) = parmQ2min(nset,0)
105 parmQ2max(nset,i) = parmQ2max(nset,0)
106 enddo
107 endif
108 return
109
110 entry GetXminM(nset,mem,xmin)
111 xmin = parmXmin(nset,mem)
112 return
113
114 entry GetXmaxM(nset,mem,xmax)
115 xmax = parmXmax(nset,mem)
116 return
117
118 entry GetQ2minM(nset,mem,q2min)
119 q2min = parmQ2min(nset,mem)
120 return
121
122 entry GetQ2maxM(nset,mem,q2max)
123 q2max = parmQ2max(nset,mem)
124 return
125
126 entry GetMinMaxM(nset,mem,xmin,xmax,q2min,q2max)
127 xmin = parmXmin(nset,mem)
128 xmax = parmXmax(nset,mem)
129 q2min = parmQ2min(nset,mem)
130 q2max = parmQ2max(nset,mem)
131 return
132
133 end
134
135
136