Multi-Parton Program

Here we provide code for a somewhat more elaborate program for interpolating parton distributions from tables. The tables are in files, such as cteq3m.ptn or mrsap.ptn. The user's program must first call the subroutine OPENPARTON, as in
        CALL OPENPARTON(cteq3m.ptn,1)
        CALL OPENPARTON(cteq3d.ptn,2)
        CALL OPENPARTON(cteq5m.ptn,3)
        CALL OPENPARTON(cteq5d.ptn,4)
Then calls to the function
PARTON(NPARTON,X,MU,J)
will produce the parton distribution for partons of type NPARTON from set J at momentum fraction X and scale MU.
Beware, an early version of this code had a BUG .


Davison E. Soper and Parvez Anandam
Institute of Theoretical Science
University of Oregon
Eugene OR 97403. USA.
soper@bovine.uoregon.edu
email: soper@bovine.uoregon.edu