Instantiation of AliConfig removed
[u/mrichter/AliRoot.git] / ALIROOT / pydata.F
1  
2 C*********************************************************************
3  
4       BLOCK DATA PYDATA
5  
6 C...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  
20 C...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  
35 C...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  
85 C...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  
165 C...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  
283 C...Cross sections and slope offsets.
284       DATA SIGT/294*0./
285  
286       END