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