[u/mrichter/AliRoot.git] / PDF / tpdf / testpdf.F

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