LHAPDF/lhapdf5.3.1/wrapdog.f

   1       subroutine DOGevolvep0(xin,qin,p2in,ip2in,pdf)
   2       include 'parmsetup.inc'
   3       real*8 xin,qin,q2in,p2in,pdf(-6:6),xval(45),qcdl4,qcdl5
   4       real*8 upv,dnv,usea,dsea,str,chm,bot,top,glu
   5       character*16 name(nmxset)
   6       integer nmem(nmxset),ndef(nmxset),mmem
   7       common/NAME/name,nmem,ndef,mmem
   8       integer ns
   9
  10       save
  11
  12       call DOPHO1(xin,qin,upv,dnv,usea,dsea,str,chm,bot,glu)
  13
  14       pdf(-6)= 0.0d0
  15       pdf(6)= 0.0d0
  16       pdf(-5)= bot
  17       pdf(5 )= bot
  18       pdf(-4)= chm
  19       pdf(4 )= chm
  20       pdf(-3)= str
  21       pdf(3 )= str
  22       pdf(-2)= usea
  23       pdf(2 )= upv
  24       pdf(-1)= dsea
  25       pdf(1 )= dnv
  26       pdf(0 )= glu
  27
  28       return
  29 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
  30       entry DOGevolvep1(xin,qin,p2in,ip2in,pdf)
  31
  32       call DOPHO2(xin,qin,upv,dnv,usea,dsea,str,chm,bot,glu)
  33
  34       pdf(-6)= 0.0d0
  35       pdf(6)= 0.0d0
  36       pdf(-5)= bot
  37       pdf(5 )= bot
  38       pdf(-4)= chm
  39       pdf(4 )= chm
  40       pdf(-3)= str
  41       pdf(3 )= str
  42       pdf(-2)= usea
  43       pdf(2 )= upv
  44       pdf(-1)= dsea
  45       pdf(1 )= dnv
  46       pdf(0 )= glu
  47
  48       return
  49 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
  50       entry DOGread(nset)
  51       read(1,*)nmem(nset),ndef(nset)
  52       return
  53 c
  54 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
  55       entry DOGalfa(alfas,qalfa)
  56         call getnset(iset)
  57         call GetOrderAsM(iset,iord)
  58         call Getlam4M(iset,imem,qcdl4)
  59         call Getlam5M(iset,imem,qcdl5)
  60         call aspdflib(alfas,Qalfa,iord,qcdl5)
  61       return
  62 c
  63 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
  64       entry DOGinit(Eorder,Q2fit)
  65       return
  66 c
  67 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
  68       entry DOGpdf(mem)
  69       imem = mem
  70       return
  71 c
  72  1000 format(5e13.5)
  73       end
  74 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
  75        SUBROUTINE DOPHO1(DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL)
  76 C********************************************************************
  77 C*                                                                  *
  78 C*    Parametrization of parton distribution functions              *
  79 C*    in the photon (LO analysis) - asymptotic solution of AP eq.!  *
  80 C*                                                                  *
  81 C* authors:  D.Duke and H.Owens (DO)                                *
  82 C*           /Phys.Rev. D26 (1982) 1600/                            *
  83 C*                                                                  *
  84 C* Prepared by:                                                     *
  85 C*             Krzysztof Charchula, DESY                            *
  86 C*             bitnet: F1PCHA@DHHDESY3                              *
  87 C*             decnet: 13313::CHARCHULA                             *
  88 C*                                                                  *
  89 C* Modified by:                                                     *
  90 C*             H. Plothow-Besch/CERN-PPE                            *
  91 C*                                                                  *
  92 C********************************************************************
  93 C
  94       implicit real*8 (a-h,o-z)
  95        double precision
  96      +        CQ(5),
  97      +        DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL
  98       PARAMETER (ALPEM=7.29927D-3, PI=3.141592D0)
  99       PARAMETER (ALAM=0.2D0)
 100       DATA CQ/0.33333D0,0.66666D0,0.33333D0,0.66666D0,0.33333D0/
 101 C
 102        Q2 = DQ*DQ
 103        ALAM2=ALAM**2
 104        FQ=ALPEM/(2.*PI)*LOG(Q2/ALAM2)
 105 C
 106 C...gluons
 107        POMG=0.194*(1.-DX)**1.03/(DX**0.97)
 108        DGL=POMG*FQ
 109 C
 110 C...quarks
 111         POM1=(1.81-1.67*DX+2.16*DX**2)
 112         POM2=DX**0.7/(1.-0.4*LOG(1.-DX))
 113         POM3=38.D-4*(1.-DX)**1.82/(DX**1.18)
 114           DDB=(CQ(1)**2*POM1*POM2+POM3)*FQ
 115           DDV=DDB
 116           DUB=(CQ(2)**2*POM1*POM2+POM3)*FQ
 117           DUV=DUB
 118           DSB=(CQ(3)**2*POM1*POM2+POM3)*FQ
 119           DCB=(CQ(4)**2*POM1*POM2+POM3)*FQ
 120           DBB=(CQ(5)**2*POM1*POM2+POM3)*FQ
 121        RETURN
 122        END
 123 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
 124        SUBROUTINE DOPHO2(DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL)
 125 C********************************************************************
 126 C*                                                                  *
 127 C*    Parametrization of parton distribution functions              *
 128 C*    in the photon (LO analysis) - asymptotic solution of AP eq.!  *
 129 C*                                                                  *
 130 C* authors:  D.Duke and H.Owens (DO)                                *
 131 C*           /Phys.Rev. D26 (1982) 1600/                            *
 132 C*                                                                  *
 133 C* Prepared by:                                                     *
 134 C*             Krzysztof Charchula, DESY                            *
 135 C*             bitnet: F1PCHA@DHHDESY3                              *
 136 C*             decnet: 13313::CHARCHULA                             *
 137 C*                                                                  *
 138 C* Modified by:                                                     *
 139 C*             H. Plothow-Besch/CERN-PPE                            *
 140 C*                                                                  *
 141 C********************************************************************
 142 C
 143       implicit real*8 (a-h,o-z)
 144       double precision
 145      +        CQ(5),
 146      +        DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL
 147       PARAMETER (ALPEM=7.29927D-3,PI=3.141592D0)
 148       PARAMETER (ALAM=0.4D0)
 149       DATA CQ/0.33333D0,0.66666D0,0.33333D0,0.66666D0,0.33333D0/
 150 C
 151        Q2 = DQ*DQ
 152        ALAM2=ALAM**2
 153        FQ=ALPEM/(2.*PI)*LOG(Q2/ALAM2)
 154 C
 155 C...gluons
 156        POMG=0.194*(1.-DX)**1.03/(DX**0.97)
 157        DGL=POMG*FQ
 158 C
 159 C...quarks
 160         POM1=(1.81-1.67*DX+2.16*DX**2)
 161         POM2=DX**0.7/(1.-0.4*LOG(1.-DX))
 162         POM3=38.D-4*(1.-DX)**1.82/(DX**1.18)
 163           DDB=(CQ(1)**2*POM1*POM2+POM3)*FQ
 164           DDV=DDB
 165           DUB=(CQ(2)**2*POM1*POM2+POM3)*FQ
 166           DUV=DUB
 167           DSB=(CQ(3)**2*POM1*POM2+POM3)*FQ
 168           DCB=(CQ(4)**2*POM1*POM2+POM3)*FQ
 169           DBB=(CQ(5)**2*POM1*POM2+POM3)*FQ
 170        RETURN
 171        END
 172 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc