Back to home page

sPhenix code displayed by LXR

 
 

    


File indexing completed on 2025-08-05 08:21:11

0001  
0002 C*********************************************************************
0003  
0004 C...PYGRVS
0005 C...Auxiliary for the GRV 94 parton distribution functions
0006 C...for s, c and b sea.
0007 C...Authors: M. Glueck, E. Reya and A. Vogt.
0008  
0009       FUNCTION PYGRVS (X, S, STH, AL, BE, AK, AG, B, D, E, ES)
0010  
0011 C...Double precision declaration.
0012       IMPLICIT DOUBLE PRECISION (A - Z)
0013  
0014 C...Evaluation.
0015       IF(S.LE.STH) THEN
0016         PYGRVS = 0D0
0017       ELSE
0018         DX = SQRT (X)
0019         LX = LOG (1D0/X)
0020         PYGRVS = (S - STH)**AL / LX**AK * (1D0+ AG*DX + B*X) *
0021      &     (1D0- X)**D * EXP (-E + SQRT (ES * S**BE * LX))
0022       ENDIF
0023  
0024       RETURN
0025       END