| 21886bb6 |
1 | #include "pdf/pilot.h" |
| 2 | c --------------------------------------------------------------------- |
| 3 | double precision function faux5MI(ifl,x,q) |
| 4 | c auxiliary function for parametrization of CTEQ5MI (J. Pumplin 9/99). |
| 5 | c --------------------------------------------------------------------- |
| 6 | #include "pdf/impdp.inc" |
| 7 | integer ifl |
| 8 | |
| 9 | parameter (nex=8, nlf=2) |
| 10 | dimension am(0:nex,0:nlf,-5:2) |
| 11 | dimension alfvec(-5:2), qmavec(-5:2) |
| 12 | dimension mexvec(-5:2), mlfvec(-5:2) |
| 13 | dimension ut1vec(-5:2), ut2vec(-5:2) |
| 14 | dimension af(0:nex) |
| 15 | |
| 16 | data mexvec( 2) / 8 / |
| 17 | data mlfvec( 2) / 2 / |
| 18 | data ut1vec( 2) / 0.5141718E+01 / |
| 19 | data ut2vec( 2) / -0.1346944E+01 / |
| 20 | data alfvec( 2) / 0.5260555E+00 / |
| 21 | data qmavec( 2) / 0.0000000E+00 / |
| 22 | data (am( 0,k, 2),k=0, 2) |
| 23 | & / 0.4289071E+01, -0.2536870E+01, -0.1259948E+01 / |
| 24 | data (am( 1,k, 2),k=0, 2) |
| 25 | & / 0.9839410E+00, 0.4168426E-01, -0.5018952E-01 / |
| 26 | data (am( 2,k, 2),k=0, 2) |
| 27 | & / -0.1651961E+02, 0.9246261E+01, 0.5996400E+01 / |
| 28 | data (am( 3,k, 2),k=0, 2) |
| 29 | & / -0.2077936E+02, 0.9786469E+01, 0.7656465E+01 / |
| 30 | data (am( 4,k, 2),k=0, 2) |
| 31 | & / 0.3054926E+02, 0.1889536E+01, 0.1380541E+01 / |
| 32 | data (am( 5,k, 2),k=0, 2) |
| 33 | & / 0.3084695E+02, -0.1212303E+02, -0.1053551E+02 / |
| 34 | data (am( 6,k, 2),k=0, 2) |
| 35 | & / -0.1426778E+02, 0.6239537E+01, 0.5254819E+01 / |
| 36 | data (am( 7,k, 2),k=0, 2) |
| 37 | & / -0.1909811E+02, 0.3695678E+01, 0.5495729E+01 / |
| 38 | data (am( 8,k, 2),k=0, 2) |
| 39 | & / 0.1889751E-01, 0.5027193E-02, 0.6624896E-03 / |
| 40 | |
| 41 | data mexvec( 1) / 8 / |
| 42 | data mlfvec( 1) / 2 / |
| 43 | data ut1vec( 1) / 0.4138426E+01 / |
| 44 | data ut2vec( 1) / -0.3221374E+01 / |
| 45 | data alfvec( 1) / 0.4960962E+00 / |
| 46 | data qmavec( 1) / 0.0000000E+00 / |
| 47 | data (am( 0,k, 1),k=0, 2) |
| 48 | & / 0.1332497E+01, -0.3703718E+00, 0.1288638E+00 / |
| 49 | data (am( 1,k, 1),k=0, 2) |
| 50 | & / 0.7544687E+00, 0.3255075E-01, -0.4706680E-01 / |
| 51 | data (am( 2,k, 1),k=0, 2) |
| 52 | & / -0.7638814E+00, 0.5008313E+00, -0.9237374E-01 / |
| 53 | data (am( 3,k, 1),k=0, 2) |
| 54 | & / -0.3689889E+00, -0.1055098E+01, -0.4645065E+00 / |
| 55 | data (am( 4,k, 1),k=0, 2) |
| 56 | & / 0.3991610E+02, 0.1979881E+01, 0.1775814E+01 / |
| 57 | data (am( 5,k, 1),k=0, 2) |
| 58 | & / 0.6201080E+01, 0.2046288E+01, 0.3804571E+00 / |
| 59 | data (am( 6,k, 1),k=0, 2) |
| 60 | & / -0.8027900E+00, -0.7011688E+00, -0.8049612E+00 / |
| 61 | data (am( 7,k, 1),k=0, 2) |
| 62 | & / -0.8631305E+01, -0.3981200E+01, 0.6970153E+00 / |
| 63 | data (am( 8,k, 1),k=0, 2) |
| 64 | & / 0.2371230E-01, 0.5372683E-02, 0.1118701E-02 / |
| 65 | |
| 66 | data mexvec( 0) / 8 / |
| 67 | data mlfvec( 0) / 2 / |
| 68 | data ut1vec( 0) / -0.1026789E+01 / |
| 69 | data ut2vec( 0) / -0.9051707E+01 / |
| 70 | data alfvec( 0) / 0.9462977E+00 / |
| 71 | data qmavec( 0) / 0.0000000E+00 / |
| 72 | data (am( 0,k, 0),k=0, 2) |
| 73 | & / 0.1191990E+03, -0.8548739E+00, -0.1963040E+01 / |
| 74 | data (am( 1,k, 0),k=0, 2) |
| 75 | & / -0.9449972E+02, 0.1074771E+01, 0.2056055E+01 / |
| 76 | data (am( 2,k, 0),k=0, 2) |
| 77 | & / 0.3701064E+01, -0.1167947E-02, 0.1933573E+00 / |
| 78 | data (am( 3,k, 0),k=0, 2) |
| 79 | & / 0.1171345E+03, -0.1064540E+01, -0.1875312E+01 / |
| 80 | data (am( 4,k, 0),k=0, 2) |
| 81 | & / -0.1014453E+03, -0.5707427E+00, 0.4511242E-01 / |
| 82 | data (am( 5,k, 0),k=0, 2) |
| 83 | & / 0.6365168E+01, 0.1275354E+01, -0.4964081E+00 / |
| 84 | data (am( 6,k, 0),k=0, 2) |
| 85 | & / -0.3370693E+01, -0.1122020E+01, 0.5947751E-01 / |
| 86 | data (am( 7,k, 0),k=0, 2) |
| 87 | & / -0.5327270E+01, -0.9293556E+00, 0.6629940E+00 / |
| 88 | data (am( 8,k, 0),k=0, 2) |
| 89 | & / 0.2437513E-01, 0.1600939E-02, 0.6855336E-03 / |
| 90 | |
| 91 | data mexvec(-1) / 8 / |
| 92 | data mlfvec(-1) / 2 / |
| 93 | data ut1vec(-1) / 0.5243571E+01 / |
| 94 | data ut2vec(-1) / -0.2870513E+01 / |
| 95 | data alfvec(-1) / 0.6701448E+00 / |
| 96 | data qmavec(-1) / 0.0000000E+00 / |
| 97 | data (am( 0,k,-1),k=0, 2) |
| 98 | & / 0.2428863E+02, 0.1907035E+01, -0.4606457E+00 / |
| 99 | data (am( 1,k,-1),k=0, 2) |
| 100 | & / 0.2006810E+01, -0.1265915E+00, 0.7153556E-02 / |
| 101 | data (am( 2,k,-1),k=0, 2) |
| 102 | & / -0.1884546E+02, -0.2339471E+01, 0.5740679E+01 / |
| 103 | data (am( 3,k,-1),k=0, 2) |
| 104 | & / -0.2527892E+02, -0.2044124E+01, 0.1280470E+02 / |
| 105 | data (am( 4,k,-1),k=0, 2) |
| 106 | & / -0.1013824E+03, -0.1594199E+01, 0.2216401E+00 / |
| 107 | data (am( 5,k,-1),k=0, 2) |
| 108 | & / 0.8070930E+02, 0.1792072E+01, -0.2164364E+02 / |
| 109 | data (am( 6,k,-1),k=0, 2) |
| 110 | & / -0.4641050E+02, 0.1977338E+00, 0.1273014E+02 / |
| 111 | data (am( 7,k,-1),k=0, 2) |
| 112 | & / -0.3910568E+02, 0.1719632E+01, 0.1086525E+02 / |
| 113 | data (am( 8,k,-1),k=0, 2) |
| 114 | & / -0.1185496E+01, -0.1905847E+00, -0.8744118E-03 / |
| 115 | |
| 116 | data mexvec(-2) / 7 / |
| 117 | data mlfvec(-2) / 2 / |
| 118 | data ut1vec(-2) / 0.4782210E+01 / |
| 119 | data ut2vec(-2) / -0.1976856E+02 / |
| 120 | data alfvec(-2) / 0.7558374E+00 / |
| 121 | data qmavec(-2) / 0.0000000E+00 / |
| 122 | data (am( 0,k,-2),k=0, 2) |
| 123 | & / -0.6216935E+00, 0.2369963E+00, -0.7909949E-02 / |
| 124 | data (am( 1,k,-2),k=0, 2) |
| 125 | & / 0.1245440E+01, -0.1031510E+00, 0.4916523E-02 / |
| 126 | data (am( 2,k,-2),k=0, 2) |
| 127 | & / -0.7060824E+01, -0.3875283E-01, 0.1784981E+00 / |
| 128 | data (am( 3,k,-2),k=0, 2) |
| 129 | & / -0.7430595E+01, 0.1964572E+00, -0.1284999E+00 / |
| 130 | data (am( 4,k,-2),k=0, 2) |
| 131 | & / -0.6897810E+01, 0.2620543E+01, 0.8012553E-02 / |
| 132 | data (am( 5,k,-2),k=0, 2) |
| 133 | & / 0.1507713E+02, 0.2340307E-01, 0.2482535E+01 / |
| 134 | data (am( 6,k,-2),k=0, 2) |
| 135 | & / -0.1815341E+01, -0.1538698E+01, -0.2014208E+01 / |
| 136 | data (am( 7,k,-2),k=0, 2) |
| 137 | & / -0.2571932E+02, 0.2903941E+00, -0.2848206E+01 / |
| 138 | |
| 139 | data mexvec(-3) / 7 / |
| 140 | data mlfvec(-3) / 2 / |
| 141 | data ut1vec(-3) / 0.4518239E+01 / |
| 142 | data ut2vec(-3) / -0.2690590E+01 / |
| 143 | data alfvec(-3) / 0.6124079E+00 / |
| 144 | data qmavec(-3) / 0.0000000E+00 / |
| 145 | data (am( 0,k,-3),k=0, 2) |
| 146 | & / -0.2734458E+01, -0.7245673E+00, -0.6351374E+00 / |
| 147 | data (am( 1,k,-3),k=0, 2) |
| 148 | & / 0.2927174E+01, 0.4822709E+00, -0.1088787E-01 / |
| 149 | data (am( 2,k,-3),k=0, 2) |
| 150 | & / -0.1771017E+02, -0.1416635E+01, 0.8467622E+01 / |
| 151 | data (am( 3,k,-3),k=0, 2) |
| 152 | & / -0.4972782E+02, -0.3348547E+01, 0.1767061E+02 / |
| 153 | data (am( 4,k,-3),k=0, 2) |
| 154 | & / -0.7102770E+01, -0.3205337E+01, 0.4101704E+00 / |
| 155 | data (am( 5,k,-3),k=0, 2) |
| 156 | & / 0.7169698E+02, -0.2205985E+01, -0.2463931E+02 / |
| 157 | data (am( 6,k,-3),k=0, 2) |
| 158 | & / -0.4090347E+02, 0.2103486E+01, 0.1416507E+02 / |
| 159 | data (am( 7,k,-3),k=0, 2) |
| 160 | & / -0.2952639E+02, 0.5376136E+01, 0.7825585E+01 / |
| 161 | |
| 162 | data mexvec(-4) / 7 / |
| 163 | data mlfvec(-4) / 2 / |
| 164 | data ut1vec(-4) / 0.2783230E+01 / |
| 165 | data ut2vec(-4) / -0.1746328E+01 / |
| 166 | data alfvec(-4) / 0.1115653E+01 / |
| 167 | data qmavec(-4) / 0.1300000E+01 / |
| 168 | data (am( 0,k,-4),k=0, 2) |
| 169 | & / -0.1743872E+01, -0.1128921E+01, -0.2841969E+00 / |
| 170 | data (am( 1,k,-4),k=0, 2) |
| 171 | & / 0.3345755E+01, 0.3187765E+00, 0.1378124E+00 / |
| 172 | data (am( 2,k,-4),k=0, 2) |
| 173 | & / -0.2037615E+02, 0.4121687E+01, 0.2236520E+00 / |
| 174 | data (am( 3,k,-4),k=0, 2) |
| 175 | & / -0.4703104E+02, 0.5353087E+01, -0.1455347E+01 / |
| 176 | data (am( 4,k,-4),k=0, 2) |
| 177 | & / -0.1060230E+02, -0.1551122E+01, -0.1078863E+01 / |
| 178 | data (am( 5,k,-4),k=0, 2) |
| 179 | & / 0.5088892E+02, -0.8197304E+01, 0.8083451E+01 / |
| 180 | data (am( 6,k,-4),k=0, 2) |
| 181 | & / -0.2819070E+02, 0.4554086E+01, -0.5890995E+01 / |
| 182 | data (am( 7,k,-4),k=0, 2) |
| 183 | & / -0.1098238E+02, 0.2590096E+01, -0.8062879E+01 / |
| 184 | |
| 185 | data mexvec(-5) / 6 / |
| 186 | data mlfvec(-5) / 2 / |
| 187 | data ut1vec(-5) / 0.1619654E+02 / |
| 188 | data ut2vec(-5) / -0.3367346E+01 / |
| 189 | data alfvec(-5) / 0.5109891E-02 / |
| 190 | data qmavec(-5) / 0.4500000E+01 / |
| 191 | data (am( 0,k,-5),k=0, 2) |
| 192 | & / -0.6800138E+01, 0.2493627E+01, -0.1075724E+01 / |
| 193 | data (am( 1,k,-5),k=0, 2) |
| 194 | & / 0.3036555E+01, 0.3324733E+00, 0.2008298E+00 / |
| 195 | data (am( 2,k,-5),k=0, 2) |
| 196 | & / -0.5203879E+01, -0.8493476E+01, -0.4523208E+01 / |
| 197 | data (am( 3,k,-5),k=0, 2) |
| 198 | & / -0.1524239E+01, -0.3411912E+01, -0.1771867E+02 / |
| 199 | data (am( 4,k,-5),k=0, 2) |
| 200 | & / -0.1099444E+02, 0.1320930E+01, -0.2353831E+01 / |
| 201 | data (am( 5,k,-5),k=0, 2) |
| 202 | & / 0.1699299E+02, -0.3565802E+02, 0.3566872E+02 / |
| 203 | data (am( 6,k,-5),k=0, 2) |
| 204 | & / -0.1465793E+02, 0.2703365E+02, -0.2176372E+02 / |
| 205 | |
| 206 | if(q .le. qmavec(ifl)) then |
| 207 | faux5MI = 0.d0 |
| 208 | return |
| 209 | endif |
| 210 | |
| 211 | if(x .ge. 1.d0) then |
| 212 | faux5MI = 0.d0 |
| 213 | return |
| 214 | endif |
| 215 | |
| 216 | tmp = log(q/alfvec(ifl)) |
| 217 | if(tmp .le. 0.d0) then |
| 218 | faux5MI = 0.d0 |
| 219 | return |
| 220 | endif |
| 221 | |
| 222 | sb = log(tmp) |
| 223 | sb1 = sb - 1.2d0 |
| 224 | sb2 = sb1*sb1 |
| 225 | |
| 226 | do i = 0, nex |
| 227 | af(i) = 0.d0 |
| 228 | sbx = 1.d0 |
| 229 | do k = 0, mlfvec(ifl) |
| 230 | af(i) = af(i) + sbx*am(i,k,ifl) |
| 231 | sbx = sb1*sbx |
| 232 | enddo |
| 233 | enddo |
| 234 | |
| 235 | y = -log(x) |
| 236 | u = log(x/0.00001d0) |
| 237 | |
| 238 | part1 = af(1)*y**(1.d0+0.01d0*af(4))*(1.d0+ af(8)*u) |
| 239 | part2 = af(0)*(1.d0 - x) + af(3)*x |
| 240 | part3 = x*(1.d0-x)*(af(5)+af(6)*(1.d0-x)+af(7)*x*(1.d0-x)) |
| 241 | part4 = ut1vec(ifl)*log(1.d0-x) + |
| 242 | & AF(2)*log(1.d0+exp(ut2vec(ifl))-x) |
| 243 | |
| 244 | faux5MI = exp(log(x) + part1 + part2 + part3 + part4) |
| 245 | |
| 246 | c include threshold factor... |
| 247 | faux5MI = faux5MI * (1.d0 - qmavec(ifl)/q) |
| 248 | |
| 249 | return |
| 250 | end |
git://git.uio.no / u / mrichter / AliRoot.git / blame