]>
Commit | Line | Data |
---|---|---|
3820ca8e | 1 | |
2 | CDECK ID>, HWSFBR. | |
3 | ||
4 | *CMZ :- -15/07/92 14.08.45 by Mike Seymour | |
5 | ||
6 | *-- Author : Bryan Webber | |
7 | ||
8 | C----------------------------------------------------------------------- | |
9 | ||
10 | SUBROUTINE HWSFBR(X,QQ,FORCED,ID,IW,ID1,ID2,IW1,IW2,Z) | |
11 | ||
12 | C----------------------------------------------------------------------- | |
13 | ||
14 | C FINDS BRANCHING (ID1->ID+ID2) AND Z=X/X1 IN BACKWARD | |
15 | ||
16 | C EVOLUTION AT ENERGY FRACTION X AND SCALE QQ | |
17 | ||
18 | C | |
19 | ||
20 | C FORCED=.TRUE. FORCES SPLITTING OF NON-VALENCE PARTON | |
21 | ||
22 | C | |
23 | ||
24 | C IW,IW1,IW2 ARE COLOUR CONNECTION WORDS | |
25 | ||
26 | C | |
27 | ||
28 | C ID1.LT.0 ON RETURN MEANS NO PHASE SPACE | |
29 | ||
30 | C ID1.EQ.0 ON RETURN FLAGS REJECTED BRANCHINGS | |
31 | ||
32 | C----------------------------------------------------------------------- | |
33 | ||
34 | INCLUDE 'HERWIG61.INC' | |
35 | ||
36 | DOUBLE PRECISION HWBVMC,HWR,HWUALF,HWUAEM,QP,X,QQ,Z,WQG,WQV, | |
37 | ||
38 | & WQP,XQV,ZMIN,ZMAX,YMIN,YMAX,DELY,YY,PSUM,EZ,WQN,WR,ZR,WZ,ZZ,AZ, | |
39 | ||
40 | & PVAL,EY,DIST(13),PROB(13,100),PPHO | |
41 | ||
42 | INTEGER ID,IW,ID1,ID2,IW1,IW2,NZ,IDHAD,IP,IZ | |
43 | ||
44 | LOGICAL HWRLOG,HWSVAL,FORCED,NONF,NONV,PHOTPR | |
45 | ||
46 | EXTERNAL HWBVMC,HWR,HWUALF,HWUAEM,HWRLOG,HWSVAL | |
47 | ||
48 | ID1=-1 | |
49 | ||
50 | QP=HWBVMC(ID) | |
51 | ||
52 | WQG=1.-QG/QQ | |
53 | ||
54 | WQV=1.-QV/QQ | |
55 | ||
56 | WQP=1.-QP/QQ | |
57 | ||
58 | XQV=X/WQV | |
59 | ||
60 | NONV=.NOT.HWSVAL(ID) | |
61 | ||
62 | NONF=.NOT.FORCED | |
63 | ||
64 | 5 IF (ID.EQ.13) THEN | |
65 | ||
66 | ZMIN=X | |
67 | ||
68 | IF (NONF) THEN | |
69 | ||
70 | ZMAX=WQG | |
71 | ||
72 | ELSE | |
73 | ||
74 | ZMAX=WQV | |
75 | ||
76 | ENDIF | |
77 | ||
78 | ELSE | |
79 | ||
80 | IF (NONV) THEN | |
81 | ||
82 | ZMIN=XQV | |
83 | ||
84 | IF (NONF) THEN | |
85 | ||
86 | ZMAX=WQG | |
87 | ||
88 | ELSE | |
89 | ||
90 | ZMAX=WQP | |
91 | ||
92 | ENDIF | |
93 | ||
94 | ELSE | |
95 | ||
96 | ZMIN=X | |
97 | ||
98 | ZMAX=MAX(WQG,WQP) | |
99 | ||
100 | ENDIF | |
101 | ||
102 | ENDIF | |
103 | ||
104 | IF (ZMIN.GE.ZMAX) RETURN | |
105 | ||
106 | ID1=0 | |
107 | ||
108 | C---INTERPOLATION VARIABLE IS Y=LN(Z/(1-Z)) | |
109 | ||
110 | YMIN=LOG(ZMIN/(1.-ZMIN)) | |
111 | ||
112 | YMAX=LOG(ZMAX/(1.-ZMAX)) | |
113 | ||
114 | DELY=YMAX-YMIN | |
115 | ||
116 | NZ=MIN(INT(ZBINM*DELY)+1,NZBIN) | |
117 | ||
118 | DELY=(YMAX-YMIN)/FLOAT(NZ) | |
119 | ||
120 | YY=YMIN+0.5*DELY | |
121 | ||
122 | PSUM=0. | |
123 | ||
124 | IDHAD=IDHW(INHAD) | |
125 | ||
126 | C---SET UP TABLES FOR CHOOSING BRANCHING | |
127 | ||
128 | DO 40 IZ=1,NZ | |
129 | ||
130 | EZ=EXP(YY) | |
131 | ||
132 | WR=1.+EZ | |
133 | ||
134 | ZR=WR/EZ | |
135 | ||
136 | WZ=1./WR | |
137 | ||
138 | ZZ=WZ*EZ | |
139 | ||
140 | AZ=WZ*ZZ*HWUALF(5-2*SUDORD,MAX(WZ*QQ,QG)) | |
141 | ||
142 | CALL HWSFUN(X*ZR,QQ,IDHAD,NSTRU,DIST,JNHAD) | |
143 | ||
144 | IF (ID.NE.13) THEN | |
145 | ||
146 | C---SPLITTING INTO QUARK | |
147 | ||
148 | DO 10 IP=1,ID-1 | |
149 | ||
150 | 10 PROB(IP,IZ)=PSUM | |
151 | ||
152 | IF (NONF) PSUM=PSUM+DIST(ID)*AZ*CFFAC*(1.+ZZ*ZZ)*WR | |
153 | ||
154 | DO 20 IP=ID,12 | |
155 | ||
156 | 20 PROB(IP,IZ)=PSUM | |
157 | ||
158 | PSUM=PSUM+DIST(13)*AZ*0.5*(ZZ*ZZ+WZ*WZ) | |
159 | ||
160 | PROB(13,IZ)=PSUM | |
161 | ||
162 | ELSE | |
163 | ||
164 | C---SPLITTING INTO GLUON | |
165 | ||
166 | DO 30 IP=1,12 | |
167 | ||
168 | PSUM=PSUM+DIST(IP)*AZ*CFFAC*(1.+WZ*WZ)*ZR | |
169 | ||
170 | 30 PROB(IP,IZ)=PSUM | |
171 | ||
172 | IF (NONF) PSUM=PSUM+DIST(13)*AZ*2.*CAFAC*(WZ*ZR+ZZ*WR+WZ*ZZ) | |
173 | ||
174 | PROB(13,IZ)=PSUM | |
175 | ||
176 | ENDIF | |
177 | ||
178 | 40 YY=YY+DELY | |
179 | ||
180 | 50 PHOTPR=IDHAD.EQ.59.AND.ID.NE.13 | |
181 | ||
182 | IF (PHOTPR) THEN | |
183 | ||
184 | C---ALLOW ANOMALOUS PHOTON SPLITTING | |
185 | ||
186 | PPHO=HWUAEM(-QQ*QQ)*CAFAC*(ZMIN**2+(1.-ZMIN)**2) | |
187 | ||
188 | & *ICHRG(ID)**2/(18.*PIFAC) | |
189 | ||
190 | IF (PPHO.GT.(PPHO+PSUM*DELY)*HWR()) THEN | |
191 | ||
192 | C---ANOMALOUS PHOTON SPLITTING OCCURRED | |
193 | ||
194 | ID1=59 | |
195 | ||
196 | RETURN | |
197 | ||
198 | ENDIF | |
199 | ||
200 | ENDIF | |
201 | ||
202 | IF (PSUM.LE.ZERO) RETURN | |
203 | ||
204 | C---CHOOSE Z | |
205 | ||
206 | PVAL=PSUM*HWR() | |
207 | ||
208 | DO 60 IZ=1,NZ | |
209 | ||
210 | IF (PROB(13,IZ).GT.PVAL) GOTO 70 | |
211 | ||
212 | 60 CONTINUE | |
213 | ||
214 | IZ=NZ | |
215 | ||
216 | 70 EY=EXP(YMIN+DELY*(FLOAT(IZ)-HWR())) | |
217 | ||
218 | ZZ=EY/(1.+EY) | |
219 | ||
220 | C---CHOOSE BRANCHING | |
221 | ||
222 | DO 80 IP=1,13 | |
223 | ||
224 | IF (PROB(IP,IZ).GT.PVAL) GOTO 90 | |
225 | ||
226 | 80 CONTINUE | |
227 | ||
228 | IP=13 | |
229 | ||
230 | C---CHECK THAT Z IS INSIDE PHASE SPACE (RETURN IF NOT) | |
231 | ||
232 | 90 CONTINUE | |
233 | ||
234 | IF (ID.NE.13) THEN | |
235 | ||
236 | IF (IP.EQ.ID) THEN | |
237 | ||
238 | IF ((NONV.AND.ZZ*WQP.LT.XQV).OR.ZZ.GT.WQG) THEN | |
239 | ||
240 | IF (PHOTPR) GOTO 50 | |
241 | ||
242 | RETURN | |
243 | ||
244 | ENDIF | |
245 | ||
246 | ELSE | |
247 | ||
248 | IF (ZZ.LT.XQV.OR.ZZ.GT.WQP) THEN | |
249 | ||
250 | IF (PHOTPR) GOTO 50 | |
251 | ||
252 | RETURN | |
253 | ||
254 | ENDIF | |
255 | ||
256 | ENDIF | |
257 | ||
258 | ELSE | |
259 | ||
260 | IF (IP.EQ.ID) THEN | |
261 | ||
262 | IF (ZZ.LT.XQV.OR.ZZ.GT.WQG) RETURN | |
263 | ||
264 | ELSEIF (.NOT.HWSVAL(IP)) THEN | |
265 | ||
266 | WQN=1.-HWBVMC(IP)/QQ | |
267 | ||
268 | IF (ZZ*WQN.LT.XQV.OR.ZZ.GT.WQN) RETURN | |
269 | ||
270 | ENDIF | |
271 | ||
272 | ENDIF | |
273 | ||
274 | C---EVERYTHING OK: LABEL NEW BRANCHES | |
275 | ||
276 | Z=ZZ | |
277 | ||
278 | ID1=IP | |
279 | ||
280 | IW1=IW*2 | |
281 | ||
282 | IW2=IW1+1 | |
283 | ||
284 | IF (ID.LE.6) THEN | |
285 | ||
286 | IF (ID1.EQ.13) THEN | |
287 | ||
288 | ID2=ID+6 | |
289 | ||
290 | ELSE | |
291 | ||
292 | ID2=13 | |
293 | ||
294 | IW2=IW1 | |
295 | ||
296 | ENDIF | |
297 | ||
298 | ELSE IF (ID.NE.13) THEN | |
299 | ||
300 | IF (ID1.EQ.13) THEN | |
301 | ||
302 | ID2=ID-6 | |
303 | ||
304 | IW2=IW1 | |
305 | ||
306 | ELSE | |
307 | ||
308 | ID2=13 | |
309 | ||
310 | ENDIF | |
311 | ||
312 | ELSE | |
313 | ||
314 | ID2=ID1 | |
315 | ||
316 | IF (ID1.EQ.13) THEN | |
317 | ||
318 | IF (HWRLOG(HALF)) IW2=IW1 | |
319 | ||
320 | ELSE IF (ID1.GT.6) THEN | |
321 | ||
322 | IW2=IW1 | |
323 | ||
324 | END IF | |
325 | ||
326 | END IF | |
327 | ||
328 | IF (IW2.EQ.IW1) IW1=IW1+1 | |
329 | ||
330 | 999 END | |
331 | ||
332 | CDECK ID>, HWSFUN. | |
333 | ||
334 | *CMZ :- -02/05/91 11.30.51 by Federico Carminati | |
335 | ||
336 | *-- Author : Miscellaneous, combined by Bryan Webber | |
337 | ||
338 | C----------------------------------------------------------------------- | |
339 | ||
340 | SUBROUTINE HWSFUN(X,SCALE,IDHAD,NSET,DIST,IBEAM) | |
341 | ||
342 | C----------------------------------------------------------------------- | |
343 | ||
344 | C NUCLEON AND PION STRUCTURE FUNCTIONS DIST=X*QRK(X,Q=SCALE) | |
345 | ||
346 | C | |
347 | ||
348 | C IDHAD = TYPE OF HADRON: | |
349 | ||
350 | C 73=P 91=PBAR 75=N 93=NBAR 38=PI+ 30=PI- 59=PHOTON | |
351 | ||
352 | C | |
353 | ||
354 | C NEW SPECIAL CODES: | |
355 | ||
356 | C 71=`REMNANT PHOTON' 72=`REMNANT NUCLEON' | |
357 | ||
358 | C | |
359 | ||
360 | C NSET = STRUCTURE FUNCTION SET | |
361 | ||
362 | C = 1,2 FOR DUKE+OWENS SETS 1,2 (SOFT/HARD GLUE) | |
363 | ||
364 | C = 3,4 FOR EICHTEN ET AL SETS 1,2 (NUCLEON ONLY) | |
365 | ||
366 | C = 5 FOR OWENS SET 1.1 (PREPRINT FSU-HEP-910606) | |
367 | ||
368 | C | |
369 | ||
370 | C FOR PHOTON DREES+GRASSIE IS USED | |
371 | ||
372 | C | |
373 | ||
374 | C N.B. IF IBEAM.GT.0.AND.MODPDF(IBEAM).GE.0 THEN NSET IS | |
375 | ||
376 | C IGNORED AND CERN PDFLIB WITH AUTHOR GROUP=AUTPDF(IBEAM) AND | |
377 | ||
378 | C SET=MODPDF(IBEAM) IS USED. FOR COMPATABILITY WITH VERSIONS 3 | |
379 | ||
380 | C AND EARLIER, AUTPDF SHOULD BE SET TO 'MODE' | |
381 | ||
382 | C NOTE THAT NO CONSISTENCY CHECK IS MADE, FOR EXAMPLE THAT THE | |
383 | ||
384 | C REQUESTED SET FOR A PHOTON IS ACTUALLY A PHOTON SET | |
385 | ||
386 | C | |
387 | ||
388 | C IF (ISPAC.GT.0) SCALE IS REPLACED BY MAX(SCALE,QSPAC) | |
389 | ||
390 | C | |
391 | ||
392 | C FOR PHOTON, IF (PHOMAS.GT.0) THEN QUARK DISTRIBUTIONS ARE | |
393 | ||
394 | C SUPPRESSED BY LOG((Q**2+PHOMAS**2)/(P**2+PHOMAS**2)) | |
395 | ||
396 | C L = -------------------------------------- , | |
397 | ||
398 | C LOG((Q**2+PHOMAS**2)/( PHOMAS**2)) | |
399 | ||
400 | C WHILE GLUON DISTRIBUTIONS ARE SUPPRESSED BY L**2, | |
401 | ||
402 | C WHERE Q=SCALE AND P=VIRTUALITY OF THE PHOTON | |
403 | ||
404 | C | |
405 | ||
406 | C DUKE+OWENS = D.W.DUKE AND J.F.OWENS, PHYS. REV. D30 (1984) 49 (P/N) | |
407 | ||
408 | C + J.F.OWENS, PHYS. REV. D30 (1984) 943 (PI+/-) | |
409 | ||
410 | C WITH EXTRA SIGNIFICANT FIGURES VIA ED BERGER | |
411 | ||
412 | C WARNING....MOMENTUM SUM RULE BADLY VIOLATED ABOVE 1 TEV | |
413 | ||
414 | C DUKE+OWENS SETS 1,2 OBSOLETE. SET 1 UPDATED TO OWENS 1.1 (1991) | |
415 | ||
416 | C PION NOT RELIABLE ABOVE SCALE = 50 GEV | |
417 | ||
418 | C | |
419 | ||
420 | C EICHTEN ET AL = E.EICHTEN,I.HINCHLIFFE,K.LANE AND C.QUIGG, | |
421 | ||
422 | C REV. MOD. PHYS. 56 (1984) 579 | |
423 | ||
424 | C REVISED AS IN REV. MOD. PHYS. 58 (1986) 1065 | |
425 | ||
426 | C RELIABLE RANGE : SQRT(5)GEV < SCALE < 10TEV, 1E-4 < X < 1 | |
427 | ||
428 | C | |
429 | ||
430 | C DREES+GRASSIE = M.DREES & K.GRASSIE, ZEIT. PHYS. C28 (1985) 451 | |
431 | ||
432 | C MODIFIED IN M.DREES & C.S.KIM, DESY 91-039 | |
433 | ||
434 | C AND C.S.KIM, DTP/91/16 FOR HEAVY QUARKS | |
435 | ||
436 | C | |
437 | ||
438 | C FOR CERN PDFLIB DETAILS SEE PDFLIB DOC Q ON CERNVM OR | |
439 | ||
440 | C CERN_ROOT:[DOC]PDFLIB.TXT ON VXCERN | |
441 | ||
442 | C----------------------------------------------------------------------- | |
443 | ||
444 | INCLUDE 'HERWIG61.INC' | |
445 | ||
446 | DOUBLE PRECISION HWSGAM,X,SCALE,XOLD,QOLD,XMWN,QSCA,SS,SMIN,S,T, | |
447 | ||
448 | & TMIN,TMAX,VX,AA,VT,WT,UPV,DNV,SEA,STR,CHM,BTM,TOP,GLU,WX,XQSUM, | |
449 | ||
450 | & DMIN,TPMIN,TPMAX,DIST(13),G(2),Q0(5),QL(5),F(5),A(6,5), | |
451 | ||
452 | & B(3,6,5,4),XQ(6),TX(6),TT(6),TB(6),NEHLQ(8,2),CEHLQ(6,6,2,8,2), | |
453 | ||
454 | & BB(4,6,5),VAL(20),USEA,DSEA,TOTAL,SCALEF,FAC,TBMIN(2),TTMIN(2) | |
455 | ||
456 | REAL HWSDGG,HWSDGQ,XSP,Q2,P2,W2,EMB2,EMC2,ALAM2,XPGA(-6:6),F2GM, | |
457 | ||
458 | & XPVMD,XPANL,XPANH,XPBEH,XPDIR | |
459 | ||
460 | COMMON/SASCOM/XPVMD(-6:6),XPANL(-6:6),XPANH(-6:6),XPBEH(-6:6), | |
461 | ||
462 | & XPDIR(-6:6) | |
463 | ||
464 | LOGICAL PDFWRX(2,2),PDFWRQ(2,2) | |
465 | ||
466 | DOUBLE PRECISION PDFXMN,PDFXMX,PDFQMN,PDFQMX | |
467 | ||
468 | COMMON /W50513/PDFXMN,PDFXMX,PDFQMN,PDFQMX | |
469 | ||
470 | INTEGER IDHAD,NSET,IBEAM,IOLD,NOLD,IP,I,J,K,NX,IT,IX,IFL,NFL, | |
471 | ||
472 | & MPDF,IHAD,ISET,IOP1,IOP2,IP2 | |
473 | ||
474 | CHARACTER*20 PARM(20) | |
475 | ||
476 | EXTERNAL HWSGAM,HWSDGG,HWSDGQ | |
477 | ||
478 | SAVE QOLD,IOLD,NOLD,XOLD,SS,S,T,TMIN,TMAX,G,A,TX,TT,TB,IP,NX | |
479 | ||
480 | DATA PDFWRX,PDFWRQ/8*.TRUE./ | |
481 | ||
482 | DATA (((B(I,J,K,1),I=1,3),J=1,6),K=1,5)/ | |
483 | ||
484 | &3.,0.,0.,.419,.004383,-.007412, | |
485 | ||
486 | &3.46,.72432,-.065998,4.4,-4.8644,1.3274, | |
487 | ||
488 | &6*0.,1., | |
489 | ||
490 | &0.,0.,.763,-.23696,.025836,4.,.62664,-.019163, | |
491 | ||
492 | &0.,-.42068,.032809,6*0.,1.265,-1.1323,.29268, | |
493 | ||
494 | &0.,-.37162,-.028977,8.05,1.5877,-.15291, | |
495 | ||
496 | &0.,6.3059,-.27342,0.,-10.543,-3.1674, | |
497 | ||
498 | &0.,14.698,9.798,0.,.13479,-.074693, | |
499 | ||
500 | &-.0355,-.22237,-.057685,6.3494,3.2649,-.90945, | |
501 | ||
502 | &0.,-3.0331,1.5042,0.,17.431,-11.255, | |
503 | ||
504 | &0.,-17.861,15.571,1.564,-1.7112,.63751, | |
505 | ||
506 | &0.,-.94892,.32505,6.,1.4345,-1.0485, | |
507 | ||
508 | &9.,-7.1858,.25494,0.,-16.457,10.947, | |
509 | ||
510 | &0.,15.261,-10.085/ | |
511 | ||
512 | DATA (((B(I,J,K,2),I=1,3),J=1,6),K=1,5)/ | |
513 | ||
514 | &3.,0.,0.,.3743,.013946,-.00031695, | |
515 | ||
516 | &3.329,.75343,-.076125,6.032,-6.2153,1.5561, | |
517 | ||
518 | &6*0.,1.,0., | |
519 | ||
520 | &0.,.7608,-.2317,.023232,3.83,.62746,-.019155, | |
521 | ||
522 | &0.,-.41843,.035972,6*0.,1.6714,-1.9168,.58175, | |
523 | ||
524 | &0.,-.27307,-.16392,9.145,.53045,-.76271, | |
525 | ||
526 | &0.,15.665,-2.8341,0.,-100.63,44.658, | |
527 | ||
528 | &0.,223.24,-116.76,0.,.067368,-.030574, | |
529 | ||
530 | &-.11989,-.23293,-.023273,3.5087,3.6554,-.45313, | |
531 | ||
532 | &0.,-.47369,.35793,0.,9.5041,-5.4303, | |
533 | ||
534 | &0.,-16.563,15.524,.8789,-.97093,.43388, | |
535 | ||
536 | &0.,-1.1612,.4759,4.,1.2271,-.25369, | |
537 | ||
538 | &9.,-5.6354,-.81747,0.,-7.5438,5.5034, | |
539 | ||
540 | &0.,-.59649,.12611/ | |
541 | ||
542 | DATA (((B(I,J,K,3),I=1,3),J=1,6),K=1,5)/ | |
543 | ||
544 | &1.,0.,0.,0.4,-0.06212,-0.007109,0.7,0.6478,0.01335,27*0., | |
545 | ||
546 | &0.9,-0.2428,0.1386,0.,-0.2120,0.003671,5.0,0.8673,0.04747, | |
547 | ||
548 | &0.,1.266,-2.215,0.,2.382,0.3482,3*0., | |
549 | ||
550 | &0.,0.07928,-0.06134,-0.02212,-0.3785,-0.1088,2.894,9.433, | |
551 | ||
552 | &-10.852,0.,5.248,-7.187,0.,8.388,-11.61,3*0., | |
553 | ||
554 | &0.888,-1.802,1.812,0.,-1.576,1.20,3.11,-0.1317,0.5068, | |
555 | ||
556 | &6.0,2.801,-12.16,0.,-17.28,20.49,3*0./ | |
557 | ||
558 | DATA (((B(I,J,K,4),I=1,3),J=1,6),K=1,5)/ | |
559 | ||
560 | &1.,0.,0.,0.4,-0.05909,-0.006524,0.628,0.6436,0.01451,27*0., | |
561 | ||
562 | &0.90,-0.1417,-0.1740,0.,-0.1697,-0.09623,5.0,-2.474,1.575, | |
563 | ||
564 | &0.,-2.534,1.378,0.,0.5621,-0.2701,3*0., | |
565 | ||
566 | &0.,0.06229,-0.04099,-0.0882,-0.2892,-0.1082,1.924,0.2424, | |
567 | ||
568 | &2.036,0.,-4.463,5.209,0.,-0.8367,-0.04840,3*0., | |
569 | ||
570 | &0.794,-0.9144,0.5966,0.,-1.237,0.6582,2.89,0.5966,-0.2550, | |
571 | ||
572 | &6.0,-3.671,-2.304,0.,-8.191,7.758,3*0./ | |
573 | ||
574 | C---COEFFTS FOR NEW OWENS 1.1 SET | |
575 | ||
576 | DATA BB/3.,3*0.,.665,-.1097,-.002442,0., | |
577 | ||
578 | &3.614,.8395,-.02186,0.,.8673,-1.6637,.342,0., | |
579 | ||
580 | &0.,1.1049,-.2369,5*0.,1.,3*0., | |
581 | ||
582 | &.8388,-.2092,.02657,0.,4.667,.7951,.1081,0., | |
583 | ||
584 | &0.,-1.0232,.05799,0.,0.,.8616,.153,5*0., | |
585 | ||
586 | &.909,-.4023,.006305,0., | |
587 | ||
588 | &0.,-.3823,.02766,0.,7.278,-.7904,.8108,0., | |
589 | ||
590 | &0.,-1.6629,.5719,0.,0.,-.01333,.5299,0., | |
591 | ||
592 | &0.,.1211,-.1739,0.,0.,.09469,-.07066,.01236, | |
593 | ||
594 | &-.1447,-.402,.1533,-.06479,6.7599,1.6596,.6798,-.8525, | |
595 | ||
596 | &0.,-4.4559,3.3756,-.9468, | |
597 | ||
598 | &0.,7.862,-3.6591,.03672,0.,-.2472,-.751,.0487, | |
599 | ||
600 | &3.017,-4.7347,3.3594,-.9443,0.,-.9342,.5454,-.1668, | |
601 | ||
602 | &5.304,1.4654,-1.4292,.7569,0.,-3.9141,2.8445,-.8411, | |
603 | ||
604 | &0.,9.0176,-10.426,4.0983,0.,-5.9602,7.515,-2.7329/ | |
605 | ||
606 | C...THE FOLLOWING DATA LINES ARE COEFFICIENTS NEEDED IN THE | |
607 | ||
608 | C...EICHTEN, HINCHLIFFE, LANE, QUIGG PROTON STRUCTURE FUNCTION | |
609 | ||
610 | C...POWERS OF 1-X IN DIFFERENT CASES | |
611 | ||
612 | DATA NEHLQ/3,4,7,5,7,7,7,7,3,4,7,6,7,7,7,7/ | |
613 | ||
614 | C...EXPANSION COEFFICIENTS FOR UP VALENCE QUARK DISTRIBUTION | |
615 | ||
616 | DATA (((CEHLQ(IX,IT,NX,1,1),IX=1,6),IT=1,6),NX=1,2)/ | |
617 | ||
618 | 1 7.677E-01,-2.087E-01,-3.303E-01,-2.517E-02,-1.570E-02,-1.000E-04, | |
619 | ||
620 | 2-5.326E-01,-2.661E-01, 3.201E-01, 1.192E-01, 2.434E-02, 7.620E-03, | |
621 | ||
622 | 3 2.162E-01, 1.881E-01,-8.375E-02,-6.515E-02,-1.743E-02,-5.040E-03, | |
623 | ||
624 | 4-9.211E-02,-9.952E-02, 1.373E-02, 2.506E-02, 8.770E-03, 2.550E-03, | |
625 | ||
626 | 5 3.670E-02, 4.409E-02, 9.600E-04,-7.960E-03,-3.420E-03,-1.050E-03, | |
627 | ||
628 | 6-1.549E-02,-2.026E-02,-3.060E-03, 2.220E-03, 1.240E-03, 4.100E-04, | |
629 | ||
630 | 1 2.395E-01, 2.905E-01, 9.778E-02, 2.149E-02, 3.440E-03, 5.000E-04, | |
631 | ||
632 | 2 1.751E-02,-6.090E-03,-2.687E-02,-1.916E-02,-7.970E-03,-2.750E-03, | |
633 | ||
634 | 3-5.760E-03,-5.040E-03, 1.080E-03, 2.490E-03, 1.530E-03, 7.500E-04, | |
635 | ||
636 | 4 1.740E-03, 1.960E-03, 3.000E-04,-3.400E-04,-2.900E-04,-1.800E-04, | |
637 | ||
638 | 5-5.300E-04,-6.400E-04,-1.700E-04, 4.000E-05, 6.000E-05, 4.000E-05, | |
639 | ||
640 | 6 1.700E-04, 2.200E-04, 8.000E-05, 1.000E-05,-1.000E-05,-1.000E-05/ | |
641 | ||
642 | DATA (((CEHLQ(IX,IT,NX,1,2),IX=1,6),IT=1,6),NX=1,2)/ | |
643 | ||
644 | 1 7.237E-01,-2.189E-01,-2.995E-01,-1.909E-02,-1.477E-02, 2.500E-04, | |
645 | ||
646 | 2-5.314E-01,-2.425E-01, 3.283E-01, 1.119E-01, 2.223E-02, 7.070E-03, | |
647 | ||
648 | 3 2.289E-01, 1.890E-01,-9.859E-02,-6.900E-02,-1.747E-02,-5.080E-03, | |
649 | ||
650 | 4-1.041E-01,-1.084E-01, 2.108E-02, 2.975E-02, 9.830E-03, 2.830E-03, | |
651 | ||
652 | 5 4.394E-02, 5.116E-02,-1.410E-03,-1.055E-02,-4.230E-03,-1.270E-03, | |
653 | ||
654 | 6-1.991E-02,-2.539E-02,-2.780E-03, 3.430E-03, 1.720E-03, 5.500E-04, | |
655 | ||
656 | 1 2.410E-01, 2.884E-01, 9.369E-02, 1.900E-02, 2.530E-03, 2.400E-04, | |
657 | ||
658 | 2 1.765E-02,-9.220E-03,-3.037E-02,-2.085E-02,-8.440E-03,-2.810E-03, | |
659 | ||
660 | 3-6.450E-03,-5.260E-03, 1.720E-03, 3.110E-03, 1.830E-03, 8.700E-04, | |
661 | ||
662 | 4 2.120E-03, 2.320E-03, 2.600E-04,-4.900E-04,-3.900E-04,-2.300E-04, | |
663 | ||
664 | 5-6.900E-04,-8.200E-04,-2.000E-04, 7.000E-05, 9.000E-05, 6.000E-05, | |
665 | ||
666 | 6 2.400E-04, 3.100E-04, 1.100E-04, 0.000E+00,-2.000E-05,-2.000E-05/ | |
667 | ||
668 | C...EXPANSION COEFFICIENTS FOR DOWN VALENCE QUARK DISTRIBUTION | |
669 | ||
670 | DATA (((CEHLQ(IX,IT,NX,2,1),IX=1,6),IT=1,6),NX=1,2)/ | |
671 | ||
672 | 1 3.813E-01,-8.090E-02,-1.634E-01,-2.185E-02,-8.430E-03,-6.200E-04, | |
673 | ||
674 | 2-2.948E-01,-1.435E-01, 1.665E-01, 6.638E-02, 1.473E-02, 4.080E-03, | |
675 | ||
676 | 3 1.252E-01, 1.042E-01,-4.722E-02,-3.683E-02,-1.038E-02,-2.860E-03, | |
677 | ||
678 | 4-5.478E-02,-5.678E-02, 8.900E-03, 1.484E-02, 5.340E-03, 1.520E-03, | |
679 | ||
680 | 5 2.220E-02, 2.567E-02,-3.000E-05,-4.970E-03,-2.160E-03,-6.500E-04, | |
681 | ||
682 | 6-9.530E-03,-1.204E-02,-1.510E-03, 1.510E-03, 8.300E-04, 2.700E-04, | |
683 | ||
684 | 1 1.261E-01, 1.354E-01, 3.958E-02, 8.240E-03, 1.660E-03, 4.500E-04, | |
685 | ||
686 | 2 3.890E-03,-1.159E-02,-1.625E-02,-9.610E-03,-3.710E-03,-1.260E-03, | |
687 | ||
688 | 3-1.910E-03,-5.600E-04, 1.590E-03, 1.590E-03, 8.400E-04, 3.900E-04, | |
689 | ||
690 | 4 6.400E-04, 4.900E-04,-1.500E-04,-2.900E-04,-1.800E-04,-1.000E-04, | |
691 | ||
692 | 5-2.000E-04,-1.900E-04, 0.000E+00, 6.000E-05, 4.000E-05, 3.000E-05, | |
693 | ||
694 | 6 7.000E-05, 8.000E-05, 2.000E-05,-1.000E-05,-1.000E-05,-1.000E-05/ | |
695 | ||
696 | DATA (((CEHLQ(IX,IT,NX,2,2),IX=1,6),IT=1,6),NX=1,2)/ | |
697 | ||
698 | 1 3.578E-01,-8.622E-02,-1.480E-01,-1.840E-02,-7.820E-03,-4.500E-04, | |
699 | ||
700 | 2-2.925E-01,-1.304E-01, 1.696E-01, 6.243E-02, 1.353E-02, 3.750E-03, | |
701 | ||
702 | 3 1.318E-01, 1.041E-01,-5.486E-02,-3.872E-02,-1.038E-02,-2.850E-03, | |
703 | ||
704 | 4-6.162E-02,-6.143E-02, 1.303E-02, 1.740E-02, 5.940E-03, 1.670E-03, | |
705 | ||
706 | 5 2.643E-02, 2.957E-02,-1.490E-03,-6.450E-03,-2.630E-03,-7.700E-04, | |
707 | ||
708 | 6-1.218E-02,-1.497E-02,-1.260E-03, 2.240E-03, 1.120E-03, 3.500E-04, | |
709 | ||
710 | 1 1.263E-01, 1.334E-01, 3.732E-02, 7.070E-03, 1.260E-03, 3.400E-04, | |
711 | ||
712 | 2 3.660E-03,-1.357E-02,-1.795E-02,-1.031E-02,-3.880E-03,-1.280E-03, | |
713 | ||
714 | 3-2.100E-03,-3.600E-04, 2.050E-03, 1.920E-03, 9.800E-04, 4.400E-04, | |
715 | ||
716 | 4 7.700E-04, 5.400E-04,-2.400E-04,-3.900E-04,-2.400E-04,-1.300E-04, | |
717 | ||
718 | 5-2.600E-04,-2.300E-04, 2.000E-05, 9.000E-05, 6.000E-05, 4.000E-05, | |
719 | ||
720 | 6 9.000E-05, 1.000E-04, 2.000E-05,-2.000E-05,-2.000E-05,-1.000E-05/ | |
721 | ||
722 | C...EXPANSION COEFFICIENTS FOR UP AND DOWN SEA QUARK DISTRIBUTIONS | |
723 | ||
724 | DATA (((CEHLQ(IX,IT,NX,3,1),IX=1,6),IT=1,6),NX=1,2)/ | |
725 | ||
726 | 1 6.870E-02,-6.861E-02, 2.973E-02,-5.400E-03, 3.780E-03,-9.700E-04, | |
727 | ||
728 | 2-1.802E-02, 1.400E-04, 6.490E-03,-8.540E-03, 1.220E-03,-1.750E-03, | |
729 | ||
730 | 3-4.650E-03, 1.480E-03,-5.930E-03, 6.000E-04,-1.030E-03,-8.000E-05, | |
731 | ||
732 | 4 6.440E-03, 2.570E-03, 2.830E-03, 1.150E-03, 7.100E-04, 3.300E-04, | |
733 | ||
734 | 5-3.930E-03,-2.540E-03,-1.160E-03,-7.700E-04,-3.600E-04,-1.900E-04, | |
735 | ||
736 | 6 2.340E-03, 1.930E-03, 5.300E-04, 3.700E-04, 1.600E-04, 9.000E-05, | |
737 | ||
738 | 1 1.014E+00,-1.106E+00, 3.374E-01,-7.444E-02, 8.850E-03,-8.700E-04, | |
739 | ||
740 | 2 9.233E-01,-1.285E+00, 4.475E-01,-9.786E-02, 1.419E-02,-1.120E-03, | |
741 | ||
742 | 3 4.888E-02,-1.271E-01, 8.606E-02,-2.608E-02, 4.780E-03,-6.000E-04, | |
743 | ||
744 | 4-2.691E-02, 4.887E-02,-1.771E-02, 1.620E-03, 2.500E-04,-6.000E-05, | |
745 | ||
746 | 5 7.040E-03,-1.113E-02, 1.590E-03, 7.000E-04,-2.000E-04, 0.000E+00, | |
747 | ||
748 | 6-1.710E-03, 2.290E-03, 3.800E-04,-3.500E-04, 4.000E-05, 1.000E-05/ | |
749 | ||
750 | DATA (((CEHLQ(IX,IT,NX,3,2),IX=1,6),IT=1,6),NX=1,2)/ | |
751 | ||
752 | 1 1.008E-01,-7.100E-02, 1.973E-02,-5.710E-03, 2.930E-03,-9.900E-04, | |
753 | ||
754 | 2-5.271E-02,-1.823E-02, 1.792E-02,-6.580E-03, 1.750E-03,-1.550E-03, | |
755 | ||
756 | 3 1.220E-02, 1.763E-02,-8.690E-03,-8.800E-04,-1.160E-03,-2.100E-04, | |
757 | ||
758 | 4-1.190E-03,-7.180E-03, 2.360E-03, 1.890E-03, 7.700E-04, 4.100E-04, | |
759 | ||
760 | 5-9.100E-04, 2.040E-03,-3.100E-04,-1.050E-03,-4.000E-04,-2.400E-04, | |
761 | ||
762 | 6 1.190E-03,-1.700E-04,-2.000E-04, 4.200E-04, 1.700E-04, 1.000E-04, | |
763 | ||
764 | 1 1.081E+00,-1.189E+00, 3.868E-01,-8.617E-02, 1.115E-02,-1.180E-03, | |
765 | ||
766 | 2 9.917E-01,-1.396E+00, 4.998E-01,-1.159E-01, 1.674E-02,-1.720E-03, | |
767 | ||
768 | 3 5.099E-02,-1.338E-01, 9.173E-02,-2.885E-02, 5.890E-03,-6.500E-04, | |
769 | ||
770 | 4-3.178E-02, 5.703E-02,-2.070E-02, 2.440E-03, 1.100E-04,-9.000E-05, | |
771 | ||
772 | 5 8.970E-03,-1.392E-02, 2.050E-03, 6.500E-04,-2.300E-04, 2.000E-05, | |
773 | ||
774 | 6-2.340E-03, 3.010E-03, 5.000E-04,-3.900E-04, 6.000E-05, 1.000E-05/ | |
775 | ||
776 | C...EXPANSION COEFFICIENTS FOR GLUON DISTRIBUTION | |
777 | ||
778 | DATA (((CEHLQ(IX,IT,NX,4,1),IX=1,6),IT=1,6),NX=1,2)/ | |
779 | ||
780 | 1 9.482E-01,-9.578E-01, 1.009E-01,-1.051E-01, 3.456E-02,-3.054E-02, | |
781 | ||
782 | 2-9.627E-01, 5.379E-01, 3.368E-01,-9.525E-02, 1.488E-02,-2.051E-02, | |
783 | ||
784 | 3 4.300E-01,-8.306E-02,-3.372E-01, 4.902E-02,-9.160E-03, 1.041E-02, | |
785 | ||
786 | 4-1.925E-01,-1.790E-02, 2.183E-01, 7.490E-03, 4.140E-03,-1.860E-03, | |
787 | ||
788 | 5 8.183E-02, 1.926E-02,-1.072E-01,-1.944E-02,-2.770E-03,-5.200E-04, | |
789 | ||
790 | 6-3.884E-02,-1.234E-02, 5.410E-02, 1.879E-02, 3.350E-03, 1.040E-03, | |
791 | ||
792 | 1 2.948E+01,-3.902E+01, 1.464E+01,-3.335E+00, 5.054E-01,-5.915E-02, | |
793 | ||
794 | 2 2.559E+01,-3.955E+01, 1.661E+01,-4.299E+00, 6.904E-01,-8.243E-02, | |
795 | ||
796 | 3-1.663E+00, 1.176E+00, 1.118E+00,-7.099E-01, 1.948E-01,-2.404E-02, | |
797 | ||
798 | 4-2.168E-01, 8.170E-01,-7.169E-01, 1.851E-01,-1.924E-02,-3.250E-03, | |
799 | ||
800 | 5 2.088E-01,-4.355E-01, 2.239E-01,-2.446E-02,-3.620E-03, 1.910E-03, | |
801 | ||
802 | 6-9.097E-02, 1.601E-01,-5.681E-02,-2.500E-03, 2.580E-03,-4.700E-04/ | |
803 | ||
804 | DATA (((CEHLQ(IX,IT,NX,4,2),IX=1,6),IT=1,6),NX=1,2)/ | |
805 | ||
806 | 1 2.367E+00, 4.453E-01, 3.660E-01, 9.467E-02, 1.341E-01, 1.661E-02, | |
807 | ||
808 | 2-3.170E+00,-1.795E+00, 3.313E-02,-2.874E-01,-9.827E-02,-7.119E-02, | |
809 | ||
810 | 3 1.823E+00, 1.457E+00,-2.465E-01, 3.739E-02, 6.090E-03, 1.814E-02, | |
811 | ||
812 | 4-1.033E+00,-9.827E-01, 2.136E-01, 1.169E-01, 5.001E-02, 1.684E-02, | |
813 | ||
814 | 5 5.133E-01, 5.259E-01,-1.173E-01,-1.139E-01,-4.988E-02,-2.021E-02, | |
815 | ||
816 | 6-2.881E-01,-3.145E-01, 5.667E-02, 9.161E-02, 4.568E-02, 1.951E-02, | |
817 | ||
818 | 1 3.036E+01,-4.062E+01, 1.578E+01,-3.699E+00, 6.020E-01,-7.031E-02, | |
819 | ||
820 | 2 2.700E+01,-4.167E+01, 1.770E+01,-4.804E+00, 7.862E-01,-1.060E-01, | |
821 | ||
822 | 3-1.909E+00, 1.357E+00, 1.127E+00,-7.181E-01, 2.232E-01,-2.481E-02, | |
823 | ||
824 | 4-2.488E-01, 9.781E-01,-8.127E-01, 2.094E-01,-2.997E-02,-4.710E-03, | |
825 | ||
826 | 5 2.506E-01,-5.427E-01, 2.672E-01,-3.103E-02,-1.800E-03, 2.870E-03, | |
827 | ||
828 | 6-1.128E-01, 2.087E-01,-6.972E-02,-2.480E-03, 2.630E-03,-8.400E-04/ | |
829 | ||
830 | C...EXPANSION COEFFICIENTS FOR STRANGE SEA QUARK DISTRIBUTION | |
831 | ||
832 | DATA (((CEHLQ(IX,IT,NX,5,1),IX=1,6),IT=1,6),NX=1,2)/ | |
833 | ||
834 | 1 4.968E-02,-4.173E-02, 2.102E-02,-3.270E-03, 3.240E-03,-6.700E-04, | |
835 | ||
836 | 2-6.150E-03,-1.294E-02, 6.740E-03,-6.890E-03, 9.000E-04,-1.510E-03, | |
837 | ||
838 | 3-8.580E-03, 5.050E-03,-4.900E-03,-1.600E-04,-9.400E-04,-1.500E-04, | |
839 | ||
840 | 4 7.840E-03, 1.510E-03, 2.220E-03, 1.400E-03, 7.000E-04, 3.500E-04, | |
841 | ||
842 | 5-4.410E-03,-2.220E-03,-8.900E-04,-8.500E-04,-3.600E-04,-2.000E-04, | |
843 | ||
844 | 6 2.520E-03, 1.840E-03, 4.100E-04, 3.900E-04, 1.600E-04, 9.000E-05, | |
845 | ||
846 | 1 9.235E-01,-1.085E+00, 3.464E-01,-7.210E-02, 9.140E-03,-9.100E-04, | |
847 | ||
848 | 2 9.315E-01,-1.274E+00, 4.512E-01,-9.775E-02, 1.380E-02,-1.310E-03, | |
849 | ||
850 | 3 4.739E-02,-1.296E-01, 8.482E-02,-2.642E-02, 4.760E-03,-5.700E-04, | |
851 | ||
852 | 4-2.653E-02, 4.953E-02,-1.735E-02, 1.750E-03, 2.800E-04,-6.000E-05, | |
853 | ||
854 | 5 6.940E-03,-1.132E-02, 1.480E-03, 6.500E-04,-2.100E-04, 0.000E+00, | |
855 | ||
856 | 6-1.680E-03, 2.340E-03, 4.200E-04,-3.400E-04, 5.000E-05, 1.000E-05/ | |
857 | ||
858 | DATA (((CEHLQ(IX,IT,NX,5,2),IX=1,6),IT=1,6),NX=1,2)/ | |
859 | ||
860 | 1 6.478E-02,-4.537E-02, 1.643E-02,-3.490E-03, 2.710E-03,-6.700E-04, | |
861 | ||
862 | 2-2.223E-02,-2.126E-02, 1.247E-02,-6.290E-03, 1.120E-03,-1.440E-03, | |
863 | ||
864 | 3-1.340E-03, 1.362E-02,-6.130E-03,-7.900E-04,-9.000E-04,-2.000E-04, | |
865 | ||
866 | 4 5.080E-03,-3.610E-03, 1.700E-03, 1.830E-03, 6.800E-04, 4.000E-04, | |
867 | ||
868 | 5-3.580E-03, 6.000E-05,-2.600E-04,-1.050E-03,-3.800E-04,-2.300E-04, | |
869 | ||
870 | 6 2.420E-03, 9.300E-04,-1.000E-04, 4.500E-04, 1.700E-04, 1.100E-04, | |
871 | ||
872 | 1 9.868E-01,-1.171E+00, 3.940E-01,-8.459E-02, 1.124E-02,-1.250E-03, | |
873 | ||
874 | 2 1.001E+00,-1.383E+00, 5.044E-01,-1.152E-01, 1.658E-02,-1.830E-03, | |
875 | ||
876 | 3 4.928E-02,-1.368E-01, 9.021E-02,-2.935E-02, 5.800E-03,-6.600E-04, | |
877 | ||
878 | 4-3.133E-02, 5.785E-02,-2.023E-02, 2.630E-03, 1.600E-04,-8.000E-05, | |
879 | ||
880 | 5 8.840E-03,-1.416E-02, 1.900E-03, 5.800E-04,-2.500E-04, 1.000E-05, | |
881 | ||
882 | 6-2.300E-03, 3.080E-03, 5.500E-04,-3.700E-04, 7.000E-05, 1.000E-05/ | |
883 | ||
884 | C...EXPANSION COEFFICIENTS FOR CHARM SEA QUARK DISTRIBUTION | |
885 | ||
886 | DATA (((CEHLQ(IX,IT,NX,6,1),IX=1,6),IT=1,6),NX=1,2)/ | |
887 | ||
888 | 1 9.270E-03,-1.817E-02, 9.590E-03,-6.390E-03, 1.690E-03,-1.540E-03, | |
889 | ||
890 | 2 5.710E-03,-1.188E-02, 6.090E-03,-4.650E-03, 1.240E-03,-1.310E-03, | |
891 | ||
892 | 3-3.960E-03, 7.100E-03,-3.590E-03, 1.840E-03,-3.900E-04, 3.400E-04, | |
893 | ||
894 | 4 1.120E-03,-1.960E-03, 1.120E-03,-4.800E-04, 1.000E-04,-4.000E-05, | |
895 | ||
896 | 5 4.000E-05,-3.000E-05,-1.800E-04, 9.000E-05,-5.000E-05,-2.000E-05, | |
897 | ||
898 | 6-4.200E-04, 7.300E-04,-1.600E-04, 5.000E-05, 5.000E-05, 5.000E-05, | |
899 | ||
900 | 1 8.098E-01,-1.042E+00, 3.398E-01,-6.824E-02, 8.760E-03,-9.000E-04, | |
901 | ||
902 | 2 8.961E-01,-1.217E+00, 4.339E-01,-9.287E-02, 1.304E-02,-1.290E-03, | |
903 | ||
904 | 3 3.058E-02,-1.040E-01, 7.604E-02,-2.415E-02, 4.600E-03,-5.000E-04, | |
905 | ||
906 | 4-2.451E-02, 4.432E-02,-1.651E-02, 1.430E-03, 1.200E-04,-1.000E-04, | |
907 | ||
908 | 5 1.122E-02,-1.457E-02, 2.680E-03, 5.800E-04,-1.200E-04, 3.000E-05, | |
909 | ||
910 | 6-7.730E-03, 7.330E-03,-7.600E-04,-2.400E-04, 1.000E-05, 0.000E+00/ | |
911 | ||
912 | DATA (((CEHLQ(IX,IT,NX,6,2),IX=1,6),IT=1,6),NX=1,2)/ | |
913 | ||
914 | 1 9.980E-03,-1.945E-02, 1.055E-02,-6.870E-03, 1.860E-03,-1.560E-03, | |
915 | ||
916 | 2 5.700E-03,-1.203E-02, 6.250E-03,-4.860E-03, 1.310E-03,-1.370E-03, | |
917 | ||
918 | 3-4.490E-03, 7.990E-03,-4.170E-03, 2.050E-03,-4.400E-04, 3.300E-04, | |
919 | ||
920 | 4 1.470E-03,-2.480E-03, 1.460E-03,-5.700E-04, 1.200E-04,-1.000E-05, | |
921 | ||
922 | 5-9.000E-05, 1.500E-04,-3.200E-04, 1.200E-04,-6.000E-05,-4.000E-05, | |
923 | ||
924 | 6-4.200E-04, 7.600E-04,-1.400E-04, 4.000E-05, 7.000E-05, 5.000E-05, | |
925 | ||
926 | 1 8.698E-01,-1.131E+00, 3.836E-01,-8.111E-02, 1.048E-02,-1.300E-03, | |
927 | ||
928 | 2 9.626E-01,-1.321E+00, 4.854E-01,-1.091E-01, 1.583E-02,-1.700E-03, | |
929 | ||
930 | 3 3.057E-02,-1.088E-01, 8.022E-02,-2.676E-02, 5.590E-03,-5.600E-04, | |
931 | ||
932 | 4-2.845E-02, 5.164E-02,-1.918E-02, 2.210E-03,-4.000E-05,-1.500E-04, | |
933 | ||
934 | 5 1.311E-02,-1.751E-02, 3.310E-03, 5.100E-04,-1.200E-04, 5.000E-05, | |
935 | ||
936 | 6-8.590E-03, 8.380E-03,-9.200E-04,-2.600E-04, 1.000E-05,-1.000E-05/ | |
937 | ||
938 | C...EXPANSION COEFFICIENTS FOR BOTTOM SEA QUARK DISTRIBUTION | |
939 | ||
940 | DATA (((CEHLQ(IX,IT,NX,7,1),IX=1,6),IT=1,6),NX=1,2)/ | |
941 | ||
942 | 1 9.010E-03,-1.401E-02, 7.150E-03,-4.130E-03, 1.260E-03,-1.040E-03, | |
943 | ||
944 | 2 6.280E-03,-9.320E-03, 4.780E-03,-2.890E-03, 9.100E-04,-8.200E-04, | |
945 | ||
946 | 3-2.930E-03, 4.090E-03,-1.890E-03, 7.600E-04,-2.300E-04, 1.400E-04, | |
947 | ||
948 | 4 3.900E-04,-1.200E-03, 4.400E-04,-2.500E-04, 2.000E-05,-2.000E-05, | |
949 | ||
950 | 5 2.600E-04, 1.400E-04,-8.000E-05, 1.000E-04, 1.000E-05, 1.000E-05, | |
951 | ||
952 | 6-2.600E-04, 3.200E-04, 1.000E-05,-1.000E-05, 1.000E-05,-1.000E-05, | |
953 | ||
954 | 1 8.029E-01,-1.075E+00, 3.792E-01,-7.843E-02, 1.007E-02,-1.090E-03, | |
955 | ||
956 | 2 7.903E-01,-1.099E+00, 4.153E-01,-9.301E-02, 1.317E-02,-1.410E-03, | |
957 | ||
958 | 3-1.704E-02,-1.130E-02, 2.882E-02,-1.341E-02, 3.040E-03,-3.600E-04, | |
959 | ||
960 | 4-7.200E-04, 7.230E-03,-5.160E-03, 1.080E-03,-5.000E-05,-4.000E-05, | |
961 | ||
962 | 5 3.050E-03,-4.610E-03, 1.660E-03,-1.300E-04,-1.000E-05, 1.000E-05, | |
963 | ||
964 | 6-4.360E-03, 5.230E-03,-1.610E-03, 2.000E-04,-2.000E-05, 0.000E+00/ | |
965 | ||
966 | DATA (((CEHLQ(IX,IT,NX,7,2),IX=1,6),IT=1,6),NX=1,2)/ | |
967 | ||
968 | 1 8.980E-03,-1.459E-02, 7.510E-03,-4.410E-03, 1.310E-03,-1.070E-03, | |
969 | ||
970 | 2 5.970E-03,-9.440E-03, 4.800E-03,-3.020E-03, 9.100E-04,-8.500E-04, | |
971 | ||
972 | 3-3.050E-03, 4.440E-03,-2.100E-03, 8.500E-04,-2.400E-04, 1.400E-04, | |
973 | ||
974 | 4 5.300E-04,-1.300E-03, 5.600E-04,-2.700E-04, 3.000E-05,-2.000E-05, | |
975 | ||
976 | 5 2.000E-04, 1.400E-04,-1.100E-04, 1.000E-04, 0.000E+00, 0.000E+00, | |
977 | ||
978 | 6-2.600E-04, 3.200E-04, 0.000E+00,-3.000E-05, 1.000E-05,-1.000E-05, | |
979 | ||
980 | 1 8.672E-01,-1.174E+00, 4.265E-01,-9.252E-02, 1.244E-02,-1.460E-03, | |
981 | ||
982 | 2 8.500E-01,-1.194E+00, 4.630E-01,-1.083E-01, 1.614E-02,-1.830E-03, | |
983 | ||
984 | 3-2.241E-02,-5.630E-03, 2.815E-02,-1.425E-02, 3.520E-03,-4.300E-04, | |
985 | ||
986 | 4-7.300E-04, 8.030E-03,-5.780E-03, 1.380E-03,-1.300E-04,-4.000E-05, | |
987 | ||
988 | 5 3.460E-03,-5.380E-03, 1.960E-03,-2.100E-04, 1.000E-05, 1.000E-05, | |
989 | ||
990 | 6-4.850E-03, 5.950E-03,-1.890E-03, 2.600E-04,-3.000E-05, 0.000E+00/ | |
991 | ||
992 | C...EXPANSION COEFFICIENTS FOR TOP SEA QUARK DISTRIBUTION | |
993 | ||
994 | DATA (((CEHLQ(IX,IT,NX,8,1),IX=1,6),IT=1,6),NX=1,2)/ | |
995 | ||
996 | 1 4.410E-03,-7.480E-03, 3.770E-03,-2.580E-03, 7.300E-04,-7.100E-04, | |
997 | ||
998 | 2 3.840E-03,-6.050E-03, 3.030E-03,-2.030E-03, 5.800E-04,-5.900E-04, | |
999 | ||
1000 | 3-8.800E-04, 1.660E-03,-7.500E-04, 4.700E-04,-1.000E-04, 1.000E-04, | |
1001 | ||
1002 | 4-8.000E-05,-1.500E-04, 1.200E-04,-9.000E-05, 3.000E-05, 0.000E+00, | |
1003 | ||
1004 | 5 1.300E-04,-2.200E-04,-2.000E-05,-2.000E-05,-2.000E-05,-2.000E-05, | |
1005 | ||
1006 | 6-7.000E-05, 1.900E-04,-4.000E-05, 2.000E-05, 0.000E+00, 0.000E+00, | |
1007 | ||
1008 | 1 6.623E-01,-9.248E-01, 3.519E-01,-7.930E-02, 1.110E-02,-1.180E-03, | |
1009 | ||
1010 | 2 6.380E-01,-9.062E-01, 3.582E-01,-8.479E-02, 1.265E-02,-1.390E-03, | |
1011 | ||
1012 | 3-2.581E-02, 2.125E-02, 4.190E-03,-4.980E-03, 1.490E-03,-2.100E-04, | |
1013 | ||
1014 | 4 7.100E-04, 5.300E-04,-1.270E-03, 3.900E-04,-5.000E-05,-1.000E-05, | |
1015 | ||
1016 | 5 3.850E-03,-5.060E-03, 1.860E-03,-3.500E-04, 4.000E-05, 0.000E+00, | |
1017 | ||
1018 | 6-3.530E-03, 4.460E-03,-1.500E-03, 2.700E-04,-3.000E-05, 0.000E+00/ | |
1019 | ||
1020 | DATA (((CEHLQ(IX,IT,NX,8,2),IX=1,6),IT=1,6),NX=1,2)/ | |
1021 | ||
1022 | 1 4.260E-03,-7.530E-03, 3.830E-03,-2.680E-03, 7.600E-04,-7.300E-04, | |
1023 | ||
1024 | 2 3.640E-03,-6.050E-03, 3.030E-03,-2.090E-03, 5.900E-04,-6.000E-04, | |
1025 | ||
1026 | 3-9.200E-04, 1.710E-03,-8.200E-04, 5.000E-04,-1.200E-04, 1.000E-04, | |
1027 | ||
1028 | 4-5.000E-05,-1.600E-04, 1.300E-04,-9.000E-05, 3.000E-05, 0.000E+00, | |
1029 | ||
1030 | 5 1.300E-04,-2.100E-04,-1.000E-05,-2.000E-05,-2.000E-05,-1.000E-05, | |
1031 | ||
1032 | 6-8.000E-05, 1.800E-04,-5.000E-05, 2.000E-05, 0.000E+00, 0.000E+00, | |
1033 | ||
1034 | 1 7.146E-01,-1.007E+00, 3.932E-01,-9.246E-02, 1.366E-02,-1.540E-03, | |
1035 | ||
1036 | 2 6.856E-01,-9.828E-01, 3.977E-01,-9.795E-02, 1.540E-02,-1.790E-03, | |
1037 | ||
1038 | 3-3.053E-02, 2.758E-02, 2.150E-03,-4.880E-03, 1.640E-03,-2.500E-04, | |
1039 | ||
1040 | 4 9.200E-04, 4.200E-04,-1.340E-03, 4.600E-04,-8.000E-05,-1.000E-05, | |
1041 | ||
1042 | 5 4.230E-03,-5.660E-03, 2.140E-03,-4.300E-04, 6.000E-05, 0.000E+00, | |
1043 | ||
1044 | 6-3.890E-03, 5.000E-03,-1.740E-03, 3.300E-04,-4.000E-05, 0.000E+00/ | |
1045 | ||
1046 | DATA TBMIN,TTMIN/8.1905,7.4474,11.5528,10.8097/ | |
1047 | ||
1048 | DATA XOLD,QOLD,IOLD,NOLD/-1.,0.,0,0/ | |
1049 | ||
1050 | DATA DMIN,Q0,QL/0.,2*2.,2*2.236,2.,.2,.4,.2,.29,.177/ | |
1051 | ||
1052 | IF (X.LE.ZERO) CALL HWWARN('HWSFUN',100,*999) | |
1053 | ||
1054 | XMWN=ONE-X | |
1055 | ||
1056 | IF (XMWN.LE.ZERO) THEN | |
1057 | ||
1058 | DO 1 I=1,13 | |
1059 | ||
1060 | DIST(I)=0 | |
1061 | ||
1062 | 1 CONTINUE | |
1063 | ||
1064 | RETURN | |
1065 | ||
1066 | ENDIF | |
1067 | ||
1068 | C---FREEZE THE SCALE IF REQUIRED | |
1069 | ||
1070 | SCALEF=SCALE | |
1071 | ||
1072 | IF (ISPAC.GT.0) SCALEF=MAX(SCALEF,QSPAC) | |
1073 | ||
1074 | C---CHECK IF PDFLIB REQUESTED | |
1075 | ||
1076 | IF (IBEAM.EQ.1.OR.IBEAM.EQ.2) THEN | |
1077 | ||
1078 | MPDF=MODPDF(IBEAM) | |
1079 | ||
1080 | ELSE | |
1081 | ||
1082 | MPDF=-1 | |
1083 | ||
1084 | ENDIF | |
1085 | ||
1086 | QSCA=ABS(SCALEF) | |
1087 | ||
1088 | IF (IDHAD.EQ.59.OR.IDHAD.EQ.71) THEN | |
1089 | ||
1090 | IF (MPDF.GE.0) THEN | |
1091 | ||
1092 | C---USE PDFLIB PHOTON STRUCTURE FUNCTIONS | |
1093 | ||
1094 | PARM(1)=AUTPDF(IBEAM) | |
1095 | ||
1096 | VAL(1)=FLOAT(MPDF) | |
1097 | ||
1098 | C---FIX TO CALL SCHULER-SJOSTRAND CODE | |
1099 | ||
1100 | IF (AUTPDF(IBEAM).EQ.'SaSph') THEN | |
1101 | ||
1102 | XSP=X | |
1103 | ||
1104 | IF ( XSP.LE.ZERO) CALL HWWARN('HWSFUN',102,*999) | |
1105 | ||
1106 | IF (ONE-XSP.LE.ZERO) CALL HWWARN('HWSFUN',103,*999) | |
1107 | ||
1108 | Q2=QSCA**2 | |
1109 | ||
1110 | ISET=MOD(MODPDF(IBEAM),10) | |
1111 | ||
1112 | IOP1=MOD(MODPDF(IBEAM)/10,2) | |
1113 | ||
1114 | IOP2=MOD(MODPDF(IBEAM)/20,2) | |
1115 | ||
1116 | IP2=MODPDF(IBEAM)/100 | |
1117 | ||
1118 | IF (IOP2.EQ.0) THEN | |
1119 | ||
1120 | P2=0. | |
1121 | ||
1122 | ELSE | |
1123 | ||
1124 | IHAD=IBEAM | |
1125 | ||
1126 | IF (JDAHEP(1,IHAD).NE.0) IHAD=JDAHEP(1,IHAD) | |
1127 | ||
1128 | P2=PHEP(5,IHAD)**2 | |
1129 | ||
1130 | ENDIF | |
1131 | ||
1132 | CALL SASGAM(ISET,XSP,Q2,P2,IP2,F2GM,XPGA) | |
1133 | ||
1134 | IF (IOP1.EQ.1 .AND. ISTAT.LT.10) THEN | |
1135 | ||
1136 | DO 5 I=-6,6 | |
1137 | ||
1138 | 5 XPGA(I)=XPVMD(I)+XPANL(I)+XPBEH(I)+XPDIR(I) | |
1139 | ||
1140 | ENDIF | |
1141 | ||
1142 | UPV=XPGA(2) | |
1143 | ||
1144 | DNV=XPGA(1) | |
1145 | ||
1146 | USEA=XPGA(2) | |
1147 | ||
1148 | DSEA=XPGA(1) | |
1149 | ||
1150 | STR=XPGA(3) | |
1151 | ||
1152 | CHM=XPGA(4) | |
1153 | ||
1154 | BTM=XPGA(5) | |
1155 | ||
1156 | TOP=XPGA(6) | |
1157 | ||
1158 | GLU=XPGA(0) | |
1159 | ||
1160 | ELSE | |
1161 | ||
1162 | CALL PDFSET(PARM,VAL) | |
1163 | ||
1164 | IF (X.LT.PDFXMN.AND.PDFWRX(IBEAM,1) .OR. | |
1165 | ||
1166 | & X.GT.PDFXMX.AND.PDFWRX(IBEAM,2)) THEN | |
1167 | ||
1168 | CALL HWWARN('HWSFUN',2,*999) | |
1169 | ||
1170 | WRITE (6,'(2A)') ' WARNING: PDFLIB CALLED WITH X', | |
1171 | ||
1172 | & ' OUTSIDE ALLOWED RANGE!' | |
1173 | ||
1174 | WRITE (6,'(1P,3(A,E9.3))') ' X VALUE=',X, | |
1175 | ||
1176 | & ', MINIMUM=',PDFXMN,', MAXIMUM=',PDFXMX | |
1177 | ||
1178 | WRITE (6,'(A)') ' NO FURTHER WARNINGS WILL BE ISSUED' | |
1179 | ||
1180 | IF (X.LT.PDFXMN) PDFWRX(IBEAM,1)=.FALSE. | |
1181 | ||
1182 | IF (X.GT.PDFXMX) PDFWRX(IBEAM,2)=.FALSE. | |
1183 | ||
1184 | ENDIF | |
1185 | ||
1186 | IF (QSCA**2.LT.PDFQMN.AND.PDFWRQ(IBEAM,1) .OR. | |
1187 | ||
1188 | & QSCA**2.GT.PDFQMX.AND.PDFWRQ(IBEAM,2)) THEN | |
1189 | ||
1190 | CALL HWWARN('HWSFUN',3,*999) | |
1191 | ||
1192 | WRITE (6,'(2A)') ' WARNING: PDFLIB CALLED WITH Q', | |
1193 | ||
1194 | & ' OUTSIDE ALLOWED RANGE!' | |
1195 | ||
1196 | WRITE (6,'(1P,3(A,E9.3))') ' Q VALUE=',QSCA, | |
1197 | ||
1198 | & ', MINIMUM=',SQRT(PDFQMN),', MAXIMUM=',SQRT(PDFQMX) | |
1199 | ||
1200 | WRITE (6,'(A)') ' NO FURTHER WARNINGS WILL BE ISSUED' | |
1201 | ||
1202 | IF (QSCA**2.LT.PDFQMN) PDFWRQ(IBEAM,1)=.FALSE. | |
1203 | ||
1204 | IF (QSCA**2.GT.PDFQMN) PDFWRQ(IBEAM,2)=.FALSE. | |
1205 | ||
1206 | ENDIF | |
1207 | ||
1208 | CALL STRUCTM(X,QSCA,UPV,DNV,USEA,DSEA,STR,CHM,BTM,TOP,GLU) | |
1209 | ||
1210 | ENDIF | |
1211 | ||
1212 | DIST(1)=DSEA | |
1213 | ||
1214 | DIST(2)=USEA | |
1215 | ||
1216 | DIST(7)=DSEA | |
1217 | ||
1218 | DIST(8)=USEA | |
1219 | ||
1220 | ELSE | |
1221 | ||
1222 | XSP=X | |
1223 | ||
1224 | IF ( XSP.LE.ZERO) CALL HWWARN('HWSFUN',102,*999) | |
1225 | ||
1226 | IF (ONE-XSP.LE.ZERO) CALL HWWARN('HWSFUN',103,*999) | |
1227 | ||
1228 | Q2=SCALEF**2 | |
1229 | ||
1230 | W2=Q2*(1-X)/X | |
1231 | ||
1232 | EMC2=4*RMASS(4)**2 | |
1233 | ||
1234 | EMB2=4*RMASS(5)**2 | |
1235 | ||
1236 | ALAM2=0.160 | |
1237 | ||
1238 | NFL=3 | |
1239 | ||
1240 | IF (Q2.GT.50.) NFL=4 | |
1241 | ||
1242 | IF (Q2.GT.500.) NFL=5 | |
1243 | ||
1244 | STR=HWSDGQ(XSP,Q2,NFL,1) | |
1245 | ||
1246 | CHM=HWSDGQ(XSP,Q2,NFL,2) | |
1247 | ||
1248 | GLU=HWSDGG(XSP,Q2,NFL) | |
1249 | ||
1250 | DIST(1)=STR | |
1251 | ||
1252 | DIST(2)=CHM | |
1253 | ||
1254 | DIST(7)=STR | |
1255 | ||
1256 | DIST(8)=CHM | |
1257 | ||
1258 | IF (W2.GT.EMB2) THEN | |
1259 | ||
1260 | BTM=STR | |
1261 | ||
1262 | IF (W2*ALAM2.LT.Q2*EMB2) | |
1263 | ||
1264 | & BTM=BTM*LOG(W2/EMB2)/LOG(Q2/ALAM2) | |
1265 | ||
1266 | ELSE | |
1267 | ||
1268 | BTM=0. | |
1269 | ||
1270 | ENDIF | |
1271 | ||
1272 | IF (W2.GT.EMC2) THEN | |
1273 | ||
1274 | IF (W2*ALAM2.LT.Q2*EMC2) | |
1275 | ||
1276 | & CHM=CHM*LOG(W2/EMC2)/LOG(Q2/ALAM2) | |
1277 | ||
1278 | ELSE | |
1279 | ||
1280 | CHM=0. | |
1281 | ||
1282 | ENDIF | |
1283 | ||
1284 | TOP=0. | |
1285 | ||
1286 | ENDIF | |
1287 | ||
1288 | C---INCLUDE SUPPRESSION FROM PHOTON VIRTUALITY IF NECESSARY | |
1289 | ||
1290 | IF (PHOMAS.GT.ZERO.AND.(IBEAM.EQ.1.OR.IBEAM.EQ.2)) THEN | |
1291 | ||
1292 | IHAD=IBEAM | |
1293 | ||
1294 | IF (JDAHEP(1,IHAD).NE.0) IHAD=JDAHEP(1,IHAD) | |
1295 | ||
1296 | IF (IDHW(IHAD).EQ.59) THEN | |
1297 | ||
1298 | FAC=LOG((QSCA**2+PHOMAS**2)/(PHEP(5,IHAD)**2+PHOMAS**2))/ | |
1299 | ||
1300 | $ LOG((QSCA**2+PHOMAS**2)/( PHOMAS**2)) | |
1301 | ||
1302 | IF (FAC.LT.ZERO) FAC=ZERO | |
1303 | ||
1304 | DIST(1)=DIST(1)*FAC | |
1305 | ||
1306 | DIST(2)=DIST(2)*FAC | |
1307 | ||
1308 | DIST(7)=DIST(7)*FAC | |
1309 | ||
1310 | DIST(8)=DIST(8)*FAC | |
1311 | ||
1312 | STR=STR*FAC | |
1313 | ||
1314 | CHM=CHM*FAC | |
1315 | ||
1316 | BTM=BTM*FAC | |
1317 | ||
1318 | TOP=TOP*FAC | |
1319 | ||
1320 | GLU=GLU*FAC**2 | |
1321 | ||
1322 | ELSE | |
1323 | ||
1324 | CALL HWWARN('HWSFUN',1,*999) | |
1325 | ||
1326 | ENDIF | |
1327 | ||
1328 | ENDIF | |
1329 | ||
1330 | GOTO 900 | |
1331 | ||
1332 | ENDIF | |
1333 | ||
1334 | IF (MPDF.GE.0) THEN | |
1335 | ||
1336 | C---USE PDFLIB NUCLEON STRUCTURE FUNCTIONS | |
1337 | ||
1338 | PARM(1)=AUTPDF(IBEAM) | |
1339 | ||
1340 | VAL(1)=FLOAT(MPDF) | |
1341 | ||
1342 | CALL PDFSET(PARM,VAL) | |
1343 | ||
1344 | IF (X.LT.PDFXMN.AND.PDFWRX(IBEAM,1) .OR. | |
1345 | ||
1346 | & X.GT.PDFXMX.AND.PDFWRX(IBEAM,2)) THEN | |
1347 | ||
1348 | CALL HWWARN('HWSFUN',4,*999) | |
1349 | ||
1350 | WRITE (6,'(2A)') ' WARNING: PDFLIB CALLED WITH X', | |
1351 | ||
1352 | & ' OUTSIDE ALLOWED RANGE!' | |
1353 | ||
1354 | WRITE (6,'(1P,3(A,E9.3))') ' X VALUE=',X, | |
1355 | ||
1356 | & ', MINIMUM=',PDFXMN,', MAXIMUM=',PDFXMX | |
1357 | ||
1358 | WRITE (6,'(A)') ' NO FURTHER WARNINGS WILL BE ISSUED' | |
1359 | ||
1360 | IF (X.LT.PDFXMN) PDFWRX(IBEAM,1)=.FALSE. | |
1361 | ||
1362 | IF (X.GT.PDFXMX) PDFWRX(IBEAM,2)=.FALSE. | |
1363 | ||
1364 | ENDIF | |
1365 | ||
1366 | IF (QSCA**2.LT.PDFQMN.AND.PDFWRQ(IBEAM,1) .OR. | |
1367 | ||
1368 | & QSCA**2.GT.PDFQMX.AND.PDFWRQ(IBEAM,2)) THEN | |
1369 | ||
1370 | CALL HWWARN('HWSFUN',5,*999) | |
1371 | ||
1372 | WRITE (6,'(2A)') ' WARNING: PDFLIB CALLED WITH Q', | |
1373 | ||
1374 | & ' OUTSIDE ALLOWED RANGE!' | |
1375 | ||
1376 | WRITE (6,'(1P,3(A,E9.3))') ' Q VALUE=',QSCA, | |
1377 | ||
1378 | & ', MINIMUM=',SQRT(PDFQMN),', MAXIMUM=',SQRT(PDFQMX) | |
1379 | ||
1380 | WRITE (6,'(A)') ' NO FURTHER WARNINGS WILL BE ISSUED' | |
1381 | ||
1382 | IF (QSCA**2.LT.PDFQMN) PDFWRQ(IBEAM,1)=.FALSE. | |
1383 | ||
1384 | IF (QSCA**2.GT.PDFQMN) PDFWRQ(IBEAM,2)=.FALSE. | |
1385 | ||
1386 | ENDIF | |
1387 | ||
1388 | CALL STRUCTM(X,QSCA,UPV,DNV,USEA,DSEA,STR,CHM,BTM,TOP,GLU) | |
1389 | ||
1390 | ELSE | |
1391 | ||
1392 | IF (QSCA.LT.Q0(NSET)) QSCA=Q0(NSET) | |
1393 | ||
1394 | IF (QSCA.NE.QOLD.OR.IDHAD.NE.IOLD.OR.NSET.NE.NOLD) THEN | |
1395 | ||
1396 | C---INITIALIZE | |
1397 | ||
1398 | IF (NSET.LT.1.OR.NSET.GT.5) CALL HWWARN('HWSFUN',400,*999) | |
1399 | ||
1400 | QOLD=QSCA | |
1401 | ||
1402 | IOLD=IDHAD | |
1403 | ||
1404 | NOLD=NSET | |
1405 | ||
1406 | SS=LOG(QSCA/QL(NSET)) | |
1407 | ||
1408 | SMIN=LOG(Q0(NSET)/QL(NSET)) | |
1409 | ||
1410 | IF (NSET.LT.3.OR.NSET.EQ.5) THEN | |
1411 | ||
1412 | S=LOG(SS/SMIN) | |
1413 | ||
1414 | ELSE | |
1415 | ||
1416 | T=2.*SS | |
1417 | ||
1418 | TMIN=2.*SMIN | |
1419 | ||
1420 | TMAX=2.*LOG(1.E4/QL(NSET)) | |
1421 | ||
1422 | ENDIF | |
1423 | ||
1424 | IF (IDHAD.GE.72) THEN | |
1425 | ||
1426 | IF (NSET.LT.3) THEN | |
1427 | ||
1428 | IP=NSET | |
1429 | ||
1430 | DO 10 I=1,5 | |
1431 | ||
1432 | DO 10 J=1,6 | |
1433 | ||
1434 | 10 A(J,I)=B(1,J,I,IP)+S*(B(2,J,I,IP)+S*B(3,J,I,IP)) | |
1435 | ||
1436 | DO 20 K=1,2 | |
1437 | ||
1438 | AA=ONE+A(2,K)+A(3,K) | |
1439 | ||
1440 | 20 G(K)=HWSGAM(AA)/((ONE+A(2,K)*A(4,K)/AA)*HWSGAM(A(2,K)) | |
1441 | ||
1442 | & *HWSGAM(ONE+A(3,K))) | |
1443 | ||
1444 | ELSEIF (NSET.EQ.5) THEN | |
1445 | ||
1446 | DO 21 I=1,5 | |
1447 | ||
1448 | DO 21 J=1,6 | |
1449 | ||
1450 | 21 A(J,I)=BB(1,J,I)+S*(BB(2,J,I)+S*(BB(3,J,I)+S*BB(4,J,I))) | |
1451 | ||
1452 | DO 22 K=1,2 | |
1453 | ||
1454 | AA=ONE+A(2,K)+A(3,K) | |
1455 | ||
1456 | 22 G(K)=HWSGAM(AA)/((ONE+A(2,K)/AA*(A(4,K)+ | |
1457 | ||
1458 | & (ONE+A(2,K))/(ONE+AA)*A(5,K)))*HWSGAM(A(2,K)) | |
1459 | ||
1460 | & *HWSGAM(ONE+A(3,K))) | |
1461 | ||
1462 | ELSE | |
1463 | ||
1464 | IP=NSET-2 | |
1465 | ||
1466 | VT=MAX(-ONE,MIN(ONE,(2.*T-TMAX-TMIN)/(TMAX-TMIN))) | |
1467 | ||
1468 | WT=VT*VT | |
1469 | ||
1470 | C...CHEBYSHEV POLYNOMIALS FOR T EXPANSION | |
1471 | ||
1472 | TT(1)=1. | |
1473 | ||
1474 | TT(2)=VT | |
1475 | ||
1476 | TT(3)= 2.*WT- 1. | |
1477 | ||
1478 | TT(4)= (4.*WT- 3.)*VT | |
1479 | ||
1480 | TT(5)= (8.*WT- 8.)*WT+1. | |
1481 | ||
1482 | TT(6)=((16.*WT-20.)*WT+5.)*VT | |
1483 | ||
1484 | ENDIF | |
1485 | ||
1486 | ELSEIF (NSET.LT.3) THEN | |
1487 | ||
1488 | IP=NSET+2 | |
1489 | ||
1490 | DO 30 I=1,5 | |
1491 | ||
1492 | DO 30 J=1,6 | |
1493 | ||
1494 | 30 A(J,I)=B(1,J,I,IP)+S*(B(2,J,I,IP)+S*B(3,J,I,IP)) | |
1495 | ||
1496 | AA=ONE+A(2,1)+A(3,1) | |
1497 | ||
1498 | G(1)=HWSGAM(AA)/(HWSGAM(A(2,1))*HWSGAM(ONE+A(3,1))) | |
1499 | ||
1500 | G(2)=0. | |
1501 | ||
1502 | ENDIF | |
1503 | ||
1504 | ENDIF | |
1505 | ||
1506 | C | |
1507 | ||
1508 | IF (NSET.LT.3.OR.NSET.EQ.5) THEN | |
1509 | ||
1510 | DO 50 I=1,5 | |
1511 | ||
1512 | 50 F(I)=A(1,I)*X**A(2,I)*XMWN**A(3,I)*(ONE+X* | |
1513 | ||
1514 | & (A(4,I)+X*(A(5,I) + X*A(6,I)))) | |
1515 | ||
1516 | F(1)=F(1)*G(1) | |
1517 | ||
1518 | F(2)=F(2)*G(2) | |
1519 | ||
1520 | UPV=F(1)-F(2) | |
1521 | ||
1522 | DNV=F(2) | |
1523 | ||
1524 | SEA=F(3)/6. | |
1525 | ||
1526 | STR=SEA | |
1527 | ||
1528 | CHM=F(4) | |
1529 | ||
1530 | BTM=0. | |
1531 | ||
1532 | TOP=0. | |
1533 | ||
1534 | GLU=F(5) | |
1535 | ||
1536 | ELSE | |
1537 | ||
1538 | IF (X.NE.XOLD) THEN | |
1539 | ||
1540 | XOLD=X | |
1541 | ||
1542 | IF (X.GT.0.1) THEN | |
1543 | ||
1544 | NX=1 | |
1545 | ||
1546 | VX=(2.*X-1.1)/0.9 | |
1547 | ||
1548 | ELSE | |
1549 | ||
1550 | NX=2 | |
1551 | ||
1552 | VX=MAX(-ONE,(2.*LOG(X)+11.51293)/6.90776) | |
1553 | ||
1554 | ENDIF | |
1555 | ||
1556 | WX=VX*VX | |
1557 | ||
1558 | TX(1)=1. | |
1559 | ||
1560 | TX(2)=VX | |
1561 | ||
1562 | TX(3)= 2.*WX- 1. | |
1563 | ||
1564 | TX(4)= (4.*WX- 3.)*VX | |
1565 | ||
1566 | TX(5)= (8.*WX- 8.)*WX+1. | |
1567 | ||
1568 | TX(6)=((16.*WX-20.)*WX+5.)*VX | |
1569 | ||
1570 | ENDIF | |
1571 | ||
1572 | C...CALCULATE STRUCTURE FUNCTIONS | |
1573 | ||
1574 | DO 120 IFL=1,6 | |
1575 | ||
1576 | XQSUM=0. | |
1577 | ||
1578 | DO 110 IT=1,6 | |
1579 | ||
1580 | DO 110 IX=1,6 | |
1581 | ||
1582 | 110 XQSUM=XQSUM+CEHLQ(IX,IT,NX,IFL,IP)*TX(IX)*TT(IT) | |
1583 | ||
1584 | 120 XQ(IFL)=XQSUM*XMWN**NEHLQ(IFL,IP) | |
1585 | ||
1586 | UPV=XQ(1) | |
1587 | ||
1588 | DNV=XQ(2) | |
1589 | ||
1590 | STR=XQ(5) | |
1591 | ||
1592 | CHM=XQ(6) | |
1593 | ||
1594 | SEA=XQ(3) | |
1595 | ||
1596 | GLU=XQ(4) | |
1597 | ||
1598 | C...SPECIAL EXPANSION FOR BOTTOM (THRESHOLD EFFECTS) | |
1599 | ||
1600 | IF (NFLAV.LT.5.OR.T.LE.TBMIN(IP)) THEN | |
1601 | ||
1602 | BTM=0. | |
1603 | ||
1604 | ELSE | |
1605 | ||
1606 | VT=MAX(-ONE,MIN(ONE,(2.*T-TMAX-TBMIN(IP))/(TMAX-TBMIN(IP)))) | |
1607 | ||
1608 | WT=VT*VT | |
1609 | ||
1610 | TB(1)=1. | |
1611 | ||
1612 | TB(2)=VT | |
1613 | ||
1614 | TB(3)= 2.*WT- 1. | |
1615 | ||
1616 | TB(4)= (4.*WT- 3.)*VT | |
1617 | ||
1618 | TB(5)= (8.*WT- 8.)*WT+1. | |
1619 | ||
1620 | TB(6)=((16.*WT-20.)*WT+5.)*VT | |
1621 | ||
1622 | XQSUM=0. | |
1623 | ||
1624 | DO 130 IT=1,6 | |
1625 | ||
1626 | DO 130 IX=1,6 | |
1627 | ||
1628 | 130 XQSUM=XQSUM+CEHLQ(IX,IT,NX,7,IP)*TX(IX)*TB(IT) | |
1629 | ||
1630 | BTM=XQSUM*XMWN**NEHLQ(7,IP) | |
1631 | ||
1632 | ENDIF | |
1633 | ||
1634 | C...SPECIAL EXPANSION FOR TOP (THRESHOLD EFFECTS) | |
1635 | ||
1636 | TPMIN=TTMIN(IP)+TMTOP | |
1637 | ||
1638 | C---TMTOP=2.*LOG(TOPMAS/30.) | |
1639 | ||
1640 | TPMAX=TMAX+TMTOP | |
1641 | ||
1642 | IF (NFLAV.LT.6.OR.T.LE.TPMIN) THEN | |
1643 | ||
1644 | TOP=0. | |
1645 | ||
1646 | ELSE | |
1647 | ||
1648 | VT=MAX(-ONE,MIN(ONE,(2.*T-TPMAX-TPMIN)/(TPMAX-TPMIN))) | |
1649 | ||
1650 | WT=VT*VT | |
1651 | ||
1652 | TB(1)=1. | |
1653 | ||
1654 | TB(2)=VT | |
1655 | ||
1656 | TB(3)= 2.*WT- 1. | |
1657 | ||
1658 | TB(4)= (4.*WT- 3.)*VT | |
1659 | ||
1660 | TB(5)= (8.*WT- 8.)*WT+1. | |
1661 | ||
1662 | TB(6)=((16.*WT-20.)*WT+5.)*VT | |
1663 | ||
1664 | XQSUM=0. | |
1665 | ||
1666 | DO 150 IT=1,6 | |
1667 | ||
1668 | DO 150 IX=1,6 | |
1669 | ||
1670 | 150 XQSUM=XQSUM+CEHLQ(IX,IT,NX,8,IP)*TX(IX)*TB(IT) | |
1671 | ||
1672 | TOP=XQSUM*XMWN**NEHLQ(8,IP) | |
1673 | ||
1674 | ENDIF | |
1675 | ||
1676 | ENDIF | |
1677 | ||
1678 | ENDIF | |
1679 | ||
1680 | IF (MPDF.LT.0) THEN | |
1681 | ||
1682 | USEA=SEA | |
1683 | ||
1684 | DSEA=USEA | |
1685 | ||
1686 | ENDIF | |
1687 | ||
1688 | IF (IDHAD.EQ.73.OR.IDHAD.EQ.72) THEN | |
1689 | ||
1690 | DIST(1)=DSEA+DNV | |
1691 | ||
1692 | DIST(2)=USEA+UPV | |
1693 | ||
1694 | DIST(7)=DSEA | |
1695 | ||
1696 | DIST(8)=USEA | |
1697 | ||
1698 | ELSEIF (IDHAD.EQ.91) THEN | |
1699 | ||
1700 | DIST(1)=DSEA | |
1701 | ||
1702 | DIST(2)=USEA | |
1703 | ||
1704 | DIST(7)=DSEA+DNV | |
1705 | ||
1706 | DIST(8)=USEA+UPV | |
1707 | ||
1708 | ELSEIF (IDHAD.EQ.75) THEN | |
1709 | ||
1710 | DIST(1)=USEA+UPV | |
1711 | ||
1712 | DIST(2)=DSEA+DNV | |
1713 | ||
1714 | DIST(7)=USEA | |
1715 | ||
1716 | DIST(8)=DSEA | |
1717 | ||
1718 | ELSEIF (IDHAD.EQ.93) THEN | |
1719 | ||
1720 | DIST(1)=USEA | |
1721 | ||
1722 | DIST(2)=DSEA | |
1723 | ||
1724 | DIST(7)=USEA+UPV | |
1725 | ||
1726 | DIST(8)=DSEA+DNV | |
1727 | ||
1728 | ELSEIF (IDHAD.EQ.38) THEN | |
1729 | ||
1730 | DIST(1)=USEA | |
1731 | ||
1732 | DIST(2)=USEA+UPV | |
1733 | ||
1734 | DIST(7)=USEA+UPV | |
1735 | ||
1736 | DIST(8)=USEA | |
1737 | ||
1738 | ELSEIF (IDHAD.EQ.30) THEN | |
1739 | ||
1740 | DIST(1)=USEA+UPV | |
1741 | ||
1742 | DIST(2)=USEA | |
1743 | ||
1744 | DIST(7)=USEA | |
1745 | ||
1746 | DIST(8)=USEA+UPV | |
1747 | ||
1748 | ELSE | |
1749 | ||
1750 | PRINT *,' CALLED HWSFUN FOR IDHAD =',IDHAD | |
1751 | ||
1752 | CALL HWWARN('HWSFUN',400,*999) | |
1753 | ||
1754 | ENDIF | |
1755 | ||
1756 | 900 DIST(3)=STR | |
1757 | ||
1758 | DIST(4)=CHM | |
1759 | ||
1760 | DIST(5)=BTM | |
1761 | ||
1762 | DIST(6)=TOP | |
1763 | ||
1764 | DIST(9)=STR | |
1765 | ||
1766 | DIST(10)=CHM | |
1767 | ||
1768 | DIST(11)=BTM | |
1769 | ||
1770 | DIST(12)=TOP | |
1771 | ||
1772 | DIST(13)=GLU | |
1773 | ||
1774 | DO 901 I=1,13 | |
1775 | ||
1776 | IF (DIST(I).LT.DMIN) DIST(I)=DMIN | |
1777 | ||
1778 | 901 CONTINUE | |
1779 | ||
1780 | C---FOR REMNANT NUCLEONS SWITCH OFF VALENCE QUARKS, | |
1781 | ||
1782 | C WHILE MAINTAINING MOMENTUM SUM RULE | |
1783 | ||
1784 | IF (IDHAD.EQ.72) THEN | |
1785 | ||
1786 | TOTAL=0 | |
1787 | ||
1788 | DO 910 I=1,13 | |
1789 | ||
1790 | TOTAL=TOTAL+DIST(I) | |
1791 | ||
1792 | 910 CONTINUE | |
1793 | ||
1794 | DIST(1)=DIST(1)-DNV | |
1795 | ||
1796 | DIST(2)=DIST(2)-UPV | |
1797 | ||
1798 | IF (TOTAL.GT.DNV+UPV) THEN | |
1799 | ||
1800 | DO 920 I=1,13 | |
1801 | ||
1802 | DIST(I)=DIST(I)*TOTAL/(TOTAL-DNV-UPV) | |
1803 | ||
1804 | 920 CONTINUE | |
1805 | ||
1806 | ENDIF | |
1807 | ||
1808 | ENDIF | |
1809 | ||
1810 | 999 END | |
1811 | ||
1812 | CDECK ID>, HWSGAM. | |
1813 | ||
1814 | *CMZ :- -26/04/91 11.11.56 by Bryan Webber | |
1815 | ||
1816 | *-- Author : Adapted by Bryan Webber | |
1817 | ||
1818 | C----------------------------------------------------------------------- | |
1819 | ||
1820 | FUNCTION HWSGAM(ZINPUT) | |
1821 | ||
1822 | C----------------------------------------------------------------------- | |
1823 | ||
1824 | C Gamma function computed by eq. 6.1.40, Abramowitz. | |
1825 | ||
1826 | C B(M) = B2m/(2m *(2m-1)) where B2m is the 2m'th Bernoulli number. | |
1827 | ||
1828 | C HLNTPI = .5*LOG(2.*PI) | |
1829 | ||
1830 | C----------------------------------------------------------------------- | |
1831 | ||
1832 | DOUBLE PRECISION HWSGAM,ZINPUT,B(10),HLNTPI,Z,SHIFT,G,T,RECZSQ | |
1833 | ||
1834 | INTEGER I | |
1835 | ||
1836 | DATA B/ | |
1837 | ||
1838 | 1 0.83333333333333333333D-01, -0.27777777777777777778D-02, | |
1839 | ||
1840 | 1 0.79365079365079365079D-03, -0.59523809523809523810D-03, | |
1841 | ||
1842 | 1 0.84175084175084175084D-03, -0.19175269175269175269D-02, | |
1843 | ||
1844 | 1 0.64102564102564102564D-02, -0.29550653594771241830D-01, | |
1845 | ||
1846 | 1 0.17964437236883057316D0 , -1.3924322169059011164D0 / | |
1847 | ||
1848 | DATA HLNTPI/0.91893853320467274178D0/ | |
1849 | ||
1850 | C | |
1851 | ||
1852 | C Shift argument to large value ( > 20 ) | |
1853 | ||
1854 | C | |
1855 | ||
1856 | Z=ZINPUT | |
1857 | ||
1858 | SHIFT=1. | |
1859 | ||
1860 | 10 IF (Z.LT.20.D0) THEN | |
1861 | ||
1862 | SHIFT = SHIFT*Z | |
1863 | ||
1864 | Z = Z + 1.D0 | |
1865 | ||
1866 | GOTO 10 | |
1867 | ||
1868 | ENDIF | |
1869 | ||
1870 | C | |
1871 | ||
1872 | C Compute asymptotic formula | |
1873 | ||
1874 | C | |
1875 | ||
1876 | G = (Z-.5D0)*LOG(Z) - Z + HLNTPI | |
1877 | ||
1878 | T = 1.D0/Z | |
1879 | ||
1880 | RECZSQ = T**2 | |
1881 | ||
1882 | DO 20 I = 1,10 | |
1883 | ||
1884 | G = G + B(I)*T | |
1885 | ||
1886 | T = T*RECZSQ | |
1887 | ||
1888 | 20 CONTINUE | |
1889 | ||
1890 | HWSGAM = EXP(G)/SHIFT | |
1891 | ||
1892 | END |