Corrected path to mapping libraries.
[u/mrichter/AliRoot.git] / PDF / tpdf / tpdfupv.F
CommitLineData
21886bb6 1*
2* $Id$
3*
4* $Log$
5* Revision 1.1.1.2 1996/10/30 08:32:57 cernlib
6* Version 7.04
7*
8* Revision 1.1.1.1 1996/04/12 15:30:18 plothow
9* Version 7.01
10*
11*
12#include "pdf/pilot.h"
13 PROGRAM PDFUPV
14C
15 PARAMETER (LFLPRT = 2)
16#include "pdf/w5051p2.inc"
17#include "pdf/w5051p7.inc"
18#include "pdf/w50510.inc"
19#include "pdf/w50511.inc"
20#include "pdf/w505110.inc"
21#include "pdf/w505120.inc"
22#include "pdf/w505122.inc"
23C define HBOOK settings
24 PARAMETER (NHBMEM = 500000)
25 COMMON /PAWC/ HMEM(NHBMEM)
26 PARAMETER (NB=900, ID= 0)
27C define DOUBLE PRECISION variables for calling sequences
28#include "pdf/expdp.inc"
29 + DX,DSCALE,DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL,
30 + ALF,ALPHAS2
31 REAL X, SCALE, UPV, DNV, USEA, DSEA, STR, CHM, BOT, TOP, GL
32#include "pdf/w50516.inc"
33 CHARACTER*20 PARM(NCHDIM)
34#include "pdf/expdp.inc"
35 + VAL(NCHDIM)
36C
37 DATA SCALE/80.140/
38 DATA X00/0.002/, DX0/0.001/
39C
40 N6 = L6
41 IFLPRT = LFLPRT
42C
43 CALL HLIMIT(NHBMEM)
44C
45 DSCALE=SCALE
46C
47C First call to PDFSET to initialize COMMON/W505120/
48 PARM(1) = 'Init0'
49 VAL(1) = 0.D0
50 CALL PDFSET(PARM,VAL)
51C loop over all existing sets of Nucleon structure functions (SF)
52 NPTYPE = 1
53 NSF = 0
54 DO IGR = 1,NGRMAX
55 IF(NPGSMX(NPTYPE,IGR).EQ.0) GOTO 30
56 DO ISET = 1,NPGSMX(NPTYPE,IGR)
57 WRITE (N6,*) ' =============================================
58 +================================================================='
59 NSF = NSF + 1
60C book histograms for each set of SF separately
61 CALL HBOOK1(ID+NSF,'U Valence$',NB,0.0015,0.9015,0.)
62C force label printing for each set of SF (not only 1st)
63 FIRST = .TRUE.
64C select and set parameters
65 PARM(1) = 'Nptype'
66 VAL(1) = NPTYPE
67 PARM(2) = 'Ngroup'
68 VAL(2) = IGR
69 PARM(3) = 'Nset'
70 VAL(3) = ISET
71 CALL PDFSET(PARM,VAL)
72C loop over all x bins
73 DO 10 I=1,NB
74 X= X00 + (I-1)*DX0
75 DX=X
76 IF(X.LT.0.0015 .OR. X.GT.0.9015) GOTO 10
77 CALL STRUCTM(DX,DSCALE,
78 + DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL)
79 UPV=DUPV
80 IF(X.GT.0.499 .AND. X.LE.0.500) WRITE(N6,1000) X,SCALE,UPV
81 1000 FORMAT(/,' X= ',F6.4,' Q= ',F6.3,' UPV= ',F8.4)
82 CALL HF1(ID+NSF,X,UPV)
83 10 CONTINUE
84C get alpha(s) for selected set of SF at Q = SCALE
85 ALF = ALPHAS2(DSCALE)
86 WRITE(N6,4000) NPTYPE,NGROUP,NSET,SFNAME(NPTYPE,NGROUP,NSET),
87 + MODECR(NPTYPE,NGROUP,NSET),NSF
88 WRITE(N6,4001) DSCALE,ALF
89 4000 FORMAT(/,' Nptype = ',I1,' Ngroup = ',I1,
90 + ' Nset = ',I2,' Name = "',A8,
91 + '" CrMode = ',I3,' HBId = ',I3)
92 4001 FORMAT(1H ,' Scale = ',F8.4,' alpha(s) = ',F6.4)
93C get error summary for each set of SF
94 CALL PDFSTA
95 20 CONTINUE
96 ENDDO
97 30 CONTINUE
98 ENDDO
99C
100 STOP
101 END