[u/mrichter/AliRoot.git] / PDF / spdf / rxqbar.F

*
* $Id$
*
* $Log$
* Revision 1.1.1.2  1996/10/30 08:29:17  cernlib
* Version 7.04
*
* Revision 1.1.1.1  1996/04/12 15:29:29  plothow
* Version 7.01
*
*
#include "pdf/pilot.h"
      SUBROUTINE RXQBAR (X,QQ,XQBAR)
C.......................................................................
C.
C.    SEA
C.
C.......................................................................
#include "pdf/impdp.inc"
      DATA XLAM /0.50D00/, Q0/2.1D0/, Q00/1.8D0/
      DATA ONED/1.D0/, THRED/3.D0/
      IF(QQ.LT.Q0) QQ=Q0
      S=LOG(QQ/XLAM**2)/LOG(Q00/XLAM**2)
      S=LOG(S)
      AS=0.1467D00-0.1211D00*S+0.0274D00*S**2
      ETAS=7.0D00+0.0217D00*S+0.0037D00*S**2
      APS=0.1853D00*S-0.0608D00*S**2
      ETAPS=9.5041D00+1.0165D00*S-0.1049D00*S**2
      BS=0.1682D00*S+0.4473D00*S**2
      CS=25.8997D00+3.9572D00*S+1.6331D00*S**2
      XQBAR=AS*((ONED-X)**ETAS)+APS*((ONED-X)**ETAPS)+BS*EXP(-CS*X)
      RETURN
      END
Commit	Line	Data
21886bb6	1	*
	2	* $Id$
	3	*
	4	* $Log$
	5	* Revision 1.1.1.2 1996/10/30 08:29:17 cernlib
	6	* Version 7.04
	7	*
	8	* Revision 1.1.1.1 1996/04/12 15:29:29 plothow
	9	* Version 7.01
	10	*
	11	*
	12	#include "pdf/pilot.h"
	13	SUBROUTINE RXQBAR (X,QQ,XQBAR)
	14	C.......................................................................
	15	C.
	16	C. SEA
	17	C.
	18	C.......................................................................
	19	#include "pdf/impdp.inc"
	20	DATA XLAM /0.50D00/, Q0/2.1D0/, Q00/1.8D0/
	21	DATA ONED/1.D0/, THRED/3.D0/
	22	IF(QQ.LT.Q0) QQ=Q0
	23	S=LOG(QQ/XLAM2)/LOG(Q00/XLAM2)
	24	S=LOG(S)
	25	AS=0.1467D00-0.1211D00S+0.0274D00S**2
	26	ETAS=7.0D00+0.0217D00S+0.0037D00S**2
	27	APS=0.1853D00S-0.0608D00S**2
	28	ETAPS=9.5041D00+1.0165D00S-0.1049D00S**2
	29	BS=0.1682D00S+0.4473D00S**2
	30	CS=25.8997D00+3.9572D00S+1.6331D00S**2
	31	XQBAR=AS((ONED-X)ETAS)+APS((ONED-X)*ETAPS)+BSEXP(-CS*X)
	32	RETURN
	33	END