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

*
* $Id$
*
* $Log$
* Revision 1.1.1.2  1996/10/30 08:28:40  cernlib
* Version 7.04
*
* Revision 1.2  1996/10/28 15:02:20  pdflib
* PDFLIB new version 7.04
*
* Revision 1.1.1.1  1996/04/12 15:29:23  plothow
* Version 7.01
*
*
#include "pdf/pilot.h"
       SUBROUTINE GRVGALO (ZX,ZQ,ZUV,ZDV,ZUB,ZDB,ZSB,ZCB,ZBB,ZGL)
* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
*                                                                 *
*      G R V - P H O T O N - P A R A M E T R I Z A T I O N S      *
*                                                                 *
*                 FOR A DETAILED EXPLANATION SEE :                *
*              M. GLUECK, E.REYA, A.VOGT: DO-TH 91/31             *
*                                                                 *
*    THE OUTPUT IS ALWAYS   1./ ALPHA(EM) * X * PARTON DENSITY    *
*    output modified by HPB to be always    X * PARTON DENSITY    *
*                                                                 *
*   THE PARAMETRIZATIONS ARE FITTED TO THE PARTON DISTRIBUTIONS   *
*   FOR Q ** 2 BETWEEN MU ** 2 (=  0.25 / 0.30  GEV ** 2  IN LO   *
*   / HO) AND  1.E6 GEV ** 2  AND FOR X BETWEEN  1.E-5  AND  1.   *
*                                                                 *
*              HEAVY QUARK THRESHOLDS  Q(H) = M(H) :              *
*         M(C)  =  1.5,  M(B)  =  4.5,  M(T)  =  100  GEV         *
*                                                                 *
*      CORRESPONDING LAMBDA(F) VALUES FOR F ACTIVE FLAVOURS :     *
*      LO :   LAMBDA(3)  =  0.232,   LAMBDA(4)  =  0.200,         *
*             LAMBDA(5)  =  0.153,   LAMBDA(6)  =  0.082  GEV     *
*      HO :   LAMBDA(3)  =  0.248,   LAMBDA(4)  =  0.200,         *
*             LAMBDA(5)  =  0.131,   LAMBDA(6)  =  0.053  GEV     *
*                                                                 *
*      HO DISTRIBUTIONS REFER TO THE DIS(GAMMA) SCHEME, SEE :     *
*              M. GLUECK, E.REYA, A.VOGT: DO-TH 91/26             *
*                                                                 *
* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
C
       IMPLICIT REAL (A - Y)
#include "pdf/expdp.inc"
     +        ZX,ZQ,ZUV,ZDV,ZUB,ZDB,ZSB,ZCB,ZBB,ZGL
       REAL  X, Q
       DATA ALPHEM/7.29927D-3/
       X = ZX
       Q = ZQ
       MU2  = 0.25
       LAM2 = 0.232 * 0.232
       Q2 = Q*Q
       S  = ALOG (ALOG(Q2/LAM2) / ALOG(MU2/LAM2))
       SS = SQRT (S)
       S2 = S * S
C...X * U = X * UBAR :
       AL =  1.717
       BE =  0.641
       AK =  0.500 - 0.176 * S
       BK = 15.00  - 5.687 * SS - 0.552 * S2
       AG =  0.235 + 0.046 * SS
       BG =  0.082 - 0.051 * S  + 0.168 * S2
       C  =   0.0  + 0.459 * S
       D  =  0.354 - 0.061 * S
       E  =  4.899 + 1.678 * S
       ES =  2.046 + 1.389 * S
       UL =  GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
       ZUV = UL * ALPHEM
       ZUB = ZUV
C...X * D = X * DBAR :
       AL =  1.549
       BE =  0.782
       AK =  0.496 + 0.026 * S
       BK =  0.685 - 0.580 * SS + 0.608 * S2
       AG =  0.233 + 0.302 * S
       BG =   0.0  - 0.818 * S  + 0.198 * S2
       C  =  0.114 + 0.154 * S
       D  =  0.405 - 0.195 * S  + 0.046 * S2
       E  =  4.807 + 1.226 * S
       ES =  2.166 + 0.664 * S
       DL  =  GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
       ZDV = DL * ALPHEM
       ZDB = ZDV
C...X * G :
       AL =  0.676
       BE =  1.089
       AK =  0.462 - 0.524 * SS
       BK =  5.451              - 0.804 * S2
       AG =  0.535 - 0.504 * SS + 0.288 * S2
       BG =  0.364 - 0.520 * S
       C  = -0.323              + 0.115 * S2
       D  =  0.233 + 0.790 * S  - 0.139 * S2
       E  =  0.893 + 1.968 * S
       ES =  3.432 + 0.392 * S
       GL =  GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
       ZGL = GL * ALPHEM
C...X * S = X * SBAR :
       SF =   0.0
       AL =  1.609
       BE =  0.962
       AK =  0.470              - 0.099 * S2
       BK =  3.246
       AG =  0.121 - 0.068 * SS
       BG = -0.090 + 0.074 * S
       C  =  0.062 + 0.034 * S
       D  =   0.0  + 0.226 * S  - 0.060 * S2
       E  =  4.288 + 1.707 * S
       ES =  2.122 + 0.656 * S
       SL =  GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES)
       ZSB = SL * ALPHEM
C...X * C = X * CBAR :
       SF =  0.888
       AL =  0.970
       BE =  0.545
       AK =  1.254 - 0.251 * S
       BK =  3.932              - 0.327 * S2
       AG =  0.658 + 0.202 * S
       BG = -0.699
       C  =  0.965
       D  =   0.0  + 0.141 * S  - 0.027 * S2
       E  =  4.911 + 0.969 * S
       ES =  2.796 + 0.952 * S
       CL =  GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES)
       ZCB = CL * ALPHEM
C...X * B = X * BBAR :
       SF =  1.351
       AL =  1.016
       BE =  0.338
       AK =  1.961 - 0.370 * S
       BK =  0.923 + 0.119 * S
       AG =  0.815 + 0.207 * S
       BG = -2.275
       C  =  1.480
       D  = -0.223 + 0.173 * S
       E  =  5.426 + 0.623 * S
       ES =  3.819 + 0.901 * S
       BL =  GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES)
       ZBB = BL * ALPHEM
C
       RETURN
       END
Commit	Line	Data
fe4da5cc	1	*
	2	* $Id$
	3	*
	4	* $Log$
	5	* Revision 1.1.1.2 1996/10/30 08:28:40 cernlib
	6	* Version 7.04
	7	*
	8	* Revision 1.2 1996/10/28 15:02:20 pdflib
	9	* PDFLIB new version 7.04
	10	*
	11	* Revision 1.1.1.1 1996/04/12 15:29:23 plothow
	12	* Version 7.01
	13	*
	14	*
	15	#include "pdf/pilot.h"
	16	SUBROUTINE GRVGALO (ZX,ZQ,ZUV,ZDV,ZUB,ZDB,ZSB,ZCB,ZBB,ZGL)
	17	* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
	18	* *
	19	* G R V - P H O T O N - P A R A M E T R I Z A T I O N S *
	20	* *
	21	* FOR A DETAILED EXPLANATION SEE : *
	22	* M. GLUECK, E.REYA, A.VOGT: DO-TH 91/31 *
	23	* *
	24	* THE OUTPUT IS ALWAYS 1./ ALPHA(EM) * X * PARTON DENSITY *
	25	* output modified by HPB to be always X * PARTON DENSITY *
	26	* *
	27	* THE PARAMETRIZATIONS ARE FITTED TO THE PARTON DISTRIBUTIONS *
	28	* FOR Q 2 BETWEEN MU 2 (= 0.25 / 0.30 GEV ** 2 IN LO *
	29	* / HO) AND 1.E6 GEV ** 2 AND FOR X BETWEEN 1.E-5 AND 1. *
	30	* *
	31	* HEAVY QUARK THRESHOLDS Q(H) = M(H) : *
	32	* M(C) = 1.5, M(B) = 4.5, M(T) = 100 GEV *
	33	* *
	34	* CORRESPONDING LAMBDA(F) VALUES FOR F ACTIVE FLAVOURS : *
	35	* LO : LAMBDA(3) = 0.232, LAMBDA(4) = 0.200, *
	36	* LAMBDA(5) = 0.153, LAMBDA(6) = 0.082 GEV *
	37	* HO : LAMBDA(3) = 0.248, LAMBDA(4) = 0.200, *
	38	* LAMBDA(5) = 0.131, LAMBDA(6) = 0.053 GEV *
	39	* *
	40	* HO DISTRIBUTIONS REFER TO THE DIS(GAMMA) SCHEME, SEE : *
	41	* M. GLUECK, E.REYA, A.VOGT: DO-TH 91/26 *
	42	* *
	43	* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
	44	C
	45	IMPLICIT REAL (A - Y)
	46	#include "pdf/expdp.inc"
	47	+ ZX,ZQ,ZUV,ZDV,ZUB,ZDB,ZSB,ZCB,ZBB,ZGL
	48	REAL X, Q
	49	DATA ALPHEM/7.29927D-3/
	50	X = ZX
	51	Q = ZQ
	52	MU2 = 0.25
	53	LAM2 = 0.232 * 0.232
	54	Q2 = Q*Q
	55	S = ALOG (ALOG(Q2/LAM2) / ALOG(MU2/LAM2))
	56	SS = SQRT (S)
	57	S2 = S * S
	58	C...X * U = X * UBAR :
	59	AL = 1.717
	60	BE = 0.641
	61	AK = 0.500 - 0.176 * S
	62	BK = 15.00 - 5.687 * SS - 0.552 * S2
	63	AG = 0.235 + 0.046 * SS
	64	BG = 0.082 - 0.051 * S + 0.168 * S2
65	C = 0.0 + 0.459 * S
66	D = 0.354 - 0.061 * S
67	E = 4.899 + 1.678 * S
68	ES = 2.046 + 1.389 * S
69	UL = GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
70	ZUV = UL * ALPHEM
71	ZUB = ZUV
72	C...X * D = X * DBAR :
73	AL = 1.549
74	BE = 0.782
75	AK = 0.496 + 0.026 * S
76	BK = 0.685 - 0.580 * SS + 0.608 * S2
77	AG = 0.233 + 0.302 * S
78	BG = 0.0 - 0.818 * S + 0.198 * S2
79	C = 0.114 + 0.154 * S
80	D = 0.405 - 0.195 * S + 0.046 * S2
81	E = 4.807 + 1.226 * S
82	ES = 2.166 + 0.664 * S
83	DL = GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
84	ZDV = DL * ALPHEM
85	ZDB = ZDV
86	C...X * G :
87	AL = 0.676
88	BE = 1.089
89	AK = 0.462 - 0.524 * SS
90	BK = 5.451 - 0.804 * S2
91	AG = 0.535 - 0.504 * SS + 0.288 * S2
92	BG = 0.364 - 0.520 * S
93	C = -0.323 + 0.115 * S2
94	D = 0.233 + 0.790 * S - 0.139 * S2
95	E = 0.893 + 1.968 * S
96	ES = 3.432 + 0.392 * S
97	GL = GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
98	ZGL = GL * ALPHEM
99	C...X * S = X * SBAR :
100	SF = 0.0
101	AL = 1.609
102	BE = 0.962
103	AK = 0.470 - 0.099 * S2
104	BK = 3.246
105	AG = 0.121 - 0.068 * SS
106	BG = -0.090 + 0.074 * S
107	C = 0.062 + 0.034 * S
108	D = 0.0 + 0.226 * S - 0.060 * S2
109	E = 4.288 + 1.707 * S
110	ES = 2.122 + 0.656 * S
111	SL = GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES)
112	ZSB = SL * ALPHEM
113	C...X * C = X * CBAR :
114	SF = 0.888
115	AL = 0.970
116	BE = 0.545
117	AK = 1.254 - 0.251 * S
118	BK = 3.932 - 0.327 * S2
119	AG = 0.658 + 0.202 * S
120	BG = -0.699
121	C = 0.965
122	D = 0.0 + 0.141 * S - 0.027 * S2
123	E = 4.911 + 0.969 * S
124	ES = 2.796 + 0.952 * S
125	CL = GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES)
126	ZCB = CL * ALPHEM
127	C...X * B = X * BBAR :
128	SF = 1.351
129	AL = 1.016
130	BE = 0.338
131	AK = 1.961 - 0.370 * S
132	BK = 0.923 + 0.119 * S
133	AG = 0.815 + 0.207 * S
134	BG = -2.275
135	C = 1.480
136	D = -0.223 + 0.173 * S
137	E = 5.426 + 0.623 * S
138	ES = 3.819 + 0.901 * S
139	BL = GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES)
140	ZBB = BL * ALPHEM
141	C
142	RETURN
143	END