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