]>
Commit | Line | Data |
---|---|---|
21886bb6 | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.1.1.2 1996/10/30 08:30:53 cernlib | |
6 | * Version 7.04 | |
7 | * | |
8 | * Revision 1.1.1.1 1996/04/12 15:29:48 plothow | |
9 | * Version 7.01 | |
10 | * | |
11 | * | |
12 | #include "pdf/pilot.h" | |
13 | SUBROUTINE SMRSPI(X,SCALE,MODE,UPV,DNV,SEA,STR,CHM,BOT,GLU) | |
14 | C***************************************************************C | |
15 | C C | |
16 | C PION DISTRIBUTION FUNCTIONS EXTRACTED FROM DRELL-YAN AND C | |
17 | C PROMPT PHOTON PION BEAM DATA, USING (NUCLEAR TARGET C | |
18 | C CORRECTED) HMRS(B) DISTRIBUTIONS FOR THE NUCLEON. C | |
19 | C MAIN UNCERTAINTY IS IN SEA, HENCE FOLLOWING THREE SETS... C | |
20 | C MODE=2 IS THE "BEST FIT" SET. C | |
21 | C C | |
22 | C MODE 1 CORRESPONDS TO A 10% MOMENTUM SEA (0.11213) C | |
23 | C MODE 2 CORRESPONDS TO A 15% MOMENTUM SEA (0.16119) C | |
24 | C MODE 3 CORRESPONDS TO A 20% MOMENTUM SEA (0.23785) C | |
25 | C C | |
26 | C (NUMBER IN BRACKETS IS LAST ENTRY IN FIRST ROW OF GRID) C | |
27 | C C | |
28 | C -*- C | |
29 | C C | |
30 | C REFERENCE: A.D. MARTIN. R.G. ROBERTS. W.J. STIRLING C | |
31 | C AND P.J. SUTTON, PREPRINT RAL-91-058 (1991). C | |
32 | C -*- C | |
33 | C C | |
34 | C FOR THESE SETS.... C | |
35 | C * LAMBDA(MSBAR,NF=4) = 190 MEV C | |
36 | C * UV=DV, STR=USEA=DSEA, BTM=0 C | |
37 | C C | |
38 | C NOTE THAT X TIMES THE PARTON DISTRIBUTION FUNCTION C | |
39 | C IS RETURNED I.E. G(X) = GLU/X ETC, AND THAT "SEA" C | |
40 | C IS THE QUARK SEA / 6 I.E. UBAR(X)=DBAR(X) C | |
41 | C = SEA/X FOR A PION. IF IN DOUBT, CHECK THE C | |
42 | C MOMENTUM SUM RULE! NOTE ALSO THAT SCALE=Q IN GEV C | |
43 | C C | |
44 | C -*- C | |
45 | C C | |
46 | C THE RANGE OF APPLICABILITY IS CURRENTLY: C | |
47 | C 10**-5 < X < 1 AND 5 < Q**2 < 1.31 * 10**6 C | |
48 | C HIGHER Q**2 VALUES CAN BE SUPPLIED ON REQUEST C | |
49 | C - PROBLEMS, COMMENTS ETC TO WJS@UK.AC.DUR.HEP C | |
50 | C C | |
51 | C C | |
52 | C***************************************************************C | |
53 | #include "pdf/impdp.inc" | |
54 | IF(MODE.EQ.1) CALL SMRSP31(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GLU) | |
55 | IF(MODE.EQ.2) CALL SMRSP32(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GLU) | |
56 | IF(MODE.EQ.3) CALL SMRSP33(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GLU) | |
57 | RETURN | |
58 | END |