PDF/tpdf/testgamo.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.2  1996/10/30 08:32:51  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 TESTGAMO
  14 C
  15 #include "pdf/impdp.inc"
  16       PARAMETER (LFLPRT = 2)
  17 C+SEQ, W5051P2.
  18       PARAMETER (MODEMX = 281)
  19 #include "pdf/w5051p7.inc"
  20 #include "pdf/w50510.inc"
  21 #include "pdf/w50511.inc"
  22 #include "pdf/w50512.inc"
  23 #include "pdf/w50513.inc"
  24 #include "pdf/w50516.inc"
  25 #include "pdf/w50517.inc"
  26 C
  27 #include "pdf/expdp.inc"
  28      +       DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DTB,DGL
  29       REAL       GLUCK
  30       COMMON/W5051I8/GLUCK
  31 #include "pdf/expdp.inc"
  32      +               XI(32),WI(32),XX(33)
  33       COMMON/W5051I9/XI,WI,XX,NTERMS
  34 #include "pdf/expdp.inc"
  35      +       XV(3),ZQ(5)
  36       DATA XV/0.05D0,0.50D0,0.95D0/
  37       DATA ZQ/10.D0,100.D0,1000.D0,1.D+4,5.5D+5/
  38       DATA FLAV/4.D0/
  39       DATA ALQED/137.D0/
  40       DATA ICALL/0/
  41       CHARACTER*20 PARM(20)
  42 #include "pdf/expdp.inc"
  43      +       VAL(20)
  44 C
  45       N6 = L6
  46       IFLPRT = LFLPRT
  47 C
  48       IF(GLUCK.LT.0.01) GLUCK = 0.
  49       PI=4.D0*ATAN(1.D0)
  50       PI2=PI**2
  51       CF=4.D0/3.D0
  52       NF=INT(FLAV+0.1D0)
  53       FLAVOR=FLAV
  54       B0=11.D0-2.D0*NF/3.D0
  55       B1=102.D0-38.D0*NF/3.D0
  56       CCOEG=2.D0/9.D0
  57       IF(NF.EQ.4) CCOEG=5.D0/18.D0
  58       COEG=2.D0*FLAVOR*CCOEG
  59       CALL WATE32
  60 C
  61       DO 300 ISTRUC=272,272
  62       ICALL = 0
  63       FIRST=.TRUE.
  64       VAL(1)  = ISTRUC
  65       PARM(1) = 'Mode'
  66       IF(ISTRUC.GE.201 .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          CALL PDFSET(PARM,VAL)
  75          Q02 = Q2MIN
  76          OWLAM = QCDL4
  77          OWLAM2=OWLAM**2
  78          IORD = LO - 1
  79 *
  80       DO 100 J=3,3
  81       Q2=Q02
  82       IF(J.GE.2) Q2=ZQ(J-1)
  83       ALQ2= LOG(Q2/OWLAM2)
  84       ALFPI= 2.D0 /(B0*ALQ2+B1* LOG(ALQ2)/B0)
  85       DO 70 K=1,3
  86       X=XV(K)
  87       DX = X
  88       Q = SQRT(Q2)
  89       DQ = Q
  90       CALL STRUCTM(DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DTB,DGL)
  91       UB =DUB
  92       DB =DDB
  93       SB =DSB
  94       CB =DCB
  95       GL =DGL
  96       F2=(4.D0*  DUB   +  DDB   +  DSB)/9.D0*2.D0
  97       F2C=F2+(4.D0/9.D0)*DCB*2.D0
  98       IF(ICALL.EQ.0) THEN
  99         ICALL = 1
 100         WRITE (N6,1000) Q02,IORD,OWLAM,GLUCK,FLAV
 101  1000   FORMAT(1X,' Q02= ',F7.3,'   IORD= ',I2,'   LAMBDA=',F7.4,
 102      +               '   GLUCK= ',F5.2,'   NF = ',F5.2)
 103         ENDIF
 104       IF (K.EQ.1) THEN
 105         WRITE (N6,1001) Q2
 106  1001   FORMAT(/,'  Q2=  ',E10.3,' GEV**2')
 107         WRITE(N6,1002)
 108  1002   FORMAT(/,4X,' X',6X,'X*UB  ',6X,'X*DB  ',6X,'X*SB  ',6X,'X*CB  '
 109      +                  ,6X,'X*GL',8X,' F2 ',6X,' F2/ALPHA'
 110      +                            ,5X,' F2C ',5X,' F2C/ALPHA')
 111       ENDIF
 112       IF (GLUCK.LT..5) GOTO 5
 113       F2=DUV
 114  5    IF(IORD)13,13,14
 115 14    CONTINUE
 116       AL=0.5D0*ALFPI
 117       AL1= LOG(1.D0-X)
 118       XQ=F2+F2*AL*CF*(-9.D0-2.D0*PI2/3.D0+AL1*(-3.D0+2.D0*AL1))
 119 C     WRITE(N6,1111) X,Q,F2,AL1,XQ
 120 C1111 FORMAT(1X,' X,F2 = ',5E13.5)
 121       DO 81 I=1,NTERMS
 122       Y=0.5D0*(1.D0-X)*XI(I)+0.5D0*(1.D0+X)
 123       XY=X/Y
 124       AL1= LOG(1.D0-Y)
 125       DX = XY
 126       CALL STRUCTM(DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DTB,DGL)
 127       XQQ=(4.D0*  DUB   +  DDB   +  DSB)/9.D0*2.D0
 128       IF (GLUCK.LT..5) GOTO 6
 129       XQQ=DUV
 130  6    C22=CF*(6.D0+4.D0*Y-2.D0*(1.D0+Y*Y)/
 131      +          (1.D0-Y)* LOG(Y)-2.D0*(1.D0+Y)*AL1)
 132       C23=CF*(-3.D0+4.D0*AL1)/(1.D0-Y)
 133       CG2=(-1.D0+8.D0*Y*(1.D0-Y)+(1.D0-2.D0*Y+2.D0*Y*Y)*
 134      +                   LOG(1.D0/Y-1.D0))*COEG
 135       XQ=XQ +0.5D0*(1.D0-X)*WI(I)*AL*(C22*XQQ+C23*(XQQ- F2))
 136       XQ=XQ +0.5D0*(1.D0-X)*WI(I)*AL*CG2*DGL
 137 C     WRITE(N6,1112) Y,AL1,C22,C23,CG2,XQQ,XQ
 138 C1112 FORMAT(1X,' Y,C22 = ',7E13.5)
 139  81   CONTINUE
 140       F2=3.D0*3.D0*(2.D0/27.D0)*(1.D0/(2.D0*PI*ALQED))*2.D0
 141       F2=XQ+F2*(8.D0*X*(1.D0-X)-1.D0+(1.D0-2.D0*X+2.D0*X*X)*
 142      +            LOG(1.D0/X-1.D0))*X
 143       F2C=F2+(4.D0/9.D0)*CB*2.D0
 144    13 F2N=F2*ALQED
 145       F2CN=F2C*ALQED
 146       WRITE(N6,1003)X,UB,DB,SB,CB,GL,F2,F2N,F2C,F2CN
 147  1003 FORMAT(1X,' ',F6.3,9E12.4)
 148    70 CONTINUE
 149   100 CONTINUE
 150       WRITE(N6,*) '  '
 151   300 CONTINUE
 152       STOP
 153       END