Corrected path to mapping libraries.
[u/mrichter/AliRoot.git] / PDF / tpdf / testpdgo.F
CommitLineData
21886bb6 1*
2* $Id$
3*
4* $Log$
5* Revision 1.1.1.2 1996/10/30 08:32:56 cernlib
6* Version 7.04
7*
8* Revision 1.1.1.1 1996/04/12 15:30:18 plothow
9* Version 7.01
10*
11*
12#include "pdf/pilot.h"
13 PROGRAM TESTPDGO
14C
15 PARAMETER (LFLPRT = 2)
16C+SEQ, W5051P2.
17 PARAMETER (MODEMX = 281)
18#include "pdf/w5051p7.inc"
19#include "pdf/w50510.inc"
20#include "pdf/w50511.inc"
21#include "pdf/w50512.inc"
22#include "pdf/w50515.inc"
23#include "pdf/w50516.inc"
24#include "pdf/w50517.inc"
25C
26#include "pdf/expdp.inc"
27 + DX,DQ,DXPDF(-6:6),
28 + ALF,ALPHAS2
29 REAL X, Q, UPV, DNV, USEA, DSEA, SSEA, CSEA, BOT, TOP, GLUE
30 DIMENSION XX(3),Q2(3)
31 CHARACTER*20 PARM(20)
32#include "pdf/expdp.inc"
33 + VAL(20)
34C
35 DATA NXX,XX/3, .002,.140,.700/
36 DATA NQ2,Q2/3, 1.30,24.0,33.5/
37C
38 N6 = L6
39 IFLPRT = LFLPRT
40C
41 CH1=4./9.
42 CH2=1./9.
43C
44 SCAL=SQRT(Q2(2))
45C
46 DO 300 ISTRUC=-3,MODEMX
47 FIRST=.TRUE.
48 IF(ISTRUC.LE.-2) IFLSET = 0
49C... test different MODE writings
50 IF(ISTRUC .LT. 0) GOTO 310
51 PARM(1) = 'MODE'
52 IF(ISTRUC.GT.20) PARM(1) = 'Mode'
53 IF(ISTRUC.GT.30) PARM(1) = 'mode'
54 IF(ISTRUC.GT.40) PARM(1) = 'MOde'
55 IF(ISTRUC.GT.50) PARM(1) = 'moDE'
56 VAL(1) = ISTRUC
57 IF(ISTRUC.EQ. 19 ) GOTO 300
58 IF(ISTRUC.GE. 48 .AND. ISTRUC.LE. 49) GOTO 300
59 IF(ISTRUC.GE. 61 .AND. ISTRUC.LE. 69) GOTO 300
60 IF(ISTRUC.EQ. 74 ) GOTO 300
61 IF(ISTRUC.GE. 76 .AND. ISTRUC.LE. 79) GOTO 300
62 IF(ISTRUC.GE. 81 .AND. ISTRUC.LE.105) GOTO 300
63 IF(ISTRUC.GE.108 .AND. ISTRUC.LE.129) GOTO 300
64 IF(ISTRUC.GE.133 .AND. ISTRUC.LE.171) GOTO 300
65 IF(ISTRUC.GE.174 .AND. ISTRUC.LE.179) GOTO 300
66 IF(ISTRUC.GE.183 .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
74C... test different MODE writings and INTEGER settings
75 310 CONTINUE
76 IF(ISTRUC.GE.0 .AND. ISTRUC.LE.MODEMX) THEN
77 CALL PDFSET(PARM,VAL)
78 ELSEIF (ISTRUC.EQ.-3) THEN
79 GOTO 5
80 ELSEIF (ISTRUC.EQ.-2) THEN
81 PARM(1) = 'MODE'
82 VAL(1) = 45.D0
83 CALL PDFSET(PARM,VAL)
84 ELSEIF (ISTRUC.EQ.-1) THEN
85 PARM(1) = 'mode'
86 VAL(1) = 45.D0
87 CALL PDFSET(PARM,VAL)
88 ENDIF
89C
90 DQ = SCAL
91 ALF = ALPHAS2(DQ)
92 WRITE (N6,2000) MODE,QCDL4,QCDL5,DQ,ALF
93 2000 FORMAT(/,' Mode = ',I3,' QCDL4, QCDL5 = ',2F7.4,
94 + ' SCALE = ',F7.4,' alpha(s) = ',F8.5)
95C
96 5 CONTINUE
97 DO 200 IX=1,NXX
98 X=XX(IX)
99 Q=SQRT(Q2(IX))
100 DX=X
101 DQ=Q
102 CALL PFTOPDG(DX,DQ,DXPDF)
103 DNV=DXPDF(1)
104 UPV=DXPDF(2)
105 DSEA=DXPDF(-1)
106 USEA=DXPDF(-2)
107 SSEA=DXPDF(3)
108 CSEA=DXPDF(4)
109 GLUE=DXPDF(0)
110C
111 FFPU = UPV
112 FFPUB = USEA
113 FFPD = DNV
114 FFPDB = DSEA
115 FFPS = SSEA
116 FFPSB = SSEA
117 FFPC = CSEA
118 FFPCB = CSEA
119 FFNU = FFPD
120 FFNUB = FFPDB
121 FFND = FFPU
122 FFNDB = FFPUB
123 FFNS = FFPS
124 FFNSB = FFPSB
125 FFNC = FFPC
126 FFNCB = FFPCB
127 FFPGL = GLUE
128 SFFP=CH1*(FFPU+FFPUB+FFPC+FFPCB)+CH2*(FFPD+FFPDB+FFPS+FFPSB)
129 SFFN=CH1*(FFNU+FFNUB+FFNC+FFNCB)+CH2*(FFND+FFNDB+FFNS+FFNSB)
130 FACNP = 0.
131 IF (SFFP .NE. 0.) FACNP=SFFN/SFFP
132 UP = FFPU
133 UV = FFPU - FFPUB
134 US = FFPUB
135 DP = FFPD
136 DV = FFPD - FFPDB
137 DS = FFPDB
138 GL = FFPGL
139 F2P = SFFP
140C
141 WRITE(N6,1000)X,Q2(IX),FACNP,UP,UV,US,DP,DV,GL,F2P
142 1000 FORMAT(1X,' X:',F6.3,2X,'Q2:',F5.1,
143 1 3X,'N/P:',F7.4,2X,'UP:',F6.3,
144 1 1X,'UV:',F6.3,1X,'US:',F6.3,1X,'DP:',F6.3,1X,'DV:',F6.3,
145 2 1X,'GL:',F6.3,2X,'F2P:',F7.4)
146C
147 200 CONTINUE
148C
149 IF(ISTRUC.EQ.-3) THEN
150 DQ = SCAL
151 ALF = ALPHAS2(DQ)
152 WRITE (N6,2000) MODE,QCDL4,QCDL5,DQ,ALF
153 ENDIF
154 CALL PDFSTA
155C
156 300 CONTINUE
157C
158 STOP
159 END