[u/mrichter/AliRoot.git] / LHAPDF / lhapdf5.2.2 / wrapdog.f

      subroutine DOGevolvep0(xin,qin,p2in,ip2in,pdf)
      include 'parmsetup.inc'
      real*8 xin,qin,q2in,p2in,pdf(-6:6),xval(45),qcdl4,qcdl5
      real*8 upv,dnv,usea,dsea,str,chm,bot,top,glu
      character*16 name(nmxset)
      integer nmem(nmxset),ndef(nmxset),mmem
      common/NAME/name,nmem,ndef,mmem
      integer ns
      
      save 

      call DOPHO1(xin,qin,upv,dnv,usea,dsea,str,chm,bot,glu)

      pdf(-6)= 0.0d0
      pdf(6)= 0.0d0
      pdf(-5)= bot
      pdf(5 )= bot
      pdf(-4)= chm
      pdf(4 )= chm
      pdf(-3)= str
      pdf(3 )= str
      pdf(-2)= usea
      pdf(2 )= upv
      pdf(-1)= dsea
      pdf(1 )= dnv
      pdf(0 )= glu
      
      return
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
      entry DOGevolvep1(xin,qin,p2in,ip2in,pdf)
      
      call DOPHO2(xin,qin,upv,dnv,usea,dsea,str,chm,bot,glu)

      pdf(-6)= 0.0d0
      pdf(6)= 0.0d0
      pdf(-5)= bot
      pdf(5 )= bot
      pdf(-4)= chm
      pdf(4 )= chm
      pdf(-3)= str
      pdf(3 )= str
      pdf(-2)= usea
      pdf(2 )= upv
      pdf(-1)= dsea
      pdf(1 )= dnv
      pdf(0 )= glu
      
      return
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
      entry DOGread(nset)
      read(1,*)nmem(nset),ndef(nset)
      return
c
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
      entry DOGalfa(alfas,qalfa)
        call getnset(iset)
	call GetOrderAsM(iset,iord)
        call Getlam4M(iset,imem,qcdl4)
        call Getlam5M(iset,imem,qcdl5)
        call aspdflib(alfas,Qalfa,iord,qcdl5)
      return
c
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
      entry DOGinit(Eorder,Q2fit)
      return
c
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
      entry DOGpdf(mem)
      imem = mem
      return
c
 1000 format(5e13.5)
      end
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
       SUBROUTINE DOPHO1(DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL)
C********************************************************************
C*                                                                  *
C*    Parametrization of parton distribution functions              *
C*    in the photon (LO analysis) - asymptotic solution of AP eq.!  *
C*                                                                  *
C* authors:  D.Duke and H.Owens (DO)                                *
C*           /Phys.Rev. D26 (1982) 1600/                            *
C*                                                                  *
C* Prepared by:                                                     *
C*             Krzysztof Charchula, DESY                            *
C*             bitnet: F1PCHA@DHHDESY3                              *
C*             decnet: 13313::CHARCHULA                             *
C*                                                                  *
C* Modified by:                                                     *
C*             H. Plothow-Besch/CERN-PPE                            *
C*                                                                  *
C********************************************************************
C
      implicit real*8 (a-h,o-z)
       double precision
     +        CQ(5),
     +        DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL
      PARAMETER (ALPEM=7.29927D-3, PI=3.141592D0)
      PARAMETER (ALAM=0.2D0)
      DATA CQ/0.33333D0,0.66666D0,0.33333D0,0.66666D0,0.33333D0/
C
       Q2 = DQ*DQ
       ALAM2=ALAM**2
       FQ=ALPEM/(2.*PI)*LOG(Q2/ALAM2)
C
C...gluons
       POMG=0.194*(1.-DX)**1.03/(DX**0.97)
       DGL=POMG*FQ
C
C...quarks
        POM1=(1.81-1.67*DX+2.16*DX**2)
        POM2=DX**0.7/(1.-0.4*LOG(1.-DX))
        POM3=38.D-4*(1.-DX)**1.82/(DX**1.18)
          DDB=(CQ(1)**2*POM1*POM2+POM3)*FQ
          DDV=DDB
          DUB=(CQ(2)**2*POM1*POM2+POM3)*FQ
          DUV=DUB
          DSB=(CQ(3)**2*POM1*POM2+POM3)*FQ
          DCB=(CQ(4)**2*POM1*POM2+POM3)*FQ
          DBB=(CQ(5)**2*POM1*POM2+POM3)*FQ
       RETURN
       END
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
       SUBROUTINE DOPHO2(DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL)
C********************************************************************
C*                                                                  *
C*    Parametrization of parton distribution functions              *
C*    in the photon (LO analysis) - asymptotic solution of AP eq.!  *
C*                                                                  *
C* authors:  D.Duke and H.Owens (DO)                                *
C*           /Phys.Rev. D26 (1982) 1600/                            *
C*                                                                  *
C* Prepared by:                                                     *
C*             Krzysztof Charchula, DESY                            *
C*             bitnet: F1PCHA@DHHDESY3                              *
C*             decnet: 13313::CHARCHULA                             *
C*                                                                  *
C* Modified by:                                                     *
C*             H. Plothow-Besch/CERN-PPE                            *
C*                                                                  *
C********************************************************************
C
      implicit real*8 (a-h,o-z)
      double precision
     +        CQ(5),
     +        DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL
      PARAMETER (ALPEM=7.29927D-3,PI=3.141592D0)
      PARAMETER (ALAM=0.4D0)
      DATA CQ/0.33333D0,0.66666D0,0.33333D0,0.66666D0,0.33333D0/
C
       Q2 = DQ*DQ
       ALAM2=ALAM**2
       FQ=ALPEM/(2.*PI)*LOG(Q2/ALAM2)
C
C...gluons
       POMG=0.194*(1.-DX)**1.03/(DX**0.97)
       DGL=POMG*FQ
C
C...quarks
        POM1=(1.81-1.67*DX+2.16*DX**2)
        POM2=DX**0.7/(1.-0.4*LOG(1.-DX))
        POM3=38.D-4*(1.-DX)**1.82/(DX**1.18)
          DDB=(CQ(1)**2*POM1*POM2+POM3)*FQ
          DDV=DDB
          DUB=(CQ(2)**2*POM1*POM2+POM3)*FQ
          DUV=DUB
          DSB=(CQ(3)**2*POM1*POM2+POM3)*FQ
          DCB=(CQ(4)**2*POM1*POM2+POM3)*FQ
          DBB=(CQ(5)**2*POM1*POM2+POM3)*FQ
       RETURN
       END
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
Commit	Line	Data
3c5d1739	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.67DX+2.16DX**2)
112	POM2=DX*0.7/(1.-0.4LOG(1.-DX))
113	POM3=38.D-4(1.-DX)1.82/(DX*1.18)
114	DDB=(CQ(1)*2POM1POM2+POM3)FQ
115	DDV=DDB
116	DUB=(CQ(2)*2POM1POM2+POM3)FQ
117	DUV=DUB
118	DSB=(CQ(3)*2POM1POM2+POM3)FQ
119	DCB=(CQ(4)*2POM1POM2+POM3)FQ
120	DBB=(CQ(5)*2POM1POM2+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.67DX+2.16DX**2)
161	POM2=DX*0.7/(1.-0.4LOG(1.-DX))
162	POM3=38.D-4(1.-DX)1.82/(DX*1.18)
163	DDB=(CQ(1)*2POM1POM2+POM3)FQ
164	DDV=DDB
165	DUB=(CQ(2)*2POM1POM2+POM3)FQ
166	DUV=DUB
167	DSB=(CQ(3)*2POM1POM2+POM3)FQ
168	DCB=(CQ(4)*2POM1POM2+POM3)FQ
169	DBB=(CQ(5)*2POM1POM2+POM3)FQ
170	RETURN
171	END
172	ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc