Adding some QCD diffractive states to the PDG list

[u/mrichter/AliRoot.git] / PDF / tpdf / testpdf.F

diff --git a/PDF/tpdf/testpdf.F b/PDF/tpdf/testpdf.F
index 8a10d801cb27484303d3c4961812c9c8dd7db3d9..97b2a931828324897f12f35c5c62b29f612b517d 100644 (file)
--- a/PDF/tpdf/testpdf.F
+++ b/PDF/tpdf/testpdf.F
@@ -2,8 +2,11 @@
 * $Id$
 *
 * $Log$
 * $Id$
 *
 * $Log$

-* Revision 1.1.1.4  1997/07/13 13:25:32  cernlib
-* Import version 7.09
+* Revision 1.1.1.5  2000/04/19 13:30:02  mclareni
+* Import version 8.04
+*
+* Revision 1.5  2000/04/09 10:06:56  plothow
+* new version 8.00 - MRS 98 and 99, CTEQ5 and GRV98 added

 *
 * Revision 1.4  1997/07/02 15:03:12  pdflib
 * potential bug removed in GRSGALO
 *
 * Revision 1.4  1997/07/02 15:03:12  pdflib
 * potential bug removed in GRSGALO


@@ -38,13 +41,13 @@ C
      +       DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL,
      +       DQ2,DP2,ALF,ALPHAS2
       REAL    X, Q, UPV, DNV, USEA, DSEA, SSEA, CSEA, BOT, TOP, GLUE
      +       DX,DQ,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL,
      +       DQ2,DP2,ALF,ALPHAS2
       REAL    X, Q, UPV, DNV, USEA, DSEA, SSEA, CSEA, BOT, TOP, GLUE

-      DIMENSION XX(3),Q2(3)
+      DIMENSION XX(5),Q2(5)

       CHARACTER*20 PARM(NCHDIM)
 #include "pdf/expdp.inc"
      +       VAL(NCHDIM)
 C
       CHARACTER*20 PARM(NCHDIM)
 #include "pdf/expdp.inc"
      +       VAL(NCHDIM)
 C

-      DATA NXX,XX/3, .0005,.140,.700/
-      DATA NQ2,Q2/3, 1.30,24.0,33.5/
+      DATA NXX,XX/5, .0005,.140,.700,0.01,0.05/
+      DATA NQ2,Q2/5, 1.30,24.0,33.5,200.0,200.0/

       DATA P2/5./
       DATA IP2/0/
 C
       DATA P2/5./
       DATA IP2/0/
 C


@@ -68,7 +71,7 @@ C Test new PDFLIB versions (4.0 or bigger)
             DO ISET=-9,NPGSMX(KPTY,KGRO)
                IF((KPTY*KGRO).NE.1 .AND. ISET.LE.0) GOTO 300
                WRITE (N6,*) ' ==========================================
             DO ISET=-9,NPGSMX(KPTY,KGRO)
                IF((KPTY*KGRO).NE.1 .AND. ISET.LE.0) GOTO 300
                WRITE (N6,*) ' ==========================================

-     +================================================================='
+     +============================================'

                FIRST=.TRUE.
                IF(ISET.LE.-8) IFLSET = 0
                IF(ISET.LT.0) GOTO 50
                FIRST=.TRUE.
                IF(ISET.LE.-8) IFLSET = 0
                IF(ISET.LT.0) GOTO 50


@@ -149,8 +152,8 @@ C
                WRITE(N6,4000) NPTYPE,NGROUP,NSET,
      +         SFNAME(NPTYPE,NGROUP,NSET),MODECR(NPTYPE,NGROUP,NSET)
                WRITE(N6,4001) QCDL4,QCDL5,DQ,ALF
                WRITE(N6,4000) NPTYPE,NGROUP,NSET,
      +         SFNAME(NPTYPE,NGROUP,NSET),MODECR(NPTYPE,NGROUP,NSET)
                WRITE(N6,4001) QCDL4,QCDL5,DQ,ALF

- 4000          FORMAT(/,'  Nptype = ',I1,'  Ngroup = ',I1,
-     +                  '  Nset = ',I2,'  Name  = "',A8,
+ 4000          FORMAT(/,'  Nptype = ',I1,'  Ngroup = ',I2,
+     +                  '  Nset = ',I3,'  Name  = "',A8,

      +                  '"  CrMode = ',I3)
  4001          FORMAT(1H ,' QCDL4, QCDL5 = ',2F7.4,
      +                    '    Scale = ',F8.4,'  alpha(s) = ',F7.5)
      +                  '"  CrMode = ',I3)
  4001          FORMAT(1H ,' QCDL4, QCDL5 = ',2F7.4,
      +                    '    Scale = ',F8.4,'  alpha(s) = ',F7.5)