VBF: coefmom.f

File coefmom.f, 72.0 KB (added by trac, 12 years ago)
Line 
1
2 PROGRAM MAIN
3C THIS PROGRAM COMPUTES THE MOMENTS OF THE COEFFICIENT FUNCTIONS
4C THE COEFFICIENT FUNCTIONS ARE PRESENTED IN THE DIS AND THE MSBAR SCHEME .
5
6 IMPLICIT REAL*8(A-H,O-Z)
7 DOUBLE PRECISION LAMBDA,NF
8 COMMON/PAR/NF,Q2,FS,RS,DN
9 COMMON/PAR1/LAMBDA,Q02
10 COMMON/DIS/NDIS
11 COMMON/CHOICE/MODE
12 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
13 EXTERNAL F30R,F31R,F31S,F32R,F32S
14 EXTERNAL F20R,F21R,F21S,F22R,F22S
15 EXTERNAL F1G1,F1G2,F1Q2
16csm
17 external FG21R,FG22R,FPS22R,FL1R,FL2R,FPSL2R,FLG1R,FLG2R
18
19C CHOOSE THE MOMENT DN=N M^(N) = \int_0^1 X^N C_{i,k}(X,Q^2/\MU^2)
20
21 write(*,*)
22 # 'enter number of flavours nf'
23 read (*,*) NF
24 write(*,*)
25 # 'enter Mellin moment N'
26 read (*,*) DNNN
27
28 DN = DNNN - 1.d0
29csm DN=5.0D0
30
31C CHOOSE FA=1.0D0 IF ONE TAKES THE ZEROTH MOMENT DN=0.0D0
32C CHOOSE FA=2.0D0 IF ONE TAKES DN > 0.0D0
33
34 FA=2.0D0
35 KEY=2
36 EPS=1.0D-10
37 MAX=10000
38C HERE FOLLOW THE COLOUR FACTORS CF, CF^2, CA.CF, CF.T_f
39
40 CF=4.0D0/3.0D0
41C CF=0.0D0
42 CA=3.0D0
43C CA=0.0D0
44 CF2=16.0D0/9.0D0
45C CF2=0.0D0
46 CACF=4.0D0
47C CACF=0.0D0
48 CFTF=2.0D0/3.0D0
49C CFTF=0.0D0
50
51C CHOOSE DIS OR DISP SCHEME
52C DISP SCHEME IS A SCHEME IN WHICH THE MOMENTUM CONSERVATION SUM RULE
53C IS SATISFIED
54
55C DISP SCHEME : NDIS=1
56C DIS SCHEME : NDIS=2
57
58 NDIS=0
59
60 PI=3.14159265359D0
61
62C-------------------------------------------------------------------
63C HERE FOLLOWS THE VIRTUAL PHOTON MOMENTUM IN GEV
64 Q=DSQRT(10.0D0)
65csm WRITE(*,21) Q
66 21 FORMAT(' ','THE VIRTUAL PHOTON MOMENTUM=',D20.10)
67 Q2=Q*Q
68C HERE FOLLOWS THE FACTORIZATION SCALE IN GEV \MU=FS
69 FS=Q
70csm WRITE(*,23) FS
71 23 FORMAT(' ','THE FACTORIZATION SCALE=',D20.10)
72C HERE FOLLOWS THE RENORMALIZATION SCALE IN GEV
73 RS=Q
74csm WRITE(*,24) RS
75 24 FORMAT(' ','THE RENORMALIZATION SCALE=',D20.10)
76
77C HERE FOLLOWS THE NUMBER OF FLAVOURS
78csm NF=4.0D0
79 WRITE(*,25) NF
80 25 FORMAT(' ','nf = ',g20.2)
81
82 DNN = DN + 1.0d0
83 WRITE(*,325) DNN
84 325 FORMAT(' ','N = ',g20.2)
85
86 R20R=DAIND(0.0D0,1.0D0,F20R,EPS,KEY,MAX,KOUNT,EST)
87 WRITE(*,29) R20R
88 29 FORMAT(' ','c2qq0 = ',D20.10)
89 R21R=DAIND(0.0D0,1.0D0,F21R,EPS,KEY,MAX,KOUNT,EST)
90 IF(FA.GT.1.1D0) THEN
91 R21S=DAIND(0.0D0,1.0D0,F21S,EPS,KEY,MAX,KOUNT,EST)
92 ELSE
93 R21S=0.0D0
94 ENDIF
95 F2QR1=R21R
96 F2QL1=R21S
97 F2QVS1=COR1(0.0D0)
98 F2Q1=F2QR1+F2QL1+F2QVS1
99 WRITE(*,33) F2Q1
100 33 FORMAT(' ','c2qq1 = ',D20.10)
101csm beg added
102 RG21R=DAIND(0.0D0,1.0D0,FG21R,EPS,KEY,MAX,KOUNT,EST)
103 RG21S=0.0D0
104 F2GR1=RG21R
105 F2GL1=RG21S
106 F2G1=F2GR1+F2GL1
107 WRITE(*,333) F2G1
108 333 FORMAT(' ','c2gg1 = ',D20.10)
109csm end added
110 35 FORMAT(' ','c3qq1 = ',D20.10)
111 R22R=DAIND(0.0D0,1.0D0,F22R,EPS,KEY,MAX,KOUNT,EST)
112 IF(FA.GT.1.1D0) THEN
113 R22S=DAIND(0.0D0,1.0D0,F22S,EPS,KEY,MAX,KOUNT,EST)
114 ELSE
115 R22S=0.0D0
116 ENDIF
117 F2QR2=R22R
118 F2QL2=R22S
119 F2QVS2=COR2(0.0D0)
120 F2Q2=F2QR2+F2QL2+F2QVS2
121csm WRITE(*,*) F2QR2,F2QL2, F2QVS2
122 WRITE(*,37) F2Q2
123 37 FORMAT(' ','c2qq2 = ',D20.10)
124csm beg added
125 RPS22R=DAIND(0.0D0,1.0D0,FPS22R,EPS,KEY,MAX,KOUNT,EST)
126 RPS22S=0.0D0
127 F2PSR2=RPS22R
128 F2PSL2=RPS22S
129 F2PS2=F2PSR2+F2PSL2
130 WRITE(*,337) F2PS2
131 337 FORMAT(' ','c2qqps2 = ',D20.10)
132 RG22R=DAIND(0.0D0,1.0D0,FG22R,EPS,KEY,MAX,KOUNT,EST)
133 RG22S=0.0D0
134 F2GR2=RG22R
135 F2GL2=RG22S
136 F2G2=F2GR2+F2GL2
137 WRITE(*,338) F2G2
138 338 FORMAT(' ','c2gg2 = ',D20.10)
139
140 RL1R=DAIND(0.0D0,1.0D0,FL1R,EPS,KEY,MAX,KOUNT,EST)
141 FLQR1=RL1R
142 FLQL1=0.d0
143 FLQ1=FLQR1+FLQL1
144 WRITE(*,331) FLQ1
145 331 FORMAT(' ','cLqq1 = ',D20.10)
146
147 RL2R=DAIND(0.0D0,1.0D0,FL2R,EPS,KEY,MAX,KOUNT,EST)
148 FLQR2=RL2R
149 FLQL2=0.d0
150 FLQ2=FLQR2+FLQL2
151 WRITE(*,332) FLQ2
152 332 FORMAT(' ','cLqq2 = ',D20.10)
153
154 RPSL2R=DAIND(0.0D0,1.0D0,FPSL2R,EPS,KEY,MAX,KOUNT,EST)
155 RPSL2S=0.0D0
156 FLPSR2=RPSL2R
157 FLPSL2=RPSL2S
158 FLPS2=FLPSR2+FLPSL2
159 WRITE(*,334) FLPS2
160 334 FORMAT(' ','cLqqps2 = ',D20.10)
161
162 RGL1R=DAIND(0.0D0,1.0D0,FLG1R,EPS,KEY,MAX,KOUNT,EST)
163 RGL1S=0.0D0
164 FLGR1=RGL1R
165 FLGL1=RGL1S
166 FLG1=FLGR1+FLGL1
167 WRITE(*,335) FLG1
168 335 FORMAT(' ','cLgg1 = ',D20.10)
169
170 RGL2R=DAIND(0.0D0,1.0D0,FLG2R,EPS,KEY,MAX,KOUNT,EST)
171 RGL2S=0.0D0
172 FLGR2=RGL2R
173 FLGL2=RGL2S
174 FLG2=FLGR2+FLGL2
175 WRITE(*,336) FLG2
176 336 FORMAT(' ','cLgg2 = ',D20.10)
177
178csm end added
179
180 WRITE(*,325) DN
181 R30R=DAIND(0.0D0,1.0D0,F30R,EPS,KEY,MAX,KOUNT,EST)
182 WRITE(*,31) R30R
183 31 FORMAT(' ','c3qq0 = ',D20.10)
184 R31R=DAIND(0.0D0,1.0D0,F31R,EPS,KEY,MAX,KOUNT,EST)
185 IF(FA.GT.1.1D0) THEN
186 R31S=DAIND(0.0D0,1.0D0,F31S,EPS,KEY,MAX,KOUNT,EST)
187 ELSE
188 R31S=0.0D0
189 ENDIF
190 F3QR1=R31R
191 F3QL1=R31S
192 F3QVS1=COR1(0.0D0)
193 F3Q1=F3QR1+F3QL1+F3QVS1
194 WRITE(*,35) F3Q1
195 R32R=DAIND(0.0D0,1.0D0,F32R,EPS,KEY,MAX,KOUNT,EST)
196 IF(FA.GT.1.1D0) THEN
197 R32S=DAIND(0.0D0,1.0D0,F32S,EPS,KEY,MAX,KOUNT,EST)
198 ELSE
199 R32S=0.0D0
200 ENDIF
201 F3QR2=R32R
202 F3QL2=R32S
203 F3QVS2=COR2(0.0D0)
204 F3Q2=F3QR2+F3QL2+F3QVS2
205 WRITE(*,39) F3Q2
206 39 FORMAT(' ','c3qq2 = ',D20.10)
207*
208* csm polarized coefficient functions
209*
210 R1G1=DAIND(0.0D0,1.0D0,F1G1,EPS,KEY,MAX,KOUNT,EST)
211 G1G1=R1G1
212csm WRITE(*,41) G1G1
213 41 FORMAT(' ','THE ORDER ALPHAS CONTRIBUTION TO G1 =',D20.10)
214 R1G2=DAIND(0.0D0,1.0D0,F1G2,EPS,KEY,MAX,KOUNT,EST)
215 G1G2=R1G2
216csm WRITE(*,43) G1G2
217 43 FORMAT(' ','ORDER ALPHAS**2 GLUON CONTRIBUTION TO G1 =',D20.10)
218 R1Q2=DAIND(0.0D0,1.0D0,F1Q2,EPS,KEY,MAX,KOUNT,EST)
219 G1Q2=R1Q2
220csm WRITE(*,44) G1Q2
221 44 FORMAT(' ','ORDER ALPHAS**2 QUARK CONTRIBUTION TO G1 =',D20.10)
222*
223*
224 END
225
226
227C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
228C V + Q -> Q + G (SINGLET AND NONSINGLET) , F2(X,Q^2) , MSBAR SCHEME .
229
230 DOUBLE PRECISION FUNCTION C2NS1(Y)
231 IMPLICIT REAL*8(A-H,O-Z)
232 DOUBLE PRECISION NF
233 COMMON/PAR/NF,Q2,FS,RS,DN
234 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
235 FS2=FS*FS
236 DL1=DLOG(Q2/FS2)
237 DL=DLOG(Y)
238 DLM=DLOG(1.0D0-Y)
239 DM=1.0D0/(1.0D0-Y)
240 A=2.0D0*(1.0D0+Y)*(DL-DLM-DL1)-4.0D0*DM*DL+6.0D0+4.0D0*Y
241 C2NS1=CF*A
242 RETURN
243 END
244
245
246C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
247C V + Q -> Q + G (SINGLET AND NONSINGLET) , F3(X,Q^2) , MSBAR SCHEME .
248
249
250 DOUBLE PRECISION FUNCTION C3NS1(Y)
251 IMPLICIT REAL*8(A-H,O-Z)
252 DOUBLE PRECISION NF
253 COMMON/PAR/NF,Q2,FS,RS,DN
254 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
255 FS2=FS*FS
256 DL1=DLOG(Q2/FS2)
257 DL=DLOG(Y)
258 DLM=DLOG(1.0D0-Y)
259 DM=1.0D0/(1.0D0-Y)
260 A=2.0D0*(1.0D0+Y)*(DL-DLM-DL1)-4.0D0*DM*DL+4.0D0+2.0D0*Y
261 C3NS1=CF*A
262 RETURN
263 END
264
265C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
266C V + Q -> Q + G (SINGLET AND NONSINGLET) , F2(X,Q^2) AND F3(X,Q^2) ,
267C MSBAR SCHEME .THIS IS THE SOFT GLUON CONTRIBUTION.
268
269
270
271 DOUBLE PRECISION FUNCTION SOFT1(Y)
272 IMPLICIT REAL*8(A-H,O-Z)
273 DOUBLE PRECISION NF
274 COMMON/PAR/NF,Q2,FS,RS,DN
275 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
276 FS2=FS*FS
277 DL1=DLOG(Q2/FS2)
278 DLM=DLOG(1.0D0-Y)
279 DM=1.0D0/(1.0D0-Y)
280 A=4.0D0*(DLM+DL1)-3.0D0
281 SOFT1=CF*A*DM
282 RETURN
283 END
284
285C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
286C V + Q -> Q + G (SINGLET AND NONSINGLET) , F2(X,Q^2)
287C MSBAR SCHEME . THIS IS THE SOFT PLUS VIRTUAL GLUON CONTRIBUTION.
288
289
290
291 DOUBLE PRECISION FUNCTION BACOR1(Z)
292 IMPLICIT REAL*8(A-H,O-Z)
293 DOUBLE PRECISION NF
294 COMMON/PAR/NF,Q2,FS,RS,DN
295 COMMON/CHOICE/MODE
296 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
297 ZETA2=1.6449340668D0
298 FS2=FS*FS
299 DL1=DLOG(Q2/FS2)
300 CALL STRUCT(Z,FS,MODE,UV,DV,SS,ST,CH,BOT,GL)
301 SIGMA=UV+DV+4.0D0*SS+2.0D0*(ST+CH)
302 FA=SIGMA
303 DLM=DLOG(1.0D0-Z)
304 C1=2.0D0*DLM**2-3.0D0*DLM-9.0D0-4.0D0*ZETA2
305 C2=4.0D0*DLM+3.0D0
306 BACOR1=4.0D0*(C1+C2*DL1)*FA/3.0D0
307 RETURN
308 END
309
310 DOUBLE PRECISION FUNCTION COR1(Z)
311 IMPLICIT REAL*8(A-H,O-Z)
312 DOUBLE PRECISION NF
313 COMMON/PAR/NF,Q2,FS,RS,DN
314 COMMON/CHOICE/MODE
315 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
316 ZETA2=1.6449340668D0
317 FS2=FS*FS
318 DL1=DLOG(Q2/FS2)
319 DLM=DLOG(1.0D0-Z)
320 C1=2.0D0*DLM**2-3.0D0*DLM-9.0D0-4.0D0*ZETA2
321 C2=4.0D0*DLM+3.0D0
322 COR1=CF*(C1+C2*DL1)
323 RETURN
324 END
325
326C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
327C V + Q -> Q + G (NONSINGLET) , F3(X,Q^2) , MSBAR SCHEME .
328C THIS IS THE SOFT PLUS VIRTUAL GLUON CONTRIBUTION.
329
330
331
332 DOUBLE PRECISION FUNCTION COR3(Z)
333 IMPLICIT REAL*8(A-H,O-Z)
334 DOUBLE PRECISION NF
335 COMMON/PAR/NF,Q2,FS,RS,DN
336 COMMON/CHOICE/MODE
337 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
338 ZETA2=1.6449340668D0
339 FS2=FS*FS
340 DL1=DLOG(Q2/FS2)
341 CALL STRUCT(Z,FS,MODE,UV,DV,SS,ST,CH,BOT,GL)
342 FA=UV+DV
343 DLM=DLOG(1.0D0-Z)
344 C1=2.0D0*DLM**2-3.0D0*DLM-9.0D0-4.0D0*ZETA2
345 C2=4.0D0*DLM+3.0D0
346 COR3=CF*(C1+C2*DL1)*FA
347 RETURN
348 END
349
350C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
351C V + Q -> Q + G (NONSINGLET) , F3(X,Q^2) , DIS SCHEME .
352
353
354 DOUBLE PRECISION FUNCTION D3NS1(Y)
355 IMPLICIT REAL*8(A-H,O-Z)
356 DOUBLE PRECISION NF
357 COMMON/PAR/NF,Q2,FS,RS,DN
358 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
359 FS2=FS*FS
360 DL1=DLOG(Q2/FS2)
361 A=-2.0D0*(1.0D0+Y)*DL1-2.0D0-2.0D0*Y
362 D3NS1=CF*A
363 RETURN
364 END
365
366
367C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
368C V + Q -> Q + G (SINGLET AND NONSINGLET) , F2(X,Q^2) AND F3(X,Q^2) ,
369C DIS SCHEME .THIS IS THE SOFT GLUON CONTRIBUTION.
370
371
372
373 DOUBLE PRECISION FUNCTION DSOFT1(Y)
374 IMPLICIT REAL*8(A-H,O-Z)
375 DOUBLE PRECISION NF
376 COMMON/PAR/NF,Q2,FS,RS,DN
377 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
378 FS2=FS*FS
379 DL1=DLOG(Q2/FS2)
380 DM=1.0D0/(1.0D0-Y)
381 A=4.0D0*DL1
382 DSOFT1=CF*A*DM
383 RETURN
384 END
385
386
387C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
388C V + Q -> Q + G (SINGLET AND NONSINGLET) , F2(X,Q^2)
389C DIS SCHEME . THIS IS THE SOFT PLUS VIRTUAL GLUON CONTRIBUTION.
390
391
392
393 DOUBLE PRECISION FUNCTION DCOR1(Z)
394 IMPLICIT REAL*8(A-H,O-Z)
395 DOUBLE PRECISION NF
396 COMMON/PAR/NF,Q2,FS,RS,DN
397 COMMON/CHOICE/MODE
398 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
399 FS2=FS*FS
400 DL1=DLOG(Q2/FS2)
401 CALL STRUCT(Z,FS,MODE,UV,DV,SS,ST,CH,BOT,GL)
402 SIGMA=UV+DV+4.0D0*SS+2.0D0*(ST+CH)
403 FA=SIGMA
404 DLM=DLOG(1.0D0-Z)
405 A=(4.0D0*DLM+3.0D0)*DL1
406 DCOR1=CF*A*FA
407 RETURN
408 END
409
410
411C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
412C V + Q -> Q + G (NONSINGLET) , F3(X,Q^2) , DIS SCHEME .
413C THIS IS THE SOFT PLUS VIRTUAL GLUON CONTRIBUTION.
414
415
416 DOUBLE PRECISION FUNCTION DCOR3(Z)
417 IMPLICIT REAL*8(A-H,O-Z)
418 DOUBLE PRECISION NF
419 COMMON/PAR/NF,Q2,FS,RS,DN
420 COMMON/CHOICE/MODE
421 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
422 FS2=FS*FS
423 DL1=DLOG(Q2/FS2)
424 CALL STRUCT(Z,FS,MODE,UV,DV,SS,ST,CH,BOT,GL)
425 FA=UV+DV
426 DLM=DLOG(1.0D0-Z)
427 A=(4.0D0*DLM+3.0D0)*DL1
428 DCOR3=CF*A*FA
429 RETURN
430 END
431
432
433C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
434C V + G -> Q + QB , F2(X,Q^2) , MSBAR SCHEME .
435
436
437 DOUBLE PRECISION FUNCTION C2G1(Y)
438 IMPLICIT REAL*8(A-H,O-Z)
439 DOUBLE PRECISION NF
440 COMMON/PAR/NF,Q2,FS,RS,DN
441 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
442 FS2=FS*FS
443 DL1=DLOG(Q2/FS2)
444 DL=DLOG(Y)
445 DLM=DLOG(1.0D0-Y)
446 A1=(2.0D0-4.0D0*Y+4.0D0*Y*Y)*(DLM+DL1-DL)+16.0D0*Y*(1.0D0-Y)
447 1-2.0D0
448 C2G1=NF*A1
449 RETURN
450 END
451
452
453C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESSES :
454C V + Q -> Q + G + G ; V + Q -> Q + Q + QB WITH NONIDENTICAL QUARKS
455C IN THE FINAL STATE (NONSINGLET) , F2(X,Q^2) , MSBAR SCHEME .
456
457 DOUBLE PRECISION FUNCTION C2NSP(Y)
458 IMPLICIT REAL*8(A-H,O-Z)
459 COMPLEX*16 WGPLG
460 DOUBLE PRECISION NF
461 COMMON/PAR/NF,Q2,FS,RS,DN
462 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
463 FS2=FS*FS
464 RS2=RS*RS
465 DL1=DLOG(Q2/FS2)
466 DL2=DL1*DL1
467 EL1=DLOG(FS2/RS2)
468 S121MX=WGPLG(1,2,1.0D0-Y)
469 S211MX=WGPLG(2,1,1.0D0-Y)
470 S21MX=WGPLG(2,1,-Y)
471 S111MX=WGPLG(1,1,1.0D0-Y)
472 S11MX=WGPLG(1,1,-Y)
473 ZETA2=1.6449340668D0
474 ZETA3=1.2020569031D0
475 DM=1.0D0/(1.0D0-Y)
476 DL=DLOG(Y)
477 DLP=DLOG(1.0D0+Y)
478 DLM=DLOG(1.0D0-Y)
479 A1=DL**3*(3.0D0+3.0D0*Y-8.0D0*DM/3.0D0)+
480 1DL**2*DLM*(-16.0D0-16.0D0*Y+24.0D0*DM)+
481 2DL**2*(16.0D0+36.0D0*Y*Y*Y/5.0D0+24*Y*Y+
482 314.0D0*Y-3.0D0*DM)+
483 4DL*DLM**2*(16.0D0+16.0D0*Y-28.0D0*DM)
484 A2=DL*DLM*(-40.0D0-56.0D0*Y+24.0D0*DM)+
485 2DL*DLP*(40.0D0-72.0D0*Y*Y*Y/5.0D0-48.0D0*Y*Y+8.0D0*Y)+
486 38.0D0*DL*DLP/Y/Y/5.0D0-8.0D0*DL/Y/5.0D0+
487 4ZETA2*DL*(-28.0D0-28.0D0*Y+48.0D0*DM)+
488 5DL*S111MX*(-4.0D0-4.0D0*Y)+
489 6DL*S11MX*(16.0D0+64.0D0*Y-48.0D0*DM)
490 A3=DL*(-146.0D0/5.0D0+72.0D0*Y*Y/5.0D0-356.0D0*Y/5.0D0+
491 161.0D0*DM)+
492 2DLM**3*(-4.0D0-4.0D0*Y)+
493 3DLM**2*(14.0D0+18.0D0*Y)+
494 4ZETA2*DLM*(24.0D0-24.0D0*Y)+
495 5DLM*S111MX*(8.0D0+8.0D0*Y-8.0D0*DM)+
496 6DLM*(-32.0D0+84.0D0*Y)+
497 7ZETA2*(-28.0D0-72.0D0*Y*Y*Y/5.0D0-48.0D0*Y*Y-12.0D0*Y)
498 A4=ZETA3*(-32.0D0-32.0D0*Y)+
499 1S111MX*(-8.0D0-24.0D0*Y-12.0D0*DM)+
500 2S211MX*(-16.0D0-16.0D0*Y+24.0D0*DM)+
501 3S121MX*(4.0D0+52.0D0*Y-24.0D0*DM)+
502 4S11MX*(40.0D0-72.0D0*Y*Y*Y/5.0D0-48.0D0*Y*Y+8.0D0*Y)+
503 58.0D0*S11MX/Y/Y/5.0D0+8.0D0/Y/5.0D0+
504 6S21MX*(-32.0D0-128.0D0*Y+96.0D0*DM)+72.0D0*ZETA3*DM+
505 772.0D0*Y*Y/5.0D0+38.0D0/5.0D0-543.0D0*Y/5.0D0
506 A5=DL*(6.0D0+6.0D0*Y-8.0D0*DM)+
507 1DLM*(-8.0D0-8.0D0*Y)-
508 210.0D0-2.0D0*Y
509 A6=DL**2*(-8.0D0-8.0D0*Y+8.0D0*DM)+
510 1DL*DLM*(28.0D0+28.0D0*Y-48.0D0*DM)+
511 2DL*(-8.0D0-24.0D0*Y-12.0D0*DM)+
512 3DLM**2*(-12.0D0-12.0D0*Y)+
513 4DLM*(16.0D0+24.0D0*Y)+
514 5ZETA2*(16.0D0+16.0D0*Y)+
515 6S111MX*(4.0D0+4.0D0*Y)+
516 722.0D0+44.0D0*Y
517 FCF=CF2*(A1+A2+A3+A4+A5*DL2+A6*DL1)
518 B1=DL**3*(1.0D0+Y-2.0D0*DM)+
519 1DL**2*DLM*(-2.0D0-2.0D0*Y+4.0D0*DM)+
520 2DL**2*(43.0D0/6.0D0-18.0D0*Y*Y*Y/5.0D0-12.0D0*Y*Y+
521 367.0D0*Y/6.0D0-55.0D0*DM/3.0D0)+
522 4DL*DLM*(-32.0D0/3.0D0-32.0D0*Y/3.0D0+88.0D0*DM/3.0D0)+
523 5DL*DLP*(-20.0D0+36.0D0*Y*Y*Y/5.0D0+24.0D0*Y*Y-4.0D0*Y)
524 6-4.0D0*DL*DLP/Y/Y/5.0D0+4.0D0*DL/Y/5.0D0
525 B2=DL*S111MX*(-4.0D0-4.0D0*Y+8.0D0*DM)+
526 1DL*S11MX*(-8.0D0-32.0D0*Y+24.0D0*DM)+
527 2DL*(604.0D0/15.0D0-36.0D0*Y*Y/5.0D0+1474.0D0*Y/15.0D0-
528 3239.0D0*DM/3.0D0)+
529 4DLM**2*(11.0D0/3.0D0+11.0D0*Y/3.0D0)+
530 5ZETA2*DLM*(24.0D0*Y)+
531 6DLM*S111MX*(-4.0D0-4.0D0*Y+8.0D0*DM)
532 B3=DLM*(16.0D0/9.0D0-740.0D0*Y/9.0D0)+
533 1ZETA2*(-28.0D0/3.0D0+36.0D0*Y*Y*Y/5.0D0+24.0D0*Y*Y
534 2-52.0D0*Y/3.0D0)+ZETA3*(-2.0D0-2.0D0*Y)+
535 3S111MX*(-10.0D0/3.0D0-10.0D0*Y/3.0D0+44.0D0*DM/3.0D0)+
536 4S211MX*(12.0D0+12.0D0*Y-24.0D0*DM)+
537 5S121MX*(-8.0D0-32.0D0*Y+24.0D0*DM)
538 B4=S11MX*(-20.0D0+36.0D0*Y*Y*Y/5.0D0+24.0D0*Y*Y-4.0D0*Y)
539 1-4.0D0*S11MX/Y/Y/5.0D0-4.0D0/Y/5.0D0+
540 2S21MX*(16.0D0+64.0D0*Y-48.0D0*DM)-36.0D0*ZETA3*DM+
541 31522.0D0/135.0D0-36.0D0*Y*Y/5.0D0+18733.0D0*Y/135.0D0
542 B5=11.0D0/3.0D0+11.0D0*Y/3.0D0
543 B6=DL**2*(-2.0D0-2.0D0*Y+4.0D0*DM)+
544 1DL*(-32.0D0/3.0D0-32.0D0*Y/3.0D0+88.0D0*DM/3.0D0)+
545 2DLM*(22.0D0/3.0D0+22.0D0*Y/3.0D0)+
546 3ZETA2*(4.0D0+4.0D0*Y)-
547 492.0D0/9.0D0-506.0D0*Y/9.0D0
548 B7=(22.0D0/3.0D0+22.0D0*Y/3.0D0)*(DLM-DL)+
549 144.0D0*DL*DM/3.0D0-22.0D0-44.0D0*Y/3.0D0
550 FCA=CACF*(B1+B2+B3+B4+B5*DL2+2.0D0*B5*DL1*EL1+B6*DL1+
551 1B7*EL1)
552 C1=DL**2*(-5.0D0/3.0D0-5.0D0*Y/3.0D0+10.0D0*DM/3.0D0)+
553 1DL*DLM*(8.0D0/3.0D0+8.0D0*Y/3.0D0-16.0D0*DM/3.0D0)+
554 2DL*(-26.0D0/3.0D0-38.0D0*Y/3.0D0+38.0D0*DM/3.0D0)+
555 3DLM**2*(-2.0D0/3.0D0-2.0D0*Y/3.0D0)+
556 4DLM*(32.0D0/9.0D0+68.0D0*Y/9.0D0)+
557 5ZETA2*(4.0D0/3.0D0+4.0D0*Y/3.0D0)+
558 6S111MX*(4.0D0/3.0D0+4.0D0*Y/3.0D0-8.0D0*DM/3.0D0)
559 7-158.0D0/27.0D0-488.0D0*Y/27.0D0
560 C2=-2.0D0/3.0D0-2.0D0*Y/3.0D0
561 C3=DL*(8.0D0/3.0D0+8.0D0*Y/3.0D0-16.0D0*DM/3.0D0)+
562 1DLM*(-4.0D0/3.0D0-4.0D0*Y/3.0D0)+
563 232.0D0/9.0D0+68.0D0*Y/9.0D0
564 C4=(-4.0D0/3.0D0-4.0D0*Y/3.0D0)*(DLM-DL)-
565 18.0D0*DL*DM/3.0D0+4.0D0+8.0D0*Y/3.0D0
566 FNF=2.0D0*NF*CFTF*(C1+C2*DL2+2.0D0*C2*DL1*EL1+C3*DL1+C4*EL1)
567 C2NSP=FCF+FCA+FNF
568 RETURN
569 END
570
571C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESSES :
572C V + Q -> Q + G + G ; V + Q -> Q + Q + QB WITH NONIDENTICAL QUARKS
573C IN THE FINAL STATE (NONSINGLET) , F3(X,Q^2) , MS BAR SCHEME .
574
575 DOUBLE PRECISION FUNCTION C3NSP(Y)
576 IMPLICIT REAL*8(A-H,O-Z)
577 COMPLEX*16 WGPLG
578 DOUBLE PRECISION NF
579 COMMON/PAR/NF,Q2,FS,RS,DN
580 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
581 FS2=FS*FS
582 RS2=RS*RS
583 DL1=DLOG(Q2/FS2)
584 DL2=DL1*DL1
585 EL1=DLOG(FS2/RS2)
586 S121MX=WGPLG(1,2,1.0D0-Y)
587 S211MX=WGPLG(2,1,1.0D0-Y)
588 S21MX=WGPLG(2,1,-Y)
589 S111MX=WGPLG(1,1,1.0D0-Y)
590 S11MX=WGPLG(1,1,-Y)
591 ZETA2=1.6449340668D0
592 ZETA3=1.2020569031D0
593 DM=1.0D0/(1.0D0-Y)
594 DL=DLOG(Y)
595 DLP=DLOG(1.0D0+Y)
596 DLM=DLOG(1.0D0-Y)
597 A1=DL**3*(3.0D0+3.0D0*Y-8.0D0*DM/3.0D0)+
598 1DL**2*DLM*(-16.0D0-16.0D0*Y+24.0D0*DM)+
599 2DL**2*(16.0D0+4*Y*Y+22.0D0*Y-3.0D0*DM)+
600 3DL*DLM**2*(16.0D0+16.0D0*Y-28.0D0*DM)
601 A2=DL*DLM*(-28.0D0-44.0D0*Y+24.0D0*DM)+
602 1DL*DLP*(-8.0D0-8.0D0*Y*Y-8.0D0*Y-8.0D0/Y)+
603 2ZETA2*DL*(-28.0D0-28.0D0*Y+48.0D0*DM)+
604 3DL*S111MX*(-4.0D0-4.0D0*Y)+
605 4DL*S11MX*(32.0D0+16.0D0*Y-48.0D0*DM)
606 A3=DL*(-26.0D0-32.0D0*Y+61.0D0*DM)+
607 1DLM**3*(-4.0D0-4.0D0*Y)+
608 2DLM**2*(10.0D0+14.0D0*Y)+
609 3ZETA2*DLM*(8.0D0+24.0D0*Y)+
610 4DLM*S111MX*(8.0D0+8.0D0*Y-8.0D0*DM)+
611 5DLM*(10.0D0+6.0D0*Y)+
612 6ZETA2*(-28.0D0-8.0D0*Y*Y-28.0D0*Y)
613 A4=ZETA3*(-32.0D0-32.0D0*Y)+
614 1S111MX*(-4.0D0-20.0D0*Y-12.0D0*DM)+
615 2S211MX*(-16.0D0-16.0D0*Y+24.0D0*DM)+
616 3S121MX*(20.0D0+4.0D0*Y-24.0D0*DM)+
617 4S11MX*(-8.0D0-8.0D0*Y*Y-8.0D0*Y-8.0D0/Y)+
618 6S21MX*(-64.0D0-32.0D0*Y+96.0D0*DM)+72.0D0*ZETA3*DM-
619 741.0D0-10.0D0*Y
620 A5=DL*(6.0D0+6.0D0*Y-8.0D0*DM)+
621 1DLM*(-8.0D0-8.0D0*Y)-
622 210.0D0-2.0D0*Y
623 A6=DL**2*(-8.0D0-8.0D0*Y+8.0D0*DM)+
624 1DL*DLM*(28.0D0+28.0D0*Y-48.0D0*DM)+
625 2DL*(-4.0D0-20.0D0*Y-12.0D0*DM)+
626 3DLM**2*(-12.0D0-12.0D0*Y)+
627 4DLM*(8.0D0+16.0D0*Y)+
628 5ZETA2*(16.0D0+16.0D0*Y)+
629 6S111MX*(4.0D0+4.0D0*Y)+
630 716.0D0+38.0D0*Y
631 FCF=CF2*(A1+A2+A3+A4+A5*DL2+A6*DL1)
632 B1=DL**3*(1.0D0+Y-2.0D0*DM)+
633 1DL**2*DLM*(-2.0D0-2.0D0*Y+4.0D0*DM)+
634 2DL**2*(31.0D0/6.0D0-2.0D0*Y*Y+31.0D0*Y/6.0D0-
635 355.0D0*DM/3.0D0)+
636 4DL*DLM*(-32.0D0/3.0D0-32.0D0*Y/3.0D0+88.0D0*DM/3.0D0)+
637 5DL*DLP*(4.0D0+4.0D0*Y*Y+4.0D0*Y+4.0D0/Y)
638 B2=DL*S111MX*(-4.0D0-4.0D0*Y+8.0D0*DM)+
639 1DL*S11MX*(-16.0D0-8.0D0*Y+24.0D0*DM)+
640 2DL*(32.0D0+66.0D0*Y-239.0D0*DM/3.0D0)+
641 3DLM**2*(11.0D0/3.0D0+11.0D0*Y/3.0D0)+
642 4ZETA2*DLM*(8.0D0)+
643 5DLM*S111MX*(-4.0D0-4.0D0*Y+8.0D0*DM)
644 B3=DLM*(-134.0D0/9.0D0-314.0D0*Y/9.0D0)+
645 1ZETA2*(-16.0D0/3.0D0+4.0D0*Y*Y-16.0D0*Y/3.0D0)+
646 2ZETA3*(-2.0D0-2.0D0*Y)+
647 3S111MX*(-10.0D0/3.0D0-10.0D0*Y/3.0D0+44.0D0*DM/3.0D0)+
648 4S211MX*(12.0D0+12.0D0*Y-24.0D0*DM)+
649 5S121MX*(-16.0D0-8.0D0*Y+24.0D0*DM)
650 B4=S11MX*(4.0D0+4.0D0*Y*Y+4.0D0*Y+4.0D0/Y)+
651 1S21MX*(32.0D0+16.0D0*Y-48.0D0*DM)-36.0D0*ZETA3*DM+
652 2725.0D0/27.0D0+1820.0D0*Y/27.0D0
653 B5=11.0D0/3.0D0+11.0D0*Y/3.0D0
654 B6=DL**2*(-2.0D0-2.0D0*Y+4.0D0*DM)+
655 1DL*(-32.0D0/3.0D0-32.0D0*Y/3.0D0+88.0D0*DM/3.0D0)+
656 2DLM*(22.0D0/3.0D0+22.0D0*Y/3.0D0)+
657 3ZETA2*(4.0D0+4.0D0*Y)-
658 426.0D0/9.0D0-440.0D0*Y/9.0D0
659 B7=(22.0D0/3.0D0+22.0D0*Y/3.0D0)*(DLM-DL)+
660 144.0D0*DL*DM/3.0D0-44.0D0/3.0D0-22.0D0*Y/3.0D0
661 FCA=CACF*(B1+B2+B3+B4+B5*DL2+2.0D0*B5*DL1*EL1+B6*DL1+
662 1B7*EL1)
663 C1=DL**2*(-5.0D0/3.0D0-5.0D0*Y/3.0D0+10.0D0*DM/3.0D0)+
664 1DL*DLM*(8.0D0/3.0D0+8.0D0*Y/3.0D0-16.0D0*DM/3.0D0)+
665 2DL*(-6.0D0-10.0D0*Y+38.0D0*DM/3.0D0)+
666 3DLM**2*(-2.0D0/3.0D0-2.0D0*Y/3.0D0)+
667 4DLM*(20.0D0/9.0D0+56.0D0*Y/9.0D0)+
668 5ZETA2*(4.0D0/3.0D0+4.0D0*Y/3.0D0)+
669 6S111MX*(4.0D0/3.0D0+4.0D0*Y/3.0D0-8.0D0*DM/3.0D0)
670 7-116.0D0/27.0D0-302.0D0*Y/27.0D0
671 C2=-2.0D0/3.0D0-2.0D0*Y/3.0D0
672 C3=DL*(8.0D0/3.0D0+8.0D0*Y/3.0D0-16.0D0*DM/3.0D0)+
673 1DLM*(-4.0D0/3.0D0-4.0D0*Y/3.0D0)+
674 220.0D0/9.0D0+56.0D0*Y/9.0D0
675 C4=(-4.0D0/3.0D0-4.0D0*Y/3.0D0)*(DLM-DL)-
676 18.0D0*DL*DM/3.0D0+8.0D0/3.0D0+4.0D0*Y/3.0D0
677 FNF=2.0D0*NF*CFTF*(C1+C2*DL2+2.0D0*C2*DL1*EL1+C3*DL1+C4*EL1)
678 C3NSP=FCF+FCA+FNF
679C C3NSP=FNF
680 RETURN
681 END
682
683C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESSES :
684C V + Q -> Q + G + G ; V + Q -> Q + Q + QB WITH NONIDENTICAL QUARKS
685C IN THE FINAL STATE (NONSINGLET) ,F2(X,Q^2) AND F3(X,Q^2) ,
686C MSBAR SCHEME .THIS IS THE SOFT GLUON CONTRIBUTION.
687
688 DOUBLE PRECISION FUNCTION SOFT2(Y)
689 IMPLICIT REAL*8(A-H,O-Z)
690 DOUBLE PRECISION NF
691 COMMON/PAR/NF,Q2,FS,RS,DN
692 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
693 FS2=FS*FS
694 RS2=RS*RS
695 DL1=DLOG(Q2/FS2)
696 DL2=DL1*DL1
697 EL1=DLOG(FS2/RS2)
698 ZETA2=1.6449340668D0
699 ZETA3=1.2020569031D0
700 DLM=DLOG(1.0D0-Y)
701 DM=1.0D0/(1.0D0-Y)
702 FCF=8.0D0*DLM**3-18.0D0*DLM**2-DLM*(32.0D0*ZETA2
703 1+27.0D0)-8.0D0*ZETA3+36.0D0*ZETA2+25.5D0+
704 2(16.0D0*DLM+12.0D0)*DL2+(24.0D0*DLM**2-12.0D0*DLM-32.0D0*
705 3ZETA2-45.0D0)*DL1
706 FCA=-22.0D0*DLM**2/3.0D0+DLM*(-8.0D0*ZETA2+367.0D0/9.0D0)+
707 140.0D0*ZETA3+44.0D0*ZETA2/3.0D0-3155.0D0/54.0D0-
708 222.0D0*(DL2+2.0D0*DL1*EL1)/3.0D0+(-44.0D0*DLM/3.0D0-8.0D0*
709 3ZETA2+367.0D0/9.0D0)*DL1+(-44.0D0*DLM/3.0D0+11.0D0)*EL1
710 FNF=4.0D0*DLM**2/3.0D0-58.0D0*DLM/9.0D0-8.0D0*ZETA2/3.0D0
711 1+247.0D0/27.0D0+4.0D0*(DL2+2.0D0*DL1*EL1)/3.0D0+
712 2(8.0D0*DLM/3.0D0-58.0D0/9.0D0)*DL1+(8.0D0*DLM/3.0D0-2.0D0)
713 3*EL1
714 SOFT2=DM*(CF2*FCF+CACF*FCA+2.0D0*NF*CFTF*FNF)
715 RETURN
716 END
717
718C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESSES :
719C V + Q -> Q + G + G ; V + Q -> Q + Q + QB WITH NONIDENTICAL QUARKS
720C IN THE FINAL STATE (NONSINGLET) ,F2(X,Q^2) , MS BAR SCHEME .
721C THIS IS THE SOFT PLUS VIRTUAL GLUON CONTRIBUTION.
722
723
724 DOUBLE PRECISION FUNCTION BACOR2(Z)
725 IMPLICIT REAL*8(A-H,O-Z)
726 DOUBLE PRECISION NF
727 COMMON/PAR/NF,Q2,FS,RS,DN
728 COMMON/CHOICE/MODE
729 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
730 ZETA2=1.6449340668D0
731 ZETA3=1.2020569031D0
732 FS2=FS*FS
733 RS2=RS*RS
734 DL1=DLOG(Q2/FS2)
735 DL2=DL1*DL1
736 EL1=DLOG(FS2/RS2)
737 CALL STRUCT(Z,FS,MODE,UV,DV,SS,ST,CH,BOT,GL)
738 SIGMA=UV+DV+4.0D0*SS+2.0D0*(ST+CH)
739 FA=SIGMA
740 DLM=DLOG(1.0D0-Z)
741 A1=2.0D0*DLM**4-6.0D0*DLM**3-(16.0D0*ZETA2+13.5D0)*DLM**2+
742 1(-8.0D0*ZETA3+36.0D0*ZETA2+25.5D0)*DLM+6.0D0*ZETA2**2-
743 278.0D0*ZETA3+69.0D0*ZETA2+331/8.0D0
744 A2=8.0D0*DLM**2+12.0D0*DLM+4.5D0-8.0D0*ZETA2
745 A3=8.0D0*DLM**3-6.0D0*DLM**2+(45.0D0-32.0D0*ZETA2)*DLM-25.5D0
746 1-12.0D0*ZETA2+40.0D0*ZETA3
747 FCF=CF2*(A1+A2*DL2+A3*DL1)
748 B1=-22.0D0*DLM**3/9.0D0+(-4.0D0*ZETA2+367.0D0/18.0D0)*
749 1DLM**2+(40.0D0*ZETA3+44.0D0*ZETA2/3.0D0-3155.0D0/54.0D0)*DLM+
750 271.0D0*ZETA2**2/5.0D0+140.0D0*ZETA3/3.0D0-251.0D0*ZETA2/3.0D0-
751 35464.0D0/72.0D0
752 B2=-22.0D0*DLM/3.0D0-5.5D0
753 B3=-22.0D0*DLM**2/3.0D0+(367.0D0/9.0D0-8.0D0*ZETA2)*DLM+
754 1215.0D0/6.0D0+88.0D0*ZETA2/3.0D0-12.0D0*ZETA3
755 B4=-22.0D0*DLM**2/3.0D0+11.0D0*DLM+33.0D0+44.0D0*ZETA2/3.0D0
756 FCA=CACF*(B1+B2*DL2+2.0D0*B2*DL1*EL1+B3*DL1+B4*EL1)
757 C1=4.0D0*DLM**3/9.0D0-29.0D0*DLM**2/9.0D0+(-8.0D0*ZETA2/3.0D0+
758 1247.0D0/27.0D0)*DLM+4.0D0*ZETA3/3.0D0+38.0D0*ZETA2/3.0D0+
759 2457.0D0/36.0D0
760 C2=4.0D0*DLM/3.0D0+1.0D0
761 C3=4.0D0*DLM**2/3.0D0-58.0D0*DLM/9.0D0-19.0D0/3.0D0-
762 116.0D0*ZETA2/3.0D0
763 C4=4.0D0*DLM**2/3.0D0-2.0D0*DLM-6.0D0-8.0D0*ZETA2/3.0D0
764 FNF=2.0D0*NF*CFTF*(C1+C2*DL2+2.0D0*C2*DL1*EL1+C3*DL1+C4*EL1)
765 BACOR2=(FCF+FCA+FNF)*FA
766 RETURN
767 END
768
769 DOUBLE PRECISION FUNCTION COR2(Z)
770 IMPLICIT REAL*8(A-H,O-Z)
771 DOUBLE PRECISION NF
772 COMMON/PAR/NF,Q2,FS,RS,DN
773 COMMON/CHOICE/MODE
774 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
775 ZETA2=1.6449340668D0
776 ZETA3=1.2020569031D0
777 FS2=FS*FS
778 RS2=RS*RS
779 DL1=DLOG(Q2/FS2)
780 DL2=DL1*DL1
781 EL1=DLOG(FS2/RS2)
782 DLM=DLOG(1.0D0-Z)
783 A1=2.0D0*DLM**4-6.0D0*DLM**3-(16.0D0*ZETA2+13.5D0)*DLM**2+
784 1(-8.0D0*ZETA3+36.0D0*ZETA2+25.5D0)*DLM+6.0D0*ZETA2**2-
785 278.0D0*ZETA3+69.0D0*ZETA2+331/8.0D0
786 A2=8.0D0*DLM**2+12.0D0*DLM+4.5D0-8.0D0*ZETA2
787 A3=8.0D0*DLM**3-6.0D0*DLM**2+(45.0D0-32.0D0*ZETA2)*DLM-25.5D0
788 1-12.0D0*ZETA2+40.0D0*ZETA3
789 FCF=CF2*(A1+A2*DL2+A3*DL1)
790 B1=-22.0D0*DLM**3/9.0D0+(-4.0D0*ZETA2+367.0D0/18.0D0)*
791 1DLM**2+(40.0D0*ZETA3+44.0D0*ZETA2/3.0D0-3155.0D0/54.0D0)*DLM+
792 271.0D0*ZETA2**2/5.0D0+140.0D0*ZETA3/3.0D0-251.0D0*ZETA2/3.0D0-
793 35465.0D0/72.0D0
794 B2=-22.0D0*DLM/3.0D0-5.5D0
795 B3=-22.0D0*DLM**2/3.0D0+(367.0D0/9.0D0-8.0D0*ZETA2)*DLM+
796 1215.0D0/6.0D0+88.0D0*ZETA2/3.0D0-12.0D0*ZETA3
797 B4=-22.0D0*DLM**2/3.0D0+11.0D0*DLM+33.0D0+44.0D0*ZETA2/3.0D0
798 FCA=CACF*(B1+B2*DL2+2.0D0*B2*DL1*EL1+B3*DL1+B4*EL1)
799 C1=4.0D0*DLM**3/9.0D0-29.0D0*DLM**2/9.0D0+(-8.0D0*ZETA2/3.0D0+
800 1247.0D0/27.0D0)*DLM+4.0D0*ZETA3/3.0D0+38.0D0*ZETA2/3.0D0+
801 2457.0D0/36.0D0
802 C2=4.0D0*DLM/3.0D0+1.0D0
803 C3=4.0D0*DLM**2/3.0D0-58.0D0*DLM/9.0D0-19.0D0/3.0D0-
804 116.0D0*ZETA2/3.0D0
805 C4=4.0D0*DLM**2/3.0D0-2.0D0*DLM-6.0D0-8.0D0*ZETA2/3.0D0
806 FNF=2.0D0*NF*CFTF*(C1+C2*DL2+2.0D0*C2*DL1*EL1+C3*DL1+C4*EL1)
807 COR2=FCF+FCA+FNF
808C COR2=FNF
809 RETURN
810 END
811
812
813C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESSES :
814C V + Q -> Q + G + G ; V + Q -> Q + Q + QB WITH NONIDENTICAL QUARKS
815C IN THE FINAL STATE (NONSINGLET) ,F3(X,Q^2) , MS BAR SCHEME .
816C THIS IS THE SOFT PLUS VIRTUAL GLUON CONTRIBUTION.
817
818
819 DOUBLE PRECISION FUNCTION COR4(Z)
820 IMPLICIT REAL*8(A-H,O-Z)
821 DOUBLE PRECISION NF
822 COMMON/PAR/NF,Q2,FS,RS,DN
823 COMMON/CHOICE/MODE
824 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
825 ZETA2=1.6449340668D0
826 ZETA3=1.2020569031D0
827 FS2=FS*FS
828 RS2=RS*RS
829 DL1=DLOG(Q2/FS2)
830 DL2=DL1*DL1
831 EL1=DLOG(FS2/RS2)
832 DLM=DLOG(1.0D0-Z)
833 A1=2.0D0*DLM**4-6.0D0*DLM**3-(16.0D0*ZETA2+13.5D0)*DLM**2+
834 1(-8.0D0*ZETA3+36.0D0*ZETA2+25.5D0)*DLM+6.0D0*ZETA2**2-
835 278.0D0*ZETA3+69.0D0*ZETA2+331/8.0D0
836 A2=8.0D0*DLM**2+12.0D0*DLM+4.5D0-8.0D0*ZETA2
837 A3=8.0D0*DLM**3-6.0D0*DLM**2-(45.0D0+32.0D0*ZETA2)*DLM-25.5D0
838 1-12.0D0*ZETA2+40.0D0*ZETA3
839 FCF=CF2*(A1+A2*DL2+A3*DL1)
840 B1=-22.0D0*DLM**3/9.0D0+(-4.0D0*ZETA2+367.0D0/18.0D0)*
841 1DLM**2+(40.0D0*ZETA3+44.0D0*ZETA2/3.0D0-3155.0D0/54.0D0)*DLM+
842 271.0D0*ZETA2**2/5.0D0+140.0D0*ZETA3/3.0D0-251.0D0*ZETA2/3.0D0-
843 35465.0D0/72.0D0
844 B2=-22.0D0*DLM/3.0D0-5.5D0
845 B3=-22.0D0*DLM**2/3.0D0+(367.0D0/9.0D0-8.0D0*ZETA2)*DLM+
846 1215.0D0/6.0D0+88.0D0*ZETA2/3.0D0-12.0D0*ZETA3
847 B4=-22.0D0*DLM**2/3.0D0+11.0D0*DLM+33.0D0+44.0D0*ZETA2/3.0D0
848 FCA=CACF*(B1+B2*DL2+2.0D0*B2*DL1*EL1+B3*DL1+B4*EL1)
849 C1=4.0D0*DLM**3/9.0D0-29.0D0*DLM**2/9.0D0+(-8.0D0*ZETA2/3.0D0+
850 1247.0D0/27.0D0)*DLM+4.0D0*ZETA3/3.0D0+38.0D0*ZETA2/3.0D0+
851 2457.0D0/36.0D0
852 C2=4.0D0*DLM/3.0D0+1.0D0
853 C3=4.0D0*DLM**2/3.0D0-58.0D0*DLM/9.0D0-19.0D0/3.0D0-
854 116.0D0*ZETA2/3.0D0
855 C4=4.0D0*DLM**2/3.0D0-2.0D0*DLM-6.0D0-8.0D0*ZETA2/3.0D0
856 FNF=2.0D0*NF*CFTF*(C1+C2*DL2+2.0D0*C2*DL1*EL1+C3*DL1+C4*EL1)
857 COR4=(FCF+FCA+FNF)
858 RETURN
859 END
860
861C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESSES :
862C V + Q -> Q + G + G ; V + Q -> Q + Q + QB WITH NONIDENTICAL QUARKS
863C IN THE FINAL STATE (NONSINGLET) , F3(X,Q^2) , DIS SCHEME .
864
865 DOUBLE PRECISION FUNCTION D3NSP(Y)
866 IMPLICIT REAL*8(A-H,O-Z)
867 COMPLEX*16 WGPLG
868 DOUBLE PRECISION NF
869 COMMON/PAR/NF,Q2,FS,RS,DN
870 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
871 FS2=FS*FS
872 RS2=RS*RS
873 DL1=DLOG(Q2/FS2)
874 DL2=DL1*DL1
875 EL1=DLOG(FS2/RS2)
876 S121MX=WGPLG(1,2,1.0D0-Y)
877 S21MX=WGPLG(2,1,-Y)
878 S111MX=WGPLG(1,1,1.0D0-Y)
879 S11MX=WGPLG(1,1,-Y)
880 ZETA2=1.6449340668D0
881 DM=1.0D0/(1.0D0-Y)
882 DL=DLOG(Y)
883 DLP=DLOG(1.0D0+Y)
884 DLM=DLOG(1.0D0-Y)
885 A1=DL**2*(2-36.0D0*Y*Y*Y/5.0D0-20*Y*Y+10.0D0*Y)+
886 1DL*DLM*(8.0D0+8.0D0*Y)+
887 2DL*DLP*(-48.0D0+72.0D0*Y*Y*Y/5.0D0+40.0D0*Y*Y-16.0D0*Y-
888 38.0D0/Y)-8.0D0*DL*DLP/Y/Y/5.0D0+8.0D0*DL/Y/5.0D0+
889 4DL*S11MX*(16.0D0-48.0D0*Y)+
890 5DL*(-34.0D0/5.0D0-72.0D0*Y*Y/5.0D0+186.0D0*Y/5.0D0)+
891 6ZETA2*DLM*(-16.0D0+48.0D0*Y)
892 A2=DLM*(36.0D0-84.0D0*Y)+
893 1ZETA2*(-8.0D0+72.0D0*Y*Y*Y/5.0D0+40.0D0*Y*Y-24.0D0*Y)+
894 2S111MX*(8.0D0+8.0D0*Y)+
895 3S121MX*(16.0D0-48.0D0*Y)+
896 4S11MX*(-48.0D0+72.0D0*Y*Y*Y/5.0D0+40.0D0*Y*Y-16.0D0*Y-
897 58.0D0/Y)-8.0D0*S11MX/Y/Y/5.0D0-8.0D0/Y/5.0D0+
898 6S21MX*(-32.0D0+96.0D0*Y)-
899 7283.0D0/5.0D0-72.0D0*Y*Y/5.0D0+353.0D0*Y/5.0D0
900 A3=DL*(6.0D0+6.0D0*Y-8.0D0*DM)+
901 1DLM*(-8.0D0-8.0D0*Y)-
902 210.0D0-2.0D0*Y
903 A4=DL**2*(-2.0D0-2.0D0*Y)+
904 1DL*DLM*(8.0D0+8.0D0*Y-16.0D0*DM)+
905 2DL*(4.0D0-4.0D0*Y-12.0D0*DM)+
906 3DLM*(-8.0D0-8.0D0*Y)-
907 426.0D0+14.0D0*Y
908 FCF=CF2*(A1+A2+A3*DL2+A4*DL1)
909 B1=DL**2*(-2.0D0+18.0D0*Y*Y*Y/5.0D0+10.0D0*Y*Y-6.0D0*Y)+
910 1DL*DLP*(24.0D0-36.0D0*Y*Y*Y/5.0D0-20.0D0*Y*Y+8.0D0*Y+
911 24.0D0/Y)+4.0D0*DL*DLP/Y/Y/5.0D0-4.0D0*DL/Y/5.0D0+
912 3DL*S11MX*(-8.0D0+24.0D0*Y)+
913 4DL*(-124.0D0/15.0D0+36.0D0*Y*Y/5.0D0-484.0D0*Y/15.0D0)+
914 5ZETA2*DLM*(8.0D0-24.0D0*Y)
915 B2=DLM*(-50.0D0/3.0D0+142.0D0*Y/3.0D0)+
916 1ZETA2*(4.0D0-36.0D0*Y*Y*Y/5.0D0-20.0D0*Y*Y+12.0D0*Y)+
917 2S121MX*(-8.0D0+24.0D0*Y)+
918 3S11MX*(24.0D0-36.0D0*Y*Y*Y/5.0D0-20.0D0*Y*Y+8.0D0*Y+
919 44.0D0/Y)+4.0D0*S11MX/Y/Y/5.0D0+4.0D0/Y/5.0D0+
920 5S21MX*(16.0D0-48.0D0*Y)+
921 6701.0D0/45.0D0+36.0D0*Y*Y/5.0D0-3211.0D0*Y/45.0D0
922 B3=11.0D0/3.0D0+11.0D0*Y/3.0D0
923 B4=DL**2*(-2.0D0-2.0D0*Y+4.0D0*DM)+
924 1DL*(-32.0D0/3.0D0-32.0D0*Y/3.0D0+88.0D0*DM/3.0D0)+
925 2DLM*(22.0D0/3.0D0+22.0D0*Y/3.0D0)+
926 3ZETA2*(4.0D0+4.0D0*Y)-
927 426.0D0/9.0D0-440.0D0*Y/9.0D0
928 B5=22.0D0/3.0D0+22.0D0*Y/3.0D0
929 FCA=CACF*(B1+B2+B3*DL2+2.0D0*B3*DL1*EL1+B4*DL1+
930 1B5*EL1)
931 C1=DL*(8.0D0/3.0D0+8.0D0*Y/3.0D0)+
932 1DLM*(-4.0D0/3.0D0-4.0D0*Y/3.0D0)+
933 214.0D0/9.0D0+62.0D0*Y/9.0D0
934 C2=-2.0D0/3.0D0-2.0D0*Y/3.0D0
935 C3=DL*(8.0D0/3.0D0+8.0D0*Y/3.0D0-16.0D0*DM/3.0D0)+
936 1DLM*(-4.0D0/3.0D0-4.0D0*Y/3.0D0)+
937 220.0D0/9.0D0+56.0D0*Y/9.0D0
938 C4=-4.0D0/3.0D0-4.0D0*Y/3.0
939 FNF=2.0D0*NF*CFTF*(C1+C2*DL2+2.0D0*C2*DL1*EL1+C3*DL1+C4*EL1)
940 D3NSP=FCF+FCA+FNF
941 RETURN
942 END
943
944C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESSES :
945C V + Q -> Q + G + G ; V + Q -> Q + Q + QB WITH NONIDENTICAL QUARKS
946C IN THE FINAL STATE (NONSINGLET) ,F2(X,Q^2) AND F3(X,Q^2) ,
947C DIS SCHEME .THIS IS THE SOFT GLUON CONTRIBUTION.
948
949 DOUBLE PRECISION FUNCTION DSOFT2(Y)
950 IMPLICIT REAL*8(A-H,O-Z)
951 DOUBLE PRECISION NF
952 COMMON/PAR/NF,Q2,FS,RS,DN
953 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
954 FS2=FS*FS
955 RS2=RS*RS
956 DL1=DLOG(Q2/FS2)
957 DL2=DL1*DL1
958 EL1=DLOG(FS2/RS2)
959 DLM=DLOG(1.0D0-Y)
960 DM=1.0D0/(1.0D0-Y)
961 FCF=(16.0D0*DLM+12.0D0)*DL2
962 FCA=22.0D0*(DL2+2.0D0*DL1*EL1)/3.0D0
963 FNF=4.0D0*(DL2+2.0D0*DL1*EL1)/3.0D0
964 DSOFT2=DM*(CF2*FCF+CACF*FCA+2.0D0*NF*CFTF*FNF)
965 RETURN
966 END
967
968C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESSES :
969C V + Q -> Q + G + G ; V + Q -> Q + Q + QB WITH NONIDENTICAL QUARKS
970C IN THE FINAL STATE (NONSINGLET) , F2(X,Q^2) , DIS SCHEME .
971C THIS IS THE SOFT PLUS VIRTUAL GLUON CONTRIBUTION.
972
973
974 DOUBLE PRECISION FUNCTION DCOR2(Z)
975 IMPLICIT REAL*8(A-H,O-Z)
976 DOUBLE PRECISION NF
977 COMMON/PAR/NF,Q2,FS,RS,DN
978 COMMON/CHOICE/MODE
979 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
980 ZETA2=1.6449340668D0
981 FS2=FS*FS
982 RS2=RS*RS
983 DL1=DLOG(Q2/FS2)
984 DL2=DL1*DL1
985 EL1=DLOG(FS2/RS2)
986 CALL STRUCT(Z,FS,MODE,UV,DV,SS,ST,CH,BOT,GL)
987 SIGMA=UV+DV+4.0D0*SS+2.0D0*(ST+CH)
988 FA=SIGMA
989 DLM=DLOG(1.0D0-Z)
990 A1=8.0D0*DLM**2+12.0D0*DLM+4.5D0-8.0D0*ZETA2
991 FCF=CF2*A1*DL2
992 B1=-22.0D0*DLM/3.0D0-5.5D0
993 FCA=CACF*(B1*DL2+2.0D0*B1*DL1*EL1)
994 C1=4.0D0*DLM/3.0D0+1.0D0
995 FNF=2.0D0*NF*CFTF*(C1*DL2+2.0D0*C1*DL1*EL1)
996 DCOR2=(FCF+FCA+FNF)*FA
997 RETURN
998 END
999
1000C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESSES :
1001C V + Q -> Q + G + G ; V + Q -> Q + Q + QB WITH NONIDENTICAL QUARKS
1002C IN THE FINAL STATE (NONSINGLET) , F3(X,Q^2) , DIS SCHEME .
1003C THIS IS THE SOFT PLUS VIRTUAL GLUON CONTRIBUTION.
1004
1005
1006 DOUBLE PRECISION FUNCTION DCOR4(Z)
1007 IMPLICIT REAL*8(A-H,O-Z)
1008 DOUBLE PRECISION NF
1009 COMMON/PAR/NF,Q2,FS,RS,DN
1010 COMMON/CHOICE/MODE
1011 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1012 ZETA2=1.6449340668D0
1013 FS2=FS*FS
1014 RS2=RS*RS
1015 DL1=DLOG(Q2/FS2)
1016 DL2=DL1*DL1
1017 EL1=DLOG(FS2/RS2)
1018 CALL STRUCT(Z,FS,MODE,UV,DV,SS,ST,CH,BOT,GL)
1019 FA=UV+DV
1020 DLM=DLOG(1.0D0-Z)
1021 A1=8.0D0*DLM**2+12.0D0*DLM+4.5D0-8.0D0*ZETA2
1022 FCF=CF2*A1*DL2
1023 B1=-22.0D0*DLM/3.0D0-5.5D0
1024 FCA=CACF*(B1*DL2+2.0D0*B1*DL1*EL1)
1025 C1=4.0D0*DLM/3.0D0+1.0D0
1026 FNF=2.0D0*NF*CFTF*(C1*DL2+2.0D0*C1*DL1*EL1)
1027 DCOR4=(FCF+FCA+FNF)*FA
1028 RETURN
1029 END
1030
1031
1032
1033
1034C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS:
1035C V + Q -> Q + Q + QB WITH IDENTICAL QUARKS IN THE FINAL STATE
1036C (NONSINGLET) , F2(X,Q^2) , MSBAR SCHEME .
1037
1038
1039 DOUBLE PRECISION FUNCTION C2NSM(Y)
1040 IMPLICIT REAL*8(A-H,O-Z)
1041 COMPLEX*16 WGPLG
1042 DOUBLE PRECISION NF
1043 COMMON/PAR/NF,Q2,FS,RS,DN
1044 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1045 FS2=FS*FS
1046 DL1=DLOG(Q2/FS2)
1047 S121MX=WGPLG(1,2,1.0D0-Y)
1048 S12MX=WGPLG(1,2,-Y)
1049 S211MX=WGPLG(2,1,1.0D0-Y)
1050 S21MX=WGPLG(2,1,-Y)
1051 Z=(1.0D0-Y)/(1.0D0+Y)
1052 S21Z=WGPLG(2,1,Z)
1053 S21MZ=WGPLG(2,1,-Z)
1054 S111MX=WGPLG(1,1,1.0D0-Y)
1055 S11MX=WGPLG(1,1,-Y)
1056 ZETA2=1.6449340668D0
1057 ZETA3=1.2020569031D0
1058 DP=1.0D0/(1.0D0+Y)
1059 DL=DLOG(Y)
1060 DLP=DLOG(1.0D0+Y)
1061 DLM=DLOG(1.0D0-Y)
1062 B1=DL**3*(2.0D0-2.0D0*Y-4.0D0*DP)+
1063 1DL**2*DLM*(-4.0D0+4.0D0*Y+8.0D0*DP)+
1064 2DL**2*DLP*(-16.0D0+40.0D0*Y+40.0D0*DP)+
1065 3DL**2*(-4.0D0+36.0D0*Y*Y*Y/5.0D0-24.0D0*Y*Y-16.0D0*Y)+
1066 4DL*DLM*DLP*(16.0D0-16.0D0*Y-32.0D0*DP)
1067 B2=DL*DLM*(8.0D0+8.0D0*Y)+
1068 1DL*DLP**2*(-48.0D0*Y-16.0D0*DP)+
1069 2DL*DLP*(32.0D0-72.0D0*Y*Y*Y/5.0D0+48.0D0*Y*Y+32.0D0*Y)+
1070 38.0D0*DL*DLP/Y/Y/5.0D0-8.0D0*DL/Y/5.0D0+
1071 4DL*S111MX*(-8.0D0+8.0D0*Y+16.0D0*DP)+
1072 5DL*S11MX*(-16.0D0+16.0D0*Y+32.0D0*DP)+
1073 6DL*(14.0D0/5.0D0+72.0D0*Y*Y/5.0D0-146.0D0*Y/5.0D0
1074 7-16.0D0*DP)+
1075 8ZETA2*DLM*(8.0D0-8.0D0*Y-16.0D0*DP)
1076 B3=DLM*S11MX*(16.0D0-16.0D0*Y-32.0D0*DP)+
1077 1DLM*(16.0D0-16.0D0*Y)+
1078 2ZETA2*DLP*(-48.0D0*Y-16.0D0*DP)+
1079 3DLP*S11MX*(-96.0D0*Y-32.0D0*DP)+
1080 4S21MZ*(16.0D0-16.0D0*Y-32.0D0*DP)+
1081 5S21Z*(-16.0D0+16.0D0*Y+32.0D0*DP)+
1082 6ZETA2*(-4.0D0-72.0D0*Y*Y*Y/5.0D0+48.0D0*Y*Y+20.0D0*Y)+
1083 7ZETA3*(48.0D0*Y+16.0D0*DP)
1084 B4=S111MX*(8.0D0+8.0D0*Y)+
1085 1S211MX*(16.0D0-16.0D0*Y-32.0D0*DP)+
1086 2S121MX*(-8.0D0+8.0D0*Y+16.0D0*DP)+
1087 3S11MX*(32.0D0-72.0D0*Y*Y*Y/5.0D0+48.0D0*Y*Y+32.0D0*Y)+
1088 48.0D0*S11MX/Y/Y/5.0D0+8.0D0/Y/5.0D0+
1089 5S21MX*(48.0D0*Y+16.0D0*DP)+
1090 6S12MX*(-96.0D0*Y-32.0D0*DP)+
1091 772.0D0*Y*Y/5.0D0-162.0D0/5.0D0+82.0D0*Y/5.0D0
1092 B5=DL**2*(-4.0D0+4.0D0*Y+8.0D0*DP)+
1093 1DL*DLP*(16.0D0-16.0D0*Y-32.0D0*DP)+
1094 2DL*(8.0D0+8.0D0*Y)+
1095 3ZETA2*(8.0D0-8.0D0*Y-16.0D0*DP)+
1096 4S11MX*(16.0D0-16.0D0*Y-32.0D0*DP)+
1097 516.0D0-16.0D0*Y
1098 COLOR=CF2-0.5D0*CACF
1099 C2NSM=COLOR*(B1+B2+B3+B4+B5*DL1)
1100 RETURN
1101 END
1102
1103C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS:
1104C V + Q -> Q + Q + QB WITH IDENTICAL QUARKS IN THE FINAL STATE
1105C (NONSINGLET) , F3(X,Q^2) , MSBAR SCHEME .
1106
1107
1108 DOUBLE PRECISION FUNCTION C3NSM(Y)
1109 IMPLICIT REAL*8(A-H,O-Z)
1110 COMPLEX*16 WGPLG
1111 DOUBLE PRECISION NF
1112 COMMON/PAR/NF,Q2,FS,RS,DN
1113 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1114 FS2=FS*FS
1115 DL1=DLOG(Q2/FS2)
1116 S121MX=WGPLG(1,2,1.0D0-Y)
1117 S12MX=WGPLG(1,2,-Y)
1118 S211MX=WGPLG(2,1,1.0D0-Y)
1119 S21MX=WGPLG(2,1,-Y)
1120 Z=(1.0D0-Y)/(1.0D0+Y)
1121 S21Z=WGPLG(2,1,Z)
1122 S21MZ=WGPLG(2,1,-Z)
1123 S111MX=WGPLG(1,1,1.0D0-Y)
1124 S11MX=WGPLG(1,1,-Y)
1125 ZETA2=1.6449340668D0
1126 ZETA3=1.2020569031D0
1127 DP=1.0D0/(1.0D0+Y)
1128 DL=DLOG(Y)
1129 DLP=DLOG(1.0D0+Y)
1130 DLM=DLOG(1.0D0-Y)
1131 B1=DL**3*(2.0D0-2.0D0*Y-4.0D0*DP)+
1132 1DL**2*DLM*(-4.0D0+4.0D0*Y+8.0D0*DP)+
1133 2DL**2*DLP*(-24.0D0+16.0D0*Y+40.0D0*DP)+
1134 3DL**2*(-8.0D0-4.0D0*Y*Y-4.0D0*Y)+
1135 4DL*DLM*DLP*(16.0D0-16.0D0*Y-32.0D0*DP)
1136 B2=DL*DLM*(8.0D0+8.0D0*Y)+
1137 1DL*DLP**2*(16.0D0-16.0D0*DP)+
1138 2DL*DLP*(8.0D0*Y*Y+8.0D0/Y)+
1139 3DL*S111MX*(-8.0D0+8.0D0*Y+16.0D0*DP)+
1140 4DL*S11MX*(-16.0D0+16.0D0*Y+32.0D0*DP)+
1141 5DL*(14.0D0+30.0D0*Y-16.0D0*DP)+
1142 6ZETA2*DLM*(8.0D0-8.0D0*Y-16.0D0*DP)
1143 B3=DLM*S11MX*(16.0D0-16.0D0*Y-32.0D0*DP)+
1144 1DLM*(16.0D0-16.0D0*Y)+
1145 2ZETA2*DLP*(16.0D0-16.0D0*DP)+
1146 3DLP*S11MX*(32.0D0-32.0D0*DP)+
1147 4S21MZ*(16.0D0-16.0D0*Y-32.0D0*DP)+
1148 5S21Z*(-16.0D0+16.0D0*Y+32.0D0*DP)+
1149 6ZETA2*(4.0D0+8.0D0*Y*Y-4.0D0*Y)+
1150 7ZETA3*(-16.0D0+16.0D0*DP)
1151 B4=S111MX*(8.0D0+8.0D0*Y)+
1152 1S211MX*(16.0D0-16.0D0*Y-32.0D0*DP)+
1153 2S121MX*(-8.0D0+8.0D0*Y+16.0D0*DP)+
1154 3S11MX*(8.0D0*Y*Y+8.0D0/Y)+
1155 4S21MX*(-16.0D0+16.0D0*DP)+
1156 5S12MX*(32.0D0-32.0D0*DP)+
1157 630.0D0-30.0D0*Y
1158 B5=DL**2*(-4.0D0+4.0D0*Y+8.0D0*DP)+
1159 1DL*DLP*(16.0D0-16.0D0*Y-32.0D0*DP)+
1160 2DL*(8.0D0+8.0D0*Y)+
1161 3ZETA2*(8.0D0-8.0D0*Y-16.0D0*DP)+
1162 4S11MX*(16.0D0-16.0D0*Y-32.0D0*DP)+
1163 516.0D0-16.0D0*Y
1164 C3NSM=(CF2-0.5D0*CACF)*(B1+B2+B3+B4+B5*DL1)
1165 RETURN
1166 END
1167
1168C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS:
1169C V + Q -> Q + Q + QB WITH IDENTICAL QUARKS IN THE FINAL STATE
1170C (NONSINGLET) , F3(X,Q^2) , DIS SCHEME .
1171
1172
1173 DOUBLE PRECISION FUNCTION D3NSM(Y)
1174 IMPLICIT REAL*8(A-H,O-Z)
1175 COMPLEX*16 WGPLG
1176 DOUBLE PRECISION NF
1177 COMMON/PAR/NF,Q2,FS,RS,DN
1178 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1179 FS2=FS*FS
1180 DL1=DLOG(Q2/FS2)
1181 S12MX=WGPLG(1,2,-Y)
1182 S21MX=WGPLG(2,1,-Y)
1183 S11MX=WGPLG(1,1,-Y)
1184 ZETA2=1.6449340668D0
1185 ZETA3=1.2020569031D0
1186 DP=1.0D0/(1.0D0+Y)
1187 DL=DLOG(Y)
1188 DLP=DLOG(1.0D0+Y)
1189 B1=DL**2*DLP*(-8.0D0-24.0D0*Y)+
1190 1DL**2*(-4.0D0-36.0D0*Y*Y*Y/5.0D0+20.0D0*Y*Y+12.0D0*Y)+
1191 2DL*DLP**2*(16.0D0+48.0D0*Y)+
1192 3DL*DLP*(-32.0D0+72.0D0*Y*Y*Y/5.0D0-40.0D0*Y*Y-32.0D0*Y+
1193 48.0D0/Y)-8.0D0*DL*DLP/Y/Y/5.0D0+8.0D0*DL/Y/5.0D0+
1194 5DL*(56.0D0/5.0D0-72.0D0*Y*Y/5.0D0+296.0D0*Y/5.0D0)+
1195 6ZETA2*DLP*(16.0D0+48.0D0*Y)
1196 B2=DLP*S11MX*(32.0D0+96.0D0*Y)+
1197 1ZETA2*(8.0D0+72.0D0*Y*Y*Y/5.0D0-40.0D0*Y*Y-24.0D0*Y)+
1198 2ZETA3*(-16.0D0-48.0D0*Y)+
1199 3S11MX*(-32.0D0+72.0D0*Y*Y*Y/5.0D0-40.0D0*Y*Y-32.0D0*Y+
1200 48.0D0/Y)-8.0D0*S11MX/Y/Y/5.0D0-8.0D0/Y/5.0D0+
1201 5S21MX*(-16.0D0-48.0D0*Y)+
1202 6S12MX*(32.0D0+96.0D0*Y)+
1203 7312.0D0/5.0D0-72.0D0*Y*Y/5.0D0-232.0D0*Y/5.0D0
1204 B3=DL**2*(-4.0D0+4.0D0*Y+8.0D0*DP)+
1205 1DL*DLP*(16.0D0-16.0D0*Y-32.0D0*DP)+
1206 2DL*(8.0D0+8.0D0*Y)+
1207 3ZETA2*(8.0D0-8.0D0*Y-16.0D0*DP)+
1208 4S11MX*(16.0D0-16.0D0*Y-32.0D0*DP)+
1209 516.0D0-16.0D0*Y
1210 D3NSM=(CF2-0.5D0*CACF)*(B1+B2+B3*DL1)
1211 RETURN
1212 END
1213
1214
1215C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
1216C V + Q -> Q + QB + Q (SINGLET) , F2(X,Q^2) , MSBAR SCHEME .
1217
1218
1219 DOUBLE PRECISION FUNCTION C2S(Y)
1220 IMPLICIT REAL*8(A-H,O-Z)
1221 COMPLEX*16 WGPLG
1222 DOUBLE PRECISION NF
1223 COMMON/PAR/NF,Q2,FS,RS,DN
1224 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1225 FS2=FS*FS
1226 DL1=DLOG(Q2/FS2)
1227 DL2=DL1*DL1
1228 S211MX=WGPLG(2,1,1.0D0-Y)
1229 S111MX=WGPLG(1,1,1.0D0-Y)
1230 S11MX=WGPLG(1,1,-Y)
1231 ZETA2=1.6449340668D0
1232 DL=DLOG(Y)
1233 DLP=DLOG(1.0D0+Y)
1234 DLM=DLOG(1.0D0-Y)
1235 A1=DL**3*(10.0D0/3.0D0+10.0D0*Y/3.0D0)+
1236 1DL**2*DLM*(-8.0D0-8.0D0*Y)+
1237 2DL**2*(-1.0D0-32.0D0*Y*Y/3.0D0+15.0D0*Y)+
1238 3DL*DLM**2*(4.0D0+4.0D0*Y)+
1239 4DL*DLM*(16.0D0*Y*Y)+
1240 5DL*DLP*(-16.0D0-16.0D0*Y*Y/3.0D0-16.0D0*Y-16.0D0/Y/3.0D0)
1241 A2=DL*ZETA2*(-8.0D0-8.0D0*Y)+
1242 1DL*S111MX*(-8.0D0-8.0D0*Y)+
1243 2DL*(56.0D0-128.0D0*Y*Y/9.0D0-88.0D0*Y/3.0D0)+
1244 3DLM**2*(2.0D0-8.0D0*Y*Y/3.0D0-2.0D0*Y+8.0D0/Y/3.0D0)+
1245 5DLM*S111MX*(8.0D0+8.0D0*Y)+
1246 6DLM*(-104.0D0/3.0D0-32.0D0*Y*Y/9.0D0+80.0D0*Y/3.0D0
1247 7+104.0D0/Y/9.0D0)
1248 A3=ZETA2*(-4.0D0+16.0D0*Y*Y/3.0D0-12.0D0*Y-32.0D0/Y/3.0D0)+
1249 1S111MX*(4.0D0+32.0D0*Y*Y/3.0D0-4.0D0*Y+16.0D0/Y/3.0D0)+
1250 2S211MX*(-8.0D0-8.0D0*Y)+
1251 3S11MX*(-16.0D0-16.0D0*Y*Y/3.0D0-16.0D0*Y-16.0D0/Y/3.0D0)+
1252 4158.0D0/9.0D0+448.0D0*Y*Y/27.0D0-422.0D0*Y/9.0D0+
1253 5344.0D0/Y/27.0D0
1254C A3=A3+32.0D0*ZETA2/Y/3.0D0-344.0D0/Y/27.0D0
1255 A4=DL*(4.0D0+4.0D0*Y)+
1256 12.0D0-2.0D0*Y-8.0D0*Y*Y/3.0D0+8.0D0/Y/3.0D0
1257 A5=DL**2*(-8.0D0-8.0D0*Y)+
1258 1DL*DLM*(8.0D0+8.0D0*Y)+
1259 2DL*(16.0D0*Y*Y)+
1260 3DLM*(4.0D0-4.0D0*Y-16.0D0*Y*Y/3.0D0+16.0D0/Y/3.0D0)+
1261 4S111MX*(8.0D0+8.0D0*Y)-104.0D0/3.0D0+80.0D0*Y/3.0D0-
1262 532.0D0*Y*Y/9.0D0+104.0D0/Y/9.0D0
1263 C2S=2.0D0*NF*CFTF*(A1+A2+A3+A4*DL2+A5*DL1)
1264 RETURN
1265 END
1266
1267
1268C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
1269C V + G -> Q + QB + G , F2(X,Q^2) , MSBAR SCHEME .
1270
1271
1272 DOUBLE PRECISION FUNCTION C2G2(Y)
1273 IMPLICIT REAL*8(A-H,O-Z)
1274 COMPLEX*16 WGPLG
1275 DOUBLE PRECISION NF
1276 COMMON/PAR/NF,Q2,FS,RS,DN
1277 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1278 FS2=FS*FS
1279 RS2=RS*RS
1280 DL1=DLOG(Q2/FS2)
1281 DL2=DL1*DL1
1282 EL1=DLOG(FS2/RS2)
1283 S121MX=WGPLG(1,2,1.0D0-Y)
1284 S12MX=WGPLG(1,2,-Y)
1285 S211MX=WGPLG(2,1,1.0D0-Y)
1286 S21MX=WGPLG(2,1,-Y)
1287 Z=(1.0D0-Y)/(1.0D0+Y)
1288 S21Z=WGPLG(2,1,Z)
1289 S21MZ=WGPLG(2,1,-Z)
1290 S111MX=WGPLG(1,1,1.0D0-Y)
1291 S11MX=WGPLG(1,1,-Y)
1292 ZETA2=1.6449340668D0
1293 ZETA3=1.2020569031D0
1294 DL=DLOG(Y)
1295 DLP=DLOG(1.0D0+Y)
1296 DLM=DLOG(1.0D0-Y)
1297 A1=DL**3*(-5.0D0/3.0D0-20.0D0*Y*Y/3.0D0+10.0D0*Y/3.0D0)+
1298 1DL**2*DLM*(8.0D0+24.0D0*Y*Y-16.0D0*Y)+
1299 2DL**2*DLP*(8.0D0+8.0D0*Y*Y+16.0D0*Y)+
1300 3DL**2*(-1.5D0-48.0D0*Y*Y*Y/5.0D0-36.0D0*Y*Y+22.0D0*Y/3.0D0)+
1301 4DL*DLM**2*(-8.0D0-20.0D0*Y*Y+16.0D0*Y)+
1302 5DL*DLM*(16.0D0+72.0D0*Y*Y-56.0D0*Y)+
1303 6DL*DLP**2*(-16.0D0-16.0D0*Y*Y-32.0D0*Y)+
1304 7DL*DLP*(48.0D0+96.0D0*Y*Y*Y/5.0D0+64.0D0*Y/3.0D0)
1305 A2=8.0D0*DL*DLP/Y/Y/15.0D0-8.0D0*DL/Y/15.0D0+
1306 1DL*ZETA2*(12.0D0+32.0D0*Y*Y-24.0D0*Y)+
1307 2DL*S111MX*(4.0D0+16.0D0*Y*Y-8.0D0*Y)+
1308 3DL*S11MX*(-16.0D0-16.0D0*Y*Y+32.0D0*Y)+
1309 4DL*(-236.0D0/15.0D0-216.0D0*Y*Y/5.0D0+113.0D0*Y/5.0D0)+
1310 5DLM**3*(10.0D0/3.0D0+20.0D0*Y*Y/3.0D0-20.0D0*Y/3.0D0)+
1311 6DLM**2*(-13.0D0-36.0D0*Y*Y+40.0D0*Y)+
1312 7ZETA2*DLM*(-16.0D0*Y*Y)
1313 A3=DLM*S111MX*(8.0D0+8.0D0*Y*Y-16.0D0*Y)+
1314 1DLM*(14.0D0+24.0D0*Y*Y-40.0D0*Y)+
1315 2ZETA2*DLP*(-16.0D0-16.0D0*Y*Y-32.0D0*Y)+
1316 3DLP*S11MX*(-32.0D0-32.0D0*Y*Y-64.0D0*Y)+
1317 4ZETA2*(26.0D0+96.0D0*Y*Y*Y/5.0D0+72.0D0*Y*Y-176.0D0*Y/3.0D0)+
1318 5ZETA3*(52.0D0+88.0D0*Y*Y-40.0D0*Y)
1319 A4=S111MX*(-10.0D0+24.0D0*Y)+
1320 1S211MX*(-16.0D0-24.0D0*Y*Y+32.0D0*Y)+
1321 2S121MX*(-4.0D0+8.0D0*Y*Y+8.0D0*Y)+
1322 3S11MX*(48.0D0+96.0D0*Y*Y*Y/5.0D0+64.0D0*Y/3.0D0)+
1323 48.0D0*S11MX/Y/Y/15.0D0+8.0D0/Y/15.0D0+
1324 5S21MX*(48.0D0+48.0D0*Y*Y-32.0D0*Y)+
1325 6S12MX*(-32.0D0-32.0D0*Y*Y-64.0D0*Y)-
1326 7647.0D0/15.0D0-36.0D0*Y*Y/5.0D0+239.0D0*Y/5.0D0
1327 A5=DL*(-2.0D0+4.0D0*Y-8.0D0*Y*Y)+
1328 1DLM*(4.0D0-8.0D0*Y+8.0D0*Y*Y)-
1329 21.0D0+4.0D0*Y
1330 A6=DL**2*(4.0D0-8.0D0*Y+16.0D0*Y*Y)+
1331 1DL*DLM*(-12.0D0+24.0D0*Y-32.0D0*Y*Y)+
1332 2DL*(4.0D0-24.0D0*Y+40.0D0*Y*Y)+
1333 3DLM**2*(8.0D0-16.0D0*Y+16.0D0*Y*Y)+
1334 4DLM*(-14.0D0+48.0D0*Y-40.0D0*Y*Y)+
1335 5ZETA2*(-16.0D0+32.0D0*Y-32.0D0*Y*Y)+
1336 6S111MX*(4.0D0-8.0D0*Y)+
1337 718.0D0-34.0D0*Y+8.0D0*Y*Y
1338 FCF=CF*(A1+A2+A3+A4+A5*DL2+A6*DL1)
1339 B1=DL**3*(10.0D0/3.0D0+28.0D0*Y/3.0D0)+
1340 1DL**2*DLM*(-4.0D0+8.0D0*Y*Y-32.0D0*Y)+
1341 2DL**2*DLP*(4.0D0+12.0D0*Y*Y+8.0D0*Y)+
1342 3DL**2*(-1.0D0-194.0D0*Y*Y/3.0D0+88.0D0*Y)+
1343 4DL*DLM**2*(-8.0D0*Y*Y+24.0D0*Y)+
1344 5DL*DLM*DLP*(-8.0D0-16.0D0*Y*Y-16.0D0*Y)+
1345 6DL*DLM*(8.0D0+148.0D0*Y*Y-144.0D0*Y)+
1346 8DL*DLP**2*(8.0D0+8.0D0*Y*Y+16.0D0*Y)
1347 B2=DL*DLP*(-24.0D0+80.0D0*Y*Y/3.0D0-16.0D0/Y/3.0D0)+
1348 1DL*ZETA2*(16.0D0*Y*Y-48.0D0*Y)+
1349 2DL*S111MX*(-8.0D0-16.0D0*Y)+
1350 3DL*S11MX*(8.0D0+8.0D0*Y*Y+16.0D0*Y)+
1351 4DL*(58.0D0-2090.0D0*Y*Y/9.0D0+584.0D0*Y/3.0D0)+
1352 5DLM**3*(2.0D0/3.0D0+4.0D0*Y*Y/3.0D0-4.0D0*Y/3.0D0)+
1353 6DLM**2*(-2.0D0-122.0D0*Y*Y/3.0D0+36.0D0*Y+8.0D0/Y/3.0D0)+
1354 7ZETA2*DLM*(-20.0D0-32.0D0*Y*Y+24.0D0*Y)
1355 B3=DLM*S111MX*(4.0D0-8.0D0*Y*Y+40.0D0*Y)+
1356 1DLM*S11MX*(-8.0D0-16.0D0*Y*Y-16.0D0*Y)+
1357 2DLM*(-62.0D0/3.0D0+1570.0D0*Y*Y/9.0D0-454.0D0*Y/3.0D0
1358 3+104.0D0/Y/9.0D0)+
1359 4ZETA2*DLP*(8.0D0+8.0D0*Y*Y+16.0D0*Y)+
1360 5DLP*S11MX*(16.0D0+16.0D0*Y*Y+32.0D0*Y)+
1361 6S21MZ*(-8.0D0-16.0D0*Y*Y-16.0D0*Y)+
1362 7S21Z*(8.0D0+16.0D0*Y*Y+16.0D0*Y)
1363 B4=ZETA2*(4.0D0+268.0D0*Y*Y/3.0D0-80.0D0*Y-32.0D0/Y/3.0D0)+
1364 1ZETA3*(-10.0D0-12.0D0*Y*Y-12.0D0*Y)+
1365 2S111MX*(4.0D0+176.0D0*Y*Y/3.0D0-64.0D0*Y+16.0D0/Y/3.0D0)+
1366 3S211MX*(-4.0D0+8.0D0*Y*Y-72.0D0*Y)+
1367 4S121MX*(4.0D0+8.0D0*Y)+
1368 5S11MX*(-24.0D0+80.0D0*Y*Y/3.0D0-16.0D0/Y/3.0D0)+
1369 6S21MX*(-8.0D0+8.0D0*Y*Y-16.0D0*Y)+
1370 7S12MX*(16.0D0+16.0D0*Y*Y+32.0D0*Y)+344.0D0/Y/27.0D0+
1371 8239.0D0/9.0D0-4493.0D0*Y*Y/27.0D0+1072.0D0*Y/9.0D0
1372C B4=B4+16.0D0*ZETA2/Y/3.0D0-344.0D0/Y/27.0D0
1373 B5=DL*(4.0D0+16.0D0*Y)+
1374 1DLM*(4.0D0-8.0D0*Y+8.0D0*Y*Y)+2.0D0+16.0D0*Y-
1375 262.0D0*Y*Y/3.0D0+8.0D0/Y/3.0D0
1376 B6=-22.0D0/3.0D0+44.0D0*Y/3.0D0-44.0D0*Y*Y/3.0D0
1377 B7=DL**2*(-8.0D0-24.0D0*Y)+
1378 1DL*DLM*(48.0D0*Y-16.0D0*Y*Y)+
1379 2DL*DLP*(-8.0D0-16.0D0*Y-16.0D0*Y*Y)+
1380 3DL*(-96.0D0*Y+100.0D0*Y*Y)+
1381 4DLM**2*(4.0D0-8.0D0*Y+8.0D0*Y*Y)+
1382 5DLM*(-4.0D0+80.0D0*Y-268.0D0*Y*Y/3.0D0+16.0D0/Y/3.0D0)+
1383 6ZETA2*(-8.0D0-16.0D0*Y*Y)+
1384 7S111MX*(8.0D0+32.0D0*Y)+
1385 8S11MX*(-8.0D0-16.0D0*Y-16.0D0*Y*Y)-110.0D0/3.0D0-
1386 9184.0D0*Y/3.0D0+814.0D0*Y*Y/9.0D0+104.0D0/Y/9.0D0
1387 B8=22.0D0*(1.0D0-2.0D0*Y+2.0D0*Y*Y)*(DL-DLM)/3.0D0+
1388 122.0D0/3.0D0-176.0D0*Y/3.0D0+176.0D0*Y*Y/3.0D0
1389 FCA=CA*(B1+B2+B3+B4+B5*DL2+B6*DL1*EL1+B7*DL1+B8*EL1)
1390 C1=4.0D0/3.0D0-8.0D0*Y/3.0D0+8.0D0*Y*Y/3.0D0
1391 C2=4.0D0*(1.0D0-2.0D0*Y+2.0D0*Y*Y)*(DLM-DL)/3.0D0-
1392 14.0D0/3.0D0+32.0D0*Y/3.0D0-32.0D0*Y*Y/3.0D0
1393 FNF=NF*(C1*DL1*EL1+C2*EL1)
1394 C2G2=NF*(FCF+FCA+FNF)
1395 RETURN
1396 END
1397
1398C THIS IS THE LONGITUDINAL WILSON COEFFICIENT CORRESPONDING
1399C TO THE PROCESS : V + Q -> Q + G (SINGLET AND NONSINGLET)
1400C FL(X,Q^2 , MSBAR SCHEME .
1401
1402 DOUBLE PRECISION FUNCTION CLNS1(Y)
1403 IMPLICIT REAL*8(A-H,O-Z)
1404 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1405 A1=4.0D0*Y
1406 CLNS1=CF*A1
1407 RETURN
1408 END
1409
1410
1411C THIS IS THE LONGITUDINAL WILSON COEFFICIENT CORRESPONDING
1412C TO THE PROCESS : V + G -> Q + QB
1413C FL(X,Q^2) , MSBAR SCHEME .
1414
1415 DOUBLE PRECISION FUNCTION CLG1(Y)
1416 IMPLICIT REAL*8(A-H,O-Z)
1417 DOUBLE PRECISION NF
1418 COMMON/PAR/NF,Q2,FS,RS,DN
1419 A1=8.0D0*Y-8.0D0*Y*Y
1420 CLG1=NF*A1
1421 RETURN
1422 END
1423
1424
1425C THIS IS THE LONGITUDINAL WILSON COEFFICIENT CORRESPONDING TO THE PROCESSES :
1426C V + Q -> Q + G + G ; V + Q -> Q + Q + QB WITH NONIDENTICAL QUARKS
1427C IN THE FINAL STATE (NONSINGLET) , FL(X,Q^2) , MSBAR SCHEME .
1428
1429
1430
1431 DOUBLE PRECISION FUNCTION CLNSP(Y)
1432 IMPLICIT REAL*8(A-H,O-Z)
1433 COMPLEX*16 WGPLG
1434 DOUBLE PRECISION NF
1435 COMMON/PAR/NF,Q2,FS,RS,DN
1436 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1437 FS2=FS*FS
1438 RS2=RS*RS
1439 DL1=DLOG(Q2/FS2)
1440 CL1=DLOG(Q2/RS2)
1441 S121MX=WGPLG(1,2,1.0D0-Y)
1442 S21MX=WGPLG(2,1,-Y)
1443 S111MX=WGPLG(1,1,1.0D0-Y)
1444 S11MX=WGPLG(1,1,-Y)
1445 ZETA2=1.6449340668D0
1446 DL=DLOG(Y)
1447 DLP=DLOG(1.0D0+Y)
1448 DLM=DLOG(1.0D0-Y)
1449 A1=DL**2*(24.0D0*Y*Y*Y/5.0D0+16.0D0*Y*Y)+
1450 1DL*DLM*(-24.0D0*Y)+
1451 2DL*DLP*(32.0D0-48.0D0*Y*Y*Y/5.0D0-32.0D0*Y*Y+16.0D0/Y)+
1452 332.0D0*DL*DLP/Y/Y/5.0D0-32.0D0*DL/Y/5.0D0+
1453 4DL*S11MX*(32.0D0*Y)+
1454 5DL*(-24.0D0/5.0D0+48.0D0*Y*Y/5.0D0+16.0D0*Y/5.0D0)+
1455 6DLM**2*(8.0D0*Y)+
1456 7ZETA2*DLM*(-32.0D0*Y)
1457 A2=DLM*(8.0D0+28.0D0*Y)+
1458 1ZETA2*(-48.0D0*Y*Y*Y/5.0D0-32.0D0*Y*Y)+
1459 2S111MX*(-8.0D0*Y)+
1460 3S121MX*(32.0D0*Y)+
1461 4S11MX*(32.0D0-48.0D0*Y*Y*Y/5.0D0-32.0D0*Y*Y+16.0D0/Y)+
1462 532.0D0*S11MX/Y/Y/5.0D0+32.0D0/Y/5.0D0+
1463 6S21MX*(-64.0D0*Y)+252.0D0/5.0D0+
1464 748.0D0*Y*Y/5.0D0-502.0D0*Y/5.0D0
1465 A3=DL*(-8.0D0*Y)+
1466 1DLM*(16.0D0*Y)+
1467 28.0D0+4.0D0*Y
1468 FCF=CF2*(A1+A2+A3*DL1)
1469 B1=DL**2*(-12.0D0*Y*Y*Y/5.0D0-8.0D0*Y*Y+4.0D0*Y)+
1470 1DL*DLP*(-16.0D0+24.0D0*Y*Y*Y/5.0D0+16.0D0*Y*Y-8.0D0/Y)
1471 2-16.0D0*DL*DLP/Y/Y/5.0D0+16.0D0*DL/Y/5.0D0+
1472 3DL*S11MX*(-16.0D0*Y)+
1473 4DL*(-8.0D0/5.0D0-24.0D0*Y*Y/5.0D0+416.0D0*Y/15.0D0)+
1474 5ZETA2*DLM*(16.0D0*Y)+
1475 6DLM*(-92.0D0*Y/3.0D0)
1476 B2=ZETA2*(24.0D0*Y*Y*Y/5.0D0+16.0D0*Y*Y-8.0D0*Y)+
1477 1S121MX*(-16.0D0*Y)+
1478 2S11MX*(-16.0D0+24.0D0*Y*Y*Y/5.0D0+16.0D0*Y*Y-8.0D0/Y)
1479 3-16.0D0*S11MX/Y/Y/5.0D0-16.0D0/Y/5.0D0+
1480 4S21MX*(32.0D0*Y)-508.0D0/15.0D0
1481 5-24.0D0*Y*Y/5.0D0+4034.0D0*Y/45.0D0
1482 B3=-44.0D0*Y/3.0D0
1483 FCA=CACF*(B1+B2+B3*CL1)
1484 C1=DL*(-16.0D0*Y/3.0)+
1485 1DLM*(8.0D0*Y/3.0D0)+
1486 28.0D0/3.0D0-100.0D0*Y/9.0D0
1487 C2=8.0D0*Y/3.0D0
1488 FNF=2.0D0*NF*CFTF*(C1+C2*CL1)
1489 CLNSP=FCF+FCA+FNF
1490 RETURN
1491 END
1492
1493C THIS IS THE LONGITUDINAL WILSON COEFFICIENT CORRESPONDING TO THE PROCESS:
1494C V + Q -> Q + Q + QB WITH IDENTICAL QUARKS IN THE FINAL STATE
1495C (NONSINGLET) , FL(X,Q^2) , MSBAR AND DIS SCHEME .
1496
1497
1498 DOUBLE PRECISION FUNCTION CLNSM(Y)
1499 IMPLICIT REAL*8(A-H,O-Z)
1500 COMPLEX*16 WGPLG
1501 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1502 S12MX=WGPLG(1,2,-Y)
1503 S21MX=WGPLG(2,1,-Y)
1504 S11MX=WGPLG(1,1,-Y)
1505 ZETA2=1.6449340668D0
1506 ZETA3=1.2020569031D0
1507 DL=DLOG(Y)
1508 DLP=DLOG(1.0D0+Y)
1509 A1=DL**2*DLP*(16.0D0*Y)+
1510 1DL**2*(24.0D0*Y*Y*Y/5.0D0-16.0D0*Y*Y-8.0D0*Y)+
1511 2DL*DLP**2*(-32.0D0*Y)+
1512 3DL*DLP*(32.0D0-48.0D0*Y*Y*Y/5.0D0+32.0D0*Y*Y+32.0D0*Y-16.0D0/Y)+
1513 432.0D0*DL*DLP/Y/Y/5.0D0-32.0D0*DL/Y/5.0D0+
1514 5DL*(16.0D0/5.0D0+48.0D0*Y*Y/5.0D0-224.0D0*Y/5.0D0)+
1515 6ZETA2*DLP*(-32.0D0*Y)+
1516 7DLP*S11MX*(-64.0D0*Y)
1517 A2=ZETA2*(-48.0D0*Y*Y*Y/5.0D0+32.0D0*Y*Y+16.0D0*Y)+
1518 1ZETA3*(32.0D0*Y)+
1519 2S11MX*(32.0D0-48.0D0*Y*Y*Y/5.0D0+32.0D0*Y*Y+32.0D0*Y-16.0D0/Y)+
1520 332.0D0*S11MX/Y/Y/5.0D0+32.0D0/Y/5.0D0+
1521 4S21MX*(32.0D0*Y)+
1522 5S12MX*(-64.0D0*Y)-208.0D0/5.0D0+
1523 648.0D0*Y*Y/5.0D0+128.0D0*Y/5.0D0
1524 CLNSM=(CF2-0.5D0*CACF)*(A1+A2)
1525 RETURN
1526 END
1527
1528C THIS IS THE LONGITUDINAL WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
1529C V + Q -> Q + QB + Q (SINGLET) , FL(X,Q^2) , MSBAR SCHEME .
1530
1531
1532 DOUBLE PRECISION FUNCTION CLS(Y)
1533 IMPLICIT REAL*8(A-H,O-Z)
1534 COMPLEX*16 WGPLG
1535 DOUBLE PRECISION NF
1536 COMMON/PAR/NF,Q2,FS,RS,DN
1537 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1538 S111MX=WGPLG(1,1,1.0D0-Y)
1539 DL=DLOG(Y)
1540 DLM=DLOG(1.0D0-Y)
1541 FS2=FS*FS
1542 DL1=DLOG(Q2/FS2)
1543 A1=DL**2*(16.0D0*Y)+
1544 1DL*DLM*(-16.0D0*Y)+
1545 2DL*(16.0D0-32.0D0*Y*Y-16.0D0*Y)+
1546 3DLM*(-16.0D0+32.0D0*Y*Y/3.0D0+16.0D0/Y/3.0D0)+
1547 4S111MX*(-16.0D0*Y)+16.0D0/3.0D0+
1548 5160.0D0*Y*Y/9.0D0-64.0D0*Y/3.0D0-16.0D0/Y/9.0D0
1549 A2=DL*(-16.0D0*Y)-
1550 116.0D0+32.0D0*Y*Y/3.0D0+16.0D0/Y/3.0D0
1551 CLS=2.0D0*NF*CFTF*(A1+A2*DL1)
1552 RETURN
1553 END
1554
1555C THIS IS THE LONGITUDINAL WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
1556C V + G -> Q + QB + G , FL(X,Q^2) , MSBAR SCHEME .
1557
1558
1559 DOUBLE PRECISION FUNCTION CLG2(Y)
1560 IMPLICIT REAL*8(A-H,O-Z)
1561 COMPLEX*16 WGPLG
1562 DOUBLE PRECISION NF
1563 COMMON/PAR/NF,Q2,FS,RS,DN
1564 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1565 FS2=FS*FS
1566 RS2=RS*RS
1567 DL1=DLOG(Q2/FS2)
1568 EL1=DLOG(FS2/RS2)
1569 S111MX=WGPLG(1,1,1.0D0-Y)
1570 S11MX=WGPLG(1,1,-Y)
1571 ZETA2=1.6449340668D0
1572 DL=DLOG(Y)
1573 DLP=DLOG(1.0D0+Y)
1574 DLM=DLOG(1.0D0-Y)
1575 A1=DL**2*(-32.0D0*Y*Y*Y/5.0D0-32.0D0*Y/3.0D0)+
1576 1DL*DLM*(16.0D0*Y)+
1577 2DL*DLP*(64.0D0*Y*Y*Y/5.0D0-32.0D0*Y/3.0D0)+
1578 332.0D0*DL*DLP/Y/Y/15.0D0-32.0D0*DL/Y/15.0D0+
1579 4DL*(-104.0D0/15.0D0+96.0D0*Y*Y/5.0D0-208.0D0*Y/5.0D0)+
1580 5DLM*(8.0D0-32.0D0*Y*Y+24.0D0*Y)
1581 A2=ZETA2*(64.0D0*Y*Y*Y/5.0D0-32.0D0*Y/3.0D0)+
1582 1S111MX*(16.0D0*Y)+
1583 2S11MX*(64.0D0*Y*Y*Y/5.0D0-32.0D0*Y/3.0D0)+
1584 332.0D0*S11MX/Y/Y/15.0D0+32.0D0/Y/15.0D0
1585 4-128.0D0/15.0D0+336.0D0*Y*Y/5.0D0-304.0D0*Y/5.0D0
1586 A3=DL*(16.0D0*Y)+
1587 18.0D0+8.0D0*Y-16.0D0*Y*Y
1588 FCF=CF*(A1+A2+A3*DL1)
1589 B1=DL**2*(48.0D0*Y)+
1590 1DL*DLM*(32.0D0*Y*Y-96.0D0*Y)+
1591 2DL*DLP*(32.0D0*Y*Y+32.0D0*Y)+
1592 3DL*(16.0D0-208.0D0*Y*Y+128.0D0*Y)+
1593 4DLM**2*(-16.0D0*Y*Y+16.0D0*Y)+
1594 5DLM*(-16.0D0+464.0D0*Y*Y/3.0D0-144.0D0*Y+16.0D0/Y/3.0D0)
1595 B2=ZETA2*(32.0D0*Y*Y)+
1596 1S111MX*(-64.0D0*Y)+
1597 2S11MX*(32.0D0*Y*Y+32.0D0*Y)+16.0D0/3.0D0
1598 3-848.0D0*Y*Y/9.0D0+272.0D0*Y/3.0D0-16.0D0/Y/9.0D0
1599 B3=DL*(-64.0D0*Y)+
1600 1DLM*(32.0D0*Y-32.0D0*Y*Y)-
1601 216.0D0-80.0D0*Y+272.0D0*Y*Y/3.0D0+16.0D0/Y/3.0D0
1602 B4=-88.0D0*Y/3.0D0+88.0D0*Y*Y/3.0D0
1603 FCA=CA*(B1+B2+B3*DL1+B4*EL1)
1604 C1=16.0D0*Y/3.0D0-16.0D0*Y*Y/3.0D0
1605 FNF=NF*C1*EL1
1606 CLG2=NF*(FCF+FCA+FNF)
1607 RETURN
1608 END
1609
1610C THIS IS THE LONGITUDINAL WILSON COEFFICIENT CORRESPONDING TO THE PROCESSES :
1611C V + Q -> Q + G + G ; V + Q -> Q + Q + QB WITH NONIDENTICAL QUARKS
1612C IN THE FINAL STATE (NONSINGLET) , FL(X,Q^2) , DIS SCHEME .
1613
1614
1615
1616 DOUBLE PRECISION FUNCTION DLNSP(Y)
1617 IMPLICIT REAL*8(A-H,O-Z)
1618 COMPLEX*16 WGPLG
1619 DOUBLE PRECISION NF
1620 COMMON/PAR/NF,Q2,FS,RS,DN
1621 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1622 FS2=FS*FS
1623 RS2=RS*RS
1624 DL1=DLOG(Q2/FS2)
1625 CL1=DLOG(Q2/RS2)
1626 S121MX=WGPLG(1,2,1.0D0-Y)
1627 S21MX=WGPLG(2,1,-Y)
1628 S111MX=WGPLG(1,1,1.0D0-Y)
1629 S11MX=WGPLG(1,1,-Y)
1630 ZETA2=1.6449340668D0
1631 DL=DLOG(Y)
1632 DLP=DLOG(1.0D0+Y)
1633 DLM=DLOG(1.0D0-Y)
1634 A1=DL**2*(24.0D0*Y*Y*Y/5.0D0+16.0D0*Y*Y-4.0D0*Y)+
1635 1DL*DLM*(-16.0D0*Y)+
1636 2DL*DLP*(32.0D0-48.0D0*Y*Y*Y/5.0D0-32.0D0*Y*Y+16.0D0/Y)+
1637 332.0D0*DL*DLP/Y/Y/5.0D0-32.0D0*DL/Y/5.0D0+
1638 4DL*S11MX*(32.0D0*Y)+
1639 5DL*(16.0D0/5.0D0+48.0D0*Y*Y/5.0D0-4.0D0*Y/5.0D0)+
1640 7ZETA2*DLM*(-32.0D0*Y)
1641 A2=DLM*(48.0D0*Y)+
1642 1ZETA2*(-48.0D0*Y*Y*Y/5.0D0-32.0D0*Y*Y+16.0D0*Y)+
1643 2S111MX*(-16.0D0*Y)+
1644 3S121MX*(32.0D0*Y)+
1645 4S11MX*(32.0D0-48.0D0*Y*Y*Y/5.0D0-32.0D0*Y*Y+16.0D0/Y)+
1646 532.0D0*S11MX/Y/Y/5.0D0+32.0D0/Y/5.0D0+
1647 6S21MX*(-64.0D0*Y)+232.0D0/5.0D0+
1648 748.0D0*Y*Y/5.0D0-302.0D0*Y/5.0D0
1649 A3=DL*(-8.0D0*Y)+
1650 1DLM*(16.0D0*Y)+
1651 28.0D0+4.0D0*Y
1652 FCF=CF2*(A1+A2+A3*DL1)
1653 B1=DL**2*(-12.0D0*Y*Y*Y/5.0D0-8.0D0*Y*Y+4.0D0*Y)+
1654 1DL*DLP*(-16.0D0+24.0D0*Y*Y*Y/5.0D0+16.0D0*Y*Y-8.0D0/Y)
1655 2-16.0D0*DL*DLP/Y/Y/5.0D0+16.0D0*DL/Y/5.0D0+
1656 3DL*S11MX*(-16.0D0*Y)+
1657 4DL*(-8.0D0/5.0D0-24.0D0*Y*Y/5.0D0+416.0D0*Y/15.0D0)+
1658 5ZETA2*DLM*(16.0D0*Y)+
1659 6DLM*(-92.0D0*Y/3.0D0)
1660 B2=ZETA2*(24.0D0*Y*Y*Y/5.0D0+16.0D0*Y*Y-8.0D0*Y)+
1661 1S121MX*(-16.0D0*Y)+
1662 2S11MX*(-16.0D0+24.0D0*Y*Y*Y/5.0D0+16.0D0*Y*Y-8.0D0/Y)
1663 3-16.0D0*S11MX/Y/Y/5.0D0-16.0D0/Y/5.0D0+
1664 4S21MX*(32.0D0*Y)-508.0D0/15.0D0
1665 5-24.0D0*Y*Y/5.0D0+4034.0D0*Y/45.0D0
1666 B3=-44.0D0*Y/3.0D0
1667 FCA=CACF*(B1+B2+B3*CL1)
1668 C1=DL*(-16.0D0*Y/3.0)+
1669 1DLM*(8.0D0*Y/3.0D0)+
1670 28.0D0/3.0D0-100.0D0*Y/9.0D0
1671 C2=8.0D0*Y/3.0D0
1672 FNF=2.0D0*NF*CFTF*(C1+C2*CL1)
1673 DLNSP=FCF+FCA+FNF
1674 RETURN
1675 END
1676
1677
1678C THIS IS THE LONGITUDINAL WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
1679C V + Q -> Q + QB + Q (SINGLET) , FL(X,Q^2) , DIS SCHEME .
1680
1681
1682
1683 DOUBLE PRECISION FUNCTION DLS(Y)
1684 IMPLICIT REAL*8(A-H,O-Z)
1685 COMPLEX*16 WGPLG
1686 DOUBLE PRECISION NF
1687 COMMON/PAR/NF,Q2,FS,RS,DN
1688 COMMON/DIS/NDIS
1689 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1690 S111MX=WGPLG(1,1,1.0D0-Y)
1691 ZETA2=1.6449340668D0
1692 DL=DLOG(Y)
1693 DLM=DLOG(1.0D0-Y)
1694 FS2=FS*FS
1695 DL1=DLOG(Q2/FS2)
1696 IF (NDIS.EQ.1) THEN
1697 A1=DL**2*(-16.0D0*Y*Y+24.0D0*Y)+
1698 1DL*DLM*(32.0D0*Y*Y-32.0D0*Y)+
1699 2DL*(8.0D0-80.0D0*Y*Y+8.0D0*Y)+
1700 3DLM**2*(-16.0D0*Y*Y+16.0D0*Y)+
1701 4DLM*(-8.0D0+176.0D0*Y*Y/3.0D0-56.0D0*Y+16.0D0/Y/3.0D0)+
1702 5ZETA2*(32.0D0*Y*Y-32.0D0*Y)+16.0D0/3.0D0+
1703 5736.0D0*Y*Y/9.0D0-256.0D0*Y/3.0D0-16.0D0/Y/9.0D0
1704 A2=DL*(-16.0D0*Y)-
1705 116.0D0+32.0D0*Y*Y/3.0D0+16.0D0/Y/3.0D0
1706 ELSE
1707 A1=DL**2*(16.0D0*Y)+
1708 1DL*DLM*(-16.0D0*Y)+
1709 2DL*(16.0D0-32.0D0*Y*Y-16.0D0*Y)+
1710 3DLM*(-16.0D0+32.0D0*Y*Y/3.0D0+16.0D0/Y/3.0D0)+
1711 4S111MX*(-16.0D0*Y)+16.0D0/3.0D0+
1712 5160.0D0*Y*Y/9.0D0-64.0D0*Y/3.0D0-16.0D0/Y/9.0D0
1713 A2=DL*(-16.0D0*Y)-
1714 116.0D0+32.0D0*Y*Y/3.0D0+16.0D0/Y/3.0D0
1715 ENDIF
1716 DLS=2.0D0*NF*CFTF*(A1+A2*DL1)
1717 RETURN
1718 END
1719
1720
1721
1722
1723C THIS IS THE LONGITUDINAL WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
1724C V + G -> Q + QB + G , FL(X,Q^2) , DIS SCHEME .
1725
1726
1727 DOUBLE PRECISION FUNCTION DLG2(Y)
1728 IMPLICIT REAL*8(A-H,O-Z)
1729 COMPLEX*16 WGPLG
1730 DOUBLE PRECISION NF
1731 COMMON/PAR/NF,Q2,FS,RS,DN
1732 COMMON/DIS/NDIS
1733 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1734 FS2=FS*FS
1735 RS2=RS*RS
1736 DL1=DLOG(Q2/FS2)
1737 EL1=DLOG(FS2/RS2)
1738 S111MX=WGPLG(1,1,1.0D0-Y)
1739 S11MX=WGPLG(1,1,-Y)
1740 ZETA2=1.6449340668D0
1741 DL=DLOG(Y)
1742 DLP=DLOG(1.0D0+Y)
1743 DLM=DLOG(1.0D0-Y)
1744 A1=DL**2*(-32.0D0*Y*Y*Y/5.0D0-8.0D0*Y/3.0D0)+
1745 1DL*DLP*(64.0D0*Y*Y*Y/5.0D0-32.0D0*Y/3.0D0)+
1746 232.0D0*DL*DLP/Y/Y/15.0D0-32.0D0*DL/Y/15.0D0+
1747 3DL*(16.0D0/15.0D0+16.0D0*Y*Y/5.0D0+72.0D0*Y/5.0D0)+
1748 4DLM*(-16.0D0*Y*Y+16.0D0*Y)
1749 A2=ZETA2*(64.0D0*Y*Y*Y/5.0D0-32.0D0*Y/3.0D0)+
1750 1S11MX*(64.0D0*Y*Y*Y/5.0D0-32.0D0*Y/3.0D0)+
1751 232.0D0*S11MX/Y/Y/15.0D0+32.0D0/Y/15.0D0
1752 3+112.0D0/15.0D0+16.0D0*Y*Y/5.0D0-64.0D0*Y/5.0D0
1753 A3=DL*(16.0D0*Y)+
1754 18.0D0+8.0D0*Y-16.0D0*Y*Y
1755 FCF=CF*(A1+A2+A3*DL1)
1756 IF (NDIS.EQ.1) THEN
1757 B1=DL**2*(48.0D0*Y)+
1758 1DL*DLM*(32.0D0*Y*Y-96.0D0*Y)+
1759 2DL*DLP*(32.0D0*Y*Y+32.0D0*Y)+
1760 3DL*(16.0D0-208.0D0*Y*Y+128.0D0*Y)+
1761 4DLM**2*(-16.0D0*Y*Y+16.0D0*Y)+
1762 5DLM*(-16.0D0+464.0D0*Y*Y/3.0D0-144.0D0*Y+16.0D0/Y/3.0D0)
1763 B2=ZETA2*(32.0D0*Y*Y)+
1764 1S111MX*(-64.0D0*Y)+
1765 2S11MX*(32.0D0*Y*Y+32.0D0*Y)+16.0D0/3.0D0
1766 3-848.0D0*Y*Y/9.0D0+272.0D0*Y/3.0D0-16.0D0/Y/9.0D0
1767 B3=DL*(-64.0D0*Y)+
1768 1DLM*(32.0D0*Y-32.0D0*Y*Y)-
1769 216.0D0-80.0D0*Y+272.0D0*Y*Y/3.0D0+16.0D0/Y/3.0D0
1770 B4=-88.0D0*Y/3.0D0+88.0D0*Y*Y/3.0D0
1771 FCA=CA*(B1+B2+B3*DL1+B4*EL1)
1772 C1=DL**2*(-16.0D0*Y-16.0D0*Y*Y)+
1773 1DL*DLM*(32.0D0*Y+32.0D0*Y*Y)+
1774 2DL*(-8.0D0-144.0D0*Y-72.0D0*Y*Y)+
1775 3DLM*(8.0D0+48.0D0*Y-56.0D0*Y*Y)+
1776 4S111MX*(32.0D0*Y+32.0D0*Y*Y)-
1777 520.0D0-248.0D0*Y+268.0D0*Y*Y
1778 C2=16.0D0*Y/3.0D0-16.0D0*Y*Y/3.0D0
1779 FNF=NF*(C1+C2*EL1)
1780 ELSE
1781 B1=DL**2*(48.0D0*Y)+
1782 1DL*DLM*(32.0D0*Y*Y-96.0D0*Y)+
1783 2DL*DLP*(32.0D0*Y*Y+32.0D0*Y)+
1784 3DL*(16.0D0-208.0D0*Y*Y+128.0D0*Y)+
1785 4DLM**2*(-16.0D0*Y*Y+16.0D0*Y)+
1786 5DLM*(-16.0D0+464.0D0*Y*Y/3.0D0-144.0D0*Y+16.0D0/Y/3.0D0)
1787 B2=ZETA2*(32.0D0*Y*Y)+
1788 1S111MX*(-64.0D0*Y)+
1789 2S11MX*(32.0D0*Y*Y+32.0D0*Y)+16.0D0/3.0D0
1790 3-848.0D0*Y*Y/9.0D0+272.0D0*Y/3.0D0-16.0D0/Y/9.0D0
1791 B3=DL*(-64.0D0*Y)+
1792 1DLM*(32.0D0*Y-32.0D0*Y*Y)-
1793 216.0D0-80.0D0*Y+272.0D0*Y*Y/3.0D0+16.0D0/Y/3.0D0
1794 B4=-88.0D0*Y/3.0D0+88.0D0*Y*Y/3.0D0
1795 CA=3.0D0
1796 FCA=CA*(B1+B2+B3*DL1+B4*EL1)
1797 C1=16.0D0*Y/3.0D0-16.0D0*Y*Y/3.0D0
1798 FNF=NF*C1*EL1
1799 ENDIF
1800 DLG2=NF*(FCF+FCA+FNF)
1801 RETURN
1802 END
1803
1804
1805C THIS IS THE POLARIZED WILSON COEFFICIENT CORRESPONDING TO
1806C THE PROCESS : V + G -> Q + QB , G1(X,Q^2) , MSBAR SCHEME .
1807
1808
1809 REAL*8 FUNCTION POLG1(Y)
1810 IMPLICIT REAL*8(A-H,O-Z)
1811 DOUBLE PRECISION NF
1812 COMMON/PAR/NF,Q2,FS,RS,DN
1813 COMMON/DIS/NDIS
1814 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1815 FS2=FS*FS
1816 DL1=DLOG(Q2/FS2)
1817 DL=DLOG(Y)
1818 DLM=DLOG(1.0D0-Y)
1819 A1=(8.0D0*Y-4.0D0)*(DL1+DLM-DL)+4.0D0*(3.0D0-4.0D0*Y)
1820C THE 0.5 BELOW IS DUE TO TF=1/2
1821 POLG1=0.5D0*NF*A1
1822 RETURN
1823 END
1824
1825C THIS IS THE POLARIZED WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
1826C V + Q -> Q + QB + Q (SINGLET) , G1(X,Q^2) , MSBAR SCHEME .
1827
1828
1829 REAL*8 FUNCTION POLS(Y)
1830 IMPLICIT REAL*8(A-H,O-Z)
1831 COMPLEX*16 WGPLG
1832 DOUBLE PRECISION NF
1833 COMMON/PAR/NF,Q2,FS,RS,DN
1834 COMMON/DIS/NDIS
1835 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1836 FS2=FS*FS
1837 DL1=DLOG(Q2/FS2)
1838 DL2=DL1*DL1
1839 S211MX=WGPLG(2,1,1.0D0-Y)
1840 S111MX=WGPLG(1,1,1.0D0-Y)
1841 S11MX=WGPLG(1,1,-Y)
1842 ZETA2=1.6449340668D0
1843 DL=DLOG(Y)
1844 DLP=DLOG(1.0D0+Y)
1845 DLM=DLOG(1.0D0-Y)
1846 A1=DL**3*(20.0D0/3.0D0+20.0D0*Y/3.0D0)+
1847 1DL**2*DLM*(-16.0D0-16.0D0*Y)+
1848 2DL**2*(50.0D0+16.0D0*Y*Y/3.0D0-10.0D0*Y)+
1849 3DL*DLM**2*(8.0D0+8.0D0*Y)+
1850 4DL*DLM*(-64.0D0+32.0D0*Y)+
1851 5DL*DLP*(-32.0D0-32.0D0*Y*Y/3.0D0-32.0D0*Y-32.0D0/Y/3.0D0)
1852 A2=DL*ZETA2*(-16.0D0-16.0D0*Y)+
1853 1DL*S111MX*(-16.0D0-16.0D0*Y)+
1854 2DL*(476.0D0/3.0D0-52.0D0*Y/3.0D0)+
1855 3DLM**2*(20.0D0-20.0D0*Y)+
1856 5DLM*S111MX*(16.0D0+16.0D0*Y)+
1857 6DLM*(-88.0D0+88.0D0*Y)
1858 A3=ZETA2*(-72.0D0-32.0D0*Y*Y/3.0D0+40.0D0*Y)+
1859 1S111MX*(-24.0D0-8.0D0*Y)+
1860 2S211MX*(-16.0D0-16.0D0*Y)+
1861 3S11MX*(-32.0D0-32.0D0*Y*Y/3.0D0-32.0D0*Y-32.0D0/Y/3.0D0)+
1862 4760.0D0/3.0D0-760.0D0*Y/3.0D0
1863 A4=DL*(8.0D0+8.0D0*Y)+
1864 120.0D0-20.0D0*Y
1865 A5=DL**2*(-16.0D0-16.0D0*Y)+
1866 1DL*DLM*(16.0D0+16.0D0*Y)+
1867 2DL*(-64.0D0+32.0D0*Y)+
1868 3DLM*(40.0D0-40.0D0*Y)+
1869 4S111MX*(16.0D0+16.0D0*Y)-
1870 588.0D0+88.0D0*Y
1871 POLS=NF*CFTF*(A1+A2+A3+A4*DL2+A5*DL1)
1872 RETURN
1873 END
1874
1875C THIS IS THE POLARIZED WILSON COEFFICIENT CORRESPONDING TO
1876C THE PROCESS : V + G -> Q + QB + G , G1(X,Q^2) , MSBAR SCHEME .
1877
1878
1879 REAL*8 FUNCTION POLG2(Y)
1880 IMPLICIT REAL*8(A-H,O-Z)
1881 COMPLEX*16 WGPLG
1882 DOUBLE PRECISION NF
1883 COMMON/PAR/NF,Q2,FS,RS,DN
1884 COMMON/DIS/NDIS
1885 COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1886 FS2=FS*FS
1887 RS2=RS*RS
1888 DL1=DLOG(Q2/FS2)
1889 DL2=DL1*DL1
1890 EL1=DLOG(FS2/RS2)
1891 S121MX=WGPLG(1,2,1.0D0-Y)
1892 S12MX=WGPLG(1,2,-Y)
1893 S211MX=WGPLG(2,1,1.0D0-Y)
1894 S21MX=WGPLG(2,1,-Y)
1895 Z=(1.0D0-Y)/(1.0D0+Y)
1896 S21Z=WGPLG(2,1,Z)
1897 S21MZ=WGPLG(2,1,-Z)
1898 S111MX=WGPLG(1,1,1.0D0-Y)
1899 S11MX=WGPLG(1,1,-Y)
1900 ZETA2=1.6449340668D0
1901 ZETA3=1.2020569031D0
1902 DL=DLOG(Y)
1903 DLP=DLOG(1.0D0+Y)
1904 DLM=DLOG(1.0D0-Y)
1905 A1=DL**3*(10.0D0/3.0D0-20.0D0*Y/3.0D0)+
1906 1DL**2*DLM*(-16.0D0+32.0D0*Y)+
1907 2DL**2*DLP*(16.0D0+16.0D0*Y*Y+32.0D0*Y)+
1908 3DL**2*(41.0D0-104.0D0*Y*Y/3.0D0-16.0D0*Y)+
1909 4DL*DLM**2*(16.0D0-32.0D0*Y)+
1910 5DL*DLM*(-88.0D0+96.0D0*Y)+
1911 6DL*DLP**2*(-32.0D0-32.0D0*Y*Y-64.0D0*Y)+
1912 7DL*DLP*(208.0D0*Y*Y/3.0D0+64.0D0*Y+64.0D0/Y/3.0D0)
1913 A2=DL*ZETA2*(-24.0D0+48.0D0*Y)+
1914 2DL*S111MX*(-8.0D0+16.0D0*Y)+
1915 3DL*S11MX*(32.0D0+32.0D0*Y*Y-64.0D0*Y)+
1916 4DL*(416.0D0/3.0D0-16.0D0*Y*Y-274.0D0*Y/3.0D0)+
1917 5DLM**3*(-20.0D0/3.0D0+40.0D0*Y/3.0D0)+
1918 6DLM**2*(54.0D0-72.0D0*Y)+
1919 7ZETA2*DLM*(-32.0D0*Y*Y)
1920 A3=DLM*S111MX*(-16.0D0+32.0D0*Y)+
1921 1DLM*(-124.0D0+16.0D0*Y*Y+104.0D0*Y)+
1922 2ZETA2*DLP*(-32.0D0-32.0D0*Y*Y-64.0D0*Y)+
1923 3DLP*S11MX*(-64.0D0-64.0D0*Y*Y-128.0D0*Y)+
1924 4ZETA2*(-108.0D0+208.0D0*Y*Y/3.0D0+144.0D0*Y)+
1925 5ZETA3*(-40.0D0+32.0D0*Y*Y+208.0D0*Y)
1926 A4=S111MX*(20.0D0-48.0D0*Y)+
1927 1S211MX*(32.0D0-64.0D0*Y)+
1928 2S121MX*(8.0D0+32.0D0*Y*Y-16.0D0*Y)+
1929 3S11MX*(208.0D0*Y*Y/3.0D0+64.0D0*Y+64.0D0/Y/3.0D0)+
1930 4S21MX*(-32.0D0-32.0D0*Y*Y+192.0D0*Y)+
1931 5S12MX*(-64.0D0-64.0D0*Y*Y-128.0D0*Y)+
1932 6466.0D0/3.0D0-478.0D0*Y/3.0D0
1933 A5=DL*(4.0D0-8.0D0*Y)+
1934 1DLM*(-8.0D0+16.0D0*Y)+
1935 26.0D0
1936 A6=DL**2*(-8.0D0+16.0D0*Y)+
1937 1DL*DLM*(24.0D0-48.0D0*Y)+
1938 2DL*(-48.0D0+32.0D0*Y)+
1939 3DLM**2*(-16.0D0+32.0D0*Y)+
1940 4DLM*(68.0D0-80.0D0*Y)+
1941 5ZETA2*(32.0D0-64.0D0*Y)+
1942 6S111MX*(-8.0D0+16.0D0*Y)-
1943 768.0D0+52.0D0*Y
1944 FCF=CF*(A1+A2+A3+A4+A5*DL2+A6*DL1)
1945 B1=DL**3*(28.0D0/3.0D0+40.0D0*Y/3.0D0)+
1946 1DL**2*DLM*(-32.0D0-16.0D0*Y)+
1947 2DL**2*DLP*(8.0D0-8.0D0*Y*Y+16.0D0*Y)+
1948 3DL**2*(90.0D0+88.0D0*Y*Y/3.0D0-32.0D0*Y)+
1949 4DL*DLM**2*(24.0D0)+
1950 5DL*DLM*DLP*(-16.0D0-32.0D0*Y)+
1951 6DL*DLM*(-136.0D0+128.0D0*Y)+
1952 8DL*DLP**2*(16.0D0+16.0D0*Y*Y+32.0D0*Y)
1953 B2=DL*DLP*(-32.0D0-176.0D0*Y*Y/3.0D0-64.0D0*Y-32.0D0/Y/3.0D0)+
1954 1DL*ZETA2*(-48.0D0)+
1955 2DL*S111MX*(-16.0D0-32.0D0*Y)+
1956 3DL*S11MX*(16.0D0-16.0D0*Y*Y+32.0D0*Y)+
1957 4DL*(944.0D0/3.0D0+8.0D0*Y*Y-208.0D0*Y/3.0D0)+
1958 5DLM**3*(-4.0D0/3.0D0+8.0D0*Y/3.0D0)+
1959 6DLM**2*(48.0D0-56.0D0*Y)+
1960 7ZETA2*DLM*(24.0D0+16.0D0*Y*Y-80.0D0*Y)
1961 B3=DLM*S111MX*(40.0D0+16.0D0*Y)+
1962 1DLM*S11MX*(-16.0D0-32.0D0*Y)+
1963 2DLM*(-176.0D0-8.0D0*Y*Y+212.0D0*Y)+
1964 3ZETA2*DLP*(16.0D0+16.0D0*Y*Y+32.0D0*Y)+
1965 4DLP*S11MX*(32.0D0+32.0D0*Y*Y+64.0D0*Y)+
1966 5S21MZ*(-16.0D0-32.0D0*Y)+
1967 6S21Z*(16.0D0+32.0D0*Y)
1968 B4=ZETA2*(-144.0D0-176.0D0*Y*Y/3.0D0+128.0D0*Y)+
1969 1ZETA3*(-12.0D0-16.0D0*Y*Y-40.0D0*Y)+
1970 2S111MX*(-24.0D0)+
1971 3S211MX*(-72.0D0-16.0D0*Y)+
1972 4S121MX*(8.0D0-16.0D0*Y*Y+16.0D0*Y)+
1973 5S11MX*(-32.0D0-176.0D0*Y*Y/3.0D0-64.0D0*Y-32.0D0/Y/3.0D0)+
1974 6S21MX*(-16.0D0+16.0D0*Y*Y-32.0D0*Y)+
1975 7S12MX*(32.0D0+32.0D0*Y*Y+64.0D0*Y)+
1976 81420.0D0/3.0D0-1468.0D0*Y/3.0D0
1977 B5=DL*(16.0D0+16.0D0*Y)+
1978 1DLM*(-8.0D0+16.0D0*Y)+
1979 248.0D0-48.0D0*Y
1980 B6=44.0D0/3.0D0-88.0D0*Y/3.0D0
1981 B7=DL**2*(-24.0D0-32.0D0*Y)+
1982 1DL*DLM*(48.0D0)+
1983 2DL*DLP*(-16.0D0-32.0D0*Y)+
1984 3DL*(-120.0D0+96.0D0*Y)+
1985 4DLM**2*(-8.0D0+16.0D0*Y)+
1986 5DLM*(112.0D0-128.0D0*Y)+
1987 6ZETA2*(-32.0D0*Y)+
1988 7S111MX*(32.0D0+32.0D0*Y)+
1989 8S11MX*(-16.0D0-32.0D0*Y)-
1990 9160.0D0+168.0D0*Y
1991 B8=44.0D0*(2.0D0*Y-1.0D0)*(DL-DLM)/3.0D0+
1992 1176.0D0*Y/3.0D0-44.0D0
1993 FCA=CA*(B1+B2+B3+B4+B5*DL2+B6*DL1*EL1+B7*DL1+B8*EL1)
1994 C1=16.0D0*Y/3.0D0-8.0D0/3.0D0
1995 C2=8.0D0*(2.0D0*Y-1.0D0)*(DLM-DL)/3.0D0-
1996 132.0D0*Y/3.0D0+8.0D0
1997 FNF=NF*(C1*DL1*EL1+C2*EL1)
1998C THE 0.5 BELOW IS DUE TO TF=1/2
1999 POLG2=0.5D0*NF*(FCF+FCA+FNF)
2000 RETURN
2001 END
2002
2003C ZEROTH ORDER ALPHAS CONTRIBUTION TO F3, F2 : V + Q -> Q + G
2004
2005
2006 DOUBLE PRECISION FUNCTION F30R(X)
2007 IMPLICIT REAL*8(A-H,O-Z)
2008 DOUBLE PRECISION NF
2009 COMMON/PAR/NF,Q2,FS,RS,DN
2010 COMMON/CHOICE/MODE
2011 F30R=1.0D0
2012 RETURN
2013 END
2014
2015
2016 DOUBLE PRECISION FUNCTION F20R(X)
2017 IMPLICIT REAL*8(A-H,O-Z)
2018 DOUBLE PRECISION NF
2019 COMMON/PAR/NF,Q2,FS,RS,DN
2020 COMMON/CHOICE/MODE
2021 F20R=1.0D0
2022 RETURN
2023 END
2024
2025
2026C ORDER ALPHAS CONTRIBUTION TO F3, F2 : V + Q -> Q + G
2027C REGULAR PART , MSBAR SCHEME .
2028
2029
2030 DOUBLE PRECISION FUNCTION F31R(X)
2031 IMPLICIT REAL*8(A-H,O-Z)
2032 DOUBLE PRECISION NF
2033 COMMON/PAR/NF,Q2,FS,RS,DN
2034 COMMON/CHOICE/MODE
2035 FA=X**(DN-1.d0)
2036 F31R=FA*C3NS1(X)
2037 RETURN
2038 END
2039
2040 DOUBLE PRECISION FUNCTION F21R(X)
2041 IMPLICIT REAL*8(A-H,O-Z)
2042 DOUBLE PRECISION NF
2043 COMMON/PAR/NF,Q2,FS,RS,DN
2044 COMMON/CHOICE/MODE
2045 FA=X**DN
2046 F21R=FA*C2NS1(X)
2047 RETURN
2048 END
2049
2050 DOUBLE PRECISION FUNCTION FG21R(X)
2051 IMPLICIT REAL*8(A-H,O-Z)
2052 DOUBLE PRECISION NF
2053 COMMON/PAR/NF,Q2,FS,RS,DN
2054 COMMON/CHOICE/MODE
2055 FA=X**DN
2056csm gluon coeff function
2057 FG21R=FA*C2G1(X)
2058 RETURN
2059 END
2060
2061 DOUBLE PRECISION FUNCTION FG22R(X)
2062 IMPLICIT REAL*8(A-H,O-Z)
2063 DOUBLE PRECISION NF
2064 COMMON/PAR/NF,Q2,FS,RS,DN
2065 COMMON/CHOICE/MODE
2066 FA=X**DN
2067csm gluon coeff function
2068 FG22R=FA*C2G2(X)
2069 RETURN
2070 END
2071
2072
2073 DOUBLE PRECISION FUNCTION FPS22R(X)
2074 IMPLICIT REAL*8(A-H,O-Z)
2075 DOUBLE PRECISION NF
2076 COMMON/PAR/NF,Q2,FS,RS,DN
2077 COMMON/CHOICE/MODE
2078 FA=X**DN
2079csm pure singlet coeff function
2080 FPS22R=FA*C2S(X)
2081 RETURN
2082 END
2083
2084 DOUBLE PRECISION FUNCTION FL1R(X)
2085 IMPLICIT REAL*8(A-H,O-Z)
2086 DOUBLE PRECISION NF
2087 COMMON/PAR/NF,Q2,FS,RS,DN
2088 COMMON/CHOICE/MODE
2089 FA=X**DN
2090 FL1R=FA*CLNS1(X)
2091 RETURN
2092 END
2093
2094 DOUBLE PRECISION FUNCTION FL2R(X)
2095 IMPLICIT REAL*8(A-H,O-Z)
2096 DOUBLE PRECISION NF
2097 COMMON/PAR/NF,Q2,FS,RS,DN
2098 COMMON/CHOICE/MODE
2099 FA=X**DN
2100 FL2R=FA*(CLNSP(X)+CLNSM(X))
2101 RETURN
2102 END
2103
2104 DOUBLE PRECISION FUNCTION FPSL2R(X)
2105 IMPLICIT REAL*8(A-H,O-Z)
2106 DOUBLE PRECISION NF
2107 COMMON/PAR/NF,Q2,FS,RS,DN
2108 COMMON/CHOICE/MODE
2109 FA=X**DN
2110csm pure singlet coeff function
2111 FPSL2R=FA*CLS(X)
2112 RETURN
2113 END
2114
2115
2116 DOUBLE PRECISION FUNCTION FLG1R(X)
2117 IMPLICIT REAL*8(A-H,O-Z)
2118 DOUBLE PRECISION NF
2119 COMMON/PAR/NF,Q2,FS,RS,DN
2120 COMMON/CHOICE/MODE
2121 FA=X**DN
2122 FLG1R=FA*CLG1(X)
2123 RETURN
2124 END
2125
2126
2127 DOUBLE PRECISION FUNCTION FLG2R(X)
2128 IMPLICIT REAL*8(A-H,O-Z)
2129 DOUBLE PRECISION NF
2130 COMMON/PAR/NF,Q2,FS,RS,DN
2131 COMMON/CHOICE/MODE
2132 FA=X**DN
2133 FLG2R=FA*CLG2(X)
2134 RETURN
2135 END
2136
2137
2138C ORDER ALPHAS CONTRIBUTION TO F3, F2 : V + Q -> Q + G
2139C SINGULAR PART , MSBAR SCHEME .
2140
2141 DOUBLE PRECISION FUNCTION F31S(X)
2142 IMPLICIT REAL*8(A-H,O-Z)
2143 DOUBLE PRECISION NF
2144 COMMON/PAR/NF,Q2,FS,RS,DN
2145 COMMON/CHOICE/MODE
2146 FA=X**(DN-1.d0)
2147 F31S=(FA-1.0D0)*SOFT1(X)
2148 RETURN
2149 END
2150
2151 DOUBLE PRECISION FUNCTION F21S(X)
2152 IMPLICIT REAL*8(A-H,O-Z)
2153 DOUBLE PRECISION NF
2154 COMMON/PAR/NF,Q2,FS,RS,DN
2155 COMMON/CHOICE/MODE
2156 FA=X**DN
2157 F21S=(FA-1.0D0)*SOFT1(X)
2158 RETURN
2159 END
2160
2161
2162
2163C ORDER ALPHAS**2 CORRECTION TO F3, F2 : V + Q -> Q + G + G ;
2164C V + Q -> Q + Q + QB (ID + NONID)
2165C REGULAR PART , MSBAR SCHEME .
2166
2167
2168 DOUBLE PRECISION FUNCTION F32R(X)
2169 IMPLICIT REAL*8(A-H,O-Z)
2170 DOUBLE PRECISION NF
2171 COMMON/PAR/NF,Q2,FS,RS,DN
2172 COMMON/CHOICE/MODE
2173 FA=X**(DN-1.d0)
2174 F32R=FA*(C3NSP(X)-C3NSM(X))
2175 RETURN
2176 END
2177
2178 DOUBLE PRECISION FUNCTION F22R(X)
2179 IMPLICIT REAL*8(A-H,O-Z)
2180 DOUBLE PRECISION NF
2181 COMMON/PAR/NF,Q2,FS,RS,DN
2182 COMMON/CHOICE/MODE
2183 FA=X**DN
2184C THIS IS THE COMBINATION FOR F2 IN CHARGED CURRENT INTERACTIONS
2185C F22R=FA*(C2NSP(X)-C2NSM(X))
2186C THIS IS THE COMBINATION FOR F2 IN NEUTRAL CURRENT INTERACTIONS
2187 F22R=FA*(C2NSP(X)+C2NSM(X))
2188 RETURN
2189 END
2190
2191
2192C ORDER ALPHAS**2 CORRECTION TO F3, F2 : V + Q -> Q + G
2193C SINGULAR PART , MSBAR SCHEME .
2194
2195 DOUBLE PRECISION FUNCTION F32S(X)
2196 IMPLICIT REAL*8(A-H,O-Z)
2197 DOUBLE PRECISION NF
2198 COMMON/PAR/NF,Q2,FS,RS,DN
2199 COMMON/CHOICE/MODE
2200 FA=X**(DN-1.d0)
2201 F32S=(FA-1.0D0)*SOFT2(X)
2202 RETURN
2203 END
2204
2205
2206 DOUBLE PRECISION FUNCTION F22S(X)
2207 IMPLICIT REAL*8(A-H,O-Z)
2208 DOUBLE PRECISION NF
2209 COMMON/PAR/NF,Q2,FS,RS,DN
2210 COMMON/CHOICE/MODE
2211 FA=X**DN
2212 F22S=(FA-1.0D0)*SOFT2(X)
2213 RETURN
2214 END
2215
2216
2217C ORDER ALPHAS CONTRIBUTION TO G1 : V + G -> Q + QB
2218C MSBAR SCHEME .
2219
2220
2221 DOUBLE PRECISION FUNCTION F1G1(X)
2222 IMPLICIT REAL*8(A-H,O-Z)
2223 DOUBLE PRECISION NF
2224 COMMON/PAR/NF,Q2,FS,RS,DN
2225 COMMON/CHOICE/MODE
2226 FA=X**DN
2227 F1G1=FA*POLG1(X)
2228 RETURN
2229 END
2230
2231
2232C ORDER ALPHAS**2 CORRECTION TO G1, : V + G -> G + Q + QB ;
2233C MSBAR SCHEME .
2234
2235
2236 DOUBLE PRECISION FUNCTION F1G2(X)
2237 IMPLICIT REAL*8(A-H,O-Z)
2238 DOUBLE PRECISION NF
2239 COMMON/PAR/NF,Q2,FS,RS,DN
2240 COMMON/CHOICE/MODE
2241 FA=X**DN
2242 F1G2=FA*POLG2(X)
2243 RETURN
2244 END
2245
2246
2247C ORDER ALPHAS**2 CORRECTION TO G1, : V + Q -> Q + Q + QB ;
2248C MSBAR SCHEME .
2249
2250 DOUBLE PRECISION FUNCTION F1Q2(X)
2251 IMPLICIT REAL*8(A-H,O-Z)
2252 DOUBLE PRECISION NF
2253 COMMON/PAR/NF,Q2,FS,RS,DN
2254 COMMON/CHOICE/MODE
2255 FA=X**DN
2256 F1Q2=FA*POLS(X)
2257 RETURN
2258 END
2259
2260
2261
2262
2263C HERE FOLLOWS THE TWO LOOP CORRECTED RUNNING COUPLING CONSTANT
2264
2265 DOUBLE PRECISION FUNCTION ALFAS1(RS,LAMBDA)
2266 IMPLICIT REAL*8(A-H,O-Z)
2267 DOUBLE PRECISION LAMBDA,NF
2268 COMMON/PAR/NF,Q2,FSS,RSS,DN
2269 ALFAS1=ALF(RS,LAMBDA,NF)
2270 RETURN
2271 END
2272
2273
2274 DOUBLE PRECISION FUNCTION ALFAS2(RS,LAMBDA)
2275 IMPLICIT REAL*8(A-H,O-Z)
2276 DOUBLE PRECISION LAMBDA,MB,MC
2277 MB=4.75D0
2278 MC=1.5D0
2279 A=1.0D0
2280 AMB=A*MB
2281 AMC=A*MC
2282 IF (RS.LT.AMB) GO TO 1
2283 ALF1=ALF(RS,LAMBDA,5.0D0)
2284 GO TO 3
2285 1 IF (RS.LT.AMC) GO TO 2
2286 ALF1I=1.0D0/ALF(RS,LAMBDA,4.0D0)
2287 ALF2I=1.0D0/ALF(AMB,LAMBDA,5.0D0)
2288 ALF3I=1.0D0/ALF(AMB,LAMBDA,4.0D0)
2289 ALF1=1.0D0/(ALF1I+ALF2I-ALF3I)
2290 GO TO 3
2291 2 ALF1I=1.0D0/ALF(RS,LAMBDA,3.0D0)
2292 ALF2I=1.0D0/ALF(AMC,LAMBDA,4.0D0)
2293 ALF3I=1.0D0/ALF(AMB,LAMBDA,5.0D0)
2294 ALF4I=1.0D0/ALF(AMB,LAMBDA,4.0D0)
2295 ALF5I=1.0D0/ALF(AMC,LAMBDA,3.0D0)
2296 ALF1=1.0D0/(ALF1I+ALF2I+ALF3I-ALF4I-ALF5I)
2297 3 ALFAS2=ALF1
2298 RETURN
2299 END
2300
2301
2302 DOUBLE PRECISION FUNCTION ALF(RS,LAMBDA,NF)
2303 IMPLICIT REAL*8(A-H,O-Z)
2304 DOUBLE PRECISION LAMBDA,LAMBD2,NF
2305 PI=3.14159265359D0
2306 LAMBD2=LAMBDA*LAMBDA
2307 BF=(33.0D0-2.0D0*NF)/12.0D0/PI
2308 BFP=(153.0D0-19.0D0*NF)/(66.0D0-4.0D0*NF)/PI
2309 RS2=RS*RS
2310 DL1=DLOG(RS2/LAMBD2)
2311 DL2=DLOG(DL1)
2312 ALF=(1.0D0-(BFP*DL2)/(BF*DL1))/DL1/BF
2313C ALF=1.0D0/DL1/BF
2314 RETURN
2315 END
2316 SUBROUTINE STRUCT(Z,FS,MODE,UV,DV,SS,ST,CH,BOT,GL)
2317 IMPLICIT REAL*8 (A-H,O-Z)
2318 UV=8*Z*(1-Z)**3
2319 DV=5*Z*(1-Z)**4
2320 GL=21*Z*(1-Z)**5
2321 SS=9*Z*(1-Z)**8
2322 ST=SS
2323 CH=0.0D0
2324 BOT=0.0D0
2325 RETURN
2326 END
2327 INCLUDE 'daind.f'
2328 INCLUDE 'daind1.f'
2329 INCLUDE 'wgplg.f'