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