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