]>
Commit | Line | Data |
---|---|---|
1 | * | |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.1.1.4 1997/07/13 13:25:32 cernlib | |
6 | * Import version 7.09 | |
7 | * | |
8 | * Revision 1.4 1997/07/02 15:03:12 pdflib | |
9 | * potential bug removed in GRSGALO | |
10 | * | |
11 | * Revision 1.3 1997/04/08 15:53:28 pdflib | |
12 | * Lambda_QCD for CTEQ4 sets corrected | |
13 | * | |
14 | * Revision 1.2 1996/10/28 16:10:15 pdflib | |
15 | * PDFLIB new version 7.04 | |
16 | * | |
17 | * Revision 1.1.1.1 1996/04/12 15:30:17 plothow | |
18 | * Version 7.01 | |
19 | * | |
20 | * | |
21 | #include "pdf/pilot.h" | |
22 | PROGRAM TESTPDF | |
23 | C | |
24 | PARAMETER (LFLPRT = 2) | |
25 | #include "pdf/w5051p2.inc" | |
26 | #include "pdf/w5051p7.inc" | |
27 | #include "pdf/w50510.inc" | |
28 | #include "pdf/w50511.inc" | |
29 | #include "pdf/w50512.inc" | |
30 | #include "pdf/w50515.inc" | |
31 | #include "pdf/w50516.inc" | |
32 | #include "pdf/w50517.inc" | |
33 | #include "pdf/w505110.inc" | |
34 | #include "pdf/w505120.inc" | |
35 | #include "pdf/w505122.inc" | |
36 | C | |
37 | #include "pdf/expdp.inc" | |
38 | + DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL, | |
39 | + DQ2,DP2,ALF,ALPHAS2 | |
40 | REAL X, Q, UPV, DNV, USEA, DSEA, SSEA, CSEA, BOT, TOP, GLUE | |
41 | DIMENSION XX(3),Q2(3) | |
42 | CHARACTER*20 PARM(NCHDIM) | |
43 | #include "pdf/expdp.inc" | |
44 | + VAL(NCHDIM) | |
45 | C | |
46 | DATA NXX,XX/3, .0005,.140,.700/ | |
47 | DATA NQ2,Q2/3, 1.30,24.0,33.5/ | |
48 | DATA P2/5./ | |
49 | DATA IP2/0/ | |
50 | C | |
51 | N6 = L6 | |
52 | IFLPRT = LFLPRT | |
53 | C | |
54 | CH1=4./9. | |
55 | CH2=1./9. | |
56 | C | |
57 | SCAL=SQRT(Q2(2)) | |
58 | C | |
59 | C First call to PDFSET to initialize COMMON/W505120/ | |
60 | PARM(1) = 'Init0' | |
61 | VAL(1) = 0.D0 | |
62 | CALL PDFSET(PARM,VAL) | |
63 | C. | |
64 | C Test new PDFLIB versions (4.0 or bigger) | |
65 | DO KPTY = 1,NPTYMX | |
66 | DO KGRO = 1,NGRMAX | |
67 | IF(NPGSMX(KPTY,KGRO) .EQ. 0) GOTO 400 | |
68 | DO ISET=-9,NPGSMX(KPTY,KGRO) | |
69 | IF((KPTY*KGRO).NE.1 .AND. ISET.LE.0) GOTO 300 | |
70 | WRITE (N6,*) ' ========================================== | |
71 | +=================================================================' | |
72 | FIRST=.TRUE. | |
73 | IF(ISET.LE.-8) IFLSET = 0 | |
74 | IF(ISET.LT.0) GOTO 50 | |
75 | PARM(1) = 'Nptype' | |
76 | VAL(1) = KPTY | |
77 | PARM(2) = 'Ngroup' | |
78 | VAL(2) = KGRO | |
79 | PARM(3) = 'Nset' | |
80 | VAL(3) = ISET | |
81 | C... test different NSET writings | |
82 | IF(ISET.GT. 5) PARM(3) = 'NSET' | |
83 | IF(ISET.GT.10) PARM(3) = 'nset' | |
84 | IF(ISET.GT.15) PARM(3) = 'NSet' | |
85 | IF(ISET.GT.20) PARM(3) = 'nsET' | |
86 | C | |
87 | C... test different NSET writings and INVALID parameter settings | |
88 | 50 CONTINUE | |
89 | IF(ISET.GE.0) THEN | |
90 | CALL PDFSET(PARM,VAL) | |
91 | ELSEIF (ISET.EQ.-9) THEN | |
92 | GOTO 100 | |
93 | ELSEIF (ISET.EQ.-8) THEN | |
94 | PARM(1) = 'dflm' | |
95 | VAL(1) = 8.D0 | |
96 | CALL PDFSET(PARM,VAL) | |
97 | ELSEIF (ISET.EQ.-7) THEN | |
98 | PARM(1) = 'GRV' | |
99 | VAL(1) = 3.D0 | |
100 | CALL PDFSET(PARM,VAL) | |
101 | ELSEIF (ISET.EQ.-6) THEN | |
102 | PARM(1) = 'abfkw' | |
103 | VAL(1) = 1.D0 | |
104 | CALL PDFSET(PARM,VAL) | |
105 | ELSEIF (ISET.EQ.-5) THEN | |
106 | PARM(1) = 'GRVph' | |
107 | VAL(1) = 2.D0 | |
108 | CALL PDFSET(PARM,VAL) | |
109 | ELSEIF (ISET.EQ.-4) THEN | |
110 | PARM(1) = 'Nptype' | |
111 | VAL(1) = -1 | |
112 | PARM(2) = 'Ngroup' | |
113 | VAL(2) = KGRO | |
114 | PARM(3) = 'Nset' | |
115 | VAL(3) = ISET | |
116 | CALL PDFSET(PARM,VAL) | |
117 | ELSEIF (ISET.EQ.-3) THEN | |
118 | PARM(1) = 'Nptype' | |
119 | VAL(1) = -1 | |
120 | PARM(2) = 'Ngroup' | |
121 | VAL(2) = KGRO | |
122 | PARM(3) = 'Nset' | |
123 | VAL(3) = 0 | |
124 | PARM(4) = 'QCDL4' | |
125 | VAL(4) = 0.1900D0 | |
126 | PARM(5) = 'QCDL5' | |
127 | VAL(5) = 0.1230D0 | |
128 | CALL PDFSET(PARM,VAL) | |
129 | ELSEIF (ISET.EQ.-2) THEN | |
130 | PARM(1) = 'Nptype' | |
131 | VAL(1) = KPTY | |
132 | PARM(2) = 'Ngroup' | |
133 | VAL(2) = -1 | |
134 | PARM(3) = 'Nset' | |
135 | VAL(3) = ISET | |
136 | CALL PDFSET(PARM,VAL) | |
137 | ELSEIF (ISET.EQ.-1) THEN | |
138 | PARM(1) = 'Nptype' | |
139 | VAL(1) = KPTY | |
140 | PARM(2) = 'Ngroup' | |
141 | VAL(2) = KGRO | |
142 | PARM(3) = 'Nset' | |
143 | VAL(3) = 80 | |
144 | CALL PDFSET(PARM,VAL) | |
145 | ENDIF | |
146 | C | |
147 | DQ = SCAL | |
148 | ALF = ALPHAS2(DQ) | |
149 | WRITE(N6,4000) NPTYPE,NGROUP,NSET, | |
150 | + SFNAME(NPTYPE,NGROUP,NSET),MODECR(NPTYPE,NGROUP,NSET) | |
151 | WRITE(N6,4001) QCDL4,QCDL5,DQ,ALF | |
152 | 4000 FORMAT(/,' Nptype = ',I1,' Ngroup = ',I1, | |
153 | + ' Nset = ',I2,' Name = "',A8, | |
154 | + '" CrMode = ',I3) | |
155 | 4001 FORMAT(1H ,' QCDL4, QCDL5 = ',2F7.4, | |
156 | + ' Scale = ',F8.4,' alpha(s) = ',F7.5) | |
157 | C | |
158 | 100 CONTINUE | |
159 | DO 200 IX=1,NXX | |
160 | X=XX(IX) | |
161 | Q=SQRT(Q2(IX)) | |
162 | DX=X | |
163 | DQ=Q | |
164 | IF(KPTY.EQ.3) THEN | |
165 | DQ2=Q2(IX) | |
166 | DP2=P2 | |
167 | IP=IP2 | |
168 | CALL STRUCTP(DX,DQ2,DP2,IP, | |
169 | + DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL) | |
170 | ELSE | |
171 | CALL STRUCTM(DX,DQ, | |
172 | + DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL) | |
173 | ENDIF | |
174 | UPV=DUPV | |
175 | DNV=DDNV | |
176 | USEA=DUSEA | |
177 | DSEA=DDSEA | |
178 | SSEA=DSTR | |
179 | CSEA=DCHM | |
180 | GLUE=DGL | |
181 | C | |
182 | FFPU = UPV + USEA | |
183 | FFPUB = USEA | |
184 | FFPD = DNV + DSEA | |
185 | FFPDB = DSEA | |
186 | FFPS = SSEA | |
187 | FFPSB = SSEA | |
188 | FFPC = CSEA | |
189 | FFPCB = CSEA | |
190 | FFNU = FFPD | |
191 | FFNUB = FFPDB | |
192 | FFND = FFPU | |
193 | FFNDB = FFPUB | |
194 | FFNS = FFPS | |
195 | FFNSB = FFPSB | |
196 | FFNC = FFPC | |
197 | FFNCB = FFPCB | |
198 | FFPGL = GLUE | |
199 | SFFP=CH1*(FFPU+FFPUB+FFPC+FFPCB)+CH2*(FFPD+FFPDB+FFPS+FFPSB) | |
200 | SFFN=CH1*(FFNU+FFNUB+FFNC+FFNCB)+CH2*(FFND+FFNDB+FFNS+FFNSB) | |
201 | FACNP = 0. | |
202 | IF (SFFP .NE. 0.) FACNP=SFFN/SFFP | |
203 | UP = FFPU | |
204 | UV = FFPU - FFPUB | |
205 | US = FFPUB | |
206 | DP = FFPD | |
207 | DV = FFPD - FFPDB | |
208 | DS = FFPDB | |
209 | GL = FFPGL | |
210 | F2P = SFFP | |
211 | C | |
212 | WRITE(N6,3000)X,Q2(IX),FACNP,UP,UV,US,DP,DV,GL,F2P | |
213 | 3000 FORMAT(1X,' X:',F6.3,2X,'Q2:',F5.1, | |
214 | 1 3X,'N/P:',F7.4,2X,'UP:',F6.3, | |
215 | 1 1X,'UV:',F6.3,1X,'US:',F6.3,1X,'DP:',F6.3,1X,'DV:',F6.3, | |
216 | 2 1X,'GL:',F6.3,2X,'F2P:',F7.4) | |
217 | C | |
218 | 200 CONTINUE | |
219 | C | |
220 | IF(ISET.EQ.-9) THEN | |
221 | DQ = SCAL | |
222 | ALF = ALPHAS2(DQ) | |
223 | WRITE(N6,4000) NPTYPE,NGROUP,NSET, | |
224 | + SFNAME(NPTYPE,NGROUP,NSET),MODECR(NPTYPE,NGROUP,NSET) | |
225 | WRITE(N6,4001) QCDL4,QCDL5,DQ,ALF | |
226 | ENDIF | |
227 | CALL PDFSTA | |
228 | C | |
229 | 300 CONTINUE | |
230 | ENDDO | |
231 | 400 CONTINUE | |
232 | ENDDO | |
233 | ENDDO | |
234 | C | |
235 | STOP | |
236 | END |