5 * Revision 1.1.1.2 1996/10/30 08:32:51 cernlib
8 * Revision 1.1.1.1 1996/04/12 15:30:17 plothow
12 #include "pdf/pilot.h"
15 #include "pdf/impdp.inc"
16 PARAMETER (LFLPRT = 2)
18 PARAMETER (MODEMX = 281)
19 #include "pdf/w5051p7.inc"
20 #include "pdf/w50510.inc"
21 #include "pdf/w50511.inc"
22 #include "pdf/w50512.inc"
23 #include "pdf/w50513.inc"
24 #include "pdf/w50516.inc"
25 #include "pdf/w50517.inc"
27 #include "pdf/expdp.inc"
28 + DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DTB,DGL
31 #include "pdf/expdp.inc"
32 + XI(32),WI(32),XX(33)
33 COMMON/W5051I9/XI,WI,XX,NTERMS
34 #include "pdf/expdp.inc"
36 DATA XV/0.05D0,0.50D0,0.95D0/
37 DATA ZQ/10.D0,100.D0,1000.D0,1.D+4,5.5D+5/
42 #include "pdf/expdp.inc"
48 IF(GLUCK.LT.0.01) GLUCK = 0.
55 B1=102.D0-38.D0*NF/3.D0
57 IF(NF.EQ.4) CCOEG=5.D0/18.D0
58 COEG=2.D0*FLAVOR*CCOEG
66 IF(ISTRUC.GE.201 .AND. ISTRUC.LE.205) GOTO 300
67 IF(ISTRUC.GE.208 .AND. ISTRUC.LE.209) GOTO 300
68 IF(ISTRUC.GE.214 .AND. ISTRUC.LE.219) GOTO 300
69 IF(ISTRUC.GE.223 .AND. ISTRUC.LE.229) GOTO 300
70 IF(ISTRUC.GE.233 .AND. ISTRUC.LE.270) GOTO 300
71 IF(ISTRUC.GE.274 .AND. ISTRUC.LE.279) GOTO 300
72 IF(ISTRUC.GE.282 ) GOTO 300
84 ALFPI= 2.D0 /(B0*ALQ2+B1* LOG(ALQ2)/B0)
90 CALL STRUCTM(DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DTB,DGL)
96 F2=(4.D0* DUB + DDB + DSB)/9.D0*2.D0
97 F2C=F2+(4.D0/9.D0)*DCB*2.D0
100 WRITE (N6,1000) Q02,IORD,OWLAM,GLUCK,FLAV
101 1000 FORMAT(1X,' Q02= ',F7.3,' IORD= ',I2,' LAMBDA=',F7.4,
102 + ' GLUCK= ',F5.2,' NF = ',F5.2)
106 1001 FORMAT(/,' Q2= ',E10.3,' GEV**2')
108 1002 FORMAT(/,4X,' X',6X,'X*UB ',6X,'X*DB ',6X,'X*SB ',6X,'X*CB '
109 + ,6X,'X*GL',8X,' F2 ',6X,' F2/ALPHA'
110 + ,5X,' F2C ',5X,' F2C/ALPHA')
112 IF (GLUCK.LT..5) GOTO 5
118 XQ=F2+F2*AL*CF*(-9.D0-2.D0*PI2/3.D0+AL1*(-3.D0+2.D0*AL1))
119 C WRITE(N6,1111) X,Q,F2,AL1,XQ
120 C1111 FORMAT(1X,' X,F2 = ',5E13.5)
122 Y=0.5D0*(1.D0-X)*XI(I)+0.5D0*(1.D0+X)
126 CALL STRUCTM(DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DTB,DGL)
127 XQQ=(4.D0* DUB + DDB + DSB)/9.D0*2.D0
128 IF (GLUCK.LT..5) GOTO 6
130 6 C22=CF*(6.D0+4.D0*Y-2.D0*(1.D0+Y*Y)/
131 + (1.D0-Y)* LOG(Y)-2.D0*(1.D0+Y)*AL1)
132 C23=CF*(-3.D0+4.D0*AL1)/(1.D0-Y)
133 CG2=(-1.D0+8.D0*Y*(1.D0-Y)+(1.D0-2.D0*Y+2.D0*Y*Y)*
134 + LOG(1.D0/Y-1.D0))*COEG
135 XQ=XQ +0.5D0*(1.D0-X)*WI(I)*AL*(C22*XQQ+C23*(XQQ- F2))
136 XQ=XQ +0.5D0*(1.D0-X)*WI(I)*AL*CG2*DGL
137 C WRITE(N6,1112) Y,AL1,C22,C23,CG2,XQQ,XQ
138 C1112 FORMAT(1X,' Y,C22 = ',7E13.5)
140 F2=3.D0*3.D0*(2.D0/27.D0)*(1.D0/(2.D0*PI*ALQED))*2.D0
141 F2=XQ+F2*(8.D0*X*(1.D0-X)-1.D0+(1.D0-2.D0*X+2.D0*X*X)*
142 + LOG(1.D0/X-1.D0))*X
143 F2C=F2+(4.D0/9.D0)*CB*2.D0
146 WRITE(N6,1003)X,UB,DB,SB,CB,GL,F2,F2N,F2C,F2CN
147 1003 FORMAT(1X,' ',F6.3,9E12.4)