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