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