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