]> git.uio.no Git - u/mrichter/AliRoot.git/blame - ALIROOT/pydata.F
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
These functions now in minicernlib
[u/mrichter/AliRoot.git] / ALIROOT / pydata.F
CommitLineData
fe4da5cc 1
2C*********************************************************************
3
4 BLOCK DATA PYDATA
5
6C...Give sensible default values to all status codes and parameters.
7 COMMON/PYSUBS/MSEL,MSUB(200),KFIN(2,-40:40),CKIN(200)
8 COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
9 COMMON/PYINT1/MINT(400),VINT(400)
10 COMMON/PYINT2/ISET(200),KFPR(200,2),COEF(200,20),ICOL(40,4,2)
11 COMMON/PYINT3/XSFX(2,-40:40),ISIG(1000,3),SIGH(1000)
12 COMMON/PYINT4/WIDP(21:40,0:40),WIDE(21:40,0:40),WIDS(21:40,3)
13 COMMON/PYINT5/NGEN(0:200,3),XSEC(0:200,3)
14 COMMON/PYINT6/PROC(0:200)
15 CHARACTER PROC*28
16 COMMON/PYINT7/SIGT(0:6,0:6,0:5)
17 SAVE /PYSUBS/,/PYPARS/,/PYINT1/,/PYINT2/,/PYINT3/,/PYINT4/,
18 &/PYINT5/,/PYINT6/,/PYINT7/
19
20C...Default values for allowed processes and kinematics constraints.
21 DATA MSEL/1/
22 DATA MSUB/200*0/
23 DATA ((KFIN(I,J),J=-40,40),I=1,2)/16*0,4*1,4*0,6*1,5*0,5*1,0,
24 &5*1,5*0,6*1,4*0,4*1,16*0,16*0,4*1,4*0,6*1,5*0,5*1,0,5*1,5*0,
25 &6*1,4*0,4*1,16*0/
26 DATA CKIN/
27 & 2.0, -1.0, 0.0, -1.0, 1.0, 1.0, -10., 10., -10., 10.,
28 1 -10., 10., -10., 10., -10., 10., -1.0, 1.0, -1.0, 1.0,
29 2 0.0, 1.0, 0.0, 1.0, -1.0, 1.0, -1.0, 1.0, 0., 0.,
30 3 2.0, -1.0, 0., 0., 0.0, -1.0, 0.0, -1.0, 4.0, -1.0,
31 4 12.0, -1.0, 12.0, -1.0, 12.0, -1.0, 12.0, -1.0, 0., 0.,
32 5 0.0, -1.0, 0.0, -1.0, 0.0, -1.0, 0., 0., 0., 0.,
33 6 140*0./
34
35C...Default values for main switches and parameters. Reset information.
36 DATA (MSTP(I),I=1,100)/
37 & 3, 1, 2, 0, 0, 0, 0, 0, 0, 0,
38 1 1, 0, 1, 0, 5, 0, 0, 0, 0, 0,
39 2 1, 0, 1, 0, 0, 0, 0, 0, 0, 1,
40 3 1, 2, 0, 1, 0, 2, 1, 5, 2, 0,
41 4 1, 1, 3, 7, 3, 1, 1, 2, 2, 0,
42 5 9, 1, 1, 1, 5, 1, 1, 6, 1, 0,
43 6 1, 3, 2, 2, 1, 1, 2, 0, 0, 0,
44 7 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
45 8 1, 1, 100, 0, 0, 0, 0, 0, 0, 0,
46 9 1, 4, 1, 2, 0, 0, 0, 0, 0, 0/
47 DATA (MSTP(I),I=101,200)/
48 & 3, 1, 0, 0, 0, 0, 0, 0, 0, 0,
49 1 1, 1, 1, 0, 0, 0, 0, 0, 0, 0,
50 2 0, 1, 2, 1, 1, 20, 0, 0, 10, 0,
51 3 0, 4, 0, 1, 0, 0, 0, 0, 0, 0,
52 4 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
53 5 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
54 6 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
55 7 0, 2, 0, 0, 0, 0, 0, 0, 0, 0,
56 8 5, 720, 1995, 11, 29, 408, 0, 0, 0, 0,
57 9 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
58 DATA (PARP(I),I=1,100)/
59 & 0.25, 10., 0., 0., 0., 0., 0., 0., 0., 0.,
60 1 0., 0., 1.0, 0.01, 0.6, 1.0, 1.0, 0., 0., 0.,
61 2 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
62 3 1.5, 2.0, 0.075, 1.0, 0.2, 0., 2.0, 0.70, 0.006, 0.,
63 4 0.02, 2.0, 0.10, 1000., 2054., 123., 246., 0., 0., 0.,
64 5 1.0, 0., 0., 0., 0., 0., 0., 0., 0., 0.,
65 6 0.25, 1.0, 0.25, 1.0, 2.0, 1E-3, 4.0, 1E-3, 0., 0.,
66 7 4.0, 0.25, 0., 0., 0., 0., 0., 0., 0., 0.,
67 8 1.40, 1.55, 0.5, 0.2, 0.33, 0.66, 0.7, 0.5, 0., 0.,
68 9 0.44, 0.20, 2.0, 1.0, 0., 3.0, 1.0, 0.75, 0.44, 2.0/
69 DATA (PARP(I),I=101,200)/
70 & 0.5, 0.28, 1.0, 0., 0., 0., 0., 0., 0., 0.,
71 1 2.0, 0., 0., 0., 0., 0., 0., 0., 0., 0.,
72 2 1.0, 0.4, 0., 0., 0., 0., 0., 0., 0., 0.,
73 3 0.01, 0., 0., 0., 0., 0., 0., 0., 0., 0.,
74 4 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
75 5 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
76 6 2.20, 23.6, 18.4, 11.5, 0., 0., 0., 0., 0., 0.,
77 7 0., 0., 0., 1.0, 0., 0., 0., 0., 0., 0.,
78 8 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
79 9 0., 0., 0., 0., 0., 0., 0., 0., 0., 0./
80 DATA MSTI/200*0/
81 DATA PARI/200*0./
82 DATA MINT/400*0/
83 DATA VINT/400*0./
84
85C...Constants for the generation of the various processes.
86 DATA (ISET(I),I=1,100)/
87 & 1, 1, 1, -1, 3, -1, -1, 3, -2, 2,
88 1 2, 2, 2, 2, 2, 2, -1, 2, 2, 2,
89 2 -1, 2, 2, 2, 2, 2, -1, 2, 2, 2,
90 3 2, -1, 2, 2, 2, 2, -1, -1, -1, -1,
91 4 -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
92 5 -1, -1, 2, 2, -1, -1, -1, 2, -1, -1,
93 6 -1, -1, -1, -1, -1, -1, -1, 2, 2, 2,
94 7 4, 4, 4, -1, -1, 4, 4, -1, -1, 2,
95 8 2, 2, 2, 2, 2, 2, 2, 2, 2, -2,
96 9 0, 0, 0, 0, 0, 9, -2, -2, -2, -2/
97 DATA (ISET(I),I=101,200)/
98 & -1, 1, 1, -2, -2, -2, -2, -2, -2, 2,
99 1 2, 2, 2, 2, 2, -1, -1, -1, -2, -2,
100 2 5, 5, 5, 5, -2, -2, -2, -2, -2, -2,
101 3 6, -2, -2, -2, -2, -2, -2, -2, -2, -2,
102 4 1, 1, 1, 1, 1, -2, 1, 1, 1, -2,
103 5 1, 1, 1, -2, -2, 1, 1, 1, -2, -2,
104 6 2, 2, 2, 2, 2, 2, 2, 2, -2, -2,
105 7 2, 2, 5, 5, -2, 2, 2, 5, 5, -2,
106 8 5, 5, -2, -2, -2, 5, 5, -2, -2, -2,
107 9 -2, -2, -2, -2, -2, -2, -2, -2, -2, -2/
108 DATA ((KFPR(I,J),J=1,2),I=1,50)/
109 & 23, 0, 24, 0, 25, 0, 24, 0, 25, 0,
110 & 24, 0, 23, 0, 25, 0, 0, 0, 0, 0,
111 1 0, 0, 0, 0, 21, 21, 21, 22, 21, 23,
112 1 21, 24, 21, 25, 22, 22, 22, 23, 22, 24,
113 2 22, 25, 23, 23, 23, 24, 23, 25, 24, 24,
114 2 24, 25, 25, 25, 0, 21, 0, 22, 0, 23,
115 3 0, 24, 0, 25, 0, 21, 0, 22, 0, 23,
116 3 0, 24, 0, 25, 0, 21, 0, 22, 0, 23,
117 4 0, 24, 0, 25, 0, 21, 0, 22, 0, 23,
118 4 0, 24, 0, 25, 0, 21, 0, 22, 0, 23/
119 DATA ((KFPR(I,J),J=1,2),I=51,100)/
120 5 0, 24, 0, 25, 0, 0, 0, 0, 0, 0,
121 5 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
122 6 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
123 6 0, 0, 0, 0, 21, 21, 24, 24, 23, 24,
124 7 23, 23, 24, 24, 23, 24, 23, 25, 22, 22,
125 7 23, 23, 24, 24, 24, 25, 25, 25, 0, 211,
126 8 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
127 8 443, 21,10441, 21,20443, 21, 445, 21, 0, 0,
128 9 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
129 9 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
130 DATA ((KFPR(I,J),J=1,2),I=101,150)/
131 & 23, 0, 25, 0, 25, 0, 0, 0, 0, 0,
132 & 0, 0, 0, 0, 0, 0, 0, 0, 22, 25,
133 1 21, 25, 0, 25, 21, 25, 22, 22, 21, 22,
134 1 22, 23, 23, 23, 24, 24, 0, 0, 0, 0,
135 2 25, 6, 25, 6, 25, 0, 25, 0, 0, 0,
136 2 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
137 3 23, 5, 0, 0, 0, 0, 0, 0, 0, 0,
138 3 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
139 4 32, 0, 34, 0, 37, 0, 40, 0, 39, 0,
140 4 0, 0, 7, 0, 8, 0, 38, 0, 0, 0/
141 DATA ((KFPR(I,J),J=1,2),I=151,200)/
142 5 35, 0, 35, 0, 35, 0, 0, 0, 0, 0,
143 5 36, 0, 36, 0, 36, 0, 0, 0, 0, 0,
144 6 6, 37, 39, 0, 39, 39, 39, 39, 11, 0,
145 6 11, 0, 0, 7, 0, 8, 0, 0, 0, 0,
146 7 23, 35, 24, 35, 35, 0, 35, 0, 0, 0,
147 7 23, 36, 24, 36, 36, 0, 36, 0, 0, 0,
148 8 35, 6, 35, 6, 0, 0, 0, 0, 0, 0,
149 8 36, 6, 36, 6, 0, 0, 0, 0, 0, 0,
150 9 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
151 9 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
152 DATA COEF/4000*0./
153 DATA (((ICOL(I,J,K),K=1,2),J=1,4),I=1,40)/
154 1 4,0,3,0,2,0,1,0,3,0,4,0,1,0,2,0,2,0,0,1,4,0,0,3,3,0,0,4,1,0,0,2,
155 2 3,0,0,4,1,4,3,2,4,0,0,3,4,2,1,3,2,0,4,1,4,0,2,3,4,0,3,4,2,0,1,2,
156 3 3,2,1,0,1,4,3,0,4,3,3,0,2,1,1,0,3,2,1,4,1,0,0,2,2,4,3,1,2,0,0,1,
157 4 3,2,1,4,1,4,3,2,4,2,1,3,4,2,1,3,3,4,4,3,1,2,2,1,2,0,3,1,2,0,0,0,
158 5 4,2,1,0,0,0,1,0,3,0,0,3,1,2,0,0,4,0,0,4,0,0,1,2,2,0,0,1,4,4,3,3,
159 6 2,2,1,1,4,4,3,3,3,3,4,4,1,1,2,2,3,2,1,3,1,2,0,0,4,2,1,4,0,0,1,2,
160 7 4,0,0,0,4,0,1,3,0,0,3,0,2,4,3,0,3,4,0,0,1,0,0,1,0,0,3,4,2,0,0,2,
161 8 3,0,0,0,1,0,0,0,0,0,3,0,2,0,0,0,2,0,3,1,2,0,0,0,3,2,1,0,1,0,0,0,
162 9 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
163 & 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0/
164
165C...Character constants: name of processes.
166 DATA PROC(0)/ 'All included subprocesses '/
167 DATA (PROC(I),I=1,20)/
168 1'f + f~ -> gamma*/Z0 ', 'f + f~'' -> W+/- ',
169 2'f + f~ -> H0 ', 'gamma + W+/- -> W+/- ',
170 3'Z0 + Z0 -> H0 ', 'Z0 + W+/- -> W+/- ',
171 4' ', 'W+ + W- -> H0 ',
172 5' ', 'f + f'' -> f + f'' (QFD) ',
173 6'f + f'' -> f + f'' (QCD) ','f + f~ -> f'' + f~'' ',
174 7'f + f~ -> g + g ', 'f + f~ -> g + gamma ',
175 8'f + f~ -> g + Z0 ', 'f + f~'' -> g + W+/- ',
176 9'f + f~ -> g + H0 ', 'f + f~ -> gamma + gamma ',
177 &'f + f~ -> gamma + Z0 ', 'f + f~'' -> gamma + W+/- '/
178 DATA (PROC(I),I=21,40)/
179 1'f + f~ -> gamma + H0 ', 'f + f~ -> Z0 + Z0 ',
180 2'f + f~'' -> Z0 + W+/- ', 'f + f~ -> Z0 + H0 ',
181 3'f + f~ -> W+ + W- ', 'f + f~'' -> W+/- + H0 ',
182 4'f + f~ -> H0 + H0 ', 'f + g -> f + g ',
183 5'f + g -> f + gamma ', 'f + g -> f + Z0 ',
184 6'f + g -> f'' + W+/- ', 'f + g -> f + H0 ',
185 7'f + gamma -> f + g ', 'f + gamma -> f + gamma ',
186 8'f + gamma -> f + Z0 ', 'f + gamma -> f'' + W+/- ',
187 9'f + gamma -> f + H0 ', 'f + Z0 -> f + g ',
188 &'f + Z0 -> f + gamma ', 'f + Z0 -> f + Z0 '/
189 DATA (PROC(I),I=41,60)/
190 1'f + Z0 -> f'' + W+/- ', 'f + Z0 -> f + H0 ',
191 2'f + W+/- -> f'' + g ', 'f + W+/- -> f'' + gamma ',
192 3'f + W+/- -> f'' + Z0 ', 'f + W+/- -> f'' + W+/- ',
193 4'f + W+/- -> f'' + H0 ', 'f + H0 -> f + g ',
194 5'f + H0 -> f + gamma ', 'f + H0 -> f + Z0 ',
195 6'f + H0 -> f'' + W+/- ', 'f + H0 -> f + H0 ',
196 7'g + g -> f + f~ ', 'g + gamma -> f + f~ ',
197 8'g + Z0 -> f + f~ ', 'g + W+/- -> f + f~'' ',
198 9'g + H0 -> f + f~ ', 'gamma + gamma -> f + f~ ',
199 &'gamma + Z0 -> f + f~ ', 'gamma + W+/- -> f + f~'' '/
200 DATA (PROC(I),I=61,80)/
201 1'gamma + H0 -> f + f~ ', 'Z0 + Z0 -> f + f~ ',
202 2'Z0 + W+/- -> f + f~'' ', 'Z0 + H0 -> f + f~ ',
203 3'W+ + W- -> f + f~ ', 'W+/- + H0 -> f + f~'' ',
204 4'H0 + H0 -> f + f~ ', 'g + g -> g + g ',
205 5'gamma + gamma -> W+ + W- ', 'gamma + W+/- -> Z0 + W+/- ',
206 6'Z0 + Z0 -> Z0 + Z0 ', 'Z0 + Z0 -> W+ + W- ',
207 7'Z0 + W+/- -> Z0 + W+/- ', 'Z0 + Z0 -> Z0 + H0 ',
208 8'W+ + W- -> gamma + gamma ', 'W+ + W- -> Z0 + Z0 ',
209 9'W+/- + W+/- -> W+/- + W+/- ', 'W+/- + H0 -> W+/- + H0 ',
210 &'H0 + H0 -> H0 + H0 ', 'q + gamma -> q'' + pi+/- '/
211 DATA (PROC(I),I=81,100)/
212 1'q + q~ -> Q + Q~, massive ', 'g + g -> Q + Q~, massive ',
213 2'f + q -> f'' + Q, massive ', 'g + gamma -> Q + Q~, massive',
214 3'gamma + gamma -> F + F~, mas', 'g + g -> J/Psi + g ',
215 4'g + g -> chi_0c + g ', 'g + g -> chi_1c + g ',
216 5'g + g -> chi_2c + g ', ' ',
217 6'Elastic scattering ', 'Single diffractive (XB) ',
218 7'Single diffractive (AX) ', 'Double diffractive ',
219 8'Low-pT scattering ', 'Semihard QCD 2 -> 2 ',
220 9' ', ' ',
221 &' ', ' '/
222 DATA (PROC(I),I=101,120)/
223 1'g + g -> gamma*/Z0 ', 'g + g -> H0 ',
224 2'gamma + gamma -> H0 ', ' ',
225 3' ', ' ',
226 4' ', ' ',
227 5' ', 'f + f~ -> gamma + H0 ',
228 6'f + f~ -> g + H0 ', 'q + g -> q + H0 ',
229 7'g + g -> g + H0 ', 'g + g -> gamma + gamma ',
230 8'g + g -> g + gamma ', 'g + g -> gamma + Z0 ',
231 9'g + g -> Z0 + Z0 ', 'g + g -> W+ + W- ',
232 &' ', ' '/
233 DATA (PROC(I),I=121,140)/
234 1'g + g -> Q + Q~ + H0 ', 'q + q~ -> Q + Q~ + H0 ',
235 2'f + f'' -> f + f'' + H0 ',
236 2'f + f'' -> f" + f"'' + H0 ',
237 3' ', ' ',
238 4' ', ' ',
239 5' ', ' ',
240 6'g + g -> Z0 + q + q~ ', ' ',
241 7' ', ' ',
242 8' ', ' ',
243 9' ', ' ',
244 &' ', ' '/
245 DATA (PROC(I),I=141,160)/
246 1'f + f~ -> gamma*/Z0/Z''0 ', 'f + f~'' -> W''+/- ',
247 2'f + f~'' -> H+/- ', 'f + f~'' -> R ',
248 3'q + l -> LQ ', ' ',
249 4'd + g -> d* ', 'u + g -> u* ',
250 5'g + g -> eta_techni ', ' ',
251 6'f + f~ -> H''0 ', 'g + g -> H''0 ',
252 7'gamma + gamma -> H''0 ', ' ',
253 8' ', 'f + f~ -> A0 ',
254 9'g + g -> A0 ', 'gamma + gamma -> A0 ',
255 &' ', ' '/
256 DATA (PROC(I),I=161,180)/
257 1'f + g -> f'' + H+/- ', 'q + g -> LQ + l~ ',
258 2'g + g -> LQ + LQ~ ', 'q + q~ -> LQ + LQ~ ',
259 3'f + f~ -> f'' + f~'' (gamma/Z)',
260 3'f +f~'' -> f" + f~"'' (W) ',
261 4'q + q'' -> q" + d* ', 'q + q'' -> q" + u* ',
262 5' ', ' ',
263 6'f + f~ -> Z0 + H''0 ', 'f + f~'' -> W+/- + H''0 ',
264 7'f + f'' -> f + f'' + H''0 ',
265 7'f + f'' -> f" + f"'' + H''0 ',
266 8' ', 'f + f~ -> Z0 + A0 ',
267 9'f + f~'' -> W+/- + A0 ',
268 9'f + f'' -> f + f'' + A0 ',
269 &'f + f'' -> f" + f"'' + A0 ',
270 &' '/
271 DATA (PROC(I),I=181,200)/
272 1'g + g -> Q + Q~ + H''0 ', 'q + q~ -> Q + Q~ + H''0 ',
273 2' ', ' ',
274 3' ', 'g + g -> Q + Q~ + A0 ',
275 4'q + q~ -> Q + Q~ + A0 ', ' ',
276 5' ', ' ',
277 6' ', ' ',
278 7' ', ' ',
279 8' ', ' ',
280 9' ', ' ',
281 &' ', ' '/
282
283C...Cross sections and slope offsets.
284 DATA SIGT/294*0./
285
286 END