Corrected path to mapping libraries.
[u/mrichter/AliRoot.git] / PDF / tpdf / testgamo.F
CommitLineData
21886bb6 1*
2* $Id$
3*
4* $Log$
5* Revision 1.1.1.2 1996/10/30 08:32:51 cernlib
6* Version 7.04
7*
8* Revision 1.1.1.1 1996/04/12 15:30:17 plothow
9* Version 7.01
10*
11*
12#include "pdf/pilot.h"
13 PROGRAM TESTGAMO
14C
15#include "pdf/impdp.inc"
16 PARAMETER (LFLPRT = 2)
17C+SEQ, W5051P2.
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"
26C
27#include "pdf/expdp.inc"
28 + DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DTB,DGL
29 REAL GLUCK
30 COMMON/W5051I8/GLUCK
31#include "pdf/expdp.inc"
32 + XI(32),WI(32),XX(33)
33 COMMON/W5051I9/XI,WI,XX,NTERMS
34#include "pdf/expdp.inc"
35 + XV(3),ZQ(5)
36 DATA XV/0.05D0,0.50D0,0.95D0/
37 DATA ZQ/10.D0,100.D0,1000.D0,1.D+4,5.5D+5/
38 DATA FLAV/4.D0/
39 DATA ALQED/137.D0/
40 DATA ICALL/0/
41 CHARACTER*20 PARM(20)
42#include "pdf/expdp.inc"
43 + VAL(20)
44C
45 N6 = L6
46 IFLPRT = LFLPRT
47C
48 IF(GLUCK.LT.0.01) GLUCK = 0.
49 PI=4.D0*ATAN(1.D0)
50 PI2=PI**2
51 CF=4.D0/3.D0
52 NF=INT(FLAV+0.1D0)
53 FLAVOR=FLAV
54 B0=11.D0-2.D0*NF/3.D0
55 B1=102.D0-38.D0*NF/3.D0
56 CCOEG=2.D0/9.D0
57 IF(NF.EQ.4) CCOEG=5.D0/18.D0
58 COEG=2.D0*FLAVOR*CCOEG
59 CALL WATE32
60C
61 DO 300 ISTRUC=272,272
62 ICALL = 0
63 FIRST=.TRUE.
64 VAL(1) = ISTRUC
65 PARM(1) = 'Mode'
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
73C
74 CALL PDFSET(PARM,VAL)
75 Q02 = Q2MIN
76 OWLAM = QCDL4
77 OWLAM2=OWLAM**2
78 IORD = LO - 1
79*
80 DO 100 J=3,3
81 Q2=Q02
82 IF(J.GE.2) Q2=ZQ(J-1)
83 ALQ2= LOG(Q2/OWLAM2)
84 ALFPI= 2.D0 /(B0*ALQ2+B1* LOG(ALQ2)/B0)
85 DO 70 K=1,3
86 X=XV(K)
87 DX = X
88 Q = SQRT(Q2)
89 DQ = Q
90 CALL STRUCTM(DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DTB,DGL)
91 UB =DUB
92 DB =DDB
93 SB =DSB
94 CB =DCB
95 GL =DGL
96 F2=(4.D0* DUB + DDB + DSB)/9.D0*2.D0
97 F2C=F2+(4.D0/9.D0)*DCB*2.D0
98 IF(ICALL.EQ.0) THEN
99 ICALL = 1
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)
103 ENDIF
104 IF (K.EQ.1) THEN
105 WRITE (N6,1001) Q2
106 1001 FORMAT(/,' Q2= ',E10.3,' GEV**2')
107 WRITE(N6,1002)
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')
111 ENDIF
112 IF (GLUCK.LT..5) GOTO 5
113 F2=DUV
114 5 IF(IORD)13,13,14
11514 CONTINUE
116 AL=0.5D0*ALFPI
117 AL1= LOG(1.D0-X)
118 XQ=F2+F2*AL*CF*(-9.D0-2.D0*PI2/3.D0+AL1*(-3.D0+2.D0*AL1))
119C WRITE(N6,1111) X,Q,F2,AL1,XQ
120C1111 FORMAT(1X,' X,F2 = ',5E13.5)
121 DO 81 I=1,NTERMS
122 Y=0.5D0*(1.D0-X)*XI(I)+0.5D0*(1.D0+X)
123 XY=X/Y
124 AL1= LOG(1.D0-Y)
125 DX = XY
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
129 XQQ=DUV
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
137C WRITE(N6,1112) Y,AL1,C22,C23,CG2,XQQ,XQ
138C1112 FORMAT(1X,' Y,C22 = ',7E13.5)
139 81 CONTINUE
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
144 13 F2N=F2*ALQED
145 F2CN=F2C*ALQED
146 WRITE(N6,1003)X,UB,DB,SB,CB,GL,F2,F2N,F2C,F2CN
147 1003 FORMAT(1X,' ',F6.3,9E12.4)
148 70 CONTINUE
149 100 CONTINUE
150 WRITE(N6,*) ' '
151 300 CONTINUE
152 STOP
153 END