5 * Revision 1.1.1.2 1996/10/30 08:32:54 cernlib
8 * Revision 1.1.1.1 1996/04/12 15:30:17 plothow
12 #include "pdf/pilot.h"
15 PARAMETER (LFLPRT = 2)
16 #include "pdf/w5051p2.inc"
17 #include "pdf/w5051p7.inc"
18 #include "pdf/w50510.inc"
19 #include "pdf/w50511.inc"
20 #include "pdf/w50512.inc"
21 #include "pdf/w50515.inc"
22 #include "pdf/w50516.inc"
23 #include "pdf/w50517.inc"
24 #include "pdf/w505110.inc"
25 #include "pdf/w505120.inc"
26 #include "pdf/w505122.inc"
28 #include "pdf/expdp.inc"
31 REAL X, Q, UPV, DNV, USEA, DSEA, SSEA, CSEA, BOT, TOP, GLUE
33 CHARACTER*20 PARM(NCHDIM)
34 #include "pdf/expdp.inc"
37 DATA NXX,XX/3, .002,.140,.700/
38 DATA NQ2,Q2/3, 1.30,24.0,33.5/
48 C First call to PDFSET to initialize COMMON/W505120/
53 C Test new PDFLIB versions (4.0 or bigger)
56 IF(NPGSMX(KPTY,KGRO) .EQ. 0) GOTO 400
57 DO ISET=-9,NPGSMX(KPTY,KGRO)
58 IF((KPTY*KGRO).NE.1 .AND. ISET.LE.0) GOTO 300
59 WRITE (N6,*) ' ==========================================
60 +================================================================='
62 IF(ISET.LE.-8) IFLSET = 0
70 C... test different NSET writings
71 IF(ISET.GT. 5) PARM(3) = 'NSET'
72 IF(ISET.GT.10) PARM(3) = 'nset'
73 IF(ISET.GT.15) PARM(3) = 'NSet'
74 IF(ISET.GT.20) PARM(3) = 'nsET'
76 C... test different NSET writings and INVALID parameter settings
80 ELSEIF (ISET.EQ.-9) THEN
82 ELSEIF (ISET.EQ.-8) THEN
86 ELSEIF (ISET.EQ.-7) THEN
90 ELSEIF (ISET.EQ.-6) THEN
94 ELSEIF (ISET.EQ.-5) THEN
98 ELSEIF (ISET.EQ.-4) THEN
105 CALL PDFSET(PARM,VAL)
106 ELSEIF (ISET.EQ.-3) THEN
113 CALL PDFSET(PARM,VAL)
114 ELSEIF (ISET.EQ.-2) THEN
121 CALL PDFSET(PARM,VAL)
122 ELSEIF (ISET.EQ.-1) THEN
129 CALL PDFSET(PARM,VAL)
134 WRITE(N6,4000) NPTYPE,NGROUP,NSET,
135 + SFNAME(NPTYPE,NGROUP,NSET),MODECR(NPTYPE,NGROUP,NSET)
136 WRITE(N6,4001) QCDL4,QCDL5,DQ,ALF
137 4000 FORMAT(/,' Nptype = ',I1,' Ngroup = ',I1,
138 + ' Nset = ',I2,' Name = "',A8,
140 4001 FORMAT(1H ,' QCDL4, QCDL5 = ',2F7.4,
141 + ' Scale = ',F8.4,' alpha(s) = ',F7.5)
149 CALL PFTOPDG(DX,DQ,DXPDF)
175 SFFP=CH1*(FFPU+FFPUB+FFPC+FFPCB)+CH2*(FFPD+FFPDB+FFPS+FFPSB)
176 SFFN=CH1*(FFNU+FFNUB+FFNC+FFNCB)+CH2*(FFND+FFNDB+FFNS+FFNSB)
178 IF (SFFP .NE. 0.) FACNP=SFFN/SFFP
188 WRITE(N6,3000)X,Q2(IX),FACNP,UP,UV,US,DP,DV,GL,F2P
189 3000 FORMAT(1X,' X:',F6.3,2X,'Q2:',F5.1,
190 1 3X,'N/P:',F7.4,2X,'UP:',F6.3,
191 1 1X,'UV:',F6.3,1X,'US:',F6.3,1X,'DP:',F6.3,1X,'DV:',F6.3,
192 2 1X,'GL:',F6.3,2X,'F2P:',F7.4)
199 WRITE(N6,4000) NPTYPE,NGROUP,NSET,
200 + SFNAME(NPTYPE,NGROUP,NSET),MODECR(NPTYPE,NGROUP,NSET)
201 WRITE(N6,4001) QCDL4,QCDL5,DQ,ALF