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