]>
Commit | Line | Data |
---|---|---|
21886bb6 | 1 | #include "pdf/pilot.h" |
2 | PROGRAM TEST_NPDF | |
3 | C | |
4 | PARAMETER (LFLPRT = 2) | |
5 | #include "pdf/w5051p2.inc" | |
6 | #include "pdf/w5051p7.inc" | |
7 | #include "pdf/w50510.inc" | |
8 | #include "pdf/w50511.inc" | |
9 | #include "pdf/w50511a.inc" | |
10 | #include "pdf/w50512.inc" | |
11 | #include "pdf/w50515.inc" | |
12 | #include "pdf/w50516.inc" | |
13 | #include "pdf/w50517.inc" | |
14 | #include "pdf/w505110.inc" | |
15 | #include "pdf/w505120.inc" | |
16 | #include "pdf/w505122.inc" | |
17 | C | |
18 | #include "pdf/expdp.inc" | |
19 | + DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL, | |
20 | + DUV,DDV,DUS,DDS,DST,DCH,DB,DT,DG, | |
21 | + RUV,RDV,RUB,RDB,RS,RC,RB,RT,RG, | |
22 | + DQ2,DP2,DANO,ALF,ALPHAS2 | |
23 | REAL X, Q, UPV, DNV, USEA, DSEA, SSEA, CSEA, BOT, TOP, GLUE | |
24 | REAL P2, ANO | |
25 | DIMENSION XX(6),Q2(6),ANO(3) | |
26 | CHARACTER*20 PARM(NCHDIM) | |
27 | #include "pdf/expdp.inc" | |
28 | + VAL(NCHDIM) | |
29 | C | |
30 | DATA NXX,XX/6, .0001,.001,.01,.1,.3,.700/ | |
31 | DATA NQ2,Q2/6, 2.25,10000.,39.9,2.25,39.9,39.9/ | |
32 | DATA P2/5./ | |
33 | DATA IP2/0/ | |
34 | DATA ANO/12.,117.,208./ | |
35 | C | |
36 | N6 = L6 | |
37 | IFLPRT = LFLPRT | |
38 | C | |
39 | CH1=4./9. | |
40 | CH2=1./9. | |
41 | C | |
42 | SCAL=SQRT(Q2(2)) | |
43 | C | |
44 | C First call to PDFSET to initialize COMMON/W505120/ | |
45 | PARM(1) = 'Init0' | |
46 | VAL(1) = 0.D0 | |
47 | CALL PDFSET(PARM,VAL) | |
48 | C. | |
49 | C Test new PDFLIB versions (4.0 or bigger) | |
50 | KPTY=1 | |
51 | DO KGRO = 1,NGRMAX | |
52 | IF(NPGSMX(KPTY,KGRO) .EQ. 0) GOTO 400 | |
53 | DO ISET=0,NPGSMX(KPTY,KGRO) | |
54 | IF((KPTY*KGRO).NE.1 .AND. ISET.LE.0) GOTO 300 | |
55 | WRITE (N6,*) ' ========================================== | |
56 | +============================================' | |
57 | FIRST=.TRUE. | |
58 | PARM(1) = 'Nptype' | |
59 | VAL(1) = KPTY | |
60 | PARM(2) = 'Ngroup' | |
61 | VAL(2) = KGRO | |
62 | PARM(3) = 'Nset' | |
63 | VAL(3) = ISET | |
64 | C now do the NPDF's | |
65 | PARM(4) = 'NAtype' | |
66 | VAL(4) = 4 | |
67 | PARM(5) = 'NAgroup' | |
68 | VAL(5) = 1 | |
69 | PARM(6) = 'NAset' | |
70 | VAL(6) = 1 | |
71 | C... test different NSET writings | |
72 | IF(ISET.GT. 5) PARM(3) = 'NSET' | |
73 | IF(ISET.GT.10) PARM(3) = 'nset' | |
74 | IF(ISET.GT.15) PARM(3) = 'NSet' | |
75 | IF(ISET.GT.20) PARM(3) = 'nsET' | |
76 | C | |
77 | CALL PDFSET(PARM,VAL) | |
78 | C | |
79 | DQ = SCAL | |
80 | ALF = ALPHAS2(DQ) | |
81 | WRITE(N6,4000) NPTYPE,NGROUP,NSET, | |
82 | + SFNAME(NPTYPE,NGROUP,NSET),MODECR(NPTYPE,NGROUP,NSET) | |
83 | WRITE(N6,4001) QCDL4,QCDL5,DQ,ALF | |
84 | 4000 FORMAT(/,' Nptype = ',I1,' Ngroup = ',I2, | |
85 | + ' Nset = ',I3,' Name = "',A8, | |
86 | + '" CrMode = ',I3) | |
87 | 4001 FORMAT(1H ,' QCDL4, QCDL5 = ',2F7.4, | |
88 | + ' Scale = ',F8.4,' alpha(s) = ',F7.5) | |
89 | WRITE(N6,4002) NATYPE,NAGROUP,NASET, | |
90 | + SFNAME(NATYPE,NAGROUP,NASET) | |
91 | 4002 FORMAT(/,' NAtype = ',I1,' NAgroup = ',I2, | |
92 | + ' NAset = ',I3,' Name = "',A8,'"') | |
93 | C | |
94 | 100 CONTINUE | |
95 | DO IA=1,3 | |
96 | IAFL=0 | |
97 | DO 200 IX=1,NXX | |
98 | X=XX(IX) | |
99 | Q=SQRT(Q2(IX)) | |
100 | DX=X | |
101 | DQ=Q | |
102 | DANO=ANO(IA) | |
103 | CALL STRUCTA(DX,DQ,DANO, | |
104 | + DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL) | |
105 | UPV=DUPV | |
106 | DNV=DDNV | |
107 | USEA=DUSEA | |
108 | DSEA=DDSEA | |
109 | SSEA=DSTR | |
110 | CSEA=DCHM | |
111 | GLUE=DGL | |
112 | C | |
113 | FFPU = UPV + USEA | |
114 | FFPUB = USEA | |
115 | FFPD = DNV + DSEA | |
116 | FFPDB = DSEA | |
117 | FFPS = SSEA | |
118 | FFPSB = SSEA | |
119 | FFPC = CSEA | |
120 | FFPCB = CSEA | |
121 | FFNU = FFPD | |
122 | FFNUB = FFPDB | |
123 | FFND = FFPU | |
124 | FFNDB = FFPUB | |
125 | FFNS = FFPS | |
126 | FFNSB = FFPSB | |
127 | FFNC = FFPC | |
128 | FFNCB = FFPCB | |
129 | FFPGL = GLUE | |
130 | SFFP=CH1*(FFPU+FFPUB+FFPC+FFPCB)+CH2*(FFPD+FFPDB+FFPS+FFPSB) | |
131 | SFFN=CH1*(FFNU+FFNUB+FFNC+FFNCB)+CH2*(FFND+FFNDB+FFNS+FFNSB) | |
132 | FACNP = 0. | |
133 | IF (SFFP .NE. 0.) FACNP=SFFN/SFFP | |
134 | UP = FFPU | |
135 | UV = FFPU - FFPUB | |
136 | US = FFPUB | |
137 | DP = FFPD | |
138 | DV = FFPD - FFPDB | |
139 | DS = FFPDB | |
140 | GL = FFPGL | |
141 | F2P = SFFP | |
142 | C | |
143 | IF (IAFL.EQ.0) WRITE(N6,2999) ANO(IA) | |
144 | 2999 FORMAT(1X,' A:',F5.0) | |
145 | IAFL=1 | |
146 | IF(Q2(IX).LT.100.) THEN | |
147 | WRITE(N6,3000) X,Q2(IX),FACNP,UP,UV,US,DP,DV,GL,F2P | |
148 | ELSE | |
149 | WRITE(N6,3001) X,Q2(IX),FACNP,UP,UV,US,DP,DV,GL,F2P | |
150 | ENDIF | |
151 | 3000 FORMAT(1X,' X:',F6.3,2X,'Q2:',F6.1, | |
152 | 1 3X,'N/P:',F7.4,2X,'UP:',F6.3, | |
153 | 1 1X,'UV:',F6.3,1X,'US:',F6.3,1X,'DP:',F6.3,1X,'DV:',F6.3, | |
154 | 2 1X,'GL:',F6.3,2X,'F2P:',F7.4) | |
155 | 3001 FORMAT(1X,' X:',F6.3,2X,'Q2:',F6.0, | |
156 | 1 3X,'N/P:',F7.4,2X,'UP:',F6.3, | |
157 | 1 1X,'UV:',F6.3,1X,'US:',F6.3,1X,'DP:',F6.3,1X,'DV:',F6.3, | |
158 | 2 1X,'GL:',F6.3,2X,'F2P:',F7.4) | |
159 | C | |
160 | 200 CONTINUE | |
161 | WRITE(N6,*) ' ' | |
162 | ENDDO | |
163 | C | |
164 | CALL PDFSTA | |
165 | C | |
166 | 300 CONTINUE | |
167 | ENDDO | |
168 | 400 CONTINUE | |
169 | ENDDO | |
170 | C ENDDO | |
171 | C | |
172 | STOP | |
173 | END |