]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.1.1.2 1996/10/30 08:27:50 cernlib | |
6 | * Version 7.04 | |
7 | * | |
8 | * Revision 1.1.1.1 1996/04/12 15:29:10 plothow | |
9 | * Version 7.01 | |
10 | * | |
11 | * | |
12 | #include "pdf/pilot.h" | |
13 | C | |
14 | C | |
15 | SUBROUTINE Ctq3Pds(Iset, Pdf, X, Q, Irt) | |
16 | ||
17 | C This function returns the CTEQ parton distributions xf^Iset_Iprtn/proton | |
18 | C --- the Momentum density in array form | |
19 | c | |
20 | C (Iset, X, Q): explained in header comment lines; | |
21 | ||
22 | C Irt : return error code -- cumulated over flavors: | |
23 | C see module Ctq3Pd for explanation on individual flavors. | |
24 | C Pdf (Iparton); | |
25 | C Iparton = -6, -5, ...0, 1, 2 ... 6 | |
26 | C has the same meaning as explained in the header comment lines. | |
27 | ||
28 | C Implicit Double Precision (A-H, O-Z) | |
29 | C+SEQ, IMPDP. | |
30 | Dimension Pdf (-6:6) | |
31 | ||
32 | Irt=0 | |
33 | do 10 I=-6,2 | |
34 | if(I.le.0) then | |
35 | Pdf(I) = Ctq3Pd(Iset,I,X,Q,Irt1) | |
36 | Pdf(-I)= Pdf(I) | |
37 | else | |
38 | Pdf(I) = Ctq3Pd(Iset,I,X,Q,Irt1) + Pdf(-I) | |
39 | endif | |
40 | Irt=Irt+Irt1 | |
41 | 10 Continue | |
42 | ||
43 | Return | |
44 | C ************************* | |
45 | End |