PDF/spdf/dgpho4.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.4  2000/09/18 10:02:35  hristov
   6 * Makefile added to PDF8
   7 *
   8 * Revision 1.1.1.2  1996/10/30 08:28:05  cernlib
   9 * Version 7.04
  10 *
  11 * Revision 1.1.1.1  1996/04/12 15:29:17  plothow
  12 * Version 7.01
  13 *
  14 *
  15 #include "pdf/pilot.h"
  16        SUBROUTINE DGPHO4(DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL)
  17 C********************************************************************
  18 C*                                                                  *
  19 C*    Parametrization of parton distribution functions              *
  20 C*    in the photon (LO analysis) - full  solution of AP eq.!       *
  21 C*                                                                  *
  22 C* authors:  M.Drees and K.Grassie (DG)                             *
  23 C*          /Z. Phys. C28 (1985) 451/                               *
  24 C*                                                                  *
  25 C* Prepared by:                                                     *
  26 C*             Krzysztof Charchula, DESY                            *
  27 C*             bitnet: F1PCHA@DHHDESY3                              *
  28 C*             decnet: 13313::CHARCHULA                             *
  29 C*                                                                  *
  30 C* Modified by:                                                     *
  31 C*             H. Plothow-Besch/CERN-PPE                            *
  32 C*                                                                  *
  33 C********************************************************************
  34 C
  35 #include "pdf/impdp.inc"
  36 #if defined(CERNLIB_DOUBLE)
  37       DOUBLE PRECISION
  38 #endif
  39 #if defined(CERNLIB_SINGLE)
  40       REAL
  41 #endif
  42      +        A(3,4,3),AT(3),
  43      +        B(5,4,2,3),BT(5,2),XQPOM(2),E(2),
  44      +        DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL
  45       PARAMETER (ALPEM=7.29927D-3, PI=3.141592D0)
  46       PARAMETER (ALAM=0.4D0)
  47 C...comments
  48 C...--------------------------------------------------
  49 C...        with nf=5 (valid for 200< Q2 <1D4 GeV2)
  50 C...--------------------------------------------------
  51 C
  52 C...initialization of gluon parameters array for DG
  53         DATA (((A(I,J,K),I=1,3),J=1,4),K=1,3)/
  54      >    -0.20700, -0.19870,  5.1190,
  55      >     0.61580,  0.62570, -0.2752,
  56      >     1.07400,  8.35200, -6.9930,
  57      >     0.00000,  5.02400,  2.2980,
  58      >     0.8926D-2,0.0509,  -0.2313,
  59      >     0.65940,  0.27740,  0.1382,
  60      >     0.47660, -0.39060,  6.5420,
  61      >     0.01975, -0.32120,  0.5162,
  62      >     0.03197, -0.618D-2,-0.1216,
  63      >     1.01800,  0.94760,  0.9047,
  64      >     0.24610, -0.60940,  2.6530,
  65      >     0.02707, -0.01067,  0.2003D-2/
  66 C
  67 C...initialization of quark parameters array for DG
  68         DATA (((B(I,J,K,1),I=1,5),J=1,4),K=1,2)/
  69      >     2.2850,   6.0730,  -0.4202,   -0.0808,  0.0553,
  70      >    -0.0153,  -0.8132,   0.0178,    0.6346,  1.1360,
  71      >     1.33D3, -41.310,    0.9216,    1.2080,  0.9512,
  72      >     4.2190,   3.1650,   0.1800,    0.2030,  0.0116,
  73      >    16.690,    0.1760,  -0.0208,   -0.0168, -0.1986,
  74      >    -0.7916,   0.0479,   0.3386D-2, 1.3530,  1.1000,
  75      >     1.0990D3, 1.0470,   4.8530,    1.4260,  1.1360,
  76      >     4.4280,   0.0250,   0.8404,    1.2390, -0.2779/
  77         DATA (((B(I,J,K,2),I=1,5),J=1,4),K=1,2)/
  78      >    -0.3711,  -0.1717,   0.08766,  -0.8915, -0.1816,
  79      >     1.0610,   0.7815,   0.02197,   0.2857,  0.5866,
  80      >     4.7580,   1.5350,   0.10960,   2.9730,  2.4210,
  81      >    -0.0150,   0.7067D-2,0.20400,   0.1185,  0.4059,
  82      >    -0.1207,  25.000,   -0.01230,  -0.0919,  0.02015,
  83      >     1.0710,  -1.6480,   1.16200,   0.7912,  0.9869,
  84      >     1.9770,  -0.01563,  0.48240,   0.6397, -0.07036,
  85      >    -0.8625D-2,6.4380,  -0.01100,   2.3270,  0.01694/
  86         DATA (((B(I,J,K,3),I=1,5),J=1,4),K=1,2)/
  87      >    15.8,      2.742,    0.02917,  -0.0342, -0.02302,
  88      >    -0.9464,  -0.7332,   0.04657,   0.7196,  0.9229,
  89      >    -0.5,      0.7148,   0.1785,    0.7338,  0.5873,
  90      >    -0.2118,   3.287,    0.04811,   0.08139,-0.79D-4,
  91      >     6.734,   59.88,    -0.3226D-2,-0.03321, 0.1059,
  92      >    -1.008,   -2.983,    0.8432,    0.9475,  0.6954,
  93      >    -0.08594,  4.48,     0.3616,   -0.3198, -0.6663,
  94      >     0.07625,  0.9686,   0.1383D-2, 0.02132, 0.3683/
  95 C
  96 C...specification of sets
  97          NFL=5
  98 C
  99 C...calculations
 100        Q2 = DQ*DQ
 101        ALAM2=ALAM**2
 102        T=LOG(Q2/ALAM2)
 103        LF=NFL-2
 104 C
 105 C...gluons
 106         DO 11 I=1,3
 107           AT(I)=A(I,1,LF)*T**A(I,2,LF)+A(I,3,LF)*T**(-A(I,4,LF))
 108  11     CONTINUE
 109         POMG=AT(1)*DX**AT(2)*(1.D0-DX)**AT(3)
 110         DGL=POMG*ALPEM
 111 C
 112 C...quarks
 113         E(1)=1.D0
 114         E(2)=55.D0/6.D0
 115         DO 13 J=1,2
 116           DO 15 I=1,5
 117             BTP=B(I,1,J,LF)*T**B(I,2,J,LF)
 118             BT(I,J)=BTP+B(I,3,J,LF)*T**(-B(I,4,J,LF))
 119  15       CONTINUE
 120  13     CONTINUE
 121 C
 122 C...singlet & non-singlet combinations
 123         DO 17 J=1,2
 124           POM1=DX*(DX*DX+(1.D0-DX)**2)/(BT(1,J)-BT(2,J)*LOG(1.D0-DX))
 125           POM2=BT(3,J)*DX**BT(4,J)*(1.D0-DX)**BT(5,J)
 126           XQPOM(J)=E(J)*POM1+POM2
 127  17     CONTINUE
 128 C
 129 C...quarks flavours
 130       DUB=ALPEM*1.D0/10.D0*(XQPOM(2)+15.D0/2.D0*XQPOM(1))
 131       DUV=DUB
 132       DCB=DUB
 133       DDB=ALPEM*1.D0/10.D0*(XQPOM(2)-5.D0*XQPOM(1))
 134       DDV=DDB
 135       DSB=DDB
 136       DCB=DUB
 137       DBB=DDB
 138 C
 139       RETURN
 140       END