[u/mrichter/AliRoot.git] / LHAPDF / lhapdf5.3.1 / wrapgsg96.f

      subroutine GSG96evolvep0(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
      real*8 SIG,QNS,GL
      common/gsgdat/SIG(78,11,3),QNS(78,11,3),GL(78,11,3)
      character*16 name(nmxset)
      integer nmem(nmxset),ndef(nmxset),mmem
      common/NAME/name,nmem,ndef,mmem
      integer nset
     
      save 
      
      iimem = 2
      call GS96HL(iimem,xin,qin,upv,dnv,str,chm,bot,glu)
      usea = upv
      dsea = dnv

      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 GSG96evolvep1(xin,qin,p2in,ip2in,pdf)

      iimem = 1
      call GS96HL(iimem,xin,qin,upv,dnv,str,chm,bot,glu)
      usea = upv
      dsea = dnv

      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 GSG96read(nset)
      read(1,*)nmem(nset),ndef(nset)
      
      do m = 1,2
      do k = 1,78
        read(1,*)(sig(k,j,m),j=1,4)
        read(1,*)(sig(k,j,m),j=5,8)
        read(1,*)(sig(k,j,m),j=9,11)
      enddo
      do k = 1,78
        read(1,*)(qns(k,j,m),j=1,4)
        read(1,*)(qns(k,j,m),j=5,8)
        read(1,*)(qns(k,j,m),j=9,11)
      enddo
      do k = 1,78
        read(1,*)(gl(k,j,m),j=1,4)
        read(1,*)(gl(k,j,m),j=5,8)
        read(1,*)(gl(k,j,m),j=9,11)
      enddo
      enddo
      return
c
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
      entry GSG96alfa(alfas,qalfa)
       call getnset(iset)
       call getnmem(iset,imem)
       call GetOrderAsM(iset,iord)
       call Getlam4M(iset,imem,qcdl4)
       call Getlam5M(iset,imem,qcdl5)
       call aspdflib(alfas,Qalfa,iord,qcdl5)
      return
c
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
      entry GSG96init(Eorder,Q2fit)
      return
c
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
      entry GSG96pdf(mem)
      call getnset(iset)
      call setnmem(iset,mem)
c      imem = mem
      return
c
 1000 format(5e13.5)
      end
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c      SUBROUTINE GS96HO(X,Q,U,D,S,C,B,G)
      SUBROUTINE GS96HL(iset,X,Q,U,D,S,C,B,G)
      implicit real*8 (a-h,o-z)
      PARAMETER(NP=78,NQ=11,NARG=2)
*****************************************************************
* Subroutine returns the parton distributions in the photon in  *
* next-to-leading order. u,d etc. gives the actual distributions*
* not x times the distributions; Q2 means Q^2. The distributions*
* are valid for 5.0e^-4< x < 1.0 and 5.3 GeV^2 < Q^2 < 1.0e^8.  *
* if higher Q^2 or lower x is required, these may be obtained   *
* from the authors on request.                                  * 
* Lionel Gordon April 1996 : Gordon@hep.anl.gov                 *
* John Storrow             :johns@a3.ph.man.ac.uk               *
*****************************************************************
c      DIMENSION SIG(NP,NQ),QNS(NP,NQ),GL(NP,NQ),Y(NP)
      DIMENSION Y(NP)
      DIMENSION XT(NARG),NA(NARG),A(NP+NQ),QT(NQ)
      common/gsgdat/SIG(78,11,3),QNS(78,11,3),GL(78,11,3)
      EXTERNAL GS2XCOR
c      SAVE SIG,QNS,GL,Y,ICALL
      SAVE Y,ICALL
      DATA QT /3.0D0,20.0D0,50.0D0,100.0D0,500.0D0,1.0D3,1.0D4,1.0D5,
     + 1.0D6,1.0D7,1.0D8/
******************************************************************
        q2=q*q
*if x is out of range
       IF((X.LT.5.0D-4).OR.(X.GT.1.0D0)) GOTO 90
*******************************************************************
       IF (ICALL.NE.1) THEN
* get the x coordinates
       CALL GS2XCOR(Y,NP)
        ICALL=1
        END IF
*
      DO 30 IX=1,NP
        A(IX)=Y(IX)
   30 CONTINUE
      DO 40 IQ=1,NQ
        A(NP+IQ)=QT(IQ)
   40 CONTINUE
*
      NA(1)=NP
      NA(2)=NQ  
       XT(1)=X
       XT(2)=Q2
      XSIG=DBFINT(2,XT,NA,A,SIG(1,1,iset))
      XQNS=DBFINT(2,XT,NA,A,QNS(1,1,iset))
        G =DBFINT(2,XT,NA,A,GL(1,1,iset))
*
      IF (Q2.LT.50.0D0) THEN 
C Use three flavour evolution. 
       U=(XSIG+9.0D0*XQNS)/6.0D0
       D=(XSIG-4.5D0*XQNS)/6.0D0
       S=D
       C=0.0D0
       B=0.0D0
*      
      ELSE IF((Q2.GT.50.0).AND.(Q2.LT.250.0)) THEN
C Use four flavour evolution 
      U=(XSIG+6.0D0*XQNS)/8.0D0
      D=(XSIG-6.0D0*XQNS)/8.0D0
      S=D
      C=U
      B=0.0D0
      ELSE
C Use five flavour evolution
      U=(XSIG+7.5D0*XQNS)/10.0D0
      D=(XSIG-5.0D0*XQNS)/10.0D0
      S=D
      C=U
      B=D
      ENDIF
      U=X*U
      US=U
      D=X*D
      DS=D
      S=X*S
      C=X*C
      B=X*B
      G=X*G
 90   RETURN
      END
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
      SUBROUTINE GS2XCOR(Y,NP)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION Y(NP)
      N=1
      DO 10 IX=1,20,2
         Y(N)=DBLE(IX)/2000.0D0
      N=N+1
   10 CONTINUE
      DO 20 IX=30,200,10
         Y(N)=DBLE(IX)/2000.0D0
      N=N+1
   20 CONTINUE
      DO 30 IX=240,1600,40
         Y(N)=REAL(IX)/2000.0D0
      N=N+1
   30 CONTINUE
      DO 40 IX=1625,1980,25
         Y(N)=REAL(IX)/2000.0D0
      N=N+1
   40 CONTINUE
      RETURN
      END
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
Commit	Line	Data
4e9e3152	1	subroutine GSG96evolvep0(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	real*8 SIG,QNS,GL
	6	common/gsgdat/SIG(78,11,3),QNS(78,11,3),GL(78,11,3)
	7	character*16 name(nmxset)
	8	integer nmem(nmxset),ndef(nmxset),mmem
	9	common/NAME/name,nmem,ndef,mmem
	10	integer nset
	11
	12	save
	13
	14	iimem = 2
	15	call GS96HL(iimem,xin,qin,upv,dnv,str,chm,bot,glu)
	16	usea = upv
	17	dsea = dnv
	18
	19	pdf(-6)= 0.0d0
	20	pdf(6)= 0.0d0
	21	pdf(-5)= bot
	22	pdf(5 )= bot
	23	pdf(-4)= chm
	24	pdf(4 )= chm
	25	pdf(-3)= str
	26	pdf(3 )= str
	27	pdf(-2)= usea
	28	pdf(2 )= upv
	29	pdf(-1)= dsea
	30	pdf(1 )= dnv
	31	pdf(0 )= glu
	32
	33	return
	34	ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
	35	entry GSG96evolvep1(xin,qin,p2in,ip2in,pdf)
	36
	37	iimem = 1
	38	call GS96HL(iimem,xin,qin,upv,dnv,str,chm,bot,glu)
	39	usea = upv
	40	dsea = dnv
	41
	42	pdf(-6)= 0.0d0
	43	pdf(6)= 0.0d0
	44	pdf(-5)= bot
	45	pdf(5 )= bot
	46	pdf(-4)= chm
	47	pdf(4 )= chm
	48	pdf(-3)= str
	49	pdf(3 )= str
	50	pdf(-2)= usea
	51	pdf(2 )= upv
	52	pdf(-1)= dsea
	53	pdf(1 )= dnv
	54	pdf(0 )= glu
	55
	56	return
	57	ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
	58	entry GSG96read(nset)
	59	read(1,*)nmem(nset),ndef(nset)
	60
	61	do m = 1,2
	62	do k = 1,78
	63	read(1,*)(sig(k,j,m),j=1,4)
	64	read(1,*)(sig(k,j,m),j=5,8)
65	read(1,*)(sig(k,j,m),j=9,11)
66	enddo
67	do k = 1,78
68	read(1,*)(qns(k,j,m),j=1,4)
69	read(1,*)(qns(k,j,m),j=5,8)
70	read(1,*)(qns(k,j,m),j=9,11)
71	enddo
72	do k = 1,78
73	read(1,*)(gl(k,j,m),j=1,4)
74	read(1,*)(gl(k,j,m),j=5,8)
75	read(1,*)(gl(k,j,m),j=9,11)
76	enddo
77	enddo
78	return
79	c
80	ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
81	entry GSG96alfa(alfas,qalfa)
82	call getnset(iset)
83	call getnmem(iset,imem)
84	call GetOrderAsM(iset,iord)
85	call Getlam4M(iset,imem,qcdl4)
86	call Getlam5M(iset,imem,qcdl5)
87	call aspdflib(alfas,Qalfa,iord,qcdl5)
88	return
89	c
90	ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
91	entry GSG96init(Eorder,Q2fit)
92	return
93	c
94	ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
95	entry GSG96pdf(mem)
96	call getnset(iset)
97	call setnmem(iset,mem)
98	c imem = mem
99	return
100	c
101	1000 format(5e13.5)
102	end
103	ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
104	c SUBROUTINE GS96HO(X,Q,U,D,S,C,B,G)
105	SUBROUTINE GS96HL(iset,X,Q,U,D,S,C,B,G)
106	implicit real*8 (a-h,o-z)
107	PARAMETER(NP=78,NQ=11,NARG=2)
108	*****************************************************************
109	* Subroutine returns the parton distributions in the photon in *
110	* next-to-leading order. u,d etc. gives the actual distributions*
111	* not x times the distributions; Q2 means Q^2. The distributions*
112	* are valid for 5.0e^-4< x < 1.0 and 5.3 GeV^2 < Q^2 < 1.0e^8. *
113	* if higher Q^2 or lower x is required, these may be obtained *
114	* from the authors on request. *
115	* Lionel Gordon April 1996 : Gordon@hep.anl.gov *
116	* John Storrow :johns@a3.ph.man.ac.uk *
117	*****************************************************************
118	c DIMENSION SIG(NP,NQ),QNS(NP,NQ),GL(NP,NQ),Y(NP)
119	DIMENSION Y(NP)
120	DIMENSION XT(NARG),NA(NARG),A(NP+NQ),QT(NQ)
121	common/gsgdat/SIG(78,11,3),QNS(78,11,3),GL(78,11,3)
122	EXTERNAL GS2XCOR
123	c SAVE SIG,QNS,GL,Y,ICALL
124	SAVE Y,ICALL
125	DATA QT /3.0D0,20.0D0,50.0D0,100.0D0,500.0D0,1.0D3,1.0D4,1.0D5,
126	+ 1.0D6,1.0D7,1.0D8/
127	******************************************************************
128	q2=q*q
129	*if x is out of range
130	IF((X.LT.5.0D-4).OR.(X.GT.1.0D0)) GOTO 90
131	*******************************************************************
132	IF (ICALL.NE.1) THEN
133	* get the x coordinates
134	CALL GS2XCOR(Y,NP)
135	ICALL=1
136	END IF
137	*
138	DO 30 IX=1,NP
139	A(IX)=Y(IX)
140	30 CONTINUE
141	DO 40 IQ=1,NQ
142	A(NP+IQ)=QT(IQ)
143	40 CONTINUE
144	*
145	NA(1)=NP
146	NA(2)=NQ
147	XT(1)=X
148	XT(2)=Q2
149	XSIG=DBFINT(2,XT,NA,A,SIG(1,1,iset))
150	XQNS=DBFINT(2,XT,NA,A,QNS(1,1,iset))
151	G =DBFINT(2,XT,NA,A,GL(1,1,iset))
152	*
153	IF (Q2.LT.50.0D0) THEN
154	C Use three flavour evolution.
155	U=(XSIG+9.0D0*XQNS)/6.0D0
156	D=(XSIG-4.5D0*XQNS)/6.0D0
157	S=D
158	C=0.0D0
159	B=0.0D0
160	*
161	ELSE IF((Q2.GT.50.0).AND.(Q2.LT.250.0)) THEN
162	C Use four flavour evolution
163	U=(XSIG+6.0D0*XQNS)/8.0D0
164	D=(XSIG-6.0D0*XQNS)/8.0D0
165	S=D
166	C=U
167	B=0.0D0
168	ELSE
169	C Use five flavour evolution
170	U=(XSIG+7.5D0*XQNS)/10.0D0
171	D=(XSIG-5.0D0*XQNS)/10.0D0
172	S=D
173	C=U
174	B=D
175	ENDIF
176	U=X*U
177	US=U
178	D=X*D
179	DS=D
180	S=X*S
181	C=X*C
182	B=X*B
183	G=X*G
184	90 RETURN
185	END
186	ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
187	SUBROUTINE GS2XCOR(Y,NP)
188	IMPLICIT REAL*8(A-H,O-Z)
189	DIMENSION Y(NP)
190	N=1
191	DO 10 IX=1,20,2
192	Y(N)=DBLE(IX)/2000.0D0
193	N=N+1
194	10 CONTINUE
195	DO 20 IX=30,200,10
196	Y(N)=DBLE(IX)/2000.0D0
197	N=N+1
198	20 CONTINUE
199	DO 30 IX=240,1600,40
200	Y(N)=REAL(IX)/2000.0D0
201	N=N+1
202	30 CONTINUE
203	DO 40 IX=1625,1980,25
204	Y(N)=REAL(IX)/2000.0D0
205	N=N+1
206	40 CONTINUE
207	RETURN
208	END
209	ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc