PDF/spdf/valdis.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.2  1996/10/30 08:31:58  cernlib
   6 * Version 7.04
   7 *
   8 * Revision 1.1.1.1  1996/04/12 15:30:07  plothow
   9 * Version 7.01
  10 *
  11 *
  12 #include "pdf/pilot.h"
  13 C----------------------------------------------------------------------
  14       FUNCTION VALDIS(XX)
  15       COMMON/W5051I2/ Q2L
  16       PARAMETER (Q02=0.5,ALAM2=0.045*0.045)
  17       PARAMETER (AL0=0.7, BET0=3.21, ALPRM=-0.19, BETPRM=0.82)
  18
  19       DATA ZERO/0./,HALF/0.5/,ONE/1./,TWO/2./,THRE/3./
  20       DATA IRAZ/0/, Q2OLD/-1./
  21       SAVE AV, ALPHA, BETA
  22
  23       IF(Q2L.NE.Q2OLD) THEN
  24           Q2OLD = Q2L
  25           S     = LOG( LOG( MAX( Q02,Q2L)/ALAM2)/ LOG( Q02/ALAM2))
  26           ALPHA = AL0  + ALPRM *S
  27           BETA  = BET0 + BETPRM*S
  28           AV    = THRE*GAMMA(ALPHA+BETA+ONE)/
  29      1                                 (GAMMA(ALPHA)*GAMMA(BETA+ONE))
  30       ENDIF
  31       VALDIS=AV * XX**ALPHA * (ONE-XX)**BETA
  32       RETURN
  33       END