PDF/tpdf/testpdfo.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.2  1996/10/30 08:32:53  cernlib
   6 * Version 7.04
   7 *
   8 * Revision 1.1.1.1  1996/04/12 15:30:17  plothow
   9 * Version 7.01
  10 *
  11 *
  12 #include "pdf/pilot.h"
  13       PROGRAM TESTPDFO
  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,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL,
  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                WRITE (N6,*) ' ==========================================
  77      +================================================================='
  78       IF(ISTRUC.GE.0 .AND. ISTRUC.LE.MODEMX) THEN
  79          CALL PDFSET(PARM,VAL)
  80       ELSEIF (ISTRUC.EQ.-3) THEN
  81          GOTO 5
  82       ELSEIF (ISTRUC.EQ.-2) THEN
  83          PARM(1) = 'MODE'
  84          VAL(1)  = 45.D0
  85          CALL PDFSET(PARM,VAL)
  86       ELSEIF (ISTRUC.EQ.-1) THEN
  87          PARM(1) = 'mode'
  88          VAL(1)  = 45.D0
  89          CALL PDFSET(PARM,VAL)
  90       ENDIF
  91 C
  92       DQ = SCAL
  93       ALF = ALPHAS2(DQ)
  94       WRITE (N6,2000) MODE,QCDL4,QCDL5,DQ,ALF
  95  2000 FORMAT(/,'  Mode = ',I3,'  QCDL4, QCDL5 = ',2F7.4,
  96      +         '  SCALE = ',F7.4,'  alpha(s) = ',F8.5)
  97 C
  98    5  CONTINUE
  99       DO 200 IX=1,NXX
 100       X=XX(IX)
 101       Q=SQRT(Q2(IX))
 102         DX=X
 103         DQ=Q
 104        CALL STRUCTM(DX,DQ,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       WRITE(N6,1000)X,Q2(IX),FACNP,UP,UV,US,DP,DV,GL,F2P
 144  1000 FORMAT(1X,' X:',F6.3,2X,'Q2:',F5.1,
 145      1    3X,'N/P:',F7.4,2X,'UP:',F6.3,
 146      1    1X,'UV:',F6.3,1X,'US:',F6.3,1X,'DP:',F6.3,1X,'DV:',F6.3,
 147      2    1X,'GL:',F6.3,2X,'F2P:',F7.4)
 148 C
 149  200  CONTINUE
 150 C
 151       IF(ISTRUC.EQ.-3) THEN
 152         DQ = SCAL
 153         ALF = ALPHAS2(DQ)
 154         WRITE (N6,2000) MODE,QCDL4,QCDL5,DQ,ALF
 155       ENDIF
 156       CALL PDFSTA
 157 C
 158  300  CONTINUE
 159 C
 160       STOP
 161       END