replaces AliPHOSAliEnFile: references to AliEn removed
[u/mrichter/AliRoot.git] / PYTHIA6 / pythia6203.f
1 C*********************************************************************
2 C*********************************************************************
3 C*                                                                  **
4 C*                                                   August 2001    **
5 C*                                                                  **
6 C*                       The Lund Monte Carlo                       **
7 C*                                                                  **
8 C*                        PYTHIA version 6.2                        **
9 C*                                                                  **
10 C*                        Torbjorn Sjostrand                        **
11 C*                 Department of Theoretical Physics                **
12 C*                         Lund University                          **
13 C*               Solvegatan 14A, S-223 62 Lund, Sweden              **
14 C*                    phone +46 - 46 - 222 48 16                    **
15 C*                    E-mail torbjorn@thep.lu.se                    **
16 C*                                                                  **
17 C*                  SUSY and Technicolor parts by                   **
18 C*                         Stephen Mrenna                           **
19 C*                   Physics Department, UC Davis                   **
20 C*             One Shields Avenue, Davis, CA 95616, USA             **
21 C*                   phone + 1 - 530 - 752 - 2661                   **
22 C*                E-mail mrenna@physics.ucdavis.edu                 **
23 C*                                                                  **
24 C*                  PYTHIA 7 efforts coordinated by                 **
25 C*                          Leif Lonnblad                           **
26 C*                 Department of Theoretical Physics                **
27 C*                         Lund University                          **
28 C*               Solvegatan 14A, S-223 62 Lund, Sweden              **
29 C*                    phone +46 - 46 - 222 77 80                    **
30 C*                      E-mail leif@thep.lu.se                      **
31 C*                                                                  **
32 C*         Several parts are written by Hans-Uno Bengtsson          **
33 C*          PYSHOW is written together with Mats Bengtsson          **
34 C*               PYMAEL is written by Emanuel Norrbin               **
35 C*     advanced popcorn baryon production written by Patrik Eden    **
36 C*    code for virtual photons mainly written by Christer Friberg   **
37 C*    code for low-mass strings mainly written by Emanuel Norrbin   **
38 C*        Bose-Einstein code mainly written by Leif Lonnblad        **
39 C*           Lepton number violation code by Peter Skands           **
40 C*      CTEQ  parton distributions are by the CTEQ collaboration    **
41 C*      GRV 94 parton distributions are by Glueck, Reya and Vogt    **
42 C*   SaS photon parton distributions together with Gerhard Schuler  **
43 C*     g + g and q + qbar -> t + tbar + H code by Zoltan Kunszt     **
44 C*         MSSM Higgs mass calculation code by M. Carena,           **
45 C*           J.R. Espinosa, M. Quiros and C.E.M. Wagner             **
46 C*         PYGAUS adapted from CERN library (K.S. Kolbig)           **
47 C*                                                                  **
48 C*   The latest program version and documentation is found on WWW   **
49 C*            http://www.thep.lu.se/~torbjorn/Pythia.html           **
50 C*                                                                  **
51 C*              Copyright Torbjorn Sjostrand, Lund 2001             **
52 C*                                                                  **
53 C*********************************************************************
54 C*********************************************************************
55 C                                                                    *
56 C  List of subprograms in order of appearance, with main purpose     *
57 C  (S = subroutine, F = function, B = block data)                    *
58 C                                                                    *
59 C  B   PYDATA   to contain all default values                        *
60 C  S   PYTEST   to test the proper functioning of the package        *
61 C  S   PYHEPC   to convert between /PYJETS/ and /HEPEVT/ records     *
62 C                                                                    *
63 C  S   PYINIT   to administer the initialization procedure           *
64 C  S   PYEVNT   to administer the generation of an event             *
65 C  S   PYSTAT   to print cross-section and other information         *
66 C  S   PYINRE   to initialize treatment of resonances                *
67 C  S   PYINBM   to read in beam, target and frame choices            *
68 C  S   PYINKI   to initialize kinematics of incoming particles       *
69 C  S   PYINPR   to set up the selection of included processes        *
70 C  S   PYXTOT   to give total, elastic and diffractive cross-sect.   *
71 C  S   PYMAXI   to find differential cross-section maxima            *
72 C  S   PYPILE   to select multiplicity of pileup events              *
73 C  S   PYSAVE   to save alternatives for gamma-p and gamma-gamma     *
74 C  S   PYGAGA   to handle lepton -> lepton + gamma branchings        *
75 C  S   PYRAND   to select subprocess and kinematics for event        *
76 C  S   PYSCAT   to set up kinematics and colour flow of event        *
77 C  S   PYSSPA   to simulate initial state spacelike showers          *
78 C  S   PYMEMX   auxiliary to PYSSPA for ME correction maximum        *
79 C  S   PYMEWT   auxiliary to PYSSPA for matrix element correction    *
80 C  S   PYADSH   to administrate sequential final-state showers       *
81 C  S   PYRESD   to perform resonance decays                          *
82 C  S   PYMULT   to generate multiple interactions                    *
83 C  S   PYREMN   to add on target remnants                            *
84 C  S   PYDIFF   to set up kinematics for diffractive events          *
85 C  S   PYDISG   to set up kinematics, remnant and showers for DIS    *
86 C  S   PYDOCU   to compute cross-sections and handle documentation   *
87 C  S   PYFRAM   to perform boosts between different frames           *
88 C  S   PYWIDT   to calculate full and partial widths of resonances   *
89 C  S   PYOFSH   to calculate partial width into off-shell channels   *
90 C  S   PYRECO   to handle colour reconnection in W+W- events         *
91 C  S   PYKLIM   to calculate borders of allowed kinematical region   *
92 C  S   PYKMAP   to construct value of kinematical variable           *
93 C  S   PYSIGH   to calculate differential cross-sections             *
94 C  S   PYPDFU   to evaluate parton distributions                     *
95 C  S   PYPDFL   to evaluate parton distributions at low x and Q^2    *
96 C  S   PYPDEL   to evaluate electron parton distributions            *
97 C  S   PYPDGA   to evaluate photon parton distributions (generic)    *
98 C  S   PYGGAM   to evaluate photon parton distributions (SaS sets)   *
99 C  S   PYGVMD   to evaluate VMD part of photon parton distributions  *
100 C  S   PYGANO   to evaluate anomalous part of photon pdf's           *
101 C  S   PYGBEH   to evaluate Bethe-Heitler part of photon pdf's       *
102 C  S   PYGDIR   to evaluate direct contribution to photon pdf's      *
103 C  S   PYPDPI   to evaluate pion parton distributions                *
104 C  S   PYPDPR   to evaluate proton parton distributions              *
105 C  F   PYCTEQ   to evaluate the CTEQ 3 proton parton distributions   *
106 C  S   PYGRVL   to evaluate the GRV 94L proton parton distributions  *
107 C  S   PYGRVM   to evaluate the GRV 94M proton parton distributions  *
108 C  S   PYGRVD   to evaluate the GRV 94D proton parton distributions  *
109 C  F   PYGRVV   auxiliary to the PYGRV* routines                     *
110 C  F   PYGRVW   auxiliary to the PYGRV* routines                     *
111 C  F   PYGRVS   auxiliary to the PYGRV* routines                     *
112 C  F   PYCT5L   to evaluate the CTEQ 5L proton parton distributions  *
113 C  F   PYCT5M   to evaluate the CTEQ 5M1 proton parton distributions *
114 C  S   PYPDPO   to evaluate old proton parton distributions          *
115 C  F   PYHFTH   to evaluate threshold factor for heavy flavour       *
116 C  S   PYSPLI   to find flavours left in hadron when one removed     *
117 C  F   PYGAMM   to evaluate ordinary Gamma function Gamma(x)         *
118 C  S   PYWAUX   to evaluate auxiliary functions W1(s) and W2(s)      *
119 C  S   PYI3AU   to evaluate auxiliary function I3(s,t,u,v)           *
120 C  F   PYSPEN   to evaluate Spence (dilogarithm) function Sp(x)      *
121 C  S   PYQQBH   to evaluate matrix element for g + g -> Q + Qbar + H *
122 C                                                                    *
123 C  S   PYMSIN   to initialize the supersymmetry simulation           *
124 C  S   PYAPPS   to determine MSSM parameters from SUGRA input        *
125 C  F   PYRNMQ   to determine running quark masses                    *
126 C  F   PYRNMT   to determine running top mass                        *
127 C  S   PYTHRG   to calculate sfermion third-gen. mass eigenstates    *
128 C  S   PYINOM   to calculate neutralino/chargino mass eigenstates    *
129 C  F   PYRNM3   to determine running M3, gluino mass                 *
130 C  S   PYEIG4   to calculate eigenvalues and -vectors in 4*4 matrix  *
131 C  S   PYHGGM   to determine Higgs mass spectrum                     *
132 C  S   PYSUBH   to determine Higgs masses in the MSSM                *
133 C  S   PYPOLE   to determine Higgs masses in the MSSM                *
134 C  S   PYRGHM   auxiliary to PYPOLE                                  *
135 C  S   PYGFXX   auxiliary to PYRGHM                                  *
136 C  F   PYFINT   auxiliary to PYPOLE                                  *
137 C  F   PYFISB   auxiliary to PYFINT                                  *
138 C  S   PYSFDC   to calculate sfermion decay partial widths           *
139 C  S   PYGLUI   to calculate gluino decay partial widths             *
140 C  S   PYTBBN   to calculate 3-body decay of gluino to neutralino    *
141 C  S   PYTBBC   to calculate 3-body decay of gluino to chargino      *
142 C  S   PYNJDC   to calculate neutralino decay partial widths         *
143 C  S   PYCJDC   to calculate chargino decay partial widths           *
144 C  F   PYXXZ6   auxiliary for ino 3-body decays                      *
145 C  F   PYXXGA   auxiliary for ino -> ino + gamma decay               *
146 C  F   PYX2XG   auxiliary for ino -> ino + gauge boson decay         *
147 C  F   PYX2XH   auxiliary for ino -> ino + Higgs decay               *
148 C  S   PYHEXT   to calculate non-SM Higgs decay partial widths       *
149 C  F   PYH2XX   auxiliary for H -> ino + ino decay                   *
150 C  F   PYGAUS   to perform Gaussian integration                      *
151 C  F   PYSIMP   to perform Simpson integration                       *
152 C  F   PYLAMF   to evaluate the lambda kinematics function           *
153 C  S   PYTBDY   to perform 3-body decay of gauginos                  *
154 C  S   PYTECM   to calculate techni_rho/omega masses                 *
155 C  S   PYEICG   to calculate eigenvalues of a 4*4 complex matrix     *
156 C  S   PYCMQR   auxiliary to PYEICG                                  *
157 C  S   PYCMQ2   auxiliary to PYEICG                                  *
158 C  S   PYCDIV   auxiliary to PYCMQR                                  *
159 C  S   PYCSRT   auxiliary to PYCMQR                                  *
160 C  S   PYTHAG   auxiliary to PYCMQR                                  *
161 C  S   PYCBAL   auxiliary to PYEICG                                  *
162 C  S   PYCBA2   auxiliary to PYEICG                                  *
163 C  S   PYCRTH   auxiliary to PYEICG                                  *
164 C  S   PYLDCM   auxiliary to PYSIGH, for technicolor in QCD 2 -> 2   *
165 C  S   PYBKSB   auxiliary to PYSIGH, for technicolor in QCD 2 -> 2   *
166 C  S   PYWIDX   to calculate decay widths from within PYWIDT         *
167 C  S   PYRVSF   to calculate R-violating sfermion decay widths       *
168 C  S   PYRVNE   to calculate R-violating neutralino decay widths     *
169 C  S   PYRVCH   to calculate R-violating chargino decay widths       *
170 C  F   PYRVSB   auxiliary to PYRVSF                                  *
171 C  S   PYRVGW   to calculate R-Violating 3-body widths               *
172 C  F   PYRVI1   auxiliary to PYRVGW, to do PS integration for res.   *
173 C  F   PYRVI2   auxiliary to PYRVGW, to do PS integration for LR-int.*
174 C  F   PYRVI3   auxiliary to PYRVGW, to do PS X integral for int.    *
175 C  F   PYRVG1   auxiliary to PYRVI1, general matrix element, res.    *
176 C  F   PYRVG2   auxiliary to PYRVI2, general matrix element, LR-int. *
177 C  F   PYRVG3   auxiliary to PYRVI3, to do PS Y integral for int.    *
178 C  F   PYRVG4   auxiliary to PYRVG3, general matrix element, int.    *
179 C  F   PYRVR    auxiliary to PYRVG1, Breit-Wigner                    *
180 C  F   PYRVS    auxiliary to PYRVG2 & PYRVG4                         *
181 C                                                                    *
182 C  S   PY1ENT   to fill one entry (= parton or particle)             *
183 C  S   PY2ENT   to fill two entries                                  *
184 C  S   PY3ENT   to fill three entries                                *
185 C  S   PY4ENT   to fill four entries                                 *
186 C  S   PY2FRM   to interface to generic two-fermion generator        *
187 C  S   PY4FRM   to interface to generic four-fermion generator       *
188 C  S   PY6FRM   to interface to generic six-fermion generator        *
189 C  S   PY4JET   to generate a shower from a given 4-parton config    *
190 C  S   PY4JTW   to evaluate the weight od a shower history for above *
191 C  S   PY4JTS   to set up the parton configuration for above         *
192 C  S   PYJOIN   to connect entries with colour flow information      *
193 C  S   PYGIVE   to fill (or query) commonblock variables             *
194 C  S   PYEXEC   to administrate fragmentation and decay chain        *
195 C  S   PYPREP   to rearrange showered partons along strings          *
196 C  S   PYSTRF   to do string fragmentation of jet system             *
197 C  S   PYINDF   to do independent fragmentation of one or many jets  *
198 C  S   PYDECY   to do the decay of a particle                        *
199 C  S   PYDCYK   to select parton and hadron flavours in decays       *
200 C  S   PYKFDI   to select parton and hadron flavours in fragm        *
201 C  S   PYNMES   to select number of popcorn mesons                   *
202 C  S   PYKFIN   to calculate falvour prod. ratios from input params. *
203 C  S   PYPTDI   to select transverse momenta in fragm                *
204 C  S   PYZDIS   to select longitudinal scaling variable in fragm     *
205 C  S   PYSHOW   to do timelike parton shower evolution               *
206 C  F   PYMAEL   auxiliary to PYSHOW, with gluon emission ME's        *
207 C  S   PYBOEI   to include Bose-Einstein effects (crudely)           *
208 C  S   PYBESQ   auxiliary to PYBOEI                                  *
209 C  F   PYMASS   to give the mass of a particle or parton             *
210 C  F   PYMRUN   to give the running MSbar mass of a quark            *
211 C  S   PYNAME   to give the name of a particle or parton             *
212 C  F   PYCHGE   to give three times the electric charge              *
213 C  F   PYCOMP   to compress standard KF flavour code to internal KC  *
214 C  S   PYERRM   to write error messages and abort faulty run         *
215 C  F   PYALEM   to give the alpha_electromagnetic value              *
216 C  F   PYALPS   to give the alpha_strong value                       *
217 C  F   PYANGL   to give the angle from known x and y components      *
218 C  F   PYR      to provide a random number generator                 *
219 C  S   PYRGET   to save the state of the random number generator     *
220 C  S   PYRSET   to set the state of the random number generator      *
221 C  S   PYROBO   to rotate and/or boost an event                      *
222 C  S   PYEDIT   to remove unwanted entries from record               *
223 C  S   PYLIST   to list event record or particle data                *
224 C  S   PYLOGO   to write a logo                                      *
225 C  S   PYUPDA   to update particle data                              *
226 C  F   PYK      to provide integer-valued event information          *
227 C  F   PYP      to provide real-valued event information             *
228 C  S   PYSPHE   to perform sphericity analysis                       *
229 C  S   PYTHRU   to perform thrust analysis                           *
230 C  S   PYCLUS   to perform three-dimensional cluster analysis        *
231 C  S   PYCELL   to perform cluster analysis in (eta, phi, E_T)       *
232 C  S   PYJMAS   to give high and low jet mass of event               *
233 C  S   PYFOWO   to give Fox-Wolfram moments                          *
234 C  S   PYTABU   to analyze events, with tabular output               *
235 C                                                                    *
236 C  S   PYEEVT   to administrate the generation of an e+e- event      *
237 C  S   PYXTEE   to give the total cross-section at given CM energy   *
238 C  S   PYRADK   to generate initial state photon radiation           *
239 C  S   PYXKFL   to select flavour of primary qqbar pair              *
240 C  S   PYXJET   to select (matrix element) jet multiplicity          *
241 C  S   PYX3JT   to select kinematics of three-jet event              *
242 C  S   PYX4JT   to select kinematics of four-jet event               *
243 C  S   PYXDIF   to select angular orientation of event               *
244 C  S   PYONIA   to perform generation of onium decay to gluons       *
245 C                                                                    *
246 C  S   PYBOOK   to book a histogram                                  *
247 C  S   PYFILL   to fill an entry in a histogram                      *
248 C  S   PYFACT   to multiply histogram contents by a factor           *
249 C  S   PYOPER   to perform operations between histograms             *
250 C  S   PYHIST   to print and reset all histograms                    *
251 C  S   PYPLOT   to print a single histogram                          *
252 C  S   PYNULL   to reset contents of a single histogram              *
253 C  S   PYDUMP   to dump histogram contents onto a file               *
254 C                                                                    *
255 C  S   PYKCUT   dummy routine for user kinematical cuts              *
256 C  S   PYEVWT   dummy routine for weighting events                   *
257 C  S   UPINIT   dummy routine to initialize user processes           *
258 C  S   UPEVNT   dummy routine to generate a user process event       *
259 C  S   PDFSET   dummy routine to be removed when using PDFLIB        *
260 C  S   STRUCTM  dummy routine to be removed when using PDFLIB        *
261 C  S   STRUCTP  dummy routine to be removed when using PDFLIB        *
262 C  S   PYTAUD   dummy routine for interface to tau decay libraries   *
263 C  S   PYTIME   dummy routine for giving date and time               *
264 C                                                                    *
265 C*********************************************************************
266  
267 C...PYDATA
268 C...Default values for switches and parameters,
269 C...and particle, decay and process data.
270  
271       BLOCK DATA PYDATA
272  
273 C...Double precision and integer declarations.
274       IMPLICIT DOUBLE PRECISION(A-H, O-Z)
275       IMPLICIT INTEGER(I-N)
276 C PH      INTEGER PYK,PYCHGE,PYCOMP
277 C...Commonblocks.
278       COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
279       COMMON/PYDAT2/KCHG(500,4),PMAS(500,4),PARF(2000),VCKM(4,4)
280       COMMON/PYDAT3/MDCY(500,3),MDME(8000,2),BRAT(8000),KFDP(8000,5)
281       COMMON/PYDAT4/CHAF(500,2)
282       CHARACTER CHAF*16
283       COMMON/PYDATR/MRPY(6),RRPY(100)
284       COMMON/PYSUBS/MSEL,MSELPD,MSUB(500),KFIN(2,-40:40),CKIN(200)
285       COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
286       COMMON/PYINT1/MINT(400),VINT(400)
287       COMMON/PYINT2/ISET(500),KFPR(500,2),COEF(500,20),ICOL(40,4,2)
288       COMMON/PYINT3/XSFX(2,-40:40),ISIG(1000,3),SIGH(1000)
289       COMMON/PYINT4/MWID(500),WIDS(500,5)
290       COMMON/PYINT5/NGENPD,NGEN(0:500,3),XSEC(0:500,3)
291       COMMON/PYINT6/PROC(0:500)
292       CHARACTER PROC*28
293       COMMON/PYINT7/SIGT(0:6,0:6,0:5)
294       COMMON/PYMSSM/IMSS(0:99),RMSS(0:99)
295       COMMON/PYSSMT/ZMIX(4,4),UMIX(2,2),VMIX(2,2),SMZ(4),SMW(2),
296      &SFMIX(16,4),ZMIXI(4,4),UMIXI(2,2),VMIXI(2,2)
297       COMMON/PYMSRV/RVLAM(3,3,3), RVLAMP(3,3,3), RVLAMB(3,3,3)
298       COMMON/PYBINS/IHIST(4),INDX(1000),BIN(20000)
299       SAVE /PYDAT1/,/PYDAT2/,/PYDAT3/,/PYDAT4/,/PYDATR/,/PYSUBS/,
300      &/PYPARS/,/PYINT1/,/PYINT2/,/PYINT3/,/PYINT4/,/PYINT5/,
301      &/PYINT6/,/PYINT7/,/PYMSSM/,/PYSSMT/,/PYMSRV/,/PYBINS/
302  
303 C...PYDAT1, containing status codes and most parameters.
304       DATA MSTU/
305      &   0,    0,    0, 4000,10000,  500, 8000,    0,    0,    2,
306      1   6,    1,    1,    0,    0,    1,    0,    0,    0,    0,
307      2   2,   10,    0,    0,    1,   10,    0,    0,    0,    0,
308      3   0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
309      4   2,    2,    1,    4,    2,    1,    1,    0,    0,    0,
310      5  25,   24,    0,    1,    0,    0,    0,    0,    0,    0,
311      6   0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
312      7  30*0,
313      1   1,    0,    0,    0,    0,    0,    0,    0,    0,    0,
314      2   1,    5,    3,    5,    0,    0,    0,    0,    0,    0,
315      &  80*0/
316       DATA (PARU(I),I=1,100)/
317      &  3.141592653589793D0, 6.283185307179586D0,
318      &  0.197327D0, 5.06773D0, 0.389380D0, 2.56819D0,  4*0D0,
319      1  0.001D0, 0.09D0, 0.01D0, 0D0, 0D0, 0D0, 0D0, 0D0, 0D0, 0D0,
320      2  0D0,   0D0,   0D0,   0D0,  0D0,  0D0,  0D0,  0D0,  0D0,  0D0,
321      3  0D0,   0D0,   0D0,   0D0,  0D0,  0D0,  0D0,  0D0,  0D0,  0D0,
322      4  2.0D0,  1.0D0, 0.25D0,  2.5D0, 0.05D0,
323      4  0D0,   0D0, 0.0001D0, 0D0,   0D0,
324      5  2.5D0,1.5D0,7.0D0,1.0D0,0.5D0,2.0D0,3.2D0, 0D0, 0D0, 0D0,
325      6  40*0D0/
326       DATA (PARU(I),I=101,200)/
327      &  0.00729735D0, 0.232D0, 0.007764D0, 1.0D0, 1.16639D-5,
328      &  0D0, 0D0, 0D0, 0D0,  0D0,
329      1  0.20D0, 0.25D0, 1.0D0, 4.0D0, 10D0, 0D0, 0D0,  0D0, 0D0, 0D0,
330      2 -0.693D0, -1.0D0, 0.387D0, 1.0D0, -0.08D0,
331      2 -1.0D0,  1.0D0,  1.0D0,  1.0D0,  0D0,
332      3  1.0D0,-1.0D0, 1.0D0,-1.0D0, 1.0D0,  0D0,  0D0, 0D0, 0D0, 0D0,
333      4  5.0D0, 1.0D0, 1.0D0,  0D0, 1.0D0, 1.0D0,  0D0, 0D0, 0D0, 0D0,
334      5  1.0D0, 0D0, 0D0, 0D0, 1000D0, 1.0D0, 1.0D0, 1.0D0, 1.0D0,0D0,
335      6  1.0D0, 1.0D0, 1.0D0, 1.0D0, 1.0D0,  0D0,  0D0, 0D0, 0D0, 0D0,
336      7  1.0D0, 1.0D0, 1.0D0, 1.0D0, 1.0D0, 1.0D0, 1.0D0, 0D0,0D0,0D0,
337      8  1.0D0, 1.0D0, 1.0D0, 0.0D0, 0.0D0, 1.0D0, 1.0D0, 0D0,0D0,0D0,
338      9  0D0,  0D0,  0D0,  0D0, 1.0D0,  0D0,  0D0, 0D0, 0D0, 0D0/
339       DATA MSTJ/
340      &  1,    3,    0,    0,    0,    0,    0,    0,    0,    0,
341      1  4,    2,    0,    1,    0,    2,    2,    0,    0,    0,
342      2  2,    1,    1,    2,    1,    2,    2,    0,    0,    0,
343      3  0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
344      4  2,    2,    4,    2,    5,    3,    3,    0,    0,    3,
345      5  0,    3,    0,    2,    0,    0,    1,    0,    0,    0,
346      6  40*0,
347      &  5,    2,    7,    5,    1,    1,    0,    2,    0,    2,
348      1  0,    0,    0,    0,    1,    1,    0,    0,    0,    0,
349      2  80*0/
350       DATA PARJ/
351      &  0.10D0, 0.30D0, 0.40D0, 0.05D0, 0.50D0,
352      &  0.50D0, 0.50D0,   0.6D0,   1.2D0,   0.6D0,
353      1  0.50D0,0.60D0,0.75D0, 0D0, 0D0, 0D0, 0D0, 1.0D0, 1.0D0, 0D0,
354      2  0.36D0, 1.0D0,0.01D0, 2.0D0,1.0D0,0.4D0, 0D0, 0D0, 0D0, 0D0,
355      3  0.10D0, 1.0D0, 0.8D0, 1.5D0,0D0,2.0D0,0.2D0, 0D0,0.08D0,0D0,
356      4  0.3D0, 0.58D0, 0.5D0, 0.9D0,0.5D0,1.0D0,1.0D0, 0D0, 0D0,0D0,
357      5  0.77D0, 0.77D0, 0.77D0, -0.05D0, -0.005D0,
358      5  0D0, 0D0, 0D0, 1.0D0, 0D0,
359      6  4.5D0, 0.7D0, 0D0,0.003D0, 0.5D0, 0.5D0, 0D0, 0D0, 0D0, 0D0,
360      7  10D0, 1000D0, 100D0, 1000D0, 0D0, 0.7D0,10D0, 0D0,0D0,0.5D0,
361      8  0.29D0, 1.0D0, 1.0D0,  0D0,  10D0, 10D0, 0D0, 0D0, 0D0,1D-4,
362      9  0.02D0, 1.0D0, 0.2D0,  0D0,  0D0,  0D0,  0D0, 0D0, 0D0, 0D0,
363      &  0D0,  0D0,  0D0,  0D0,   0D0,   0D0,  0D0,  0D0,  0D0,  0D0,
364      1  0D0,  0D0,  0D0,  0D0,   0D0,   0D0,  0D0,  0D0,  0D0,  0D0,
365      2  1.0D0, 0.25D0,91.187D0,2.489D0, 0.01D0,
366      2  2.0D0,  1.0D0, 0.25D0,0.002D0,   0D0,
367      3  0D0, 0D0, 0D0, 0D0, 0.01D0, 0.99D0, 0D0, 0D0,  0.2D0,   0D0,
368      4  10*0D0,
369      5  10*0D0,
370      6  10*0D0,
371      7  0D0, 0D0, 0D0, 0D0, 0D0, 0D0, 0D0, 0D0, 0D0, -0.693D0,
372      8 -1.0D0, 0.387D0, 1.0D0, -0.08D0, -1.0D0,
373      8  1.0D0,  1.0D0, -0.693D0, -1.0D0, 0.387D0,
374      9  1.0D0, -0.08D0, -1.0D0,   1.0D0, 1.0D0,
375      9  5*0D0/
376  
377 C...PYDAT2, with particle data and flavour treatment parameters.
378       DATA (KCHG(I,1),I=   1, 500)/-1,2,-1,2,-1,2,-1,2,2*0,-3,0,-3,0,
379      &-3,0,-3,6*0,3,9*0,3,2*0,3,4*0,-1,41*0,2,-1,20*0,3*3,7*0,3*3,3*0,
380      &3*3,3*0,3*3,6*0,3*3,3*0,3*3,4*0,-2,-3,2*1,2*0,4,2*3,6,2*-2,2*-3,
381      &0,2*1,2*0,2*3,-2,2*-3,2*0,-3,2*1,2*0,3,0,2*4,2*3,2*6,3,2*1,2*0,
382      &2*3,2*0,4,2*3,2*6,2*3,6,2*-2,2*-3,0,-3,0,2*1,2*0,2*3,0,3,2*-2,
383      &2*-3,2*0,2*-3,0,2*1,2*0,2*3,2*0,2*3,-2,2*-3,2*0,2*-3,2*0,-3,2*0,
384      &2*3,4*0,2*3,2*0,2*3,2*0,2*3,4*0,2*3,2*0,2*3,3*0,3,2*0,3,0,3,0,3,
385      &2*0,3,0,3,3*0,-1,2,-1,2,-1,2,-3,0,-3,0,-3,4*0,3,2*0,3,0,-1,2,-1,
386      &2,-1,2,-3,0,-3,0,-3,2*0,3,3*0,3,8*0,-1,2,-3,6*0,3,2*6,0,3,4*0,3,
387      &139*0/
388       DATA (KCHG(I,2),I=   1, 500)/8*1,12*0,2,20*0,1,107*0,-1,0,2*-1,
389      &2*0,-1,3*0,2*-1,3*0,2*-1,4*0,-1,5*0,2*-1,4*0,2*-1,5*0,2*-1,6*0,
390      &-1,7*0,2*-1,5*0,2*-1,6*0,2*-1,7*0,2*-1,8*0,-1,56*0,6*1,6*0,2,7*0,
391      &6*1,9*0,2,3*0,2,0,5*2,2*1,156*0/
392       DATA (KCHG(I,3),I=   1, 500)/8*1,2*0,8*1,5*0,1,9*0,1,2*0,1,3*0,
393      &2*1,39*0,1,0,2*1,20*0,3*1,4*0,6*1,3*0,9*1,3*0,12*1,4*0,100*1,2*0,
394      &2*1,2*0,4*1,2*0,6*1,2*0,8*1,3*0,1,0,2*1,0,3*1,0,4*1,3*0,12*1,3*0,
395      &1,2*0,1,0,12*1,0,1,3*0,1,8*0,4*1,5*0,3*1,0,1,3*0,2*1,139*0/
396       DATA (KCHG(I,4),I=   1, 290)/1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,
397      &16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,
398      &37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,
399      &58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,
400      &79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,
401      &100,110,111,113,115,130,211,213,215,221,223,225,310,311,313,315,
402      &321,323,325,331,333,335,411,413,415,421,423,425,431,433,435,441,
403      &443,445,511,513,515,521,523,525,531,533,535,541,543,545,551,553,
404      &555,990,1103,1114,2101,2103,2112,2114,2203,2212,2214,2224,3101,
405      &3103,3112,3114,3122,3201,3203,3212,3214,3222,3224,3303,3312,3314,
406      &3322,3324,3334,4101,4103,4112,4114,4122,4132,4201,4203,4212,4214,
407      &4222,4224,4232,4301,4303,4312,4314,4322,4324,4332,4334,4403,4412,
408      &4414,4422,4424,4432,4434,4444,5101,5103,5112,5114,5122,5132,5142,
409      &5201,5203,5212,5214,5222,5224,5232,5242,5301,5303,5312,5314,5322,
410      &5324,5332,5334,5342,5401,5403,5412,5414,5422,5424,5432,5434,5442,
411      &5444,5503,5512,5514,5522,5524,5532,5534,5542,5544,5554,10111,
412      &10113,10211,10213,10221,10223,10311,10313,10321,10323,10331,
413      &10333,10411,10413,10421,10423,10431,10433,10441,10443,10511,
414      &10513,10521,10523,10531,10533,10541,10543,10551,10553,20113,
415      &20213,20223,20313,20323,20333,20413,20423,20433,20443,20513/
416       DATA (KCHG(I,4),I= 291, 500)/20523,20533,20543,20553,100443,
417      &100553,1000001,1000002,1000003,1000004,1000005,1000006,1000011,
418      &1000012,1000013,1000014,1000015,1000016,1000021,1000022,1000023,
419      &1000024,1000025,1000035,1000037,1000039,2000001,2000002,2000003,
420      &2000004,2000005,2000006,2000011,2000012,2000013,2000014,2000015,
421      &2000016,3000111,3000211,3000221,3000331,3000113,3000213,3000223,
422      &3100021,3100111,3200111,3100113,3200113,3300113,3400113,4000001,
423      &4000002,4000011,4000012,5000039,9900012,9900014,9900016,9900023,
424      &9900024,9900041,9900042,9900110,9900210,9900220,9900330,9900440,
425      &9902110,9902210,139*0/
426       DATA (PMAS(I,1),I=   1, 217)/2*0.33D0,0.5D0,1.5D0,4.8D0,175D0,
427      &2*400D0,2*0D0,0.00051D0,0D0,0.10566D0,0D0,1.777D0,0D0,400D0,
428      &5*0D0,91.188D0,80.45D0,115D0,6*0D0,500D0,900D0,500D0,3*300D0,
429      &3*0D0,5000D0,200D0,40*0D0,1D0,2D0,5D0,16*0D0,0.13498D0,0.7685D0,
430      &1.318D0,0.49767D0,0.13957D0,0.7669D0,1.318D0,0.54745D0,0.78194D0,
431      &1.275D0,2*0.49767D0,0.8961D0,1.432D0,0.4936D0,0.8916D0,1.425D0,
432      &0.95777D0,1.0194D0,1.525D0,1.8693D0,2.01D0,2.46D0,1.8645D0,
433      &2.0067D0,2.46D0,1.9685D0,2.1124D0,2.5735D0,2.9798D0,3.09688D0,
434      &3.5562D0,5.2792D0,5.3248D0,5.83D0,5.2789D0,5.3248D0,5.83D0,
435      &5.3693D0,5.4163D0,6.07D0,6.594D0,6.602D0,7.35D0,9.4D0,9.4603D0,
436      &9.9132D0,0D0,0.77133D0,1.234D0,0.57933D0,0.77133D0,0.93957D0,
437      &1.233D0,0.77133D0,0.93827D0,1.232D0,1.231D0,0.80473D0,0.92953D0,
438      &1.19744D0,1.3872D0,1.11568D0,0.80473D0,0.92953D0,1.19255D0,
439      &1.3837D0,1.18937D0,1.3828D0,1.09361D0,1.3213D0,1.535D0,1.3149D0,
440      &1.5318D0,1.67245D0,1.96908D0,2.00808D0,2.4521D0,2.5D0,2.2849D0,
441      &2.4703D0,1.96908D0,2.00808D0,2.4535D0,2.5D0,2.4529D0,2.5D0,
442      &2.4656D0,2.15432D0,2.17967D0,2.55D0,2.63D0,2.55D0,2.63D0,2.704D0,
443      &2.8D0,3.27531D0,3.59798D0,3.65648D0,3.59798D0,3.65648D0,
444      &3.78663D0,3.82466D0,4.91594D0,5.38897D0,5.40145D0,5.8D0,5.81D0,
445      &5.641D0,5.84D0,7.00575D0,5.38897D0,5.40145D0,5.8D0,5.81D0,5.8D0/
446       DATA (PMAS(I,1),I= 218, 500)/5.81D0,5.84D0,7.00575D0,5.56725D0,
447      &5.57536D0,5.96D0,5.97D0,5.96D0,5.97D0,6.12D0,6.13D0,7.19099D0,
448      &6.67143D0,6.67397D0,7.03724D0,7.0485D0,7.03724D0,7.0485D0,
449      &7.21101D0,7.219D0,8.30945D0,8.31325D0,10.07354D0,10.42272D0,
450      &10.44144D0,10.42272D0,10.44144D0,10.60209D0,10.61426D0,
451      &11.70767D0,11.71147D0,15.11061D0,0.9835D0,1.231D0,0.9835D0,
452      &1.231D0,1D0,1.17D0,1.429D0,1.29D0,1.429D0,1.29D0,2*1.4D0,2.272D0,
453      &2.424D0,2.272D0,2.424D0,2.5D0,2.536D0,3.4151D0,3.46D0,5.68D0,
454      &5.73D0,5.68D0,5.73D0,5.92D0,5.97D0,7.25D0,7.3D0,9.8598D0,9.875D0,
455      &2*1.23D0,1.282D0,2*1.402D0,1.427D0,2*2.372D0,2.56D0,3.5106D0,
456      &2*5.78D0,6.02D0,7.3D0,9.8919D0,3.686D0,10.0233D0,32*500D0,
457      &3*110D0,350D0,3*210D0,500D0,125D0,250D0,400D0,2*350D0,300D0,
458      &4*400D0,1000D0,3*500D0,1200D0,750D0,2*200D0,7*0D0,139*0D0/
459       DATA (PMAS(I,2),I=   1, 500)/5*0D0,1.39816D0,16*0D0,2.47813D0,
460      &2.07115D0,0.00367D0,6*0D0,14.54029D0,0D0,16.66099D0,8.38842D0,
461      &3.3752D0,4.17669D0,3*0D0,417.29147D0,0.39162D0,60*0D0,0.151D0,
462      &0.107D0,2*0D0,0.149D0,0.107D0,0D0,0.00843D0,0.185D0,2*0D0,
463      &0.0505D0,0.109D0,0D0,0.0498D0,0.098D0,0.0002D0,0.00443D0,0.076D0,
464      &2*0D0,0.023D0,2*0D0,0.023D0,2*0D0,0.015D0,0.0013D0,0D0,0.002D0,
465      &2*0D0,0.02D0,2*0D0,0.02D0,2*0D0,0.02D0,2*0D0,0.02D0,5*0D0,0.12D0,
466      &3*0D0,0.12D0,2*0D0,2*0.12D0,3*0D0,0.0394D0,4*0D0,0.036D0,0D0,
467      &0.0358D0,2*0D0,0.0099D0,0D0,0.0091D0,74*0D0,0.06D0,0.142D0,
468      &0.06D0,0.142D0,0D0,0.36D0,0.287D0,0.09D0,0.287D0,0.09D0,0.25D0,
469      &0.08D0,0.05D0,0.02D0,0.05D0,0.02D0,0.05D0,0D0,0.014D0,0.01D0,
470      &8*0.05D0,0D0,0.01D0,2*0.4D0,0.025D0,2*0.174D0,0.053D0,3*0.05D0,
471      &0.0009D0,4*0.05D0,3*0D0,19*1D0,0D0,7*1D0,0D0,1D0,0D0,1D0,0D0,
472      &0.02911D0,0.01741D0,0.04536D0,0.09511D0,0.8686D0,0.62395D0,
473      &0.19192D0,123.27638D0,0.02296D0,0.18886D0,23.26819D0,2.86306D0,
474      &0D0,3.45903D0,2.59359D0,2.59687D0,0.42896D0,0.41912D0,0.14153D0,
475      &2*0.00098D0,0.00097D0,26.7245D0,21.74916D0,0.88159D0,0.88001D0,
476      &7*0D0,139*0D0/
477       DATA (PMAS(I,3),I=   1, 500)/5*0D0,13.98156D0,16*0D0,24.78129D0,
478      &20.71149D0,0.03669D0,6*0D0,145.40294D0,0D0,166.60993D0,
479      &83.88423D0,33.75195D0,41.76694D0,3*0D0,4172.91467D0,3.91621D0,
480      &60*0D0,0.4D0,0.25D0,2*0D0,0.4D0,0.25D0,0D0,0.1D0,0.17D0,2*0D0,
481      &0.2D0,0.12D0,0D0,0.2D0,0.12D0,0.002D0,0.015D0,0.2D0,2*0D0,0.12D0,
482      &2*0D0,0.12D0,2*0D0,0.05D0,0.005D0,0D0,0.01D0,2*0D0,0.05D0,2*0D0,
483      &0.05D0,2*0D0,0.05D0,2*0D0,0.05D0,5*0D0,0.14D0,3*0D0,0.14D0,2*0D0,
484      &2*0.14D0,3*0D0,0.04D0,4*0D0,0.035D0,0D0,0.035D0,2*0D0,0.05D0,0D0,
485      &0.05D0,74*0D0,0.05D0,0.25D0,0.05D0,0.25D0,0D0,0.2D0,0.4D0,
486      &0.005D0,0.4D0,0.01D0,0.35D0,0.001D0,0.1D0,0.08D0,0.1D0,0.08D0,
487      &0.1D0,0D0,0.05D0,0.02D0,6*0.1D0,0.05D0,0.1D0,0D0,0.02D0,2*0.3D0,
488      &0.05D0,2*0.3D0,0.02D0,2*0.1D0,0.03D0,0.001D0,4*0.1D0,3*0D0,
489      &19*10D0,0.00001D0,7*10D0,0.00001D0,10D0,0.00001D0,10D0,0.00001D0,
490      &0.29108D0,0.17412D0,0.45362D0,0.95114D0,8.68604D0,6.23946D0,
491      &1.91923D0,450D0,0.22959D0,1.88863D0,232.68185D0,28.63059D0,0D0,
492      &34.59032D0,25.93594D0,25.96873D0,4.28961D0,4.19124D0,1.41528D0,
493      &0.00977D0,0.00976D0,0.00973D0,267.24501D0,217.49162D0,8.81592D0,
494      &8.80013D0,7*0D0,139*0D0/
495       DATA (PMAS(I,4),I=   1, 500)/12*0D0,658654D0,0D0,0.0872D0,68*0D0,
496      &0.1D0,0.387D0,16*0D0,0.00003D0,2*0D0,15500D0,7804.5D0,5*0D0,
497      &26.762D0,3*0D0,3709D0,5*0D0,0.317D0,2*0D0,0.1244D0,2*0D0,0.14D0,
498      &5*0D0,0.468D0,2*0D0,0.462D0,2*0D0,0.483D0,2*0D0,0.15D0,18*0D0,
499      &44.34D0,0D0,78.88D0,4*0D0,23.96D0,2*0D0,49.1D0,0D0,87.1D0,0D0,
500      &24.6D0,4*0D0,0.0618D0,0.029D0,6*0D0,0.106D0,6*0D0,0.019D0,2*0D0,
501      &7*0.1D0,4*0D0,0.342D0,2*0.387D0,6*0D0,2*0.387D0,6*0D0,0.387D0,
502      &0D0,0.387D0,2*0D0,8*0.387D0,0D0,9*0.387D0,112*0D0,139*0D0/
503       DATA PARF/
504      &  0.5D0,0.25D0, 0.5D0,0.25D0, 1D0, 0.5D0,  0D0,  0D0,  0D0, 0D0,
505      1  0.5D0,  0D0, 0.5D0,  0D0,  1D0,  1D0,  0D0,  0D0,  0D0, 0D0,
506      2  0.5D0,  0D0, 0.5D0,  0D0,  1D0,  1D0,  0D0,  0D0,  0D0, 0D0,
507      3  0.5D0,  0D0, 0.5D0,  0D0,  1D0,  1D0,  0D0,  0D0,  0D0, 0D0,
508      4  0.5D0,  0D0, 0.5D0,  0D0,  1D0,  1D0,  0D0,  0D0,  0D0, 0D0,
509      5  0.5D0,  0D0, 0.5D0,  0D0,  1D0,  1D0,  0D0,  0D0,  0D0, 0D0,
510      6  0.75D0, 0.5D0, 0D0,0.1667D0,0.0833D0,0.1667D0,0D0,0D0,0D0, 0D0,
511      7  0D0,  0D0,  1D0,0.3333D0,0.6667D0,0.3333D0,0D0,0D0,0D0, 0D0,
512      8  0D0,  0D0,  0D0,  0D0,  0D0,  0D0,  0D0,  0D0,  0D0, 0D0,
513      9  0.0099D0, 0.0056D0, 0.199D0, 1.35D0, 4.5D0, 165D0,  4*0D0,
514      & 0.325D0,0.325D0,0.5D0,1.6D0, 5.0D0,  0D0,  0D0,  0D0,  0D0, 0D0,
515      1 0D0,0.11D0,0.16D0,0.048D0,0.50D0,0.45D0,0.55D0,0.60D0,0D0,0D0,
516      2 0.2D0, 0.1D0,  0D0,  0D0,  0D0,  0D0,  0D0,  0D0,  0D0, 0D0,
517      3 60*0D0,
518      4 0.2D0,  0.5D0,  8*0D0,
519      5 1800*0D0/
520       DATA ((VCKM(I,J),J=1,4),I=1,4)/
521      &  0.95113D0,  0.04884D0,  0.00003D0,  0.00000D0,
522      &  0.04884D0,  0.94940D0,  0.00176D0,  0.00000D0,
523      &  0.00003D0,  0.00176D0,  0.99821D0,  0.00000D0,
524      &  0.00000D0,  0.00000D0,  0.00000D0,  1.00000D0/
525  
526 C...PYDAT3, with particle decay parameters and data.
527       DATA (MDCY(I,1),I=   1, 500)/5*0,3*1,6*0,1,0,1,5*0,3*1,6*0,1,0,
528      &4*1,3*0,2*1,40*0,3*1,16*0,3*1,2*0,9*1,0,32*1,2*0,1,3*0,1,2*0,2*1,
529      &2*0,3*1,2*0,4*1,0,5*1,2*0,4*1,2*0,5*1,2*0,6*1,0,7*1,2*0,5*1,2*0,
530      &6*1,2*0,7*1,2*0,8*1,0,75*1,0,7*1,0,1,0,1,0,26*1,146*0/
531       DATA (MDCY(I,2),I=   1, 351)/1,9,17,25,33,41,56,66,2*0,76,80,82,
532      &87,89,143,145,150,2*0,153,162,174,190,210,6*0,289,0,311,334,420,
533      &503,3*0,530,539,40*0,540,541,545,16*0,554,556,561,570,579,581,
534      &583,590,598,604,613,615,617,620,630,636,639,650,656,667,673,736,
535      &739,747,808,810,818,851,853,857,858,861,863,899,900,908,944,945,
536      &953,992,993,997,1028,1029,1033,1034,1043,2*0,1045,3*0,1046,2*0,
537      &1049,1052,2*0,1053,1055,1058,2*0,1062,1063,1066,1069,0,1072,1077,
538      &1079,1082,1084,2*0,1088,1089,1090,1166,2*0,1170,1171,1172,1173,
539      &1174,2*0,1178,1179,1181,1182,1184,1188,0,1189,1193,1197,1201,
540      &1205,1209,1213,2*0,1217,1218,1219,1236,1245,2*0,1254,1255,1256,
541      &1257,1258,1267,2*0,1276,1277,1278,1279,1280,1289,1290,2*0,1299,
542      &1308,1317,1326,1335,1344,1353,1362,0,1371,1380,1389,1398,1407,
543      &1416,1425,1434,1443,1452,1453,1454,1455,1456,1461,1464,1466,1471,
544      &1473,1478,1485,1489,1491,1493,1495,1497,1499,1501,1503,1504,1506,
545      &1508,1510,1512,1514,1516,1518,1520,1522,1523,1525,1527,1541,1543,
546      &1545,1549,1551,1553,1555,1557,1559,1561,1563,1565,1567,1578,1592,
547      &1631,1652,1691,1712,1751,1775,1806,1832,1864,1890,1922,1948,2009,
548      &2160,2406,2615,2877,3155,0,3388,3431,3456,3499,3524,3567,3592,0,
549      &3628,0,3664,0,3700,3708,3716,3724,3727,3751,3777,3801,3807,3814,
550      &3821,3828,3834,3840,3849,3853,3857,3860,3862,3883,3905,3927,3949/
551       DATA (MDCY(I,2),I= 352, 500)/3964,3976,3983,146*0/
552       DATA (MDCY(I,3),I=   1, 500)/5*8,15,2*10,2*0,4,2,5,2,54,2,5,3,
553      &2*0,9,12,16,20,79,6*0,22,0,23,86,83,27,3*0,9,1,40*0,1,4,9,16*0,2,
554      &5,2*9,2*2,7,8,6,9,2*2,3,10,6,3,11,6,11,6,63,3,8,61,2,8,33,2,4,1,
555      &3,2,36,1,8,36,1,8,39,1,4,31,1,4,1,9,2,2*0,1,3*0,3,2*0,3,1,2*0,2,
556      &3,4,2*0,1,3*3,0,5,2,3,2,4,2*0,2*1,76,4,2*0,4*1,4,2*0,1,2,1,2,4,1,
557      &0,7*4,2*0,2*1,17,2*9,2*0,4*1,2*9,2*0,4*1,9,1,9,2*0,8*9,0,9*9,4*1,
558      &5,3,2,5,2,5,7,4,7*2,1,9*2,1,2*2,14,2*2,4,9*2,11,14,39,21,39,21,
559      &39,24,31,26,32,26,32,26,61,151,246,209,262,278,233,0,43,25,43,25,
560      &43,25,36,0,36,0,36,0,3*8,3,24,26,24,6,3*7,2*6,9,2*4,3,2,21,3*22,
561      &15,12,2*7,146*0/
562       DATA (MDME(I,1),I=   1,8000)/6*1,-1,7*1,-1,7*1,-1,7*1,-1,7*1,-1,
563      &7*1,-1,1,7*-1,8*1,2*-1,8*1,2*-1,73*1,-1,2*1,-1,5*1,0,2*-1,6*1,0,
564      &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,
565      &2*-1,3*1,-1,5*1,62*-1,6*1,2*-1,6*1,8*-1,3*1,-1,3*1,-1,3*1,5*-1,
566      &3*1,4*-1,6*1,2*-1,3*1,-1,12*1,62*-1,6*1,2*-1,3*1,-1,9*1,62*-1,
567      &3*1,-1,3*1,-1,1,18*-1,4*1,2*-1,2*1,-1,1225*1,2*-1,248*1,2*-1,
568      &1725*1,2*-1,6*1,2*-1,9*1,-1,3*1,-1,3*1,5*-1,3*1,-1,14*1,2*-1,6*1,
569      &2*-1,67*1,2*-1,6*1,2*-1,4*1,-1,107*1,4011*0/
570       DATA (MDME(I,2),I=   1,8000)/43*102,4*0,102,0,6*53,3*102,4*0,102,
571      &2*0,3*102,4*0,102,2*0,6*102,42,6*102,2*42,2*0,8*41,2*0,36*41,
572      &8*102,0,102,0,102,2*0,21*102,8*32,8*0,16*32,4*0,8*32,9*0,62*53,
573      &8*32,14*0,16*32,7*0,8*32,16*0,62*53,8*32,13*0,62*53,4*32,5*0,
574      &18*53,6*32,4*0,12,2*42,2*11,9*42,0,2,3,15*0,4*42,5*0,3,12*0,2,
575      &3*0,1,0,3,16*0,2*3,15*0,2*42,2*3,18*0,2*3,3*0,1,11*0,22*42,41*0,
576      &2*3,9*0,16*42,45*0,3,10*0,10*42,20*0,2*13,6*0,12,2*0,12,0,12,
577      &14*42,16*0,48,3*13,2*42,9*0,14*42,16*0,48,3*13,2*42,9*0,14*42,
578      &19*0,48,3*13,2*42,6*0,2*11,28*42,5*0,32,3*0,4*32,2*4,0,32,45*0,
579      &14*42,52*0,10*13,2*42,2*11,4*0,2*42,2*11,6*0,2*42,2*11,0,2*42,
580      &2*11,2*42,2*11,2*42,2*11,2*42,2*11,2*42,2*11,2*42,2*11,2*42,2*11,
581      &2*0,3*42,8*0,48,3*13,20*42,4*0,18*42,4*0,9*42,0,162*42,50*0,2*12,
582      &17*0,2*32,33*0,12,9*0,32,2*0,12,11*0,4*32,2*4,5*0,2108*53,4*32,
583      &3*0,6*32,3*0,4*32,3*0,4*32,8*0,8*32,14*0,16*32,12*0,8*32,8*0,
584      &46*32,3*53,12*0,8*32,13*0,66*51,6*32,9*0,9*32,4028*0/
585       DATA (BRAT(I)  ,I=   1, 346)/43*0D0,0.00003D0,0.001765D0,
586      &0.998205D0,35*0D0,1D0,6*0D0,0.1783D0,0.1735D0,0.1131D0,0.2494D0,
587      &0.003D0,0.09D0,0.0027D0,0.01D0,0.0014D0,0.0012D0,2*0.00025D0,
588      &0.0071D0,0.012D0,0.0004D0,0.00075D0,0.00006D0,2*0.00078D0,
589      &0.0034D0,0.08D0,0.011D0,0.0191D0,0.00006D0,0.005D0,0.0133D0,
590      &0.0067D0,0.0005D0,0.0035D0,0.0006D0,0.0015D0,0.00021D0,0.0002D0,
591      &0.00075D0,0.0001D0,0.0002D0,0.0011D0,3*0.0002D0,0.00022D0,
592      &0.0004D0,0.0001D0,2*0.00205D0,2*0.00069D0,0.00025D0,0.00051D0,
593      &0.00025D0,35*0D0,0.153995D0,0.11942D0,0.153984D0,0.119259D0,
594      &0.152272D0,3*0D0,0.033576D0,0.066806D0,0.033576D0,0.066806D0,
595      &0.0335D0,0.066806D0,2*0D0,0.321369D0,0.016494D0,2*0D0,0.016502D0,
596      &0.320615D0,2*0D0,0.00001D0,0.000591D0,6*0D0,2*0.108166D0,
597      &0.108087D0,0D0,0.000001D0,0D0,0.000349D0,0.048707D0,0.768308D0,
598      &4*0D0,0.000227D0,0.064048D0,0D0,0.040621D0,0.002043D0,0.000615D0,
599      &0.006981D0,0.068099D0,62*0D0,0.145835D0,0.113276D0,0.145835D0,
600      &0.113271D0,0.145781D0,0.049002D0,2*0D0,0.032025D0,0.063642D0,
601      &0.032025D0,0.063642D0,0.032022D0,0.063642D0,8*0D0,0.251225D0,
602      &0.0129D0,0.000006D0,0D0,0.0129D0,0.250764D0,0.00038D0,0D0,
603      &0.000008D0,0.000465D0,0.215418D0,5*0D0,2*0.085312D0,0.08531D0,
604      &7*0D0,0.000049D0,0.000774D0,5*0D0,0.000074D0,0D0,0.000417D0/
605       DATA (BRAT(I)  ,I= 347, 651)/0.000015D0,0.000061D0,0.30671D0,
606      &0.689011D0,0D0,0.002889D0,69*0D0,0.000001D0,0.000121D0,
607      &0.001924D0,4*0D0,0.000001D0,0.000184D0,0D0,0.003106D0,0.000015D0,
608      &0.000003D0,2*0D0,0.994646D0,66*0D0,0.000021D0,0.090135D0,2*0D0,
609      &0.000013D0,0.003714D0,0D0,0.906117D0,18*0D0,3*0.215119D0,
610      &0.214724D0,2*0D0,0.06996D0,0.069959D0,0D0,2*1D0,2*0.08D0,0.76D0,
611      &0.08D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,
612      &0.005D0,0.988D0,0.012D0,0.998739D0,0.00079D0,0.00038D0,
613      &0.000046D0,0.000045D0,2*0.34725D0,0.144D0,0.104D0,0.0245D0,
614      &2*0.01225D0,0.0028D0,0.0057D0,0.2112D0,0.1256D0,2*0.1939D0,
615      &2*0.1359D0,0.002D0,0.001D0,0.0006D0,0.999877D0,0.000123D0,
616      &0.99955D0,0.00045D0,2*0.34725D0,0.144D0,0.104D0,0.049D0,0.0028D0,
617      &0.0057D0,0.3923D0,0.321D0,0.2317D0,0.0478D0,0.0049D0,0.0013D0,
618      &0.0003D0,0.0007D0,0.89D0,0.08693D0,0.0221D0,0.00083D0,
619      &2*0.00007D0,0.564D0,0.282D0,0.072D0,0.028D0,0.023D0,2*0.0115D0,
620      &0.005D0,0.003D0,0.6861D0,0.3139D0,2*0.5D0,0.665D0,0.333D0,
621      &0.002D0,0.333D0,0.166D0,0.168D0,0.084D0,0.087D0,0.043D0,0.059D0,
622      &2*0.029D0,0.002D0,0.6352D0,0.2116D0,0.0559D0,0.0173D0,0.0482D0,
623      &0.0318D0,0.666D0,0.333D0,0.001D0,0.332D0,0.166D0,0.168D0,0.084D0,
624      &0.086D0,0.043D0,0.059D0,2*0.029D0,2*0.002D0,0.437D0,0.208D0/
625       DATA (BRAT(I)  ,I= 652, 823)/0.302D0,0.0302D0,0.0212D0,0.0016D0,
626      &0.48947D0,0.34D0,3*0.043D0,0.027D0,0.0126D0,0.0013D0,0.0003D0,
627      &0.00025D0,0.00008D0,0.444D0,2*0.222D0,0.104D0,2*0.004D0,0.07D0,
628      &0.065D0,2*0.005D0,2*0.011D0,5*0.001D0,0.07D0,0.065D0,2*0.005D0,
629      &2*0.011D0,5*0.001D0,0.026D0,0.019D0,0.066D0,0.041D0,0.045D0,
630      &0.076D0,0.0073D0,2*0.0047D0,0.026D0,0.001D0,0.0006D0,0.0066D0,
631      &0.005D0,2*0.003D0,2*0.0006D0,2*0.001D0,0.006D0,0.005D0,0.012D0,
632      &0.0057D0,0.067D0,0.008D0,0.0022D0,0.027D0,0.004D0,0.019D0,
633      &0.012D0,0.002D0,0.009D0,0.0218D0,0.001D0,0.022D0,0.087D0,0.001D0,
634      &0.0019D0,0.0015D0,0.0028D0,0.683D0,0.306D0,0.011D0,0.3D0,0.15D0,
635      &0.16D0,0.08D0,0.13D0,0.06D0,0.08D0,0.04D0,0.034D0,0.027D0,
636      &2*0.002D0,2*0.004D0,2*0.002D0,0.034D0,0.027D0,2*0.002D0,
637      &2*0.004D0,2*0.002D0,0.0365D0,0.045D0,0.073D0,0.062D0,3*0.021D0,
638      &0.0061D0,0.015D0,0.025D0,0.0088D0,0.074D0,0.0109D0,0.0041D0,
639      &0.002D0,0.0035D0,0.0011D0,0.001D0,0.0027D0,2*0.0016D0,0.0018D0,
640      &0.011D0,0.0063D0,0.0052D0,0.018D0,0.016D0,0.0034D0,0.0036D0,
641      &0.0009D0,0.0006D0,0.015D0,0.0923D0,0.018D0,0.022D0,0.0077D0,
642      &0.009D0,0.0075D0,0.024D0,0.0085D0,0.067D0,0.0511D0,0.017D0,
643      &0.0004D0,0.0028D0,0.619D0,0.381D0,0.3D0,0.15D0,0.16D0,0.08D0,
644      &0.13D0,0.06D0,0.08D0,0.04D0,0.01D0,2*0.02D0,0.03D0,2*0.005D0/
645       DATA (BRAT(I)  ,I= 824, 991)/2*0.02D0,0.03D0,2*0.005D0,0.015D0,
646      &0.037D0,0.028D0,0.079D0,0.095D0,0.052D0,0.0078D0,4*0.001D0,
647      &0.028D0,0.033D0,0.026D0,0.05D0,0.01D0,4*0.005D0,0.25D0,0.0952D0,
648      &0.94D0,0.06D0,2*0.4D0,2*0.1D0,1D0,0.0602D0,0.0601D0,0.8797D0,
649      &0.135D0,0.865D0,0.02D0,0.055D0,2*0.005D0,0.008D0,0.012D0,0.02D0,
650      &0.055D0,2*0.005D0,0.008D0,0.012D0,0.01D0,0.03D0,0.0035D0,0.011D0,
651      &0.0055D0,0.0042D0,0.009D0,0.018D0,0.015D0,0.0185D0,0.0135D0,
652      &0.025D0,0.0004D0,0.0007D0,0.0008D0,0.0014D0,0.0019D0,0.0025D0,
653      &0.4291D0,0.08D0,0.07D0,0.02D0,0.015D0,0.005D0,1D0,0.3D0,0.15D0,
654      &0.16D0,0.08D0,0.13D0,0.06D0,0.08D0,0.04D0,0.02D0,0.055D0,
655      &2*0.005D0,0.008D0,0.012D0,0.02D0,0.055D0,2*0.005D0,0.008D0,
656      &0.012D0,0.01D0,0.03D0,0.0035D0,0.011D0,0.0055D0,0.0042D0,0.009D0,
657      &0.018D0,0.015D0,0.0185D0,0.0135D0,0.025D0,0.0004D0,0.0007D0,
658      &0.0008D0,0.0014D0,0.0019D0,0.0025D0,0.4291D0,0.08D0,0.07D0,
659      &0.02D0,0.015D0,0.005D0,1D0,0.3D0,0.15D0,0.16D0,0.08D0,0.13D0,
660      &0.06D0,0.08D0,0.04D0,0.02D0,0.055D0,2*0.005D0,0.008D0,0.012D0,
661      &0.02D0,0.055D0,2*0.005D0,0.008D0,0.012D0,0.01D0,0.03D0,0.0035D0,
662      &0.011D0,0.0055D0,0.0042D0,0.009D0,0.018D0,0.015D0,0.0185D0,
663      &0.0135D0,0.025D0,2*0.0002D0,0.0007D0,2*0.0004D0,0.0014D0,0.001D0,
664      &0.0009D0,0.0025D0,0.4291D0,0.08D0,0.07D0,0.02D0,0.015D0,0.005D0/
665       DATA (BRAT(I)  ,I= 992,1183)/1D0,2*0.3D0,2*0.2D0,0.047D0,0.122D0,
666      &0.006D0,0.012D0,0.035D0,0.012D0,0.035D0,0.003D0,0.007D0,0.15D0,
667      &0.037D0,0.008D0,0.002D0,0.05D0,0.015D0,0.003D0,0.001D0,0.014D0,
668      &0.042D0,0.014D0,0.042D0,0.24D0,0.065D0,0.012D0,0.003D0,0.001D0,
669      &0.002D0,0.001D0,0.002D0,0.014D0,0.003D0,1D0,2*0.3D0,2*0.2D0,1D0,
670      &0.0252D0,0.0248D0,0.0267D0,0.015D0,0.045D0,0.015D0,0.045D0,
671      &0.7743D0,0.029D0,0.22D0,0.78D0,1D0,0.331D0,0.663D0,0.006D0,
672      &0.663D0,0.331D0,0.006D0,1D0,0.999D0,0.001D0,0.88D0,2*0.06D0,
673      &0.639D0,0.358D0,0.002D0,0.001D0,1D0,0.88D0,2*0.06D0,0.516D0,
674      &0.483D0,0.001D0,0.88D0,2*0.06D0,0.9988D0,0.0001D0,0.0006D0,
675      &0.0004D0,0.0001D0,0.667D0,0.333D0,0.9954D0,0.0011D0,0.0035D0,
676      &0.333D0,0.667D0,0.676D0,0.234D0,0.085D0,0.005D0,2*1D0,0.018D0,
677      &2*0.005D0,0.003D0,0.002D0,2*0.006D0,0.018D0,2*0.005D0,0.003D0,
678      &0.002D0,2*0.006D0,0.0066D0,0.025D0,0.016D0,0.0088D0,2*0.005D0,
679      &0.0058D0,0.005D0,0.0055D0,4*0.004D0,2*0.002D0,2*0.004D0,0.003D0,
680      &0.002D0,2*0.003D0,3*0.002D0,2*0.001D0,0.002D0,2*0.001D0,
681      &2*0.002D0,0.0013D0,0.0018D0,5*0.001D0,4*0.003D0,2*0.005D0,
682      &2*0.002D0,2*0.001D0,2*0.002D0,2*0.001D0,0.2432D0,0.057D0,
683      &2*0.035D0,0.15D0,2*0.075D0,0.03D0,2*0.015D0,2*0.08D0,0.76D0,
684      &0.08D0,4*1D0,2*0.08D0,0.76D0,0.08D0,1D0,2*0.5D0,1D0,2*0.5D0/
685       DATA (BRAT(I)  ,I=1184,1377)/2*0.08D0,0.76D0,0.08D0,1D0,2*0.08D0,
686      &0.76D0,3*0.08D0,0.76D0,3*0.08D0,0.76D0,3*0.08D0,0.76D0,3*0.08D0,
687      &0.76D0,3*0.08D0,0.76D0,3*0.08D0,0.76D0,0.08D0,2*1D0,2*0.105D0,
688      &0.04D0,0.0077D0,0.02D0,0.0235D0,0.0285D0,0.0435D0,0.0011D0,
689      &0.0022D0,0.0044D0,0.4291D0,0.08D0,0.07D0,0.02D0,0.015D0,0.005D0,
690      &2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,
691      &2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,
692      &4*1D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,
693      &0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,
694      &0.005D0,4*1D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,
695      &0.015D0,0.005D0,1D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,
696      &0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,
697      &0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,
698      &0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,
699      &0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,
700      &0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,
701      &0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,
702      &0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,
703      &0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0,
704      &0.015D0,0.005D0,2*0.105D0,0.04D0,0.5D0,0.08D0,0.14D0,0.01D0/
705       DATA (BRAT(I)  ,I=1378,1580)/0.015D0,0.005D0,2*0.105D0,0.04D0,
706      &0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,2*0.105D0,0.04D0,
707      &0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,2*0.105D0,0.04D0,
708      &0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,2*0.105D0,0.04D0,
709      &0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,2*0.105D0,0.04D0,
710      &0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,2*0.105D0,0.04D0,
711      &0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,2*0.105D0,0.04D0,
712      &0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,2*0.105D0,0.04D0,
713      &0.5D0,0.08D0,0.14D0,0.01D0,0.015D0,0.005D0,4*1D0,0.52D0,0.26D0,
714      &0.11D0,2*0.055D0,0.333D0,0.334D0,0.333D0,0.667D0,0.333D0,0.28D0,
715      &0.14D0,0.313D0,0.157D0,0.11D0,0.667D0,0.333D0,0.28D0,0.14D0,
716      &0.313D0,0.157D0,0.11D0,0.36D0,0.18D0,0.03D0,2*0.015D0,2*0.2D0,
717      &4*0.25D0,0.667D0,0.333D0,0.667D0,0.333D0,0.667D0,0.333D0,0.667D0,
718      &0.333D0,4*0.5D0,0.007D0,0.993D0,1D0,0.667D0,0.333D0,0.667D0,
719      &0.333D0,0.667D0,0.333D0,0.667D0,0.333D0,8*0.5D0,0.02D0,0.98D0,
720      &1D0,4*0.5D0,3*0.146D0,3*0.05D0,0.15D0,2*0.05D0,4*0.024D0,0.066D0,
721      &0.667D0,0.333D0,0.667D0,0.333D0,4*0.25D0,0.667D0,0.333D0,0.667D0,
722      &0.333D0,2*0.5D0,0.273D0,0.727D0,0.667D0,0.333D0,0.667D0,0.333D0,
723      &4*0.5D0,0.35D0,0.65D0,2*0.0083D0,0.1866D0,0.324D0,0.184D0,
724      &0.027D0,0.001D0,0.093D0,0.087D0,0.078D0,0.0028D0,3*0.014D0/
725       DATA (BRAT(I)  ,I=1581,3853)/0.008D0,0.024D0,0.008D0,0.024D0,
726      &0.425D0,0.02D0,0.185D0,0.088D0,0.043D0,0.067D0,0.066D0,2108*0D0,
727      &0.017431D0,0.054048D0,0.857694D0,2*0D0,0.00025D0,0.070578D0,0D0,
728      &0.022748D0,0.026576D0,0.359486D0,0.561581D0,2*0D0,0.000104D0,
729      &0.029504D0,0.011185D0,0.034681D0,0.550354D0,2*0D0,0.00016D0,
730      &0.045287D0,0.358333D0,0.445781D0,0D0,0.554219D0,0.144051D0,
731      &2*0.351902D0,0D0,0.082107D0,0.029566D0,0.001511D0,0.000726D0,
732      &0.004518D0,0.006522D0,0.004518D0,0.006522D0,0.004513D0,3*0D0,
733      &0.002908D0,0.000973D0,0.002908D0,0.000973D0,0.002908D0,
734      &0.000973D0,2*0D0,0.143982D0,0.489888D0,0.1951D0,0D0,0.114302D0,
735      &0.008426D0,0.014868D0,0.000763D0,2*0D0,0.000763D0,0.01484D0,
736      &0.000003D0,2*0D0,0.000027D0,0.001945D0,5*0D0,3*0.00503D0,0D0,
737      &0.133776D0,0.003284D0,0.37169D0,0.006838D0,2*0.030954D0,
738      &0.00163D0,0D0,0.047224D0,0.073737D0,0.047224D0,0.073732D0,
739      &0.047179D0,3*0D0,0.034761D0,0.009166D0,0.034761D0,0.009166D0,
740      &0.034759D0,0.009166D0,2*0D0,4*0.009069D0,0.510147D0,0.453576D0,
741      &6*0D0,1D0,6*0D0,1D0,4*0.001128D0,0.571047D0,0.382288D0,
742      &0.042153D0,4*0.016597D0,0.93361D0,0D0,4*0.016597D0,0.93361D0,0D0,
743      &4*0.05515D0,0.34469D0,0D0,0.228998D0,0.164208D0,0.041503D0,
744      &0.850973D0,0.005411D0,0.045025D0,0.098591D0,0.849898D0/
745       DATA (BRAT(I)  ,I=3854,3984)/0.021617D0,0.030018D0,0.098466D0,
746      &0.294448D0,0.10945D0,0.596102D0,0.389906D0,0.610094D0,3*0.0633D0,
747      &0.063299D0,0.063295D0,0.056281D0,2*0D0,6*0.020495D0,2*0D0,
748      &0.327919D0,0.04099D0,0.045236D0,0.090112D0,0D0,0.19874D0,
749      &0.010204D0,0.000003D0,0.010205D0,0.198356D0,0.000151D0,
750      &0.000006D0,0.000367D0,0.081967D0,0.19874D0,0.010204D0,0.000003D0,
751      &0.010205D0,0.198356D0,0.000151D0,0.000006D0,0.000367D0,
752      &0.081967D0,4*0D0,0.198776D0,0.010206D0,0.000003D0,0.010207D0,
753      &0.19839D0,0.000151D0,0.000006D0,0.000367D0,0.081893D0,0.198776D0,
754      &0.010206D0,0.000003D0,0.010207D0,0.19839D0,0.000151D0,0.000006D0,
755      &0.000367D0,0.081893D0,4*0D0,0.199344D0,0.010234D0,0.000003D0,
756      &0.010236D0,0.198928D0,0.000149D0,0.000006D0,0.000368D0,
757      &0.080733D0,0.199344D0,0.010234D0,0.000003D0,0.010236D0,
758      &0.198928D0,0.000149D0,0.000006D0,0.000368D0,0.080733D0,4*0D0,
759      &0.184738D0,0.104588D0,0.184738D0,0.104587D0,0.184731D0,0.09582D0,
760      &0.022902D0,0.008429D0,0.015602D0,0.022902D0,0.008429D0,
761      &0.015602D0,0.022902D0,0.008429D0,0.015602D0,0.28959D0,0.01487D0,
762      &0.000008D0,0.01487D0,0.289061D0,0.000492D0,0.000009D0,0.000536D0,
763      &0.27911D0,2*0.037151D0,0.03715D0,0.090266D0,2*0.001805D0,
764      &0.090266D0,0.001805D0,0.812263D0,0.00179D0,0.090428D0,0.001809D0/
765       DATA (BRAT(I)  ,I=3985,8000)/0.001808D0,0.090428D0,0.001808D0,
766      &0.81372D0,0D0,4011*0D0/
767       DATA (KFDP(I,1),I=   1, 377)/21,22,23,4*-24,25,21,22,23,4*24,25,
768      &21,22,23,4*-24,25,21,22,23,4*24,25,21,22,23,4*-24,25,21,22,23,
769      &4*24,25,37,1000022,1000023,1000025,1000035,1000021,1000039,21,22,
770      &23,4*-24,25,2*-37,21,22,23,4*24,25,2*37,22,23,-24,25,23,24,-12,
771      &22,23,-24,25,23,24,-12,-14,48*16,22,23,-24,25,23,24,22,23,-24,25,
772      &-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,
773      &3,4,5,6,7,8,11,12,13,14,15,16,17,18,4*-1,4*-3,4*-5,4*-7,-11,-13,
774      &-15,-17,1,2,3,4,5,6,7,8,11,13,15,17,21,2*22,23,24,1000022,
775      &2*1000023,3*1000025,4*1000035,2*1000024,2*1000037,1000001,
776      &2000001,1000001,-1000001,1000002,2000002,1000002,-1000002,
777      &1000003,2000003,1000003,-1000003,1000004,2000004,1000004,
778      &-1000004,1000005,2000005,1000005,-1000005,1000006,2000006,
779      &1000006,-1000006,1000011,2000011,1000011,-1000011,1000012,
780      &2000012,1000012,-1000012,1000013,2000013,1000013,-1000013,
781      &1000014,2000014,1000014,-1000014,1000015,2000015,1000015,
782      &-1000015,1000016,2000016,1000016,-1000016,1,2,3,4,5,6,7,8,11,12,
783      &13,14,15,16,17,18,24,37,2*23,25,35,4*-1,4*-3,4*-5,4*-7,-11,-13,
784      &-15,-17,3*24,1,2,3,4,5,6,7,8,11,13,15,17,21,2*22,23,24,23,25,24,
785      &37,23,25,36,1000022,2*1000023,3*1000025,4*1000035,2*1000024,
786      &2*1000037,1000001,2000001,1000001,-1000001,1000002,2000002/
787       DATA (KFDP(I,1),I= 378, 580)/1000002,-1000002,1000003,2000003,
788      &1000003,-1000003,1000004,2000004,1000004,-1000004,1000005,
789      &2000005,1000005,-1000005,1000006,2000006,1000006,-1000006,
790      &1000011,2000011,1000011,-1000011,1000012,2000012,1000012,
791      &-1000012,1000013,2000013,1000013,-1000013,1000014,2000014,
792      &1000014,-1000014,1000015,2000015,1000015,-1000015,1000016,
793      &2000016,1000016,-1000016,1,2,3,4,5,6,7,8,11,13,15,17,21,2*22,23,
794      &24,23,25,24,37,1000022,2*1000023,3*1000025,4*1000035,2*1000024,
795      &2*1000037,1000001,2000001,1000001,-1000001,1000002,2000002,
796      &1000002,-1000002,1000003,2000003,1000003,-1000003,1000004,
797      &2000004,1000004,-1000004,1000005,2000005,1000005,-1000005,
798      &1000006,2000006,1000006,-1000006,1000011,2000011,1000011,
799      &-1000011,1000012,2000012,1000012,-1000012,1000013,2000013,
800      &1000013,-1000013,1000014,2000014,1000014,-1000014,1000015,
801      &2000015,1000015,-1000015,1000016,2000016,1000016,-1000016,-1,-3,
802      &-5,-7,-11,-13,-15,-17,24,2*1000022,2*1000023,2*1000025,2*1000035,
803      &1000006,2000006,1000006,2000006,-1000001,-1000003,-1000011,
804      &-1000013,-1000015,-2000015,1,2,3,4,5,6,11,13,15,2,82,-11,-13,2*2,
805      &-12,-14,-16,2*-2,2*-4,-2,-4,2*22,211,111,221,13,11,213,-213,221,
806      &223,321,130,310,111,331,111,211,-12,12,-14,14,211,111,22,-13,-11/
807       DATA (KFDP(I,1),I= 581, 992)/2*211,213,113,221,223,321,211,331,
808      &22,111,211,2*22,211,22,111,211,22,211,221,111,11,211,111,2*211,
809      &321,130,310,221,111,211,111,130,310,321,2*311,321,311,323,313,
810      &323,313,321,3*311,-13,3*211,12,14,311,2*321,311,321,313,323,313,
811      &323,311,4*321,211,111,3*22,111,321,130,-213,113,213,211,22,111,
812      &11,13,211,321,130,310,221,211,111,11*-11,11*-13,-311,-313,-311,
813      &-313,-20313,2*-311,-313,-311,-313,2*111,2*221,2*331,2*113,2*223,
814      &2*333,-311,-313,2*-321,211,-311,-321,333,-311,-313,-321,211,
815      &2*-321,2*-311,-321,211,113,421,2*411,421,411,423,413,423,413,421,
816      &411,8*-11,8*-13,-321,-323,-321,-323,-311,2*-313,-311,-313,2*-311,
817      &-321,-10323,-321,-323,-321,-311,2*-313,211,111,333,3*-321,-311,
818      &-313,-321,-313,310,333,211,2*-321,-311,-313,-311,211,-321,3*-311,
819      &211,113,321,2*421,411,421,413,423,413,423,411,421,-15,5*-11,
820      &5*-13,221,331,333,221,331,333,10221,211,213,211,213,321,323,321,
821      &323,2212,221,331,333,221,2*2,2*431,421,411,423,413,82,11,13,82,
822      &443,82,6*12,6*14,2*16,3*-411,3*-413,2*-411,2*-413,2*441,2*443,
823      &2*20443,2*2,2*4,2,4,511,521,511,523,513,523,513,521,511,6*12,
824      &6*14,2*16,3*-421,3*-423,2*-421,2*-423,2*441,2*443,2*20443,2*2,
825      &2*4,2,4,521,511,521,513,523,513,523,511,521,6*12,6*14,2*16,
826      &3*-431,3*-433,2*-431,2*-433,3*441,3*443,3*20443,2*2,2*4,2,4,531/
827       DATA (KFDP(I,1),I= 993,1402)/521,511,523,513,16,2*4,2*12,2*14,
828      &2*16,4*2,4*4,2*-11,2*-13,2*-1,2*-3,2*-11,2*-13,2*-1,541,511,521,
829      &513,523,21,11,13,15,1,2,3,4,21,22,553,21,2112,2212,2*2112,2212,
830      &2112,2*2212,2112,-12,3122,3212,3112,2212,2*2112,-12,2*3122,3222,
831      &3112,2212,2112,2212,3122,3222,3212,3122,3112,-12,-14,-12,3322,
832      &3312,2*3122,3212,3322,3312,3122,3322,3312,-12,2*4122,7*-11,7*-13,
833      &2*2224,2*2212,2*2214,2*3122,2*3212,2*3214,5*3222,4*3224,2*3322,
834      &3324,2*2224,7*2212,5*2214,2*2112,2*2114,2*3122,2*3212,2*3214,
835      &2*3222,2*3224,4*2,3,2*2,1,2*2,-11,-13,2*2,4*4122,-11,-13,2*2,
836      &3*4132,3*4232,-11,-13,2*2,4332,-11,-13,2*2,-11,-13,2*2,-11,-13,
837      &2*2,-11,-13,2*2,-11,-13,2*2,-11,-13,2*2,-11,-13,2*2,2*5122,-12,
838      &-14,-16,5*4122,441,443,20443,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,
839      &2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,4*5122,-12,-14,-16,2*-2,
840      &2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,2*5132,2*5232,-12,-14,-16,
841      &2*-2,2*-4,-2,-4,5332,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,
842      &2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,
843      &2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,
844      &-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,
845      &-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,
846      &2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2/
847       DATA (KFDP(I,1),I=1403,1708)/2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,
848      &-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,-12,
849      &-14,-16,2*-2,2*-4,-2,-4,-12,-14,-16,2*-2,2*-4,-2,-4,221,223,221,
850      &223,211,111,321,130,310,213,113,-213,321,311,321,311,323,313,
851      &2*311,321,311,321,313,323,321,211,111,321,130,310,2*211,313,-313,
852      &323,-323,421,411,423,413,411,421,413,423,411,421,423,413,443,
853      &2*82,521,511,523,513,511,521,513,523,521,511,523,513,511,521,513,
854      &523,553,2*21,213,-213,113,213,10211,10111,-10211,2*221,213,2*113,
855      &-213,2*321,2*311,113,323,2*313,323,313,-313,323,-323,423,2*413,
856      &2*423,413,443,82,523,2*513,2*523,2*513,523,553,21,11,13,82,4*443,
857      &10441,20443,445,441,11,13,15,1,2,3,4,21,22,2*553,10551,20553,555,
858      &1000039,-1000024,-1000037,1000022,1000023,1000025,1000035,
859      &1000002,2000002,1000002,2000002,1000021,3*-12,3*-14,3*-16,12,11,
860      &12,11,12,11,14,13,14,13,14,13,16,15,16,15,16,15,1000039,1000024,
861      &1000037,1000022,1000023,1000025,1000035,1000001,2000001,1000001,
862      &2000001,1000021,3*-11,3*-13,3*-15,1000039,-1000024,-1000037,
863      &1000022,1000023,1000025,1000035,1000004,2000004,1000004,2000004,
864      &1000021,3*-12,3*-14,3*-16,12,11,12,11,12,11,14,13,14,13,14,13,16,
865      &15,16,15,16,15,1000039,1000024,1000037,1000022,1000023,1000025,
866      &1000035,1000003,2000003,1000003,2000003,1000021,3*-11,3*-13/
867       DATA (KFDP(I,1),I=1709,1966)/3*-15,1000039,-1000024,-1000037,
868      &1000022,1000023,1000025,1000035,1000006,2000006,1000006,2000006,
869      &1000021,3*-12,3*-14,3*-16,12,11,12,11,12,11,14,13,14,13,14,13,16,
870      &15,16,15,16,15,1000039,1000024,1000037,1000022,1000023,1000025,
871      &1000035,1000005,2000005,1000005,2000005,1000021,1000022,1000016,
872      &-1000015,3*-11,3*-13,3*-15,1000039,-1000024,-1000037,1000022,
873      &1000023,1000025,1000035,1000012,2000012,1000012,2*12,2*14,2*16,
874      &3*-14,3*-16,3*-2,3*-4,3*-6,1000039,1000024,1000037,1000022,
875      &1000023,1000025,1000035,1000011,2000011,1000011,2000011,3*-13,
876      &3*-15,3*-1,3*-3,3*-5,1000039,-1000024,-1000037,1000022,1000023,
877      &1000025,1000035,1000014,2000014,1000014,2000014,2*12,2*14,2*16,
878      &3*-12,3*-16,3*-2,3*-4,3*-6,1000039,1000024,1000037,1000022,
879      &1000023,1000025,1000035,1000013,2000013,1000013,2000013,3*-11,
880      &3*-15,3*-1,3*-3,3*-5,1000039,-1000024,-1000037,1000022,1000023,
881      &1000025,1000035,1000016,2000016,1000016,2000016,2*12,2*14,2*16,
882      &3*-12,3*-14,3*-2,3*-4,3*-6,1000039,1000024,1000037,1000022,
883      &1000023,1000025,1000035,1000015,2000015,1000015,2000015,3*-11,
884      &3*-13,3*-1,3*-3,3*-5,1000039,1000001,-1000001,2000001,-2000001,
885      &1000002,-1000002,2000002,-2000002,1000003,-1000003,2000003,
886      &-2000003,1000004,-1000004,2000004,-2000004,1000005,-1000005/
887       DATA (KFDP(I,1),I=1967,2235)/2000005,-2000005,1000006,-1000006,
888      &2000006,-2000006,6*1000022,6*1000023,6*1000025,6*1000035,1000024,
889      &-1000024,1000024,-1000024,1000024,-1000024,1000037,-1000037,
890      &1000037,-1000037,1000037,-1000037,5*1000039,4,1,-12,12,-12,12,
891      &-12,12,-12,12,-12,12,-12,12,-14,14,-14,14,-14,14,-14,14,-14,14,
892      &-14,14,-16,16,-16,16,-16,16,-16,16,-16,16,-16,16,-12,12,-11,11,
893      &-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,
894      &-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-14,14,-13,13,
895      &-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,
896      &-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-16,16,-15,15,
897      &-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,
898      &-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,5*1000039,
899      &16*1000022,1000024,-1000024,1000024,-1000024,1000024,-1000024,
900      &1000024,-1000024,1000024,-1000024,1000024,-1000024,1000037,
901      &-1000037,1000037,-1000037,1000037,-1000037,1000037,-1000037,
902      &1000037,-1000037,1000037,-1000037,1000024,-1000024,1000037,
903      &-1000037,1000001,-1000001,2000001,-2000001,1000002,-1000002,
904      &2000002,-2000002,1000003,-1000003,2000003,-2000003,1000004,
905      &-1000004,2000004,-2000004,1000005,-1000005,2000005,-2000005,
906      &1000006,-1000006,2000006,-2000006,1000011,-1000011,2000011/
907       DATA (KFDP(I,1),I=2236,2523)/-2000011,1000012,-1000012,2000012,
908      &-2000012,1000013,-1000013,2000013,-2000013,1000014,-1000014,
909      &2000014,-2000014,1000015,-1000015,2000015,-2000015,1000016,
910      &-1000016,2000016,-2000016,5*1000021,-12,12,-12,12,-12,12,-12,12,
911      &-12,12,-12,12,-14,14,-14,14,-14,14,-14,14,-14,14,-14,14,-16,16,
912      &-16,16,-16,16,-16,16,-16,16,-16,16,-12,12,-11,11,-12,12,-11,11,
913      &-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,
914      &-11,11,-12,12,-11,11,-12,12,-11,11,-14,14,-13,13,-14,14,-13,13,
915      &-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,
916      &-13,13,-14,14,-13,13,-14,14,-13,13,-16,16,-15,15,-16,16,-15,15,
917      &-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,
918      &-15,15,-16,16,-15,15,-16,16,-15,15,2*1000039,6*1000022,6*1000023,
919      &6*1000025,6*1000035,1000022,1000023,1000025,1000035,1000002,
920      &2000002,-1000001,-2000001,1000004,2000004,-1000003,-2000003,
921      &1000006,2000006,-1000005,-2000005,1000012,2000012,-1000011,
922      &-2000011,1000014,2000014,-1000013,-2000013,1000016,2000016,
923      &-1000015,-2000015,2*1000021,-12,12,-11,-12,12,-11,-12,12,-11,-12,
924      &12,-11,-12,12,-11,-12,12,-11,-14,-13,-14,-13,-14,-13,-14,14,-13,
925      &-14,14,-13,-14,14,-13,-16,-15,-16,-15,-16,-15,-16,-15,-16,-15,
926      &-16,-15,-12,2*-11,12,-12,2*-11,12,-12,2*-11,12,-12,2*-11,12,-12/
927       DATA (KFDP(I,1),I=2524,2794)/2*-11,12,-12,2*-11,12,-12,2*-11,12,
928      &-12,2*-11,12,-12,2*-11,12,-14,2*-13,14,-14,2*-13,14,-14,2*-13,14,
929      &-14,2*-13,14,-14,2*-13,14,-14,2*-13,14,-14,2*-13,14,-14,2*-13,14,
930      &-14,2*-13,14,-16,2*-15,16,-16,2*-15,16,-16,2*-15,16,-16,2*-15,16,
931      &-16,2*-15,16,-16,2*-15,16,-16,2*-15,16,-16,2*-15,16,-16,2*-15,16,
932      &5*1000039,16*1000022,16*1000023,1000024,-1000024,1000024,
933      &-1000024,1000024,-1000024,1000024,-1000024,1000024,-1000024,
934      &1000024,-1000024,1000037,-1000037,1000037,-1000037,1000037,
935      &-1000037,1000037,-1000037,1000037,-1000037,1000037,-1000037,
936      &1000024,-1000024,1000037,-1000037,1000001,-1000001,2000001,
937      &-2000001,1000002,-1000002,2000002,-2000002,1000003,-1000003,
938      &2000003,-2000003,1000004,-1000004,2000004,-2000004,1000005,
939      &-1000005,2000005,-2000005,1000006,-1000006,2000006,-2000006,
940      &1000011,-1000011,2000011,-2000011,1000012,-1000012,2000012,
941      &-2000012,1000013,-1000013,2000013,-2000013,1000014,-1000014,
942      &2000014,-2000014,1000015,-1000015,2000015,-2000015,1000016,
943      &-1000016,2000016,-2000016,5*1000021,-12,12,-12,12,-12,12,-12,12,
944      &-12,12,-12,12,-14,14,-14,14,-14,14,-14,14,-14,14,-14,14,-16,16,
945      &-16,16,-16,16,-16,16,-16,16,-16,16,-12,12,-11,11,-12,12,-11,11,
946      &-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12/
947       DATA (KFDP(I,1),I=2795,3070)/-11,11,-12,12,-11,11,-12,12,-11,11,
948      &-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,
949      &-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,
950      &-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,
951      &-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,
952      &5*1000039,16*1000022,16*1000023,16*1000025,1000024,-1000024,
953      &1000024,-1000024,1000024,-1000024,1000024,-1000024,1000024,
954      &-1000024,1000024,-1000024,1000037,-1000037,1000037,-1000037,
955      &1000037,-1000037,1000037,-1000037,1000037,-1000037,1000037,
956      &-1000037,1000024,-1000024,1000037,-1000037,1000001,-1000001,
957      &2000001,-2000001,1000002,-1000002,2000002,-2000002,1000003,
958      &-1000003,2000003,-2000003,1000004,-1000004,2000004,-2000004,
959      &1000005,-1000005,2000005,-2000005,1000006,-1000006,2000006,
960      &-2000006,1000011,-1000011,2000011,-2000011,1000012,-1000012,
961      &2000012,-2000012,1000013,-1000013,2000013,-2000013,1000014,
962      &-1000014,2000014,-2000014,1000015,-1000015,2000015,-2000015,
963      &1000016,-1000016,2000016,-2000016,5*1000021,-12,12,-12,12,-12,12,
964      &-12,12,-12,12,-12,12,-14,14,-14,14,-14,14,-14,14,-14,14,-14,14,
965      &-16,16,-16,16,-16,16,-16,16,-16,16,-16,16,-12,12,-11,11,-12,12,
966      &-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11,-12,12,-11,11/
967       DATA (KFDP(I,1),I=3071,3398)/-12,12,-11,11,-12,12,-11,11,-12,12,
968      &-11,11,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,
969      &-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,-13,13,-14,14,
970      &-13,13,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,
971      &-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,-15,15,-16,16,
972      &-15,15,2*1000039,15*1000024,6*1000022,6*1000023,6*1000025,
973      &6*1000035,1000022,1000023,1000025,1000035,1000002,2000002,
974      &-1000001,-2000001,1000004,2000004,-1000003,-2000003,1000006,
975      &2000006,-1000005,-2000005,1000012,2000012,-1000011,-2000011,
976      &1000014,2000014,-1000013,-2000013,1000016,2000016,-1000015,
977      &-2000015,2*1000021,-12,12,-11,-12,12,-11,-12,12,-11,-12,12,-11,
978      &-12,12,-11,-12,12,-11,-14,14,-13,-14,14,-13,-14,14,-13,-14,14,
979      &-13,-14,14,-13,-14,14,-13,-16,16,-15,-16,16,-15,-16,16,-15,-16,
980      &16,-15,-16,16,-15,-16,16,-15,-12,2*-11,12,-12,2*-11,12,-12,2*-11,
981      &12,-12,2*-11,12,-12,2*-11,12,-12,2*-11,12,-12,2*-11,12,-12,2*-11,
982      &12,-12,2*-11,12,-14,2*-13,14,-14,2*-13,14,-14,2*-13,14,-14,2*-13,
983      &14,-14,2*-13,14,-14,2*-13,14,-14,2*-13,14,-14,2*-13,14,-14,2*-13,
984      &14,-16,2*-15,16,-16,2*-15,16,-16,2*-15,16,-16,2*-15,16,-16,2*-15,
985      &16,-16,2*-15,16,-16,2*-15,16,-16,2*-15,16,-16,2*-15,16,1000039,
986      &-1000024,-1000037,1000022,1000023,1000025,1000035,4*1000001/
987       DATA (KFDP(I,1),I=3399,3676)/1000002,2000002,1000002,2000002,
988      &1000021,3*-12,3*-14,3*-16,12,11,12,11,12,11,14,13,14,13,14,13,16,
989      &15,16,15,16,15,1000039,1000024,1000037,1000022,1000023,1000025,
990      &1000035,4*1000002,1000001,2000001,1000001,2000001,1000021,3*-11,
991      &3*-13,3*-15,1000039,-1000024,-1000037,1000022,1000023,1000025,
992      &1000035,4*1000003,1000004,2000004,1000004,2000004,1000021,3*-12,
993      &3*-14,3*-16,12,11,12,11,12,11,14,13,14,13,14,13,16,15,16,15,16,
994      &15,1000039,1000024,1000037,1000022,1000023,1000025,1000035,
995      &4*1000004,1000003,2000003,1000003,2000003,1000021,3*-11,3*-13,
996      &3*-15,1000039,-1000024,-1000037,1000022,1000023,1000025,1000035,
997      &4*1000005,1000006,2000006,1000006,2000006,1000021,3*-12,3*-14,
998      &3*-16,12,11,12,11,12,11,14,13,14,13,14,13,16,15,16,15,16,15,
999      &1000039,1000024,1000037,1000022,1000023,1000025,1000035,
1000      &4*1000006,1000005,2000005,1000005,2000005,1000021,3*-11,3*-13,
1001      &3*-15,1000039,-1000024,-1000037,1000022,1000023,1000025,1000035,
1002      &4*1000011,1000012,2000012,1000012,2000012,2*12,2*14,2*16,3*-14,
1003      &3*-16,3*-2,3*-4,3*-6,1000039,-1000024,-1000037,1000022,1000023,
1004      &1000025,1000035,4*1000013,1000014,2000014,1000014,2000014,2*12,
1005      &2*14,2*16,3*-12,3*-16,3*-2,3*-4,3*-6,1000039,-1000024,-1000037,
1006      &1000022,1000023,1000025,1000035,4*1000015,1000016,2000016/
1007       DATA (KFDP(I,1),I=3677,8000)/1000016,2000016,2*12,2*14,2*16,
1008      &3*-12,3*-14,3*-2,3*-4,3*-6,3,4,5,6,11,13,15,21,2*4,2,4,24,-11,
1009      &-13,-15,3,4,5,6,11,13,15,21,5,6,21,2*24,2*3000211,2*22,2*23,1,2,
1010      &3,4,5,6,7,8,11,12,13,14,15,16,17,18,2*24,3*3000211,24,4*-1,4*-3,
1011      &4*-5,4*-7,-11,-13,-15,-17,22,23,22,23,24,3000211,24,3000211,1,2,
1012      &3,4,5,6,7,8,11,12,13,14,15,16,17,18,1,2,3,4,5,6,1,2,3,4,5,6,21,1,
1013      &2,3,4,5,6,21,1,2,3,4,5,6,21,1,2,3,4,5,6,1,2,3,4,5,6,1,2,3,4,5,6,
1014      &21,3100111,3200111,21,22,23,-24,21,22,23,24,22,23,-24,23,24,1,2,
1015      &3,4,5,6,7,8,11,12,13,14,15,16,17,18,21,22,23,24,0,9*11,9*-11,
1016      &2*11,2*-11,9*13,9*-13,2*13,2*-13,9*15,9*-15,2*15,2*-15,1,2,3,4,5,
1017      &6,11,12,9900012,13,14,9900014,15,16,9900016,3*-1,3*-3,3*-5,-11,
1018      &-13,-15,3*-11,2*-13,-15,24,3*-11,2*-13,-15,9900024,4011*0/
1019       DATA (KFDP(I,2),I=   1, 339)/3*1,2,4,6,8,1,3*2,1,3,5,7,2,3*3,2,4,
1020      &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,6*1000006,3*7,
1021      &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,
1022      &13,11,13,-211,-213,-211,-213,-211,-213,-211,-213,2*-211,-321,
1023      &-323,-321,2*-323,3*-321,4*-211,-213,-211,-213,-211,-213,-211,
1024      &-213,-211,-213,3*-211,-213,4*-211,-323,-321,2*-211,2*-321,3*-211,
1025      &2*15,16,15,16,15,2*17,18,17,2*18,2*17,-1,-2,-3,-4,-5,-6,-7,-8,21,
1026      &-1,-2,-3,-4,-5,-6,-7,-8,-11,-13,-15,-17,-1,-2,-3,-4,-5,-6,-7,-8,
1027      &-11,-12,-13,-14,-15,-16,-17,-18,2,4,6,8,2,4,6,8,2,4,6,8,2,4,6,8,
1028      &12,14,16,18,-1,-2,-3,-4,-5,-6,-7,-8,-11,-13,-15,-17,21,22,2*23,
1029      &-24,2*1000022,1000023,1000022,1000023,1000025,1000022,1000023,
1030      &1000025,1000035,-1000024,-1000037,-1000024,-1000037,-1000001,
1031      &2*-2000001,2000001,-1000002,2*-2000002,2000002,-1000003,
1032      &2*-2000003,2000003,-1000004,2*-2000004,2000004,-1000005,
1033      &2*-2000005,2000005,-1000006,2*-2000006,2000006,-1000011,
1034      &2*-2000011,2000011,-1000012,2*-2000012,2000012,-1000013,
1035      &2*-2000013,2000013,-1000014,2*-2000014,2000014,-1000015,
1036      &2*-2000015,2000015,-1000016,2*-2000016,2000016,-1,-2,-3,-4,-5,-6,
1037      &-7,-8,-11,-12,-13,-14,-15,-16,-17,-18,-24,-37,22,25,2*36,2,4,6,8,
1038      &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/
1039       DATA (KFDP(I,2),I= 340, 533)/-7,-8,-11,-13,-15,-17,21,22,2*23,
1040      &-24,2*25,-37,-24,3*36,2*1000022,1000023,1000022,1000023,1000025,
1041      &1000022,1000023,1000025,1000035,-1000024,-1000037,-1000024,
1042      &-1000037,-1000001,2*-2000001,2000001,-1000002,2*-2000002,2000002,
1043      &-1000003,2*-2000003,2000003,-1000004,2*-2000004,2000004,-1000005,
1044      &2*-2000005,2000005,-1000006,2*-2000006,2000006,-1000011,
1045      &2*-2000011,2000011,-1000012,2*-2000012,2000012,-1000013,
1046      &2*-2000013,2000013,-1000014,2*-2000014,2000014,-1000015,
1047      &2*-2000015,2000015,-1000016,2*-2000016,2000016,-1,-2,-3,-4,-5,-6,
1048      &-7,-8,-11,-13,-15,-17,21,22,2*23,-24,2*25,-37,-24,2*1000022,
1049      &1000023,1000022,1000023,1000025,1000022,1000023,1000025,1000035,
1050      &-1000024,-1000037,-1000024,-1000037,-1000001,2*-2000001,2000001,
1051      &-1000002,2*-2000002,2000002,-1000003,2*-2000003,2000003,-1000004,
1052      &2*-2000004,2000004,-1000005,2*-2000005,2000005,-1000006,
1053      &2*-2000006,2000006,-1000011,2*-2000011,2000011,-1000012,
1054      &2*-2000012,2000012,-1000013,2*-2000013,2000013,-1000014,
1055      &2*-2000014,2000014,-1000015,2*-2000015,2000015,-1000016,
1056      &2*-2000016,2000016,2,4,6,8,12,14,16,18,25,1000024,1000037,
1057      &1000024,1000037,1000024,1000037,1000024,1000037,2*-1000005,
1058      &2*-2000005,1000002,1000004,1000012,1000014,2*1000016,-3,-4,-5,-6/
1059       DATA (KFDP(I,2),I= 534, 938)/-7,-8,-13,-15,-17,11,-82,12,14,-1,
1060      &-3,11,13,15,1,4,3,4,1,3,22,11,-211,2*22,-13,-11,-211,211,111,211,
1061      &-321,130,310,22,2*111,-211,11,-11,13,-13,-211,111,22,14,12,111,
1062      &22,111,3*211,-311,22,211,22,111,-211,211,11,-211,13,22,-211,111,
1063      &-211,22,111,-11,-211,111,2*-211,-321,130,310,221,111,-211,111,
1064      &2*0,-211,111,22,-211,111,-211,111,-211,211,-213,113,223,221,14,
1065      &111,211,111,-11,-13,211,111,22,211,111,211,111,2*211,213,113,223,
1066      &221,22,-211,111,113,223,22,111,-321,310,211,111,2*-211,221,22,
1067      &-11,-13,-211,-321,130,310,221,-211,111,11*12,11*14,2*211,2*213,
1068      &211,20213,2*321,2*323,211,213,211,213,211,213,211,213,211,213,
1069      &211,213,3*211,213,211,2*321,8*211,2*113,3*211,111,22,211,111,211,
1070      &111,4*211,8*12,8*14,2*211,2*213,2*111,221,2*113,223,333,20213,
1071      &211,2*321,323,2*311,313,-211,111,113,2*211,321,2*211,311,321,310,
1072      &211,-211,4*211,321,4*211,113,2*211,-321,111,22,-211,111,-211,111,
1073      &-211,211,-211,211,16,5*12,5*14,3*211,3*213,211,2*111,2*113,
1074      &2*-311,2*-313,-2112,3*321,323,2*-1,22,111,321,311,321,311,-82,
1075      &-11,-13,-82,22,-82,6*-11,6*-13,2*-15,211,213,20213,211,213,20213,
1076      &431,433,431,433,311,313,311,313,311,313,-1,-4,-3,-4,-1,-3,22,
1077      &-211,111,-211,111,-211,211,-211,211,6*-11,6*-13,2*-15,211,213,
1078      &20213,211,213,20213,431,433,431,433,321,323,321,323,321,323,-1/
1079       DATA (KFDP(I,2),I= 939,1352)/-4,-3,-4,-1,-3,22,211,111,211,111,
1080      &4*211,6*-11,6*-13,2*-15,211,213,20213,211,213,20213,431,433,431,
1081      &433,221,331,333,221,331,333,221,331,333,-1,-4,-3,-4,-1,-3,22,
1082      &-321,-311,-321,-311,-15,-3,-1,2*-11,2*-13,2*-15,-1,-4,-3,-4,-3,
1083      &-4,-1,-4,2*12,2*14,2,3,2,3,2*12,2*14,2,1,22,411,421,411,421,21,
1084      &-11,-13,-15,-1,-2,-3,-4,2*21,22,21,2*-211,111,22,111,211,22,211,
1085      &-211,11,2*-211,111,-211,111,22,11,22,111,-211,211,111,211,22,211,
1086      &111,211,-211,22,11,13,11,-211,2*111,2*22,111,211,-321,-211,111,
1087      &11,2*-211,7*12,7*14,-321,-323,-311,-313,-311,-313,211,213,211,
1088      &213,211,213,111,221,331,113,223,111,221,113,223,321,323,321,-211,
1089      &-213,111,221,331,113,223,333,10221,111,221,331,113,223,211,213,
1090      &211,213,321,323,321,323,321,323,311,313,311,313,2*-1,-3,-1,2203,
1091      &3201,3203,2203,2101,2103,12,14,-1,-3,2*111,2*211,12,14,-1,-3,22,
1092      &111,2*22,111,22,12,14,-1,-3,22,12,14,-1,-3,12,14,-1,-3,12,14,-1,
1093      &-3,12,14,-1,-3,12,14,-1,-3,12,14,-1,-3,12,14,-1,-3,2*-211,11,13,
1094      &15,-211,-213,-20213,-431,-433,3*3122,1,4,3,4,1,3,11,13,15,1,4,3,
1095      &4,1,3,11,13,15,1,4,3,4,1,3,2*111,2*211,11,13,15,1,4,3,4,1,3,11,
1096      &13,15,1,4,3,4,1,3,4*22,11,13,15,1,4,3,4,1,3,22,11,13,15,1,4,3,4,
1097      &1,3,11,13,15,1,4,3,4,1,3,11,13,15,1,4,3,4,1,3,11,13,15,1,4,3,4,1,
1098      &3,11,13,15,1,4,3,4,1,3,11,13,15,1,4,3,4,1,3,11,13,15,1,4,3,4,1,3/
1099       DATA (KFDP(I,2),I=1353,1822)/11,13,15,1,4,3,4,1,3,11,13,15,1,4,3,
1100      &4,1,3,11,13,15,1,4,3,4,1,3,11,13,15,1,4,3,4,1,3,11,13,15,1,4,3,4,
1101      &1,3,11,13,15,1,4,3,4,1,3,11,13,15,1,4,3,4,1,3,11,13,15,1,4,3,4,1,
1102      &3,11,13,15,1,4,3,4,1,3,11,13,15,1,4,3,4,1,3,11,13,15,1,4,3,4,1,3,
1103      &2*111,2*211,-211,111,-321,130,310,-211,111,211,-211,111,-213,113,
1104      &-211,111,223,211,111,213,113,211,111,223,-211,111,-321,130,310,
1105      &2*-211,-311,311,-321,321,211,111,211,111,-211,111,-211,111,311,
1106      &2*321,311,22,2*-82,-211,111,-211,111,211,111,211,111,-321,-311,
1107      &-321,-311,411,421,411,421,22,2*21,-211,2*211,111,-211,111,2*211,
1108      &111,-211,211,111,211,-321,2*-311,-321,22,-211,111,211,111,-311,
1109      &311,-321,321,211,111,-211,111,321,311,22,-82,-211,111,211,111,
1110      &-321,-311,411,421,22,21,-11,-13,-82,211,111,221,111,4*22,-11,-13,
1111      &-15,-1,-2,-3,-4,2*21,211,111,3*22,1,2*2,4*1,2*-24,2*-37,2*1,3,5,
1112      &1,3,5,1,3,5,1,2,3,4,5,6,1,2,3,4,5,6,1,2,3,4,5,6,2,2*1,4*2,2*24,
1113      &2*37,2,1,3,5,1,3,5,1,3,5,3,2*4,4*3,2*-24,2*-37,3,1,3,5,1,3,5,1,3,
1114      &5,1,2,3,4,5,6,1,2,3,4,5,6,1,2,3,4,5,6,4,2*3,4*4,2*24,2*37,4,1,3,
1115      &5,1,3,5,1,3,2*5,2*6,4*5,2*-24,2*-37,5,1,3,5,1,3,5,1,3,5,1,2,3,4,
1116      &5,6,1,2,3,4,5,6,1,2,3,4,5,2*6,2*5,4*6,2*24,2*37,6,4,-15,16,1,3,5,
1117      &1,3,5,1,3,5,11,2*12,4*11,2*-24,-37,13,15,11,15,11,13,11,13,15,11,
1118      &13,15,1,3,5,1,3,5,1,3,5,12,2*11,4*12,2*24,2*37,11,13,15,11,13,15/
1119       DATA (KFDP(I,2),I=1823,2288)/1,3,5,1,3,5,1,3,5,13,2*14,4*13,
1120      &2*-24,2*-37,13,15,11,15,11,13,11,13,15,11,13,15,1,3,5,1,3,5,1,3,
1121      &5,14,2*13,4*14,2*24,2*37,11,13,15,11,13,15,1,3,5,1,3,5,1,3,5,15,
1122      &2*16,4*15,2*-24,2*-37,13,15,11,15,11,13,11,13,15,11,13,15,1,3,5,
1123      &1,3,5,1,3,5,16,2*15,4*16,2*24,2*37,11,13,15,11,13,15,1,3,5,1,3,5,
1124      &1,3,5,21,-1,1,-1,1,-2,2,-2,2,-3,3,-3,3,-4,4,-4,4,-5,5,-5,5,-6,6,
1125      &-6,6,1,3,5,2,4,6,1,3,5,2,4,6,1,3,5,2,4,6,1,3,5,2,4,6,1,-1,3,-3,5,
1126      &-5,1,-1,3,-3,5,-5,22,23,25,35,36,-1,-3,-13,13,-13,13,-13,13,-15,
1127      &15,-15,15,-15,15,-11,11,-11,11,-11,11,-15,15,-15,15,-15,15,-11,
1128      &11,-11,11,-11,11,-13,13,-13,13,-13,13,-1,1,-2,2,-1,1,-2,2,-1,1,
1129      &-2,2,-3,3,-4,4,-3,3,-4,4,-3,3,-4,4,-5,5,-6,6,-5,5,-6,6,-5,5,-6,6,
1130      &-1,1,-2,2,-1,1,-2,2,-1,1,-2,2,-3,3,-4,4,-3,3,-4,4,-3,3,-4,4,-5,5,
1131      &-6,6,-5,5,-6,6,-5,5,-6,6,-1,1,-2,2,-1,1,-2,2,-1,1,-2,2,-3,3,-4,4,
1132      &-3,3,-4,4,-3,3,-4,4,-5,5,-6,6,-5,5,-6,6,-5,5,-6,6,22,23,25,35,36,
1133      &22,23,11,13,15,12,14,16,1,3,5,2,4,25,35,36,-24,24,11,-11,13,-13,
1134      &15,-15,1,-1,3,-3,-24,24,11,-11,13,-13,15,-15,1,-1,3,-3,-37,37,
1135      &-37,37,-1,1,-1,1,-2,2,-2,2,-3,3,-3,3,-4,4,-4,4,-5,5,-5,5,-6,6,-6,
1136      &6,-11,11,-11,11,-12,12,-12,12,-13,13,-13,13,-14,14,-14,14,-15,15,
1137      &-15,15,-16,16,-16,16,1,3,5,2,4,-13,13,-13,13,-13,13,-15,15,-15,
1138      &15,-15,15,-11,11,-11,11,-11,11,-15,15,-15,15,-15,15,-11,11,-11/
1139       DATA (KFDP(I,2),I=2289,2743)/11,-11,11,-13,13,-13,13,-13,13,-1,1,
1140      &-2,2,-1,1,-2,2,-1,1,-2,2,-3,3,-4,4,-3,3,-4,4,-3,3,-4,4,-5,5,-6,6,
1141      &-5,5,-6,6,-5,5,-6,6,-1,1,-2,2,-1,1,-2,2,-1,1,-2,2,-3,3,-4,4,-3,3,
1142      &-4,4,-3,3,-4,4,-5,5,-6,6,-5,5,-6,6,-5,5,-6,6,-1,1,-2,2,-1,1,-2,2,
1143      &-1,1,-2,2,-3,3,-4,4,-3,3,-4,4,-3,3,-4,4,-5,5,-6,6,-5,5,-6,6,-5,5,
1144      &-6,6,24,37,24,-11,-13,-15,-1,-3,24,-11,-13,-15,-1,-3,24,-11,-13,
1145      &-15,-1,-3,24,-11,-13,-15,-1,-3,4*37,2*-1,2*2,2*-3,2*4,2*-5,2*6,
1146      &2*-11,2*12,2*-13,2*14,2*-15,2*16,-1,-3,-13,14,2*-13,14,2*-13,14,
1147      &-13,-15,16,2*-15,16,2*-15,16,-15,6*-11,-15,16,2*-15,16,2*-15,16,
1148      &-15,6*-11,6*-13,-1,-2,-1,2,-1,-2,-1,2,-1,-2,-1,2,-3,-4,-3,4,-3,
1149      &-4,-3,4,-3,-4,-3,4,-5,-6,-5,6,-5,-6,-5,6,-5,-6,-5,6,-1,-2,-1,2,
1150      &-1,-2,-1,2,-1,-2,-1,2,-3,-4,-3,4,-3,-4,-3,4,-3,-4,-3,4,-5,-6,-5,
1151      &6,-5,-6,-5,6,-5,-6,-5,6,-1,-2,-1,2,-1,-2,-1,2,-1,-2,-1,2,-3,-4,
1152      &-3,4,-3,-4,-3,4,-3,-4,-3,4,-5,-6,-5,6,-5,-6,-5,6,-5,-6,-5,6,22,
1153      &23,25,35,36,22,23,11,13,15,12,14,16,1,3,5,2,4,25,35,36,22,23,11,
1154      &13,15,12,14,16,1,3,5,2,4,25,35,36,-24,24,11,-11,13,-13,15,-15,1,
1155      &-1,3,-3,-24,24,11,-11,13,-13,15,-15,1,-1,3,-3,-37,37,-37,37,-1,1,
1156      &-1,1,-2,2,-2,2,-3,3,-3,3,-4,4,-4,4,-5,5,-5,5,-6,6,-6,6,-11,11,
1157      &-11,11,-12,12,-12,12,-13,13,-13,13,-14,14,-14,14,-15,15,-15,15,
1158      &-16,16,-16,16,1,3,5,2,4,-13,13,-13,13,-13,13,-15,15,-15,15,-15/
1159       DATA (KFDP(I,2),I=2744,3191)/15,-11,11,-11,11,-11,11,-15,15,-15,
1160      &15,-15,15,-11,11,-11,11,-11,11,-13,13,-13,13,-13,13,-1,1,-2,2,-1,
1161      &1,-2,2,-1,1,-2,2,-3,3,-4,4,-3,3,-4,4,-3,3,-4,4,-5,5,-6,6,-5,5,-6,
1162      &6,-5,5,-6,6,-1,1,-2,2,-1,1,-2,2,-1,1,-2,2,-3,3,-4,4,-3,3,-4,4,-3,
1163      &3,-4,4,-5,5,-6,6,-5,5,-6,6,-5,5,-6,6,-1,1,-2,2,-1,1,-2,2,-1,1,-2,
1164      &2,-3,3,-4,4,-3,3,-4,4,-3,3,-4,4,-5,5,-6,6,-5,5,-6,6,-5,5,-6,6,22,
1165      &23,25,35,36,22,23,11,13,15,12,14,16,1,3,5,2,4,25,35,36,22,23,11,
1166      &13,15,12,14,16,1,3,5,2,4,25,35,36,22,23,11,13,15,12,14,16,1,3,5,
1167      &2,4,25,35,36,-24,24,11,-11,13,-13,15,-15,1,-1,3,-3,-24,24,11,-11,
1168      &13,-13,15,-15,1,-1,3,-3,-37,37,-37,37,-1,1,-1,1,-2,2,-2,2,-3,3,
1169      &-3,3,-4,4,-4,4,-5,5,-5,5,-6,6,-6,6,-11,11,-11,11,-12,12,-12,12,
1170      &-13,13,-13,13,-14,14,-14,14,-15,15,-15,15,-16,16,-16,16,1,3,5,2,
1171      &4,-13,13,-13,13,-13,13,-15,15,-15,15,-15,15,-11,11,-11,11,-11,11,
1172      &-15,15,-15,15,-15,15,-11,11,-11,11,-11,11,-13,13,-13,13,-13,13,
1173      &-1,1,-2,2,-1,1,-2,2,-1,1,-2,2,-3,3,-4,4,-3,3,-4,4,-3,3,-4,4,-5,5,
1174      &-6,6,-5,5,-6,6,-5,5,-6,6,-1,1,-2,2,-1,1,-2,2,-1,1,-2,2,-3,3,-4,4,
1175      &-3,3,-4,4,-3,3,-4,4,-5,5,-6,6,-5,5,-6,6,-5,5,-6,6,-1,1,-2,2,-1,1,
1176      &-2,2,-1,1,-2,2,-3,3,-4,4,-3,3,-4,4,-3,3,-4,4,-5,5,-6,6,-5,5,-6,6,
1177      &-5,5,-6,6,24,37,23,11,13,15,12,14,16,1,3,5,2,4,25,35,36,24,-11,
1178      &-13,-15,-1,-3,24,-11,-13,-15,-1,-3,24,-11,-13,-15,-1,-3,24,-11/
1179       DATA (KFDP(I,2),I=3192,3692)/-13,-15,-1,-3,4*37,2*-1,2*2,2*-3,
1180      &2*4,2*-5,2*6,2*-11,2*12,2*-13,2*14,2*-15,2*16,-1,-3,-13,14,2*-13,
1181      &14,2*-13,14,-13,-15,16,2*-15,16,2*-15,16,-15,-11,12,2*-11,12,
1182      &2*-11,12,-11,-15,16,2*-15,16,2*-15,16,-15,-11,12,2*-11,12,2*-11,
1183      &12,-11,-13,14,2*-13,14,2*-13,14,-13,-1,-2,-1,2,-1,-2,-1,2,-1,-2,
1184      &-1,2,-3,-4,-3,4,-3,-4,-3,4,-3,-4,-3,4,-5,-6,-5,6,-5,-6,-5,6,-5,
1185      &-6,-5,6,-1,-2,-1,2,-1,-2,-1,2,-1,-2,-1,2,-3,-4,-3,4,-3,-4,-3,4,
1186      &-3,-4,-3,4,-5,-6,-5,6,-5,-6,-5,6,-5,-6,-5,6,-1,-2,-1,2,-1,-2,-1,
1187      &2,-1,-2,-1,2,-3,-4,-3,4,-3,-4,-3,4,-3,-4,-3,4,-5,-6,-5,6,-5,-6,
1188      &-5,6,-5,-6,-5,6,1,2*2,4*1,23,25,35,36,2*-24,2*-37,2*1,3,5,1,3,5,
1189      &1,3,5,1,2,3,4,5,6,1,2,3,4,5,6,1,2,3,4,5,6,2,2*1,4*2,23,25,35,36,
1190      &2*24,2*37,2,1,3,5,1,3,5,1,3,5,3,2*4,4*3,23,25,35,36,2*-24,2*-37,
1191      &3,1,3,5,1,3,5,1,3,5,1,2,3,4,5,6,1,2,3,4,5,6,1,2,3,4,5,6,4,2*3,
1192      &4*4,23,25,35,36,2*24,2*37,4,1,3,5,1,3,5,1,3,2*5,2*6,4*5,23,25,35,
1193      &36,2*-24,2*-37,5,1,3,5,1,3,5,1,3,5,1,2,3,4,5,6,1,2,3,4,5,6,1,2,3,
1194      &4,5,2*6,2*5,4*6,23,25,35,36,2*24,2*37,6,1,3,5,1,3,5,1,3,5,11,
1195      &2*12,4*11,23,25,35,36,2*-24,2*-37,13,15,11,15,11,13,11,13,15,11,
1196      &13,15,1,3,5,1,3,5,1,3,5,13,2*14,4*13,23,25,35,36,2*-24,2*-37,13,
1197      &15,11,15,11,13,11,13,15,11,13,15,1,3,5,1,3,5,1,3,5,15,2*16,4*15,
1198      &23,25,35,36,2*-24,2*-37,13,15,11,15,11,13,11,13,15,11,13,15,1,3/
1199       DATA (KFDP(I,2),I=3693,8000)/5,1,3,5,1,3,5,-3,-4,-5,-6,-11,-13,
1200      &-15,21,-1,-3,2*-5,5,12,14,16,-3,-4,-5,-6,-11,-13,-15,21,-5,-6,21,
1201      &-24,-3000211,-24,-3000211,3000111,3000221,3000111,3000221,-1,-2,
1202      &-3,-4,-5,-6,-7,-8,-11,-12,-13,-14,-15,-16,-17,-18,23,3000111,23,
1203      &3000111,22,3000221,2,4,6,8,2,4,6,8,2,4,6,8,2,4,6,8,12,14,16,18,
1204      &2*3000111,2*3000221,-3000211,2*-24,-3000211,-1,-2,-3,-4,-5,-6,-7,
1205      &-8,-11,-12,-13,-14,-15,-16,-17,-18,-1,-2,-3,-4,-5,-6,-1,-2,-3,-4,
1206      &-5,-6,21,-1,-2,-3,-4,-5,-6,21,-1,-2,-3,-4,-5,-6,21,-1,-2,-3,-4,
1207      &-5,-6,-1,-2,-3,-4,-5,-6,-1,-2,-3,-4,-5,-6,3*21,3*1,4*2,1,2*11,
1208      &2*12,11,-1,-2,-3,-4,-5,-6,-7,-8,-11,-12,-13,-14,-15,-16,-17,-18,
1209      &21,22,23,-24,0,3*-1,3*-3,3*-5,3*1,3*3,3*5,2*-13,2*15,3*-1,3*-3,
1210      &3*-5,3*1,3*3,3*5,2*-11,2*15,3*-1,3*-3,3*-5,3*1,3*3,3*5,2*-11,
1211      &2*13,-1,-2,-3,-4,-5,-6,-11,-12,9900012,-13,-14,9900014,-15,-16,
1212      &9900016,2,4,6,2,4,6,2,4,6,9900012,9900014,9900016,-11,-13,-15,
1213      &-13,2*-15,24,-11,-13,-15,-13,2*-15,9900024,4011*0/
1214       DATA (KFDP(I,3),I=   1,1021)/81*0,14,6*0,2*16,2*0,6*111,310,130,
1215      &2*0,3*111,310,130,321,113,211,223,221,2*113,2*211,2*223,2*221,
1216      &2*113,221,2*113,2*213,-213,113,2*111,310,130,310,130,2*310,130,
1217      &402*0,4*3,4*4,1,4,3,2*2,0,-11,8*0,-211,5*0,2*111,211,-211,211,
1218      &-211,10*0,111,4*0,2*111,-211,-11,11,-13,22,111,3*0,22,3*0,111,
1219      &211,4*0,111,11*0,111,-211,6*0,-211,3*111,7*0,111,-211,5*0,2*221,
1220      &3*0,111,5*0,111,11*0,-311,-313,-311,-321,-313,-323,111,221,331,
1221      &113,223,-311,-313,-311,-321,-313,-323,111,221,331,113,223,22*0,
1222      &111,113,2*211,-211,-311,211,111,3*211,-211,7*211,7*0,111,-211,
1223      &111,-211,-321,-323,-311,-321,-313,-323,-211,-213,-321,-323,-311,
1224      &-321,-313,-323,-211,-213,22*0,111,113,-311,2*-211,211,-211,310,
1225      &-211,2*111,211,2*-211,-321,-211,2*211,-211,111,-211,2*211,6*0,
1226      &111,-211,111,-211,0,221,331,333,321,311,221,331,333,321,311,20*0,
1227      &3,13*0,-411,-413,-10413,-10411,-20413,-415,-411,-413,-10413,
1228      &-10411,-20413,-415,-411,-413,16*0,-4,-1,-4,-3,2*-2,5*0,111,-211,
1229      &111,-211,-421,-423,-10423,-10421,-20423,-425,-421,-423,-10423,
1230      &-10421,-20423,-425,-421,-423,16*0,-4,-1,-4,-3,2*-2,5*0,111,-211,
1231      &111,-211,-431,-433,-10433,-10431,-20433,-435,-431,-433,-10433,
1232      &-10431,-20433,-435,-431,-433,19*0,-4,-1,-4,-3,2*-2,8*0,441,443,
1233      &441,443,441,443,-4,-1,-4,-3,-4,-3,-4,-1,531,533,531,533,3,2,3,2/
1234       DATA (KFDP(I,3),I=1022,2197)/511,513,511,513,1,2,13*0,2*21,11*0,
1235      &2112,6*0,2212,12*0,2*3122,3212,10*0,3322,2*0,3122,3212,3214,2112,
1236      &2114,2212,2112,3122,3212,3214,2112,2114,2212,2112,52*0,3*3,1,6*0,
1237      &4*3,4*0,4*3,6*0,4*3,0,28*3,2*0,3*4122,8*0,4,1,4,3,2*2,4*4,1,4,3,
1238      &2*2,4*4,1,4,3,2*2,4*0,4*4,1,4,3,2*2,4*4,1,4,3,2*2,4*0,4*4,1,4,3,
1239      &2*2,0,4*4,1,4,3,2*2,4*4,1,4,3,2*2,4*4,1,4,3,2*2,4*4,1,4,3,2*2,
1240      &4*4,1,4,3,2*2,4*4,1,4,3,2*2,4*4,1,4,3,2*2,4*4,1,4,3,2*2,4*4,1,4,
1241      &3,2*2,4*4,1,4,3,2*2,4*4,1,4,3,2*2,4*4,1,4,3,2*2,4*4,1,4,3,2*2,
1242      &4*4,1,4,3,2*2,4*4,1,4,3,2*2,4*4,1,4,3,2*2,4*4,1,4,3,2*2,4*4,1,4,
1243      &3,2*2,31*0,211,111,45*0,-211,2*111,-211,3*111,-211,111,211,30*0,
1244      &-211,111,13*0,2*21,-211,111,175*0,2*5,207*0,-1,-3,-5,-2,-4,-6,-1,
1245      &-3,-5,-2,-4,-6,-1,-3,-5,-2,-4,-6,-1,-3,-5,-2,-4,-6,-2,2,-4,4,-6,
1246      &6,-2,2,-4,4,-6,6,5*0,11,12,11,-11,13,-13,15,-15,11,-11,13,-13,15,
1247      &-15,11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,11,-11,13,-13,15,
1248      &-15,11,-11,13,-13,15,-15,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,
1249      &3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,
1250      &3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,
1251      &5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,
1252      &5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,7*0,-11,-13,-15,-12,-14,-16,
1253      &-1,-3,-5,-2,-4,5*0,-12,12,-14,14,-16,16,-2,2,-4,4,2*0,-12,12,-14/
1254       DATA (KFDP(I,3),I=2198,2789)/14,-16,16,-2,2,-4,4,52*0,-1,-3,-5,
1255      &-2,-4,11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,11,-11,13,-13,15,
1256      &-15,11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,11,-11,13,-13,15,
1257      &-15,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,
1258      &1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,
1259      &1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,
1260      &3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,
1261      &3,-3,5,-5,5,-5,3*0,12,14,16,2,4,0,12,14,16,2,4,0,12,14,16,2,4,0,
1262      &12,14,16,2,4,28*0,2,4,12,-11,11,14,-13,13,16,-15,15,12,-11,11,14,
1263      &-13,13,16,-15,15,12,11,14,13,16,15,12,-11,11,14,-13,13,16,-15,15,
1264      &12,11,14,13,16,15,12,11,14,13,16,15,2*2,1,-1,2*4,3,-3,2*6,5,-5,
1265      &2*2,1,-1,2*4,3,-3,2*6,5,-5,2*2,1,-1,2*4,3,-3,2*6,5,-5,2*2,1,-1,
1266      &2*4,3,-3,2*6,5,-5,2*2,1,-1,2*4,3,-3,2*6,5,-5,2*2,1,-1,2*4,3,-3,
1267      &2*6,5,-5,2*2,1,-1,2*4,3,-3,2*6,5,-5,2*2,1,-1,2*4,3,-3,2*6,5,-5,
1268      &2*2,1,-1,2*4,3,-3,2*6,5,-5,7*0,-11,-13,-15,-12,-14,-16,-1,-3,-5,
1269      &-2,-4,5*0,-11,-13,-15,-12,-14,-16,-1,-3,-5,-2,-4,5*0,-12,12,-14,
1270      &14,-16,16,-2,2,-4,4,2*0,-12,12,-14,14,-16,16,-2,2,-4,4,52*0,-1,
1271      &-3,-5,-2,-4,11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,11,-11,13,
1272      &-13,15,-15,11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,11,-11,13,
1273      &-13,15,-15,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5/
1274       DATA (KFDP(I,3),I=2790,3335)/-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,
1275      &-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,
1276      &-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,
1277      &-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,7*0,-11,-13,
1278      &-15,-12,-14,-16,-1,-3,-5,-2,-4,5*0,-11,-13,-15,-12,-14,-16,-1,-3,
1279      &-5,-2,-4,5*0,-11,-13,-15,-12,-14,-16,-1,-3,-5,-2,-4,5*0,-12,12,
1280      &-14,14,-16,16,-2,2,-4,4,2*0,-12,12,-14,14,-16,16,-2,2,-4,4,52*0,
1281      &-1,-3,-5,-2,-4,11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,11,-11,
1282      &13,-13,15,-15,11,-11,13,-13,15,-15,11,-11,13,-13,15,-15,11,-11,
1283      &13,-13,15,-15,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,
1284      &5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,
1285      &5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,
1286      &1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,1,-1,3,-3,3,-3,5,-5,5,-5,1,-1,
1287      &1,-1,3,-3,3,-3,5,-5,5,-5,3*0,-11,-13,-15,-12,-14,-16,-1,-3,-5,-2,
1288      &-4,4*0,12,14,16,2,4,0,12,14,16,2,4,0,12,14,16,2,4,0,12,14,16,2,4,
1289      &28*0,2,4,12,-11,11,14,-13,13,16,-15,15,12,-11,11,14,-13,13,16,
1290      &-15,15,12,-11,11,14,-13,13,16,-15,15,12,-11,11,14,-13,13,16,-15,
1291      &15,12,-11,11,14,-13,13,16,-15,15,12,-11,11,14,-13,13,16,-15,15,
1292      &2*2,1,-1,2*4,3,-3,2*6,5,-5,2*2,1,-1,2*4,3,-3,2*6,5,-5,2*2,1,-1,
1293      &2*4,3,-3,2*6,5,-5,2*2,1,-1,2*4,3,-3,2*6,5,-5,2*2,1,-1,2*4,3,-3/
1294       DATA (KFDP(I,3),I=3336,8000)/2*6,5,-5,2*2,1,-1,2*4,3,-3,2*6,5,-5,
1295      &2*2,1,-1,2*4,3,-3,2*6,5,-5,2*2,1,-1,2*4,3,-3,2*6,5,-5,2*2,1,-1,
1296      &2*4,3,-3,2*6,5,-5,324*0,-5,170*0,2,4,6,2,4,6,2,4,6,-2,-4,-6,-2,
1297      &-4,-6,-2,-4,-6,2*9900014,2*9900016,2,4,6,2,4,6,2,4,6,-2,-4,-6,-2,
1298      &-4,-6,-2,-4,-6,2*9900012,2*9900016,2,4,6,2,4,6,2,4,6,-2,-4,-6,-2,
1299      &-4,-6,-2,-4,-6,2*9900012,2*9900014,4052*0/
1300       DATA (KFDP(I,4),I=   1,8000)/94*0,4*111,6*0,111,2*0,-211,0,-211,
1301      &3*0,111,2*-211,0,111,0,2*111,113,221,2*111,-213,-211,211,113,
1302      &6*111,310,2*130,402*0,13*81,41*0,-11,10*0,111,-211,4*0,111,62*0,
1303      &111,211,111,211,7*0,111,211,111,211,35*0,2*-211,2*111,211,111,
1304      &-211,2*211,2*-211,13*0,-211,111,-211,111,4*0,-211,111,-211,111,
1305      &34*0,111,-211,3*111,3*-211,2*111,3*-211,14*0,-321,-311,3*0,-321,
1306      &-311,20*0,-3,43*0,6*1,39*0,6*2,42*0,6*3,14*0,8*4,4*0,4*-5,4*0,
1307      &2*-5,67*0,-211,111,5*0,-211,111,52*0,2101,2103,2*2101,6*0,4*81,
1308      &4*0,4*81,6*0,4*81,0,28*81,13*0,6*2101,18*81,4*0,18*81,4*0,9*81,0,
1309      &162*81,31*0,-211,111,6516*0/
1310       DATA (KFDP(I,5),I=   1,8000)/96*0,2*111,17*0,111,7*0,2*111,0,
1311      &3*111,0,111,597*0,-211,2*111,-211,111,-211,111,65*0,111,-211,
1312      &3*111,-211,111,7193*0/
1313  
1314 C...PYDAT4, with particle names (character strings).
1315       DATA (CHAF(I,1),I=   1, 202)/'d','u','s','c','b','t','b''','t''', 
1316      &2*' ','e-','nu_e','mu-','nu_mu','tau-','nu_tau','tau''-',         
1317      &'nu''_tau',2*' ','g','gamma','Z0','W+','h0',6*' ','Z''0','Z"0',   
1318      &'W''+','H0','A0','H+',' ','Graviton',' ','R0','LQ_ue',38*' ',     
1319      &'specflav','rndmflav','phasespa','c-hadron','b-hadron',4*' ',     
1320      &'system','cluster','string','indep.','CMshower','SPHEaxis',       
1321      &'THRUaxis','CLUSjet','CELLjet','table',' ','reggeon','pi0',       
1322      &'rho0','a_20','K_L0','pi+','rho+','a_2+','eta','omega','f_2',     
1323      &'K_S0','K0','K*0','K*_20','K+','K*+','K*_2+','eta''','phi',       
1324      &'f''_2','D+','D*+','D*_2+','D0','D*0','D*_20','D_s+','D*_s+',     
1325      &'D*_2s+','eta_c','J/psi','chi_2c','B0','B*0','B*_20','B+','B*+',  
1326      &'B*_2+','B_s0','B*_s0','B*_2s0','B_c+','B*_c+','B*_2c+','eta_b',  
1327      &'Upsilon','chi_2b','pomeron','dd_1','Delta-','ud_0','ud_1','n0',  
1328      &'Delta0','uu_1','p+','Delta+','Delta++','sd_0','sd_1','Sigma-',   
1329      &'Sigma*-','Lambda0','su_0','su_1','Sigma0','Sigma*0','Sigma+',    
1330      &'Sigma*+','ss_1','Xi-','Xi*-','Xi0','Xi*0','Omega-','cd_0',       
1331      &'cd_1','Sigma_c0','Sigma*_c0','Lambda_c+','Xi_c0','cu_0','cu_1',  
1332      &'Sigma_c+','Sigma*_c+','Sigma_c++','Sigma*_c++','Xi_c+','cs_0',   
1333      &'cs_1','Xi''_c0','Xi*_c0','Xi''_c+','Xi*_c+','Omega_c0',          
1334      &'Omega*_c0','cc_1','Xi_cc+','Xi*_cc+','Xi_cc++','Xi*_cc++'/       
1335       DATA (CHAF(I,1),I= 203, 332)/'Omega_cc+','Omega*_cc+',            
1336      &'Omega*_ccc++','bd_0','bd_1','Sigma_b-','Sigma*_b-','Lambda_b0',  
1337      &'Xi_b-','Xi_bc0','bu_0','bu_1','Sigma_b0','Sigma*_b0','Sigma_b+', 
1338      &'Sigma*_b+','Xi_b0','Xi_bc+','bs_0','bs_1','Xi''_b-','Xi*_b-',    
1339      &'Xi''_b0','Xi*_b0','Omega_b-','Omega*_b-','Omega_bc0','bc_0',     
1340      &'bc_1','Xi''_bc0','Xi*_bc0','Xi''_bc+','Xi*_bc+','Omega''_bc0',   
1341      &'Omega*_bc0','Omega_bcc+','Omega*_bcc+','bb_1','Xi_bb-',          
1342      &'Xi*_bb-','Xi_bb0','Xi*_bb0','Omega_bb-','Omega*_bb-',            
1343      &'Omega_bbc0','Omega*_bbc0','Omega*_bbb-','a_00','b_10','a_0+',    
1344      &'b_1+','f_0','h_1','K*_00','K_10','K*_0+','K_1+','f''_0','h''_1', 
1345      &'D*_0+','D_1+','D*_00','D_10','D*_0s+','D_1s+','chi_0c','h_1c',   
1346      &'B*_00','B_10','B*_0+','B_1+','B*_0s0','B_1s0','B*_0c+','B_1c+',  
1347      &'chi_0b','h_1b','a_10','a_1+','f_1','K*_10','K*_1+','f''_1',      
1348      &'D*_1+','D*_10','D*_1s+','chi_1c','B*_10','B*_1+','B*_1s0',       
1349      &'B*_1c+','chi_1b','psi''','Upsilon''','~d_L','~u_L','~s_L',       
1350      &'~c_L','~b_1','~t_1','~e_L-','~nu_eL','~mu_L-','~nu_muL',         
1351      &'~tau_1-','~nu_tauL','~g','~chi_10','~chi_20','~chi_1+',          
1352      &'~chi_30','~chi_40','~chi_2+','~Gravitino','~d_R','~u_R','~s_R',  
1353      &'~c_R','~b_2','~t_2','~e_R-','~nu_eR','~mu_R-','~nu_muR',         
1354      &'~tau_2-','~nu_tauR','pi_tc0','pi_tc+','pi''_tc0','eta_tc0'/      
1355       DATA (CHAF(I,1),I= 333, 500)/'rho_tc0','rho_tc+','omega_tc',      
1356      &'V8_tc','pi_22_1_tc','pi_22_8_tc','rho_11_tc','rho_12_tc',        
1357      &'rho_21_tc','rho_22_tc','d*','u*','e*-','nu*_e0','Graviton*',     
1358      &'nu_Re','nu_Rmu','nu_Rtau','Z_R0','W_R+','H_L++','H_R++',         
1359      &'rho_diff0','pi_diffr+','omega_di','phi_diff','J/psi_di',         
1360      &'n_diffr0','p_diffr+',139*' '/                                    
1361       DATA (CHAF(I,2),I=   1, 205)/'dbar','ubar','sbar','cbar','bbar',  
1362      &'tbar','b''bar','t''bar',2*' ','e+','nu_ebar','mu+','nu_mubar',   
1363      &'tau+','nu_taubar','tau''+','nu''_taubar',5*' ','W-',9*' ',       
1364      &'W''-',2*' ','H-',3*' ','Rbar0','LQ_uebar',39*' ','rndmflavbar',  
1365      &' ','c-hadronbar','b-hadronbar',20*' ','pi-','rho-','a_2-',4*' ', 
1366      &'Kbar0','K*bar0','K*_2bar0','K-','K*-','K*_2-',3*' ','D-','D*-',  
1367      &'D*_2-','Dbar0','D*bar0','D*_2bar0','D_s-','D*_s-','D*_2s-',      
1368      &3*' ','Bbar0','B*bar0','B*_2bar0','B-','B*-','B*_2-','B_sbar0',   
1369      &'B*_sbar0','B*_2sbar0','B_c-','B*_c-','B*_2c-',4*' ','dd_1bar',   
1370      &'Deltabar+','ud_0bar','ud_1bar','nbar0','Deltabar0','uu_1bar',    
1371      &'pbar-','Deltabar-','Deltabar--','sd_0bar','sd_1bar','Sigmabar+', 
1372      &'Sigma*bar+','Lambdabar0','su_0bar','su_1bar','Sigmabar0',        
1373      &'Sigma*bar0','Sigmabar-','Sigma*bar-','ss_1bar','Xibar+',         
1374      &'Xi*bar+','Xibar0','Xi*bar0','Omegabar+','cd_0bar','cd_1bar',     
1375      &'Sigma_cbar0','Sigma*_cbar0','Lambda_cbar-','Xi_cbar0','cu_0bar', 
1376      &'cu_1bar','Sigma_cbar-','Sigma*_cbar-','Sigma_cbar--',            
1377      &'Sigma*_cbar--','Xi_cbar-','cs_0bar','cs_1bar','Xi''_cbar0',      
1378      &'Xi*_cbar0','Xi''_cbar-','Xi*_cbar-','Omega_cbar0',               
1379      &'Omega*_cbar0','cc_1bar','Xi_ccbar-','Xi*_ccbar-','Xi_ccbar--',   
1380      &'Xi*_ccbar--','Omega_ccbar-','Omega*_ccbar-','Omega*_cccbar-'/    
1381       DATA (CHAF(I,2),I= 206, 325)/'bd_0bar','bd_1bar','Sigma_bbar+',   
1382      &'Sigma*_bbar+','Lambda_bbar0','Xi_bbar+','Xi_bcbar0','bu_0bar',   
1383      &'bu_1bar','Sigma_bbar0','Sigma*_bbar0','Sigma_bbar-',             
1384      &'Sigma*_bbar-','Xi_bbar0','Xi_bcbar-','bs_0bar','bs_1bar',        
1385      &'Xi''_bbar+','Xi*_bbar+','Xi''_bbar0','Xi*_bbar0','Omega_bbar+',  
1386      &'Omega*_bbar+','Omega_bcbar0','bc_0bar','bc_1bar','Xi''_bcbar0',  
1387      &'Xi*_bcbar0','Xi''_bcbar-','Xi*_bcbar-','Omega''_bcba',           
1388      &'Omega*_bcbar0','Omega_bccbar-','Omega*_bccbar-','bb_1bar',       
1389      &'Xi_bbbar+','Xi*_bbbar+','Xi_bbbar0','Xi*_bbbar0','Omega_bbbar+', 
1390      &'Omega*_bbbar+','Omega_bbcbar0','Omega*_bbcbar0',                 
1391      &'Omega*_bbbbar+',2*' ','a_0-','b_1-',2*' ','K*_0bar0','K_1bar0',  
1392      &'K*_0-','K_1-',2*' ','D*_0-','D_1-','D*_0bar0','D_1bar0',         
1393      &'D*_0s-','D_1s-',2*' ','B*_0bar0','B_1bar0','B*_0-','B_1-',       
1394      &'B*_0sbar0','B_1sbar0','B*_0c-','B_1c-',3*' ','a_1-',' ',         
1395      &'K*_1bar0','K*_1-',' ','D*_1-','D*_1bar0','D*_1s-',' ',           
1396      &'B*_1bar0','B*_1-','B*_1sbar0','B*_1c-',3*' ','~d_Lbar',          
1397      &'~u_Lbar','~s_Lbar','~c_Lbar','~b_1bar','~t_1bar','~e_L+',        
1398      &'~nu_eLbar','~mu_L+','~nu_muLbar','~tau_1+','~nu_tauLbar',3*' ',  
1399      &'~chi_1-',2*' ','~chi_2-',' ','~d_Rbar','~u_Rbar','~s_Rbar',      
1400      &'~c_Rbar','~b_2bar','~t_2bar','~e_R+','~nu_eRbar','~mu_R+'/       
1401       DATA (CHAF(I,2),I= 326, 500)/'~nu_muRbar','~tau_2+',              
1402      &'~nu_tauRbar',' ','pi_tc-',3*' ','rho_tc-',8*' ','d*bar','u*bar', 
1403      &'e*bar+','nu*_ebar0',5*' ','W_R-','H_L--','H_R--',' ',            
1404      &'pi_diffr-',3*' ','n_diffrbar0','p_diffrbar-',139*' '/            
1405  
1406 C...PYDATR, with initial values for the random number generator.
1407       DATA MRPY/19780503,0,0,97,33,0/
1408  
1409 C...Default values for allowed processes and kinematics constraints.
1410       DATA MSEL/1/
1411       DATA MSUB/500*0/
1412       DATA ((KFIN(I,J),J=-40,40),I=1,2)/16*0,4*1,4*0,6*1,5*0,5*1,0,
1413      &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,
1414      &6*1,4*0,4*1,16*0/
1415       DATA CKIN/
1416      &  2.0D0, -1.0D0,  0.0D0, -1.0D0,  1.0D0,
1417      &  1.0D0,  -10D0,   10D0,  -40D0,   40D0,
1418      1  -40D0,   40D0,  -40D0,   40D0,  -40D0,
1419      1   40D0, -1.0D0,  1.0D0, -1.0D0,  1.0D0,
1420      2  0.0D0,  1.0D0,  0.0D0,  1.0D0, -1.0D0,
1421      2  1.0D0, -1.0D0,  1.0D0,    0D0,    0D0,
1422      3  2.0D0, -1.0D0,    0D0,    0D0,  0.0D0,
1423      3 -1.0D0,  0.0D0, -1.0D0,  4.0D0, -1.0D0,
1424      4 12.0D0, -1.0D0, 12.0D0, -1.0D0, 12.0D0,
1425      4 -1.0D0, 12.0D0, -1.0D0,    0D0,    0D0,
1426      5  0.0D0, -1.0D0,  0.0D0, -1.0D0,  0.0D0,
1427      5 -1.0D0,    0D0,    0D0,    0D0,    0D0,
1428      6 0.0001D0, 0.99D0, 0.0001D0, 0.99D0,    0D0,
1429      6   -1D0,    0D0,   -1D0,    0D0,   -1D0,
1430      7    0D0,   -1D0, 0.0001D0, 0.99D0, 0.0001D0,
1431      7 0.99D0,    2D0,   -1D0,    0D0,    0D0,
1432      8  120*0D0/
1433  
1434 C...Default values for main switches and parameters. Reset information.
1435       DATA (MSTP(I),I=1,100)/
1436      &  3,    1,    2,    0,    0,    0,    0,    0,    0,    0,
1437      1  1,    0,    1,   30,    0,    1,    4,    3,    4,    3,
1438      2  1,    0,    1,    0,    0,    0,    0,    0,    0,    1,
1439      3  1,    8,    0,    1,    0,    2,    1,    5,    2,    0,
1440      4  2,    1,    3,    7,    3,    1,    1,    0,    1,    0,
1441      5  7,    1,    3,    1,    5,    1,    1,    5,    1,    7,
1442      6  2,    3,    2,    2,    1,    5,    2,    1,    0,    0,
1443      7  1,    0,    0,    0,    0,    0,    0,    0,    0,    0,
1444      8  1,    1,  100,    0,    0,    2,    0,    0,    0,    0,
1445      9  1,    3,    1,    3,    0,    0,    0,    0,    0,    0/
1446       DATA (MSTP(I),I=101,200)/
1447      &  3,    1,    0,    0,    0,    0,    0,    0,    0,    0,
1448      1  1,    1,    1,    0,    0,    0,    0,    0,    0,    0,
1449      2  0,    1,    2,    1,    1,  100,    0,    0,   10,    0,
1450      3  0,    4,    0,    1,    0,    0,    0,    0,    0,    0,
1451      4  0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
1452      5  0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
1453      6  0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
1454      7  0,    2,    0,    0,    0,    0,    0,    0,    0,    0,
1455      8  6,  203, 2001,   11,   13,    0,    0,    0,    0,    0,
1456      9  0,    0,    0,    0,    0,    0,    0,    0,    0,    0/
1457       DATA (PARP(I),I=1,100)/
1458      &  0.25D0,  10D0, 8*0D0,
1459      1  0D0, 0D0, 1.0D0, 0.01D0, 0.5D0, 1.0D0, 1.0D0, 0.4D0, 2*0D0,
1460      2  10*0D0,
1461      3  1.5D0,2.0D0,0.075D0,1.0D0,0.2D0,0D0,2.0D0,0.70D0,0.006D0,0D0,
1462      4  0.02D0,2.0D0,0.10D0,1000D0,2054D0,123D0,246D0,50D0,0D0,0.054D0,
1463      5  10*0D0,
1464      6  0.25D0, 1.0D0,0.25D0, 1.0D0, 2.0D0,1D-3, 1.0D0,1D-3,2*0D0,
1465      7  4.0D0, 0.25D0, 8*0D0,
1466      8  1.90D0, 1.90D0, 0.5D0, 0.2D0, 0.33D0,
1467      8  0.66D0, 0.7D0, 0.5D0, 1000D0, 0.16D0,
1468      9  1.0D0,0.40D0,5.0D0,1.0D0,0D0,3.0D0,1.0D0,0.75D0,1.0D0,5.0D0/
1469       DATA (PARP(I),I=101,200)/
1470      &  0.5D0, 0.28D0,  1.0D0, 0.8D0, 0D0, 0D0, 0D0, 0D0, 0D0, 1D0,
1471      1  2.0D0, 3*0D0, 1.5D0, 0.5D0, 0.6D0, 2.5D0, 2.0D0, 1.0D0,
1472      2  1.0D0,  0.4D0, 8*0D0,
1473      3  0.01D0, 5*0D0, 200D0, 200D0, 0.333D0, 0.05D0,
1474      4  0.33333D0, 82D0, 1.33333D0, 4D0, 1D0,
1475      4  1D0,  .0182D0, 1D0, 0D0, 1.33333D0,
1476      5  0D0, 0D0, 0D0, 0D0, 0.3651480D0, 200D0, 0D0, 0D0, 0D0, 0D0,
1477      6  2.20D0, 23.6D0, 18.4D0, 11.5D0, 0.5D0, 0D0, 0D0, 0D0, 2*0D0,
1478      7  0D0,   0D0,   0D0,  1.0D0, 6*0D0,
1479      8  0.1D0, 0.01D0, 0.01D0, 0.01D0, 0.1D0, 0.01D0, 0.01D0, 0.01D0,
1480      8  0.3D0, 0.64D0,
1481      9  0.64D0, 5.0D0, 8*0D0/
1482       DATA MSTI/200*0/
1483       DATA PARI/200*0D0/
1484       DATA MINT/400*0/
1485       DATA VINT/400*0D0/
1486  
1487 C...Constants for the generation of the various processes.
1488       DATA (ISET(I),I=1,100)/
1489      &  1,    1,    1,   -1,    3,   -1,   -1,    3,   -2,    2,
1490      1  2,    2,    2,    2,    2,    2,   -1,    2,    2,    2,
1491      2 -1,    2,    2,    2,    2,    2,   -1,    2,    2,    2,
1492      3  2,    2,    2,    2,    2,    2,   -1,   -1,   -1,   -1,
1493      4 -1,   -1,   -1,   -1,   -1,   -1,   -1,   -1,   -1,   -1,
1494      5 -1,   -1,    2,    2,   -1,   -1,   -1,    2,   -1,   -1,
1495      6 -1,   -1,   -1,   -1,   -1,   -1,   -1,    2,    2,    2,
1496      7  4,    4,    4,   -1,   -1,    4,    4,   -1,   -1,    2,
1497      8  2,    2,    2,    2,    2,    2,    2,    2,    2,   -2,
1498      9  0,    0,    0,    0,    0,    9,   -2,   -2,    8,   -2/
1499       DATA (ISET(I),I=101,200)/
1500      & -1,    1,    1,    1,    1,    2,    2,    2,   -2,    2,
1501      1  2,    2,    2,    2,    2,   -1,   -1,   -1,   -2,   -2,
1502      2  5,    5,    5,    5,   -2,   -2,   -2,   -2,   -2,   -2,
1503      3  2,    2,    2,    2,    2,    2,    2,    2,    2,    2,
1504      4  1,    1,    1,    1,    1,    1,    1,    1,    1,   -2,
1505      5  1,    1,    1,   -2,   -2,    1,    1,    1,   -2,   -2,
1506      6  2,    2,    2,    2,    2,    2,    2,    2,    2,   -2,
1507      7  2,    2,    5,    5,   -2,    2,    2,    5,    5,   -2,
1508      8  5,    5,    2,    2,    2,    5,    5,    2,    2,    2,
1509      9  1,    1,    1,    2,    2,   -2,   -2,   -2,   -2,   -2/
1510       DATA (ISET(I),I=201,300)/
1511      &  2,    2,    2,    2,    2,    2,    2,    2,    2,    2,
1512      1  2,    2,    2,    2,   -2,    2,    2,    2,    2,    2,
1513      2  2,    2,    2,    2,    2,    2,    2,    2,    2,    2,
1514      3  2,    2,    2,    2,    2,    2,    2,    2,    2,    2,
1515      4  2,    2,    2,    2,   -1,    2,    2,    2,    2,    2,
1516      5  2,    2,    2,    2,   -1,    2,   -1,    2,    2,   -2,
1517      6  2,    2,    2,    2,    2,   -1,   -1,   -1,   -1,   -1,
1518      7  2,    2,    2,    2,    2,    2,    2,    2,    2,    2,
1519      8  2,    2,    2,    2,    2,    2,    2,    2,    2,    2,
1520      9  2,    2,    2,    2,    2,    2,    2,    2,    2,    2/
1521       DATA (ISET(I),I=301,500)/
1522      &  2,   39*-2,
1523      4  1,    1,    2,    2,    2,    2,    2,    2,    2,    2,
1524      5  5,    5,    1,    1,   -1,   -1,   -1,   -1,   -1,   -1,
1525      6  2,    2,    2,    2,    2,    2,    2,    2,   -1,    2,
1526      7  2,    2,    2,    2,    2,    2,    2,   -1,   -1,   -1,
1527      8  10*-2,
1528      9  1,    1,    2,    2,    2, 5*-2,
1529      &  100*-2/
1530       DATA ((KFPR(I,J),J=1,2),I=1,50)/
1531      &  23,    0,   24,    0,   25,    0,   24,    0,   25,    0,
1532      &  24,    0,   23,    0,   25,    0,    0,    0,    0,    0,
1533      1   0,    0,    0,    0,   21,   21,   21,   22,   21,   23,
1534      1  21,   24,   21,   25,   22,   22,   22,   23,   22,   24,
1535      2  22,   25,   23,   23,   23,   24,   23,   25,   24,   24,
1536      2  24,   25,   25,   25,    0,   21,    0,   22,    0,   23,
1537      3   0,   24,    0,   25,    0,   21,    0,   22,    0,   23,
1538      3   0,   24,    0,   25,    0,   21,    0,   22,    0,   23,
1539      4   0,   24,    0,   25,    0,   21,    0,   22,    0,   23,
1540      4   0,   24,    0,   25,    0,   21,    0,   22,    0,   23/
1541       DATA ((KFPR(I,J),J=1,2),I=51,100)/
1542      5   0,   24,    0,   25,    0,    0,    0,    0,    0,    0,
1543      5   0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
1544      6   0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
1545      6   0,    0,    0,    0,   21,   21,   24,   24,   23,   24,
1546      7  23,   23,   24,   24,   23,   24,   23,   25,   22,   22,
1547      7  23,   23,   24,   24,   24,   25,   25,   25,    0,  211,
1548      8   0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
1549      8 443,   21,10441,   21,20443,   21,  445,   21,    0,    0,
1550      9   0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
1551      9   0,    0,    0,    0,    0,    0,    0,    0,    0,    0/
1552       DATA ((KFPR(I,J),J=1,2),I=101,150)/
1553      &  23,    0,   25,    0,   25,    0,10441,    0,  445,    0,
1554      & 443,   22,  443,   21,  443,   22,    0,    0,   22,   25,
1555      1  21,   25,    0,   25,   21,   25,   22,   22,   21,   22,
1556      1  22,   23,   23,   23,   24,   24,    0,    0,    0,    0,
1557      2  25,    6,   25,    6,   25,    0,   25,    0,    0,    0,
1558      2   0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
1559      3   0,   21,    0,   21,    0,   22,    0,   22,    0,    0,
1560      3   0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
1561      4  32,    0,   34,    0,   37,    0,   41,    0,   42,    0,
1562      4 4000011, 0, 4000001, 0, 4000002, 0, 3000331, 0,   0,    0/
1563       DATA ((KFPR(I,J),J=1,2),I=151,200)/
1564      5  35,    0,   35,    0,   35,    0,    0,    0,    0,    0,
1565      5  36,    0,   36,    0,   36,    0,    0,    0,    0,    0,
1566      6   6,   37,   42,    0,   42,   42,   42,   42,   11,    0,
1567      6  11,    0, 0, 4000001, 0, 4000002, 0, 4000011,    0,    0,
1568      7  23,   35,   24,   35,   35,    0,   35,    0,    0,    0,
1569      7  23,   36,   24,   36,   36,    0,   36,    0,    0,    0,
1570      8  35,    6,   35,    6,   21,   35,    0,   35,   21,   35,
1571      8  36,    6,   36,    6,   21,   36,    0,   36,   21,   36,
1572      9  3000113, 0, 3000213, 0, 3000223, 0, 11,    0,   11,    0,
1573      9   0,    0,    0,    0,    0,    0,    0,    0,    0,    0/
1574       DATA ((KFPR(I,J),J=1,2),I=201,250)/
1575      &  1000011,   1000011,   2000011,   2000011,   1000011,
1576      &  2000011,   1000013,   1000013,   2000013,   2000013,
1577      &  1000013,   2000013,   1000015,   1000015,   2000015,
1578      &  2000015,   1000015,   2000015,   1000011,   1000012,
1579      1  1000015,   1000016,   2000015,   1000016,   1000012,
1580      1  1000012,   1000016,   1000016,         0,         0,
1581      1  1000022,   1000022,   1000023,   1000023,   1000025,
1582      1  1000025,   1000035,   1000035,   1000022,   1000023,
1583      2  1000022,   1000025,   1000022,   1000035,   1000023,
1584      2  1000025,   1000023,   1000035,   1000025,   1000035,
1585      2  1000024,   1000024,   1000037,   1000037,   1000024,
1586      2  1000037,   1000022,   1000024,   1000023,   1000024,
1587      3  1000025,   1000024,   1000035,   1000024,   1000022,
1588      3  1000037,   1000023,   1000037,   1000025,   1000037,
1589      3  1000035,   1000037,   1000021,   1000022,   1000021,
1590      3  1000023,   1000021,   1000025,   1000021,   1000035,
1591      4  1000021,   1000024,   1000021,   1000037,   1000021,
1592      4  1000021,   1000021,   1000021,         0,         0,
1593      4  1000002,   1000022,   2000002,   1000022,   1000002,
1594      4  1000023,   2000002,   1000023,   1000002,   1000025/
1595       DATA ((KFPR(I,J),J=1,2),I=251,300)/
1596      5  2000002,   1000025,   1000002,   1000035,   2000002,
1597      5  1000035,   1000001,   1000024,   2000005,   1000024,
1598      5  1000001,   1000037,   2000005,   1000037,   1000002,
1599      5  1000021,   2000002,   1000021,         0,         0,
1600      6  1000006,   1000006,   2000006,   2000006,   1000006,
1601      6  2000006,   1000006,   1000006,   2000006,   2000006,
1602      6        0,         0,         0,         0,         0,
1603      6        0,         0,         0,         0,         0,
1604      7  1000002,   1000002,   2000002,   2000002,   1000002,
1605      7  2000002,   1000002,   1000002,   2000002,   2000002,
1606      7  1000002,   2000002,   1000002,   1000002,   2000002,
1607      7  2000002,   1000002,   1000002,   2000002,   2000002,
1608      8  1000005,   1000002,   2000005,   2000002,   1000005,
1609      8  2000002,   1000005,   1000002,   2000005,   2000002,
1610      8  1000005,   2000002,   1000005,   1000005,   2000005,
1611      8  2000005,   1000005,   1000005,   2000005,   2000005,
1612      9  1000005,   1000005,   2000005,   2000005,   1000005,
1613      9  2000005,   1000005,   1000021,   2000005,   1000021,
1614      9  1000005,   2000005,        37,        25,        37,
1615      9       35,        36,        25,        36,        35/
1616       DATA ((KFPR(I,J),J=1,2),I=301,500)/
1617      &       37,        37,      78*0,
1618      4  9900041,         0,   9900042,         0,   9900041,
1619      4       11,   9900042,        11,   9900041,        13,
1620      4  9900042,        13,   9900041,        15,   9900042,
1621      4       15,   9900041,   9900041,   9900042,   9900042,
1622      5  9900041,         0,   9900042,         0,   9900023,
1623      5        0,   9900024,         0,         0,         0,
1624      5        0,         0,         0,         0,         0,
1625      5        0,         0,         0,         0,         0,
1626      6       24,        24,        24,   3000211,   3000211,
1627      6  3000211,        22,   3000111,        22,   3000221,
1628      6       23,   3000111,        23,   3000221,        24,
1629      6  3000211,         0,         0,        24,        23,
1630      7       24,   3000111,   3000211,        23,   3000211,
1631      7  3000111,        22,   3000211,        23,   3000211,
1632      7       24,   3000111,        24,   3000221,         0,
1633      7        0,         0,         0,         0,         0,
1634      8     20*0,
1635      9  5000039,         0,   5000039,         0,        21,
1636      9  5000039,         0,   5000039,        21,   5000039,
1637      9     10*0,
1638      &    200*0/
1639       DATA COEF/10000*0D0/
1640       DATA (((ICOL(I,J,K),K=1,2),J=1,4),I=1,40)/
1641      &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,
1642      &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,
1643      &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,
1644      &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,
1645      &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,
1646      &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,
1647      &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,
1648      &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,
1649      &4,4,3,3,2,2,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
1650      &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/
1651  
1652 C...Treatment of resonances.
1653       DATA (MWID(I)  ,I=   1, 500)/5*0,3*1,8*0,1,5*0,3*1,6*0,1,0,4*1,
1654      &3*0,2*1,254*0,19*2,0,7*2,0,2,0,2,0,26*1,146*0/
1655  
1656 C...Character constants: name of processes.
1657       DATA PROC(0)/                    'All included subprocesses   '/
1658       DATA (PROC(I),I=1,20)/
1659      &'f + fbar -> gamma*/Z0       ',  'f + fbar'' -> W+/-           ',
1660      &'f + fbar -> h0              ',  'gamma + W+/- -> W+/-        ',
1661      &'Z0 + Z0 -> h0               ',  'Z0 + W+/- -> W+/-           ',
1662      &'                            ',  'W+ + W- -> h0               ',
1663      &'                            ',  'f + f'' -> f + f'' (QFD)      ',
1664      1'f + f'' -> f + f'' (QCD)      ','f + fbar -> f'' + fbar''      ',
1665      1'f + fbar -> g + g           ',  'f + fbar -> g + gamma       ',
1666      1'f + fbar -> g + Z0          ',  'f + fbar'' -> g + W+/-       ',
1667      1'f + fbar -> g + h0          ',  'f + fbar -> gamma + gamma   ',
1668      1'f + fbar -> gamma + Z0      ',  'f + fbar'' -> gamma + W+/-   '/
1669       DATA (PROC(I),I=21,40)/
1670      2'f + fbar -> gamma + h0      ',  'f + fbar -> Z0 + Z0         ',
1671      2'f + fbar'' -> Z0 + W+/-      ', 'f + fbar -> Z0 + h0         ',
1672      2'f + fbar -> W+ + W-         ',  'f + fbar'' -> W+/- + h0      ',
1673      2'f + fbar -> h0 + h0         ',  'f + g -> f + g              ',
1674      2'f + g -> f + gamma          ',  'f + g -> f + Z0             ',
1675      3'f + g -> f'' + W+/-          ', 'f + g -> f + h0             ',
1676      3'f + gamma -> f + g          ',  'f + gamma -> f + gamma      ',
1677      3'f + gamma -> f + Z0         ',  'f + gamma -> f'' + W+/-      ',
1678      3'f + gamma -> f + h0         ',  'f + Z0 -> f + g             ',
1679      3'f + Z0 -> f + gamma         ',  'f + Z0 -> f + Z0            '/
1680       DATA (PROC(I),I=41,60)/
1681      4'f + Z0 -> f'' + W+/-         ', 'f + Z0 -> f + h0            ',
1682      4'f + W+/- -> f'' + g          ', 'f + W+/- -> f'' + gamma      ',
1683      4'f + W+/- -> f'' + Z0         ', 'f + W+/- -> f'' + W+/-       ',
1684      4'f + W+/- -> f'' + h0         ', 'f + h0 -> f + g             ',
1685      4'f + h0 -> f + gamma         ',  'f + h0 -> f + Z0            ',
1686      5'f + h0 -> f'' + W+/-         ', 'f + h0 -> f + h0            ',
1687      5'g + g -> f + fbar           ',  'g + gamma -> f + fbar       ',
1688      5'g + Z0 -> f + fbar          ',  'g + W+/- -> f + fbar''       ',
1689      5'g + h0 -> f + fbar          ',  'gamma + gamma -> f + fbar   ',
1690      5'gamma + Z0 -> f + fbar      ',  'gamma + W+/- -> f + fbar''   '/
1691       DATA (PROC(I),I=61,80)/
1692      6'gamma + h0 -> f + fbar      ',  'Z0 + Z0 -> f + fbar         ',
1693      6'Z0 + W+/- -> f + fbar''      ', 'Z0 + h0 -> f + fbar         ',
1694      6'W+ + W- -> f + fbar         ',  'W+/- + h0 -> f + fbar''      ',
1695      6'h0 + h0 -> f + fbar         ',  'g + g -> g + g              ',
1696      6'gamma + gamma -> W+ + W-    ',  'gamma + W+/- -> Z0 + W+/-   ',
1697      7'Z0 + Z0 -> Z0 + Z0          ',  'Z0 + Z0 -> W+ + W-          ',
1698      7'Z0 + W+/- -> Z0 + W+/-      ',  'Z0 + Z0 -> Z0 + h0          ',
1699      7'W+ + W- -> gamma + gamma    ',  'W+ + W- -> Z0 + Z0          ',
1700      7'W+/- + W+/- -> W+/- + W+/-  ',  'W+/- + h0 -> W+/- + h0      ',
1701      7'h0 + h0 -> h0 + h0          ',  'q + gamma -> q'' + pi+/-     '/
1702       DATA (PROC(I),I=81,100)/
1703      8'q + qbar -> Q + Qbar, mass  ',  'g + g -> Q + Qbar, massive  ',
1704      8'f + q -> f'' + Q, massive    ', 'g + gamma -> Q + Qbar, mass ',
1705      8'gamma + gamma -> F + Fbar, m',  'g + g -> J/Psi + g          ',
1706      8'g + g -> chi_0c + g         ',  'g + g -> chi_1c + g         ',
1707      8'g + g -> chi_2c + g         ',  '                            ',
1708      9'Elastic scattering          ',  'Single diffractive (XB)     ',
1709      9'Single diffractive (AX)     ',  'Double  diffractive         ',
1710      9'Low-pT scattering           ',  'Semihard QCD 2 -> 2         ',
1711      9'                            ',  '                            ',
1712      9'q + gamma* -> q             ',  '                            '/
1713       DATA (PROC(I),I=101,120)/
1714      &'g + g -> gamma*/Z0          ',  'g + g -> h0                 ',
1715      &'gamma + gamma -> h0         ',  'g + g -> chi_0c             ',
1716      &'g + g -> chi_2c             ',  'g + g -> J/Psi + gamma      ',
1717      &'gamma + g -> J/Psi + g      ',  'gamma+gamma -> J/Psi + gamma',
1718      &'                            ',  'f + fbar -> gamma + h0      ',
1719      1'q + qbar -> g + h0          ',  'q + g -> q + h0             ',
1720      1'g + g -> g + h0             ',  'g + g -> gamma + gamma      ',
1721      1'g + g -> g + gamma          ',  'g + g -> gamma + Z0         ',
1722      1'g + g -> Z0 + Z0            ',  'g + g -> W+ + W-            ',
1723      1'                            ',  '                            '/
1724       DATA (PROC(I),I=121,140)/
1725      2'g + g -> Q + Qbar + h0      ',  'q + qbar -> Q + Qbar + h0   ',
1726      2'f + f'' -> f + f'' + h0       ',
1727      2'f + f'' -> f" + f"'' + h0     ',
1728      2'                            ',  '                            ',
1729      2'                            ',  '                            ',
1730      2'                            ',  '                            ',
1731      3'f + gamma*_T -> f + g       ',  'f + gamma*_L -> f + g       ',
1732      3'f + gamma*_T -> f + gamma   ',  'f + gamma*_L -> f + gamma   ',
1733      3'g + gamma*_T -> f + fbar    ',  'g + gamma*_L -> f + fbar    ',
1734      3'gamma*_T+gamma*_T -> f+fbar ',  'gamma*_T+gamma*_L -> f+fbar ',
1735      3'gamma*_L+gamma*_T -> f+fbar ',  'gamma*_L+gamma*_L -> f+fbar '/
1736       DATA (PROC(I),I=141,160)/
1737      4'f + fbar -> gamma*/Z0/Z''0   ', 'f + fbar'' -> W''+/-          ',
1738      4'f + fbar'' -> H+/-           ', 'f + fbar'' -> R              ',
1739      4'q + l -> LQ                 ',  'e + gamma -> e*             ',
1740      4'd + g -> d*                 ',  'u + g -> u*                 ',
1741      4'g + g -> eta_tc             ',  '                            ',
1742      5'f + fbar -> H0              ',  'g + g -> H0                 ',
1743      5'gamma + gamma -> H0         ',  '                            ',
1744      5'                            ',  'f + fbar -> A0              ',
1745      5'g + g -> A0                 ',  'gamma + gamma -> A0         ',
1746      5'                            ',  '                            '/
1747       DATA (PROC(I),I=161,180)/
1748      6'f + g -> f'' + H+/-          ', 'q + g -> LQ + lbar          ',
1749      6'g + g -> LQ + LQbar         ',  'q + qbar -> LQ + LQbar      ',
1750      6'f + fbar -> f'' + fbar'' (g/Z)',
1751      6'f +fbar'' -> f" + fbar"'' (W) ',
1752      6'q + q'' -> q" + d*           ',  'q + q'' -> q" + u*           ',
1753      6'q + qbar -> e + e*          ',  '                            ',
1754      7'f + fbar -> Z0 + H0         ', 'f + fbar'' -> W+/- + H0      ',
1755      7'f + f'' -> f + f'' + H0       ',
1756      7'f + f'' -> f" + f"'' + H0     ',
1757      7'                            ',  'f + fbar -> Z0 + A0         ',
1758      7'f + fbar'' -> W+/- + A0      ',
1759      7'f + f'' -> f + f'' + A0       ',
1760      7'f + f'' -> f" + f"'' + A0     ',
1761      7'                            '/
1762       DATA (PROC(I),I=181,200)/
1763      8'g + g -> Q + Qbar + H0      ',  'q + qbar -> Q + Qbar + H0   ',
1764      8'q + qbar -> g + H0          ',  'q + g -> q + H0             ',
1765      8'g + g -> g + H0             ',  'g + g -> Q + Qbar + A0      ',
1766      8'q + qbar -> Q + Qbar + A0   ',  'q + qbar -> g + A0          ',
1767      8'q + g -> q + A0             ',  'g + g -> g + A0             ',
1768      9'f + fbar -> rho_tc0         ',  'f + f'' -> rho_tc+/-         ',
1769      9'f + fbar -> omega_tc0      ',  'f+fbar -> f''+fbar'' (ETC)  ',
1770      9'f+fbar'' -> f"+fbar"'' (ETC)','                          ',
1771      9'                            ',  '                            ',
1772      9'                            ',  '                            '/
1773       DATA (PROC(I),I=201,220)/
1774      &'f + fbar -> ~e_L + ~e_Lbar  ',  'f + fbar -> ~e_R + ~e_Rbar  ',
1775      &'f + fbar -> ~e_R + ~e_Lbar  ',  'f + fbar -> ~mu_L + ~mu_Lbar',
1776      &'f + fbar -> ~mu_R + ~mu_Rbar',  'f + fbar -> ~mu_L + ~mu_Rbar',
1777      &'f+fbar -> ~tau_1 + ~tau_1bar',  'f+fbar -> ~tau_2 + ~tau_2bar',
1778      &'f+fbar -> ~tau_1 + ~tau_2bar',  'q + qbar'' -> ~l_L + ~nulbar ',
1779      1'q+qbar''-> ~tau_1 + ~nutaubar', 'q+qbar''-> ~tau_2 + ~nutaubar',
1780      1'f + fbar -> ~nul + ~nulbar  ',  'f+fbar -> ~nutau + ~nutaubar',
1781      1'                            ',  'f + fbar -> ~chi1 + ~chi1   ',
1782      1'f + fbar -> ~chi2 + ~chi2   ',  'f + fbar -> ~chi3 + ~chi3   ',
1783      1'f + fbar -> ~chi4 + ~chi4   ',  'f + fbar -> ~chi1 + ~chi2   '/
1784       DATA (PROC(I),I=221,240)/
1785      2'f + fbar -> ~chi1 + ~chi3   ',  'f + fbar -> ~chi1 + ~chi4   ',
1786      2'f + fbar -> ~chi2 + ~chi3   ',  'f + fbar -> ~chi2 + ~chi4   ',
1787      2'f + fbar -> ~chi3 + ~chi4   ',  'f+fbar -> ~chi+-1 + ~chi-+1 ',
1788      2'f+fbar -> ~chi+-2 + ~chi-+2 ',  'f+fbar -> ~chi+-1 + ~chi-+2 ',
1789      2'q + qbar'' -> ~chi1 + ~chi+-1', 'q + qbar'' -> ~chi2 + ~chi+-1',
1790      3'q + qbar'' -> ~chi3 + ~chi+-1', 'q + qbar'' -> ~chi4 + ~chi+-1',
1791      3'q + qbar'' -> ~chi1 + ~chi+-2', 'q + qbar'' -> ~chi2 + ~chi+-2',
1792      3'q + qbar'' -> ~chi3 + ~chi+-2', 'q + qbar'' -> ~chi4 + ~chi+-2',
1793      3'q + qbar -> ~chi1 + ~g      ',  'q + qbar -> ~chi2 + ~g      ',
1794      3'q + qbar -> ~chi3 + ~g      ',  'q + qbar -> ~chi4 + ~g      '/
1795       DATA (PROC(I),I=241,260)/
1796      4'q + qbar'' -> ~chi+-1 + ~g   ', 'q + qbar'' -> ~chi+-2 + ~g  ',
1797      4'q + qbar -> ~g + ~g         ',  'g + g -> ~g + ~g            ',
1798      4'                            ',  'qj + g -> ~qj_L + ~chi1     ',
1799      4'qj + g -> ~qj_R + ~chi1     ',  'qj + g -> ~qj_L + ~chi2     ',
1800      4'qj + g -> ~qj_R + ~chi2     ',  'qj + g -> ~qj_L + ~chi3     ',
1801      5'qj + g -> ~qj_R + ~chi3     ',  'qj + g -> ~qj_L + ~chi4     ',
1802      5'qj + g -> ~qj_R + ~chi4     ',  'qj + g -> ~qk_L + ~chi+-1   ',
1803      5'qj + g -> ~qk_R + ~chi+-1   ',  'qj + g -> ~qk_L + ~chi+-2   ',
1804      5'qj + g -> ~qk_R + ~chi+-2   ',  'qj + g -> ~qj_L + ~g        ',
1805      5'qj + g -> ~qj_R + ~g        ',  '                            '/
1806       DATA (PROC(I),I=261,300)/
1807      6'f + fbar -> ~t_1 + ~t_1bar  ',  'f + fbar -> ~t_2 + ~t_2bar  ',
1808      6'f + fbar -> ~t_1 + ~t_2bar  ',  'g + g -> ~t_1 + ~t_1bar     ',
1809      6'g + g -> ~t_2 + ~t_2bar     ',  '                            ',
1810      6'                            ',  '                            ',
1811      6'                            ',  '                            ',
1812      7'qi + qj -> ~qi_L + ~qj_L    ',  'qi + qj -> ~qi_R + ~qj_R    ',
1813      7'qi + qj -> ~qi_L + ~qj_R    ',  'qi+qjbar -> ~qi_L + ~qj_Lbar',
1814      7'qi+qjbar -> ~qi_R + ~qj_Rbar',  'qi+qjbar -> ~qi_L + ~qj_Rbar',
1815      7'f + fbar -> ~qi_L + ~qi_Lbar',  'f + fbar -> ~qi_R + ~qi_Rbar',
1816      7'g + g -> ~qi_L + ~qi_Lbar   ',  'g + g -> ~qi_R + ~qi_Rbar   ',
1817      8'b + qj -> ~b_1 + ~qj_L      ',  'b + qj -> ~b_2 + ~qj_R      ',
1818      8'b + qj -> ~b_1 + ~qj_R      ',  'b + qjbar -> ~b_1 + ~qj_Lbar',
1819      8'b + qjbar -> ~b_2 + ~qj_Rbar',  'b + qjbar -> ~b_1 + ~qj_Rbar',
1820      8'f + fbar -> ~b_1 + ~b_1bar  ',  'f + fbar -> ~b_2 + ~b_2bar  ',
1821      8'g + g -> ~b_1 + ~b_1bar     ',  'g + g -> ~b_2 + ~b_2bar     ',
1822      9'b + b -> ~b_1 + ~b_1        ',  'b + b -> ~b_2 + ~b_2        ',
1823      9'b + b -> ~b_1 + ~b_2        ',  'b + g -> ~b_1 + ~g          ',
1824      9'b + g -> ~b_2 + ~g          ',  'b + bbar -> ~b_1 + ~b_2bar  ',
1825      9'f + fbar'' -> H+/- + h0     ',  'f + fbar -> H+/- + H0       ',
1826      9'f + fbar -> A0 + h0         ',  'f + fbar -> A0 + H0         '/
1827       DATA (PROC(I),I=301,340)/
1828      &'f + fbar -> H+ + H-         ', 39*'                          '/
1829       DATA (PROC(I),I=341,380)/
1830      4'l + l -> H_L++/--           ',  'l + l -> H_R++/--           ',
1831      4'l + gamma -> H_L++/-- e-/+  ',  'l + gamma -> H_R++/-- e-/+  ',
1832      4'l + gamma -> H_L++/-- mu-/+ ',  'l + gamma -> H_R++/-- mu-/+ ',
1833      4'l + gamma -> H_L++/-- tau-/+',  'l + gamma -> H_R++/-- tau-/+',
1834      4'f + fbar -> H_L++ + H_L--   ',  'f + fbar -> H_R++ + H_R--   ',
1835      5'f + f -> f'' + f'' + H_L++/-- ',
1836      5'f + f -> f'' + f'' + H_R++/-- ','f + fbar -> Z_R0            ',
1837      5'f + fbar'' -> W_R+/-         ',5*'                            ',
1838      6'                            ',  'f + fbar -> W_L+ W_L-       ',
1839      6'f + fbar -> W_L+/- pi_T-/+  ',  'f + fbar -> pi_T+ pi_T-     ',
1840      6'f + fbar -> gamma pi_T0     ',  'f + fbar -> gamma pi_T0''    ',
1841      6'f + fbar -> Z0 pi_T0        ',  'f + fbar -> Z0 pi_T0''       ',
1842      6'f + fbar -> W+/- pi_T-/+    ',  '                            ',
1843      7'f + fbar'' -> W_L+/- Z_L0    ', 'f + fbar'' -> W_L+/- pi_T0   ',
1844      7'f + fbar'' -> pi_T+/- Z_L0   ', 'f + fbar'' -> pi_T+/- pi_T0  ',
1845      7'f + fbar'' -> gamma pi_T+/-  ', 'f + fbar'' -> Z0 pi_T+/-     ',
1846      7'f + fbar'' -> W+/- pi_T0     ',
1847      7'f + fbar'' -> W+/- pi_T0''    ',
1848      7'                            ','                              ',
1849      7'                            '/
1850       DATA (PROC(I),I=381,500)/
1851      8 10* '                            ',
1852      9'f + fbar -> G*              ','g + g -> G*                   ',
1853      9'q + qbar -> g + G*          ','q + g -> q + G*               ',
1854      9'g + g -> g + G*             ','                              ',
1855      & 104*'                      '/
1856  
1857 C...Cross sections and slope offsets.
1858       DATA SIGT/294*0D0/
1859  
1860 C...Supersymmetry switches and parameters.
1861       DATA IMSS/0,
1862      &  0,  0,  0,  1,  0,  0,  0,  1,  0,  0,
1863      1  89*0/
1864       DATA RMSS/0D0,
1865      &  80D0,160D0,500D0,800D0,2D0,250D0,200D0,800D0,700D0,800D0,
1866      1  700D0,500D0,250D0,200D0,800D0,400D0,0D0,0.1D0,850D0,0.041D0,
1867      2   1D0,800D0,1D4,1D4,1D4,0D0,0D0,0D0,24D17,0D0,
1868      3  69*0D0/
1869 C...Initial values for R-violating SUSY couplings.
1870 C...Should not be changed here. See PYMSIN.
1871       DATA RVLAM/27*0D0/
1872       DATA RVLAMP/27*0D0/
1873       DATA RVLAMB/27*0D0/
1874  
1875 C...Data for histogramming routines.
1876       DATA IHIST/1000,20000,55,1/
1877       DATA INDX/1000*0/
1878  
1879       END
1880  
1881 C*********************************************************************
1882  
1883 C...PYTEST
1884 C...A simple program (disguised as subroutine) to run at installation
1885 C...as a check that the program works as intended.
1886  
1887       SUBROUTINE PYTEST(MTEST)
1888  
1889 C...Double precision and integer declarations.
1890       IMPLICIT DOUBLE PRECISION(A-H, O-Z)
1891       IMPLICIT INTEGER(I-N)
1892       INTEGER PYK,PYCHGE,PYCOMP
1893 C...Commonblocks.
1894       COMMON/PYJETS/N,NPAD,K(4000,5),P(4000,5),V(4000,5)
1895       COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
1896       COMMON/PYDAT2/KCHG(500,4),PMAS(500,4),PARF(2000),VCKM(4,4)
1897       COMMON/PYDAT3/MDCY(500,3),MDME(8000,2),BRAT(8000),KFDP(8000,5)
1898       COMMON/PYSUBS/MSEL,MSELPD,MSUB(500),KFIN(2,-40:40),CKIN(200)
1899       COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
1900       SAVE /PYJETS/,/PYDAT1/,/PYDAT2/,/PYDAT3/,/PYSUBS/,/PYPARS/
1901 C...Local arrays.
1902       DIMENSION PSUM(5),PINI(6),PFIN(6)
1903  
1904 C...Save defaults for values that are changed.
1905       MSTJ1=MSTJ(1)
1906       MSTJ3=MSTJ(3)
1907       MSTJ11=MSTJ(11)
1908       MSTJ42=MSTJ(42)
1909       MSTJ43=MSTJ(43)
1910       MSTJ44=MSTJ(44)
1911       PARJ17=PARJ(17)
1912       PARJ22=PARJ(22)
1913       PARJ43=PARJ(43)
1914       PARJ54=PARJ(54)
1915       MST101=MSTJ(101)
1916       MST104=MSTJ(104)
1917       MST105=MSTJ(105)
1918       MST107=MSTJ(107)
1919       MST116=MSTJ(116)
1920  
1921 C...First part: loop over simple events to be generated.
1922       IF(MTEST.GE.1) CALL PYTABU(20)
1923       NERR=0
1924       DO 180 IEV=1,500
1925  
1926 C...Reset parameter values. Switch on some nonstandard features.
1927         MSTJ(1)=1
1928         MSTJ(3)=0
1929         MSTJ(11)=1
1930         MSTJ(42)=2
1931         MSTJ(43)=4
1932         MSTJ(44)=2
1933         PARJ(17)=0.1D0
1934         PARJ(22)=1.5D0
1935         PARJ(43)=1D0
1936         PARJ(54)=-0.05D0
1937         MSTJ(101)=5
1938         MSTJ(104)=5
1939         MSTJ(105)=0
1940         MSTJ(107)=1
1941         IF(IEV.EQ.301.OR.IEV.EQ.351.OR.IEV.EQ.401) MSTJ(116)=3
1942  
1943 C...Ten events each for some single jets configurations.
1944         IF(IEV.LE.50) THEN
1945           ITY=(IEV+9)/10
1946           MSTJ(3)=-1
1947           IF(ITY.EQ.3.OR.ITY.EQ.4) MSTJ(11)=2
1948           IF(ITY.EQ.1) CALL PY1ENT(1,1,15D0,0D0,0D0)
1949           IF(ITY.EQ.2) CALL PY1ENT(1,3101,15D0,0D0,0D0)
1950           IF(ITY.EQ.3) CALL PY1ENT(1,-2203,15D0,0D0,0D0)
1951           IF(ITY.EQ.4) CALL PY1ENT(1,-4,30D0,0D0,0D0)
1952           IF(ITY.EQ.5) CALL PY1ENT(1,21,15D0,0D0,0D0)
1953  
1954 C...Ten events each for some simple jet systems; string fragmentation.
1955         ELSEIF(IEV.LE.130) THEN
1956           ITY=(IEV-41)/10
1957           IF(ITY.EQ.1) CALL PY2ENT(1,1,-1,40D0)
1958           IF(ITY.EQ.2) CALL PY2ENT(1,4,-4,30D0)
1959           IF(ITY.EQ.3) CALL PY2ENT(1,2,2103,100D0)
1960           IF(ITY.EQ.4) CALL PY2ENT(1,21,21,40D0)
1961           IF(ITY.EQ.5) CALL PY3ENT(1,2101,21,-3203,30D0,0.6D0,0.8D0)
1962           IF(ITY.EQ.6) CALL PY3ENT(1,5,21,-5,40D0,0.9D0,0.8D0)
1963           IF(ITY.EQ.7) CALL PY3ENT(1,21,21,21,60D0,0.7D0,0.5D0)
1964           IF(ITY.EQ.8) CALL PY4ENT(1,2,21,21,-2,40D0,
1965      &    0.4D0,0.64D0,0.6D0,0.12D0,0.2D0)
1966  
1967 C...Seventy events with independent fragmentation and momentum cons.
1968         ELSEIF(IEV.LE.200) THEN
1969           ITY=1+(IEV-131)/16
1970           MSTJ(2)=1+MOD(IEV-131,4)
1971           MSTJ(3)=1+MOD((IEV-131)/4,4)
1972           IF(ITY.EQ.1) CALL PY2ENT(1,4,-5,40D0)
1973           IF(ITY.EQ.2) CALL PY3ENT(1,3,21,-3,40D0,0.9D0,0.4D0)
1974           IF(ITY.EQ.3) CALL PY4ENT(1,2,21,21,-2,40D0,
1975      &    0.4D0,0.64D0,0.6D0,0.12D0,0.2D0)
1976           IF(ITY.GE.4) CALL PY4ENT(1,2,-3,3,-2,40D0,
1977      &    0.4D0,0.64D0,0.6D0,0.12D0,0.2D0)
1978  
1979 C...A hundred events with random jets (check invariant mass).
1980         ELSEIF(IEV.LE.300) THEN
1981   100     DO 110 J=1,5
1982             PSUM(J)=0D0
1983   110     CONTINUE
1984           NJET=2D0+6D0*PYR(0)
1985           DO 130 I=1,NJET
1986             KFL=21
1987             IF(I.EQ.1) KFL=INT(1D0+4D0*PYR(0))
1988             IF(I.EQ.NJET) KFL=-INT(1D0+4D0*PYR(0))
1989             EJET=5D0+20D0*PYR(0)
1990             THETA=ACOS(2D0*PYR(0)-1D0)
1991             PHI=6.2832D0*PYR(0)
1992             IF(I.LT.NJET) CALL PY1ENT(-I,KFL,EJET,THETA,PHI)
1993             IF(I.EQ.NJET) CALL PY1ENT(I,KFL,EJET,THETA,PHI)
1994             IF(I.EQ.1.OR.I.EQ.NJET) MSTJ(93)=1
1995             IF(I.EQ.1.OR.I.EQ.NJET) PSUM(5)=PSUM(5)+PYMASS(KFL)
1996             DO 120 J=1,4
1997               PSUM(J)=PSUM(J)+P(I,J)
1998   120       CONTINUE
1999   130     CONTINUE
2000           IF(PSUM(4)**2-PSUM(1)**2-PSUM(2)**2-PSUM(3)**2.LT.
2001      &    (PSUM(5)+PARJ(32))**2) GOTO 100
2002  
2003 C...Fifty e+e- continuum events with matrix elements.
2004         ELSEIF(IEV.LE.350) THEN
2005           MSTJ(101)=2
2006           CALL PYEEVT(0,40D0)
2007  
2008 C...Fifty e+e- continuum event with varying shower options.
2009         ELSEIF(IEV.LE.400) THEN
2010           MSTJ(42)=1+MOD(IEV,2)
2011           MSTJ(43)=1+MOD(IEV/2,4)
2012           MSTJ(44)=MOD(IEV/8,3)
2013           CALL PYEEVT(0,90D0)
2014  
2015 C...Fifty e+e- continuum events with coherent shower.
2016         ELSEIF(IEV.LE.450) THEN
2017           CALL PYEEVT(0,500D0)
2018  
2019 C...Fifty Upsilon decays to ggg or gammagg with coherent shower.
2020         ELSE
2021           CALL PYONIA(5,9.46D0)
2022         ENDIF
2023  
2024 C...Generate event. Find total momentum, energy and charge.
2025         DO 140 J=1,4
2026           PINI(J)=PYP(0,J)
2027   140   CONTINUE
2028         PINI(6)=PYP(0,6)
2029         CALL PYEXEC
2030         DO 150 J=1,4
2031           PFIN(J)=PYP(0,J)
2032   150   CONTINUE
2033         PFIN(6)=PYP(0,6)
2034  
2035 C...Check conservation of energy, momentum and charge;
2036 C...usually exact, but only approximate for single jets.
2037         MERR=0
2038         IF(IEV.LE.50) THEN
2039           IF((PFIN(1)-PINI(1))**2+(PFIN(2)-PINI(2))**2.GE.10D0)
2040      &    MERR=MERR+1
2041           EPZREM=PINI(4)+PINI(3)-PFIN(4)-PFIN(3)
2042           IF(EPZREM.LT.0D0.OR.EPZREM.GT.2D0*PARJ(31)) MERR=MERR+1
2043           IF(ABS(PFIN(6)-PINI(6)).GT.2.1D0) MERR=MERR+1
2044         ELSE
2045           DO 160 J=1,4
2046             IF(ABS(PFIN(J)-PINI(J)).GT.0.0001D0*PINI(4)) MERR=MERR+1
2047   160     CONTINUE
2048           IF(ABS(PFIN(6)-PINI(6)).GT.0.1D0) MERR=MERR+1
2049         ENDIF
2050         IF(MERR.NE.0) WRITE(MSTU(11),5000) (PINI(J),J=1,4),PINI(6),
2051      &  (PFIN(J),J=1,4),PFIN(6)
2052  
2053 C...Check that all KF codes are known ones, and that partons/particles
2054 C...satisfy energy-momentum-mass relation. Store particle statistics.
2055         DO 170 I=1,N
2056           IF(K(I,1).GT.20) GOTO 170
2057           IF(PYCOMP(K(I,2)).EQ.0) THEN
2058             WRITE(MSTU(11),5100) I
2059             MERR=MERR+1
2060           ENDIF
2061           PD=P(I,4)**2-P(I,1)**2-P(I,2)**2-P(I,3)**2-P(I,5)**2
2062           IF(ABS(PD).GT.MAX(0.1D0,0.001D0*P(I,4)**2).OR.P(I,4).LT.0D0)
2063      &    THEN
2064             WRITE(MSTU(11),5200) I
2065             MERR=MERR+1
2066           ENDIF
2067   170   CONTINUE
2068         IF(MTEST.GE.1) CALL PYTABU(21)
2069  
2070 C...List all erroneous events and some normal ones.
2071         IF(MERR.NE.0.OR.MSTU(24).NE.0.OR.MSTU(28).NE.0) THEN
2072           IF(MERR.GE.1) WRITE(MSTU(11),6400)
2073           CALL PYLIST(2)
2074         ELSEIF(MTEST.GE.1.AND.MOD(IEV-5,100).EQ.0) THEN
2075           CALL PYLIST(1)
2076         ENDIF
2077  
2078 C...Stop execution if too many errors.
2079         IF(MERR.NE.0) NERR=NERR+1
2080         IF(NERR.GE.10) THEN
2081           WRITE(MSTU(11),6300)
2082           CALL PYLIST(1)
2083           STOP
2084         ENDIF
2085   180 CONTINUE
2086  
2087 C...Summarize result of run.
2088       IF(MTEST.GE.1) CALL PYTABU(22)
2089  
2090 C...Reset commonblock variables changed during run.
2091       MSTJ(1)=MSTJ1
2092       MSTJ(3)=MSTJ3
2093       MSTJ(11)=MSTJ11
2094       MSTJ(42)=MSTJ42
2095       MSTJ(43)=MSTJ43
2096       MSTJ(44)=MSTJ44
2097       PARJ(17)=PARJ17
2098       PARJ(22)=PARJ22
2099       PARJ(43)=PARJ43
2100       PARJ(54)=PARJ54
2101       MSTJ(101)=MST101
2102       MSTJ(104)=MST104
2103       MSTJ(105)=MST105
2104       MSTJ(107)=MST107
2105       MSTJ(116)=MST116
2106  
2107 C...Second part: complete events of various kinds.
2108 C...Common initial values. Loop over initiating conditions.
2109       MSTP(122)=MAX(0,MIN(2,MTEST))
2110       MDCY(PYCOMP(111),1)=0
2111       DO 230 IPROC=1,8
2112  
2113 C...Reset process type, kinematics cuts, and the flags used.
2114         MSEL=0
2115         DO 190 ISUB=1,500
2116           MSUB(ISUB)=0
2117   190   CONTINUE
2118         CKIN(1)=2D0
2119         CKIN(3)=0D0
2120         MSTP(2)=1
2121         MSTP(11)=0
2122         MSTP(33)=0
2123         MSTP(81)=1
2124         MSTP(82)=1
2125         MSTP(111)=1
2126         MSTP(131)=0
2127         MSTP(133)=0
2128         PARP(131)=0.01D0
2129  
2130 C...Prompt photon production at fixed target.
2131         IF(IPROC.EQ.1) THEN
2132           PZSUM=300D0
2133           PESUM=SQRT(PZSUM**2+PYMASS(211)**2)+PYMASS(2212)
2134           PQSUM=2D0
2135           MSEL=10
2136           CKIN(3)=5D0
2137           CALL PYINIT('FIXT','pi+','p',PZSUM)
2138  
2139 C...QCD processes at ISR energies.
2140         ELSEIF(IPROC.EQ.2) THEN
2141           PESUM=63D0
2142           PZSUM=0D0
2143           PQSUM=2D0
2144           MSEL=1
2145           CKIN(3)=5D0
2146           CALL PYINIT('CMS','p','p',PESUM)
2147  
2148 C...W production + multiple interactions at CERN Collider.
2149         ELSEIF(IPROC.EQ.3) THEN
2150           PESUM=630D0
2151           PZSUM=0D0
2152           PQSUM=0D0
2153           MSEL=12
2154           CKIN(1)=20D0
2155           MSTP(82)=4
2156           MSTP(2)=2
2157           MSTP(33)=3
2158           CALL PYINIT('CMS','p','pbar',PESUM)
2159  
2160 C...W/Z gauge boson pairs + pileup events at the Tevatron.
2161         ELSEIF(IPROC.EQ.4) THEN
2162           PESUM=1800D0
2163           PZSUM=0D0
2164           PQSUM=0D0
2165           MSUB(22)=1
2166           MSUB(23)=1
2167           MSUB(25)=1
2168           CKIN(1)=200D0
2169           MSTP(111)=0
2170           MSTP(131)=1
2171           MSTP(133)=2
2172           PARP(131)=0.04D0
2173           CALL PYINIT('CMS','p','pbar',PESUM)
2174  
2175 C...Higgs production at LHC.
2176         ELSEIF(IPROC.EQ.5) THEN
2177           PESUM=15400D0
2178           PZSUM=0D0
2179           PQSUM=2D0
2180           MSUB(3)=1
2181           MSUB(102)=1
2182           MSUB(123)=1
2183           MSUB(124)=1
2184           PMAS(25,1)=300D0
2185           CKIN(1)=200D0
2186           MSTP(81)=0
2187           MSTP(111)=0
2188           CALL PYINIT('CMS','p','p',PESUM)
2189  
2190 C...Z' production at SSC.
2191         ELSEIF(IPROC.EQ.6) THEN
2192           PESUM=40000D0
2193           PZSUM=0D0
2194           PQSUM=2D0
2195           MSEL=21
2196           PMAS(32,1)=600D0
2197           CKIN(1)=400D0
2198           MSTP(81)=0
2199           MSTP(111)=0
2200           CALL PYINIT('CMS','p','p',PESUM)
2201  
2202 C...W pair production at 1 TeV e+e- collider.
2203         ELSEIF(IPROC.EQ.7) THEN
2204           PESUM=1000D0
2205           PZSUM=0D0
2206           PQSUM=0D0
2207           MSUB(25)=1
2208           MSUB(69)=1
2209           MSTP(11)=1
2210           CALL PYINIT('CMS','e+','e-',PESUM)
2211  
2212 C...Deep inelastic scattering at a LEP+LHC ep collider.
2213         ELSEIF(IPROC.EQ.8) THEN
2214           P(1,1)=0D0
2215           P(1,2)=0D0
2216           P(1,3)=8000D0
2217           P(2,1)=0D0
2218           P(2,2)=0D0
2219           P(2,3)=-80D0
2220           PESUM=8080D0
2221           PZSUM=7920D0
2222           PQSUM=0D0
2223           MSUB(10)=1
2224           CKIN(3)=50D0
2225           MSTP(111)=0
2226           CALL PYINIT('3MOM','p','e-',PESUM)
2227         ENDIF
2228  
2229 C...Generate 20 events of each required type.
2230         DO 220 IEV=1,20
2231           CALL PYEVNT
2232           PESUMM=PESUM
2233           IF(IPROC.EQ.4) PESUMM=MSTI(41)*PESUM
2234  
2235 C...Check conservation of energy/momentum/flavour.
2236           PINI(1)=0D0
2237           PINI(2)=0D0
2238           PINI(3)=PZSUM
2239           PINI(4)=PESUMM
2240           PINI(6)=PQSUM
2241           DO 200 J=1,4
2242             PFIN(J)=PYP(0,J)
2243   200     CONTINUE
2244           PFIN(6)=PYP(0,6)
2245           MERR=0
2246           DEVE=ABS(PFIN(4)-PINI(4))+ABS(PFIN(3)-PINI(3))
2247           DEVT=ABS(PFIN(1)-PINI(1))+ABS(PFIN(2)-PINI(2))
2248           DEVQ=ABS(PFIN(6)-PINI(6))
2249           IF(DEVE.GT.2D-3*PESUM.OR.DEVT.GT.MAX(0.01D0,1D-4*PESUM).OR.
2250      &    DEVQ.GT.0.1D0) MERR=1
2251           IF(MERR.NE.0) WRITE(MSTU(11),5000) (PINI(J),J=1,4),PINI(6),
2252      &    (PFIN(J),J=1,4),PFIN(6)
2253  
2254 C...Check that all KF codes are known ones, and that partons/particles
2255 C...satisfy energy-momentum-mass relation.
2256           DO 210 I=1,N
2257             IF(K(I,1).GT.20) GOTO 210
2258             IF(PYCOMP(K(I,2)).EQ.0) THEN
2259               WRITE(MSTU(11),5100) I
2260               MERR=MERR+1
2261             ENDIF
2262             PD=P(I,4)**2-P(I,1)**2-P(I,2)**2-P(I,3)**2-P(I,5)**2*
2263      &      SIGN(1D0,P(I,5))
2264             IF(ABS(PD).GT.MAX(0.1D0,0.002D0*P(I,4)**2,0.002D0*P(I,5)**2)
2265      &      .OR.(P(I,5).GE.0D0.AND.P(I,4).LT.0D0)) THEN
2266               WRITE(MSTU(11),5200) I
2267               MERR=MERR+1
2268             ENDIF
2269   210     CONTINUE
2270  
2271 C...Listing of erroneous events, and first event of each type.
2272           IF(MERR.GE.1) NERR=NERR+1
2273           IF(NERR.GE.10) THEN
2274             WRITE(MSTU(11),6300)
2275             CALL PYLIST(1)
2276             STOP
2277           ENDIF
2278           IF(MTEST.GE.1.AND.(MERR.GE.1.OR.IEV.EQ.1)) THEN
2279             IF(MERR.GE.1) WRITE(MSTU(11),6400)
2280             CALL PYLIST(1)
2281           ENDIF
2282   220   CONTINUE
2283  
2284 C...List statistics for each process type.
2285         IF(MTEST.GE.1) CALL PYSTAT(1)
2286   230 CONTINUE
2287  
2288 C...Summarize result of run.
2289       IF(NERR.EQ.0) WRITE(MSTU(11),6500)
2290       IF(NERR.GT.0) WRITE(MSTU(11),6600) NERR
2291  
2292 C...Format statements for output.
2293  5000 FORMAT(/' Momentum, energy and/or charge were not conserved ',
2294      &'in following event'/' sum of',9X,'px',11X,'py',11X,'pz',11X,
2295      &'E',8X,'charge'/' before',2X,4(1X,F12.5),1X,F8.2/' after',3X,
2296      &4(1X,F12.5),1X,F8.2)
2297  5100 FORMAT(/5X,'Entry no.',I4,' in following event not known code')
2298  5200 FORMAT(/5X,'Entry no.',I4,' in following event has faulty ',
2299      &'kinematics')
2300  6300 FORMAT(/5X,'This is the tenth error experienced! Something is ',
2301      &'wrong.'/5X,'Execution will be stopped after listing of event.')
2302  6400 FORMAT(5X,'Faulty event follows:')
2303  6500 FORMAT(//5X,'End result of PYTEST: no errors detected.')
2304  6600 FORMAT(//5X,'End result of PYTEST:',I2,' errors detected.'/
2305      &5X,'This should not have happened!')
2306  
2307       RETURN
2308       END
2309  
2310 C*********************************************************************
2311  
2312 C...PYHEPC
2313 C...Converts PYTHIA event record contents to or from
2314 C...the standard event record commonblock.
2315  
2316       SUBROUTINE PYHEPC(MCONV)
2317  
2318 C...Double precision and integer declarations.
2319       IMPLICIT DOUBLE PRECISION(A-H, O-Z)
2320       IMPLICIT INTEGER(I-N)
2321       INTEGER PYK,PYCHGE,PYCOMP
2322 C...Commonblocks.
2323       COMMON/PYJETS/N,NPAD,K(4000,5),P(4000,5),V(4000,5)
2324       COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
2325       COMMON/PYDAT2/KCHG(500,4),PMAS(500,4),PARF(2000),VCKM(4,4)
2326       SAVE /PYJETS/,/PYDAT1/,/PYDAT2/
2327 C...HEPEVT commonblock.
2328       PARAMETER (NMXHEP=4000)
2329       COMMON/HEPEVT/NEVHEP,NHEP,ISTHEP(NMXHEP),IDHEP(NMXHEP),
2330      &JMOHEP(2,NMXHEP),JDAHEP(2,NMXHEP),PHEP(5,NMXHEP),VHEP(4,NMXHEP)
2331       DOUBLE PRECISION PHEP,VHEP
2332       SAVE /HEPEVT/
2333  
2334 C...Conversion from PYTHIA to standard, the easy part.
2335       IF(MCONV.EQ.1) THEN
2336         NEVHEP=0
2337         IF(N.GT.NMXHEP) CALL PYERRM(8,
2338      &  '(PYHEPC:) no more space in /HEPEVT/')
2339         NHEP=MIN(N,NMXHEP)
2340         DO 140 I=1,NHEP
2341           ISTHEP(I)=0
2342           IF(K(I,1).GE.1.AND.K(I,1).LE.10) ISTHEP(I)=1
2343           IF(K(I,1).GE.11.AND.K(I,1).LE.20) ISTHEP(I)=2
2344           IF(K(I,1).GE.21.AND.K(I,1).LE.30) ISTHEP(I)=3
2345           IF(K(I,1).GE.31.AND.K(I,1).LE.100) ISTHEP(I)=K(I,1)
2346           IDHEP(I)=K(I,2)
2347           JMOHEP(1,I)=K(I,3)
2348           JMOHEP(2,I)=0
2349           IF(K(I,1).NE.3.AND.K(I,1).NE.13.AND.K(I,1).NE.14) THEN
2350             JDAHEP(1,I)=K(I,4)
2351             JDAHEP(2,I)=K(I,5)
2352           ELSE
2353             JDAHEP(1,I)=0
2354             JDAHEP(2,I)=0
2355           ENDIF
2356           DO 100 J=1,5
2357             PHEP(J,I)=P(I,J)
2358   100     CONTINUE
2359           DO 110 J=1,4
2360             VHEP(J,I)=V(I,J)
2361   110     CONTINUE
2362  
2363 C...Check if new event (from pileup).
2364           IF(I.EQ.1) THEN
2365             INEW=1
2366           ELSE
2367             IF(K(I,1).EQ.21.AND.K(I-1,1).NE.21) INEW=I
2368           ENDIF
2369  
2370 C...Fill in missing mother information.
2371           IF(I.GE.INEW+2.AND.K(I,1).EQ.21.AND.K(I,3).EQ.0) THEN
2372             IMO1=I-2
2373             IF(I.GE.INEW+3.AND.K(I-1,1).EQ.21.AND.K(I-1,3).EQ.0)
2374      &      IMO1=IMO1-1
2375             JMOHEP(1,I)=IMO1
2376             JMOHEP(2,I)=IMO1+1
2377           ELSEIF(K(I,2).GE.91.AND.K(I,2).LE.93) THEN
2378             I1=K(I,3)-1
2379   120       I1=I1+1
2380             IF(I1.GE.I) CALL PYERRM(8,
2381      &      '(PYHEPC:) translation of inconsistent event history')
2382             IF(I1.LT.I.AND.K(I1,1).NE.1.AND.K(I1,1).NE.11) GOTO 120
2383             KC=PYCOMP(K(I1,2))
2384             IF(I1.LT.I.AND.KC.EQ.0) GOTO 120
2385             IF(I1.LT.I.AND.KCHG(KC,2).EQ.0) GOTO 120
2386             JMOHEP(2,I)=I1
2387           ELSEIF(K(I,2).EQ.94) THEN
2388             NJET=2
2389             IF(NHEP.GE.I+3.AND.K(I+3,3).LE.I) NJET=3
2390             IF(NHEP.GE.I+4.AND.K(I+4,3).LE.I) NJET=4
2391             JMOHEP(2,I)=MOD(K(I+NJET,4)/MSTU(5),MSTU(5))
2392             IF(JMOHEP(2,I).EQ.JMOHEP(1,I)) JMOHEP(2,I)=
2393      &      MOD(K(I+1,4)/MSTU(5),MSTU(5))
2394           ENDIF
2395  
2396 C...Fill in missing daughter information.
2397           IF(K(I,2).EQ.94.AND.MSTU(16).NE.2) THEN
2398             DO 130 I1=JDAHEP(1,I),JDAHEP(2,I)
2399               I2=MOD(K(I1,4)/MSTU(5),MSTU(5))
2400               JDAHEP(1,I2)=I
2401   130       CONTINUE
2402           ENDIF
2403           IF(K(I,2).GE.91.AND.K(I,2).LE.94) GOTO 140
2404           I1=JMOHEP(1,I)
2405           IF(I1.LE.0.OR.I1.GT.NHEP) GOTO 140
2406           IF(K(I1,1).NE.13.AND.K(I1,1).NE.14) GOTO 140
2407           IF(JDAHEP(1,I1).EQ.0) THEN
2408             JDAHEP(1,I1)=I
2409           ELSE
2410             JDAHEP(2,I1)=I
2411           ENDIF
2412   140   CONTINUE
2413         DO 150 I=1,NHEP
2414           IF(K(I,1).NE.13.AND.K(I,1).NE.14) GOTO 150
2415           IF(JDAHEP(2,I).EQ.0) JDAHEP(2,I)=JDAHEP(1,I)
2416   150   CONTINUE
2417  
2418 C...Conversion from standard to PYTHIA, the easy part.
2419       ELSE
2420         IF(NHEP.GT.MSTU(4)) CALL PYERRM(8,
2421      &  '(PYHEPC:) no more space in /PYJETS/')
2422         N=MIN(NHEP,MSTU(4))
2423         NKQ=0
2424         KQSUM=0
2425         DO 180 I=1,N
2426           K(I,1)=0
2427           IF(ISTHEP(I).EQ.1) K(I,1)=1
2428           IF(ISTHEP(I).EQ.2) K(I,1)=11
2429           IF(ISTHEP(I).EQ.3) K(I,1)=21
2430           K(I,2)=IDHEP(I)
2431           K(I,3)=JMOHEP(1,I)
2432           K(I,4)=JDAHEP(1,I)
2433           K(I,5)=JDAHEP(2,I)
2434           DO 160 J=1,5
2435             P(I,J)=PHEP(J,I)
2436   160     CONTINUE
2437           DO 170 J=1,4
2438             V(I,J)=VHEP(J,I)
2439   170     CONTINUE
2440           V(I,5)=0D0
2441           IF(ISTHEP(I).EQ.2.AND.PHEP(4,I).GT.PHEP(5,I)) THEN
2442             I1=JDAHEP(1,I)
2443             IF(I1.GT.0.AND.I1.LE.NHEP) V(I,5)=(VHEP(4,I1)-VHEP(4,I))*
2444      &      PHEP(5,I)/PHEP(4,I)
2445           ENDIF
2446  
2447 C...Fill in missing information on colour connection in jet systems.
2448           IF(ISTHEP(I).EQ.1) THEN
2449             KC=PYCOMP(K(I,2))
2450             KQ=0
2451             IF(KC.NE.0) KQ=KCHG(KC,2)*ISIGN(1,K(I,2))
2452             IF(KQ.NE.0) NKQ=NKQ+1
2453             IF(KQ.NE.2) KQSUM=KQSUM+KQ
2454             IF(KQ.NE.0.AND.KQSUM.NE.0) THEN
2455               K(I,1)=2
2456             ELSEIF(KQ.EQ.2.AND.I.LT.N) THEN
2457               IF(K(I+1,2).EQ.21) K(I,1)=2
2458             ENDIF
2459           ENDIF
2460   180   CONTINUE
2461         IF(NKQ.EQ.1.OR.KQSUM.NE.0) CALL PYERRM(8,
2462      &  '(PYHEPC:) input parton configuration not colour singlet')
2463       ENDIF
2464  
2465       END
2466  
2467 C*********************************************************************
2468  
2469 C...PYINIT
2470 C...Initializes the generation procedure; finds maxima of the
2471 C...differential cross-sections to be used for weighting.
2472  
2473       SUBROUTINE PYINIT(FRAME,BEAM,TARGET,WIN)
2474  
2475 C...Double precision and integer declarations.
2476       IMPLICIT DOUBLE PRECISION(A-H, O-Z)
2477       IMPLICIT INTEGER(I-N)
2478       INTEGER PYK,PYCHGE,PYCOMP
2479 C...Commonblocks.
2480       COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
2481       COMMON/PYDAT2/KCHG(500,4),PMAS(500,4),PARF(2000),VCKM(4,4)
2482       COMMON/PYDAT3/MDCY(500,3),MDME(8000,2),BRAT(8000),KFDP(8000,5)
2483       COMMON/PYDAT4/CHAF(500,2)
2484       CHARACTER CHAF*16
2485       COMMON/PYSUBS/MSEL,MSELPD,MSUB(500),KFIN(2,-40:40),CKIN(200)
2486       COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
2487       COMMON/PYINT1/MINT(400),VINT(400)
2488       COMMON/PYINT2/ISET(500),KFPR(500,2),COEF(500,20),ICOL(40,4,2)
2489       COMMON/PYINT5/NGENPD,NGEN(0:500,3),XSEC(0:500,3)
2490       SAVE /PYDAT1/,/PYDAT2/,/PYDAT3/,/PYDAT4/,/PYSUBS/,/PYPARS/,
2491      &/PYINT1/,/PYINT2/,/PYINT5/
2492 C...Local arrays and character variables.
2493       DIMENSION ALAMIN(20),NFIN(20)
2494       CHARACTER*(*) FRAME,BEAM,TARGET
2495       CHARACTER CHFRAM*12,CHBEAM*12,CHTARG*12,CHLH(2)*6
2496  
2497 C...Interface to PDFLIB.
2498       COMMON/W50512/QCDL4,QCDL5
2499       SAVE /W50512/
2500       DOUBLE PRECISION VALUE(20),QCDL4,QCDL5
2501       CHARACTER*20 PARM(20)
2502       DATA VALUE/20*0D0/,PARM/20*' '/
2503  
2504 C...Data:Lambda and n_f values for parton distributions..
2505       DATA ALAMIN/0.177D0,0.239D0,0.247D0,0.2322D0,0.248D0,0.248D0,
2506      &0.192D0,0.326D0,2*0.2D0,0.2D0,0.2D0,0.29D0,0.2D0,0.4D0,5*0.2D0/,
2507      &NFIN/20*4/
2508       DATA CHLH/'lepton','hadron'/
2509  
2510 C...Reset MINT and VINT arrays. Write headers.
2511       MSTI(53)=0
2512       DO 100 J=1,400
2513         MINT(J)=0
2514         VINT(J)=0D0
2515   100 CONTINUE
2516       IF(MSTU(12).GE.1) CALL PYLIST(0)
2517       IF(MSTP(122).GE.1) WRITE(MSTU(11),5100)
2518  
2519 C...Call user process initialization routine.
2520       IF(FRAME(1:1).EQ.'u'.OR.FRAME(1:1).EQ.'U') THEN
2521         MSEL=0
2522         CALL UPINIT
2523         MSEL=0
2524       ENDIF
2525  
2526 C...Maximum 4 generations; set maximum number of allowed flavours.
2527       MSTP(1)=MIN(4,MSTP(1))
2528       MSTU(114)=MIN(MSTU(114),2*MSTP(1))
2529       MSTP(58)=MIN(MSTP(58),2*MSTP(1))
2530  
2531 C...Sum up Cabibbo-Kobayashi-Maskawa factors for each quark/lepton.
2532       DO 120 I=-20,20
2533         VINT(180+I)=0D0
2534         IA=IABS(I)
2535         IF(IA.GE.1.AND.IA.LE.2*MSTP(1)) THEN
2536           DO 110 J=1,MSTP(1)
2537             IB=2*J-1+MOD(IA,2)
2538             IF(IB.GE.6.AND.MSTP(9).EQ.0) GOTO 110
2539             IPM=(5-ISIGN(1,I))/2
2540             IDC=J+MDCY(IA,2)+2
2541             IF(MDME(IDC,1).EQ.1.OR.MDME(IDC,1).EQ.IPM) VINT(180+I)=
2542      &      VINT(180+I)+VCKM((IA+1)/2,(IB+1)/2)
2543   110     CONTINUE
2544         ELSEIF(IA.GE.11.AND.IA.LE.10+2*MSTP(1)) THEN
2545           VINT(180+I)=1D0
2546         ENDIF
2547   120 CONTINUE
2548  
2549 C...Initialize parton distributions: PDFLIB.
2550       IF(MSTP(52).EQ.2) THEN
2551         PARM(1)='NPTYPE'
2552         VALUE(1)=1
2553         PARM(2)='NGROUP'
2554         VALUE(2)=MSTP(51)/1000
2555         PARM(3)='NSET'
2556         VALUE(3)=MOD(MSTP(51),1000)
2557         PARM(4)='TMAS'
2558         VALUE(4)=PMAS(6,1)
2559 C.... ALICE
2560         CALL PDFSET_ALICE(PARM,VALUE)
2561         MINT(93)=1000000+MSTP(51)
2562       ENDIF
2563  
2564 C...Choose Lambda value to use in alpha-strong.
2565       MSTU(111)=MSTP(2)
2566       IF(MSTP(3).GE.2) THEN
2567         ALAM=0.2D0
2568         NF=4
2569         IF(MSTP(52).EQ.1.AND.MSTP(51).GE.1.AND.MSTP(51).LE.20) THEN
2570           ALAM=ALAMIN(MSTP(51))
2571           NF=NFIN(MSTP(51))
2572         ELSEIF(MSTP(52).EQ.2) THEN
2573           ALAM=QCDL4
2574           NF=4
2575         ENDIF
2576         PARP(1)=ALAM
2577         PARP(61)=ALAM
2578         PARP(72)=ALAM
2579         PARU(112)=ALAM
2580         MSTU(112)=NF
2581         IF(MSTP(3).EQ.3) PARJ(81)=ALAM
2582       ENDIF
2583  
2584 C...Initialize the SUSY generation: couplings, masses,
2585 C...decay modes, branching ratios, and so on.
2586       CALL PYMSIN
2587 C...Initialize widths and partial widths for resonances.
2588       CALL PYINRE
2589 C...Set Z0 mass and width for e+e- routines.
2590       PARJ(123)=PMAS(23,1)
2591       PARJ(124)=PMAS(23,2)
2592  
2593 C...Identify beam and target particles and frame of process.
2594       CHFRAM=FRAME//' '
2595       CHBEAM=BEAM//' '
2596       CHTARG=TARGET//' '
2597       CALL PYINBM(CHFRAM,CHBEAM,CHTARG,WIN)
2598       IF(MINT(65).EQ.1) GOTO 170
2599  
2600 C...For gamma-p or gamma-gamma allow many (3 or 6) alternatives.
2601 C...For e-gamma allow 2 alternatives.
2602       MINT(121)=1
2603       IF(MSTP(14).EQ.10.AND.(MSEL.EQ.1.OR.MSEL.EQ.2)) THEN
2604         IF((MINT(11).EQ.22.OR.MINT(12).EQ.22).AND.
2605      &  (IABS(MINT(11)).GT.100.OR.IABS(MINT(12)).GT.100)) MINT(121)=3
2606         IF(MINT(11).EQ.22.AND.MINT(12).EQ.22) MINT(121)=6
2607         IF((MINT(11).EQ.22.OR.MINT(12).EQ.22).AND.
2608      &  (IABS(MINT(11)).EQ.11.OR.IABS(MINT(12)).EQ.11)) MINT(121)=2
2609       ELSEIF(MSTP(14).EQ.20.AND.(MSEL.EQ.1.OR.MSEL.EQ.2)) THEN
2610         IF((MINT(11).EQ.22.OR.MINT(12).EQ.22).AND.
2611      &  (IABS(MINT(11)).GT.100.OR.IABS(MINT(12)).GT.100)) MINT(121)=3
2612         IF(MINT(11).EQ.22.AND.MINT(12).EQ.22) MINT(121)=9
2613       ELSEIF(MSTP(14).EQ.25.AND.(MSEL.EQ.1.OR.MSEL.EQ.2)) THEN
2614         IF((MINT(11).EQ.22.OR.MINT(12).EQ.22).AND.
2615      &  (IABS(MINT(11)).GT.100.OR.IABS(MINT(12)).GT.100)) MINT(121)=2
2616         IF(MINT(11).EQ.22.AND.MINT(12).EQ.22) MINT(121)=4
2617       ELSEIF(MSTP(14).EQ.30.AND.(MSEL.EQ.1.OR.MSEL.EQ.2)) THEN
2618         IF((MINT(11).EQ.22.OR.MINT(12).EQ.22).AND.
2619      &  (IABS(MINT(11)).GT.100.OR.IABS(MINT(12)).GT.100)) MINT(121)=4
2620         IF(MINT(11).EQ.22.AND.MINT(12).EQ.22) MINT(121)=13
2621       ENDIF
2622       MINT(123)=MSTP(14)
2623       IF((MSTP(14).EQ.10.OR.MSTP(14).EQ.20.OR.MSTP(14).EQ.25.OR.
2624      &MSTP(14).EQ.30).AND.MSEL.NE.1.AND.MSEL.NE.2) MINT(123)=0
2625       IF(MSTP(14).GE.11.AND.MSTP(14).LE.19) THEN
2626         IF(MSTP(14).EQ.11) MINT(123)=0
2627         IF(MSTP(14).EQ.12.OR.MSTP(14).EQ.14) MINT(123)=5
2628         IF(MSTP(14).EQ.13.OR.MSTP(14).EQ.17) MINT(123)=6
2629         IF(MSTP(14).EQ.15) MINT(123)=2
2630         IF(MSTP(14).EQ.16.OR.MSTP(14).EQ.18) MINT(123)=7
2631         IF(MSTP(14).EQ.19) MINT(123)=3
2632       ELSEIF(MSTP(14).GE.21.AND.MSTP(14).LE.24) THEN
2633         IF(MSTP(14).EQ.21) MINT(123)=0
2634         IF(MSTP(14).EQ.22.OR.MSTP(14).EQ.23) MINT(123)=4
2635         IF(MSTP(14).EQ.24) MINT(123)=1
2636       ELSEIF(MSTP(14).GE.26.AND.MSTP(14).LE.29) THEN
2637         IF(MSTP(14).EQ.26.OR.MSTP(14).EQ.28) MINT(123)=8
2638         IF(MSTP(14).EQ.27.OR.MSTP(14).EQ.29) MINT(123)=9
2639       ENDIF
2640  
2641 C...Set up kinematics of process.
2642       CALL PYINKI(0)
2643  
2644 C...Set up kinematics for photons inside leptons.
2645       IF(MINT(141).NE.0.OR.MINT(142).NE.0) CALL PYGAGA(1,WTGAGA)
2646  
2647 C...Precalculate flavour selection weights.
2648       CALL PYKFIN
2649  
2650 C...Loop over gamma-p or gamma-gamma alternatives.
2651       CKIN3=CKIN(3)
2652       MSAV48=0
2653       DO 160 IGA=1,MINT(121)
2654         CKIN(3)=CKIN3
2655         MINT(122)=IGA
2656  
2657 C...Select partonic subprocesses to be included in the simulation.
2658         CALL PYINPR
2659         MINT(101)=1
2660         MINT(102)=1
2661         MINT(103)=MINT(11)
2662         MINT(104)=MINT(12)
2663  
2664 C...Count number of subprocesses on.
2665         MINT(48)=0
2666         DO 130 ISUB=1,500
2667           IF(MINT(50).EQ.0.AND.ISUB.GE.91.AND.ISUB.LE.96.AND.
2668      &    MSUB(ISUB).EQ.1.AND.MINT(121).GT.1) THEN
2669             MSUB(ISUB)=0
2670           ELSEIF(MINT(50).EQ.0.AND.ISUB.GE.91.AND.ISUB.LE.96.AND.
2671      &    MSUB(ISUB).EQ.1) THEN
2672             WRITE(MSTU(11),5200) ISUB,CHLH(MINT(41)),CHLH(MINT(42))
2673             STOP
2674           ELSEIF(MSUB(ISUB).EQ.1.AND.ISET(ISUB).EQ.-1) THEN
2675             WRITE(MSTU(11),5300) ISUB
2676             STOP
2677           ELSEIF(MSUB(ISUB).EQ.1.AND.ISET(ISUB).LE.-2) THEN
2678             WRITE(MSTU(11),5400) ISUB
2679             STOP
2680           ELSEIF(MSUB(ISUB).EQ.1) THEN
2681             MINT(48)=MINT(48)+1
2682           ENDIF
2683   130   CONTINUE
2684
2685 C...Stop or raise warning flag if no subprocesses on.
2686         IF(MINT(121).EQ.1.AND.MINT(48).EQ.0) THEN
2687           IF(MSTP(127).NE.1) THEN
2688             WRITE(MSTU(11),5500)
2689             STOP
2690           ELSE
2691             WRITE(MSTU(11),5700)
2692             MSTI(53)=1
2693           ENDIF  
2694         ENDIF
2695         MINT(49)=MINT(48)-MSUB(91)-MSUB(92)-MSUB(93)-MSUB(94)
2696         MSAV48=MSAV48+MINT(48)
2697  
2698 C...Reset variables for cross-section calculation.
2699         DO 150 I=0,500
2700           DO 140 J=1,3
2701             NGEN(I,J)=0
2702             XSEC(I,J)=0D0
2703   140     CONTINUE
2704   150   CONTINUE
2705  
2706 C...Find parametrized total cross-sections.
2707         CALL PYXTOT
2708         VINT(318)=VINT(317)
2709  
2710 C...Maxima of differential cross-sections.
2711         IF(MSTP(121).LE.1) CALL PYMAXI
2712  
2713 C...Initialize possibility of pileup events.
2714         IF(MINT(121).GT.1) MSTP(131)=0
2715         IF(MSTP(131).NE.0) CALL PYPILE(1)
2716  
2717 C...Initialize multiple interactions with variable impact parameter.
2718         IF(MINT(50).EQ.1.AND.(MINT(49).NE.0.OR.MSTP(131).NE.0).AND.
2719      &  MSTP(82).GE.2) CALL PYMULT(1)
2720  
2721 C...Save results for gamma-p and gamma-gamma alternatives.
2722         IF(MINT(121).GT.1) CALL PYSAVE(1,IGA)
2723   160 CONTINUE
2724  
2725 C...Initialization finished.
2726       IF(MSAV48.EQ.0) THEN
2727         IF(MSTP(127).NE.1) THEN
2728           WRITE(MSTU(11),5500)
2729           STOP
2730         ELSE
2731           WRITE(MSTU(11),5700)
2732           MSTI(53)=1
2733         ENDIF  
2734       ENDIF
2735   170 IF(MSTP(122).GE.1) WRITE(MSTU(11),5600)
2736  
2737 C...Formats for initialization information.
2738  5100 FORMAT('1',18('*'),1X,'PYINIT: initialization of PYTHIA ',
2739      &'routines',1X,17('*'))
2740  5200 FORMAT(1X,'Error: process number ',I3,' not meaningful for ',A6,
2741      &'-',A6,' interactions.'/1X,'Execution stopped!')
2742  5300 FORMAT(1X,'Error: requested subprocess',I4,' not implemented.'/
2743      &1X,'Execution stopped!')
2744  5400 FORMAT(1X,'Error: requested subprocess',I4,' not existing.'/
2745      &1X,'Execution stopped!')
2746  5500 FORMAT(1X,'Error: no subprocess switched on.'/
2747      &1X,'Execution stopped.')
2748  5600 FORMAT(/1X,22('*'),1X,'PYINIT: initialization completed',1X,
2749      &22('*'))
2750  5700 FORMAT(1X,'Error: no subprocess switched on.'/
2751      &1X,'Execution will stop if you try to generate events.')
2752  
2753       RETURN
2754       END
2755  
2756 C*********************************************************************
2757  
2758 C...PYEVNT
2759 C...Administers the generation of a high-pT event via calls to
2760 C...a number of subroutines.
2761  
2762       SUBROUTINE PYEVNT
2763  
2764 C...Double precision and integer declarations.
2765       IMPLICIT DOUBLE PRECISION(A-H, O-Z)
2766       IMPLICIT INTEGER(I-N)
2767       INTEGER PYK,PYCHGE,PYCOMP
2768 C...Commonblocks.
2769       COMMON/PYJETS/N,NPAD,K(4000,5),P(4000,5),V(4000,5)
2770       COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
2771       COMMON/PYDAT2/KCHG(500,4),PMAS(500,4),PARF(2000),VCKM(4,4)
2772       COMMON/PYDAT3/MDCY(500,3),MDME(8000,2),BRAT(8000),KFDP(8000,5)
2773       COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
2774       COMMON/PYINT1/MINT(400),VINT(400)
2775       COMMON/PYINT2/ISET(500),KFPR(500,2),COEF(500,20),ICOL(40,4,2)
2776       COMMON/PYINT4/MWID(500),WIDS(500,5)
2777       COMMON/PYINT5/NGENPD,NGEN(0:500,3),XSEC(0:500,3)
2778       SAVE /PYJETS/,/PYDAT1/,/PYDAT2/,/PYDAT3/,/PYPARS/,/PYINT1/,
2779      &/PYINT2/,/PYINT4/,/PYINT5/
2780 C...Local array.
2781       DIMENSION VTX(4)
2782
2783 C...Stop if no subprocesses on.
2784       IF(MINT(121).EQ.1.AND.MSTI(53).EQ.1) THEN
2785         WRITE(MSTU(11),5100)
2786         STOP
2787       ENDIF  
2788  
2789 C...Initial values for some counters.
2790       N=0
2791       MINT(5)=MINT(5)+1
2792       MINT(7)=0
2793       MINT(8)=0
2794       MINT(83)=0
2795       MINT(84)=MSTP(126)
2796       MSTU(24)=0
2797       MSTU70=0
2798       MSTJ14=MSTJ(14)
2799  
2800 C...If variable energies: redo incoming kinematics and cross-section.
2801       MSTI(61)=0
2802       IF(MSTP(171).EQ.1) THEN
2803         CALL PYINKI(1)
2804         IF(MSTI(61).EQ.1) THEN
2805           MINT(5)=MINT(5)-1
2806           RETURN
2807         ENDIF
2808         IF(MINT(121).GT.1) CALL PYSAVE(3,1)
2809         CALL PYXTOT
2810       ENDIF
2811  
2812 C...Loop over number of pileup events; check space left.
2813       IF(MSTP(131).LE.0) THEN
2814         NPILE=1
2815       ELSE
2816         CALL PYPILE(2)
2817         NPILE=MINT(81)
2818       ENDIF
2819       DO 250 IPILE=1,NPILE
2820         IF(MINT(84)+100.GE.MSTU(4)) THEN
2821           CALL PYERRM(11,
2822      &    '(PYEVNT:) no more space in PYJETS for pileup events')
2823           IF(MSTU(21).GE.1) GOTO 260
2824         ENDIF
2825         MINT(82)=IPILE
2826  
2827 C...Generate variables of hard scattering.
2828         MINT(51)=0
2829         MSTI(52)=0
2830   100   CONTINUE
2831         IF(MINT(51).NE.0.OR.MSTU(24).NE.0) MSTI(52)=MSTI(52)+1
2832         MINT(31)=0
2833         MINT(51)=0
2834         MINT(57)=0
2835         CALL PYRAND
2836         IF(MSTI(61).EQ.1) THEN
2837           MINT(5)=MINT(5)-1
2838           RETURN
2839         ENDIF
2840         IF(MINT(51).EQ.2) RETURN
2841         ISUB=MINT(1)
2842         IF(MSTP(111).EQ.-1) GOTO 240
2843  
2844         IF((ISUB.LE.90.OR.ISUB.GE.95).AND.ISUB.NE.99) THEN
2845 C...Hard scattering (including low-pT):
2846 C...reconstruct kinematics and colour flow of hard scattering.
2847           MINT31=MINT(31)
2848   110     MINT(31)=MINT31
2849           MINT(51)=0
2850           CALL PYSCAT
2851           IF(MINT(51).EQ.1) GOTO 100
2852           IPU1=MINT(84)+1
2853           IPU2=MINT(84)+2
2854           IF(ISUB.EQ.95) GOTO 120
2855  
2856 C...Showering of initial state partons (optional).
2857           NFIN=N
2858           ALAMSV=PARJ(81)
2859           PARJ(81)=PARP(72)
2860           IF(MSTP(61).GE.1.AND.MINT(47).GE.2) CALL PYSSPA(IPU1,IPU2)
2861           PARJ(81)=ALAMSV
2862           IF(MINT(51).EQ.1) GOTO 100
2863  
2864 C...Showering of final state partons (optional).
2865           ALAMSV=PARJ(81)
2866           PARJ(81)=PARP(72)
2867           IF(MSTP(71).GE.1.AND.ISET(ISUB).GE.2.AND.ISET(ISUB).LE.10)
2868      &    THEN
2869             IPU3=MINT(84)+3
2870             IPU4=MINT(84)+4
2871             IF(ISET(ISUB).EQ.5) IPU4=-3
2872             QMAX=VINT(55)
2873             IF(ISET(ISUB).EQ.2) QMAX=SQRT(PARP(71))*VINT(55)
2874             CALL PYSHOW(IPU3,IPU4,QMAX)
2875           ELSEIF(ISET(ISUB).EQ.11) THEN
2876             CALL PYADSH(NFIN)
2877           ENDIF
2878           PARJ(81)=ALAMSV
2879  
2880 C...Decay of final state resonances.
2881           MINT(32)=0
2882           IF(MSTP(41).GE.1.AND.ISET(ISUB).LE.10) CALL PYRESD(0)
2883           IF(MINT(51).EQ.1) GOTO 100
2884           MINT(52)=N
2885  
2886 C...Multiple interactions.
2887           IF(MSTP(81).GE.1.AND.MINT(50).EQ.1) CALL PYMULT(6)
2888           MINT(53)=N
2889  
2890 C...Hadron remnants and primordial kT.
2891   120     CALL PYREMN(IPU1,IPU2)
2892           IF(MINT(51).EQ.1.AND.MINT(57).GE.1.AND.MINT(57).LE.5) GOTO 110
2893           IF(MINT(51).EQ.1) GOTO 100
2894  
2895          ELSEIF(ISUB.NE.99) THEN
2896 C...Diffractive and elastic scattering.
2897           CALL PYDIFF
2898  
2899         ELSE
2900 C...DIS scattering (photon flux external).
2901           CALL PYDISG
2902           IF(MINT(51).EQ.1) GOTO 100
2903         ENDIF
2904  
2905 C...Check that no odd resonance left undecayed.
2906         IF(MSTP(111).GE.1) THEN
2907           NFIX=N
2908           DO 130 I=MINT(84)+1,NFIX
2909             IF(K(I,1).GE.1.AND.K(I,1).LE.10.AND.K(I,2).NE.21.AND.
2910      &      K(I,2).NE.22) THEN
2911               KCA=PYCOMP(K(I,2))
2912               IF(MWID(KCA).NE.0.AND.MDCY(KCA,1).GE.1) THEN
2913                 CALL PYRESD(I)
2914                 IF(MINT(51).EQ.1) GOTO 100
2915               ENDIF
2916             ENDIF
2917   130     CONTINUE
2918         ENDIF
2919  
2920 C...Boost hadronic subsystem to overall rest frame.
2921 C..(Only relevant when photon inside lepton beam.)
2922         IF(MINT(141).NE.0.OR.MINT(142).NE.0) CALL PYGAGA(4,WTGAGA)
2923  
2924 C...Recalculate energies from momenta and masses (if desired).
2925         IF(MSTP(113).GE.1) THEN
2926           DO 140 I=MINT(83)+1,N
2927             IF(K(I,1).GT.0.AND.K(I,1).LE.10) P(I,4)=SQRT(P(I,1)**2+
2928      &      P(I,2)**2+P(I,3)**2+P(I,5)**2)
2929   140     CONTINUE
2930           NRECAL=N
2931         ENDIF
2932  
2933 C...Rearrange partons along strings, check invariant mass cuts.
2934         MSTU(28)=0
2935         IF(MSTP(111).LE.0) MSTJ(14)=-1
2936         CALL PYPREP(MINT(84)+1)
2937         MSTJ(14)=MSTJ14
2938         IF(MSTP(112).EQ.1.AND.MSTU(28).EQ.3) GOTO 100
2939         IF(MSTP(125).EQ.0.OR.MSTP(125).EQ.1) THEN
2940           DO 170 I=MINT(84)+1,N
2941             IF(K(I,2).EQ.94) THEN
2942               DO 160 I1=I+1,MIN(N,I+3)
2943                 IF(K(I1,3).EQ.I) THEN
2944                   K(I1,3)=MOD(K(I1,4)/MSTU(5),MSTU(5))
2945                   IF(K(I1,3).EQ.0) THEN
2946                     DO 150 II=MINT(84)+1,I-1
2947                         IF(K(II,2).EQ.K(I1,2)) THEN
2948                           IF(MOD(K(II,4),MSTU(5)).EQ.I1.OR.
2949      &                    MOD(K(II,5),MSTU(5)).EQ.I1) K(I1,3)=II
2950                         ENDIF
2951   150               CONTINUE
2952                     IF(K(I+1,3).EQ.0) K(I+1,3)=K(I,3)
2953                   ENDIF
2954                 ENDIF
2955   160         CONTINUE
2956             ENDIF
2957   170     CONTINUE
2958           CALL PYEDIT(12)
2959           CALL PYEDIT(14)
2960           IF(MSTP(125).EQ.0) CALL PYEDIT(15)
2961           IF(MSTP(125).EQ.0) MINT(4)=0
2962           DO 190 I=MINT(83)+1,N
2963             IF(K(I,1).EQ.11.AND.K(I,4).EQ.0.AND.K(I,5).EQ.0) THEN
2964               DO 180 I1=I+1,N
2965                 IF(K(I1,3).EQ.I.AND.K(I,4).EQ.0) K(I,4)=I1
2966                 IF(K(I1,3).EQ.I) K(I,5)=I1
2967   180         CONTINUE
2968             ENDIF
2969   190     CONTINUE
2970         ENDIF
2971  
2972 C...Introduce separators between sections in PYLIST event listing.
2973         IF(IPILE.EQ.1.AND.MSTP(125).LE.0) THEN
2974           MSTU70=1
2975           MSTU(71)=N
2976         ELSEIF(IPILE.EQ.1) THEN
2977           MSTU70=3
2978           MSTU(71)=2
2979           MSTU(72)=MINT(4)
2980           MSTU(73)=N
2981         ENDIF
2982  
2983 C...Go back to lab frame (needed for vertices, also in fragmentation).
2984         CALL PYFRAM(1)
2985  
2986 C...Set nonvanishing production vertex (optional).
2987         IF(MSTP(151).EQ.1) THEN
2988           DO 200 J=1,4
2989             VTX(J)=PARP(150+J)*SQRT(-2D0*LOG(MAX(1D-10,PYR(0))))*
2990      &      SIN(PARU(2)*PYR(0))
2991   200     CONTINUE
2992           DO 220 I=MINT(83)+1,N
2993             DO 210 J=1,4
2994               V(I,J)=V(I,J)+VTX(J)
2995   210       CONTINUE
2996   220     CONTINUE
2997         ENDIF
2998  
2999 C...Perform hadronization (if desired).
3000         IF(MSTP(111).GE.1) THEN
3001           CALL PYEXEC
3002           IF(MSTU(24).NE.0) GOTO 100
3003         ENDIF
3004         IF(MSTP(113).GE.1) THEN
3005           DO 230 I=NRECAL,N
3006             IF(P(I,5).GT.0D0) P(I,4)=SQRT(P(I,1)**2+
3007      &      P(I,2)**2+P(I,3)**2+P(I,5)**2)
3008   230     CONTINUE
3009         ENDIF
3010         IF(MSTP(125).EQ.0.OR.MSTP(125).EQ.1) CALL PYEDIT(14)
3011  
3012 C...Store event information and calculate Monte Carlo estimates of
3013 C...subprocess cross-sections.
3014   240   IF(IPILE.EQ.1) CALL PYDOCU
3015  
3016 C...Set counters for current pileup event and loop to next one.
3017         MSTI(41)=IPILE
3018         IF(IPILE.GE.2.AND.IPILE.LE.10) MSTI(40+IPILE)=ISUB
3019         IF(MSTU70.LT.10) THEN
3020           MSTU70=MSTU70+1
3021           MSTU(70+MSTU70)=N
3022         ENDIF
3023         MINT(83)=N
3024         MINT(84)=N+MSTP(126)
3025         IF(IPILE.LT.NPILE) CALL PYFRAM(2)
3026   250 CONTINUE
3027  
3028 C...Generic information on pileup events. Reconstruct missing history.
3029       IF(MSTP(131).EQ.1.AND.MSTP(133).GE.1) THEN
3030         PARI(91)=VINT(132)
3031         PARI(92)=VINT(133)
3032         PARI(93)=VINT(134)
3033         IF(MSTP(133).GE.2) PARI(93)=PARI(93)*XSEC(0,3)/VINT(131)
3034       ENDIF
3035       CALL PYEDIT(16)
3036  
3037 C...Transform to the desired coordinate frame.
3038   260 CALL PYFRAM(MSTP(124))
3039       MSTU(70)=MSTU70
3040       PARU(21)=VINT(1)
3041
3042 C...Error messages
3043  5100 FORMAT(1X,'Error: no subprocess switched on.'/
3044      &1X,'Execution stopped.')
3045  
3046       RETURN
3047       END
3048  
3049 C***********************************************************************
3050  
3051 C...PYSTAT
3052 C...Prints out information about cross-sections, decay widths, branching
3053 C...ratios, kinematical limits, status codes and parameter values.
3054  
3055       SUBROUTINE PYSTAT(MSTAT)
3056  
3057 C...Double precision and integer declarations.
3058       IMPLICIT DOUBLE PRECISION(A-H, O-Z)
3059       IMPLICIT INTEGER(I-N)
3060       INTEGER PYK,PYCHGE,PYCOMP
3061 C...Parameter statement to help give large particle numbers.
3062       PARAMETER (KSUSY1=1000000,KSUSY2=2000000,KTECHN=3000000,
3063      &KEXCIT=4000000,KDIMEN=5000000)
3064       PARAMETER (EPS=1D-3)
3065 C...Commonblocks.
3066       COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
3067       COMMON/PYDAT2/KCHG(500,4),PMAS(500,4),PARF(2000),VCKM(4,4)
3068       COMMON/PYDAT3/MDCY(500,3),MDME(8000,2),BRAT(8000),KFDP(8000,5)
3069       COMMON/PYSUBS/MSEL,MSELPD,MSUB(500),KFIN(2,-40:40),CKIN(200)
3070       COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
3071       COMMON/PYINT1/MINT(400),VINT(400)
3072       COMMON/PYINT2/ISET(500),KFPR(500,2),COEF(500,20),ICOL(40,4,2)
3073       COMMON/PYINT4/MWID(500),WIDS(500,5)
3074       COMMON/PYINT5/NGENPD,NGEN(0:500,3),XSEC(0:500,3)
3075       COMMON/PYINT6/PROC(0:500)
3076       CHARACTER PROC*28, CHTMP*16
3077       COMMON/PYMSSM/IMSS(0:99),RMSS(0:99)
3078       COMMON/PYMSRV/RVLAM(3,3,3), RVLAMP(3,3,3), RVLAMB(3,3,3)
3079       SAVE /PYDAT1/,/PYDAT2/,/PYDAT3/,/PYSUBS/,/PYPARS/,/PYINT1/,
3080      &/PYINT2/,/PYINT4/,/PYINT5/,/PYINT6/,/PYMSSM/,/PYMSRV/
3081 C...Local arrays, character variables and data.
3082       DIMENSION WDTP(0:300),WDTE(0:300,0:5),NMODES(0:20),PBRAT(10)
3083       CHARACTER PROGA(6)*28,CHAU*16,CHKF*16,CHD1*16,CHD2*16,CHD3*16,
3084      &CHIN(2)*12,STATE(-1:5)*4,CHKIN(21)*18,DISGA(2)*28,
3085      &PROGG9(13)*28,PROGG4(4)*28,PROGG2(2)*28,PROGP4(4)*28
3086       CHARACTER*24 CHD0, CHDC(10)
3087       CHARACTER*6 DNAME(3)
3088       DATA PROGA/
3089      &'VMD/hadron * VMD            ','VMD/hadron * direct         ',
3090      &'VMD/hadron * anomalous      ','direct * direct             ',
3091      &'direct * anomalous          ','anomalous * anomalous       '/
3092       DATA DISGA/'e * VMD','e * anomalous'/
3093       DATA PROGG9/
3094      &'direct * direct             ','direct * VMD                ',
3095      &'direct * anomalous          ','VMD * direct                ',
3096      &'VMD * VMD                   ','VMD * anomalous             ',
3097      &'anomal