| 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 |
| 51 | c integer iset,imem |
| 52 | c 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 |
| 72 | c |
| 73 | entry GetLam5M(nset,mem,lam5) |
| 74 | lam5 = parmQCD(nset,mem,2) |
| 75 | return |
| 76 | c |
| 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 | |
git://git.uio.no / u / mrichter / AliRoot.git / blame