Old version of PDF removed
authorhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 18 Sep 2000 09:12:41 +0000 (09:12 +0000)
committerhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 18 Sep 2000 09:12:41 +0000 (09:12 +0000)
344 files changed:
PDF/Makefile [deleted file]
PDF/doc/inform.dat [deleted file]
PDF/doc/pdfdoc.dat [deleted file]
PDF/doc/v100.dat [deleted file]
PDF/dpdf/inform.dat [deleted file]
PDF/dpdf/pdfdoc.doc [deleted file]
PDF/dpdf/v100.dat [deleted file]
PDF/pdf/expdp.inc [deleted file]
PDF/pdf/impdp.inc [deleted file]
PDF/pdf/pilot.h [deleted file]
PDF/pdf/w50510.inc [deleted file]
PDF/pdf/w50511.inc [deleted file]
PDF/pdf/w505110.inc [deleted file]
PDF/pdf/w50511c.inc [deleted file]
PDF/pdf/w50512.inc [deleted file]
PDF/pdf/w505120.inc [deleted file]
PDF/pdf/w505121.inc [deleted file]
PDF/pdf/w505122.inc [deleted file]
PDF/pdf/w50513.inc [deleted file]
PDF/pdf/w50514.inc [deleted file]
PDF/pdf/w50514w.inc [deleted file]
PDF/pdf/w50515.inc [deleted file]
PDF/pdf/w50516.inc [deleted file]
PDF/pdf/w50517.inc [deleted file]
PDF/pdf/w50518.inc [deleted file]
PDF/pdf/w50519.inc [deleted file]
PDF/pdf/w5051p0.inc [deleted file]
PDF/pdf/w5051p1.inc [deleted file]
PDF/pdf/w5051p2.inc [deleted file]
PDF/pdf/w5051p3.inc [deleted file]
PDF/pdf/w5051p7.inc [deleted file]
PDF/spdf/abfkw1.F [deleted file]
PDF/spdf/abfkw2.F [deleted file]
PDF/spdf/abfkw3.F [deleted file]
PDF/spdf/acfgp1.F [deleted file]
PDF/spdf/acfgp2.F [deleted file]
PDF/spdf/afcplu.F [deleted file]
PDF/spdf/afgetfv.F [deleted file]
PDF/spdf/afgin2.F [deleted file]
PDF/spdf/afgint.F [deleted file]
PDF/spdf/afpolin.F [deleted file]
PDF/spdf/afratin.F [deleted file]
PDF/spdf/alphas2.F [deleted file]
PDF/spdf/augetfv.F [deleted file]
PDF/spdf/aurgam.F [deleted file]
PDF/spdf/aurpi1.F [deleted file]
PDF/spdf/aurpi2.F [deleted file]
PDF/spdf/aurpi3.F [deleted file]
PDF/spdf/aurpr.F [deleted file]
PDF/spdf/bbeta.F [deleted file]
PDF/spdf/beta.F [deleted file]
PDF/spdf/bxf3.F [deleted file]
PDF/spdf/bxglue.F [deleted file]
PDF/spdf/bxqbar.F [deleted file]
PDF/spdf/ctq1opd.F [deleted file]
PDF/spdf/ctq1opf.F [deleted file]
PDF/spdf/ctq1pd.F [deleted file]
PDF/spdf/ctq1pf.F [deleted file]
PDF/spdf/ctq2opd.F [deleted file]
PDF/spdf/ctq2opf.F [deleted file]
PDF/spdf/ctq2ops.F [deleted file]
PDF/spdf/ctq2ppd.F [deleted file]
PDF/spdf/ctq2ppf.F [deleted file]
PDF/spdf/ctq2pps.F [deleted file]
PDF/spdf/ctq3pd.F [deleted file]
PDF/spdf/ctq3pds.F [deleted file]
PDF/spdf/ctq3pf.F [deleted file]
PDF/spdf/ctq4Fn.F [deleted file]
PDF/spdf/ctq4PnX.F [deleted file]
PDF/spdf/ctq4RA1.F [deleted file]
PDF/spdf/ctq4RA2.F [deleted file]
PDF/spdf/ctq4RA4.F [deleted file]
PDF/spdf/ctq4RA5.F [deleted file]
PDF/spdf/ctq4RD.F [deleted file]
PDF/spdf/ctq4RHJ.F [deleted file]
PDF/spdf/ctq4RL.F [deleted file]
PDF/spdf/ctq4RM.F [deleted file]
PDF/spdf/ctq4RlQ.F [deleted file]
PDF/spdf/dbfint.F [deleted file]
PDF/spdf/dflm1.F [deleted file]
PDF/spdf/dflm2.F [deleted file]
PDF/spdf/dflm3.F [deleted file]
PDF/spdf/dflm4.F [deleted file]
PDF/spdf/dflm5.F [deleted file]
PDF/spdf/dflm6.F [deleted file]
PDF/spdf/dflm7.F [deleted file]
PDF/spdf/dflm8.F [deleted file]
PDF/spdf/dflm9.F [deleted file]
PDF/spdf/dgpho1.F [deleted file]
PDF/spdf/dgpho2.F [deleted file]
PDF/spdf/dgpho3.F [deleted file]
PDF/spdf/dgpho4.F [deleted file]
PDF/spdf/dopho1.F [deleted file]
PDF/spdf/dopho2.F [deleted file]
PDF/spdf/dpolin.F [deleted file]
PDF/spdf/fxaver.F [deleted file]
PDF/spdf/fxg160.F [deleted file]
PDF/spdf/fxg260.F [deleted file]
PDF/spdf/fxg360.F [deleted file]
PDF/spdf/fxghar.F [deleted file]
PDF/spdf/fxgsof.F [deleted file]
PDF/spdf/fxnlla.F [deleted file]
PDF/spdf/fxvhar.F [deleted file]
PDF/spdf/fxvsof.F [deleted file]
PDF/spdf/gamfun.F [deleted file]
PDF/spdf/gammado.F [deleted file]
PDF/spdf/ghrdv.F [deleted file]
PDF/spdf/ghrgl.F [deleted file]
PDF/spdf/ghruv.F [deleted file]
PDF/spdf/ghrxi.F [deleted file]
PDF/spdf/ghrxs.F [deleted file]
PDF/spdf/gr94fv.F [deleted file]
PDF/spdf/gr94fw.F [deleted file]
PDF/spdf/gr94fws.F [deleted file]
PDF/spdf/grsf1.F [deleted file]
PDF/spdf/grsf2.F [deleted file]
PDF/spdf/grsgalo.F [deleted file]
PDF/spdf/grsvh1.F [deleted file]
PDF/spdf/grsvh2.F [deleted file]
PDF/spdf/grsvl1.F [deleted file]
PDF/spdf/grsvl2.F [deleted file]
PDF/spdf/grv94di.F [deleted file]
PDF/spdf/grv94ho.F [deleted file]
PDF/spdf/grv94lo.F [deleted file]
PDF/spdf/grvfgp.F [deleted file]
PDF/spdf/grvfqbp.F [deleted file]
PDF/spdf/grvfv.F [deleted file]
PDF/spdf/grvfvp.F [deleted file]
PDF/spdf/grvfw.F [deleted file]
PDF/spdf/grvfws.F [deleted file]
PDF/spdf/grvgah0.F [deleted file]
PDF/spdf/grvgaho.F [deleted file]
PDF/spdf/grvgalo.F [deleted file]
PDF/spdf/grvgf.F [deleted file]
PDF/spdf/grvgfs.F [deleted file]
PDF/spdf/grvho.F [deleted file]
PDF/spdf/grvlo.F [deleted file]
PDF/spdf/grvpiho.F [deleted file]
PDF/spdf/grvpilo.F [deleted file]
PDF/spdf/gs2xcor.F [deleted file]
PDF/spdf/gs96ho.F [deleted file]
PDF/spdf/gs96lo.F [deleted file]
PDF/spdf/gsxcor.F [deleted file]
PDF/spdf/hmrs1eb.F [deleted file]
PDF/spdf/hmrs2eb.F [deleted file]
PDF/spdf/hmrs3eb.F [deleted file]
PDF/spdf/hmrsgeb.F [deleted file]
PDF/spdf/kmrseb.F [deleted file]
PDF/spdf/kmrsgeb.F [deleted file]
PDF/spdf/mrsdshd.F [deleted file]
PDF/spdf/mrsdshp.F [deleted file]
PDF/spdf/mrseb.F [deleted file]
PDF/spdf/mrsebp.F [deleted file]
PDF/spdf/newdo1.F [deleted file]
PDF/spdf/pdfset.F [deleted file]
PDF/spdf/pdfsta.F [deleted file]
PDF/spdf/pdfvers.F [deleted file]
PDF/spdf/pdxmt.F [deleted file]
PDF/spdf/pdzxmt.F [deleted file]
PDF/spdf/pftopdg.F [deleted file]
PDF/spdf/phlac1.F [deleted file]
PDF/spdf/phlac2.F [deleted file]
PDF/spdf/phlac3.F [deleted file]
PDF/spdf/rxdv.F [deleted file]
PDF/spdf/rxglue.F [deleted file]
PDF/spdf/rxqbar.F [deleted file]
PDF/spdf/rxud.F [deleted file]
PDF/spdf/sasan1.F [deleted file]
PDF/spdf/sasano.F [deleted file]
PDF/spdf/sasbeh.F [deleted file]
PDF/spdf/sasdir.F [deleted file]
PDF/spdf/sasgam1.F [deleted file]
PDF/spdf/sasgam2.F [deleted file]
PDF/spdf/sasvm1.F [deleted file]
PDF/spdf/sasvmd.F [deleted file]
PDF/spdf/seadis.F [deleted file]
PDF/spdf/sfabfow.F [deleted file]
PDF/spdf/sfafg1.F [deleted file]
PDF/spdf/sfbma.F [deleted file]
PDF/spdf/sfbmb.F [deleted file]
PDF/spdf/sfctq11.F [deleted file]
PDF/spdf/sfctq12.F [deleted file]
PDF/spdf/sfctq13.F [deleted file]
PDF/spdf/sfctq14.F [deleted file]
PDF/spdf/sfctq15.F [deleted file]
PDF/spdf/sfctq21.F [deleted file]
PDF/spdf/sfctq22.F [deleted file]
PDF/spdf/sfctq23.F [deleted file]
PDF/spdf/sfctq24.F [deleted file]
PDF/spdf/sfctq25.F [deleted file]
PDF/spdf/sfctq26.F [deleted file]
PDF/spdf/sfctq31.F [deleted file]
PDF/spdf/sfctq32.F [deleted file]
PDF/spdf/sfctq33.F [deleted file]
PDF/spdf/sfctq41.F [deleted file]
PDF/spdf/sfctq42.F [deleted file]
PDF/spdf/sfctq43.F [deleted file]
PDF/spdf/sfctq44.F [deleted file]
PDF/spdf/sfctq45.F [deleted file]
PDF/spdf/sfctq46.F [deleted file]
PDF/spdf/sfctq47.F [deleted file]
PDF/spdf/sfctq48.F [deleted file]
PDF/spdf/sfctq49.F [deleted file]
PDF/spdf/sfdisa.F [deleted file]
PDF/spdf/sfdisd0.F [deleted file]
PDF/spdf/sfdisdm.F [deleted file]
PDF/spdf/sfdish.F [deleted file]
PDF/spdf/sfdiss0.F [deleted file]
PDF/spdf/sfdsqa.F [deleted file]
PDF/spdf/sfehlq1.F [deleted file]
PDF/spdf/sfehlq2.F [deleted file]
PDF/spdf/sfgrvh.F [deleted file]
PDF/spdf/sfgrvl.F [deleted file]
PDF/spdf/sfgsho.F [deleted file]
PDF/spdf/sfgslo1.F [deleted file]
PDF/spdf/sfgslo2.F [deleted file]
PDF/spdf/sfkbmr2.F [deleted file]
PDF/spdf/sfkbmr5.F [deleted file]
PDF/spdf/sflacg1.F [deleted file]
PDF/spdf/sflacg2.F [deleted file]
PDF/spdf/sflacg3.F [deleted file]
PDF/spdf/sfmfap.F [deleted file]
PDF/spdf/sfmfr1.F [deleted file]
PDF/spdf/sfmfr2.F [deleted file]
PDF/spdf/sfmrsd0.F [deleted file]
PDF/spdf/sfmrsdm.F [deleted file]
PDF/spdf/sfmrsg.F [deleted file]
PDF/spdf/sfmrsj.F [deleted file]
PDF/spdf/sfmrss0.F [deleted file]
PDF/spdf/sfms105.F [deleted file]
PDF/spdf/sfms110.F [deleted file]
PDF/spdf/sfms115.F [deleted file]
PDF/spdf/sfms120.F [deleted file]
PDF/spdf/sfms125.F [deleted file]
PDF/spdf/sfms130.F [deleted file]
PDF/spdf/sfmsap.F [deleted file]
PDF/spdf/sfmsba.F [deleted file]
PDF/spdf/sfmsbb0.F [deleted file]
PDF/spdf/sfmsbd0.F [deleted file]
PDF/spdf/sfmsbdm.F [deleted file]
PDF/spdf/sfmsbh.F [deleted file]
PDF/spdf/sfmsbs0.F [deleted file]
PDF/spdf/sfmsfa.F [deleted file]
PDF/spdf/sfmsfg.F [deleted file]
PDF/spdf/sfmsjp.F [deleted file]
PDF/spdf/sfmsqa.F [deleted file]
PDF/spdf/sfmsr1.F [deleted file]
PDF/spdf/sfmsr2.F [deleted file]
PDF/spdf/sfmsr3.F [deleted file]
PDF/spdf/sfmsr4.F [deleted file]
PDF/spdf/sfmt10.F [deleted file]
PDF/spdf/sfmt11.F [deleted file]
PDF/spdf/sfmtu1.F [deleted file]
PDF/spdf/sfmtu2.F [deleted file]
PDF/spdf/sfmtu3.F [deleted file]
PDF/spdf/sfmtu4.F [deleted file]
PDF/spdf/sfmtu5.F [deleted file]
PDF/spdf/sfmtu6.F [deleted file]
PDF/spdf/sfmtu7.F [deleted file]
PDF/spdf/sfmtu8.F [deleted file]
PDF/spdf/sfmtu9.F [deleted file]
PDF/spdf/sfsas11.F [deleted file]
PDF/spdf/sfsas12.F [deleted file]
PDF/spdf/sfsas13.F [deleted file]
PDF/spdf/sfsas14.F [deleted file]
PDF/spdf/sfsas21.F [deleted file]
PDF/spdf/sfsas22.F [deleted file]
PDF/spdf/sfsas23.F [deleted file]
PDF/spdf/sfsas24.F [deleted file]
PDF/spdf/sfwhi1.F [deleted file]
PDF/spdf/sfwhi2.F [deleted file]
PDF/spdf/sfwhi3.F [deleted file]
PDF/spdf/sfwhi4.F [deleted file]
PDF/spdf/sfwhi5.F [deleted file]
PDF/spdf/sfwhi6.F [deleted file]
PDF/spdf/smrsp31.F [deleted file]
PDF/spdf/smrsp32.F [deleted file]
PDF/spdf/smrsp33.F [deleted file]
PDF/spdf/smrspi.F [deleted file]
PDF/spdf/srctq21.F [deleted file]
PDF/spdf/srctq22.F [deleted file]
PDF/spdf/srctq23.F [deleted file]
PDF/spdf/srctq24.F [deleted file]
PDF/spdf/srctq25.F [deleted file]
PDF/spdf/srctq26.F [deleted file]
PDF/spdf/strbebc.F [deleted file]
PDF/spdf/strbep.F [deleted file]
PDF/spdf/strcbp.F [deleted file]
PDF/spdf/strcep.F [deleted file]
PDF/spdf/strdo1.F [deleted file]
PDF/spdf/strdo2.F [deleted file]
PDF/spdf/strghr.F [deleted file]
PDF/spdf/strh1b.F [deleted file]
PDF/spdf/strh1e.F [deleted file]
PDF/spdf/strh2b.F [deleted file]
PDF/spdf/strh2e.F [deleted file]
PDF/spdf/strh3b.F [deleted file]
PDF/spdf/strh3b1.F [deleted file]
PDF/spdf/strh3b3.F [deleted file]
PDF/spdf/strh3e.F [deleted file]
PDF/spdf/strh3em.F [deleted file]
PDF/spdf/strh3ep.F [deleted file]
PDF/spdf/strkb0.F [deleted file]
PDF/spdf/strkb02.F [deleted file]
PDF/spdf/strkb03.F [deleted file]
PDF/spdf/strkb04.F [deleted file]
PDF/spdf/strkb05.F [deleted file]
PDF/spdf/strkb0o.F [deleted file]
PDF/spdf/strkbm.F [deleted file]
PDF/spdf/strmrs1.F [deleted file]
PDF/spdf/strmrs2.F [deleted file]
PDF/spdf/strmrs3.F [deleted file]
PDF/spdf/strowp1.F [deleted file]
PDF/spdf/strowp2.F [deleted file]
PDF/spdf/strpriv.F [deleted file]
PDF/spdf/strucb.F [deleted file]
PDF/spdf/struce.F [deleted file]
PDF/spdf/strucor.F [deleted file]
PDF/spdf/structf.F [deleted file]
PDF/spdf/structm.F [deleted file]
PDF/spdf/structp.F [deleted file]
PDF/spdf/v3.F [deleted file]
PDF/spdf/valdis.F [deleted file]
PDF/spdf/wate32.F [deleted file]
PDF/spdf/whit1g.F [deleted file]
PDF/spdf/whit1q.F [deleted file]
PDF/spdf/whit2g.F [deleted file]
PDF/spdf/whit2q.F [deleted file]
PDF/spdf/whit3g.F [deleted file]
PDF/spdf/whit3q.F [deleted file]
PDF/spdf/whit4g.F [deleted file]
PDF/spdf/whit4q.F [deleted file]
PDF/spdf/whit5g.F [deleted file]
PDF/spdf/whit5q.F [deleted file]
PDF/spdf/whit6g.F [deleted file]
PDF/spdf/whit6q.F [deleted file]
PDF/tpdf/testgam.F [deleted file]
PDF/tpdf/testgamo.F [deleted file]
PDF/tpdf/testpdf.F [deleted file]
PDF/tpdf/testpdfo.F [deleted file]
PDF/tpdf/testpdg.F [deleted file]
PDF/tpdf/testpdgo.F [deleted file]
PDF/tpdf/tpdfupv.F [deleted file]
PDF/tpdf/tpdfupvo.F [deleted file]

diff --git a/PDF/Makefile b/PDF/Makefile
deleted file mode 100644 (file)
index 333b0a8..0000000
+++ /dev/null
@@ -1,52 +0,0 @@
-############################### PDF Makefile ##################################
-
-# Include machine specific definitions
-
-include $(ALICE_ROOT)/conf/GeneralDef
-include $(ALICE_ROOT)/conf/MachineDef.$(ALICE_TARGET)
-
-PACKAGE = PDF
-
-# C++ sources
-
-SRCS          =
-
-FSRCS  = $(wildcard spdf/*.F)
-
-FOBJS  = $(patsubst %.F,tgt_$(ALICE_TARGET)/%.o,$(FSRCS))
-
-SRCS   = $(FSRCS)
-OBJS   = $(FOBJS)
-
-# C++ compilation flags
-
-CXXFLAGS      = $(CXXOPTS) $(CLIBCXXOPTS) $(CLIBDEFS) -I.
-# C compilation flags
-
-CFLAGS      = $(COPT) $(CLIBCOPT) $(CLIBDEFS) -I.
-# FORTRAN compilation flags
-
-FFLAGS      = $(FOPT) $(CLIBFOPT) $(CLIBDEFS) -I.
-# Target
-
-SLIBRARY       = $(LIBDIR)/libpdf.$(SL)
-ALIBRARY       = $(LIBDIR)/libpdf.a
-
-default:       $(SLIBRARY)
-
-$(LIBDIR)/libpdf.$(SL):  $(OBJS)
-
-depend:                        $(SRCS)
-
-TOCLEAN                = $(OBJS)
-
-############################### General Macros ################################
-
-include $(ALICE_ROOT)/conf/GeneralMacros
-
-############################ Dependencies #####################################
-
--include tgt_$(ALICE_TARGET)/Make-depend 
diff --git a/PDF/doc/inform.dat b/PDF/doc/inform.dat
deleted file mode 100644 (file)
index 89ecfcd..0000000
+++ /dev/null
@@ -1,13 +0,0 @@
-*
-* $Id$
-*
-* $Log$
-* Revision 1.1.1.1  1996/03/08 17:48:06  mclareni
-* Pdf
-*
-*
-*  GENERAL INFORMATION
-C PDF      6.06      950315 17.05 CERN Program Library     W5051 PDF
-C          Structure Functions Package
-C          Author: H. Plothow-Besch
-C
diff --git a/PDF/doc/pdfdoc.dat b/PDF/doc/pdfdoc.dat
deleted file mode 100644 (file)
index 4caf881..0000000
+++ /dev/null
@@ -1,2408 +0,0 @@
-*
-* $Id$
-*
-* $Log$
-* Revision 1.1.1.1  1996/03/08 17:48:06  mclareni
-* Pdf
-*
-*
-1***********************************************************************
- *                                                                     *
- *   CERN                                                              *
- *                                                                     *
- *   W5051 PDFLIB                                         1995.03.15   *
- *                                                                     *
- *   ---------------------------------------------------------------   *
- *                                                                     *
- *      PDFLIB: Nucleon, Pion and Photon Parton Density Functions      *
- *                      and alpha(s) Calculations                      *
- *                                                                     *
- *                           User's Manual                             *
- *                                                                     *
- *                           Version 6.06                              *
- *                                                                     *
- *                                                                     *
- *                     H. Plothow-Besch/CERN-PPE                       *
- *                                                                     *
- *                         CH-1211 Geneva 23                           *
- *            BITNET/EARN address: PLOTHOW@CERNVM.CERN.CH              *
- *                    Tel.: +41 - 22 - 767 39 80                       *
- *                                                                     *
- *   A complete manual in PostScript format exists on a separate file  *
- *     Requests can be addressed to the CERN Program Library Office    *
- *                                                                     *
- *                   Copyright by H. Plothow-Besch                     *
- *                                                                     *
- *                            Abstract                                 *
- *                                                                     *
- *       This document describes an integrated package of Parton Den-  *
- *  sity Functions called PDFLIB which has been added to the CERN      *
- *  Program Library Pool W999 and is labelled as W5051. In this pack-  *
- *  age all the different sets of parton density functions of the      *
- *  NUCLEON, the PION and the PHOTON which are available today (about  *
- *  100 different sets) have been put together. All these sets have    *
- *  been combined in a consistent way such that they all have similar  *
- *  calling sequences and no external data files have to be read in    *
- *  anymore. A default set has been prepared, although those prefering *
- *  their own set or wanting to test a new one may do so within the    *
- *  package. The selection mode, the program parameters as well as     *
- *  the possibilities to modify the defaults and to control errors     *
- *  occured during execution are described.                            *
- *                                                                     *
- *       The package offers as well a program to calculate the strong  *
- *  coupling constant alpha(s) to first or second order (default).     *
- *  The correct Lambda(QCD) associated to the selected set of struc-   *
- *  ture functions and the number of allowed flavours with respect to  *
- *  the given Q**2 is automatically used in the caculation. For those  *
- *  structure functions where the evolution has been performed to      *
- *  leading order only, an alpha(s) value fixed to first order is re-  *
- *  turned.                                                            *
- *                                                                     *
- *       The package is written in standard FORTRAN-77 and is avai-    *
- *  lable on all CERN system areas.                                    *
- *                                                                     *
- *        Submitter:   H. Plothow-Besch/CERN-PPE                       *
- *        Language:    FORTRAN 77                                      *
- *        Library:     PDFLIB                                          *
- *                                                                     *
- ***********************************************************************
-1***********************************************************************
- *                                                                     *
- *                                                                     *
- *   _______________________________________________________________   *
- *   |                                                             |   *
- *   |  PDFLIB - User's Manual                                     |   *
- *   |                                                             |   *
- *   |  CERN Program Library entry W5051                           |   *
- *   |                                                             |   *
- *   |  Copyright CERN, Geneva 1990, 1991, 1992, 1993, 1994, 1995  |   *
- *   |                                                             |   *
- *   |  Copyright and any other appropriate legal protection of    |   *
- *   |  these computer programs and associated documentation re-   |   *
- *   |  served in all countries of the world.                      |   *
- *   |                                                             |   *
- *   |  These programs or documentation may not be reproduced      |   *
- *   |  and/or redistributed by any method without prior written   |   *
- *   |  consent of the Director-General of CERN or his delegate.   |   *
- *   |                                                             |   *
- *   |  Permission for the scientific usage of any programs des-   |   *
- *   |  cribed herein is granted apriori to those scientific in-   |   *
- *   |  stitutes associated with the CERN experimental program or  |   *
- *   |  with whom CERN has concluded a scientific collaboration    |   *
- *   |  agreement.                                                 |   *
- *   |                                                             |   *
- *   |  Commercial utilisation requires explicit a priori per-     |   *
- *   |  mission from CERN and will be subjected to payment of a    |   *
- *   |  licence fee.                                               |   *
- *   |                                                             |   *
- *   |  These statements are without prejudice to any rights that  |   *
- *   |  third parties may hold in these programs.                  |   *
- *   |                                                             |   *
- *   |                                                             |   *
- *   |  Requests for information should be addressed to:           |   *
- *   |                                                             |   *
- *   |        CERN Program Library Office                          |   *
- *   |        CERN-DD Division                                     |   *
- *   |        CH-1211 Geneva 23                                    |   *
- *   |        Switzerland                                          |   *
- *   |        Tel. +41 22 767 4951                                 |   *
- *   |        Fax. +41 22 767 7155                                 |   *
- *   |        Electronic Mail address:                             |   *
- *   |                                                             |   *
- *   |        EARN/Bitnet:   CERNLIB@CERNVM                        |   *
- *   |        DECnet:        VXCERN::CERNLIB ( node 22.190 )       |   *
- *   |        Internet:      CERNLIB@CERNVM.CERN.CH                |   *
- *   |_____________________________________________________________|   *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- ******************************** - 1 - ********************************
-1***********************************************************************
- *                                                                     *
- *      Table of Contents                                              *
- *                                                               page  *
- *                                                                     *
- *      PROGRAM SUMMARY .......................................... 3   *
- *                                                                     *
- *  0.  Major Differences compared to earlier Versions             4   *
- *      0.1  The actual Version 6.06 wrt Version 5.00 ............ 4   *
- *      0.2  Version 5.00 wrt Version 4.00 ....................... 4   *
- *      0.3  Version 4.00 wrt Version 3.00 ....................... 4   *
- *      0.4  Version 3.00 wrt Version 2.00 ....................... 6   *
- *      0.5  Version 2.00 wrt Version 1.00 ....................... 7   *
- *                                                                     *
- *  1.  Introduction                                               8   *
- *                                                                     *
- *  2.  The Structure Functions Package PDFLIB                    10   *
- *      2.1  Parameter Setting and Selection of sets ............ 10   *
- *      2.2  The Calling Sequence: PDFLIB Format ................ 12   *
- *      2.3  The Calling Sequence: PDG Format ................... 14   *
- *      2.4  Error Handling ..................................... 14   *
- *      2.5  Features ........................................... 15   *
- *      2.6  PION Structure Functions ........................... 16   *
- *      2.7  PHOTON Structure Functions ......................... 17   *
- *      2.8  How to find PDFLIB ................................. 18   *
- *                                                                     *
- *  3.  alpha(s) Calculation                                      18   *
- *                                                                     *
- *  4.  Conclusions                                               19   *
- *                                                                     *
- *  5.  Acknowledgements                                          21   *
- *                                                                     *
- *  6.  Appendix                                                  22   *
- *      A : Example of how to use PDFLIB ........................ 22   *
- *      B : List of subprograms and COMMON blocks in PDFLIB            *
- *          which are relevant to the user ...................... 25   *
- *      C : Complete list of subprograms in PDFLIB .............. 26   *
- *                                                                     *
- *  7.  References                                                32   *
- *                                                                     *
- *  8.  List of Tables                                            34   *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- ******************************** - 2 - ********************************
-1***********************************************************************
- *                                                                     *
- *  PROGRAM SUMMARY                                                    *
- *                                                                     *
- *                                                                     *
- *  Titel of Program : PDFLIB, version 6.06                            *
- *                                                                     *
- *  Catalogue number : W5051 in the CERN Computer Program Library      *
- *                                                                     *
- *  Program obtainable from : CERN Computer Program Library Office     *
- *                                                                     *
- *  Computer for which the program is designed and others on which it  *
- *  has been tested:                                                   *
- *  Computer : IBM 3090, VAX, IBM-RISC, DECS, SUN, APOLLO, HP-UX,      *
- *             Silicon Graphics, NEXT, CDC, CRAY, IBM-PC               *
- *                                                                     *
- *  Operating systems : VM/CMS, VAX/VMS, ULTRIX                        *
- *                                                                     *
- *  Programming language used : FORTRAN 77, program is available also  *
- *                              in the PATCHY and the CMZ format       *
- *                                                                     *
- *  Number of lines in distributed progam,                             *
- *  including test data, etc : 99 000                                  *
- *                                                                     *
- *  Memory required,                                                   *
- *  including testprograms and some documentation : 5 700 000 bytes    *
- *                                                                     *
- *  No. of bits in a word : 64                                         *
- *                                                                     *
- *  Peripherals used : line printer (optional)                         *
- *                                                                     *
- *  Keywords : parton density functions, structure functions, W5051,   *
- *             nucleon-, proton-, pion-, photon density parametrisa-   *
- *             tions, Altarelli-Parisi equation, strong coupling       *
- *             constant alpha(s)                                       *
- *                                                                     *
- *  Nature of physical problem :                                       *
- *  All theoretical calculations of cross sections involving the       *
- *  hadronic structure of the incoming particles, i.e. protons, pions  *
- *  or photons, use the parton density at a given longitudinal momen-  *
- *  tum X of the parton and at a given momentum transfer Q**2. There   *
- *  is a variety of parametrisations available, but each of them has   *
- *  its own structure for input and/or output.                         *
- *                                                                     *
- *  Method of solution :                                               *
- *  Collection of all available parton density parametrisations and    *
- *  application in an identical format for all of them, regardless of  *
- *  the author's convention or of the particle type.                   *
- *                                                                     *
- *  References :                                                       *
- *  [1] H. Plothow-Besch, 'PDFLIB: Structure Functions and alpha(s)    *
- *      Calculations', User's Manual - Version 1.00, W5051 PDFLIB,     *
- *      1991.03.21, CERN-PPE.                                          *
- *  [2] H. Plothow-Besch, 'PDFLIB: a library of all available parton   *
- *      density functions of the nucleon, the pion and the photon and  *
- *      the corresponding alpha(s) calculations',                      *
- *      Comp. Phys. Comm. 75 (1993) 396-416.                           *
- *                                                                     *
- *                                                                     *
- ******************************** - 3 - ********************************
-1***********************************************************************
- *                                                                     *
- *  0.  Major Differences compared to earlier Versions                 *
- *                                                                     *
- *  0.1  The actual Version 6.06 wrt Version 5.00:                     *
- *                                                                     *
- *     The following major changes compared to version 5.00 of PDFLIB  *
- *  have been implemented:                                             *
- *                                                                     *
- *  Several NUCLEON structure function sets have been added:           *
- *    - the NLL sets A' and G of Martin, Roberts and Stirling (MRS-Ap) *
- *      and (MRS-G) of February 1995 with grid data as well as in its  *
- *      parametrized version (MRS-Ap-F) and (MRS-G-Fit),               *
- *    - the version 3 of the CTEQ group of October 1994 as NLL parame- *
- *      trisation in the MS_bar as well as in the DIS renormalisation  *
- *      scheme and as LO (CTEQ3M, CTEQ3D, CTEQ3L),                     *
- *    - the new sets of GRV of November 1994 as NLL in the MS_bar as   *
- *      well as in the DIS renormalisation scheme and as LO GRV        *
- *      (GRV94-MS, GRV94DI, GRV94-LO),                                 *
- *                                                                     *
- *  a NEW group of PHOTON structure function sets have been added:     *
- *    - the six NLL sets of the WHIT group (WHIT1-G to WHIT6-G).       *
- *                                                                     *
- *  Please note that as new the default set the MRS set (G) has been   *
- *  chosen (NPTYPE = 1, NGROUP = 3, NSET = 41).                        *
- *                                                                     *
- *     More details can be found on the next pages of this updated     *
- *  User's Manual.                                                     *
- *                                                                     *
- *     The response from the Physics community has been very satis-    *
- *  factory, which is encouraging us to keep up-to-date with the       *
- *  latest developments in this area. The author will be pleased to    *
- *  learn about new parton density functions, as well as to receive    *
- *  suggestions to improve both, the usage and the documentation.      *
- *                                                                     *
- *                                                                     *
- *  0.2  Version 5.00 wrt Version 4.00:                                *
- *                                                                     *
- *     The following major changes compared to version 4.00 of PDFLIB  *
- *  have been implemented:                                             *
- *                                                                     *
- *  NUCLEON structure function sets have been added:                   *
- *    - the NLL set A of Martin, Roberts and Stirling (MRS-A) and the  *
- *      parametrized version (MRS-A-Fit)                               *
- *  PHOTON structure function sets have been added:                    *
- *    - the NLL set of Aurenche, Fontannaz and Guillet (AFG-G).        *
- *                                                                     *
- *                                                                     *
- *  0.3  Version 4.00 wrt Version 3.00:                                *
- *                                                                     *
- *     The following major changes compared to version 3.00 of PDFLIB  *
- *  have been implemented:                                             *
- *                                                                     *
- *  NUCLEON structure function sets have been updated or added:        *
- *    - the updated NLL sets of Martin, Roberts and Stirling (MRS) of  *
- *      November 1992, the sets S0', D0' and D-' in the MSbar as well  *
- *      as the sets S0', D0' and D-' in the DIS renormalisation        *
- *      scheme. Please note that the previous S0, D0 and D- sets have  *
- *                                                                     *
- ******************************** - 4 - ********************************
-1***********************************************************************
- *                                                                     *
- *      been retracted by the authors. The new MRS(H) set in both      *
- *      schemes have also been added,                                  *
- *    - the NLL and LO sets of the CTEQ collaboration (CTEQ) in its    *
- *      parametrised form of January 1993. These are the sets 1M, 1MS  *
- *      and 1ML in the MSbar, the set 1D in the DIS renormalisation    *
- *      scheme and the LO set 1L,                                      *
- *    - the new NLL sets A and B of Berger and Meng (BM) of            *
- *      February 1993 with special gluon distributions, and finally    *
- *    - the new NLL and LO sets of the CTEQ collaboration (CTEQ) in    *
- *      its parametrised form of August 1993. These are the sets 2M,   *
- *      2MS, 2MF and 2ML in the MSbar, the set 2D in the DIS renorma-  *
- *      lisation scheme and the LO set 2L. The revised CTEQ2 version,  *
- *      CTEQ2p, of November 1993 is recommended by the authors.        *
- *                                                                     *
- *  Please note that the default set has been put to the GRV set HO    *
- *  (NPTYPE = 1, NGROUP = 5, NSET = 3 / former MODE = 72).             *
- *  The MOST DRASTICAL change of the package has been, however, the    *
- *  NEW FORMAT of PARAMETER SETTINGS. Instead of using only ONE para-  *
- *  meter (MODE) to select a set of parton density functions (PDF),    *
- *  each PDF set is now identified by THREE parameters. These are      *
- *  NPTYPE, NGROUP and NSET, for PARTICLE TYPE (Nucleons, Pions or     *
- *  Photons) AUTHOR GROUP (f.ex. MRS, CTEQ, GRV, etc) and the PARAME-  *
- *  TRISATION SET within the group. Therefore the user should provide  *
- *  THREE parameters (NPTYPE, NGROUP, NSET) to the subroutine PDFSET,  *
- *  each time a PDF set shall be selected. But to make life nice and   *
- *  easy, it is also possible to select a PDF set by transferring only *
- *  ONE parameter to PDFSET as before: this parameter should then      *
- *  identify the NAME of the author group in an unequivocal way; the   *
- *  selection of the PDF set within the author group is made by the    *
- *  corresponding VALUE value (f.ex. PARM(1) = 'MRS' with VALUE(1) =   *
- *  29.0D0 would select the 29th set of the Nucleon PDFs of MRS, and   *
- *  f.ex. PARM(1) = 'GRVph' with VALUE(1) = 1.0D0 would select the     *
- *  first set of the Photon PDFs of GRV). As a suggestion, the NAME for*
- *  each PDF set can be found in the NEW COMMON block /W505110/SFNAME, *
- *  where SFNAME is a CHARACTER*8 array with dimension SFNAME(NPTYMX,  *
- *  NGRMAX,NSETMX), and NPTYMX, NGRMAX, NSETMX are defined by the      *
- *  parameter sequence, +SEQ, W5051P2. In case the user wants to se-   *
- *  lect the DEFAULT, it is sufficient to transfer as the ONLY parame- *
- *  ter PARM(1) = 'Nset' with VALUE(1) = 0.0D0. This change of format  *
- *  had been necessary to take into account the inflationary growing   *
- *  number of PDFs until now and to foresee already the possibility of *
- *  adding new PDFs in a more flexible way. Please note that in case   *
- *  of any INCONSISTANCY of parameter settings, it is always chosen    *
- *  the DEFAULT settings !
- *  Please also note that the internal COMMON block /W50511/ has been  *
- *  modified to /W50511/ NPTYPE,NGROUP,NSET,MODE,NFL,LO,TMAS.          *
- *                                                                     *
- *  Full backward compatibility with previous versions is ensured !!   *
- *  This means that old programs using PDF sets up to version 3.00 of  *
- *  PDFLIB should run without ANY change. To ensure this, four new     *
- *  COMMON blocks have been added:                                     *
- *    - /W50519/ NEWVER :  to flag the NEW/OLD version with the        *
- *      LOGICAL variable NEWVER,                                       *
- *    - /W505120/ NPGSMX(NPTYMX,NGRMAX),NSETFL(NPTYMX,NGRMAX,NSETMX) : *
- *      to provide the actual maximum number of parametrisation sets   *
- *      per particle type and author group, and to flag the retracted  *
- *                                                                     *
- ******************************** - 5 - ********************************
-1***********************************************************************
- *                                                                     *
- *      PDF sets (NSETFL(I,J,K) = 0),                                  *
- *    - /W505121/ NPTYCR(MODEMX),NGROMX(MODEMX),NSETCR(MODEMAX) :      *
- *      to provide the cross reference to the NPTYPE, NGROUP and NSET  *
- *      value (NEW format) for a given MODE value (OLD format),        *
- *    and                                                              *
- *    - /W505122/ MODECR(NPTYMX,NGRMAX,NSETMX) :  to provide the       *
- *      cross reference to the MODE value (OLD version) for a given    *
- *      NPTYPE, NGROUP and NSET value (NEW format). In case of the     *
- *      above added PDF sets NO MODE value is available; in that case  *
- *      MODECR(I,J,K) = -1 is returned.                                *
- *                                                                     *
- *                                                                     *
- *  0.4  Version 3.00 wrt Version 2.00                                 *
- *                                                                     *
- *     The following major changes compared to version 2.00 of PDFLIB  *
- *  had been implemented:                                              *
- *                                                                     *
- *  NUCLEON structure function sets had been updated or added:         *
- *    - the NLL set of Aurenche et al. (ABFOW),                        *
- *    - the updated NLL set of Kwiecinski, Martin, Roberts and         *
- *      Stirling (KMRS) set B0 with u_bar not equal to d_bar           *
- *      of April 1992,                                                 *
- *    - the new NLL sets of Martin, Roberts and Stirling (MRS)         *
- *      set S0, D0 and D- with u_bar not equal to d_bar of April 1992. *
- *  Please note that the default set had been put to the MRS set S0    *
- *  (MODE = 45).                                                       *
- *                                                                     *
- *  PION structure function sets had been added:                       *
- *    - the new LO and NLL sets of Gl"uck, Reya and Vogt (GRV-P)       *
- *      of August 1991,                                                *
- *    - the new sets 2 and 3 of Aurenche et al. (ABFKW-P).             *
- *  Please also note that the ABFKW set had changed MODE number from   *
- *  110 -> 180, as well as the SMRS sets from 120ff -> 130ff, and      *
- *  that earlier inconsistencies in the parton density definitions     *
- *  had been corrected.                                                *
- *                                                                     *
- *  For the first time also PHOTON structure function sets had been    *
- *  added to the package. These sets were the following:               *
- *    - the old Duke and Owens (DO-G) LO and NLL sets,                 *
- *    - the old LO sets 1 to 4 of Drees and Grassie (DG-G),            *
- *    - the new LO sets 1 to 3 of Abramowicz, Charchula and Levy       *
- *      (LAC-G),                                                       *
- *    - the new LO sets 1 and 2 and the NLL set of Gordon and Storrow  *
- *      (GS-G) of July 1991,                                           *
- *    - the new LO and NLL sets, as well as the leading term of the    *
- *      NLL calculation of Gl"uck, Reya and Vogt (GRV-G)               *
- *      of December 1991,                                              *
- *    - the new NLL sets of Aurenche et al. (ACFGP-G) of April 1992.   *
- *                                                                     *
- *  Some improvements in the usage of the package had been made:       *
- *    - more protections had been implemented and inconsistencies      *
- *      had been corrected to make the package more robust and         *
- *      available on ALL system areas, and finally                     *
- *    - the subroutine STRUCTM had been added. STRUCTM is basically    *
- *      the same routine as the old steering subroutine STRUCTF apart  *
- *      from ONE additional parameter in the argument list which had   *
- *                                                                     *
- ******************************** - 6 - ********************************
-1***********************************************************************
- *                                                                     *
- *      been added to transfer separately the parton densities for     *
- *      u_bar and d_bar. The argument list of the subroutine STRUCTF   *
- *      was as before; in case of STRUCTF sea = (u_bar + d_bar)/2 is   *
- *      transferred. The subroutine PFTOPDG, which is the sterring     *
- *      routine in PDG format, was unchanged; this routine always      *
- *      transferred u_bar and d_bar separately.                        *
- *                                                                     *
- *                                                                     *
- *  0.5  Version 2.00 wrt Version 1.00                                 *
- *                                                                     *
- *     The following major changes compared to version 1.00 of PDFLIB  *
- *  had been implemented:                                              *
- *                                                                     *
- *  New parton density functions of the NUCLEON had been added:        *
- *    - the old LO set of Buras and Gaemers of 1978 (BEBC) for com-    *
- *      pleteness,                                                     *
- *    - all the LO sets of Diemoz, Ferroni, Longo and Martinelli (DFLM)*
- *      including their parametrisation for the NLL average fit,       *
- *    - the new LO set 1.1 of Owens (DO 1.1) of June 1991,             *
- *    and finally                                                      *
- *    - the new LO and NLL sets of Gl"uck, Reya and Vogt (GRV) of      *
- *      May 1991.                                                      *
- *                                                                     *
- *  PION structure function sets had been added to the package:        *
- *    - the old Owens (OW-P) sets 1 and 2,                             *
- *    - the newer Aurenche et al. (ABFKW-P) set and                    *
- *    - the new Sutton, Martin, Roberts and Stirling (SMRS-P) sets 1,  *
- *      2 and 3.                                                       *
- *                                                                     *
- *  Some improvements in the usage of the package had been made:       *
- *    - COMMON block names had been changed to avoid interference      *
- *      to those defined by the user (/IMODE/ -> /W50511/,             *
- *      /ALAMB/ -> /W50512/),                                          *
- *    - three new COMMON blocks had been added:                        *
- *      a) /W50513/ XMIN,XMAX,Q2MIN,Q2MAX :                            *
- *         to provide the minimum and maximum values of X and Q**2 for *
- *         each set of structure functions,                            *
- *      b) /W50510/ IFLPRT :                                           *
- *         to provide a print flag at different levels to print        *
- *         b.1) the variables of the three COMMON blocks               *
- *              /W50511/, /W50512/ and /W50513/ at job initialisation  *
- *         b.2) error messages during job execution                    *
- *      c) /W50514W/ PDFWGT :                                          *
- *         to take into account the event weight for error counting    *
- *    - a subroutine PDFSTA had been added to allow at job termination *
- *      for a summary how often X and Q**2 limits have been exceeded,  *
- *    - the MODE number of some structure functions had been changed   *
- *      to keep consistensy,                                           *
- *    - the number of parameters which can be set through the sub-     *
- *      routine PDFSET had been extended to 20,                        *
- *    - more protections had been implemented and inconsistencies      *
- *      had been corrected to make the package more robust,            *
- *    and finally                                                      *
- *    - a subroutine PFTOPDG had been provided as an interface for     *
- *      those users who prefer the flavour code convention of the PDG  *
- *      instead of the PDFLIB format.                                  *
- *                                                                     *
- ******************************** - 7 - ********************************
-1***********************************************************************
- *                                                                     *
- *  1.  Introduction                                                   *
- *                                                                     *
- *                                                                     *
- *     The number of useful sets of parton density functions of the    *
- *  NUCLEON available today is about 70 and is still growing. The ol-  *
- *  der sets are leading order evolutions only, like                   *
- *                                                                     *
- *   -  Buras and Gaemers (BEBC) [1],                                  *
- *   -  Owens and Reya (OR) [2],                                       *
- *   -  Baier, Engels and Petersson (BEP) [3],                         *
- *   -  Gl"uck, Hoffmann and Reya (GHR) [4],                           *
- *   -  Duke and Owens (DO) sets 1 and 2 [5],                          *
- *   -  Eichten, Hinchliffe, Lane and Quigg (EHLQ) sets 1 and 2 [6],   *
- *   -  Diemoz, Ferroni, Longo and Martinelli (DFLM) with              *
- *      Lambda_4(QCD) = 200 MeV for soft and hard valence quark and    *
- *      soft and hard gluon distributions [7], and finally             *
- *   -  the new LO paramerisation of Owens set 1.1 (DO 1.1) [8].       *
- *                                                                     *
- *     Most of the sets, essentially the more recent ones, are         *
- *  next-to-leading order evolutions, like                             *
- *                                                                     *
- *   -  Diemoz, Ferroni, Longo and Martinelli (DFLM) with              *
- *      Lambda_4(QCD) = 160, 260, 360 MeV [7],                         *
- *   -  Martin, Roberts and Stirling (MRS) sets 1 to 3, sets E and B   *
- *      and sets E' and B' [9],                                        *
- *   -  Harriman, Martin, Roberts and Stirling (HMRS) sets E, E+, E-   *
- *      and B [10],                                                    *
- *   -  Kwiecinski, Martin, Roberts and Stirling (KMRS) [11] sets B0   *
- *      and B-, and for the B- set with different radii for shadowing, *
- *   -  Martin, Roberts and Stirling (MRS-B Lambda) [11] for the B0    *
- *      set also with different values of Lambda_4(QCD) = 135, 160,    *
- *      200 and 235 MeV,                                               *
- *   -  Martin, Roberts and Stirling (MRS) [12] sets S0, D0 and D-     *
- *      as well as the upgraded KMRS set B0, all having u_bar not      *
- *      equal to d_bar,                                                *
- *   -  Morfin and Tung (MT) sets 1 to 6 [13],                         *
- *   -  Gl"uck, Reya and Vogt (GRV) [14] including their new para-     *
- *      metrisation from August 1991,                                  *
- *   -  Aurenche et al. (ABFOW) [15],                                  *
- *   -  the upgraded Martin, Roberts and Stirling (MRS) [16] sets      *
- *      S0', D0' and D-' in the MSbar renormalisation scheme and MRS(H)*
- *   -  the upgraded Martin, Roberts and Stirling (MRS) [17] sets      *
- *      S0', D0' and D-' in the DIS renormalisation scheme and MRS(H), *
- *    - the sets 1M, 1MS and 1ML in the MSbar and the set 1D in the    *
- *      DIS renormalisation scheme of the CTEQ collaboration (CTEQ)    *
- *      [18] in its parametrised form,                                 *
- *    - the new sets A and B of Berger and Meng (BM) [19] with special *
- *      gluon distributions, and                                       *
- *    - the new sets 2M, 2MS, 2MF and 2ML in the MSbar and the set 2D  *
- *      in the DIS renormalisation scheme of the CTEQ collaboration    *
- *      (CTEQ) [20] in its parametrised form which superseed already   *
- *      the previous CTEQ sets 1M, 1MS, 1ML and 1D. The CTEQ collabo-  *
- *      ration provides also parametrisations to leading order, the    *
- *      sets 1L [18] and 2L [20].                                      *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- ******************************** - 8 - ********************************
-1***********************************************************************
- *                                                                     *
- *     There is also a limited number of parton density functions of   *
- *  the PION available. These are the old set to leading order evolu-  *
- *  tion of                                                            *
- *                                                                     *
- *   -  Owens (OW-P) [21], sets 1 and 2,                               *
- *                                                                     *
- *  and the newer and the very new next-to-leading order evolutions of *
- *                                                                     *
- *   -  Aurenche et al. (ABFKW-P) [22], sets 1, 2 and 3,               *
- *   -  Sutton, Martin, Roberts and Stirling (SMRS-P) [23], sets 1, 2  *
- *      and 3 and                                                      *
- *   -  Gl"uck, Reya and Vogt (GRV-P) [24], the LO and NLL sets.       *
- *                                                                     *
- *     In addition there are also a few parton density functions of    *
- *  the PHOTON available. These are the old sets to leading order evo- *
- *  lution of                                                          *
- *                                                                     *
- *   -  Duke and Owens (DO-G) (and the NLL) [25] using an asymptotic   *
- *      solution of the Altarelli-Parisi equation,                     *
- *   -  Drees and Grassie (DG-G) sets 1 to 4 [26] using a full solu-   *
- *      tion of the Altarelli-Parisi equation,                         *
- *                                                                     *
- *  the new leading order evolution of                                 *
- *                                                                     *
- *   -  Abramowicz, Charchula and Levy (LAC-G) [27] sets 1 to 3 using  *
- *      a full solution of the Altarelli-Parisi equation,              *
- *                                                                     *
- *  and the new next-to-leading order evolution of                     *
- *                                                                     *
- *   -  Gordon and Storrow (GS-G) [28], the LO sets 1 and 2 and the    *
- *      NLL set,                                                       *
- *   -  Gl"uck, Reya and Vogt (GRV-G) [29], the LO, NLL and the        *
- *      leading term of the NLL sets, and                              *
- *   -  Aurenche et al. (ACFGP-G) [30] sets 1 and 2, without and with  *
- *      massive charm.                                                 *
- *                                                                     *
- *     The different sets differ also in the renormalisation scheme    *
- *  they use. The DFLM series, the sets 1 to 6 of the MT series, the   *
- *  new MRS sets S0', D0' and D-' and the new CTEQ sets 1D and 2D are  *
- *  performed in the DIS renormalisation scheme [31]. If the LO evo-   *
- *  lutions of these authors are ignored, all the other sets are       *
- *  performed in the MS_bar renormalisation scheme. All sets are made  *
- *  using four flavours in the initial state (NF = 4). The value of    *
- *  the QCD scale factor, Lambda_4(QCD), ranges from 45 to 500 MeV.    *
- *                                                                     *
- *     We have put together all these different sets of parton density *
- *  functions in one single package [32]. We have modified the sets    *
- *  such that no external data files for the grids have to be read in  *
- *  anymore. All these structure function sets have been combined in a *
- *  consistent way such that they all appear in an IDENTICAL STRUCTURE *
- *  to the user. The selection is made via THREE parameters, NPTYPE,   *
- *  NGROUP and NSET, which identify a PDF set by its PARTICLE TYPE     *
- *  (Nucleon: NPTYPE = 1, Pion: NPTYPE = 2, Photon: NPTYPE = 3), its   *
- *  AUTHOR GROUP (f.ex. MRS: NGROUP = 3, CTEQ: NGROUP = 4,             *
- *  GRV: NGROUP = 5) and its PDF SET WITHIN THE GROUP (NSET = 1,2,3    *
- *  etc.). It is also possible to select a PDF set by transferring     *
- *                                                                     *
- ******************************** - 9 - ********************************
-1***********************************************************************
- *                                                                     *
- *  only ONE parameter to PDFSET as before: this parameter should then *
- *  identify the NAME of the author group in an unequivocal way; the   *
- *  selection of the PDF set within the author group is made by the    *
- *  corresponding VALUE value (f.ex. PARM(1) = 'MRS' with VALUE(1) =   *
- *  29.0D0 would select the 29th set of the Nucleon PDFs of MRS, and   *
- *  f.ex. PARM(1) = 'GRVph' with VALUE(1) = 1.0D0 would select the     *
- *  first set of the Photon PDFs of GRV). As a suggestion, the NAME for*
- *  each PDF set can be found in the new COMMON block /W505110/SFNAME, *
- *  where SFNAME is a CHARACTER*8 array with dimension SFNAME(NPTYMX,  *
- *  NGRMAX,NSETMX), and NPTYMX, NGRMAX, NSETMX are defined by the      *
- *  parameter sequence, +SEQ, W5051P2. The parameters should be,       *
- *  NPTYPE, NGROUP, NSET or the 'NAME' of the PDF set with their VALUE *
- *  values should be set with a call to the subroutine PDFSET at the   *
- *  initialisation phase. A simple SUBROUTINE call, which is IDENTICAL *
- *  for all applications (nucleons, pions, photons) returns the parton *
- *  densities for all partons (u, d, s, c, b, t, gluon and their anti- *
- *  quarks) at a given X value (where X is the fraction of the longi-  *
- *  tudinal momentum carried by the parton) and the Q-scale SCALE (in  *
- *  GeV). A default set has been prepared, in which case it is suffi-  *
- *  cient to transfer as the ONLY parameter PARM(1) = 'Nset' with      *
- *  VALUE(1) = 0.0D0, although those preferring their own private set  *
- *  or wanting to test a new one may do so within the package. Error   *
- *  control can be obtained by setting a print flag to obtain output   *
- *  either during execution and/or as a summary at job termination via *
- *  a call to the subroutine PDFSTA. In addition, the different parton *
- *  density sets have been modified such that the library can be used  *
- *  on all the different computer systems known today. The library has *
- *  been tested on the different system areas to ensure identical re-  *
- *  sults within the given machine precision.                          *
- *                                                                     *
- *     Within the same package a program is also provided to calculate *
- *  the strong coupling constant alpha(s) to second order (default)    *
- *  or to first order (by user's choice). The correct Lambda(QCD) as-  *
- *  sociated to the selected set of structure functions and the number *
- *  of allowed flavours with respect to the given Q**2 is automati-    *
- *  cally used in the calculation. For those structure functions where *
- *  the evolution has been performed to leading order only, the        *
- *  alpha(s) value to first order is returned.                         *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *  2.  The Structure Functions Package PDFLIB                         *
- *                                                                     *
- *                                                                     *
- *     In the following the use of and the access to the PDFLIB pack-  *
- *  age is described.                                                  *
- *                                                                     *
- *  2.1  Parameter Setting and Selection of Sets                       *
- *                                                                     *
- *                                                                     *
- *     The user may access his/her preferred set of structure func-    *
- *  tions by setting the THREE parameters, NPTYPE, NGROUP and NSET,    *
- *  which identify each set of parton density functions via a call     *
- *  to the subroutine                                                  *
- *                                                                     *
- *       CALL PDFSET(PARM,VALUE)                                       *
- *                                                                     *
- ******************************* - 10 - ********************************
-1***********************************************************************
- *                                                                     *
- *  at the initialization phase of his/her MAIN program, or use the    *
- *  package as a 'black box' with the default values (see below). The  *
- *  arguments PARM and VALUE are vectors of dimension 20 and have the  *
- *  following meaning                                                  *
- *                                                                     *
- *       PARM(I)  = character*20 variable, which defines in any order  *
- *                  the variables                                      *
- *                  'NPTYPE', 'NGROUP', 'NSET', 'MODE', 'INIT0',       *
- *                  'NFL', 'LO', 'TMAS', 'QCDL4', 'QCDL5' and          *
- *                  'XMIN', 'XMAX', 'Q2MIN', 'Q2MAX'                   *
- *       VALUE(I) = the corresponding numerical value of the variable  *
- *                  PARM(I)                                            *
- *                  (TMAS, QCDL4, QCDL5, XMIN, XMAX, Q2MIN and Q2MAX   *
- *                  are DOUBLE PRECISION variables),                   *
- *                                                                     *
- *  where                                                              *
- *                                                                     *
- *       NPTYPE = number of particle type ranging from 1 to 3          *
- *                (Nucleons: NPTYPE = 1, Pions: NPTYPE = 2 and         *
- *                 Photons: NPTYPE = 3) (Default: NPTYPE = 1)          *
- *       NGROUP = number of author group ranging from 1 to 7           *
- *                (Default: NGROUP = 5)                                *
- *       NSET   = number of a selected structure function set within   *
- *                the author group ranging from 1 to 34                *
- *                (Default: NSET = 3)                                  *
- *      (MODE   = number of a selected structure function set ranging  *
- *                from 0 to 281 - OLD format !!)                       *
- *                (Default: MODE = 45)                                 *
- *       INIT0  = in case of PARM(1) = 'INIT0' PDFSET fills as the     *
- *                only action the COMMON blocks /W505110/, /W505120/   *
- *                /W505121/ and /W505122/                              *
- *       NFL    = desired number of flavours in the alpha(s) calcula-  *
- *                tion ranging from 3 to 6                             *
- *                (Default: NFL = 5)                                   *
- *       LO     = order of alpha(s) calculation; if LO = 1, alpha(s)   *
- *                to first order only                                  *
- *                (Default: LO = 2)                                    *
- *       TMAS   = the user defined value of the top-quark mass in      *
- *                GeV/c**2 (optional)                                  *
- *                (Default: TMAS = 100.0D0)                            *
- *       QCDL4  = QCD scale, Lambda_4(QCD), in GeV for four flavours   *
- *       QCDL5  = QCD scale, Lambda_5(QCD), in GeV for five flavours   *
- *                corresponding to QCDL4                               *
- *                                                                     *
- *  and                                                                *
- *                                                                     *
- *       XMIN   = minimum allowed X value                              *
- *       XMAX   = maximum allowed X value                              *
- *       Q2MIN  = minimum allowed Q**2 value (in (GeV/c)**2)           *
- *       Q2MAX  = maximum allowed Q**2 value (in (GeV/c)**2)           *
- *                                                                     *
- *  for each set of structure functions.                               *
- *                                                                     *
- *     If the user wants to select the default settings, it is suffi-  *
- *  cient to give PARM(1) = 'Nset' with VALUE(1) = 0.0D0 as the ONLY   *
- *  parameter to PDFSET. Please note that a PDF set may also be iden-  *
- *                                                                     *
- ******************************* - 11 - ********************************
-1***********************************************************************
- *                                                                     *
- *  tified by the NAME of the author group provided in an unequivocal  *
- *  way. The selection of the PDF set within the author group is, in   *
- *  that case, given by the corresponding VALUE value. The NAME for    *
- *  each PDF set can be found in the COMMON block /W505110/SFNAME,     *
- *  where SFNAME is a CHARACTER*8 array with dimension SFNAME(NPTYMX,  *
- *  NGRMAX,NSETMX), and NPTYMX, NGRMAX, NSETMX are defined by the para-*
- *  meter sequence, +SEQ, W5051P2. Please also note that the variable  *
- *  names can be transferred in lower, upper or mixed character modes. *
- *                                                                     *
- *     Please also note that PDFSET can be called as often as the user *
- *  likes. In order to redefine the parameters to select other sets of *
- *  structure functions which the user wants to investigate, it might  *
- *  be necessary to call PDFSET in an alternating way, but it is al-   *
- *  ways mandatory to transfer either the THREE parameters, NPTYPE,    *
- *  NGROUP, NSET, or the 'NAME' of the author group with their corres- *
- *  ponding VALUE values. The subroutine PDFSET fills the internal     *
- *  COMMON blocks                                                      *
- *                                                                     *
- *       COMMON/W50511/ NPTYPE,NGROUP,NSET,MODE,NFL,LO,TMAS            *
- *       COMMON/W50512/ QCDL4,QCDL5                                    *
- *       COMMON/W50513/ XMIN,XMAX,Q2MIN,Q2MAX                          *
- *                                                                     *
- *  at the time the routine is called. Please note that in case of a   *
- *  multiple call to PDFSET with fewer parameters redefined than in a  *
- *  preceeding call, always the last parameters are kept in memory.    *
- *  All variables of the three COMMON blocks can be automatically      *
- *  printed at job initialisation by setting the print flag IFLPRT in  *
- *  the COMMON block                                                   *
- *                                                                     *
- *       COMMON/W50510/ IFLPRT                                         *
- *                                                                     *
- *  to IFLPRT = 2.                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *  2.2  The Calling Sequence: PDFLIB Format                           *
- *                                                                     *
- *                                                                     *
- *     The main steering routine for a set of structure functions is   *
- *  accessed as follows:                                               *
- *                                                                     *
- *       CALL STRUCTM(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL)    *
- *                                                                     *
- *     Please note that all variables are defined as DOUBLE PRECISION. *
- *  The user has to provide the following INPUTs:                      *
- *                                                                     *
- *       X      = X value of parton                                    *
- *       SCALE  = QCD scale in GeV                                     *
- *                                                                     *
- *     The subroutine STRUCTM returns the following OUTPUT:            *
- *                                                                     *
- *       UPV    = up valence quark                                     *
- *       DNV    = down valence quark                                   *
- *       USEA   = sea (up_bar)                                         *
- *       DSEA   = sea (down_bar)                                       *
- *       STR    = strange quark                                        *
- *                                                                     *
- ******************************* - 12 - ********************************
-1***********************************************************************
- *                                                                     *
- *       CHM    = charm quark                                          *
- *       BOT    = bottom quark                                         *
- *       TOP    = top quark                                            *
- *       GL     = gluon                                                *
- *                                                                     *
- *     This steering routine has separate arguments for u_bar and      *
- *  d_bar-quarks to take into account that some of the new sets of     *
- *  nucleon density functions provide separately the densities for     *
- *  u_bar and d_bar. In case u_bar is not given separately from d_bar  *
- *  it is set USEA = DSEA.                                             *
- *                                                                     *
- *  Please note that STRUCTM returns  X * parton distribution          *
- *  function !                                                         *
- *                                                                     *
- *     An older version of the main steering routine
- *                                                                     *
- *       CALL STRUCTF(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,TOP,GL)          *
- *                                                                     *
- *  is still available. Also here all variables are defined as         *
- *  DOUBLE PRECISION. Please note that STRUCTF has only ONE argument   *
- *  for the sea-quark, which gives u_bar = d_bar or (u_bar + d_bar)/2. *
- *  The user has to provide as above the following INPUTs:             *
- *                                                                     *
- *       X      = X value of parton                                    *
- *       SCALE  = QCD scale in GeV                                     *
- *                                                                     *
- *     The subroutine STRUCTF returns the following OUTPUT:            *
- *                                                                     *
- *       UPV    = up valence quark                                     *
- *       DNV    = down valence quark                                   *
- *       SEA    = sea (up_bar = down_bar or (up_bar + down_bar)/2)     *
- *       STR    = strange quark                                        *
- *       CHM    = charm quark                                          *
- *       BOT    = bottom quark                                         *
- *       TOP    = top quark                                            *
- *       GL     = gluon                                                *
- *                                                                     *
- *     Please note that also STRUCTF returns  X * parton distribution  *
- *  function !                                                         *
- *                                                                     *
- *     A list of available sets of NUCLEON structure functions in the  *
- *  package is given in Table 1.                                       *
- *                                                                     *
- *     The recommended set of structure functions from the different   *
- *  authors is indicated with a flash. If NSET is set to zero          *
- *  (or if any of the three parameters, NPTYPE, NGROUP and NSET, is    *
- *  undefined) a default set is selected which is a nucleon parton     *
- *  dennsity set, the set of MRS (G) with Lambda_4(QCD) = 255 MeV      *
- *  (NPTYPE = 1, NGROUP = 3, NSET = 41).                               *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- ******************************* - 13 - ********************************
-1***********************************************************************
- *                                                                     *
- *  2.3  The Calling Sequence: PDG Format                              *
- *                                                                     *
- *                                                                     *
- *     To allow the use of the package with the flavour code conven-   *
- *  tion of the Particle Data Group (PDG) an interface has been writ-  *
- *  ten which translates the PDFLIB format (see Section 2.2) into the  *
- *  PDG format. Instead of calling the subroutine STRUCTM the user     *
- *  preferring the flavour code convention of the PDG accesses the     *
- *  package via a call to the subroutine                               *
- *                                                                     *
- *       CALL PFTOPDG(X,SCALE,DXPDF).                                  *
- *                                                                     *
- *     All variables are defined as DOUBLE PRECISION as before, and    *
- *  X and SCALE are INPUTs provided by the user having the same        *
- *  meaning as above (see Section 2.2).                                *
- *                                                                     *
- *     The subroutine PFTOPDG outputs a vector DXPDF(-6:6) of which    *
- *  its variables have the following meaning:                          *
- *                                                                     *
- *       DXPDF(0)  = gluon                                             *
- *       DXPDF(1)  = sum of down valence + down_bar quark              *
- *       DXPDF(2)  = sum of up valence + up_bar quark                  *
- *       DXPDF(3)  = strange quark                                     *
- *       DXPDF(4)  = charm quark                                       *
- *       DXPDF(5)  = bottom quark                                      *
- *       DXPDF(6)  = top quark                                         *
- *  and  DXPDF(-1) to DXPDF(-6) are the corresponding antiquarks.      *
- *                                                                     *
- *     DXPDF(-1) = DXPDF(-2) stands for down_bar = up_bar or           *
- *  (down_bar + up_bar)/2 in almost all parametrisations. In case      *
- *  where u_bar is not equal to d_bar (like for the retracted MRS S0,  *
- *  D0 and D- sets, the updated MRS S0', D0' and D-' in MS_bar and DIS *
- *  renormalisation scheme as well as the upgraded KMRS set B0 and one *
- *  of the MT series) DXPDF(-1) = d_bar and DXPDF(-2) = u_bar. In all  *
- *  sets is DXPDF(3) = DXPDF(-3), DXPDF(4) = DXPDF(-4),                *
- *  DXPDF(5) = DXPDF(-5) and DXPDF(6) = DXPDF(-6) so far.              *
- *                                                                     *
- *     Please note that also PFTOPDG returns X * parton distribution   *
- *  function !                                                         *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *  2.4  Error Handling                                                *
- *                                                                     *
- *                                                                     *
- *     A note of caution should be addressed here. All structure       *
- *  functions are limited in X and Q**2 which range for most of the    *
- *  sets from                                                          *
- *                                                                     *
- *       4 - 10            < Q**2 < 10**6 - 10**8 (GeV/c)**2           *
- *       10**{-5} - 10**{-4} < X < 1.                                  *
- *                                                                     *
- *     If the user wants to extend the X or Q**2 ranges, it is pos-    *
- *  sible for some of the structure function sets, but the result      *
- *  should be looked at with great caution. If in doubt, please check  *
- *  with the authors of the structure function set. The program is     *
- *                                                                     *
- ******************************* - 14 - ********************************
-1***********************************************************************
- *                                                                     *
- *  protected against calculations of Q**2 values below Q2MIN (in this *
- *  case Q**2 = Q2MIN) and of X values in unphysical regions (X < 0 or *
- *  X > 1). The execution of the program is stopped in the later case. *
- *                                                                     *
- *     A COMMON block can be accessed                                  *
- *                                                                     *
- *       COMMON/W50513/ XMIN,XMAX,Q2MIN,Q2MAX                          *
- *                                                                     *
- *  where                                                              *
- *                                                                     *
- *       XMIN   = minimum allowed X value                              *
- *       XMAX   = maximum allowed X value                              *
- *       Q2MIN  = minimum allowed Q**2 value (in (GeV/c)**2)           *
- *       Q2MAX  = maximum allowed Q**2 value (in (GeV/c)**2)           *
- *                                                                     *
- *  is given for each set of structure functions. To control how often *
- *  the X or Q**2 ranges have been exceeded during execution of the    *
- *  user's program a call to the subroutine                            *
- *                                                                     *
- *       CALL PDFSTA                                                   *
- *                                                                     *
- *  at the termination phase of his/her MAIN program allows to print   *
- *  a summary of these errors. If the COMMON block                     *
- *                                                                     *
- *       COMMON/W50514W/ PDFWGT                                        *
- *                                                                     *
- *  where                                                              *
- *                                                                     *
- *       PDFWGT = weight                                               *
- *                                                                     *
- *  is filled by the user on an event-to-event basis the statistics    *
- *  in the error summary takes into account this weight.               *
- *                                                                     *
- *     The print flag in the COMMON/W50510/IFLPRT may be set to        *
- *  IFLPRT = 3 to print an error message EACH TIME a limit is exeeded  *
- *  during job execution. Please note that setting this flag may pro-  *
- *  duce an ENORMOUS AMOUNT of output!                                 *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *  2.5  Features                                                      *
- *                                                                     *
- *                                                                     *
- *     If NPTYPE = 1, NGROUP = 1 and NSET = 1, the user has the possi- *
- *  bility to use his/her own private set of structure functions       *
- *  (i.e. CDHS effective structure functions could be inseted here,    *
- *  or the user could implement a new set of structure functions for   *
- *  testing purposes). To insert the private set the DUMMY subroutine  *
- *                                                                     *
- *    SUBROUTINE STRPRIV(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL) *
- *                                                                     *
- *  should be activated by the user. All variables have to be defined  *
- *  as DOUBLE PRECISION. The variables have the same meaning as des-   *
- *  cribed for the subroutine STRUCTM in Section 2.2 (for explanations *
- *  see there).                                                        *
- *                                                                     *
- *                                                                     *
- ******************************* - 15 - ********************************
-1***********************************************************************
- *                                                                     *
- *     For the calculation of the strong coupling constant alpha(s)    *
- *  (see Section 3) with the user's private set of structure functions *
- *  the user should provide a value for the variables QCDL4 and QCDL5  *
- *  via the subroutine PDFSET at the initialization phase. In fact,    *
- *  because only the parameter QCDL5 is used in the alpha(s) calcul-   *
- *  ation, only this definition is mandatory.                          *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *  2.6  PION Structure Functions                                      *
- *                                                                     *
- *                                                                     *
- *     All the available sets of structure functions of the PION have  *
- *  been added to the package as well:                                 *
- *                                                                     *
- *                                                                     *
- *   -  the old Owens (OW-P) sets 1 and 2 [21],                        *
- *   -  the newer sets 1 to 3 of Aurenche et al. (ABFKW-P) [22],       *
- *   -  the very new Sutton, Martin, Roberts and Stirling (SMRS-P)     *
- *      sets 1 to 3 [23] and                                           *
- *   -  the very new Gl"uck, Reya and Vogt (GRV-P) [24] sets of LO     *
- *      and NLL calculations.                                          *
- *                                                                     *
- *     The calling sequence to the pion set of structure functions is  *
- *  kept identical to those described in Section 2.2 for the nucleon   *
- *  structure fuctions. It can be either                               *
- *                                                                     *
- *       CALL STRUCTM(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL)    *
- *                                                                     *
- *  or                                                                 *
- *                                                                     *
- *       CALL PFTOPDG(X,SCALE,DXPDF)                                   *
- *                                                                     *
- *  or the old version of the main steering routine                    *
- *                                                                     *
- *       CALL STRUCTF(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,TOP,GL)          *
- *                                                                     *
- *  where, as before, the INPUTs are to be provided by the user in     *
- *  DOUBLE PRECISION. The OUTPUT has the same definitions as before    *
- *  (see Section 2.2). In case of the pion structure functions it is   *
- *  also returned  X * parton distribution function with               *
- *  UPV = DNV, SEA = STR (USEA = DSEA = STR) and TOP = 0               *
- *  (DXPDF(1) = DXPDF(2), DXPDF(3) = DXPDF(-1) = DXPDF(-2) and         *
- *  DXPDF(6) = 0).                                                     *
- *                                                                     *
- *     A list of available sets of PION structure functions in the     *
- *  package is given in Table 2.                                       *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- ******************************* - 16 - ********************************
-1***********************************************************************
- *                                                                     *
- *  2.7  PHOTON Structure Functions                                    *
- *                                                                     *
- *                                                                     *
- *     In the previous version, all the available sets of structure    *
- *  functions of the PHOTON have been added to the package as well,    *
- *                                                                     *
- *   -  the old Duke and Owens (DO-G) LO and NLL sets [25],            *
- *   -  the old Drees and Grassie (DG-G) sets 1 to 4 [26],             *
- *   -  the new Abramowicz, Charchula and Levy (LAC-G) [27]            *
- *      sets 1 to 3,                                                   *
- *   -  the new Gordon and Storrow (GS-G) [28] sets, the LO sets 1     *
- *      and 2 and the NLL set,                                         *
- *   -  the new Gl"uck, Reya and Vogt (GRV-G) [29] sets, the LO, NLL   *
- *      and the leading term of the NLL calculations, and              *
- *   -  the new Aurenche et al. (ACFGP-G) [30] NLL sets 1 and 2,       *
- *      without and with massive charm.                                *
- *                                                                     *
- *     The calling sequence to the photon set of structure functions   *
- *  is also kept identical to those described in Section 2.2 for the   *
- *  nucleon structure fuctions. It can be either                       *
- *                                                                     *
- *       CALL STRUCTM(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL)    *
- *                                                                     *
- *  or                                                                 *
- *                                                                     *
- *       CALL PFTOPDG(X,SCALE,DXPDF)                                   *
- *                                                                     *
- *  or the old version of the main steering routine                    *
- *                                                                     *
- *       CALL STRUCTF(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,TOP,GL)          *
- *                                                                     *
- *  where, as before, the INPUTs are to be provided by the user in     *
- *  DOUBLE PRECISION. The OUTPUT has the same definitions as before    *
- *  (see Section 2.2). In case of the photon structure functions it    *
- *  is also returned  X * parton distribution function (note that the  *
- *  alpha(QED) has to be taken care of by the user !) with UPV = USEA, *
- *  DNV = DSEA and TOP = 0 (DXPDF(1) = 2 * DXPDF(-1), DXPDF(2) =       *
- *  2 * DXPDF(-2), as a result of the definitions - see Section 2.3 -, *
- *  and DXPDF(6) = 0), while all the other quarks are set identical to *
- *  their antiquark densities.                                         *
- *                                                                     *
- *     A list of available sets of PHOTON structure functions in the   *
- *  package is given in Table 3.                                       *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- ******************************* - 17 - ********************************
-1***********************************************************************
- *                                                                     *
- *  2.8  How to find PDFLIB                                            *
- *                                                                     *
- *                                                                     *
- *     At CERN the PDFLIB library and corresponding documentation are  *
- *  available on all systems in the standard areas.                    *
- *                                                                     *
- *    System    Library                     Documentation              *
- *                                                                     *
- *    CERNVM    PDFLIB TXTLIB Q             PDFLIB DOC Q               *
- *    VXCERN    CERN:[PRO.LIB]PDFLIB.OLB    CERN:[PRO.DOC]PDFLIB.PS    *
- *    Unix      /cern/pro/lib/libpdflib.a   /cern/pro/doc/pdf.ps       *
- *    systems                                                          *
- *                                                                     *
- *     The package may be accessed via the CERNLIB command:            *
- *                                                                     *
- *       CERNLIB PDFLIB                                                *
- *                                                                     *
- *                                                                     *
- *     Whenever you use this program, please give reference to the     *
- *  following papers:                                                  *
- *                                                                     *
- *       H. Plothow-Besch, 'PDFLIB: a library of all available parton  *
- *       density functions of the nucleon, the pion and the photon     *
- *       and the corresponding alpha(s) calculations',                 *
- *       Comp. Phys. Comm. 75 (1993) 396-416.                          *
- *                                                                     *
- *   and/or                                                            *
- *                                                                     *
- *       H. Plothow-Besch, 'PDFLIB: Structure Functions and            *
- *       alpha_s Calculations', User's Manual - Version 1.00,          *
- *       W5051 PDFLIB, 1991.03.21, CERN-PPE.                           *
- *                                                                     *
- *   and/or                                                            *
- *                                                                     *
- *       H. Plothow-Besch, 'PDFLIB: Nucleon, Pion and Photon Parton    *
- *       Density Functions and alpha_s Calculations',                  *
- *       User's Manual - Version 6.06, W5051 PDFLIB, 1995.03.15,       *
- *       CERN-PPE.                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *  3.  alpha(s) Calculation                                           *
- *                                                                     *
- *                                                                     *
- *     Within the same package a program is provided to calculate the  *
- *  strong coupling constant alpha(s) to second order as a function    *
- *  of Lambda(QCD) of five flavours and the desired number of fla-     *
- *  vours (NFL) for the selected set of structure functions, which     *
- *  fixes Lambda(QCD). The formula on which the calculations are       *
- *  based upon can be found in Ref. [33]. The same three parameters,   *
- *  NPTYPE, NGROUP and NSET, which select a structure function set,    *
- *  is used to steer the calculation of the alpha(s) value at a given  *
- *  scale from the Lambda(QCD) defined in the selected PDF set.        *
- *                                                                     *
- *      The value of alpha(s) is matched at the thresholds q = m_q.    *
- *  When invoked with NFL < 0, it chooses NFL as the number of fla-    *
- *                                                                     *
- ******************************* - 18 - ********************************
-1***********************************************************************
- *                                                                     *
- *  vours for which the masses are less then q. For the quark masses   *
- *  where thresholds are changed the following values have been used:  *
- *  m_charm = 1.5 GeV/c**2, m_bottom = 4.75 GeV/c**2 and               *
- *    m_top = 100 GeV/c**2.                                            *
- *                                                                     *
- *     The alpha(s) value can be obtained via a call to:               *
- *                                                                     *
- *       FUNCTION ALPHAS2(SCALE)                                       *
- *                                                                     *
- *     The user has to provide the following INPUT:                    *
- *                                                                     *
- *       SCALE  =  QCD scale in GeV                                    *
- *                                                                     *
- *     The function ALPHAS2 returns the following OUTPUT:              *
- *                                                                     *
- *       ALPHAS2 =  alpha strong to second order,                      *
- *                                                                     *
- *  if LO not equal to one. For those structure functions, for which   *
- *  the evolution is done to leading order only, alpha(s) to first     *
- *  order is returned.                                                 *
- *                                                                     *
- *     The same subroutine                                             *
- *                                                                     *
- *       SUBROUTINE PDFSET(PARM,VALUE)                                 *
- *                                                                     *
- *  as described in Section 2.1 should be used to change the default   *
- *  values of the variables NPTYPE, NGROUP, NSET, NFL, LO, TMAS,       *
- *  QCDL4, QCDL5, XMIN, XMAX, Q2MIN and Q2MAX.                         *
- *                                                                     *
- *     The internal COMMON block                                       *
- *                                                                     *
- *       COMMON/W50512/ QCDL4,QCDL5                                    *
- *                                                                     *
- *  where                                                              *
- *                                                                     *
- *       QCDL4  = QCD scale, Lambda_4(QCD), in GeV for four flavours   *
- *       QCDL5  = QCD scale, Lambda_5(QCD), in GeV for five flavours   *
- *                corresponding to QCDL4                               *
- *                                                                     *
- *  provides the actual value of Lambda_4(QCD) and Lambda_5(QCD) used  *
- *  in the alpha(s) calculation for four and five flavours, respectiv- *
- *  ly, for each set of structure functions.                           *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *  4.  Conclusions                                                    *
- *                                                                     *
- *                                                                     *
- *     In the Long Write-up as an example, the NUCLEON structure       *
- *  function distributions at the scale Q**2 = m_W**2 for u and d va-  *
- *  lence quarks as well as for the sea and gluon contributions are    *
- *  shown in Figs. 1 to 4 respectively, as a function of the parton X. *
- *  The average X values for CERN and FNAL as well as for LHC and SSC  *
- *  energies at that scale are indicated. Please, note the logarithmic *
- *  scale for the gluon distribution. The d/u ratio of the valence     *
- *  quarks is shown in Fig. 5. It can be seen from all these figures   *
- *                                                                     *
- ******************************* - 19 - ********************************
-1***********************************************************************
- *                                                                     *
- *  these figures that there is a large spread in shape for the dif-   *
- *  ferent sets in all contributions leading to different results in   *
- *  the cross section calculations. It should be pointed out that the  *
- *  spread from the different parametrisations in the SAME renormali-  *
- *  sation scheme is MUCH LARGER than the spread from parametrisations *
- *  in different renormalisation schemes. This is true for ALL parton  *
- *  densities, and in particular at small X values !                   *
- *                                                                     *
- *     A comparison to recent data is shown in Fig. 6, where the ratio *
- *  F_2(n)/F_2(p) at Q**2 = m_W**2 is displayed. This ratio mainly re- *
- *  flects the d/u ratio of valence quarks which is the dominant un-   *
- *  certainty of the W and Z boson cross section ratio in pp or pp_bar *
- *  reactions at lower energies. In fact, the older sets of structure  *
- *  functions, namely the set of Owens and Reya [2] and the set of     *
- *  Baier et al. [3] have only been implemented for completeness.      *
- *  These sets should not be used for cross section calculations of    *
- *  any hard processes at high energies anymore. From Fig. 6 we con-   *
- *  clude that also the other older sets of structure functions,       *
- *  namely the two sets of Duke and Owens [5] and the two sets of      *
- *  Eichten et al. [6] should be used with some care, because they do  *
- *  not fit the recent low energy deep-inelastic lepton-nucleon data   *
- *  from NMC and BCDMS [34]. This is not surprising because these data *
- *  were not available when the sets have been made. The new set of    *
- *  DO 1.1 superseeds the old DO sets 1 and 2 but still gives very     *
- *  limited results in the low X region, even though data are avai-    *
- *  lable there. (Please note that the last data point at the lowest   *
- *  X value shown in Fig. 6 is out of range in Q**2 for all sets of    *
- *  structure functions and that the point before last is very much at *
- *  the limit of the allowed range).                                   *
- *                                                                     *
- *     In Fig. 7 the F_2(p) distribution of the recent data from Ref.  *
- *  [35] at Q**2 = 5 GeV**2, together with the very recent data at     *
- *  Q**2 = 15 GeV**2 from the H1 and the ZEUS experiments [36] at the  *
- *  the ep-collider HERA is shown as a function of X. Overlayed on the *
- *  same Figure are a few theoretical predictions of more recent par-  *
- *  ton parametrisations [12,13,14,15,16,20] at Q**2 = 15 GeV**2.      *
- *  Please note that most of these parametrisations are NOT real pre-  *
- *  dictions because the data are fully or partially used in the fits. *
- *  Apart from Ref. [13] set B2, from Ref. [14] set HO, from Ref. [16] *
- *  set D-' and (H) and from Ref. [20] the new set 2pM, most of these  *
- *  more recent parametrisations fail to describe well enough the low  *
- *  X behaviour of the present data.                                   *
- *                                                                     *
- *     From Figs. 6 and 7 we conclude that the recent sets of parton   *
- *  densities, nameley the GRV set HO [14], but also the old MT set    *
- *  B2 [13] - both sets are real predictions -, and the updated MRS    *
- *  set (H) [16], could be preferred for all theoretical predictions   *
- *  involving structure functions, because they fit the present data   *
- *  best which, at the moment, still suffer from large uncertainties.  *
- *                                                                     *
- *     As an example, the u_bar, strange and gluon densities for the   *
- *  different sets of structure functions of the PHOTON are shown in   *
- *  Figs. 7, 8 and 9 of the Long Write-up, respectively, as a function *
- *  of X at the scale Q**2 = 100 (GeV/c)**2. In Fig. 10 a comparison   *
- *  of F_2(gamma)/alpha(QED) as a function of X at Q**2 = 100 (GeV/c)**2
- *  for the different photon structure function sets is made. The      *
- *                                                                     *
- ******************************* - 20 - ********************************
-1***********************************************************************
- *                                                                     *
- *  charm quark density has been taken into account. The theoretical   *
- *  predictions are compared with experimental data from the JADE      *
- *  experiment [37] at the same Q**2 value. Because of the large       *
- *  experimental uncertainties no distinction between the              *
- *  different parton density sets can be made yet.                     *
- *                                                                     *
- *     Coming new sets of structure functions could and should be      *
- *  easily implemented in this package. Authors of new sets are kindly *
- *  requested to provide us with the relevant information.             *
- *                                                                     *
- *     Please return any problems, questions, suggestions for im-      *
- *  provements to the author of the package                            *
- *  (e-mail address: PLOTHOW@CERNVM.CERN.CH).                          *
- *                                                                     *
- *     Requests for the Long Write-up of the User's Manual can be ad-  *
- *  dressed either to the author or to the CERN Program Library Office.*
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *  5.  Acknowledgements                                               *
- *                                                                     *
- *                                                                     *
- *     The author of the package would like to thank                   *
- *  M.A. Marquina/CERN-CN for help and support making this package     *
- *  available to the Physics community as part of the CERN Program     *
- *  Library. In addition we give many thanks to G. Folger/CERN-CN and  *
- *  I. McLaren/CERN-CN for thoroughly testing out the package on the   *
- *  different system areas.                                            *
- *                                                                     *
- *     We also would like to thank J.M. Marraffino/FNAL for the help-  *
- *  ful comments and the correction set which he supplied.             *
- *                                                                     *
- *     And finally we acknowledge K. Charchula/DESY who provided the   *
- *  source code of some of the photon structure functions which had    *
- *  previously been implemented in the package.                        *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- ******************************* - 21 - ********************************
-1***********************************************************************
- *                                                                     *
- *  Appendix A                                                         *
- *                                                                     *
- *                                                                     *
- *     As an example how to use the package, the FORTRAN code to ob-   *
- *  tain Fig. 1 of the Long Write-up which shows the parton density    *
- *  function distributions for the up valence quark as a function of   *
- *  the parton X for all sets of structure functions is given below:   *
- *                                                                     *
- *      PROGRAM PDFUPV                                                 *
- *  C  define the maximum number of PDF sets, the name and the cross reference
- *      PARAMETER (NPTYMX = 3, NGRMAX = 8, NSETMX = 41)                *
- *      COMMON /W505120/ NPGSMX(NPTYMX,NGRMAX),NSETFL(NPTYMX,NGRMAX,NSETMX)
- *      CHARACTER*8 SFNAME(NPTYMX,NGRMAX,NSETMX)                       *
- *      COMMON /W505110/ SFNAME                                        *
- *      COMMON /W505122/ MODECR(NPTYMX,NGRMAX,NSETMX)                  *
- *  C  define HBOOK settings                                           *
- *      PARAMETER (NHBMEM = 500000)                                    *
- *      COMMON /PAWC/ HMEM(NHBMEM)                                     *
- *      PARAMETER (NB=900, ID=100)                                     *
- *  C  define DOUBLE PRECISION variables for calling sequence to STRUCTM
- *      DOUBLE PRECISION DX,DSCALE,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL
- *      DOUBLE PRECISION ALF,ALPHAS2                                   *
- *      REAL    X, SCALE, UPV, DNV, USEA, DSEA, STR, CHM, BOT, TOP, GL *
- *      COMMON/W50516/ FIRST                                           *
- *      LOGICAL FIRST                                                  *
- *      CHARACTER*20 PARM(20)                                          *
- *      DOUBLE PRECISION VAL(20)                                       *
- *      DATA SCALE/80.140/                                             *
- *      DATA X00/0.002/,DX0/0.001/,XLOW/0.0015/,XUP/0.9015/            *
- *  C                                                                  *
- *      CALL HLIMIT(NHBMEM)                                            *
- *      DSCALE=SCALE                                                   *
- *  C  first call to PDFSET to initialize COMMON/W505120/              *
- *      PARM(1) = 'Init0'                                              *
- *      VAL(1)  = 0.D0                                                 *
- *      CALL PDFSET(PARM,VAL)                                          *
- *  C  loop over all existing sets of Nucleon structure functions (SF) *
- *      NPTYPE = 1                                                     *
- *      NHB = 0                                                        *
- *      DO 20 IGR = 1,NGRMAX                                           *
- *         IF(NPGSMX(NPTYPE,IGR).EQ.0) GOTO 20                         *
- *         DO ISET = 1,NPGSMX(NPTYPE,IGR)                              *
- *  C  book histograms for each set of SF separately                   *
- *         NHB = NHB + 1                                               *
- *         CALL HBOOK1(ID+NHB,'U Valence quark ',NB,XLOW,XUP,0.)       *
- *  C  force label printing for each set of SF (not only the 1st)      *
- *         FIRST = .TRUE.                                              *
- *  C  define and set parameters                                       *
- *         PARM(1) = 'Nptype'                                          *
- *         VAL(1)  = NPTYPE                                            *
- *         PARM(2) = 'Ngroup'                                          *
- *         VAL(2)  = IGR                                               *
- *         PARM(3) = 'Nset'                                            *
- *         VAL(3)  = ISET                                              *
- *         CALL PDFSET(PARM,VAL)                                       *
- *                                                                     *
- *                                                                     *
- ******************************* - 22 - ********************************
-1***********************************************************************
- *                                                                     *
- *  Example of Fig. 1, cont'd                                          *
- *                                                                     *
- *  C  loop over all X bins                                            *
- *        DO 10 I=1,NB                                                 *
- *           X = X00 + (I-1)*DX0                                       *
- *           DX = X                                                    *
- *           IF(X.LT.XLOW .OR. X.GT.XUP) GOTO 10                       *
- *           CALL STRUCTM(DX,DSCALE,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
- *           UPV=DUPV                                                  *
- *           IF(X.GT.0.499 .AND. X.LE.0.500) WRITE(6,1000) X,SCALE,UPV *
- *  1000     FORMAT(/,'  X= ',F6.4,'  Q= ',F6.3,'  UPV= ',8.4)         *
- *           CALL HF1(ID+NHB,X,UPV)                                    *
- *  10    CONTINUE                                                     *
- *  C  get alpha(s) for selected set of SF at Q = SCALE                *
- *        ALF = ALPHAS2(DSCALE)                                        *
- *        WRITE(6,4000) NPTYPE,IGR,ISET,SFNAME(NPTYPE,IGR,ISET),       *
- *   +                  MODECR(NPTYPE,IGR,ISET),NSF                    *
- *        WRITE(6,4001) DSCALE,ALF                                     *
- *  4000  FORMAT(/,'  Nptype = ',I1,'  Ngroup = ',I1,'  Nset = ',I2,   *
- *      +          '  Name = "'A8,'"  CrMode = ',I3,'  HBId = ',I3)    *
- *  4001  FORMAT(1H ,' SCALE  = ',F8.4,'   alpha(s) = ',F6             *
- *  C  get error summary for each set of SF                            *
- *        CALL PDFSTA                                                  *
- *     ENDDO                                                           *
- *  20 CONTINUE                                                        *
- *     STOP                                                            *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- ******************************* - 23 - ********************************
-1***********************************************************************
- *                                                                     *
- *     This code will produce the following output:                    *
- *                                                                     *
- *                                                                     *
- *         1. HBOOK Output :                                           *
- *                                                                     *
- *     one histogram for each set of structure functions. The histo-   *
- *  gram output has then been accessed, and each histogram has been    *
- *  superimposed on one single frame using the PAW package [38] to     *
- *  obtain Fig. 1 of the Long Write-up.                                *
- *                                                                     *
- *                                                                     *
- *         2. Print Output (Example is given for NPTYPE = 1,           *
- *  NGROUP = 3, NSET = 41 only :                                       *
- *                                                                     *
- * ___________________________________________________________________ *
- *|                                                                   |*
- *| 1*****  PDFLIB Version:  6.06  Released on  950315  at  17.05  in |*
- *| the CERN Program Library     W5051  *****                         |*
- *|  *****  Library compiled on      950315  at        2337   *****   |*
- *|                                                                   |*
- *| Nucleon PDFs:  Ngroup = 3,  Nset = 41,  for MRS Set (G) (L255-MSb)|*
- *| ------------------------------------------------------------------|*
- *|                                             Structure Functions   |*
- *|                                             -------------------   |*
- *|                                                                   |*
- *| X= 0.5000  Q= 80.140  UPV=   0.1176                               |*
- *|                                                                   |*
- *| Nptype = 1  Ngroup = 3  Nset = 41 Name = "MRS-G   "  CrMode =  -1  HBId = 60
- *| Scale  =  80.1400  alpha(s) =  0.1154                             |*
- *|                                                                   |*
- *| PDFLIB : Summary from PDFSTA                                      |*
- *| Nptype = 1  Ngroup = 3  Nset = 41  Name = "MRS-G   "  CrMode =  -1|*
- *| Nfl    = -5  LO = 2  Tmas = 180.00 GeV/c**2                       |*
- *| QCDL4  = 0.2550 GeV,  QCDL5 =   0.1708 GeV                        |*
- *| Xmin   = 0.10E-04,  Xmax =  0.99999E+00,  Q2min =  5.000 (GeV/c)**2,
- *|   Q2max =  0.27E+10 (GeV/c)**2                                    |*
- *|                                                                   |*
- *| PDFSTA:    NO errors occured                                      |*
- *|___________________________________________________________________|*
- *                                                                     *
- *                                                                     *
- *     Apart from the two top lines the print output shown above is    *
- *  repeated for each set of structure functions.                      *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- ******************************* - 24 - ********************************
-1***********************************************************************
- *                                                                     *
- *  Appendix B                                                         *
- *                                                                     *
- *                                                                     *
- *     A list of subroutines, functions, COMMON blocks and parameters  *
- *  of RELEVANCE to the user which are used in PDFLIB is given below:  *
- *                                                                     *
- *                                                                     *
- *     List of relevant subprograms and COMMON blocks in PDFLIB:       *
- *                                                                     *
- *  _________________________________________________________________  *
- *  |                                                               |  *
- *  |  (S = Subroutine, F = Function, C = COMMON, P = Parameter)    |  *
- *  |                                                               |  *
- *  | ------------------------------------------------------------- |  *
- *  |                                                               |  *
- *  |  S  PDFSET   : to set all parameters for PDFLIB               |  *
- *  |  S  STRUCTF  : to access PDFLIB with output in PDFLIB format  |  *
- *  |  S  STRUCTM  : to access PDFLIB with output in PDFLIB format  |  *
- *  |                (with u_bar not equal to d_bar)                |  *
- *  |  S  PFTOPDG  : to access PDFLIB with output in PDG format     |  *
- *  |  S  PDFSTA   : to print summary of parameters and error statistics
- *  |  F  ALPHAS2  : to calculate alpha(s) to second order QCD      |  *
- *  |                                                               |  *
- *  |  C  W50510   : IFLPRT - to set print flag                     |  *
- *  |  C  W50511   : NPTYPE,NGROUP,NSET,MODE,NFL,LO,TMAS            |  *
- *  |                - to select a structure function set, to       |  *
- *  |                define number of flavours, flag for leading    |  *
- *  |                order and value for the mass of the top quark  |  *
- *  |  C  W50512   : QCDL4,QCDL5 - contains Lambda_4(QCD) and       |  *
- *  |                Lambda_5(QCD) (in GeV)                         |  *
- *  |  C  W50513   : XMIN,XMAX,Q2MIN,Q2MAX - contains minimum and   |  *
- *  |                maximum values of X and Q**2                   |  *
- *  |  C  W50516   : FIRST - to force printing name of structure    |  *
- *  |                function set in case of multiple choices       |  *
- *  |  C  W50519   : NEWVER - flags NEW/OLD version of PDFLIB format|  *
- *  |  C  W505110  : SFNAME - contains NAME (CHARACTER*8 array with |  *
- *  |                dimension SFNAME(NPTYMX,NGRMAX,NSETMX)) for    |  *
- *  |                each set of PDFs                               |  *
- *  |  C  W505120  : NPGSMX(NPTYMX,NGRMAX),NSETFL(NPTYMX,NGRMAX,NSETMX)*
- *  |                - contains maximum number of stucture functions|  *
- *  |                per particle type and author group, and FLAG   |  *
- *  |                for retracted PDF set (NSETFL(I,J,K) = 0),     |  *
- *  |  C  W505121  : NPTYCR(MODEMX),NGROCR(MODEMX),NSETCR(MODEMX)   |  *
- *  |                - cross reference to the NPTYPE, NGROUP and    |  *
- *  |                NSET value (NEW format) for a given MODE       |  *
- *  |                value (OLD format)                             |  *
- *  |  C  W505122  : MODECR(NPTYMX,NGRMAX,NSETMX) - cross reference |  *
- *  |                to the MODE value (OLD format) for a given     |  *
- *  |                NPTYPE, NGROUP and NSET value (NEW format)     |  *
- *  |  P  W5051P2  : MODEMX - maximum number of structure function  |  *
- *  |                sets in the OLD PDFLIB format                  |  *
- *  |              : NPTYMX,NGRMAX,NSETMX - maximum number of       |  *
- *  |                structure function sets in the NEW PDFLIB      |  *
- *  |                format for particle type, author group number  |  *
- *  |                and number of sets within an author group      |  *
- *  |  P  W5051P7  : L6 - to define logical print unit number       |  *
- *  |_______________________________________________________________|  *
- ******************************* - 25 - ********************************
-1***********************************************************************
- *                                                                     *
- *  Appendix C                                                         *
- *                                                                     *
- *                                                                     *
- *     A complete list of subroutines, functions, COMMON blocks and    *
- *  parameters which are used in PDFLIB is given below:                *
- *                                                                     *
- *                                                                     *
- *     Complete list of subprograms and COMMON blocks in PDFLIB:       *
- *                       (in alphabetical order)                       *
- *                                                                     *
- *  _________________________________________________________________  *
- *  |                                                               |  *
- *  |  (S = Subroutine, F = Function, C = COMMON, P = Parameter)    |  *
- *  |                                                               |  *
- *  | ------------------------------------------------------------- |  *
- *  |                                                               |  *
- *  |  S  ABFKW1   : ABFKW SF set 1 in the pion (NP=2,NG=6,NS=1)    |  *
- *  |  S  ABFKW2   : ABFKW SF set 2 in the pion (NP=2,NG=6,NS=2)    |  *
- *  |  S  ABFKW3   : ABFKW SF set 3 in the pion (NP=2,NG=6,NS=3)    |  *
- *  |  S  ACFGP1   : ACFGP SF set 1 in the photon (NP=3,NG=6,NS=1)  |  *
- *  |  S  ACFGP2   : ACFGP SF set 1 in the photon (mc) (NP=3,NG=6,NS=2)*
- *  |  F  AFCPLU   : to get PD with massive charm for ACFGP         |  *
- *  |  F  AFGETFV  : to interpolate PDs in grid for ACFGP           |  *
- *  |  F  AFPOLIN  : to interpolate PDs in grid for ACFGP           |  *
- *  |  F  AFRATIN  : to interpolate PDs in grid for ACFGP           |  *
- *  |  F  ALPHAS2  : to calculate alpha(s) to second order QCD      |  *
- *  |  F  AUGETFV  : to interpolate logarithmically PDs in grid     |  *
- *  |  S  AURGAM   : to get the PDs in the photon for ACFGP-G       |  *
- *  |  S  AURPI1   : to get the PDs in the pion for ABFKW-P set 1   |  *
- *  |  S  AURPI2   : to get the PDs in the pion for ABFKW-P set 2   |  *
- *  |  S  AURPI3   : to get the PDs in the pion for ABFKW-P set 3   |  *
- *  |  S  AURPR    : to get the PDs in the nucleon for ABFOW        |  *
- *  |  F  BBETA    : DOUBLE PRECISION beta function using DGAMMA    |  *
- *  |  F  BETA     : SINGLE PRECISION beta function using GAMMA     |  *
- *  |  S  BXF3     : BEP SF set for up and down valence quarks      |  *
- *  |  S  BXGLUE   : BEP SF set for gluons                          |  *
- *  |  S  BXQBAR   : BEP SF set for antiquarks                      |  *
- *  |  F  CTQ1PD   : CTEQ SF sets 1M, 1MS, 1ML, 1D and 1L           |  *
- *  |  F  CTQ1PF   : Description of CTEQ sets 1M, 1MS, 1ML, 1D, 1L  |  *
- *  |  F  CTQ2OPD  : CTEQ SF sets 2M, 2MS, 2ML, 2D and 2L           |  *
- *  |  F  CTQ2OPF  : Description of CTEQ sets 2M, 2MS, 2ML, 2D, 2L  |  *
- *  |  S  CTQ2OPS  : CTEQ SF sets 2M, 2MS, 2ML, 2D and 2L           |  *
- *  |  F  CTQ2PPD  : CTEQ SF sets 2pM, 2pMS, 2pML, 2pD and 2pL      |  *
- *  |  F  CTQ2PPF  : Description of CTEQ sets 2pM, 2pMS, 2pML, 2pD, 2pL*
- *  |  S  CTQ2PPS  : CTEQ SF sets 2pM, 2pMS, 2pML, 2pD and 2pL      |  *
- *  |  F  DBFINT   : to interpolate linearly PDs in grid (double precision)
- *  |  S  DFLM1    : interface to DFLM 160 (NP=1,NG=2,NS=7)         |  *
- *  |  S  DFLM2    : interface to DFLM 260 (NP=1,NG=2,NS=8)         |  *
- *  |  S  DFLM3    : interface to DFLM 360 (NP=1,NG=2,NS=9)         |  *
- *  |  S  DFLM4    : interface to DFLM soft valence quarks (NP=1,NG=2,NS=1)
- *  |  S  DFLM5    : interface to DFLM hard valence quarks (NP=1,NG=2,NS=2)
- *  |  S  DFLM6    : interface to DFLM soft gluons (NP=1,NG=2,NS=3) |  *
- *  |  S  DFLM7    : interface to DFLM hard gluons (NP=1,NG=2,NS=4) |  *
- *  |  S  DFLM8    : interface to DFLM LO average fit (NP=1,NG=2,NS=5) *
- *  |  S  DFLM9    : interface to DFLM NLL average fit (NP=1,NG=2,NS=6)*
- *  |  S  DGPHO1   : DG SF set 1 in the photon (NP=3,NG=2,NS=1)     |  *
- *  |_______________________________________________________________|  *
- ******************************* - 26 - ********************************
-1***********************************************************************
- *  _________________________________________________________________  *
- *  |  Complete list of subprograms, cont'd                         |  *
- *  |                                                               |  *
- *  |  S  DGPHO2   : DG SF set 2 in the photon (NP=3,NG=2,NS=2)     |  *
- *  |  S  DGPHO3   : DG SF set 3 in the photon (NP=3,NG=2,NS=3)     |  *
- *  |  S  DGPHO4   : DG SF set 4 in the photon (NP=3,NG=2,NS=4)     |  *
- *  |  S  DOPHO1   : DO SF set 1 in the photon (NP=3,NG=1,NS=1)     |  *
- *  |  S  DOPHO2   : DO SF set 2 in the photon (NP=3,NG=1,NS=2)     |  *
- *  |  S  FXAVER   : DFLM SF set of LO central average fit          |  *
- *  |  S  FXG160   : DFLM SF set of NLL with QCDL4 = 160            |  *
- *  |  S  FXG260   : DFLM SF set of NLL with QCDL4 = 260            |  *
- *  |  S  FXG360   : DFLM SF set of NLL with QCDL4 = 360            |  *
- *  |  S  FXGHAR   : DFLM SF set of LO hard gluon                   |  *
- *  |  S  FXGSOF   : DFLM SF set of LO soft gluon                   |  *
- *  |  S  FXNLLA   : DFLM SF set of NLL central average fit         |  *
- *  |  S  FXVHAR   : DFLM SF set of LO hard valence                 |  *
- *  |  S  FXVSOF   : DFLM SF set of LO soft valence                 |  *
- *  |  F  GAMFUN   : to calculate GAMMA funct. in either SP or DP prec.*
- *  |  F  GAMMADO  : to calculate GAMMA funct. in SP (Hastings)     |  *
- *  |  F  GHRDV    : GHR SF set for down quarks                     |  *
- *  |  F  GHRGL    : GHR SF set for gluons                          |  *
- *  |  F  GHRUV    : GHR SF set for up quarks                       |  *
- *  |  F  GHRXI    : GHR SF set for antiquarks                      |  *
- *  |  F  GHRXS    : GHR SF set for strange quarks                  |  *
- *  |  F  GRVFGP   : GRV SF set for gluons in the pion              |  *
- *  |  F  GRVFQBP  : GRV SF set for s, c and b quarks in the pion   |  *
- *  |  F  GRVFV    : GRV SF set for up + down valence quarks        |  *
- *  |  F  GRVFVP   : GRV SF set for valence quarks in the pion      |  *
- *  |  F  GRVFW    : GRV SF set for antiquarks and gluons           |  *
- *  |  F  GRVFWS   : GRV SF set for strange, charm and bottom quarks|  *
- *  |  S  GRVGAH0  : GRV SF set lead. terms of NLL in the photon (NP=3,NG=5,NS=1
- *  |  S  GRVGAHO  : GRV SF set of NLL in the photon (NP=3,NG=5,NS=2)  *
- *  |  S  GRVGALO  : GRV SF set of LO in the photon (NP=3,NG=5,NS=3)|  *
- *  |  F  GRVGF    : GRV SF set for u and d quarks and gluons in the photon
- *  |  F  GRVGFS   : GRV SF set for s, c and b quarks in the photon |  *
- *  |  S  GRVHO    : GRV SF set of NLL (NP=1,NG=5,NS=3)             |  *
- *  |  S  GRVLO    : GRV SF set of LO (NP=1,NG=5,NS=4)              |  *
- *  |  S  GRVPIHO  : GRV SF set of NLL in the pion (NP=2,NG=5,NS=1) |  *
- *  |  S  GRVPILO  : GRV SF set of LO in the pion (NP=2,NG=5,NS=2)  |  *
- *  |  S  GSXCOR   : to get X coordinates for SFGSHO + LO sets 1, 2 |  *
- *  |  S  HMRS1EB  : Description of HMRS sets E, B (1.90-retracted) |  *
- *  |  S  HMRS2EB  : Description of HMRS sets E, B (3.90)           |  *
- *  |  S  HMRS3EB  : Description of HMRS sets E, B (4.90)           |  *
- *  |  S  HMRSGEB  : Description of HMRS sets E, E+, E- and B       |  *
- *  |  S  KMRSEB   : Description of KMRS sets E and B               |  *
- *  |  S  KMRSGEB  : Description of KMRS B0 with L=135,160,200,235  |  *
- *  |  S  MRSDSHP  : Description of MRS sets D0', S0', D-', (H)     |  *
- *  |  S  MRSEB    : Description of MRS sets E and B                |  *
- *  |  S  MRSEBP   : Description of MRS sets E' and B'              |  *
- *  |  S  NEWDO1   : DO SF set 1.1 (NP=1,NG=1,NS=10)                |  *
- *  |  S  PDFSET   : to set all parameters for PDFLIB               |  *
- *  |  S  PDFSTA   : to print summary of param. and error statistics|  *
- *  |  S  PDFVERS  : to get date and time of PDFLIB Library release |  *
- *  |  F  PDXMT    : Description of MT sets 1 to 11                 |  *
- *  |  S  PDZXMT   : MT SF sets 1 to 11                             |  *
- *  |  S  PFTOPDG  : to access PDFLIB with output in PDG format     |  *
- *  |  S  PHLAC1   : interface to LAC-G set 1 in the photon         |  *
- *  |_______________________________________________________________|  *
- ******************************* - 27 - ********************************
-1***********************************************************************
- *  _________________________________________________________________  *
- *  |  Complete list of subprograms, cont'd                         |  *
- *  |                                                               |  *
- *  |  S  PHLAC2   : interface to LAC-G set 2 in the photon         |  *
- *  |  S  PHLAC3   : interface to LAC-G set 3 in the photon         |  *
- *  |  S  RXDV     : OR SF set for down valence quarks              |  *
- *  |  S  RXGLUE   : OR SF set for gluons                           |  *
- *  |  S  RXQBAR   : OR SF set for antiquarks                       |  *
- *  |  S  RXUD     : OR SF set for the sum of up + down quarks      |  *
- *  |  S  SEADIS   : BEBC SF set for antiquarks                     |  *
- *  |  S  SFABFOW  : ABFOW SF set (NP=1,NG=6,NS=1)                  |  *
- *  |  S  SFBMA    : BM SF set A (NP=1,NG=7,NS=1)                   |  *
- *  |  S  SFBMB    : BM SF set B (NP=1,NG=7,NS=2)                   |  *
- *  |  S  SFCTQ11  : CTEQ SF set 1L (LO) (NP=1,NG=4,NS=12)          |  *
- *  |  S  SFCTQ12  : CTEQ SF set 1M (NP=1,NG=4,NS=13)               |  *
- *  |  S  SFCTQ13  : CTEQ SF set 1MS (NP=1,NG=4,NS=14)              |  *
- *  |  S  SFCTQ14  : CTEQ SF set 1ML (NP=1,NG=4,NS=15)              |  *
- *  |  S  SFCTQ15  : CTEQ SF set 1D (NP=1,NG=4,NS=16)               |  *
- *  |  S  SFCTQ21  : CTEQ SF set 2L (LO) (NP=1,NG=4,NS=17)          |  *
- *  |  S  SFCTQ22  : CTEQ SF set 2M (NP=1,NG=4,NS=18)               |  *
- *  |  S  SFCTQ23  : CTEQ SF set 2MS (NP=1,NG=4,NS=19)              |  *
- *  |  S  SFCTQ24  : CTEQ SF set 2MF (NP=1,NG=4,NS=20)              |  *
- *  |  S  SFCTQ25  : CTEQ SF set 2ML (NP=1,NG=4,NS=21)              |  *
- *  |  S  SFCTQ26  : CTEQ SF set 2D (NP=1,NG=4,NS=22)               |  *
- *  |  S  SFDISD0  : MRS SF set D0' (DIS) (NP=1,NG=3,NS=33)         |  *
- *  |  S  SFDISDM  : MRS SF set D-' (DIS) (NP=1,NG=3,NS=34)         |  *
- *  |  S  SFDISH   : MRS SF set (H) (NP=1,NG=3,NS=36)               |  *
- *  |  S  SFDISS0  : MRS SF set S0' (DIS) (NP=1,NG=3,NS=32)         |  *
- *  |  S  SFEHLQ1  : EHLQ SF set 1 (NP=1,NG=1,NS=8)                 |  *
- *  |  S  SFEHLQ2  : EHLQ SF set 2 (NP=1,NG=1,NS=9)                 |  *
- *  |  S  SFGRVH   : old GRV SF set of NLL (NP=1,NG=5,NS=1)         |  *
- *  |  S  SFGRVL   : old GRV SF set of LO (NP=1,NG=5,NS=2)          |  *
- *  |  S  SFGSHO   : GS SF set of NLL in the photon (NP=3,NG=4,NS=1)|  *
- *  |  S  SFGSLO1  : GS SF set 1 of LO in the photon (NP=3,NG=4,NS=2)  *
- *  |  S  SFGSLO2  : GS SF set 2 of LO in the photon (NP=3,NG=4,NS=3)  *
- *  |  S  SFKBMR2  : KMRS SF set B- with R2 shadowing (NP=1,NG=3,NS=19)*
- *  |  S  SFKBMR5  : KMRS SF set B- with R5 shadowing (NP=1,NG=3,NS=20)*
- *  |  S  SFLACG1  : LAC SF set 1 in the photon (NP=3,NG=3,NS=1)    |  *
- *  |  S  SFLACG2  : LAC SF set 2 in the photon (NP=3,NG=3,NS=2)    |  *
- *  |  S  SFLACG3  : LAC SF set 3 in the photon (NP=3,NG=3,NS=3)    |  *
- *  |  S  SFMRSD0  : MRS SF set D0 (NP=1,NG=3,NS=27)                |  *
- *  |  S  SFMRSDM  : MRS SF set D- (NP=1,NG=3,NS=28)                |  *
- *  |  S  SFMRSS0  : MRS SF set S0 (NP=1,NG=3,NS=26)                |  *
- *  |  S  SFMSBB0  : MRS SF set B0' (NP=1,NG=3,NS=                  |  *
- *  |  S  SFMSBD0  : MRS SF set D0' (NP=1,NG=3.NS=30)               |  *
- *  |  S  SFMSBDM  : MRS SF set D-' (NP=1,NG=3,NS=31)               |  *
- *  |  S  SFMSBH   : MRS SF set (H) (NP=1,NG=3,NS=35)               |  *
- *  |  S  SFMSBS0  : MRS SF set S0' (NP=1,NG=3,NS=29)               |  *
- *  |  S  SFMT10   : interface to MT set 6 (MS_bar)                 |  *
- *  |  S  SFMT11   : interface to MT set LO                         |  *
- *  |  S  SFMTU1   : interface to MT set S1 (DIS)                   |  *
- *  |  S  SFMTU2   : interface to MT set B1 (DIS)                   |  *
- *  |  S  SFMTU3   : interface to MT set B2 (DIS)                   |  *
- *  |  S  SFMTU4   : interface to MT set E1 (DIS)                   |  *
- *  |  S  SFMTU5   : interface to MT set 6 (u_bar + d_bar)/2 (DIS)  |  *
- *  |  S  SFMTU6   : interface to MT set S1 (MS_bar)                |  *
- *  |  S  SFMTU7   : interface to MT set B1 (MS_bar)                |  *
- *  |_______________________________________________________________|  *
- ******************************* - 28 - ********************************
-1***********************************************************************
- *  _________________________________________________________________  *
- *  |  Complete list of subprograms, cont'd                         |  *
- *  |                                                               |  *
- *  |  S  SFMTU8   : interface to MT set B2 (MS_bar)                |  *
- *  |  S  SFMTU9   : interface to MT set E1 (MS_bar)                |  *
- *  |  S  SMRSP31  : SMRS SF set 1 in the pion (NP=2,NG=3,NS=1)     |  *
- *  |  S  SMRSP32  : SMRS SF set 2 in the pion (NP=2,NG=3,NS=2)     |  *
- *  |  S  SMRSP33  : SMRS SF set 3 in the pion (NP=2,NG=3,NS=3)     |  *
- *  |  S  SMRSPI   : Description of SMRS sets 1, 2, 3 in the pion   |  *
- *  |  S  SRCTQ21  : CTEQ SF set 2pL (LO) (NP=1,NG=4,NS=23)         |  *
- *  |  S  SRCTQ22  : CTEQ SF set 2pM (NP=1,NG=4,NS=24)              |  *
- *  |  S  SRCTQ23  : CTEQ SF set 2pMS (NP=1,NG=4,NS=25)             |  *
- *  |  S  SRCTQ24  : CTEQ SF set 2pMF (NP=1,NG=4,NS=26)             |  *
- *  |  S  SRCTQ25  : CTEQ SF set 2pML (NP=1,NG=4,NS=27)             |  *
- *  |  S  SRCTQ26  : CTEQ SF set 2pD (NP=1,NG=4,NS=28)              |  *
- *  |  S  STRBEBC  : BEBC SF set (NP=1,NG=1,NS=2)                   |  *
- *  |  S  STRBEP   : BEP SF set (NP=1,NG=1,NS=4)                    |  *
- *  |  S  STRCBP   : MRS SF set B' (NP=1,NG=3,NS=7)                 |  *
- *  |  S  STRCEP   : MRS SF set E' (NP=1,NG=3,NS=6)                 |  *
- *  |  S  STRDO1   : DO SF set 1 (NP=1,NG=1,NS=6)                   |  *
- *  |  S  STRDO2   : DO SF set 2 (NP=1,NG=1,NS=7)                   |  *
- *  |  S  STRGHR   : GHR SF set (NP=1,NG=1,NS=5)                    |  *
- *  |  S  STRH1B   : HMRS SF set B (NP=1,NG=3,NS=9 - retracted)     |  *
- *  |  S  STRH1E   : HMRS SF set E (NP=1,NG=3,NS=8 - retracted)     |  *
- *  |  S  STRH2B   : HMRS SF set B (NP=1,NG=3,NS=11)                |  *
- *  |  S  STRH2E   : HMRS SF set E (NP=1,NG=3,NS=10)                |  *
- *  |  S  STRH3B   : HMRS SF set B (NP=1,NG=3,NS=15)                |  *
- *  |  S  STRH3B1  : HMRS SF set B - L=100 (NP=1,NG=3,NS=16)        |  *
- *  |  S  STRH3B3  : HMRS SF set B - L=300 (NP=1,NG=3,NS=17)        |  *
- *  |  S  STRH3E   : HMRS SF set E (NP=1,NG=3,NS=14)                |  *
- *  |  S  STRH3EM  : HMRS SF set E- (NP=1,NG=3,NS=13)               |  *
- *  |  S  STRH3EP  : HMRS SF set E+ (NP=1,NG=3,NS=12)               |  *
- *  |  S  STRKB0   : KMRS SF set B0 (NP=1,NG=3,NS=21), updated version *
- *  |  S  STRKB02  : KMRS SF set B0 - L=135 (NP=1,NG=3,NS=22)       |  *
- *  |  S  STRKB03  : KMRS SF set B0 - L=160 (NP=1,NG=3,NS=23)       |  *
- *  |  S  STRKB04  : KMRS SF set B0 - L=200 (NP=1,NG=3,NS=24)       |  *
- *  |  S  STRKB05  : KMRS SF set B0 - L=235 (NP=1,NG=3,NS=25)       |  *
- *  |  S  STRKB0O  : KMRS SF set B0 (old version)                   |  *
- *  |  S  STRKBM   : KMRS SF set B- (NP=1,NG=3,NS=18)               |  *
- *  |  S  STRMRS1  : MRS SF set 1 (NP=1,NG=3,NS=1)                  |  *
- *  |  S  STRMRS2  : MRS SF set 2 (NP=1,NG=3,NS=2)                  |  *
- *  |  S  STRMRS3  : MRS SF set 3 (NP=1,NG=3,NS=3)                  |  *
- *  |  S  STROWP1  : OW SF set 1 in the pion (NP=2,NG=1,NS=1)       |  *
- *  |  S  STROWP2  : OW SF set 2 in the pion (NP=2,NG=1,NS=2)       |  *
- *  |  S  STRPRIV  : Dummy subroutine for a new set of SF           |  *
- *  |  S  STRUCB   : MRS SF set B (NP=1,NG=3,NS=4)                  |  *
- *  |  S  STRUCE   : MRS SF set E (NP=1,NG=3,NS=5)                  |  *
- *  |  S  STRUCOR  : OR SF set (NP=1,NG=1,NS=3)                     |  *
- *  |  S  STRUCTF  : to access PDFLIB with output in PDFLIB format  |  *
- *  |  S  STRUCTM  : to access PDFLIB with output in PDFLIB format  |  *
- *  |     TESTGAM  : main program to test PDFLIB for photons        |  *
- *  |     TESTGAMO : main program to test PDFLIB for photons (OLD)  |  *
- *  |     TESTPDF  : main program to test PDFLIB in PDFLIB format   |  *
- *  |     TESTPDFO : main program to test PDFLIB in PDFLIB format (OLD)*
- *  |     TESTPDG  : main program to test PDFLIB in PDG format      |  *
- *  |     TESTPDGO : main program to test PDFLIB in PDG format (OLD)|  *
- *  |     TMAINUPV : main program to show how to use the package    |  *
- *  |_______________________________________________________________|  *
- ******************************* - 29 - ********************************
-1***********************************************************************
- *  _________________________________________________________________  *
- *  |  Complete list of subprograms, cont'd                         |  *
- *  |                                                               |  *
- *  |     TMAINUPVO: main program to show how to use the package (OLD) *
- *  |  F  V3       : BEBC SF set to calculate X * VALDIS            |  *
- *  |  S  VALDIS   : BEBC SF set for up + down valence quarks       |  *
- *  |  S  VLAMBD   : Entry in PDZXMT to get Lambda value for MT sets|  *
- *  |  S  WATE32   : defines 32 point Gaussian quadrature routine   |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |                                                               |  *
- *  |_______________________________________________________________|  *
- ******************************* - 30 - ********************************
-1***********************************************************************
- *  _________________________________________________________________  *
- *  |  Complete list of subprograms, cont'd                         |  *
- *  |                                                               |  *
- *  |  C  W50510   : IFLPRT - to set print flag                     |  *
- *  |  C  W50511   : NPTYPE,NGROUP,NSET,MODE,NFL,LO,TMAS -          |  *
- *  |                to select structure function set, to define    |  *
- *  |                number of flavours, flag for leading order     |  *
- *  |                and value for the mass of the top quark        |  *
- *  |  C  W50512   : QCDL4,QCDL5 - contains Lambda_4(QCD) and       |  *
- *  |                Lambda_5(QCD) (in GeV)                         |  *
- *  |  C  W50513   : XMIN,XMAX.Q2MIN,Q2MAX - contains minimum       |  *
- *  |                and maximum values of X and Q**2               |  *
- *  |  C  W50514   : WXMIN,WXMAX,WQ2MIN,WQ2MAX,WTXMIN,WTXMAX,WQ2MIN,|  *
- *  |                WQ2MAX - error counting for fractional and     |  *
- *  |                total number of X and Q**2                     |  *
- *  |  C  W50514W  : PDFWGT - to transmit user filled event weight  |  *
- *  |  C  W50515   : IFLSET,IFLSTA - to transmit flag for call to   |  *
- *  |                subroutine PDFSET or PDFSTA                    |  *
- *  |  C  W50516   : FIRST - to force printing name of structure    |  *
- *  |                function set in case of multiple choices       |  *
- *  |  C  W50517   : N6 - contains logical print unit number        |  *
- *  |  C  W50518   : save of COMMONs W50510, W50511, W50512, W50513 |  *
- *  |  C  W50519   : NEWVER - to flag NEW or OLD version of PDFLIB  |  *
- *  |                format                                         |  *
- *  |  C  W505110  : SFNAME - contains NAME (CHARACTER*8 array with |  *
- *  |                dimension SFNAME(NPTYMX,NGRMAX,NSETMX)) for    |  *
- *  |                each set of PDFs                               |  *
- *  |  C  W505120  : NPGSMX(NPTYMX,NGRMAX),NSETFL(NPTYMX,NGRMAX,NSETMX)*
- *  |                - to provide the actual maximum number of stuc-|  *
- *  |                ture function sets per particle type and author|  *
- *  |                group, NSETFL(I,J,K) = 0 for retracted PDF set |  *
- *  |  C  W505121  : NPTYCR(MODEMX),NGROCR(MODEMX),NSETCR(MODEMX)   |  *
- *  |                - to provide the cross reference to the        |  *
- *  |                NPTYPE, NGROUP and NSET value (NEW format) for |  *
- *  |                a given MODE value (OLD format)                |  *
- *  |  C  W505122  : MODECR(NPTYMX,NGRMAX,NSETMX) - to provide the  |  *
- *  |                cross reference to the MODE value (OLD format) |  *
- *  |                for a given NPTYPE, NGROUP and NSET value      |  *
- *  |                (NEW format)                                   |  *
- *  |  P  W5051P0  : LFLPRT - to define print level (D=0)           |  *
- *  |  P  W5051P1  : LPTYPE,LGROUP,LNSET,IFL,LORD,TM - to define    |  *
- *  |                the defaults for a PDF set given by the        |  *
- *  |                particle type, author group and number of set  |  *
- *  |                within an author group, defaults for number of |  *
- *  |                flavours, leading order flag and mass of the   |  *
- *  |                top quark (D=1, 5, 3, 5, 2, 100.D0)            |  *
- *  |  P  W5051P2  : MODEMX - maximum number of structure function  |  *
- *  |                sets in the OLD PDFLIB format                  |  *
- *  |              : NPTYMX,NGRMAX,NSETMX - maximum number of       |  *
- *  |                structure function sets in the NEW PDFLIB      |  *
- *  |                format for particle type, author group number  |  *
- *  |                and number of sets within an author group      |  *
- *  |  P  W5051P2  : MODEMX - to define maximum number of SF sets   |  *
- *  |                in PDFLIB (OLD format) (D=281)                 |  *
- *  |  P  W5051P7  : L6 - to define logical print unit number (D=6) |  *
- *  |  P  W5051P7  : L6 - to define logical print unit number       |  *
- *  |     W50511C  : Comments only - complete list of available SFs |  *
- *  |_______________________________________________________________|  *
- ******************************* - 31 - ********************************
-1***********************************************************************
- *                                                                     *
- *  References                                                         *
- *                                                                     *
- *                                                                     *
- *  [1] A.J.Buras and K.J.F. Gaemers : Nucl. Phys. B132 (1978) 249;    *
- *      BEBC WA 59 Collaboration, K. Varnell et al., Z. Phys. C36 (1987) 1.
- *                                                                     *
- *  [2] J.F.Owens and E. Reya : Phys. Rev. D17 (1978) 3003.            *
- *                                                                     *
- *  [3] R. Baier, J. Engels and B. Petersson : Z. Phys. C2 (1979) 265. *
- *                                                                     *
- *  [4] Gl"uck, E. Hoffman and E. Reya : Z. Phys. C13 (1982) 119.      *
- *                                                                     *
- *  [5] D.W. Duke and J.F.Owens : Phys. Rev. D30 (1984) 49.            *
- *                                                                     *
- *  [6] E. Eichten, I Hinchliffe, K. Lane and C. Quigg :               *
- *      Rev. Mod. Phys. 56 (1984) 579 and Rev. Mod. Phys. 58 (1985) 1065.
- *                                                                     *
- *  [7] M. Diemoz, F. Ferroni, E. Longo and G. Martinelli :            *
- *      Z. Phys. C39 (1988) 21.                                        *
- *                                                                     *
- *  [8] J.F.Owens : Phys. Lett. 266B (1991) 126.                       *
- *                                                                     *
- *  [9] A.D. Martin, R.G. Roberts and W.J. Stirling : Phys.Rev. D37 (1988) 1161;
- *      Phys. Lett. 206B (1988) 327 and Mod. Phys. Lett. A4 (1989) 1135.
- *                                                                     *
- * [10] P.N. Harriman, A.D. Martin, R.G. Roberts and W.J. Stirling :   *
- *      Phys. Rev. D42 (1990) 798 and Phys. Lett. 243B (1990) 421.     *
- *                                                                     *
- * [11] J. Kwiecinski, A.D. Martin, R.G. Roberts and W.J. Stirling :   *
- *      Phys. Rev. D42 (1990) 3645;                                    *
- *      A.D. Martin, R.G. Roberts and W.J. Stirling :                  *
- *      Phys. Rev. D43 (1991) 3648.                                    *
- *                                                                     *
- * [12] A.D. Martin, R.G. Roberts and W.J. Stirling :                  *
- *      Phys. Rev. D47 (1993) 867.                                     *
- *                                                                     *
- * [13] J. Morfin and W.K. Tung : Z. Phys. C52 (1991) 13.              *
- *                                                                     *
- * [14] M. Gl"uck, E. Reya and A. Vogt : Z. Phys. C48 (1990) 471;      *
- *      Z. Phys. C53 (1992) 127 and Phys. Lett. 306B (1993) 391.       *
- *                                                                     *
- * [15] P. Aurenche et al. : Phys. Rev. D39 (1989) 3275.               *
- *                                                                     *
- * [16] A.D. Martin, R.G. Roberts and W.J. Stirling :                  *
- *      Phys. Lett. 306B (1993) 145 and Phys. Lett. 309B (1993) 492.   *
- *                                                                     *
- * [17] A.D. Martin, R.G. Roberts and W.J. Stirling :                  *
- *      RAL Preprint, RAL-92-078 (1992) and RAL Preprint, RAL-93-077 (199
- *                                                                     *
- * [18] CTEQ collaboration, J. Botts et al. : Phys. Lett. 304B (1993) 159.
- *                                                                     *
- * [19] E.L. Berger and R. Meng : Phys. Lett. 304B (1993) 318;         *
- *      CERN-TH 6739/92; ANL-HEP-CP-92-108 and                         *
- *      E.L. Berger, R. Meng and J. Qiu : ANL-HEP-CP-92-79.            *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- ******************************* - 32 - ********************************
-1***********************************************************************
- *                                                                     *
- * [20] CTEQ collaboration, J. Botts et al. :                          *
- *      to be published                                                *
- *                                                                     *
- *                                                                     *
- * [21] J.F.Owens : Phys. Rev. D30 (1984) 943.                         *
- *                                                                     *
- * [22] P. Aurenche et al. : Phys. Lett. 233B (1989) 517.              *
- *                                                                     *
- * [23] P.J. Sutton, A.D. Martin, R.G. Roberts and W.J. Stirling :     *
- *      Phys. Rev. D45 (1992) 2349.                                    *
- *                                                                     *
- * [24] M. Gl"uck, E. Reya and A. Vogt : Z. Phys. C53 (1992) 651.      *
- *                                                                     *
- * [25] D.W. Duke and J.F.Owens : Phys. Rev. D26 (1982) 1600.          *
- *                                                                     *
- * [26] M. Drees and K. Grassie : Z. Phys. C28 (1985) 451.             *
- *                                                                     *
- * [27] H. Abramowicz, K. Charchula and A. Levy :                      *
- *      Phys. Lett. 269B (1991) 458.                                   *
- *                                                                     *
- * [28] L.E. Gordon and J.K. Storrow : Z. Phys. C56 (1992) 307.        *
- *                                                                     *
- * [29] M. Gl"uck, E. Reya and A. Vogt : Phys. Rev. D46 (1992) 1973    *
- *      and Phys. Rev. D45 (1992) 3986.                                *
- *                                                                     *
- * [30] P. Aurenche et al. : Z. Phys. C56 (1992) 589.                  *
- *                                                                     *
- * [31] G. Altarelli, R.K. Ellis and G. Martinelli :                   *
- *      Nucl. Phys. B143 (1978) 521 and B157 (1979) 461.               *
- *                                                                     *
- * [32] H. Plothow-Besch, 'PDFLIB: a library of all available parton   *
- *      density functions of the nucleon, the pion and the photon      *
- *      and the corresponding alpha(s) calculations',                  *
- *      Comp. Phys. Comm. 75 (1993) 396-416;                           *
- *      H. Plothow-Besch, 'PDFLIB: Structure Functions and alpha(s)    *
- *      Calculations', User's Manual - Version 1.00, W5051 PDFLIB,     *
- *      1991.03.21, CERN-PPE.                                          *
- *                                                                     *
- * [33] Review of Particle Properties, June 1992,                      *
- *      Phys. Rev. D45 (1992), Part. II, p. III.54.                    *
- *                                                                     *
- * [34] BCDMS Collaboration, A.C. Benvenuti et al. : Phys. Lett. 223B (1989) 485
- *      Phys. Lett. 237B (1990) 592 and Phys. Lett. 237B (1990) 599;   *
- *      EMC Collaboration, M. Arneodo et al. : Nucl. Phys. B333 (1990) 1;
- *      NMC Collaboration, D. Allasia et al. : Phys. Lett. 249B (1990) 366;
- *                         P. Amaudruz et al. : Nucl. Phys. B371 (1992) 3.
- *                                                                     *
- * [35] NMC Collaboration, P. Amaudruz et al. : Phys. Lett. 295B (1992) 159.
- *                                                                     *
- * [36] H1 Collaboration, I. Abt et al. : DESY preprint 93-117 (1993); *
- *      ZEUS Collaboration, M. Derrick et al. : DESY preprint 93-110 (1993).
- *                                                                     *
- * [37] JADE Collaboration, W. Bartel et al. : Z. Phys. C24 (1984) 231.*
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- ******************************* - 33 - ********************************
-1***********************************************************************
- *                                                                     *
- * [38] PAW, CERN Program Library Entry Q121, Long Write-up;           *
- *      R. Brun, O. Couet, C. Vandoni and P Zanarini;                  *
- *      see also to References therein.                                *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *  List of Tables                                                     *
- *                                                                     *
- *                                                                     *
- *  1  List of available sets of NUCLEON structure functions in        *
- *     PDFLIB version 6.06 ...................................... 35   *
- *                                                                     *
- *                                                                     *
- *  2  List of available sets of PION structure functions in PDFLIB    *
- *     version 6.06 ............................................. 37   *
- *                                                                     *
- *                                                                     *
- *  3  List of available sets of PHOTON structure functions in PDFLIB  *
- *     version 6.06 ............................................. 38   *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- ******************************* - 34 - ********************************
-1***********************************************************************
- *                                                                     *
- *                             Table 1                                 *
- *                                                                     *
- *       Available sets of NUCLEON structure functions in PDFLIB       *
- *                                                                     *
- * ------------------------------------------------------------------- *
- *                                                                     *
- * Nptype                                                              *
- *  Ngroup      Q**2  Name of Set   Scheme        Reference            *
- *     Nset     min                                                    *
- *         QCDL4                                                       *
- * ------------------------------------------------------------------- *
- *       0  200 0.3 ->GRV HO        MS_bar  NLL   Default              *
- * ------------------------------------------------------------------- *
- * ------------------------------------------------------------------- *
- * 1  1  1            PRIVATE to user                                  *
- * 1  1  2   45 0.5   BEBC                  LO    [Nucl.Phys. B132 (1978) 249]
- * 1  1  3  500 1.8   OR                    LO    [Phys.Rev. D17 (1978) 3003]
- * 1  1  4  470  5    BEP                   LO    [Z. Phys. C2 (1979) 265]
- * 1  1  5  400  4    GHR                   LO    [Z. Phys. C13 (1982) 119]
- * 1  1  6  200  4    DO Set 1              LO    [Phys.Rev. D30 (1984) 49]
- * 1  1  7  400  4    DO Set 2              LO                         *
- * 1  1  8  200  5    EHLQ Set 1            LO    [Rev.Mod.Phys. 56 (1984) 579;
- * 1  1  9  290  5    EHLQ Set 2            LO     Rev.Mod.Phys. 58 (1985) 1065]
- * 1  1 10  177  4    DO Set 1.1            LO    [Phys.Lett.266B (1991) 3275]
- * ------------------------------------------------------------------- *
- * 1  2  1  200 10    DFLM soft valence     LO    [Z. Phys. C39 (1988) 21
- * 1  2  2  200 10    DFLM hard valence     LO                         *
- * 1  2  3  200 10    DFLM soft gluon       LO                         *
- * 1  2  4  200 10    DFLM hard gluon       LO                         *
- * 1  2  5  200 10    DFLM central av.      LO                         *
- * 1  2  6  300 10    DFLM central av. DIS  NLL                        *
- * 1  2  7  160 10    DFLM 160         DIS  NLL   [unpublished]        *
- * 1  2  8  260 10  ->DFLM 260         DIS  NLL                        *
- * 1  2  9  360 10    DFLM 360         DIS  NLL                        *
- * ------------------------------------------------------------------- *
- * 1  3  1  107  5    MRS 1         MS_bar  NLL   [Phys.Rev. D37 (1988) 1181]
- * 1  3  2  250  5    MRS 2         MS_bar  NLL                        *
- * 1  3  3  178  5    MRS 3         MS_bar  NLL                        *
- * 1  3  4   91  5    MRS E         MS_bar  NLL   [Phys.Lett. 206B (1988) 327]
- * 1  3  5  228  5    MRS B         MS_bar  NLL                        *
- * 1  3  6   91  5    MRS E'        MS_bar  NLL   [Mod.Phys. Lett.A4 (1989) 1135
- * 1  3  7  228  5    MRS B'        MS_bar  NLL                        *
- * 1  3  8  100  5    HMRS E (1.90) MS_bar  NLL   [retracted]          *
- * 1  3  9  190  5    HMRS B (1.90) MS_bar  NLL                        *
- * 1  3 10  100  5    HMRS E (3.90) MS_bar  NLL                        *
- * 1  3 11  190  5    HMRS B (3.90) MS_bar  NLL                        *
- * 1  3 12  100  5    HMRS E+(4.90) MS_bar  NLL   [Phys.Lett. 243B (1990) 421]
- * 1  3 13  100  5    HMRS E-(4.90) MS_bar  NLL                        *
- * 1  3 14  100  5    HMRS E (4.90) MS_bar  NLL   [Phys.Rev. D4 (1990) 798]
- * 1  3 15  190  5    HMRS B (4.90) MS_bar  NLL                        *
- * 1  3 16  100  5    HMRS B (8.90) MS_bar  NLL   [unpublished]        *
- * 1  3 17  300  5    HMRS B (8.90) MS_bar  NLL                        *
- * 1  3 18  190  5    KMRS B-(7.90) MS_bar  NLL   [Phys.Rev. D42 (1990) 3645]
- * 1  3 19  190  5    KMRS B-R2 Shadow (7.90) MS_bar  NLL              *
- * 1  3 20  190  5    KMRS B-R5 Shadow (7.90) MS_bar  NLL              *
- * 1  3 21  190  5    KMRS B0-190 (7.90) MS_bar  NLL (as updated 4.92) *
- *                                                                     *
- ******************************* - 35 - ********************************
-1***********************************************************************
- *                                                                     *
- *                          Table 1 (cont'd)                           *
- *                                                                     *
- * ------------------------------------------------------------------- *
- *                                                                     *
- * Nptype                                                              *
- *  Ngroup      Q**2  Name of Set        Scheme      Reference         *
- *     Nset     min                                                    *
- *         QCDL4                                                       *
- * ------------------------------------------------------------------- *
- * 1  3 22  135  5    MRS B0-135 (10.90) MS_bar  NLL [Phys.Rev. D43 (1991) 3648]
- * 1  3 23  160  5    MRS B0-160 (10.90) MS_bar  NLL                   *
- * 1  3 24  200  5    MRS B0-200 (10.90) MS_bar  NLL                   *
- * 1  3 25  235  5    MRS B0-235 (10.90) MS_bar  NLL                   *
- * 1  3 26  215  5    MRS S0     (4.92)  MS_bar  NLL  [retracted]      *
- * 1  3 27  215  5    MRS D0     (4.92)  MS_bar  NLL                   *
- * 1  3 28  215  5    MRS D-     (4.92)  MS_bar  NLL                   *
- * 1  3 29  230  5    MRS S0'    (11.92) MS_bar  NLL  [Phys.Lett. 306B (1993)145
- * 1  3 30  230  5    MRS D0'    (11.92) MS_bar  NLL                   *
- * 1  3 31  230  5    MRS D-'    (11.92) MS_bar  NLL                   *
- * 1  3 32  230  5    MRS S0'    (2.93)  DIS     NLL  [RAL-92-078 (1992)]
- * 1  3 33  230  5    MRS D0'    (2.93)  DIS     NLL                   *
- * 1  3 34  230  5    MRS D-'    (2.93)  DIS     NLL                   *
- * 1  3 35  230  5    MRS (H)    (11.93) MS_bar  NLL                   *
- * 1  3 36  230  5    MRS (H)    (11.93) DIS     NLL                   *
- * 1  3 37  230  5    MRS (A)     (5.94  MS_bar  NLL  [RAL-94-055 (1994)]
- * 1  3 38  230  5    MRS (A) Fit (5.94) MS_bar  NLL                   *
- * 1  3 39  230  5    MRS (Ap)    (2.95  MS_bar  NLL  [RAL-94-104 (1994)]
- * 1  3 40  231  5    MRS (Ap) Fit(2.95) MS_bar  NLL                   *
- * 1  3 41  255  5  ->MRS (G)     (2.95  MS_bar  NLL  [RAL-95-021 (1995)]
- * 1  3 42  255  5    MRS (G) Fit (2.95) MS_bar  NLL                   *
- * ------------------------------------------------------------------- *
- * 1  4  1  212  4    MT S1              DIS     NLL  [Z. Phys. C52 (1991) 13]
- * 1  4  2  194  4    MT B1              DIS     NLL                   *
- * 1  4  3  191  4    MT B2              DIS     NLL                   *
- * 1  4  4  155  4    MT E1              DIS     NLL                   *
- * 1  4  5  237  4    MT 6  (1/2s)       DIS     NLL                   *
- * 1  4  6  212  4    MT S1              MS_bar  NLL                   *
- * 1  4  7  194  4    MT B1              MS_bar  NLL                   *
- * 1  4  8  191  4    MT B2              MS_bar  NLL                   *
- * 1  4  9  155  4    MT E1              MS_bar  NLL                   *
- * 1  4 10  237  4    MT 6  (1/2s)       MS_bar  NLL                   *
- * 1  4 11  144  4    MT LO                      LO                    *
- * 1  4 12  168  4    CTEQ 1L                    LO   [retracted]      *
- * 1  4 13  231  4    CTEQ 1M            MS_bar  NLL                   *
- * 1  4 14  231  4    CTEQ 1MS           MS_bar  NLL                   *
- * 1  4 15  322  4    CTEQ 1ML           MS_bar  NLL                   *
- * 1  4 16  247  4    CTEQ 1D            DIS     NLL                   *
- * 1  4 17  190  4    CTEQ 2L                    LO   [retracted]      *
- * 1  4 18  213  4    CTEQ 2M            MS_bar  NLL                   *
- * 1  4 19  208  4    CTEQ 2MS           MS_bar  NLL                   *
- * 1  4 20  208  4    CTEQ 2MF           MS_bar  NLL                   *
- * 1  4 21  322  4    CTEQ 2ML           MS_bar  NLL                   *
- * 1  4 22  235  4    CTEQ 2D            DIS     NLL                   *
- * 1  4 23  190  4    CTEQ 2pL                   LO   [to be published]*
- * 1  4 24  213  4    CTEQ 2pM           MS_bar  NLL                   *
- * ------------------------------------------------------------------- *
- *                                                                     *
- ******************************* - 36 - ********************************
-1***********************************************************************
- *                                                                     *
- *                          Table 1 (cont'd)                           *
- *                                                                     *
- * ------------------------------------------------------------------- *
- *                                                                     *
- * Nptype                                                              *
- *  Ngroup      Q**2  Name of Set        Scheme      Reference         *
- *     Nset     min                                                    *
- *         QCDL4                                                       *
- * ------------------------------------------------------------------- *
- * 1  4 25  208  4    CTEQ 2pMS          MS_bar  NLL                   *
- * 1  4 26  208  4    CTEQ 2pMF          MS_bar  NLL                   *
- * 1  4 27  322  4    CTEQ 2pML          MS_bar  NLL                   *
- * 1  4 28  235  4    CTEQ 2pD           DIS     NLL                   *
- * 1  4 29  177  4    CTEQ 3L                    LO   [MSU-HEP/41024]  *
- * 1  4 30  239  4  ->CTEQ 3M            MS_bar  NLL                   *
- * 1  4 31  247  4    CTEQ 3D            DIS     NLL                   *
- * ------------------------------------------------------------------- *
- * 1  5  1  160 0.2   old GRV HO         MS_bar  NLL  [retracted]      *
- * 1  5  2  220 0.2   old GRV LO                 LO                    *
- * 1  5  3  200 0.3   GRV HO             MS_bar  NLL  [Z. Phys. C53 (1992) 127]
- * 1  5  4  200 0.25  GRV LO                     LO                    *
- * 1  5  5  200 0.4   GRV 94 LO                  LO   [DO-TH 94/24 (1994)]
- * 1  5  6  200 0.4 ->GRV 94 HO          MS_bar  NLL                   *
- * 1  5  7  200 0.4   GRV 94 HO          DIS     NLL                   *
- * ------------------------------------------------------------------- *
- * 1  6  1  230  2    ABFOW              MS_bar  NLL  [Phys.Rev D39 (1989) 3275]
- * ------------------------------------------------------------------- *
- * 1  7  1  254  5    BM A               MS_bar  NLL  [Phys.Lett. 304B (1993)318
- * 1  7  2  254  5    BM B               MS_bar  NLL                   *
- * ------------------------------------------------------------------- *
- *                                                                     *
- *                                                                     *
- *                             Table 2                                 *
- *                                                                     *
- *         Available sets of PION structure functions in PDFLIB        *
- *                                                                     *
- * ------------------------------------------------------------------- *
- *                                                                     *
- * Nptype                                                              *
- *  Ngroup      Q**2  Name of Set   Scheme        Reference            *
- *     Nset     min                                                    *
- *         QCDL4                                                       *
- * ------------------------------------------------------------------- *
- * 2  1  1  200  4    OW-P Set 1            LO    [Phys.Rev. D30 (1984) 943]
- * 2  1  2  400  4    OW-P Set 2            LO                         *
- * ------------------------------------------------------------------- *
- * 2  3  1  190  5    SMRS-P 1      MS_bar  NLL   [Phys.Rev. D45 (1992) 2349]
- * 2  3  2  190  5  ->SMRS-P 2      MS_bar  NLL                        *
- * 2  3  3  190  5    SMRS-P 3      MS_bar  NLL                        *
- * ------------------------------------------------------------------- *
- * 2  5  1  200 0.3 ->GRV-P HO      MS_bar  NLL   [Z.Phys. C53 (1992) 651]
- * 2  5  2  200 0.25  GRV-P LO              LO                         *
- * ------------------------------------------------------------------- *
- * 2  6  1  231  2  ->ABFKW-P 1     MS_bar  NLL   [Phys.Lett.233B (189) 517]
- * 2  6  2  181  2    ABFKW-P 2     MS_bar  NLL                        *
- * 2  6  3  281  2    ABFKW-P 3     MS_bar  NLL                        *
- * ------------------------------------------------------------------- *
- ******************************* - 37 - ********************************
-1***********************************************************************
- *                                                                     *
- *                             Table 3                                 *
- *                                                                     *
- *        Available sets of PHOTON structure functions in PDFLIB       *
- *                                                                     *
- * ------------------------------------------------------------------- *
- *                                                                     *
- * Nptype                                                              *
- *  Ngroup      Q**2  Name of Set   Scheme        Reference            *
- *     Nset     min                                                    *
- *         QCDL4                                                       *
- * ------------------------------------------------------------------- *
- * 3  1  1  380 10    DO-G LO              LO     [Phys.Rev. D26 (1982) 1600]
- * 3  1  2  440 10  ->DO-G HO       MS_bar NLL                         *
- * ------------------------------------------------------------------- *
- * 3  2  1  400  1  ->DG-G Set 1           LO     [Z. Phys. C28 (1985) 451]
- * 3  2  2  400  1    DG-G Set 2           LO                          *
- * 3  2  3  400 20    DG-G Set 3           LO                          *
- * 3  2  4  400 200   DG-G Set 4           LO                          *
- * ------------------------------------------------------------------- *
- * 3  3  1  200  5    LAC-G 1              LO     [Phys.Lett.269B (1991) 458]
- * 3  3  2  200  5  ->LAC-G 2              LO                          *
- * 3  3  3  200  5    LAC-G 3              LO                          *
- * ------------------------------------------------------------------- *
- * 3  4  1  200 5.3 ->GS-G HO       MS_bar NLL    [Z.Phys. C56 (1992) 307]
- * 3  4  2  200 5.3   GS-G LO set 1        LO                          *
- * 3  4  3  200 5.3   GS-G LO set 2        LO                          *
- * ------------------------------------------------------------------- *
- * 3  5  1  200 0.3   GRV-G LHO     DIS*)  NLL    [Phys.Rev. D46 (1992) 1973;
- * 3  5  2  200 0.3 ->GRV-G HO      DIS*)  NLL     Phys.Rev. D45 (1992) 3986]
- * 3  5  3  200 0.25  GRV-G LO             LO                          *
- * ------------------------------------------------------------------- *
- * 3  6  1  200  2    ACFGP-G HO    MS_bar NLL    [Z.Phys. C56 (1992) 589]
- * 3  6  2  200  2  ->ACFGP-G HO-mc MS_bar NLL                         *
- * 3  6  3  200  2  ->AFG-G HO      MS_bar NLL                         *
- * ------------------------------------------------------------------- *
- * 3  8  1  400  4    WHIT-G 1             LO     [KEK preprint 93-160]*
- * 3  8  2  400  4  ->WHIT-G 2             LO                          *
- * 3  8  3  400  4    WHIT-G 3             LO                          *
- * 3  8  4  400  4    WHIT-G 4             LO                          *
- * 3  8  5  400  4    WHIT-G 5             LO                          *
- * 3  8  6  400  4    WHIT-G 6             LO                          *
- * ------------------------------------------------------------------- *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *  *) not standard, please consult references.                        *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- ******************************* - 38 - ********************************
-1***********************************************************************
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- ******************************* - 39 - ********************************
diff --git a/PDF/doc/v100.dat b/PDF/doc/v100.dat
deleted file mode 100644 (file)
index 9b457ff..0000000
+++ /dev/null
@@ -1,12 +0,0 @@
-*
-* $Id$
-*
-* $Log$
-* Revision 1.1.1.1  1996/03/08 17:48:06  mclareni
-* Pdf
-*
-*
-* VERSION 6.06   950315 17.05 CERN PROGRAM LIBRARY     W5051
-C
-C- Structuring of the code in PATCHY format
-C- Addition of the test
diff --git a/PDF/dpdf/inform.dat b/PDF/dpdf/inform.dat
deleted file mode 100644 (file)
index dccfe39..0000000
+++ /dev/null
@@ -1,31 +0,0 @@
-*
-* $Id$
-*
-* $Log$
-* Revision 1.1.1.6  1997/07/13 13:22:18  cernlib
-* Import version 7.09
-*
-* Revision 1.7  1997/07/02 15:02:58  pdflib
-* potential bug removed in GRSGALO
-*
-* Revision 1.6  1997/04/08 15:53:15  pdflib
-* Lambda_QCD for CTEQ4 sets corrected
-*
-* Revision 1.5  1996/12/10 10:49:21  pdflib
-* Version 7.07: consistency change for u,d-valence in CTEQ4
-*
-* Revision 1.4  1996/11/15 15:06:49  pdflib
-* Version 7.06: tabs removed and bug fixed in grsgalo
-*
-* Revision 1.3  1996/11/07 09:34:23  pdflib
-* New version 7.05: cut routine Ctq4RTb in 9 pieces
-*
-* Revision 1.1.1.1  1996/04/12 15:28:51  plothow
-* Version 7.01
-*
-*
-*  GENERAL INFORMATION
-C PDF      7.09      970702 16.05 CERN Program Library     W5051 PDF
-C          Structure Functions Package
-C          Author: H. Plothow-Besch
-C
diff --git a/PDF/dpdf/pdfdoc.doc b/PDF/dpdf/pdfdoc.doc
deleted file mode 100644 (file)
index 53a2eca..0000000
+++ /dev/null
@@ -1,2829 +0,0 @@
-1***********************************************************************
- *                                                                     *
- *   CERN                                                              *
- *                                                                     *
- *   W5051 PDFLIB                                         1997.07.02   *
- *                                                                     *
- *   ---------------------------------------------------------------   *
- *                                                                     *
- *      PDFLIB: Nucleon, Pion and Photon Parton Density Functions      *
- *                      and alpha(s) Calculations                      *
- *                                                                     *
- *                           User's Manual                             *
- *                                                                     *
- *                           Version 7.09                              *
- *                                                                     *
- *                                                                     *
- *                     H. Plothow-Besch/CERN-PPE                       *
- *                                                                     *
- *                         CH-1211 Geneva 23                           *
- *               E-mail address: plothow@mail.cern.ch                  *
- *                    Tel.: +41 - 22 - 767 39 80                       *
- *                                                                     *
- *   A complete manual in PostScript format exists on a separate file  *
- *     Requests can be addressed to the CERN Program Library Office    *
- *                                                                     *
- *                   Copyright by H. Plothow-Besch                     *
- *                                                                     *
- *                            Abstract                                 *
- *                                                                     *
- *       This document describes an integrated package of Parton Den-  *
- *  sity Functions called PDFLIB which has been added to the CERN      *
- *  Program Library Pool W999 and is labelled as W5051. In this pack-  *
- *  age all the different sets of parton density functions of the      *
- *  NUCLEON, the PION and the PHOTON which are available today (about  *
- *  100 different sets) have been put together. All these sets have    *
- *  been combined in a consistent way such that they all have similar  *
- *  calling sequences and no external data files have to be read in    *
- *  anymore. A default set has been prepared, although those prefering *
- *  their own set or wanting to test a new one may do so within the    *
- *  package. The selection mode, the program parameters as well as     *
- *  the possibilities to modify the defaults and to control errors     *
- *  occured during execution are described.                            *
- *                                                                     *
- *       The package offers as well a program to calculate the strong  *
- *  coupling constant alpha(s) to first or second order (default).     *
- *  The correct Lambda(QCD) associated to the selected set of struc-   *
- *  ture functions and the number of allowed flavours with respect to  *
- *  the given Q**2 is automatically used in the caculation. For those  *
- *  structure functions where the evolution has been performed to      *
- *  leading order only, an alpha(s) value fixed to first order is re-  *
- *  turned.                                                            *
- *                                                                     *
- *       The package is written in standard FORTRAN-77 and is avai-    *
- *  lable on all CERN system areas.                                    *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *        Submitter:   H. Plothow-Besch/CERN-PPE                       *
- *        Language:    FORTRAN 77                                      *
- *        Library:     PDFLIB                                          *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- *                                                                     *
- ***********************************************************************
-1***********************************************************************
- *                                                                     *
- *                                                                     *
- *   _______________________________________________________________   *
- *   |                                                             |   *
- *   |  PDFLIB - User's Manual                                     |   *
- *   |                                                             |   *
- *   |  CERN Computer Program Library entry W5051                  |   *
- *   |                                                             |   *
- *   |  Copyright CERN, Geneva 1990, 1991, 1992, 1993, 1994, 1995  |   *
- *   |  Copyright CERN, Geneva 1996                                |   *
- *   |                                                             |   *
- *   |  Copyright and any other appropriate legal protection of    |   *
- *   |  these computer programs and associated documentation re-   |   *
- *   |  served in all countries of the world.                      |   *
- *   |                                                             |   *
- *   |  These programs or documentation may not be reproduced      |   *
- *   |  and/or redistributed by any method without prior written   |   *
- *   |  consent of the Director-General of CERN or his delegate.   |   *
- *   |                                                             |   *
- *   |  Permission for the scientific usage of any programs des-   |   *
- *   |  cribed herein is granted apriori to those scientific in-   |   *
- *   |  stitutes associated with the CERN experimental program or  |   *
- *   |  with whom CERN has concluded a scientific collaboration    |   *
- *   |  agreement.                                                 |   *
- *   |                                                             |   *
- *   |  Commercial utilisation requires explicit a priori per-     |   *
- *   |  mission from CERN and will be subjected to payment of a    |   *
- *   |  licence fee.                                               |   *
- *   |