]> git.uio.no Git - u/mrichter/AliRoot.git/blame - HIJING/hipyset1_35/pydata_hijing.F
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Numerical stability
[u/mrichter/AliRoot.git] / HIJING / hipyset1_35 / pydata_hijing.F
CommitLineData
e74335a4 1* $Id$
2
3C*********************************************************************
4
5 BLOCK DATA PYDATA_HIJING
6
7C...Give sensible default values to all status codes and parameters.
8#include "pysubs_hijing.inc"
9#include "pypars_hijing.inc"
10#include "pyint1_hijing.inc"
11#include "pyint2_hijing.inc"
12#include "pyint3_hijing.inc"
13#include "pyint4_hijing.inc"
14#include "pyint5_hijing.inc"
15#include "pyint6_hijing.inc"
16
17C...Default values for allowed processes and kinematics constraints.
18 DATA MSEL/1/
19 DATA MSUB/200*0/
20 DATA ((KFIN(I,J),J=-40,40),I=1,2)/40*1,0,80*1,0,40*1/
21 DATA CKIN/
22 & 2.0, -1.0, 0.0, -1.0, 1.0, 1.0, -10., 10., -10., 10.,
23 1 -10., 10., -10., 10., -10., 10., -1.0, 1.0, -1.0, 1.0,
24 2 0.0, 1.0, 0.0, 1.0, -1.0, 1.0, -1.0, 1.0, 0., 0.,
25 3 2.0, -1.0, 0., 0., 0., 0., 0., 0., 0., 0.,
26 4 160*0./
27
28C...Default values for main switches and parameters. Reset information.
29 DATA (MSTP(I),I=1,100)/
30 & 3, 1, 2, 0, 0, 0, 0, 0, 0, 0,
31 1 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
32 2 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
33 3 1, 2, 0, 0, 0, 2, 0, 0, 0, 0,
34 4 1, 0, 3, 7, 1, 0, 0, 0, 0, 0,
35 5 1, 1, 20, 6, 0, 0, 0, 0, 0, 0,
36 6 1, 2, 2, 2, 1, 0, 0, 0, 0, 0,
37 7 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
38 8 1, 1, 100, 0, 0, 0, 0, 0, 0, 0,
39 9 1, 4, 0, 0, 0, 0, 0, 0, 0, 0/
40 DATA (MSTP(I),I=101,200)/
41 & 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
42 1 1, 1, 1, 0, 0, 0, 0, 0, 0, 0,
43 2 0, 1, 2, 1, 1, 20, 0, 0, 0, 0,
44 3 0, 4, 0, 1, 0, 0, 0, 0, 0, 0,
45 4 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
46 5 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
47 6 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
48 7 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
49 8 5, 3, 1989, 11, 24, 0, 0, 0, 0, 0,
50 9 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
51 DATA (PARP(I),I=1,100)/
52 & 0.25, 10., 0., 0., 0., 0., 0., 0., 0., 0.,
53 1 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
54 2 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
55 3 1.5, 2.0, 0.075, 0., 0.2, 0., 0., 0., 0., 0.,
56 4 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
57 5 1.0, 2.26, 1.E4, 1.E-4, 0., 0., 0., 0., 0., 0.,
58 6 0.25, 1.0, 0.25, 1.0, 2.0, 1.E-3, 4.0, 0., 0., 0.,
59 7 4.0, 0., 0., 0., 0., 0., 0., 0., 0., 0.,
60 8 1.6, 1.85, 0.5, 0.2, 0.33, 0.66, 0.7, 0.5, 0., 0.,
61 9 0.44, 0.44, 2.0, 1.0, 0., 3.0, 1.0, 0.75, 0., 0./
62 DATA (PARP(I),I=101,200)/
63 & -0.02, 0., 0., 0., 0., 0., 0., 0., 0., 0.,
64 1 2.0, 0., 0., 0., 0., 0., 0., 0., 0., 0.,
65 2 0.4, 0., 0., 0., 0., 0., 0., 0., 0., 0.,
66 3 0.01, 0., 0., 0., 0., 0., 0., 0., 0., 0.,
67 4 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
68 5 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
69 6 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
70 7 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
71 8 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
72 9 0., 0., 0., 0., 0., 0., 0., 0., 0., 0./
73 DATA MSTI/200*0/
74 DATA PARI/200*0./
75 DATA MINT/400*0/
76 DATA VINT/400*0./
77
78C...Constants for the generation of the various processes.
79 DATA (ISET(I),I=1,100)/
80 & 1, 1, 1, -1, 3, -1, -1, 3, -2, -2,
81 1 2, 2, 2, 2, 2, 2, -1, 2, 2, 2,
82 2 -1, 2, 2, 2, 2, 2, -1, 2, 2, 2,
83 3 2, -1, -1, -1, -1, -1, -1, -1, -1, -1,
84 4 -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
85 5 -1, -1, 2, -1, -1, -1, -1, -1, -1, -1,
86 6 -1, -1, -1, -1, -1, -1, -1, 2, -1, -1,
87 7 4, 4, 4, -1, -1, 4, 4, -1, -1, -2,
88 8 2, 2, -2, -2, -2, -2, -2, -2, -2, -2,
89 9 0, 0, 0, -1, 0, 5, -2, -2, -2, -2/
90 DATA (ISET(I),I=101,200)/
91 & -1, 1, -2, -2, -2, -2, -2, -2, -2, -2,
92 1 2, 2, 2, 2, -1, -1, -1, -2, -2, -2,
93 2 -1, -2, -2, -2, -2, -2, -2, -2, -2, -2,
94 3 -2, -2, -2, -2, -2, -2, -2, -2, -2, -2,
95 4 1, 1, 1, -2, -2, -2, -2, -2, -2, -2,
96 5 -2, -2, -2, -2, -2, -2, -2, -2, -2, -2,
97 6 2, -2, -2, -2, -2, -2, -2, -2, -2, -2,
98 7 -2, -2, -2, -2, -2, -2, -2, -2, -2, -2,
99 8 -2, -2, -2, -2, -2, -2, -2, -2, -2, -2,
100 9 -2, -2, -2, -2, -2, -2, -2, -2, -2, -2/
101 DATA ((KFPR(I,J),J=1,2),I=1,50)/
102 & 23, 0, 24, 0, 25, 0, 24, 0, 25, 0,
103 & 24, 0, 23, 0, 25, 0, 0, 0, 0, 0,
104 1 0, 0, 0, 0, 21, 21, 21, 22, 21, 23,
105 1 21, 24, 21, 25, 22, 22, 22, 23, 22, 24,
106 2 22, 25, 23, 23, 23, 24, 23, 25, 24, 24,
107 2 24, 25, 25, 25, 0, 21, 0, 22, 0, 23,
108 3 0, 24, 0, 25, 0, 21, 0, 22, 0, 23,
109 3 0, 24, 0, 25, 0, 21, 0, 22, 0, 23,
110 4 0, 24, 0, 25, 0, 21, 0, 22, 0, 23,
111 4 0, 24, 0, 25, 0, 21, 0, 22, 0, 23/
112 DATA ((KFPR(I,J),J=1,2),I=51,100)/
113 5 0, 24, 0, 25, 0, 0, 0, 0, 0, 0,
114 5 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
115 6 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
116 6 0, 0, 0, 0, 21, 21, 24, 24, 22, 24,
117 7 23, 23, 24, 24, 23, 24, 23, 25, 22, 22,
118 7 23, 23, 24, 24, 24, 25, 25, 25, 0, 0,
119 8 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
120 8 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
121 9 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
122 9 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
123 DATA ((KFPR(I,J),J=1,2),I=101,150)/
124 & 23, 0, 25, 0, 0, 0, 0, 0, 0, 0,
125 & 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
126 1 21, 25, 0, 25, 21, 25, 22, 22, 22, 23,
127 1 23, 23, 24, 24, 0, 0, 0, 0, 0, 0,
128 2 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
129 2 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
130 3 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
131 3 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
132 4 32, 0, 37, 0, 40, 0, 0, 0, 0, 0,
133 4 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
134 DATA ((KFPR(I,J),J=1,2),I=151,200)/
135 5 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
136 5 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
137 6 0, 37, 0, 0, 0, 0, 0, 0, 0, 0,
138 6 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
139 7 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
140 7 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
141 8 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
142 8 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
143 9 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
144 9 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
145 DATA COEF/4000*0./
146 DATA (((ICOL(I,J,K),K=1,2),J=1,4),I=1,40)/
147 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,
148 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,
149 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,
150 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,
151 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,
152 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,
153 7 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,
154 8 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,
155 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,
156 & 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/
157
158C...Character constants: name of processes.
159 DATA PROC(0)/ 'All included subprocesses '/
160 DATA (PROC(I),I=1,20)/
161 1'f + fb -> gamma*/Z0 ', 'f + fb'' -> W+/- ',
162 2'f + fb -> H0 ', 'gamma + W+/- -> W+/- ',
163 3'Z0 + Z0 -> H0 ', 'Z0 + W+/- -> W+/- ',
164 4' ', 'W+ + W- -> H0 ',
165 5' ', ' ',
166 6'f + f'' -> f + f'' ','f + fb -> f'' + fb'' ',
167 7'f + fb -> g + g ', 'f + fb -> g + gamma ',
168 8'f + fb -> g + Z0 ', 'f + fb'' -> g + W+/- ',
169 9'f + fb -> g + H0 ', 'f + fb -> gamma + gamma ',
170 &'f + fb -> gamma + Z0 ', 'f + fb'' -> gamma + W+/- '/
171 DATA (PROC(I),I=21,40)/
172 1'f + fb -> gamma + H0 ', 'f + fb -> Z0 + Z0 ',
173 2'f + fb'' -> Z0 + W+/- ', 'f + fb -> Z0 + H0 ',
174 3'f + fb -> W+ + W- ', 'f + fb'' -> W+/- + H0 ',
175 4'f + fb -> H0 + H0 ', 'f + g -> f + g ',
176 5'f + g -> f + gamma ', 'f + g -> f + Z0 ',
177 6'f + g -> f'' + W+/- ', 'f + g -> f + H0 ',
178 7'f + gamma -> f + g ', 'f + gamma -> f + gamma ',
179 8'f + gamma -> f + Z0 ', 'f + gamma -> f'' + W+/- ',
180 9'f + gamma -> f + H0 ', 'f + Z0 -> f + g ',
181 &'f + Z0 -> f + gamma ', 'f + Z0 -> f + Z0 '/
182 DATA (PROC(I),I=41,60)/
183 1'f + Z0 -> f'' + W+/- ', 'f + Z0 -> f + H0 ',
184 2'f + W+/- -> f'' + g ', 'f + W+/- -> f'' + gamma ',
185 3'f + W+/- -> f'' + Z0 ', 'f + W+/- -> f'' + W+/- ',
186 4'f + W+/- -> f'' + H0 ', 'f + H0 -> f + g ',
187 5'f + H0 -> f + gamma ', 'f + H0 -> f + Z0 ',
188 6'f + H0 -> f'' + W+/- ', 'f + H0 -> f + H0 ',
189 7'g + g -> f + fb ', 'g + gamma -> f + fb ',
190 8'g + Z0 -> f + fb ', 'g + W+/- -> f + fb'' ',
191 9'g + H0 -> f + fb ', 'gamma + gamma -> f + fb ',
192 &'gamma + Z0 -> f + fb ', 'gamma + W+/- -> f + fb'' '/
193 DATA (PROC(I),I=61,80)/
194 1'gamma + H0 -> f + fb ', 'Z0 + Z0 -> f + fb ',
195 2'Z0 + W+/- -> f + fb'' ', 'Z0 + H0 -> f + fb ',
196 3'W+ + W- -> f + fb ', 'W+/- + H0 -> f + fb'' ',
197 4'H0 + H0 -> f + fb ', 'g + g -> g + g ',
198 5'gamma + gamma -> W+ + W- ', 'gamma + W+/- -> gamma + W+/-',
199 6'Z0 + Z0 -> Z0 + Z0 ', 'Z0 + Z0 -> W+ + W- ',
200 7'Z0 + W+/- -> Z0 + W+/- ', 'Z0 + Z0 -> Z0 + H0 ',
201 8'W+ + W- -> gamma + gamma ', 'W+ + W- -> Z0 + Z0 ',
202 9'W+/- + W+/- -> W+/- + W+/- ', 'W+/- + H0 -> W+/- + H0 ',
203 &'H0 + H0 -> H0 + H0 ', ' '/
204 DATA (PROC(I),I=81,100)/
205 1'q + qb -> Q + QB, massive ', 'g + g -> Q + QB, massive ',
206 2' ', ' ',
207 3' ', ' ',
208 4' ', ' ',
209 5' ', ' ',
210 6'Elastic scattering ', 'Single diffractive ',
211 7'Double diffractive ', 'Central diffractive ',
212 8'Low-pT scattering ', 'Semihard QCD 2 -> 2 ',
213 9' ', ' ',
214 &' ', ' '/
215 DATA (PROC(I),I=101,120)/
216 1'g + g -> gamma*/Z0 ', 'g + g -> H0 ',
217 2' ', ' ',
218 3' ', ' ',
219 4' ', ' ',
220 5' ', ' ',
221 6'f + fb -> g + H0 ', 'q + g -> q + H0 ',
222 7'g + g -> g + H0 ', 'g + g -> gamma + gamma ',
223 8'g + g -> gamma + Z0 ', 'g + g -> Z0 + Z0 ',
224 9'g + g -> W+ + W- ', ' ',
225 &' ', ' '/
226 DATA (PROC(I),I=121,140)/
227 1'g + g -> f + fb + H0 ', ' ',
228 2' ', ' ',
229 3' ', ' ',
230 4' ', ' ',
231 5' ', ' ',
232 6' ', ' ',
233 7' ', ' ',
234 8' ', ' ',
235 9' ', ' ',
236 &' ', ' '/
237 DATA (PROC(I),I=141,160)/
238 1'f + fb -> gamma*/Z0/Z''0 ', 'f + fb'' -> H+/- ',
239 2'f + fb -> R ', ' ',
240 3' ', ' ',
241 4' ', ' ',
242 5' ', ' ',
243 6' ', ' ',
244 7' ', ' ',
245 8' ', ' ',
246 9' ', ' ',
247 &' ', ' '/
248 DATA (PROC(I),I=161,180)/
249 1'f + g -> f'' + H+/- ', ' ',
250 2' ', ' ',
251 3' ', ' ',
252 4' ', ' ',
253 5' ', ' ',
254 6' ', ' ',
255 7' ', ' ',
256 8' ', ' ',
257 9' ', ' ',
258 &' ', ' '/
259 DATA (PROC(I),I=181,200)/ 20*' '/
260
261 END