EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
leptod.F
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file leptod.F
1 
2 C **********************************************************************
3 
4  BLOCK DATA leptod
5 
6  IMPLICIT NONE
7 C...Give sensible default values to switches and parameters.
8 
9 
10 *
11 * to avoid variable conflictions, a second keep element is necessary
12 * with the same common block name (see LPTOU2)
13 *
14  COMMON /leptou/ cut(14),lst(40),parl(30),
15  & x,y,w2,q2,u
16  REAL cut,parl,x,y,w2,q2,u
17  INTEGER lst
18  SAVE /leptou/
19 
20  COMMON /linter/ pari(50),ewqc(2,2,8),qc(8),zl(2,4),zq(2,8),pq(17)
21  REAL pari,ewqc,qc,zl,zq,pq
22  SAVE /linter/
23 
24  COMMON /lflmix/ cabibo(4,4)
25  REAL cabibo
26 *DEC Station gave a warning SAVE /CABIBO/
27 
28  COMMON /loptim/ optx(4),opty(4),optq2(4),optw2(4),comfac
29  REAL optx,opty,optq2,optw2,comfac
30  SAVE /loptim/
31 
32  COMMON /lgrid/ nxx,nww,xx(31),ww(21),pqg(31,21,3),pqqb(31,21,2),
33  &qgmax(31,21,3),qqbmax(31,21,2),ycut(31,21),xtot(31,21),np
34  REAL xx,ww,pqg,pqqb,qgmax,qqbmax,ycut,xtot
35  INTEGER nxx,nww,np
36  SAVE /lgrid/
37 
38 
39  COMMON /flgrid/ nfx,nfq,xr(2),qr(2),flqt(41,16),flgt(41,16),
40  &flmt(41,16)
41  REAL xr,qr,flqt,flgt,flmt
42  INTEGER nfx,nfq
43  SAVE /flgrid/
44 
45  COMMON /pypara/ ipy(80),pypar(80),pyvar(80)
46  REAL pypar,pyvar
47  INTEGER ipy
48  SAVE /pypara/
49 
50  COMMON /lminui/ xkin(4),ukin(4),wkin(4),ain(4),bin(4),
51  &maxfin,relup,relerr,reler2,fcnmax
52  REAL xkin,ukin,wkin,ain,bin,relerr,relup,reler2,fcnmax
53  INTEGER maxfin
54  SAVE /lminui/
55 
56  COMMON /lminuc/ namkin(4),nam(30)
57  CHARACTER*10 namkin,nam
58  SAVE /lminuc/
59 
60  INTEGER imxpdf
61  parameter(imxpdf=40)
62  COMMON /pepadm/cpdfnam(2,imxpdf),ipdfnam(2,imxpdf),
63  & iplst(10),cunpol,cpol
64  CHARACTER*256 cpdfnam,cunpol,cpol
65  INTEGER iplst,ipdfnam
66  SAVE /pepadm/
67 
68 **************************************************************
69 *
70 * IPLST(1) = 0 (default) : number of PDF warnings
71 * IPLST(2) = 11 (default) : unit -1- for pdf files
72 * IPLST(3) = 12 (default) : unit -2- for pdf files
73 *
74 **************************************************************
75 
76 
77  INTEGER i
78 
79 C...LEPTOU: Cuts, basic switches and parameters.
80  DATA cut/1.e-04,1.,0.,1.,4.,1.e+08,5.,1.e+08,1.,1.e+08,1.,1.e+08,
81  &0.,3.1416/
82 C... 0 1 2 3 4 5 6 7 8 9
83  DATA lst/ 0, 1, 5, 1, 3, 1, 1, 12, 5,
84  1 1, 0, 4, 5, 4, 9, 1, 0, 2, -10,
85  2 5, 0, 0, 0, 0, 0, 0, 0, 0, 0,
86  3 0, 0, 0, 0, 1, 1, 0, 0, 1, 0,
87  4 0 /
88 C... 0 1 2 3 4 5 6 7 8 9
89  DATA parl/ 1., 1., 0.44,0.75,.2319,0., 0.5, .04, 4.,
90  1 0.0, 0.01, 0.01,0.1, 0.35,0.01,7.29735e-03,
91  & 1.16639e-05,0.044,0.03,
92  2 0.1,10*0./
93 C...Internally used variables.
94  DATA pari/50*0./
95  DATA qc/-.33333,.66667,-.33333,.66667,-.33333,.66667,
96  & -.33333,.66667/
97  DATA cabibo/.95,.05,2*0.,.05,.948,.002,2*0.,.002,.998,4*0.,1./
98  DATA optx/1.,3*0./,opty/1.,3*0./,optq2/1.,3*0./,optw2/1.,3*0./
99  DATA nxx,nww/31,21/
100  DATA pqg,pqqb,qgmax,qqbmax/6510*0./,ycut/651*0./,xtot/651*1.0e-12/
101  DATA nfx,nfq/41,16/,flqt,flgt,flmt/1968*0./
102  DATA xkin/1.,2.,3.,4./,ukin,wkin,ain,bin/16*0./,maxfin/2000/
103  DATA relup,relerr,reler2/0.1,0.05,0.05/
104  DATA namkin/' x',' ',' ',' '/
105  DATA ipy/
106  1 0, 0, 2, 2, 6, 1, 1, 6, 3, 1,
107  2 3, 1, 1, 2, 1, 1, 4, 1, 1, 1,
108  3 0, 1, 1, 1, 1, 1, 1, 0, 0, 0,
109  4 1, 2, 1, 1, 30, 33, 1, 1, 7, 0,
110  5 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
111  6 0, 0, 0, 1, 100, 0, 0, 0, 0, 0,
112  7 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
113  8 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
114  DATA (pypar(i),i=1,40)/
115  1 7.299e-03, 2.290e-01, 2.000e-01, 2.500e-01, 4.000e+00,
116  1 1.000e+00, 4.400e-01, 4.400e-01, 7.500e-02, 0.000e+00,
117  2 2.000e+00, 0.000e+00, 1.000e+00, 0.000e+00, 3.000e+00,
118  2 1.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 1.000e+00,
119  3 2.500e-01, 1.000e+00, 2.000e+00, 1.000e-03, 1.000e+00,
120  3 1.000e+00, 1.000e+00, -2.000e-02, -1.000e-02, 0.000e+00,
121  4 0.000e+00, 1.600e+00, 0.500e+00, 0.200e+00, 3.894e-01,
122  4 1.000e+00, 3.300e-01, 6.600e-01, 0.000e+00, 1.000e+00/
123  DATA (pypar(i),i=41,80)/
124  5 2.260e+00, 1.000e+04, 1.000e-04, 0.000e+00, 0.000e+00,
125  5 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00,
126  6 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00,
127  6 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00,
128  7 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00,
129  7 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00,
130  8 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00,
131  8 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00/
132  DATA pyvar/80*0./
133 
134  DATA iplst/0,11,12,7*0/
135  DATA cpdfnam /
136  1 './pdf/MRSa1.dat','./pdf/polnloA.dat',
137  2 './pdf/MRSa1.dat','./pdf/polnloB.dat',
138  3 './pdf/MRSa1.dat','./pdf/polnloC.dat',
139  4 ' ','./pdf/polparA.dat',' ','./pdf/polparB.dat',
140  6 ' ','./pdf/polparC.dat',
141  7 ' ','./pdf/STDNLO.GRID',' ','./pdf/VALNLO.GRID',
142  9 ' ','./pdf/STDLO.GRID',' ','./pdf/VALLO.GRID',
143  1 ' ','./pdf/S1NLO.GRID',
144  2 ' ','./pdf/S2NLO.GRID',
145  3 ' ','./pdf/S3NLO.GRID',
146  4 ' ','./pdf/S1LO.GRID',
147  5 ' ','./pdf/S2LO.GRID',
148  6 ' ','./pdf/S3LO.GRID',
149  7 './pdf/cteq4lq.tbl',' ',
150  8 './pdf/MRSasq.dat','./pdf/MRSalq.dat',
151 ****** grsv with grv98 a unpol refenrence |bs>*********************
152  9 './pdf/grv98lo.grid','./pdf/std2000_lo_g1.grid',
153  1 './pdf/grv98lo.grid','./pdf/val2000_lo_g1.grid',
154  2 './pdf/grv98nlo.grid','./pdf/std2000_nlo_g1.grid',
155  3 './pdf/grv98nlo.grid','./pdf/std2000_nlo_g1.grid',
156 ****** CTEQ 5 unpol ***********************************************
157  4 './pdf/cteq5l.tbl',' ',
158  5 './pdf/cteq5m.tbl',' ',
159  6 './pdf/cteq5m1.tbl',' ',
160 ****** MRS 99 unpol ***********************************************
161  7 './pdf/cor01.dat',' ',
162  8 './pdf/cor02.dat',' ',
163  9 './pdf/cor03.dat',' ',
164  1 './pdf/cor04.dat',' ',
165  2 './pdf/cor05.dat',' ',
166  3 './pdf/cor06.dat',' ',
167  4 './pdf/cor07.dat',' ',
168  5 './pdf/cor08.dat',' ',
169  6 './pdf/cor09.dat',' ',
170  7 './pdf/cor10.dat',' ',
171  8 './pdf/cor11.dat',' ',
172  9 './pdf/cor12.dat',' ',
173  1 './pdf/cteq6l.tbl',' ',
174  2 './pdf/cteq6d.tbl',' ',
175  3 './pdf/cteq6m.tbl',' ' /
176  DATA ipdfnam /
177  1 107,1,108,2,109,3,
178  4 113,4,114,5,115,6,
179  7 116,7,117,8,118,9,119,10,
180  1 124,11,125,12,126,13,127,14,128,15,129,16,
181  7 133,17,137,18,144,19,145,20,146,21,
182  1 147,22,150,23,151,24,152,25,161,26,
183  2 162,27,163,28,164,29,165,30,166,31,167,32,
184  3 168,33,169,34,170,35,171,36,172,37,173,38,174,39,175,40 /
185 
186  END