]>
Commit | Line | Data |
---|---|---|
21886bb6 | 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 |