Corrected path to mapping libraries.
[u/mrichter/AliRoot.git] / PDF / tpdf / testpdgo.F
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
14 C
15       PARAMETER (LFLPRT = 2)
16 C+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"
25 C
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)
34 C
35       DATA NXX,XX/3, .002,.140,.700/
36       DATA NQ2,Q2/3, 1.30,24.0,33.5/
37 C
38       N6 = L6
39       IFLPRT = LFLPRT
40 C
41       CH1=4./9.
42       CH2=1./9.
43 C
44       SCAL=SQRT(Q2(2))
45 C
46       DO 300 ISTRUC=-3,MODEMX
47       FIRST=.TRUE.
48       IF(ISTRUC.LE.-2) IFLSET = 0
49 C... 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
73 C
74 C... 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
89 C
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)
95 C
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)
110 C
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
140 C
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)
146 C
147  200  CONTINUE
148 C
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
155 C
156  300  CONTINUE
157 C
158       STOP
159       END