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

*
* $Id$
*
* $Log$
* Revision 1.1.1.2  1996/10/30 08:28:35  cernlib
* Version 7.04
*
* Revision 1.1.1.1  1996/04/12 15:29:22  plothow
* Version 7.01
*
*
#include "pdf/pilot.h"
*
*...NLO PARAMETRIZATION (MS(BAR)) :
*
       SUBROUTINE GRV94HO (ZX,ZQ,ZUV,ZDV,ZUB,ZDB,ZSB,ZCB,ZBB,ZTB,ZGL)
       IMPLICIT REAL (A - Y)
C+SEQ, IMPDP.
#include "pdf/expdp.inc"
     +      ZX,ZQ,ZUV,ZDV,ZDEL,ZUDB,ZUB,ZDB,ZSB,ZCB,ZBB,ZTB,ZGL
       REAL X, Q
       X = ZX
       Q = ZQ
       MU2  = 0.34
       LAM2 = 0.248   * 0.248
       Q2 = Q*Q
       ARG1 = MAX(MU2,Q2)/LAM2
       ARG2 = MU2/LAM2
       S  =  LOG ( LOG( ARG1  ) /  LOG( ARG2   ))
       DS =  SQRT (S)
       S2 = S * S
       S3 = S2 * S
*...UV :
       NU  =  1.304 + 0.863 * S
       AKU =  0.558 - 0.020 * S
       BKU =          0.183 * S
       AU  = -0.113 + 0.283 * S - 0.321 * S2
       BU  =  6.843 - 5.089 * S + 2.647 * S2 - 0.527 * S3
       CU  =  7.771 - 10.09 * S + 2.630 * S2
       DU  =  3.315 + 1.145 * S - 0.583 * S2 + 0.154 * S3
       UV  = GR94FV (X, NU, AKU, BKU, AU, BU, CU, DU)
       ZUV = UV
*...DV :
       ND  =  0.102 - 0.017 * S + 0.005 * S2
       AKD =  0.270 - 0.019 * S
       BKD =  0.260
       AD  =  2.393 + 6.228 * S - 0.881 * S2
       BD  =  46.06 + 4.673 * S - 14.98 * S2 + 1.331 * S3
       CD  =  17.83 - 53.47 * S + 21.24 * S2
       DD  =  4.081 + 0.976 * S - 0.485 * S2 + 0.152 * S3
       DV  = GR94FV (X, ND, AKD, BKD, AD, BD, CD, DD)
       ZDV = DV
*...DEL :
       NE  =  0.070 + 0.042 * S - 0.011 * S2 + 0.004 * S3
       AKE =  0.409 - 0.007 * S
       BKE =  0.782 + 0.082 * S
       AE  = -29.65 + 26.49 * S + 5.429 * S2
       BE  =  90.20 - 74.97 * S + 4.526 * S2
       CE  =  0.0
       DE  =  8.122 + 2.120 * S - 1.088 * S2 + 0.231 * S3
       DEL = GR94FV (X, NE, AKE, BKE, AE, BE, CE, DE)
       ZDEL= DEL
*...UDB :
       ALX =  0.877
       BEX =  0.561
       AKX =  0.275
       BKX =  0.0
       AGX =  0.997
       BGX =  3.210 - 1.866 * S
       CX  =  7.300
       DX  =  9.010 + 0.896 * DS + 0.222 * S2
       EX  =  3.077 + 1.446 * S
       ESX =  3.173 - 2.445 * DS + 2.207 * S
       UDB = GR94FW(X, S, ALX, BEX, AKX, BKX, AGX, BGX, CX, DX, EX, ESX)
       ZUDB= UDB
*...UB
       ZUB = (ZUDB - ZDEL)/2.D0
*...DB
       ZDB = (ZUDB + ZDEL)/2.D0
*...SB :
       ALS =  0.756
       BES =  0.216
       AKS =  1.690 + 0.650 * DS - 0.922 * S
       AS  = -4.329 + 1.131 * S
       BS  =  9.568 - 1.744 * S
       DST =  9.377 + 1.088 * DS - 1.320 * S + 0.130 * S2
       EST =  3.031 + 1.639 * S
       ESS =  5.837 + 0.815 * S
       SB  = GR94FWS (X, S, ALS, BES, AKS, AS, BS, DST, EST, ESS)
       ZSB = SB
*...GL :
       ALG =  1.014
       BEG =  1.738
       AKG =  1.724 + 0.157 * S
       BKG =  0.800 + 1.016 * S
       AG  =  7.517 - 2.547 * S
       BG  =  34.09 - 52.21 * DS + 17.47 * S
       CG  =  4.039 + 1.491 * S
       DG  =  3.404 + 0.830 * S
       EG  = -1.112 + 3.438 * S  - 0.302 * S2
       ESG =  3.256 - 0.436 * S
       GL  = GR94FW (X, S, ALG, BEG, AKG, BKG, AG, BG, CG, DG, EG, ESG)
       ZGL = GL
*...
       ZCB = 0.D0
       ZBB = 0.D0
       ZTB = 0.D0
*...
       RETURN
       END
Commit	Line	Data
21886bb6	1	*
	2	* $Id$
	3	*
	4	* $Log$
	5	* Revision 1.1.1.2 1996/10/30 08:28:35 cernlib
	6	* Version 7.04
	7	*
	8	* Revision 1.1.1.1 1996/04/12 15:29:22 plothow
	9	* Version 7.01
	10	*
	11	*
	12	#include "pdf/pilot.h"
	13	*
	14	*...NLO PARAMETRIZATION (MS(BAR)) :
	15	*
	16	SUBROUTINE GRV94HO (ZX,ZQ,ZUV,ZDV,ZUB,ZDB,ZSB,ZCB,ZBB,ZTB,ZGL)
	17	IMPLICIT REAL (A - Y)
	18	C+SEQ, IMPDP.
	19	#include "pdf/expdp.inc"
	20	+ ZX,ZQ,ZUV,ZDV,ZDEL,ZUDB,ZUB,ZDB,ZSB,ZCB,ZBB,ZTB,ZGL
	21	REAL X, Q
	22	X = ZX
	23	Q = ZQ
	24	MU2 = 0.34
	25	LAM2 = 0.248 * 0.248
	26	Q2 = Q*Q
	27	ARG1 = MAX(MU2,Q2)/LAM2
	28	ARG2 = MU2/LAM2
	29	S = LOG ( LOG( ARG1 ) / LOG( ARG2 ))
	30	DS = SQRT (S)
	31	S2 = S * S
	32	S3 = S2 * S
	33	*...UV :
	34	NU = 1.304 + 0.863 * S
	35	AKU = 0.558 - 0.020 * S
	36	BKU = 0.183 * S
	37	AU = -0.113 + 0.283 * S - 0.321 * S2
	38	BU = 6.843 - 5.089 * S + 2.647 * S2 - 0.527 * S3
	39	CU = 7.771 - 10.09 * S + 2.630 * S2
	40	DU = 3.315 + 1.145 * S - 0.583 * S2 + 0.154 * S3
	41	UV = GR94FV (X, NU, AKU, BKU, AU, BU, CU, DU)
	42	ZUV = UV
	43	*...DV :
	44	ND = 0.102 - 0.017 * S + 0.005 * S2
	45	AKD = 0.270 - 0.019 * S
	46	BKD = 0.260
	47	AD = 2.393 + 6.228 * S - 0.881 * S2
	48	BD = 46.06 + 4.673 * S - 14.98 * S2 + 1.331 * S3
	49	CD = 17.83 - 53.47 * S + 21.24 * S2
	50	DD = 4.081 + 0.976 * S - 0.485 * S2 + 0.152 * S3
	51	DV = GR94FV (X, ND, AKD, BKD, AD, BD, CD, DD)
	52	ZDV = DV
	53	*...DEL :
	54	NE = 0.070 + 0.042 * S - 0.011 * S2 + 0.004 * S3
	55	AKE = 0.409 - 0.007 * S
	56	BKE = 0.782 + 0.082 * S
	57	AE = -29.65 + 26.49 * S + 5.429 * S2
	58	BE = 90.20 - 74.97 * S + 4.526 * S2
	59	CE = 0.0
	60	DE = 8.122 + 2.120 * S - 1.088 * S2 + 0.231 * S3
	61	DEL = GR94FV (X, NE, AKE, BKE, AE, BE, CE, DE)
	62	ZDEL= DEL
	63	*...UDB :
	64	ALX = 0.877
65	BEX = 0.561
66	AKX = 0.275
67	BKX = 0.0
68	AGX = 0.997
69	BGX = 3.210 - 1.866 * S
70	CX = 7.300
71	DX = 9.010 + 0.896 * DS + 0.222 * S2
72	EX = 3.077 + 1.446 * S
73	ESX = 3.173 - 2.445 * DS + 2.207 * S
74	UDB = GR94FW(X, S, ALX, BEX, AKX, BKX, AGX, BGX, CX, DX, EX, ESX)
75	ZUDB= UDB
76	*...UB
77	ZUB = (ZUDB - ZDEL)/2.D0
78	*...DB
79	ZDB = (ZUDB + ZDEL)/2.D0
80	*...SB :
81	ALS = 0.756
82	BES = 0.216
83	AKS = 1.690 + 0.650 * DS - 0.922 * S
84	AS = -4.329 + 1.131 * S
85	BS = 9.568 - 1.744 * S
86	DST = 9.377 + 1.088 * DS - 1.320 * S + 0.130 * S2
87	EST = 3.031 + 1.639 * S
88	ESS = 5.837 + 0.815 * S
89	SB = GR94FWS (X, S, ALS, BES, AKS, AS, BS, DST, EST, ESS)
90	ZSB = SB
91	*...GL :
92	ALG = 1.014
93	BEG = 1.738
94	AKG = 1.724 + 0.157 * S
95	BKG = 0.800 + 1.016 * S
96	AG = 7.517 - 2.547 * S
97	BG = 34.09 - 52.21 * DS + 17.47 * S
98	CG = 4.039 + 1.491 * S
99	DG = 3.404 + 0.830 * S
100	EG = -1.112 + 3.438 * S - 0.302 * S2
101	ESG = 3.256 - 0.436 * S
102	GL = GR94FW (X, S, ALG, BEG, AKG, BKG, AG, BG, CG, DG, EG, ESG)
103	ZGL = GL
104	*...
105	ZCB = 0.D0
106	ZBB = 0.D0
107	ZTB = 0.D0
108	*...
109	RETURN
110	END