PDF/spdf/newdo1.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.2  1996/10/30 08:32:48  cernlib
   6 * Version 7.04
   7 *
   8 * Revision 1.2  1996/10/28 15:03:18  pdflib
   9 * PDFLIB new version 7.04
  10 *
  11 * Revision 1.1.1.1  1996/04/12 15:29:25  plothow
  12 * Version 7.01
  13 *
  14 *
  15 #include "pdf/pilot.h"
  16               SUBROUTINE NEWDO1(R8X,R8QSQ,R8U,R8D,R8SS,R8STR,R8C,R8G)
  17 C--------------------------------------------------------------
  18 C       THIS SUBROUTINE EVALUATES THE UPDATED DUKE-OWENS QUARK
  19 C       DISTRIBUTIONS , SET 1--------------U. BAUR (JUNE 9, 1991)
  20 C       adapted and corrected by---H.Plothow-Besch (Aug 14, 1991)
  21 C       BASED ON FSU-HEP-910606 PREPRINT
  22 C--------------------------------------------------------------
  23 C
  24 #include "pdf/expdp.inc"
  25      +       R8X,R8QSQ,R8U,R8D,R8SS,R8STR,R8C,R8G
  26         REAL ND,NUD,X,QSQ
  27         DATA Q0/2./, ALAM/0.177/
  28         X = R8X
  29         QSQ = R8QSQ
  30         S= LOG( LOG( MAX(Q0,QSQ)/ALAM)/ LOG(Q0/ALAM))
  31         S2=S*S
  32         S3=S2*S
  33         ONEMX = 1. - X
  34 C-------SEA DIST-----------------------------------------------
  35         AS = 0.9090 - 0.4023*S + 0.006305*S2
  36         AS1 = -0.3823*S + 0.02766*S2
  37         BS1 = 7.278 - 0.7904*S + 0.8108*S2
  38         ALS = -1.6629*S + 0.5719*S2
  39         BES = -0.01333*S + 0.5299*S2
  40         GAS = 0.1211*S - 0.1739*S2
  41 C       SS=AS*X**(AS1-1.)*ONEMX**BS1*(1.+ALS*X+BES*X*X+GAS*X**3)/6.
  42         SS = AS * X**AS1 * ONEMX**BS1 * (1.+ALS*X+BES*X*X+GAS*X*X*X)/ 6.
  43 C-------D-DIST-----------------------------------------
  44         H3 = 0.8388 - 0.2092*S + 0.02657*S2
  45 C       H4 = 4.6670 + 0.7951*S - 0.1081*S2
  46         H4 = 4.6670 + 0.7951*S + 0.1081*S2
  47         GD1 = -1.0232*S + 0.05799*S2
  48         GD2 = 0.8616*S + 0.1530*S2
  49         CALL GAMMADO(H3,G3,IER)
  50         CALL GAMMADO(1.+H3,G31,IER)
  51         CALL GAMMADO(2.+H3,G32,IER)
  52         CALL GAMMADO(1.+H4,G4,IER)
  53         CALL GAMMADO(1.+H3+H4,G34,IER)
  54         CALL GAMMADO(2.+H3+H4,G341,IER)
  55         CALL GAMMADO(3.+H3+H4,G342,IER)
  56         BD1=G3*G4/G34
  57         BD2=G31*G4/G341
  58         BD3=G32*G4/G342
  59         ND=BD1+GD1*BD2+GD2*BD3
  60 C       D=X**(H3-1.)*(1.-X)**H4*(1.+GD1*X+GD2*X**2)/ND
  61         D=X**H3*(1.-X)**H4*(1.+GD1*X+GD2*X**2)/ND
  62 C-------VALENCE U-DIST-----------------------------------------
  63         H1 = 0.6650 - 0.1097*S - 0.002442*S2
  64         H2 = 3.6140 + 0.8395*S - 0.02186*S2
  65         GUD1 = 0.8673 - 1.6637*S + 0.3420*S2
  66 C       GUD2 = 1.1049*S - 0.22369*S2
  67         GUD2 = 1.1049*S - 0.2369*S2
  68         CALL GAMMADO(H1,G1,IER)
  69         CALL GAMMADO(1.+H1,G11,IER)
  70         CALL GAMMADO(2.+H1,G21,IER)
  71         CALL GAMMADO(1.+H2,G2,IER)
  72         CALL GAMMADO(1.+H1+H2,G12,IER)
  73         CALL GAMMADO(2.+H1+H2,G121,IER)
  74         CALL GAMMADO(3.+H1+H2,G122,IER)
  75         BUD1=G1*G2/G12
  76         BUD2=G11*G2/G121
  77         BUD3=G21*G2/G122
  78         NUD=BUD1+GUD1*BUD2+GUD2*BUD3
  79 C       UPD=3*X**(H1-1.)*(1.-X)**H2*(1.+GUD1*X+GUD2*X**2)/NUD
  80         UPD=3. * X**H1 * ONEMX**H2 * (1.+GUD1*X+GUD2*X*X)/ NUD
  81 C       U=UPD-D+SS
  82 C       D=D+SS
  83         U=UPD-D
  84 C-------CHARM DIST---------------------------------------------
  85         AC = 0.09469*S - 0.07066*S2 + 0.01236*S3
  86         AC1 = -0.1447 - 0.4020*S + 0.1533*S2 - 0.06479*S3
  87         BC1 = 6.7599 + 1.6596*S + 0.6798*S2 - 0.8525*S3
  88         ALC = -4.4559*S + 3.3756*S2 - 0.9468*S3
  89         BEC = 7.862*S - 3.6591*S2 + 0.03672*S3
  90         GAC = -0.2472*S - 0.7510*S2 + 0.04870*S3
  91 C       C=AC*X**(AC1-1.)*(1.-X)**BC1*(1.+ALC*X+BEC*X*X+GAC*X**3)
  92         C=AC * X**AC1 * ONEMX**BC1 * (1.+ALC*X+BEC*X*X+GAC*X*X*X)
  93 C-------GLUE DIST----------------------------------------------
  94         AG = 3.0170 - 4.7347*S + 3.3594*S2 - 0.9443*S3
  95         AG1 = -0.9342*S + 0.5454*S2 - 0.1668*S3
  96         BG1 = 5.3040 + 1.4654*S - 1.4292*S2 + 0.7569*S3
  97         ALG = -3.9141*S + 2.8445*S2 - 0.8411*S3
  98         BEG = 9.0176*S - 10.426*S2 + 4.0983*S3
  99         GAG = -5.9602*S + 7.5150*S2 - 2.7329*S3
 100 C       G=AG*X**(AG1-1.)*(1.-X)**BG1*(1.+ALG*X+BEG*X*X+GAG*X**3)
 101         G=AG * X**AG1 * ONEMX**BG1 * (1.+ALG*X+BEG*X*X+GAG*X*X*X)
 102         R8U = U
 103         R8D = D
 104         R8SS = SS
 105         R8STR = SS
 106         R8C = C
 107         R8G = G
 108         RETURN
 109         END