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

*
* $Id$
*
* $Log$
* Revision 1.1.1.2  1996/10/30 08:28:44  cernlib
* Version 7.04
*
* Revision 1.1.1.1  1996/04/12 15:29:24  plothow
* Version 7.01
*
*
#include "pdf/pilot.h"
       SUBROUTINE GRVPIHO (ZX,ZQ,ZUV,ZDV,ZUDB,ZSB,ZCB,ZBB,ZTB,ZGL)
* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
*                                                                 *
*         G R V - P I 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/16             *
*                                                                 *
*   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.E8 GEV ** 2  AND FOR X BETWEEN  1.E-5  AND  1.   *
*   REGIONS, WHERE THE DISTRIBUTION UNDER CONSIDERATION IS NEG-   *
*   LIGIBLE, I.E. BELOW ABOUT 1.E-4, WERE EXCLUDED FROM THE FIT.  *
*                                                                 *
*              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 DISTRIBUTION REFER TO THE MS-BAR SCHEME OF BARDEEN ET AL.  *
*                                                                 *
* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
C
       IMPLICIT REAL (A - Y)
#include "pdf/expdp.inc"
     +        ZX,ZQ,ZUV,ZDV,ZUDB,ZSB,ZCB,ZBB,ZTB,ZGL
       REAL  X, Q
       X = ZX
       Q = ZQ
       MU2  = 0.3
       LAM2 = 0.248 * 0.248
       Q2 = Q*Q
       S  = ALOG (ALOG(Q2/LAM2) / ALOG(MU2/LAM2))
       DS = SQRT (S)
       S2 = S * S
C...X * VALENCE :
       NV  =  0.456 + 0.150 * DS + 0.112 * S - 0.019 * S2
       AKV =  0.505 - 0.033 * S
       AGV =  0.748 - 0.669 * DS - 0.133 * S
       DV  =  0.365 + 0.197 * DS + 0.394 * S
       VAP =  GRVFVP (X, NV, AKV, AGV, DV)
       ZUV = VAP
       ZDV = ZUV
C...X * GLUON :
       ALG =  1.096
       BEG =  1.371
       AKG =  0.437 - 0.689 * DS
       BKG = -0.631
       AGG =  1.324 - 0.441 * DS - 0.130 * S
       BGG = -0.955 + 0.259 * S
       CG  =  1.075 - 0.302 * S
       DG  =  1.158 + 1.229 * S
       EG  =   0.0  + 2.510 * S
       ESG =  2.604 + 0.165 * S
       GLP =  GRVFGP(X,S, ALG, BEG, AKG, BKG, AGG, BGG, CG, DG, EG, ESG)
       ZGL = GLP
C...X * QBAR (SU(3)-SYMMETRIC SEA) :
       SL  =   0.0
       ALS =   0.85
       BES =   0.96
       AKS = -0.350 + 0.806 * S
       AGS = -1.663
       BS  =  3.148
       DS  =  2.273 + 1.438 * S
       EST =  3.214 + 1.545 * S
       ESS =  1.341 + 1.938 * S
       QBP =  GRVFQBP (X, S, SL, ALS, BES, AKS, AGS, BS, DS, EST, ESS)
       ZUDB = QBP
       ZSB = ZUDB
C...X * CBAR = X * C :
       SC  =  0.820
       ALC =   0.98
       BEC =   0.0
       AKC =   0.0  - 0.457 * S
       AGC =   0.0
       BC  =  -1.00 +  1.40 * S
       DC  =  1.318 + 0.584 * S
       EC  =   4.45 + 1.235 * S
       ESC =  1.496 + 1.010 * S
       CBP =  GRVFQBP (X, S, SC, ALC, BEC, AKC, AGC, BC, DC, EC, ESC)
       ZCB = CBP
C...X * BBAR = X * B :
       SBO =  1.297
       ALB =   0.99
       BEB =   0.0
       AKB =   0.0  - 0.172 * S
       AGB =   0.0
       BBO =   0.0
       DB  =  1.447 + 0.485 * S
       EB  =   4.79 + 1.164 * S
       ESB =  1.724 + 2.121 * S
       BBP =  GRVFQBP (X, S, SBO, ALB, BEB, AKB, AGB, BBO, DB, EB, ESB)
       ZBB = BBP
C...X * TBAR = X * T :
       TBP = 0.
       ZTB = TBP
       RETURN
       END
Commit	Line	Data
fe4da5cc	1	*
	2	* $Id$
	3	*
	4	* $Log$
	5	* Revision 1.1.1.2 1996/10/30 08:28:44 cernlib
	6	* Version 7.04
	7	*
	8	* Revision 1.1.1.1 1996/04/12 15:29:24 plothow
	9	* Version 7.01
	10	*
	11	*
	12	#include "pdf/pilot.h"
	13	SUBROUTINE GRVPIHO (ZX,ZQ,ZUV,ZDV,ZUDB,ZSB,ZCB,ZBB,ZTB,ZGL)
	14	* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
	15	* *
	16	* G R V - P I O N - P A R A M E T R I Z A T I O N S *
	17	* *
	18	* FOR A DETAILED EXPLANATION SEE : *
	19	* M. GLUECK, E.REYA, A.VOGT: DO-TH 91/16 *
	20	* *
	21	* THE PARAMETRIZATIONS ARE FITTED TO THE PARTON DISTRIBUTIONS *
	22	* FOR Q 2 BETWEEN MU 2 (= 0.25 / 0.30 GEV ** 2 IN LO *
	23	* / HO) AND 1.E8 GEV ** 2 AND FOR X BETWEEN 1.E-5 AND 1. *
	24	* REGIONS, WHERE THE DISTRIBUTION UNDER CONSIDERATION IS NEG- *
	25	* LIGIBLE, I.E. BELOW ABOUT 1.E-4, WERE EXCLUDED FROM THE FIT. *
	26	* *
	27	* HEAVY QUARK THRESHOLDS Q(H) = M(H) : *
	28	* M(C) = 1.5, M(B) = 4.5, M(T) = 100 GEV *
	29	* *
	30	* CORRESPONDING LAMBDA(F) VALUES FOR F ACTIVE FLAVOURS : *
	31	* LO : LAMBDA(3) = 0.232, LAMBDA(4) = 0.200, *
	32	* LAMBDA(5) = 0.153, LAMBDA(6) = 0.082 GEV *
	33	* HO : LAMBDA(3) = 0.248, LAMBDA(4) = 0.200, *
	34	* LAMBDA(5) = 0.131, LAMBDA(6) = 0.053 GEV *
	35	* *
	36	* HO DISTRIBUTION REFER TO THE MS-BAR SCHEME OF BARDEEN ET AL. *
	37	* *
	38	* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
	39	C
	40	IMPLICIT REAL (A - Y)
	41	#include "pdf/expdp.inc"
	42	+ ZX,ZQ,ZUV,ZDV,ZUDB,ZSB,ZCB,ZBB,ZTB,ZGL
	43	REAL X, Q
	44	X = ZX
	45	Q = ZQ
	46	MU2 = 0.3
	47	LAM2 = 0.248 * 0.248
	48	Q2 = Q*Q
	49	S = ALOG (ALOG(Q2/LAM2) / ALOG(MU2/LAM2))
	50	DS = SQRT (S)
	51	S2 = S * S
	52	C...X * VALENCE :
	53	NV = 0.456 + 0.150 * DS + 0.112 * S - 0.019 * S2
	54	AKV = 0.505 - 0.033 * S
	55	AGV = 0.748 - 0.669 * DS - 0.133 * S
	56	DV = 0.365 + 0.197 * DS + 0.394 * S
	57	VAP = GRVFVP (X, NV, AKV, AGV, DV)
	58	ZUV = VAP
	59	ZDV = ZUV
	60	C...X * GLUON :
	61	ALG = 1.096
	62	BEG = 1.371
	63	AKG = 0.437 - 0.689 * DS
	64	BKG = -0.631
65	AGG = 1.324 - 0.441 * DS - 0.130 * S
66	BGG = -0.955 + 0.259 * S
67	CG = 1.075 - 0.302 * S
68	DG = 1.158 + 1.229 * S
69	EG = 0.0 + 2.510 * S
70	ESG = 2.604 + 0.165 * S
71	GLP = GRVFGP(X,S, ALG, BEG, AKG, BKG, AGG, BGG, CG, DG, EG, ESG)
72	ZGL = GLP
73	C...X * QBAR (SU(3)-SYMMETRIC SEA) :
74	SL = 0.0
75	ALS = 0.85
76	BES = 0.96
77	AKS = -0.350 + 0.806 * S
78	AGS = -1.663
79	BS = 3.148
80	DS = 2.273 + 1.438 * S
81	EST = 3.214 + 1.545 * S
82	ESS = 1.341 + 1.938 * S
83	QBP = GRVFQBP (X, S, SL, ALS, BES, AKS, AGS, BS, DS, EST, ESS)
84	ZUDB = QBP
85	ZSB = ZUDB
86	C...X * CBAR = X * C :
87	SC = 0.820
88	ALC = 0.98
89	BEC = 0.0
90	AKC = 0.0 - 0.457 * S
91	AGC = 0.0
92	BC = -1.00 + 1.40 * S
93	DC = 1.318 + 0.584 * S
94	EC = 4.45 + 1.235 * S
95	ESC = 1.496 + 1.010 * S
96	CBP = GRVFQBP (X, S, SC, ALC, BEC, AKC, AGC, BC, DC, EC, ESC)
97	ZCB = CBP
98	C...X * BBAR = X * B :
99	SBO = 1.297
100	ALB = 0.99
101	BEB = 0.0
102	AKB = 0.0 - 0.172 * S
103	AGB = 0.0
104	BBO = 0.0
105	DB = 1.447 + 0.485 * S
106	EB = 4.79 + 1.164 * S
107	ESB = 1.724 + 2.121 * S
108	BBP = GRVFQBP (X, S, SBO, ALB, BEB, AKB, AGB, BBO, DB, EB, ESB)
109	ZBB = BBP
110	C...X * TBAR = X * T :
111	TBP = 0.
112	ZTB = TBP
113	RETURN
114	END