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