PDF/spdf/pdfsta.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.2  1996/10/30 08:28:58  cernlib
   6 * Version 7.04
   7 *
   8 * Revision 1.1.1.1  1996/04/12 15:29:25  plothow
   9 * Version 7.01
  10 *
  11 *
  12 #include "pdf/pilot.h"
  13       SUBROUTINE PDFSTA
  14 C
  15 C **********************************************************************
  16 C *                                                                    *
  17 C *   Routine to print status of parameters in the COMMON/W50511/,     *
  18 C *                                                COMMON/W50512/,     *
  19 C *                                                COMMON/W50513/      *
  20 C *                                                COMMON/W50514/      *
  21 C *                                                COMMON/W505110/     *
  22 C *                                            and COMMON/W505122/     *
  23 C *                                                                    *
  24 C *   Author:  H. Plothow-Besch                                        *
  25 C *            CERN-PPE, CH - 1211 Geneva 23, Switzerland              *
  26 C *                                                                    *
  27 C **********************************************************************
  28 C
  29 #include "pdf/impdp.inc"
  30 C
  31 #include "pdf/w5051p2.inc"
  32 #include "pdf/w5051p7.inc"
  33 #include "pdf/w50510.inc"
  34 #include "pdf/w50511.inc"
  35 #include "pdf/w50512.inc"
  36 #include "pdf/w50513.inc"
  37 #include "pdf/w50514.inc"
  38 #include "pdf/w50515.inc"
  39 #include "pdf/w50517.inc"
  40 #include "pdf/w505110.inc"
  41 #include "pdf/w505122.inc"
  42       DATA ZEROD/0.D0/
  43 C
  44       IFLSTA = 1
  45 C.
  46 C print values
  47       IF(N6.LE.0) N6 = L6
  48       WRITE(N6,2000)
  49       WRITE(N6,2001) NPTYPE,NGROUP,NSET,SFNAME(NPTYPE,NGROUP,NSET),
  50      +               MODECR(NPTYPE,NGROUP,NSET)
  51       WRITE(N6,2002) NFL,LO,TMAS
  52       WRITE(N6,2003) QCDL4,QCDL5
  53       WRITE(N6,2004) XMIN,XMAX,Q2MIN,Q2MAX
  54  2000 FORMAT(/,'  PDFLIB : Summary from PDFSTA')
  55  2001 FORMAT(1H ,' Nptype = ',I1,'  Ngroup = ',I1,'  Nset = ',I2,
  56      +         '  Name  = "',A8,'"  CrMode = ',I3)
  57  2002 FORMAT(1H ,' Nfl    = ',I2,',  LO = ',I1,
  58      +           ',  Tmas = ',F7.2,' GeV/c**2')
  59  2003 FORMAT(1H ,' QCDL4  = ',F7.4,' GeV,  QCDL5 = ',F7.4,' GeV')
  60  2004 FORMAT(1H ,' Xmin   = ',1X,E8.2,',  Xmax = ',E11.5,
  61      +          ',  Q2min = ',F7.3,' (GeV/c)**2,  Q2max = ',E8.2,
  62      +          ' (GeV/c)**2')
  63 C.
  64       WRITE(N6,3000)
  65       IF(WXMIN .GT. ZEROD) THEN
  66          PERC = 100.D0*WXMIN/WTXMIN
  67          WRITE(N6,3001) WXMIN,WTXMIN,PERC,XMIN
  68       ENDIF
  69       IF(WXMAX .GT. ZEROD) THEN
  70          PERC = 100.D0*WXMAX/WTXMAX
  71          WRITE(N6,3002) WXMAX,WTXMAX,PERC,XMAX
  72       ENDIF
  73       IF(WQ2MIN .GT. ZEROD) THEN
  74          PERC = 100.D0*WQ2MIN/WTQ2MIN
  75          WRITE(N6,3003) WQ2MIN,WTQ2MIN,PERC,Q2MIN
  76       ENDIF
  77       IF(WQ2MAX .GT. ZEROD) THEN
  78          PERC = 100.D0*WQ2MAX/WTQ2MAX
  79          WRITE(N6,3004) WQ2MAX,WTQ2MAX,PERC,Q2MAX
  80       ENDIF
  81       IF((WXMIN+WXMAX+WQ2MIN+WQ2MAX) .EQ. ZEROD) WRITE(N6,3005)
  82  3000 FORMAT(' ')
  83  3001 FORMAT(1H ,' PDFSTA:    X values ',E9.3,' times out of ',E9.3,
  84      +         ' ( = ',E9.3,' %) smaller allowed minimum of ',E9.3)
  85  3002 FORMAT(1H ,' PDFSTA:    X values ',E9.3,' times out of ',E9.3,
  86      +         ' ( = ',E9.3,' %) bigger  allowed maximum of ',E11.5)
  87  3003 FORMAT(1H ,' PDFSTA: Q**2 values ',E9.3,' times out of ',E9.3,
  88      +         ' ( = ',E9.3,' %) smaller allowed minimum of ',E9.3,
  89      +         ' (GeV/c)**2')
  90  3004 FORMAT(1H ,' PDFSTA: Q**2 values ',E9.3,' times out of ',E9.3,
  91      +         ' ( = ',E9.3,' %) bigger  allowed maximum of ',E9.3,
  92      +         ' (GeV/c)**2')
  93  3005 FORMAT(1H ,' PDFSTA:    NO errors occured')
  94 C.
  95       RETURN
  96       END