-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathludata.f
447 lines (438 loc) · 28.3 KB
/
ludata.f
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
C*********************************************************************
BLOCK DATA LUDATA
C...Purpose: to give default values to parameters and particle and
C...decay data.
COMMON/LUDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
COMMON/LUDAT2/KCHG(500,3),PMAS(500,4),PARF(2000),VCKM(4,4)
COMMON/LUDAT3/MDCY(500,3),MDME(2000,2),BRAT(2000),KFDP(2000,5)
COMMON/LUDAT4/CHAF(500)
CHARACTER CHAF*8
COMMON/LUDATR/MRLU(6),RRLU(100)
SAVE /LUDAT1/,/LUDAT2/,/LUDAT3/,/LUDAT4/,/LUDATR/
C...LUDAT1, containing status codes and most parameters.
DATA MSTU/
& 0, 0, 0, 4000,10000, 500, 2000, 0, 0, 2,
1 6, 1, 1, 0, 1, 1, 0, 0, 0, 0,
2 2, 10, 0, 0, 1, 10, 0, 0, 0, 0,
3 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
4 2, 2, 1, 4, 2, 1, 1, 0, 0, 0,
5 25, 24, 0, 1, 0, 0, 0, 0, 0, 0,
6 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
7 30*0,
& 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1 1, 5, 3, 5, 0, 0, 0, 0, 0, 0,
2 60*0,
8 7, 410, 1997, 01, 20, 700, 0, 0, 0, 0,
9 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
DATA PARU/
& 3.1415927, 6.2831854, 0.1973, 5.068, 0.3894, 2.568, 4*0.,
1 0.001, 0.09, 0.01, 0., 0., 0., 0., 0., 0., 0.,
2 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
3 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
4 2.0, 1.0, 0.25, 2.5, 0.05, 0., 0., 0.0001, 0., 0.,
5 2.5, 1.5, 7.0, 1.0, 0.5, 2.0, 3.2, 0., 0., 0.,
6 40*0.,
& 0.00729735, 0.232, 0.007764, 1.0, 1.16639E-5, 0., 0., 0.,
& 0., 0.,
1 0.20, 0.25, 1.0, 4.0, 10., 0., 0., 0., 0., 0.,
2 -0.693, -1.0, 0.387, 1.0, -0.08, -1.0, 1.0, 1.0, 1.0, 0.,
3 1.0, -1.0, 1.0, -1.0, 1.0, 0., 0., 0., 0., 0.,
4 5.0, 1.0, 1.0, 0., 1.0, 1.0, 0., 0., 0., 0.,
5 1.0, 0., 0., 0., 1000., 1.0, 1.0, 1.0, 1.0, 0.,
6 1.0, 1.0, 1.0, 1.0, 1.0, 0., 0., 0., 0., 0.,
7 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 0., 0., 0.,
8 1.0, 1.0, 1.0, 0.0, 0.0, 1.0, 1.0, 0.0, 0.0, 0.,
9 0., 0., 0., 0., 1.0, 0., 0., 0., 0., 0./
DATA MSTJ/
& 1, 3, 0, 0, 0, 0, 0, 0, 0, 0,
1 4, 2, 0, 1, 0, 0, 0, 0, 0, 0,
2 2, 1, 1, 2, 1, 2, 2, 0, 0, 0,
3 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
4 2, 2, 4, 2, 5, 3, 3, 0, 0, 3,
5 0, 3, 0, 0, 0, 0, 0, 0, 0, 0,
6 40*0,
& 5, 2, 7, 5, 1, 1, 0, 2, 0, 2,
1 0, 0, 0, 0, 1, 1, 0, 0, 0, 0,
2 80*0/
DATA PARJ/
& 0.10, 0.30, 0.40, 0.05, 0.50, 0.50, 0.50, 0., 0., 0.,
1 0.50, 0.60, 0.75, 0., 0., 0., 0., 1.0, 1.0, 0.,
2 0.36, 1.0, 0.01, 2.0, 1.0, 0.4, 0., 0., 0., 0.,
3 0.10, 1.0, 0.8, 1.5, 0., 2.0, 0.2, 2.5, 0.6, 0.,
4 0.3, 0.58, 0.5, 0.9, 0.5, 1.0, 1.0, 1.0, 0., 0.,
5 0.77,0.77,0.77,-0.05,-0.005,-0.00001,-0.00001,-0.00001,1.0,0.,
6 4.5, 0.7, 0., 0.003, 0.5, 0.5, 0., 0., 0., 0.,
7 10., 1000., 100., 1000., 0., 0.7, 10., 0., 0., 0.,
8 0.29, 1.0, 1.0, 0., 10., 10., 0., 0., 0., 0.,
9 0.02, 1.0, 0.2, 0., 0., 0., 0., 0., 0., 0.,
& 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
1 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
2 1.0, 0.25,91.187,2.489, 0.01, 2.0, 1.0, 0.25,0.002, 0.,
3 0., 0., 0., 0., 0.01, 0.99, 0., 0., 0.2, 0.,
4 60*0./
C...LUDAT2, with particle data and flavour treatment parameters.
DATA (KCHG(I,1),I= 1, 500)/-1,2,-1,2,-1,2,-1,2,2*0,-3,0,-3,0,
&-3,0,-3,6*0,3,9*0,3,2*0,3,0,-1,44*0,2,-1,2,-1,2,3,11*0,3,0,2*3,0,
&3,0,3,0,3,10*0,3,0,2*3,0,3,0,3,0,3,10*0,3,0,2*3,0,3,0,3,0,3,10*0,
&3,0,2*3,0,3,0,3,0,3,10*0,3,0,2*3,0,3,0,3,0,3,10*0,3,0,2*3,0,3,0,
&3,0,3,70*0,3,0,3,28*0,3,2*0,3,8*0,-3,8*0,3,0,-3,0,3,-3,3*0,3,6,0,
&3,5*0,-3,0,3,-3,0,-3,4*0,-3,0,3,6,-3,0,3,-3,0,-3,0,3,6,0,3,5*0,
&-3,0,3,-3,0,-3,114*0/
DATA (KCHG(I,2),I= 1, 500)/8*1,12*0,2,16*0,2,1,50*0,-1,410*0/
DATA (KCHG(I,3),I= 1, 500)/8*1,2*0,8*1,5*0,1,9*0,1,2*0,1,0,2*1,
&41*0,1,0,7*1,10*0,10*1,10*0,10*1,10*0,10*1,10*0,10*1,10*0,10*1,
&10*0,10*1,70*0,3*1,22*0,1,5*0,1,0,2*1,6*0,1,0,2*1,6*0,2*1,0,5*1,
&0,6*1,4*0,6*1,4*0,16*1,4*0,6*1,114*0/
DATA (PMAS(I,1),I= 1, 500)/0.0099,0.0056,0.199,1.35,5.,160.,
&2*250.,2*0.,0.00051,0.,0.1057,0.,1.777,0.,250.,5*0.,91.187,80.25,
&80.,6*0.,500.,900.,500.,3*300.,350.,200.,5000.,60*0.,0.1396,
&0.4977,0.4936,1.8693,1.8645,1.9688,5.2787,5.2786,5.47972,6.594,
&0.135,0.5475,0.9578,2.9788,9.4,320.,2*500.,2*0.,0.7669,0.8961,
&0.8916,2.0101,2.0071,2.11,2*5.325,5.5068,6.602,0.7683,0.782,
&1.0194,3.0969,9.4603,320.,2*500.,2*0.,1.232,2*1.29,2*2.424,2.536,
&2*5.73,5.97,7.3,1.232,1.17,1.4,3.46,9.875,320.,2*500.,2*0.,0.983,
&2*1.429,2*2.272,2.5,2*5.68,5.92,7.25,0.9827,1.,1.4,3.4151,9.8598,
&320.,2*500.,2*0.,1.26,2*1.402,2*2.372,2.56,2*5.78,6.02,7.3,1.26,
&1.282,1.42,3.5106,9.8919,320.,2*500.,2*0.,1.318,1.432,1.425,
&2*2.46,2.61,2*5.83,6.07,7.35,1.318,1.275,1.525,3.5562,9.9132,
&320.,2*500.,2*0.,2*0.4977,8*0.,3.686,3*0.,10.0233,70*0.,1.1156,
&5*0.,2.2849,0.,2.473,2.466,6*0.,5.641,0.,2*5.84,6*0.,0.9396,
&0.9383,0.,1.1974,1.1926,1.1894,1.3213,1.3149,0.,2.4525,2.4529,
&2.4527,2*2.55,2.73,4*0.,3*5.8,2*5.96,6.12,4*0.,1.234,1.233,1.232,
&1.231,1.3872,1.3837,1.3828,1.535,1.5318,1.6724,3*2.5,2*2.63,2.8,
&4*0.,3*5.81,2*5.97,6.13,114*0./
DATA (PMAS(I,2),I= 1, 500)/22*0.,2.489,2.066,88*0.,0.0002,
&0.001,6*0.,0.149,0.0505,0.0498,7*0.,0.151,0.00843,0.0044,7*0.,
&0.155,2*0.09,2*0.02,0.,4*0.05,0.155,0.36,0.08,2*0.01,5*0.,0.057,
&2*0.287,7*0.05,0.057,0.,0.25,0.014,6*0.,0.4,2*0.174,7*0.05,0.4,
&0.024,0.06,0.0009,6*0.,0.11,0.109,0.098,2*0.019,5*0.02,0.11,
&0.185,0.076,0.002,146*0.,4*0.12,0.0394,0.036,0.0358,0.0099,
&0.0091,131*0./
DATA (PMAS(I,3),I= 1, 500)/22*0.,2*20.,88*0.,0.002,0.005,6*0.,
&0.4,2*0.2,7*0.,0.4,0.1,0.015,7*0.,0.25,0.005,0.01,2*0.08,0.,
&4*0.1,0.25,0.2,0.001,2*0.02,5*0.,0.05,2*0.4,6*0.1,2*0.05,0.,0.35,
&0.05,6*0.,3*0.3,2*0.1,0.03,4*0.1,0.3,0.05,0.02,0.001,6*0.,0.25,
&4*0.12,5*0.05,0.25,0.17,0.2,0.01,146*0.,4*0.14,0.04,2*0.035,
&2*0.05,131*0./
DATA (PMAS(I,4),I= 1, 500)/12*0.,658650.,0.,0.0914,68*0.,0.1,
&0.387,15*0.,7804.,0.,3709.,0.32,0.1259,0.135,3*0.387,0.15,110*0.,
&15500.,26.75,83*0.,78.88,5*0.,0.057,0.,0.025,0.09,6*0.,0.387,0.,
&2*0.387,9*0.,44.3,0.,23.95,49.1,86.9,6*0.,0.13,9*0.,0.387,13*0.,
&24.60001,130*0./
DATA PARF/
& 0.5, 0.25, 0.5, 0.25, 1., 0.5, 0., 0., 0., 0.,
1 0.5, 0., 0.5, 0., 1., 1., 0., 0., 0., 0.,
2 0.5, 0., 0.5, 0., 1., 1., 0., 0., 0., 0.,
3 0.5, 0., 0.5, 0., 1., 1., 0., 0., 0., 0.,
4 0.5, 0., 0.5, 0., 1., 1., 0., 0., 0., 0.,
5 0.5, 0., 0.5, 0., 1., 1., 0., 0., 0., 0.,
6 0.75, 0.5, 0., 0.1667, 0.0833, 0.1667, 0., 0., 0., 0.,
7 0., 0., 1., 0.3333, 0.6667, 0.3333, 0., 0., 0., 0.,
8 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
9 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
& 0.325, 0.325, 0.5, 1.6, 5.0, 0., 0., 0., 0., 0.,
1 0., 0.11, 0.16, 0.048, 0.50, 0.45, 0.55, 0.60, 0., 0.,
2 0.2, 0.1, 0., 0., 0., 0., 0., 0., 0., 0.,
3 1870*0./
DATA ((VCKM(I,J),J=1,4),I=1,4)/
1 0.95113, 0.04884, 0.00003, 0.00000,
2 0.04884, 0.94940, 0.00176, 0.00000,
3 0.00003, 0.00176, 0.99821, 0.00000,
4 0.00000, 0.00000, 0.00000, 1.00000/
C...LUDAT3, with particle decay parameters and data.
DATA (MDCY(I,1),I= 1, 500)/5*0,3*1,6*0,1,0,1,5*0,3*1,6*0,1,0,1,
&2*0,4*1,42*0,7*1,12*0,1,0,15*1,2*0,18*1,2*0,18*1,2*0,18*1,2*0,
&18*1,2*0,18*1,3*0,1,8*0,1,3*0,1,70*0,1,5*0,1,0,2*1,6*0,1,0,2*1,
&9*0,5*1,0,6*1,4*0,6*1,4*0,16*1,4*0,6*1,114*0/
DATA (MDCY(I,2),I= 1, 500)/1,9,17,25,33,41,50,60,2*0,70,74,76,
&81,83,124,126,132,2*0,135,144,156,172,192,6*0,209,0,231,254,274,
&292,301,304,305,42*0,314,315,319,328,331,336,338,11*0,358,359,
&361,367,430,491,524,560,596,635,666,668,675,681,682,683,684,685,
&2*0,686,688,691,694,697,699,700,701,702,703,704,708,713,721,724,
&733,734,735,2*0,736,737,742,747,749,751,753,755,757,759,761,762,
&765,769,770,771,772,773,2*0,774,775,777,779,781,783,785,787,789,
&791,793,794,799,804,806,808,809,810,2*0,811,813,815,817,819,821,
&823,825,827,829,831,833,846,850,852,854,855,856,2*0,857,863,873,
&884,892,900,904,912,920,924,928,936,945,951,953,955,956,957,2*0,
&958,966,8*0,968,3*0,979,70*0,993,5*0,997,0,1073,1074,6*0,1075,0,
&1092,1093,9*0,1094,1096,1097,1100,1101,0,1103,1104,1105,1106,
&1107,1108,4*0,1109,1110,1111,1112,1113,1114,4*0,1115,1116,1119,
&1122,1123,1126,1129,1132,1134,1136,1140,1141,1142,1143,1145,1147,
&4*0,1148,1149,1150,1151,1152,1153,114*0/
DATA (MDCY(I,3),I= 1, 500)/5*8,9,2*10,2*0,4,2,5,2,41,2,6,3,2*0,
&9,12,16,20,17,6*0,22,0,23,20,18,9,3,1,9,42*0,1,4,9,3,5,2,20,11*0,
&1,2,6,63,61,33,2*36,39,31,2,7,6,5*1,2*0,2,3*3,2,5*1,4,5,8,3,9,
&3*1,2*0,1,2*5,7*2,1,3,4,5*1,2*0,1,9*2,1,2*5,2*2,3*1,2*0,11*2,13,
&4,2*2,3*1,2*0,6,10,11,2*8,4,2*8,2*4,8,9,6,2*2,3*1,2*0,8,2,8*0,11,
&3*0,14,70*0,4,5*0,76,0,2*1,6*0,17,0,2*1,9*0,2,1,3,1,2,0,6*1,4*0,
&6*1,4*0,1,2*3,1,3*3,2*2,4,3*1,2*2,1,4*0,6*1,114*0/
DATA (MDME(I,1),I= 1,2000)/6*1,-1,7*1,-1,7*1,-1,7*1,-1,7*1,-1,
&7*1,-1,1,-1,8*1,2*-1,8*1,2*-1,61*1,-1,2*1,-1,6*1,2*-1,7*1,2*-1,
&3*1,-1,6*1,2*-1,6*1,2*-1,3*1,-1,3*1,-1,3*1,5*-1,3*1,-1,6*1,2*-1,
&3*1,-1,11*1,2*-1,6*1,8*-1,3*1,-1,3*1,-1,3*1,5*-1,3*1,4*-1,6*1,
&2*-1,3*1,-1,5*1,-1,8*1,2*-1,3*1,-1,9*1,-1,3*1,-1,9*1,2*-1,2*1,-1,
&16*1,-1,2*1,3*-1,1665*1/
DATA (MDME(I,2),I= 1,2000)/75*102,42,6*102,2*42,2*0,7*41,2*0,
&24*41,6*102,45,29*102,8*32,8*0,16*32,4*0,8*32,4*0,32,4*0,8*32,
&14*0,16*32,7*0,8*32,4*0,32,7*0,8*32,4*0,32,5*0,4*32,5*0,3*32,0,
&6*32,3*0,12,2*42,2*11,9*42,2*45,31,2*45,2*33,31,2*45,20*46,7*0,
&24*42,41*0,16*42,46*0,10*42,20*0,2*13,14*42,16*0,48,3*13,16*42,
&16*0,48,3*13,16*42,19*0,48,3*13,2*42,0,2*11,28*42,0,2,4*0,2,8*0,
&12,32,86,87,88,3,0,2*3,0,2*3,0,2*3,0,3,6*0,3,3*0,1,0,3,2*0,2*3,
&3*0,1,4*0,12,3*0,4*32,2*4,86,87,88,33*0,12,32,86,87,88,31*0,12,0,
&32,86,87,88,40*0,12,0,32,86,87,88,95*0,12,0,32,86,87,88,2*0,4*42,
&6*0,12,11*0,4*32,2*4,9*0,14*42,52*0,10*13,2*84,3*42,8*0,48,3*13,
&2*42,2*85,14*0,84,5*0,85,886*0/
DATA (BRAT(I) ,I= 1, 439)/75*0.,1.,6*0.,0.179,0.178,0.116,
&0.235,0.005,0.056,0.018,0.023,0.011,2*0.004,0.0067,0.014,2*0.002,
&2*0.001,0.0022,0.054,0.002,0.016,0.005,0.011,0.0101,5*0.006,
&0.002,2*0.001,5*0.002,6*0.,1.,29*0.,0.15394,0.11936,0.15394,
&0.11926,0.15254,3*0.,0.03368,0.06664,0.03368,0.06664,0.03368,
&0.06664,2*0.,0.3214,0.0165,2*0.,0.0165,0.3207,2*0.,0.00001,
&0.00059,6*0.,3*0.1081,3*0.,0.0003,0.048,0.8705,4*0.,0.0002,
&0.0603,0.,0.0199,0.0008,3*0.,0.143,0.111,0.143,0.111,0.143,0.085,
&2*0.,0.03,0.058,0.03,0.058,0.03,0.058,8*0.,0.25,0.01,2*0.,0.01,
&0.25,4*0.,0.24,5*0.,3*0.08,6*0.,0.01,0.08,0.82,5*0.,0.09,11*0.,
&0.01,0.08,0.82,5*0.,0.09,9*0.,1.,6*0.,0.01,0.98,0.01,1.,4*0.215,
&2*0.,2*0.07,0.,1.,2*0.08,0.76,0.08,2*0.105,0.04,0.5,0.08,0.14,
&0.01,0.015,0.005,1.,3*0.,1.,4*0.,1.,0.25,0.01,2*0.,0.01,0.25,
&4*0.,0.24,5*0.,3*0.08,0.,1.,2*0.5,0.635,0.212,0.056,0.017,0.048,
&0.032,0.07,0.065,2*0.005,2*0.011,5*0.001,0.07,0.065,2*0.005,
&2*0.011,5*0.001,0.026,0.019,0.066,0.041,0.045,0.076,0.0073,
&2*0.0047,0.026,0.001,0.0006,0.0066,0.005,2*0.003,2*0.0006,
&2*0.001,0.006,0.005,0.012,0.0057,0.067,0.008,0.0022,0.027,0.004,
&0.019,0.012,0.002,0.009,0.0218,0.001,0.022,0.087,0.001,0.0019,
&0.0015,0.0028,0.034,0.027,2*0.002,2*0.004,2*0.002,0.034,0.027/
DATA (BRAT(I) ,I= 440, 655)/2*0.002,2*0.004,2*0.002,0.0365,
&0.045,0.073,0.062,3*0.021,0.0061,0.015,0.025,0.0088,0.074,0.0109,
&0.0041,0.002,0.0035,0.0011,0.001,0.0027,2*0.0016,0.0018,0.011,
&0.0063,0.0052,0.018,0.016,0.0034,0.0036,0.0009,0.0006,0.015,
&0.0923,0.018,0.022,0.0077,0.009,0.0075,0.024,0.0085,0.067,0.0511,
&0.017,0.0004,0.0028,0.01,2*0.02,0.03,2*0.005,2*0.02,0.03,2*0.005,
&0.015,0.037,0.028,0.079,0.095,0.052,0.0078,4*0.001,0.028,0.033,
&0.026,0.05,0.01,4*0.005,0.25,0.0952,0.02,0.055,2*0.005,0.008,
&0.012,0.02,0.055,2*0.005,0.008,0.012,0.01,0.03,0.0035,0.011,
&0.0055,0.0042,0.009,0.018,0.015,0.0185,0.0135,0.025,0.0004,
&0.0007,0.0008,0.0014,0.0019,0.0025,0.4291,0.08,0.07,0.02,0.015,
&0.005,0.02,0.055,2*0.005,0.008,0.012,0.02,0.055,2*0.005,0.008,
&0.012,0.01,0.03,0.0035,0.011,0.0055,0.0042,0.009,0.018,0.015,
&0.0185,0.0135,0.025,0.0004,0.0007,0.0008,0.0014,0.0019,0.0025,
&0.4291,0.08,0.07,0.02,0.015,0.005,0.02,0.055,2*0.005,0.008,0.012,
&0.02,0.055,2*0.005,0.008,0.012,0.01,0.03,0.0035,0.011,0.0055,
&0.0042,0.009,0.018,0.015,0.0185,0.0135,0.025,2*0.0002,0.0007,
&2*0.0004,0.0014,0.001,0.0009,0.0025,0.4291,0.08,0.07,0.02,0.015,
&0.005,0.047,0.122,0.006,0.012,0.035,0.012,0.035,0.003,0.007,0.15,
&0.037,0.008,0.002,0.05,0.015,0.003,0.001,0.014,0.042,0.014,0.042/
DATA (BRAT(I) ,I= 656, 931)/0.24,0.065,0.012,0.003,0.001,0.002,
&0.001,0.002,0.014,0.003,0.988,0.012,0.389,0.319,0.2367,0.049,
&0.005,0.001,0.0003,0.441,0.206,0.3,0.03,0.022,0.001,5*1.,0.99955,
&0.00045,0.665,0.333,0.002,0.666,0.333,0.001,0.65,0.3,0.05,0.56,
&0.44,5*1.,0.99912,0.00079,0.00005,0.00004,0.888,0.085,0.021,
&2*0.003,0.49,0.344,3*0.043,0.023,0.013,0.001,0.0627,0.0597,
&0.8776,3*0.027,0.015,0.045,0.015,0.045,0.77,0.029,4*1.,0.28,0.14,
&0.313,0.157,0.11,0.28,0.14,0.313,0.157,0.11,0.667,0.333,0.667,
&0.333,2*0.5,0.667,0.333,0.667,0.333,4*0.5,1.,0.333,0.334,0.333,
&4*0.25,6*1.,0.667,0.333,0.667,0.333,0.667,0.333,0.667,0.333,
&2*0.5,0.667,0.333,0.667,0.333,4*0.5,1.,0.52,0.26,0.11,2*0.055,
&0.62,0.31,0.035,2*0.0175,0.007,0.993,0.02,0.98,3*1.,2*0.5,0.667,
&0.333,0.667,0.333,0.667,0.333,0.667,0.333,2*0.5,0.667,0.333,
&0.667,0.333,6*0.5,3*0.12,0.097,0.043,4*0.095,4*0.03,4*0.25,0.273,
&0.727,0.35,0.65,3*1.,2*0.35,0.144,0.105,0.048,0.003,0.333,0.166,
&0.168,0.084,0.087,0.043,0.059,2*0.029,0.002,0.332,0.166,0.168,
&0.084,0.086,0.043,0.059,2*0.029,2*0.002,0.3,0.15,0.16,0.08,0.13,
&0.06,0.08,0.04,0.3,0.15,0.16,0.08,0.13,0.06,0.08,0.04,2*0.3,
&2*0.2,0.3,0.15,0.16,0.08,0.13,0.06,0.08,0.04,0.3,0.15,0.16,0.08,
&0.13,0.06,0.08,0.04,2*0.3,2*0.2,2*0.3,2*0.2,2*0.35,0.144,0.105/
DATA (BRAT(I) ,I= 932,2000)/0.024,2*0.012,0.003,0.566,0.283,
&0.069,0.028,0.023,2*0.0115,0.005,0.003,0.356,2*0.178,0.28,
&2*0.004,0.135,0.865,0.22,0.78,3*1.,0.217,0.124,2*0.193,2*0.135,
&0.002,0.001,0.686,0.314,2*0.0083,0.1866,0.324,0.184,0.027,0.001,
&0.093,0.087,0.078,0.0028,3*0.014,0.008,0.024,0.008,0.024,0.425,
&0.02,0.185,0.088,0.043,0.067,0.066,0.641,0.357,2*0.001,0.018,
&2*0.005,0.003,0.002,2*0.006,0.018,2*0.005,0.003,0.002,2*0.006,
&0.0066,0.025,0.016,0.0088,2*0.005,0.0058,0.005,0.0055,4*0.004,
&2*0.002,2*0.004,0.003,0.002,2*0.003,3*0.002,2*0.001,0.002,
&2*0.001,2*0.002,0.0013,0.0018,5*0.001,4*0.003,2*0.005,2*0.002,
&2*0.001,2*0.002,2*0.001,0.2432,0.057,2*0.035,0.15,2*0.075,0.03,
&2*0.015,2*1.,2*0.105,0.04,0.0077,0.02,0.0235,0.0285,0.0435,
&0.0011,0.0022,0.0044,0.4291,0.08,0.07,0.02,0.015,0.005,2*1.,
&0.999,0.001,1.,0.516,0.483,0.001,1.,0.995,0.005,13*1.,0.331,
&0.663,0.006,0.663,0.331,0.006,1.,0.88,2*0.06,0.88,2*0.06,0.88,
&2*0.06,0.667,2*0.333,0.667,0.676,0.234,0.085,0.005,3*1.,4*0.5,
&7*1.,847*0./
DATA (KFDP(I,1),I= 1, 507)/21,22,23,4*-24,25,21,22,23,4*24,25,
&21,22,23,4*-24,25,21,22,23,4*24,25,21,22,23,4*-24,25,21,22,23,
&4*24,25,37,21,22,23,4*-24,25,2*-37,21,22,23,4*24,25,2*37,22,23,
&-24,25,23,24,-12,22,23,-24,25,23,24,-12,-14,35*16,22,23,-24,25,
&23,24,-89,22,23,-24,25,-37,23,24,37,1,2,3,4,5,6,7,8,21,1,2,3,4,5,
&6,7,8,11,13,15,17,1,2,3,4,5,6,7,8,11,12,13,14,15,16,17,18,4*-1,
&4*-3,4*-5,4*-7,-11,-13,-15,-17,1,2,3,4,5,6,7,8,11,13,15,17,21,
&2*22,23,24,1,2,3,4,5,6,7,8,11,12,13,14,15,16,17,18,24,37,2*23,25,
&35,4*-1,4*-3,4*-5,4*-7,-11,-13,-15,-17,3*24,1,2,3,4,5,6,7,8,11,
&13,15,17,21,2*22,23,24,23,25,36,1,2,3,4,5,6,7,8,11,13,15,17,21,
&2*22,23,24,23,-1,-3,-5,-7,-11,-13,-15,-17,24,5,6,21,2,1,2,3,4,5,
&6,11,13,15,82,-11,-13,2*2,-12,-14,-16,2*-2,2*-4,-2,-4,2*89,37,
&2*-89,2*5,-37,2*89,4*-1,4*-3,4*-5,4*-7,-11,-13,-15,-17,-13,130,
&310,-13,3*211,12,14,11*-11,11*-13,-311,-313,-311,-313,-20313,
&2*-311,-313,-311,-313,2*111,2*221,2*331,2*113,2*223,2*333,-311,
&-313,2*-321,211,-311,-321,333,-311,-313,-321,211,2*-321,2*-311,
&-321,211,113,8*-11,8*-13,-321,-323,-321,-323,-311,2*-313,-311,
&-313,2*-311,-321,-10323,-321,-323,-321,-311,2*-313,211,111,333,
&3*-321,-311,-313,-321,-313,310,333,211,2*-321,-311,-313,-311,211,
&-321,3*-311,211,113,321,-15,5*-11,5*-13,221,331,333,221,331,333/
DATA (KFDP(I,1),I= 508, 924)/10221,211,213,211,213,321,323,321,
&323,2212,221,331,333,221,2*2,6*12,6*14,2*16,3*-411,3*-413,2*-411,
&2*-413,2*441,2*443,2*20443,2*2,2*4,2,4,6*12,6*14,2*16,3*-421,
&3*-423,2*-421,2*-423,2*441,2*443,2*20443,2*2,2*4,2,4,6*12,6*14,
&2*16,3*-431,3*-433,2*-431,2*-433,3*441,3*443,3*20443,2*2,2*4,2,4,
&16,2*4,2*12,2*14,2*16,4*2,4*4,2*-11,2*-13,2*-1,2*-3,2*-11,2*-13,
&2*-1,3*22,111,211,2*22,211,22,211,111,3*22,111,82,21,3*0,2*211,
&321,3*311,2*321,421,2*411,2*421,431,511,521,531,541,211,111,13,
&11,211,22,211,2*111,321,130,-213,113,213,211,22,111,11,13,82,11,
&13,15,1,2,3,4,21,22,3*0,223,321,311,323,313,2*311,321,313,323,
&321,423,2*413,2*423,413,523,2*513,2*523,2*513,523,223,213,113,
&-213,313,-313,323,-323,82,21,3*0,221,321,2*311,321,421,2*411,421,
&411,421,521,2*511,2*521,2*511,521,221,211,111,321,130,310,211,
&111,321,130,310,443,82,553,21,3*0,113,213,323,2*313,323,423,
&2*413,2*423,413,523,2*513,2*523,2*513,523,213,-213,10211,10111,
&-10211,2*221,213,2*113,-213,2*321,2*311,313,-313,323,-323,443,82,
&553,21,3*0,213,113,221,223,321,211,321,311,323,313,323,313,321,
&4*311,321,313,323,313,323,311,4*321,421,411,423,413,423,413,421,
&2*411,421,413,423,413,423,411,2*421,411,423,413,521,511,523,513,
&523,513,521,2*511,521,513,523,513,523,511,2*521,511,523,513,511/
DATA (KFDP(I,1),I= 925,2000)/521,513,523,213,-213,221,223,321,
&130,310,111,211,111,2*211,321,130,310,221,111,321,130,310,221,
&211,111,443,82,553,21,3*0,111,211,-12,12,-14,14,211,111,211,111,
&11,13,82,4*443,10441,20443,445,441,11,13,15,1,2,3,4,21,22,2*553,
&10551,20553,555,2212,2*2112,-12,7*-11,7*-13,2*2224,2*2212,2*2214,
&2*3122,2*3212,2*3214,5*3222,4*3224,2*3322,3324,2*2224,7*2212,
&5*2214,2*2112,2*2114,2*3122,2*3212,2*3214,2*3222,2*3224,4*2,3,
&2*2,1,2*2,2*0,-12,-14,-16,5*4122,441,443,20443,2*-2,2*-4,-2,-4,
&2*0,2112,-12,3122,2212,2112,2212,3*3122,3*4122,4132,4232,0,
&3*5122,5132,5232,0,2112,2212,2*2112,2212,2112,2*2212,3122,3212,
&3112,3122,3222,3112,3122,3222,3212,3322,3312,3322,3312,3122,3322,
&3312,-12,3*4122,2*4132,2*4232,4332,3*5122,5132,5232,5332,847*0/
DATA (KFDP(I,2),I= 1, 476)/3*1,2,4,6,8,1,3*2,1,3,5,7,2,3*3,2,4,
&6,8,3,3*4,1,3,5,7,4,3*5,2,4,6,8,5,3*6,1,3,5,7,6,5,3*7,2,4,6,8,7,
&4,6,3*8,1,3,5,7,8,5,7,2*11,12,11,12,2*11,2*13,14,13,14,13,11,13,
&-211,-213,-211,-213,-211,-213,3*-211,-321,-323,-321,-323,3*-321,
&4*-211,-213,-211,-213,-211,-213,-211,-213,-211,-213,6*-211,2*15,
&16,15,16,15,18,2*17,18,17,2*18,2*17,-1,-2,-3,-4,-5,-6,-7,-8,21,
&-1,-2,-3,-4,-5,-6,-7,-8,-11,-13,-15,-17,-1,-2,-3,-4,-5,-6,-7,-8,
&-11,-12,-13,-14,-15,-16,-17,-18,2,4,6,8,2,4,6,8,2,4,6,8,2,4,6,8,
&12,14,16,18,-1,-2,-3,-4,-5,-6,-7,-8,-11,-13,-15,-17,21,22,2*23,
&-24,-1,-2,-3,-4,-5,-6,-7,-8,-11,-12,-13,-14,-15,-16,-17,-18,-24,
&-37,22,25,2*36,2,4,6,8,2,4,6,8,2,4,6,8,2,4,6,8,12,14,16,18,23,22,
&25,-1,-2,-3,-4,-5,-6,-7,-8,-11,-13,-15,-17,21,22,2*23,-24,2*25,
&36,-1,-2,-3,-4,-5,-6,-7,-8,-11,-13,-15,-17,21,22,2*23,-24,25,2,4,
&6,8,12,14,16,18,25,-5,-6,21,11,-3,-4,-5,-6,-7,-8,-13,-15,-17,-82,
&12,14,-1,-3,11,13,15,1,4,3,4,1,3,5,3,5,6,4,21,22,4,7,5,2,4,6,8,2,
&4,6,8,2,4,6,8,2,4,6,8,12,14,16,18,14,2*0,14,111,211,111,-11,-13,
&11*12,11*14,2*211,2*213,211,20213,2*321,2*323,211,213,211,213,
&211,213,211,213,211,213,211,213,3*211,213,211,2*321,8*211,2*113,
&2*211,8*12,8*14,2*211,2*213,2*111,221,2*113,223,333,20213,211,
&2*321,323,2*311,313,-211,111,113,2*211,321,2*211,311,321,310,211/
DATA (KFDP(I,2),I= 477, 857)/-211,4*211,321,4*211,113,2*211,-321,
&16,5*12,5*14,3*211,3*213,211,2*111,2*113,2*-311,2*-313,-2112,
&3*321,323,2*-1,6*-11,6*-13,2*-15,211,213,20213,211,213,20213,431,
&433,431,433,311,313,311,313,311,313,-1,-4,-3,-4,-1,-3,6*-11,
&6*-13,2*-15,211,213,20213,211,213,20213,431,433,431,433,321,323,
&321,323,321,323,-1,-4,-3,-4,-1,-3,6*-11,6*-13,2*-15,211,213,
&20213,211,213,20213,431,433,431,433,221,331,333,221,331,333,221,
&331,333,-1,-4,-3,-4,-1,-3,-15,-3,-1,2*-11,2*-13,2*-15,-1,-4,-3,
&-4,-3,-4,-1,-4,2*12,2*14,2,3,2,3,2*12,2*14,2,1,22,11,22,111,-211,
&211,11,-211,13,-211,111,113,223,22,111,-82,21,3*0,111,22,-211,
&111,22,211,111,22,211,111,22,111,6*22,-211,22,-13,-11,-211,111,
&-211,2*111,-321,310,211,111,2*-211,221,22,-11,-13,-82,-11,-13,
&-15,-1,-2,-3,-4,2*21,3*0,211,-213,113,-211,111,223,213,113,211,
&111,223,211,111,-211,111,321,311,-211,111,211,111,-321,-311,411,
&421,111,-211,111,211,-311,311,-321,321,-82,21,3*0,211,-211,111,
&211,111,211,111,-211,111,311,321,-211,111,211,111,-321,-311,411,
&421,111,-211,111,-321,130,310,-211,111,-321,130,310,22,-82,22,21,
&3*0,211,111,-211,111,211,111,211,111,-211,111,321,311,-211,111,
&211,111,-321,-311,411,421,-211,211,-211,111,2*211,111,-211,211,
&111,211,-321,2*-311,-321,-311,311,-321,321,22,-82,22,21,3*0,111/
DATA (KFDP(I,2),I= 858,2000)/3*211,-311,22,-211,111,-211,111,
&-211,211,-213,113,223,221,211,111,211,111,2*211,213,113,223,221,
&22,211,111,211,111,4*211,-211,111,-211,111,-211,211,-211,211,321,
&311,321,311,-211,111,-211,111,-211,211,-211,2*211,111,211,111,
&4*211,-321,-311,-321,-311,411,421,411,421,-211,211,111,211,-321,
&130,310,22,-211,111,2*-211,-321,130,310,221,111,-321,130,310,221,
&-211,111,22,-82,22,21,3*0,111,-211,11,-11,13,-13,-211,111,-211,
&111,-11,-13,-82,211,111,221,111,4*22,-11,-13,-15,-1,-2,-3,-4,
&2*21,211,111,3*22,-211,111,22,11,7*12,7*14,-321,-323,-311,-313,
&-311,-313,211,213,211,213,211,213,111,221,331,113,223,111,221,
&113,223,321,323,321,-211,-213,111,221,331,113,223,333,10221,111,
&221,331,113,223,211,213,211,213,321,323,321,323,321,323,311,313,
&311,313,2*-1,-3,-1,2203,3201,3203,2203,2101,2103,2*0,11,13,15,
&-211,-213,-20213,-431,-433,3*3122,1,4,3,4,1,3,2*0,-211,11,22,111,
&211,22,-211,111,22,-211,111,211,2*22,0,-211,111,211,2*22,0,
&2*-211,111,22,111,211,22,211,2*-211,2*111,-211,2*211,111,211,
&-211,2*111,211,-321,-211,111,11,-211,111,211,111,22,111,2*22,
&-211,111,211,3*22,847*0/
DATA (KFDP(I,3),I= 1, 944)/75*0,14,6*0,2*16,2*0,5*111,310,130,
&2*0,2*111,310,130,321,113,211,223,221,2*113,2*211,2*223,2*221,
&2*113,221,113,2*213,-213,195*0,4*3,4*4,1,4,3,2*2,10*81,25*0,-211,
&3*111,-311,-313,-311,-321,-313,-323,111,221,331,113,223,-311,
&-313,-311,-321,-313,-323,111,221,331,113,223,22*0,111,113,2*211,
&-211,-311,211,111,3*211,-211,7*211,-321,-323,-311,-321,-313,-323,
&-211,-213,-321,-323,-311,-321,-313,-323,-211,-213,22*0,111,113,
&-311,2*-211,211,-211,310,-211,2*111,211,2*-211,-321,-211,2*211,
&-211,111,-211,2*211,0,221,331,333,321,311,221,331,333,321,311,
&20*0,3,0,-411,-413,-10413,-10411,-20413,-415,-411,-413,-10413,
&-10411,-20413,-415,-411,-413,16*0,-4,-1,-4,-3,2*-2,-421,-423,
&-10423,-10421,-20423,-425,-421,-423,-10423,-10421,-20423,-425,
&-421,-423,16*0,-4,-1,-4,-3,2*-2,-431,-433,-10433,-10431,-20433,
&-435,-431,-433,-10433,-10431,-20433,-435,-431,-433,19*0,-4,-1,-4,
&-3,2*-2,3*0,441,443,441,443,441,443,-4,-1,-4,-3,-4,-3,-4,-1,531,
&533,531,533,3,2,3,2,511,513,511,513,1,2,0,-11,0,2*111,-211,-11,
&11,-13,2*221,3*0,111,27*0,111,2*0,22,111,5*0,111,12*0,2*21,103*0,
&-211,2*111,-211,3*111,-211,111,211,14*0,111,6*0,111,-211,8*0,111,
&-211,9*0,111,-211,111,-211,4*0,111,-211,111,-211,8*0,111,-211,
&111,-211,4*0,111,-211,111,-211,11*0,-211,6*0,111,211,4*0,111/
DATA (KFDP(I,3),I= 945,2000)/13*0,2*111,211,-211,211,-211,7*0,
&-211,111,13*0,2*21,-211,111,6*0,2212,3122,3212,3214,2112,2114,
&2212,2112,3122,3212,3214,2112,2114,2212,2112,52*0,3*3,1,8*0,
&3*4122,8*0,4,1,4,3,2*2,3*0,2112,43*0,3322,861*0/
DATA (KFDP(I,4),I= 1,2000)/88*0,3*111,8*0,-211,0,-211,3*0,111,
&2*-211,0,111,0,2*111,113,221,111,-213,-211,211,195*0,13*81,41*0,
&111,211,111,211,7*0,111,211,111,211,35*0,2*-211,2*111,211,111,
&-211,2*211,2*-211,2*0,-211,111,-211,111,4*0,-211,111,-211,111,
&34*0,111,-211,3*111,3*-211,2*111,3*-211,4*0,-321,-311,3*0,-321,
&-311,20*0,-3,31*0,6*1,30*0,6*2,33*0,6*3,9*0,8*4,4*0,4*-5,4*0,
&2*-5,7*0,-11,264*0,111,-211,4*0,111,57*0,-211,111,5*0,-211,111,
&52*0,2101,2103,2*2101,19*0,6*2101,909*0/
DATA (KFDP(I,5),I= 1,2000)/90*0,111,16*0,111,7*0,111,0,2*111,
&303*0,-211,2*111,-211,111,-211,111,54*0,111,-211,3*111,-211,111,
&1510*0/
C...LUDAT4, with character strings.
DATA (CHAF(I) ,I= 1, 281)/'d','u','s','c','b','t','l','h',
&2*' ','e','nu_e','mu','nu_mu','tau','nu_tau','chi','nu_chi',
&2*' ','g','gamma','Z','W','H',2*' ','reggeon','pomeron',2*' ',
&'Z''','Z"','W''','H''','A','H','eta_tech','LQ_ue','R',40*' ',
&'specflav','rndmflav','phasespa','c-hadron','b-hadron',
&'t-hadron','l-hadron','h-hadron','Wvirt','diquark','cluster',
&'string','indep.','CMshower','SPHEaxis','THRUaxis','CLUSjet',
&'CELLjet','table',' ','pi',2*'K',2*'D','D_s',2*'B','B_s','B_c',
&'pi','eta','eta''','eta_c','eta_b','eta_t','eta_l','eta_h',2*' ',
&'rho',2*'K*',2*'D*','D*_s',2*'B*','B*_s','B*_c','rho','omega',
&'phi','J/psi','Upsilon','Theta','Theta_l','Theta_h',2*' ','b_1',
&2*'K_1',2*'D_1','D_1s',2*'B_1','B_1s','B_1c','b_1','h_1','h''_1',
&'h_1c','h_1b','h_1t','h_1l','h_1h',2*' ','a_0',2*'K*_0',2*'D*_0',
&'D*_0s',2*'B*_0','B*_0s','B*_0c','a_0','f_0','f''_0','chi_0c',
&'chi_0b','chi_0t','chi_0l','chi_0h',2*' ','a_1',2*'K*_1',
&2*'D*_1','D*_1s',2*'B*_1','B*_1s','B*_1c','a_1','f_1','f''_1',
&'chi_1c','chi_1b','chi_1t','chi_1l','chi_1h',2*' ','a_2',
&2*'K*_2',2*'D*_2','D*_2s',2*'B*_2','B*_2s','B*_2c','a_2','f_2',
&'f''_2','chi_2c','chi_2b','chi_2t','chi_2l','chi_2h',2*' ','K_L',
&'K_S',8*' ','psi''',3*' ','Upsilon''',45*' ','pi_diffr'/
DATA (CHAF(I) ,I= 282, 500)/'n_diffr','p_diffr','rho_diff',
&'omega_di','phi_diff','J/psi_di',18*' ','Lambda',5*' ',
&'Lambda_c',' ',2*'Xi_c',6*' ','Lambda_b',' ',2*'Xi_b',6*' ','n',
&'p',' ',3*'Sigma',2*'Xi',' ',3*'Sigma_c',2*'Xi''_c','Omega_c',
&4*' ',3*'Sigma_b',2*'Xi''_b','Omega_b',4*' ',4*'Delta',
&3*'Sigma*',2*'Xi*','Omega',3*'Sigma*_c',2*'Xi*_c','Omega*_c',
&4*' ',3*'Sigma*_b',2*'Xi*_b','Omega*_b',114*' '/
C...LUDATR, with initial values for the random number generator.
DATA MRLU/19780503,0,0,97,33,0/
END