G4RPGNucleonInelastic.cc

Go to the documentation of this file.
00001 //
00002 // ********************************************************************
00003 // * License and Disclaimer                                           *
00004 // *                                                                  *
00005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
00006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
00007 // * conditions of the Geant4 Software License,  included in the file *
00008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
00009 // * include a list of copyright holders.                             *
00010 // *                                                                  *
00011 // * Neither the authors of this software system, nor their employing *
00012 // * institutes,nor the agencies providing financial support for this *
00013 // * work  make  any representation or  warranty, express or implied, *
00014 // * regarding  this  software system or assume any liability for its *
00015 // * use.  Please see the license in the file  LICENSE  and URL above *
00016 // * for the full disclaimer and the limitation of liability.         *
00017 // *                                                                  *
00018 // * This  code  implementation is the result of  the  scientific and *
00019 // * technical work of the GEANT4 collaboration.                      *
00020 // * By using,  copying,  modifying or  distributing the software (or *
00021 // * any work based  on the software)  you  agree  to acknowledge its *
00022 // * use  in  resulting  scientific  publications,  and indicate your *
00023 // * acceptance of all terms of the Geant4 Software license.          *
00024 // ********************************************************************
00025 //
00026 // $Id$
00027 //
00028  
00029 #include "G4RPGNucleonInelastic.hh"
00030 #include "G4SystemOfUnits.hh"
00031 #include "Randomize.hh"
00032 
00033 G4RPGNucleonInelastic::G4RPGNucleonInelastic(const G4String& modelName)
00034  :G4RPGInelastic(modelName)
00035 {
00036   SetMinEnergy( 0.0 );
00037   SetMaxEnergy( 30.*GeV );
00038 
00039   // Initialize t1_dSigma_dMult, t0_dSigma_dMult,
00040   //   nucleon-nucleon inelastic cross sections for a given multiplicity 
00041   //   for |T_z| = 1 and 0, respectively 
00042 
00043   G4int i, j, k;
00044   G4int start, stop;
00045 
00046   for (j = 0; j < 8; j++) {
00047     start = pPindex[j][0];
00048     stop = pPindex[j][1] + 1;
00049     for (k = 0; k < 30; k++) {
00050       t1_dSigma_dMult[j][k] = 0.0;
00051       for (i = start; i < stop; i++) t1_dSigma_dMult[j][k] += pPCrossSections[i][k];
00052     }
00053 
00054     start = pNindex[j][0];
00055     stop = pNindex[j][1] + 1;
00056     for (k = 0; k < 30; k++) {
00057       t0_dSigma_dMult[j][k] = 0.0;
00058       for (i = start; i < stop; i++) t0_dSigma_dMult[j][k] += pNCrossSections[i][k];
00059     }
00060   }
00061 
00062   // Initialize total cross section array
00063 
00064   for (k = 0; k < 30; k++) {
00065     pPtot[k] = 0.0;
00066     pNtot[k] = 0.0;
00067     for (j = 0; j < 8; j++) {
00068       pPtot[k] += t1_dSigma_dMult[j][k];
00069       pNtot[k] += t0_dSigma_dMult[j][k];
00070     }
00071   }
00072   
00073   //  printCrossSections();
00074 }
00075 
00076 /*
00077 void G4RPGNucleonInelastic::printCrossSections() const
00078 {
00079   G4cout << " pp total cross sections (mb) " << G4endl;
00080   for (G4int i = 0; i < 5; i++) {
00081     G4int istart = i*6;
00082     G4int istop = istart + 6;
00083     for (G4int t = istart; t < istop; t++) G4cout << pPtot[t] << "  " ;
00084     G4cout << G4endl;
00085   }
00086   G4cout << G4endl;
00087 
00088   G4cout << " pn total cross sections (mb) " << G4endl;
00089   for (G4int i = 0; i < 5; i++) {
00090     G4int istart = i*6;
00091     G4int istop = istart + 6;
00092     for (G4int t = istart; t < istop; t++) G4cout << pNtot[t] << "  " ;
00093     G4cout << G4endl;
00094   }
00095 }
00096 */
00097 
00098 
00099 G4int G4RPGNucleonInelastic::GetMultiplicityT0(G4double KE) const
00100 {
00101   G4double multint(0.);
00102   std::vector<G4double> sigma;
00103 
00104   std::pair<G4int, G4double> epair = interpolateEnergy(KE);
00105   G4int k = epair.first;
00106   G4double fraction = epair.second;
00107 
00108   for(G4int j = 0; j < 8; j++) {
00109     multint = t0_dSigma_dMult[j][k]
00110          + fraction*(t0_dSigma_dMult[j][k+1] - t0_dSigma_dMult[j][k]);
00111       sigma.push_back(multint);
00112   }
00113 
00114   return sampleFlat(sigma) + 2;
00115 }
00116 
00117 
00118 G4int G4RPGNucleonInelastic::GetMultiplicityT1(G4double KE) const
00119 {
00120   G4double multint(0.);
00121   std::vector<G4double> sigma;
00122 
00123   std::pair<G4int, G4double> epair = interpolateEnergy(KE);
00124   G4int k = epair.first;
00125   G4double fraction = epair.second;
00126 
00127   for(G4int j = 0; j < 8; j++) {
00128     multint = t1_dSigma_dMult[j][k]
00129          + fraction*(t1_dSigma_dMult[j][k+1] - t1_dSigma_dMult[j][k]);
00130       sigma.push_back(multint);
00131   }
00132 
00133   return sampleFlat(sigma) + 2;
00134 }
00135 
00136 
00137 std::vector<G4int>
00138 G4RPGNucleonInelastic::GetFSPartTypesForT0(G4int mult, G4double KE) const
00139 {
00140   G4int i;
00141   G4double sigint(0.);
00142   std::vector<G4double> sigma;
00143 
00144   std::pair<G4int, G4double> epair = interpolateEnergy(KE);
00145   G4int k = epair.first;
00146   G4double fraction = epair.second;
00147  
00148   G4int start = pNindex[mult-2][0];
00149   G4int stop = pNindex[mult-2][1];
00150   
00151   for(i = start; i < stop; i++) {
00152       sigint = pNCrossSections[i][k]
00153           + fraction*(pNCrossSections[i][k+1] - pNCrossSections[i][k]);
00154       sigma.push_back(sigint);
00155   }
00156  
00157   G4int channel = sampleFlat(sigma);
00158   
00159   std::vector<G4int> kinds;
00160 
00161   if (mult == 2) {
00162     for(i = 0; i < mult; i++) kinds.push_back(T0_2bfs[channel][i]);
00163   } else if (mult == 3) {
00164     for(i = 0; i < mult; i++) kinds.push_back(T0_3bfs[channel][i]);
00165   } else if (mult == 4) {
00166     for(i = 0; i < mult; i++) kinds.push_back(T0_4bfs[channel][i]);
00167   } else if (mult == 5) {
00168     for(i = 0; i < mult; i++) kinds.push_back(T0_5bfs[channel][i]);
00169   } else if (mult == 6) {
00170     for(i = 0; i < mult; i++) kinds.push_back(T0_6bfs[channel][i]);
00171   } else if (mult == 7) {
00172     for(i = 0; i < mult; i++) kinds.push_back(T0_7bfs[channel][i]);
00173   } else if (mult == 8) {
00174     for(i = 0; i < mult; i++) kinds.push_back(T0_8bfs[channel][i]);
00175   } else if (mult == 9) {
00176     for(i = 0; i < mult; i++) kinds.push_back(T0_9bfs[channel][i]);
00177   } else {
00178     G4cout << " Illegal multiplicity " << G4endl;
00179   }
00180   
00181   return kinds;
00182 }
00183 
00184 
00185 std::vector<G4int>
00186 G4RPGNucleonInelastic::GetFSPartTypesForT1(G4int mult, G4double KE, 
00187                                                        G4int tzindex) const
00188 {
00189   G4int i;
00190   G4double sigint(0.);
00191   std::vector<G4double> sigma;
00192 
00193   std::pair<G4int, G4double> epair = interpolateEnergy(KE);
00194   G4int k = epair.first;
00195   G4double fraction = epair.second;
00196   
00197   G4int start = pPindex[mult-2][0];
00198   G4int stop = pPindex[mult-2][1];
00199   
00200   for(i = start; i < stop; i++) {
00201       sigint = pPCrossSections[i][k]
00202           + fraction*(pPCrossSections[i][k+1] - pPCrossSections[i][k]);
00203       sigma.push_back(sigint);
00204   }
00205  
00206   G4int channel = sampleFlat(sigma);
00207  
00208   std::vector<G4int> kinds;
00209   
00210   if (mult == 2) {
00211     for(i = 0; i < mult; i++) kinds.push_back(T1_2bfs[tzindex][channel][i]);
00212   } else if (mult == 3) {
00213     for(i = 0; i < mult; i++) kinds.push_back(T1_3bfs[tzindex][channel][i]);
00214   } else if (mult == 4) {
00215     for(i = 0; i < mult; i++) kinds.push_back(T1_4bfs[tzindex][channel][i]);
00216   } else if (mult == 5) {
00217     for(i = 0; i < mult; i++) kinds.push_back(T1_5bfs[tzindex][channel][i]);
00218   } else if (mult == 6) {
00219     for(i = 0; i < mult; i++) kinds.push_back(T1_6bfs[tzindex][channel][i]);
00220   } else if (mult == 7) {
00221     for(i = 0; i < mult; i++) kinds.push_back(T1_7bfs[tzindex][channel][i]);
00222   } else if (mult == 8) {
00223     for(i = 0; i < mult; i++) kinds.push_back(T1_8bfs[tzindex][channel][i]);
00224   } else if (mult == 9) {
00225     for(i = 0; i < mult; i++) kinds.push_back(T1_9bfs[tzindex][channel][i]);
00226   } else {
00227     G4cout << " Illegal multiplicity " << G4endl;
00228   }
00229   
00230   return kinds;
00231 } 
00232 
00234 //                                                                           //
00235 //   p p and n n (|Tz| = 1) cross sections                                   //
00236 //   and final state particle types                                          //
00237 //                                                                           //
00239 
00240 // Total p p cross section as a function of kinetic energy
00241 G4double G4RPGNucleonInelastic::pPtot[30];
00242 
00243 // p p multiplicities as a function of kinetic energy
00244 G4double G4RPGNucleonInelastic::t1_dSigma_dMult[8][30];
00245 
00246 const G4int G4RPGNucleonInelastic::pPindex[8][2] =
00247  {{0, 0}, {1, 6}, {7,24}, {25,56}, {57,63}, {64,71}, {72,81}, {82,92}};  
00248 
00249 // Outgoing particle types of a given multiplicity
00250 // T1_nbfs = final state types for p p and n n
00251 
00252 const G4int G4RPGNucleonInelastic::T1_2bfs[2][1][2] =
00253   {{{pro,pro}},
00254 
00255    {{neu,neu}}};
00256 
00257 const G4int G4RPGNucleonInelastic::T1_3bfs[2][6][3] =
00258   {{{pro,pro,pi0}, {pro,neu,pip}, {pro,lam,kp}, 
00259     {pro,s0,kp},   {pro,sp,k0},   {neu,sp,kp}},
00260 
00261    {{neu,neu,pi0}, {pro,neu,pim}, {neu,lam,k0}, 
00262     {neu,s0,k0},   {neu,sm,kp},   {pro,sm,k0}}};
00263 
00264 const G4int G4RPGNucleonInelastic::T1_4bfs[2][18][4] =
00265   {{{pro,pro,pip,pim},{pro,neu,pip,pi0},{pro,pro,pi0,pi0}, 
00266     {neu,neu,pip,pip},{pro,lam,kp,pi0}, {pro,lam,k0,pip}, 
00267     {neu,lam,kp,pip}, {neu,s0,kp,pip},  {pro,s0,kp,pi0},   
00268     {pro,s0,k0,pip},  {pro,sm,kp,pip},  {pro,sp,k0,pi0}, 
00269     {neu,sp,k0,pip},  {pro,sp,kp,pim},  {neu,sp,kp,pi0},  
00270     {pro,pro,k0,k0b}, {pro,pro,kp,km},  {pro,neu,kp,k0b}},
00271 
00272    {{neu,neu,pip,pim},{pro,neu,pim,pi0},{neu,neu,pi0,pi0},
00273     {pro,pro,pim,pim},{neu,lam,k0,pi0}, {neu,lam,kp,pim},
00274     {pro,lam,k0,pim}, {pro,s0,k0,pim},  {neu,s0,k0,pi0}, 
00275     {neu,s0,kp,pim},  {neu,sp,k0,pim},  {neu,sm,kp,pi0}, 
00276     {pro,sm,kp,pim},  {neu,sm,k0,pip},  {pro,sm,k0,pi0},
00277     {neu,neu,kp,km},  {neu,neu,k0,k0b}, {pro,neu,k0,km}}};
00278 
00279 const G4int G4RPGNucleonInelastic::T1_5bfs[2][32][5] =
00280   {{{pro,pro,pip,pim,pi0},{pro,pro,pi0,pi0,pi0},{pro,neu,pip,pip,pim},
00281     {pro,neu,pip,pi0,pi0},{neu,neu,pip,pip,pi0},{pro,lam,kp,pip,pim},
00282     {pro,lam,kp,pi0,pi0}, {pro,lam,k0,pip,pi0}, {pro,s0,kp,pip,pim},
00283     {pro,s0,kp,pi0,pi0},  {pro,s0,k0,pip,pi0},  {pro,sp,k0,pip,pim}, 
00284     {pro,sp,k0,pi0,pi0},  {pro,sp,kp,pim,pi0},  {pro,sm,kp,pip,pi0},
00285     {pro,sm,k0,pip,pip},  {neu,lam,kp,pip,pi0}, {neu,lam,k0,pip,pip},
00286     {neu,s0,kp,pip,pi0},  {neu,s0,k0,pip,pip},  {neu,sp,k0,pip,pi0},
00287     {neu,sp,kp,pip,pim},  {neu,sp,kp,pi0,pi0},  {neu,sm,kp,pip,pip},
00288     {pro,pro,pip,k0,km},  {pro,pro,pim,kp,k0b}, {pro,pro,pi0,k0,k0b},
00289     {pro,pro,pi0,kp,km},  {pro,neu,pip,k0,k0b}, {pro,neu,pip,kp,km}, 
00290     {pro,neu,pi0,kp,k0b}, {neu,neu,pip,kp,k0b}},
00291 
00292    {{neu,neu,pip,pim,pi0},{neu,neu,pi0,pi0,pi0},{pro,neu,pip,pim,pim},
00293     {pro,neu,pim,pi0,pi0},{pro,pro,pim,pim,pi0},{neu,lam,k0,pip,pim},
00294     {neu,lam,k0,pi0,pi0}, {neu,lam,kp,pim,pi0}, {neu,s0,k0,pip,pim},
00295     {neu,s0,k0,pi0,pi0},  {neu,s0,kp,pim,pi0},  {neu,sm,kp,pip,pim},
00296     {neu,sm,kp,pi0,pi0},  {neu,sm,k0,pip,pi0},  {neu,sp,k0,pim,pi0},
00297     {neu,sp,kp,pim,pim},  {pro,lam,k0,pim,pi0}, {pro,lam,kp,pim,pim},
00298     {pro,s0,k0,pim,pi0},  {pro,s0,kp,pim,pim},  {pro,sm,kp,pim,pi0},
00299     {pro,sm,k0,pip,pim},  {pro,sm,k0,pi0,pi0},  {pro,sp,k0,pim,pim},
00300     {neu,neu,pim,kp,k0b}, {neu,neu,pip,k0,km},  {neu,neu,pi0,kp,km},
00301     {neu,neu,pi0,k0,k0b}, {pro,neu,pim,kp,km},  {pro,neu,pim,k0,k0b},
00302     {pro,neu,pi0,k0,km},  {pro,pro,pim,k0,km}}};
00303 
00304 const G4int G4RPGNucleonInelastic::T1_6bfs[2][7][6] =
00305   {{{pro,pro,pip,pip,pim,pim},{pro,pro,pip,pim,pi0,pi0},
00306     {pro,pro,pi0,pi0,pi0,pi0},{pro,neu,pip,pip,pim,pi0},
00307     {pro,neu,pip,pi0,pi0,pi0},{neu,neu,pip,pip,pip,pim},
00308     {neu,neu,pip,pip,pi0,pi0}},
00309 
00310    {{neu,neu,pip,pip,pim,pim},{neu,neu,pip,pim,pi0,pi0},
00311     {neu,neu,pi0,pi0,pi0,pi0},{pro,neu,pip,pim,pim,pi0},
00312     {pro,neu,pim,pi0,pi0,pi0},{pro,pro,pip,pim,pim,pim},
00313     {pro,pro,pim,pim,pi0,pi0}}};
00314 
00315 const G4int G4RPGNucleonInelastic::T1_7bfs[2][8][7] =
00316   {{{pro,pro,pip,pip,pim,pim,pi0},{pro,pro,pip,pim,pi0,pi0,pi0}, 
00317     {pro,pro,pi0,pi0,pi0,pi0,pi0},{pro,neu,pip,pip,pip,pim,pim}, 
00318     {pro,neu,pip,pip,pim,pi0,pi0},{pro,neu,pip,pi0,pi0,pi0,pi0}, 
00319     {neu,neu,pip,pip,pip,pim,pi0},{neu,neu,pip,pip,pi0,pi0,pi0}}, 
00320 
00321    {{neu,neu,pip,pip,pim,pim,pi0},{neu,neu,pip,pim,pi0,pi0,pi0},
00322     {neu,neu,pi0,pi0,pi0,pi0,pi0},{pro,neu,pip,pip,pim,pim,pim},
00323     {pro,neu,pip,pim,pim,pi0,pi0},{pro,neu,pim,pi0,pi0,pi0,pi0},
00324     {pro,pro,pip,pim,pim,pim,pi0},{pro,pro,pim,pim,pi0,pi0,pi0}}};
00325 
00326 const G4int G4RPGNucleonInelastic::T1_8bfs[2][10][8] =
00327   {{{pro,pro,pip,pip,pip,pim,pim,pim},{pro,pro,pip,pip,pim,pim,pi0,pi0}, 
00328     {pro,pro,pip,pim,pi0,pi0,pi0,pi0},{pro,pro,pi0,pi0,pi0,pi0,pi0,pi0}, 
00329     {pro,neu,pip,pip,pip,pim,pim,pi0},{pro,neu,pip,pip,pim,pi0,pi0,pi0}, 
00330     {pro,neu,pip,pi0,pi0,pi0,pi0,pi0},{neu,neu,pip,pip,pip,pip,pim,pim}, 
00331     {neu,neu,pip,pip,pip,pim,pi0,pi0},{neu,neu,pip,pip,pi0,pi0,pi0,pi0}}, 
00332 
00333    {{neu,neu,pip,pip,pip,pim,pim,pim},{neu,neu,pip,pip,pim,pim,pi0,pi0},
00334     {neu,neu,pip,pim,pi0,pi0,pi0,pi0},{neu,neu,pi0,pi0,pi0,pi0,pi0,pi0},
00335     {pro,neu,pip,pip,pim,pim,pim,pi0},{pro,neu,pip,pim,pim,pi0,pi0,pi0},
00336     {pro,neu,pim,pi0,pi0,pi0,pi0,pi0},{pro,pro,pip,pip,pim,pim,pim,pim},
00337     {pro,pro,pip,pim,pim,pim,pi0,pi0},{pro,pro,pim,pim,pi0,pi0,pi0,pi0}}};
00338 
00339 const G4int G4RPGNucleonInelastic::T1_9bfs[2][11][9] =
00340 {{{pro,pro,pip,pip,pip,pim,pim,pim,pi0},{pro,pro,pip,pip,pim,pim,pi0,pi0,pi0}, 
00341   {pro,pro,pip,pim,pi0,pi0,pi0,pi0,pi0},{pro,pro,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, 
00342   {pro,neu,pip,pip,pip,pip,pim,pim,pim},{pro,neu,pip,pip,pip,pim,pim,pi0,pi0}, 
00343   {pro,neu,pip,pip,pim,pi0,pi0,pi0,pi0},{pro,neu,pip,pi0,pi0,pi0,pi0,pi0,pi0}, 
00344   {neu,neu,pip,pip,pip,pip,pim,pim,pi0},{neu,neu,pip,pip,pip,pim,pi0,pi0,pi0}, 
00345   {neu,neu,pip,pim,pi0,pi0,pi0,pi0,pi0}}, 
00346 
00347  {{neu,neu,pip,pip,pip,pim,pim,pim,pi0},{neu,neu,pip,pip,pim,pim,pi0,pi0,pi0},
00348   {neu,neu,pip,pim,pi0,pi0,pi0,pi0,pi0},{neu,neu,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
00349   {pro,neu,pip,pip,pip,pim,pim,pim,pim},{pro,neu,pip,pip,pim,pim,pim,pi0,pi0},
00350   {pro,neu,pip,pim,pim,pi0,pi0,pi0,pi0},{pro,neu,pim,pi0,pi0,pi0,pi0,pi0,pi0},
00351   {pro,pro,pip,pip,pim,pim,pim,pim,pi0},{pro,pro,pip,pim,pim,pim,pi0,pi0,pi0},
00352   {pro,pro,pip,pim,pi0,pi0,pi0,pi0,pi0}}};
00353 
00354 //
00355 // Cross sections (in mb) for p p -> 2-9 body final states
00356 //
00357 // first index:      0: channels for mult = 2
00358 //                 1-6: channels for mult = 3
00359 //                7-24: channels for mult = 4
00360 //               25-56: channels for mult = 5
00361 //               57-63: channels for mult = 6
00362 //               64-71: channels for mult = 7
00363 //               72-81: channels for mult = 8
00364 //               82-92: channels for mult = 9
00365 //
00366 // second index: kinetic energy
00367 //
00368 
00369 const G4float G4RPGNucleonInelastic::pPCrossSections[93][30] = {
00370 //
00371 // multiplicity 2 (1 channel)
00372 //
00373 //  p p (n n)
00374  {  0.0,330.0,240.0,160.0,110.0, 85.0, 63.0, 44.0, 33.0, 28.0,
00375    25.0, 24.0, 23.0, 23.0, 26.3, 26.1, 25.0, 23.5, 21.0, 18.0,
00376    16.0, 14.3, 12.5, 11.2, 10.3,  9.6,  9.0,  8.5,  8.0,  7.7 },
00377 //
00378 // multiplicity 3 (6 channels)
00379 //
00380 //  p p pi0 (n n pi0)
00381  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00382     0.0,  0.0,  0.0,  0.0,  0.0,  1.4,  4.0,  4.3,  4.0,  4.0,
00383     3.6,  3.0,  2.8,  2.5,  1.7,  1.3,  1.1,  1.0,  0.9,  0.85 },
00384 
00385 //  p n pi+ (p n pi-)
00386  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00387     0.0,  0.0,  0.0,  0.0,  0.7,  4.5, 15.0, 19.1, 18.0, 16.0,
00388    13.0, 10.0,  8.2,  6.0,  4.3,  3.3,  2.6,  2.0,  1.65, 1.4 },
00389 
00390 //  p L K+ (n L K0)
00391  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00392     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.012,
00393     0.03, 0.06, 0.06, 0.055,0.05, 0.047,0.043,0.04, 0.037,0.033 },
00394 
00395 //  p S0 K+ (n S0 K0)
00396  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00397     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00398     0.006,0.02, 0.027,0.026,0.021,0.018,0.015,0.011,0.009,0.007 },
00399 
00400 //  p S+ K0 (n S- K+)
00401  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00402     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00403     0.013,0.025,0.03, 0.029,0.027,0.026,0.024,0.022,0.021,0.019 },
00404 
00405 //  n S+ K+ (p S- K0)
00406  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00407     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00408     0.015,0.06, 0.07, 0.065,0.05, 0.04, 0.033,0.026,0.02, 0.015 },
00409 //
00410 // multiplicity 4 (18 channels)
00411 //
00412 //  p p pi+ pi- (n n pi+ pi-)
00413  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00414     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.05, 0.6,  1.9,
00415     2.8,  3.0,  3.0,  2.8,  2.5,  2.1,  1.9,  1.6,  1.4,  1.2 },
00416 
00417 //  p n pi+ pi0 (p n pi- pi0)
00418  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00419     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.05, 0.6,  3.5,
00420     4.0,  3.9,  3.5,  3.1,  2.8,  2.4,  2.2,  1.9,  1.7,  1.5 },
00421 
00422 //  p p pi0 pi0 (n n pi0 pi0)
00423  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00424     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.02, 0.24, 0.76,
00425     1.1,  1.2,  1.2,  1.1,  1.0,  0.84, 0.76, 0.64, 0.56, 0.48 },
00426 
00427 //  n n pi+ pi+ (p p pi- pi-)
00428  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00429     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.02, 0.24, 1.4,
00430     1.6,  1.6,  1.4,  1.2,  1.1,  1.0,  0.88, 0.76, 0.68, 0.6 },
00431 
00432 //  L K+ p pi0 (L K0 n pi0)
00433  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00434     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00435     0.002,0.015,0.05, 0.06, 0.052,0.042,0.037,0.029,0.025,0.020 },
00436 
00437 //  L K0 p pi+ (L K+ n pi-)
00438  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00439     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00440     0.002,0.015,0.06, 0.086,0.09, 0.082,0.072,0.06, 0.051,0.043 },
00441 
00442 //  L K+ n pi+ (L K0 p pi-)
00443  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00444     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00445     0.002,0.01, 0.04, 0.08, 0.066,0.058,0.05, 0.04, 0.035,0.03 },
00446 
00447 //  S0 K+ n pi+ (S0 K0 p pi-)
00448  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00449     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00450     0.0,  0.003,0.02, 0.03, 0.03, 0.025,0.02, 0.015,0.011,0.01 },
00451 
00452 //  S0 K+ p pi0 (S0 K0 n pi0)
00453  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00454     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00455     0.0,  0.005,0.02, 0.025,0.022,0.02, 0.015,0.01, 0.008,0.007 },
00456 
00457 //  S0 K0 p pi+ (S0 K+ n pi-)
00458  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00459     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00460     0.0,  0.003,0.012,0.04, 0.037,0.03, 0.027,0.022,0.019,0.016 },
00461 
00462 //  S- K+ p pi+ (S+ K0 n pi-)
00463  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00464     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00465     0.0,  0.004,0.016,0.037,0.031,0.028,0.023,0.02, 0.017,0.014 },
00466 
00467 //  S+ K0 p pi0 (S- K+ n pi0)
00468  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00469     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00470     0.0,  0.003,0.007,0.021,0.021,0.019,0.017,0.014,0.012,0.01 },
00471 
00472 //  S+ K0 n pi+ (S- K+ p pi-)
00473  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00474     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00475     0.0,  0.003,0.009,0.036,0.052,0.043,0.038,0.03, 0.026,0.02 },
00476 
00477 //  S+ K+ p pi- (S- K0 n pi+)
00478  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00479     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00480     0.0,  0.004,0.012,0.038,0.037,0.03, 0.026,0.02, 0.017,0.014 },
00481 
00482 //  S+ K+ n pi0 (S- K0 p pi0)
00483  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00484     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00485     0.0,  0.004,0.012,0.038,0.037,0.03, 0.026,0.02, 0.017,0.014 },
00486 
00487 //  p p K0 K0bar (n n K+ K-)
00488  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00489     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00490     0.0,  0.001,0.003,0.009,0.02, 0.02, 0.017,0.014,0.012,0.009 },
00491 
00492 //  p p K+ K- (n n K0 K0bar)
00493  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00494     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00495     0.0,  0.001,0.003,0.009,0.02, 0.02, 0.017,0.014,0.012,0.009 },
00496 
00497 //  p n K+ K0bar (p n K0 K-)
00498  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00499     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00500     0.0,  0.002,0.007,0.029,0.024,0.02, 0.017,0.014,0.012,0.009 },
00501 //
00502 // multiplicity 5 (32 channels)
00503 //
00504 //  p p pi+ pi- pi0 (n n pi+ pi- pi0)
00505  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00506     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.01, 0.06,
00507     0.4,  1.1,  1.8,  2.4,  2.4,  2.2,  2.0,  1.7,  1.5,  1.3 },
00508 
00509 //  p p pi0 pi0 pi0 (n n pi0 pi0 pi0)
00510  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00511     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.003,0.018,
00512     0.12, 0.33, 0.54, 0.72, 0.72, 0.66, 0.6,  0.51, 0.45, 0.39 },
00513 
00514 //  p n pi+ pi+ pi- (p n pi+ pi- pi-)
00515  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00516     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.12, 0.26,
00517     0.7,  1.6,  2.4,  2.6,  2.3,  2.0,  1.8,  1.6,  1.4,  1.2 },
00518 
00519 //  p n pi+ pi0 pi0 (p n pi- pi0 pi0)
00520  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00521     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.036,0.078,
00522     0.21, 0.48, 0.72, 0.78, 0.69, 0.6,  0.54, 0.48, 0.42, 0.36 },
00523 
00524 //  n n pi+ pi+ pi0 (p p pi- pi- pi0)
00525  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00526     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.01, 0.036,
00527     0.24, 0.66, 1.08, 1.44, 1.44, 1.32, 1.2,  1.0,  0.9,  0.78 },
00528 
00529 //  p L K+ pi+ pi- (n L K0 pi+ pi-)
00530  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00531     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00532     0.0,  0.0,  0.002,0.010,0.040,0.045,0.040,0.035,0.030,0.020 },
00533 
00534 //  p L K+ pi0 pi0 (n L K0 pi0 pi0)
00535  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00536     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00537     0.0,  0.0,  0.001,0.005,0.020,0.022,0.020,0.017,0.015,0.010 },
00538 
00539 //  p L K0 pi+ pi0 (n L K+ pi- pi0)
00540  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00541     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00542     0.0,  0.0,  0.002,0.010,0.040,0.045,0.040,0.035,0.030,0.020 },
00543 
00544 //  p S0 K+ pi+ pi- (n S0 K0 pi+ pi-)
00545  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00546     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00547     0.0,  0.0,  0.001,0.007,0.030,0.035,0.030,0.028,0.021,0.017 },
00548 
00549 //  p S0 K+ pi0 pi0 (n S0 K0 pi0 pi0)
00550  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00551     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00552     0.0,  0.0,  0.001,0.004,0.015,0.017,0.015,0.014,0.011,0.009 },
00553 
00554 //  p S0 K0 pi+ pi0 (n S0 K+ pi- pi0)
00555  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00556     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00557     0.0,  0.0,  0.001,0.007,0.030,0.035,0.030,0.028,0.021,0.017 },
00558 
00559 //  p S+ K0 pi+ pi- (n S- K+ pi+ pi-)
00560  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00561     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00562     0.0,  0.001,0.002,0.004,0.018,0.040,0.033,0.028,0.021,0.017 },
00563 
00564 //  p S+ K0 pi0 pi0 (n S- K+ pi0 pi0)
00565  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00566     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00567     0.0,  0.0,  0.001,0.002,0.009,0.020,0.013,0.014,0.011,0.009 },
00568 
00569 //  p S+ K+ pi- pi0 (n S- K0 pi+ pi0)
00570  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00571     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00572     0.0,  0.001,0.002,0.004,0.018,0.040,0.033,0.028,0.021,0.017 },
00573 
00574 //  p S- K+ pi+ pi0 (n S+ K0 pi- pi0)
00575  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00576     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00577     0.0,  0.001,0.002,0.004,0.018,0.040,0.033,0.028,0.021,0.017 },
00578 
00579 //  p S- K0 pi+ pi+ (n S+ K+ pi- pi-)
00580  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00581     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00582     0.0,  0.0,  0.001,0.007,0.030,0.035,0.030,0.028,0.021,0.017 },
00583 
00584 //  n L K+ pi+ pi0 (p L K0 pi- pi0)
00585  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00586     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00587     0.0,  0.0,  0.002,0.011,0.042,0.039,0.030,0.022,0.018,0.014 },
00588 
00589 //  n L K0 pi+ pi+ (p L K+ pi- pi-)
00590  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00591     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00592     0.0,  0.0,  0.002,0.011,0.042,0.039,0.030,0.022,0.018,0.014 },
00593 
00594 //  n S0 K+ pi+ pi0 (p S0 K0 pi- pi0)
00595  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00596     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00597     0.0,  0.0,  0.001,0.005,0.021,0.020,0.015,0.011,0.009,0.007 },
00598 
00599 //  n S0 K0 pi+ pi+ (p S0 K+ pi- pi-)
00600  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00601     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00602     0.0,  0.0,  0.001,0.005,0.021,0.020,0.015,0.011,0.009,0.007 },
00603 
00604 //  n S+ K0 pi+ pi0 (p S- K+ pi- pi0)
00605  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00606     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00607     0.0,  0.0,  0.001,0.005,0.021,0.020,0.015,0.011,0.009,0.007 },
00608 
00609 //  n S+ K+ pi+ pi- (p S- K0 pi+ pi-)
00610  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00611     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00612     0.0,  0.0,  0.002,0.011,0.042,0.039,0.030,0.022,0.018,0.014 },
00613 
00614 //  n S+ K+ pi0 pi0 (p S- K0 pi0 pi0)
00615  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00616     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00617     0.0,  0.0,  0.001,0.005,0.021,0.020,0.015,0.011,0.009,0.007 },
00618 
00619 //  n S- K+ pi+ pi+ (p S+ K0 pi- pi-)
00620  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00621     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00622     0.0,  0.0,  0.001,0.005,0.021,0.020,0.015,0.011,0.009,0.007 },
00623 
00624 //  p p pi+ K0 K- (n n pi- K+ K0bar)
00625  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00626     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00627     0.0,  0.0,  0.002,0.004,0.012,0.040,0.055,0.045,0.036,0.030 },
00628 
00629 //  p p pi- K+ K0bar (n n pi+ K0 K-)
00630  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00631     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00632     0.0,  0.001,0.002,0.004,0.018,0.040,0.033,0.028,0.021,0.017 },
00633 
00634 //  p p pi0 K0 K0bar (n n pi0 K+ K-)
00635  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00636     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00637     0.0,  0.0,  0.002,0.004,0.012,0.040,0.055,0.045,0.036,0.030 },
00638 
00639 //  p p pi0 K+ K- (n n pi0 K0 K0bar)
00640  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00641     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00642     0.0,  0.0,  0.002,0.004,0.012,0.040,0.055,0.045,0.036,0.030 },
00643 
00644 //  p n pi+ K0 K0bar (p n pi- K+ K-)
00645  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00646     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00647     0.0,  0.0,  0.002,0.012,0.035,0.058,0.050,0.033,0.023,0.016 },
00648 
00649 //  p n pi+ K+ K- (p n pi- K0 K0bar)
00650  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00651     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00652     0.0,  0.0,  0.002,0.012,0.035,0.058,0.050,0.033,0.023,0.016 },
00653 
00654 //  p n pi0 K+ K0bar (p n pi0 K0 K-)
00655  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00656     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00657     0.0,  0.001,0.002,0.004,0.018,0.040,0.033,0.028,0.021,0.017 },
00658 
00659 //  n n pi+ K+ K0bar (p p pi- K0 K-)
00660  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00661     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00662     0.0,  0.001,0.002,0.004,0.018,0.040,0.033,0.028,0.021,0.017 },
00663 //
00664 // multiplicity 6 (7 channels)
00665 //
00666 //  p p pi+ pi+ pi- pi- (n n pi+ pi+ pi- pi-)
00667  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00668     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00669     0.06, 0.1,  0.18, 0.38, 0.49, 0.46, 0.43, 0.40, 0.38, 0.36 },
00670 
00671 //  p p pi+ pi- pi0 pi0 (n n pi+ pi- pi0 pi0)
00672  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00673     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00674     0.03, 0.05, 0.09, 0.19, 0.25, 0.23, 0.22, 0.2,  0.19, 0.18 },
00675 
00676 //  p p pi0 pi0 pi0 pi0 (n n pi0 pi0 pi0 pi0)
00677  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00678     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00679     0.01, 0.02, 0.05, 0.1,  0.13, 0.12, 0.11, 0.1,  0.1,  0.09 },
00680 
00681 //  p n pi+ pi+ pi- pi0 (p n pi+ pi- pi- pi0)
00682  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00683     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00684     0.06, 0.1,  0.18, 0.38, 0.49, 0.46, 0.43, 0.40, 0.38, 0.36 },
00685 
00686 //  p n pi+ pi0 pi0 pi0 (p n pi- pi0 pi0 pi0)
00687  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00688     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00689     0.03, 0.05, 0.09, 0.19, 0.25, 0.23, 0.22, 0.2,  0.19, 0.18 },
00690 
00691 //  n n pi+ pi+ pi+ pi- (p p pi+ pi- pi- pi-)
00692  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00693     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00694     0.03, 0.05, 0.09, 0.19, 0.25, 0.23, 0.22, 0.2,  0.19, 0.18 },
00695 
00696 //  n n pi+ pi+ pi0 pi0 (p p pi- pi- pi0 pi0)
00697  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00698     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00699     0.03, 0.05, 0.09, 0.19, 0.25, 0.23, 0.22, 0.2,  0.19, 0.18 },
00700 //
00701 // multiplicity 7 (8 channels)
00702 //
00703 //  p p pi+ pi+ pi- pi- pi0 (n n pi+ pi+ pi- pi- pi0)
00704  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00705     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00706     0.0,  0.0,  0.06, 0.17, 0.5,  0.7,  0.7,  0.69, 0.66, 0.62 },
00707 
00708 //  p p pi+ pi- pi0 pi0 pi0 (n n pi+ pi- pi0 pi0 pi0)
00709  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00710     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00711     0.0,  0.0,  0.036,0.096,0.30, 0.42, 0.42, 0.42, 0.40, 0.37 },
00712 
00713 //  p p pi0 pi0 pi0 pi0 pi0 (n n pi0 pi0 pi0 pi0 pi0)
00714  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00715     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00716     0.0,  0.0,  0.018,0.048,0.14, 0.20, 0.22, 0.20, 0.19, 0.18 },
00717 
00718 //  p n pi+ pi+ pi+ pi- pi- (p n pi+ pi+ pi- pi- pi-)
00719  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00720     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00721     0.0,  0.0,  0.06, 0.19, 0.31, 0.41, 0.44, 0.47, 0.45, 0.45 },
00722 
00723 //  p n pi+ pi+ pi- pi0 pi0 (p n pi+ pi- pi- pi0 pi0)
00724  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00725     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00726     0.0,  0.0,  0.036,0.12, 0.18, 0.24, 0.26, 0.23, 0.28, 0.26 },
00727 
00728 //  p n pi+ pi0 pi0 pi0 pi0 (p n pi- pi0 pi0 pi0 pi0)
00729  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00730     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00731     0.0,  0.0,  0.024,0.06, 0.08, 0.12, 0.13, 0.14, 0.13, 0.13 },
00732 
00733 //  n n pi+ pi+ pi+ pi- pi0 (p p pi+ pi- pi- pi- pi0)
00734  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00735     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00736     0.0,  0.0,  0.06, 0.17, 0.5,  0.7,  0.7,  0.69, 0.66, 0.62 },
00737 
00738 //  n n pi+ pi+ pi0 pi0 pi0 (p p pi- pi- pi0 pi0 pi0)
00739  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00740     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00741     0.0,  0.0,  0.036,0.096,0.30, 0.42, 0.42, 0.41, 0.40, 0.37 },
00742 //
00743 // multiplicity 8 (10 channels)
00744 //
00745 //  p p pi+ pi+ pi+ pi- pi- pi- (n n pi+ pi+ pi+ pi- pi- pi-)
00746  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00747     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00748     0.0,  0.003,0.01, 0.024,0.075,0.18,0.27,  0.30, 0.27, 0.24 },
00749 
00750 //  p p pi+ pi+ pi- pi- pi0 pi0 (n n pi+ pi+ pi- pi- pi0 pi0)
00751  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00752     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00753     0.0,  0.003,0.01, 0.024,0.075,0.18,0.27,  0.30, 0.27, 0.24 },
00754 
00755 //  p p pi+ pi- pi0 pi0 pi0 pi0 (n n pi+ pi- pi0 pi0 pi0 pi0)
00756  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00757     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00758     0.0,  0.003,0.006,0.015,0.045,0.12, 0.15, 0.18, 0.15, 0.15 },
00759 
00760 //  p p pi0 pi0 pi0 pi0 pi0 pi0 (n n pi0 pi0 pi0 pi0 pi0 pi0)
00761  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00762     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00763     0.0,  0.0,  0.003,0.01, 0.024,0.06, 0.09, 0.12, 0.09, 0.09 },
00764 
00765 //  p n pi+ pi+ pi+ pi- pi- pi0 (p n pi+ pi+ pi- pi- pi- pi0)
00766  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00767     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00768     0.0,  0.003,0.01, 0.024,0.075,0.18, 0.27, 0.30, 0.27, 0.24 },
00769 
00770 //  p n pi+ pi+ pi- pi0 pi0 pi0 (p n pi+ pi- pi- pi0 pi0 pi0)
00771  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00772     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00773     0.0,  0.003,0.006,0.015,0.045,0.12, 0.15, 0.18, 0.15, 0.15 },
00774 
00775 //  p n pi+ pi0 pi0 pi0 pi0 pi0 (p n pi- pi0 pi0 pi0 pi0 pi0)
00776  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00777     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00778     0.0,  0.0,  0.003,0.01, 0.027,0.06, 0.09, 0.12, 0.09, 0.09 },
00779 
00780 //  n n pi+ pi+ pi+ pi+ pi- pi- (p p pi+ pi+ pi- pi- pi- pi-)
00781  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00782     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00783     0.0,  0.003,0.01, 0.024,0.075,0.18, 0.27, 0.30, 0.27, 0.24 },
00784 
00785 //  n n pi+ pi+ pi+ pi- pi0 pi0 (p p pi+ pi- pi- pi- pi0 pi0)
00786  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00787     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00788     0.0,  0.003,0.006,0.015,0.045,0.12, 0.15, 0.18, 0.15, 0.15 },
00789 
00790 //  n n pi+ pi+ pi0 pi0 pi0 pi0 (p p pi- pi- pi0 pi0 pi0 pi0)
00791  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00792     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00793     0.0,  0.0,  0.003,0.01, 0.027,0.06, 0.09, 0.12, 0.09, 0.09 },
00794 //
00795 // multiplicity 9 (11 channels)
00796 //
00797 //  p p pi+ pi+ pi+ pi- pi- pi- pi0 (n n pi+ pi+ pi+ pi- pi- pi- pi0)
00798  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00799     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00800     0.0,  0.002,0.008,0.025,0.074,0.11, 0.14, 0.15, 0.15, 0.15 },
00801 
00802 //  p p pi+ pi+ pi- pi- pi0 pi0 pi0 (n n pi+ pi+ pi- pi- pi0 pi0 pi0)
00803  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00804     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00805     0.0,  0.001,0.006,0.020,0.058,0.091,0.11, 0.12, 0.12, 0.12 },
00806 
00807 //  p p pi+ pi- pi0 pi0 pi0 pi0 pi0 (n n pi+ pi- pi0 pi0 pi0 pi0 pi0)
00808  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00809     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00810     0.0,  0.001,0.004,0.012,0.035,0.055,0.065,0.07, 0.07, 0.07 },
00811 
00812 //  p p pi0 pi0 pi0 pi0 pi0 pi0 pi0 (n n pi0 pi0 pi0 pi0 pi0 pi0 pi0)
00813  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00814     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00815     0.0,  0.0,  0.0,  0.006,0.027,0.032,0.04, 0.042,0.042,0.042 },
00816 
00817 //  p n pi+ pi+ pi+ pi+ pi- pi- pi- (p n pi+ pi+ pi+ pi- pi- pi- pi-)
00818  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00819     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00820     0.0,  0.002,0.006,0.02, 0.06, 0.15, 0.19, 0.22, 0.22, 0.22 },
00821 
00822 //  p n pi+ pi+ pi+ pi- pi- pi0 pi0 (p n pi+ pi+ pi- pi- pi- pi0 pi0)
00823  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00824     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00825     0.0,  0.003,0.008,0.026,0.078,0.20, 0.25, 0.29, 0.29, 0.29 },
00826 
00827 //  p n pi+ pi+ pi- pi0 pi0 pi0 pi0 (p n pi+ pi- pi- pi0 pi0 pi0 pi0)
00828  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00829     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00830     0.0,  0.001,0.005,0.016,0.047,0.12, 0.15, 0.17, 0.17, 0.17 },
00831 
00832 //  p n pi+ pi0 pi0 pi0 pi0 pi0 pi0 (p n pi- pi0 pi0 pi0 pi0 pi0 pi0)
00833  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00834     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00835     0.0,  0.001,0.003,0.009,0.029,0.07, 0.094,0.10, 0.10, 0.10 },
00836 
00837 //  n n pi+ pi+ pi+ pi+ pi- pi- pi0 (p p pi+ pi+ pi- pi- pi- pi- pi0)
00838  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00839     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00840     0.0,  0.002,0.006,0.02, 0.06, 0.15, 0.19, 0.22, 0.22, 0.22 },
00841 
00842 //  n n pi+ pi+ pi+ pi- pi0 pi0 pi0 (p p pi+ pi- pi- pi- pi0 pi0 pi0)
00843  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00844     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00845     0.0,  0.001,0.005,0.015,0.047,0.12, 0.15, 0.17, 0.17, 0.17 },
00846 
00847 //  n n pi+ pi- pi0 pi0 pi0 pi0 pi0 (p p pi+ pi- pi0 pi0 pi0 pi0 pi0)
00848  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00849     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00850     0.0,  0.001,0.003,0.009,0.029,0.07, 0.094,0.10, 0.10, 0.10 }};
00851 
00853 //                                                                           //
00854 //   p n and n p (|Tz| = 0) cross sections                                   //
00855 //   and final state particle types                                          //
00856 //                                                                           //
00858 
00859 // Total p n cross section as a function of kinetic energy
00860 G4double G4RPGNucleonInelastic::pNtot[30];
00861 
00862 // p n multiplicities as a function of kinetic energy
00863 G4double G4RPGNucleonInelastic::t0_dSigma_dMult[8][30];
00864 
00865 const G4int G4RPGNucleonInelastic::pNindex[8][2] =
00866  {{0, 0}, {1,9}, {10,31}, {32,69}, {70,76}, {77,85}, {86,95}, {96,107}};  
00867 // first index:      0: channels for mult = 2
00868 //                 1-9: channels for mult = 3
00869 //               10-31: channels for mult = 4
00870 //               32-69: channels for mult = 5
00871 //               70-76: channels for mult = 6
00872 //               77-85: channels for mult = 7
00873 //               86-95: channels for mult = 8
00874 //              96-107: channels for mult = 9
00875 
00876 // Outgoing particle types of a given multiplicity
00877 // T0_nbfs = final state types for p n and n p
00878 
00879 const G4int G4RPGNucleonInelastic::T0_2bfs[1][2] =
00880   {{pro,neu}};
00881 
00882 const G4int G4RPGNucleonInelastic::T0_3bfs[9][3] =
00883   {{pro,pro,pim},{pro,neu,pi0},{neu,neu,pip},{pro,lam,k0}, 
00884    {pro,s0,k0},  {pro,sm,kp},  {neu,lam,kp}, {neu,s0,kp}, 
00885    {neu,sp,k0}};
00886 
00887 const G4int G4RPGNucleonInelastic::T0_4bfs[22][4] =
00888   {{pro,neu,pip,pim},{pro,pro,pim,pi0},{pro,neu,pi0,pi0}, 
00889    {neu,neu,pip,pi0},{pro,lam,kp,pim}, {pro,s0,kp,pim},
00890    {pro,lam,k0,pi0}, {pro,s0,k0,pi0},  {pro,sp,k0,pim},
00891    {pro,sm,kp,pi0},  {pro,sm,k0,pip},  {neu,lam,kp,pi0}, 
00892    {neu,lam,k0,pip}, {neu,sp,kp,pim},  {neu,sp,k0,pi0}, 
00893    {neu,s0,kp,pi0},  {neu,s0,k0,pip},  {neu,sm,kp,pip},
00894    {pro,neu,kp,km},  {pro,neu,k0,k0b}, {pro,pro,k0,km},
00895    {neu,neu,kp,k0b}};
00896 
00897 const G4int G4RPGNucleonInelastic::T0_5bfs[38][5] =
00898   {{pro,neu,pip,pim,pi0},{pro,neu,pi0,pi0,pi0},{pro,pro,pip,pim,pim}, 
00899    {pro,pro,pim,pi0,pi0},{neu,neu,pip,pip,pim},{neu,neu,pip,pi0,pi0}, 
00900    {pro,lam,kp,pim,pi0}, {pro,lam,k0,pip,pim}, {pro,lam,k0,pi0,pi0}, 
00901    {pro,s0,k0,pip,pim},  {pro,s0,k0,pi0,pi0},  {pro,s0,kp,pim,pi0}, 
00902    {pro,sp,kp,pim,pim},  {pro,sp,k0,pim,pi0},  {pro,sm,k0,pip,pi0}, 
00903    {pro,sm,kp,pip,pim},  {pro,sm,kp,pi0,pi0},  {neu,lam,kp,pip,pim},
00904    {neu,lam,kp,pi0,pi0}, {neu,lam,k0,pip,pi0}, {neu,s0,kp,pip,pim}, 
00905    {neu,s0,kp,pi0,pi0},  {neu,s0,k0,pip,pi0},  {neu,sp,k0,pip,pim}, 
00906    {neu,sp,k0,pi0,pi0},  {neu,sp,kp,pim,pi0},  {neu,sm,kp,pip,pi0}, 
00907    {neu,sm,k0,pip,pip},  {pro,neu,kp,km,pi0},  {pro,neu,k0,k0b,pi0}, 
00908    {pro,neu,k0,km,pip},  {pro,neu,kp,k0b,pim}, {pro,pro,k0,k0b,pim}, 
00909    {pro,pro,kp,km,pim},  {pro,pro,k0,km,pi0},  {neu,neu,kp,km,pip}, 
00910    {neu,neu,k0,k0b,pip}, {neu,neu,kp,k0b,pi0}}; 
00911 
00912 const G4int G4RPGNucleonInelastic::T0_6bfs[7][6] =
00913   {{pro,neu,pip,pip,pim,pim},{pro,neu,pip,pim,pi0,pi0}, 
00914    {pro,neu,pi0,pi0,pi0,pi0},{pro,pro,pip,pim,pim,pi0}, 
00915    {pro,pro,pim,pi0,pi0,pi0},{neu,neu,pip,pip,pim,pi0}, 
00916    {neu,neu,pip,pi0,pi0,pi0}}; 
00917 
00918 const G4int G4RPGNucleonInelastic::T0_7bfs[9][7] =
00919   {{pro,neu,pip,pip,pim,pim,pi0},{pro,neu,pip,pim,pi0,pi0,pi0}, 
00920    {pro,neu,pi0,pi0,pi0,pi0,pi0},{pro,pro,pip,pip,pim,pim,pim}, 
00921    {pro,pro,pip,pim,pim,pi0,pi0},{pro,pro,pim,pi0,pi0,pi0,pi0}, 
00922    {neu,neu,pip,pip,pip,pim,pim},{neu,neu,pip,pip,pim,pi0,pi0}, 
00923    {neu,neu,pip,pi0,pi0,pi0,pi0}}; 
00924 
00925 const G4int G4RPGNucleonInelastic::T0_8bfs[10][8] =
00926 {{pro,neu,pip,pip,pip,pim,pim,pim},{pro,neu,pip,pip,pim,pim,pi0,pi0}, 
00927  {pro,neu,pip,pim,pi0,pi0,pi0,pi0},{pro,neu,pi0,pi0,pi0,pi0,pi0,pi0}, 
00928  {pro,pro,pip,pip,pim,pim,pim,pi0},{pro,pro,pip,pim,pim,pi0,pi0,pi0}, 
00929  {pro,pro,pim,pi0,pi0,pi0,pi0,pi0},{neu,neu,pip,pip,pip,pim,pim,pi0}, 
00930  {neu,neu,pip,pip,pim,pi0,pi0,pi0},{neu,neu,pip,pi0,pi0,pi0,pi0,pi0}}; 
00931 
00932 const G4int G4RPGNucleonInelastic::T0_9bfs[12][9] =
00933 {{pro,neu,pip,pip,pip,pim,pim,pim,pi0},{pro,neu,pip,pip,pim,pim,pi0,pi0,pi0}, 
00934  {pro,neu,pip,pim,pi0,pi0,pi0,pi0,pi0},{pro,neu,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, 
00935  {pro,pro,pip,pip,pip,pim,pim,pim,pim},{pro,pro,pip,pip,pim,pim,pim,pi0,pi0}, 
00936  {pro,pro,pip,pim,pim,pi0,pi0,pi0,pi0},{pro,pro,pim,pi0,pi0,pi0,pi0,pi0,pi0}, 
00937  {neu,neu,pip,pip,pip,pip,pim,pim,pim},{neu,neu,pip,pip,pip,pim,pim,pi0,pi0}, 
00938  {neu,neu,pip,pip,pim,pi0,pi0,pi0,pi0},{neu,neu,pip,pi0,pi0,pi0,pi0,pi0,pi0}}; 
00939 
00940 //
00941 // Cross sections (in mb) for p n -> 2-9 body final states
00942 //
00943 // first index:      0: channels for mult = 2
00944 //                 1-9: channels for mult = 3
00945 //               10-31: channels for mult = 4
00946 //               32-69: channels for mult = 5
00947 //               70-76: channels for mult = 6
00948 //               77-85: channels for mult = 7
00949 //               86-95: channels for mult = 8
00950 //              96-107: channels for mult = 9
00951 //
00952 // second index: kinetic energy
00953 //
00954 const G4float G4RPGNucleonInelastic::pNCrossSections[108][30] = {
00955 //
00956 // multiplicity 2 (1 channel)
00957 //
00958 //  p n (p n)
00959  {  0.0, 46.0, 46.0, 46.0, 46.0, 46.0, 46.0, 46.0, 46.0, 46.0,
00960    44.0, 42.0, 40.0, 35.0, 31.0, 27.0, 23.0, 19.0, 17.0, 15.5,
00961    14.0, 13.0, 12.0, 11.0, 10.0,  9.5,  9.0,  8.5,  8.0,  7.7 },
00962 //
00963 // multiplicity 3 (9 channels)
00964 //
00965 //  p p pi- (n n pi+)
00966  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00967     0.0,  0.0,  0.0,  0.0,  0.25, 0.9,  1.75, 2.3,  2.8,  2.8,
00968     2.2,  1.9,  1.6,  1.35, 1.1,  0.95, 0.8,  0.7,  0.6,  0.53 },
00969 
00970 //  p n pi0 (p n pi0)
00971  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00972     0.0,  0.0,  0.0,  0.0,  1.8,  4.7,  8.3, 11.3, 12.0, 10.2,
00973     8.2,  6.0,  4.9,  3.6,  2.5,  2.0,  1.6,  1.2,  1.0,  0.08 },
00974 
00975 //  n n pi+ (p p pi-)
00976  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00977     0.0,  0.0,  0.0,  0.0,  0.95, 2.4,  4.2,  5.6,  6.1,  5.1,
00978     4.1,  3.0,  2.5,  1.8,  1.2,  1.0,  0.8,  0.6,  0.5,  0.41},
00979 
00980 //  p L K0 (n L K+)
00981  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00982     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00983     0.004,0.013,0.021,0.025,0.021,0.019,0.018,0.016,0.014,0.012},
00984 
00985 //  p S0 K0 (n S0 K+)
00986  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00987     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00988     0.003,0.010,0.016,0.020,0.016,0.015,0.014,0.013,0.011,0.010},
00989 
00990 //  p S- K+ (n S+ K0)
00991  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00992     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00993     0.004,0.013,0.021,0.025,0.021,0.019,0.018,0.016,0.014,0.012},
00994 
00995 //  n L K+ (p L K0)
00996  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00997     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00998     0.004,0.013,0.021,0.025,0.021,0.019,0.018,0.016,0.014,0.012},
00999 
01000 //  n S0 K+ (p S0 K0)
01001  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01002     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01003     0.003,0.010,0.016,0.020,0.016,0.015,0.014,0.013,0.011,0.010},
01004 
01005 //  n S+ K0 (p S- K+)
01006  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01007     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01008     0.004,0.013,0.021,0.025,0.021,0.019,0.018,0.016,0.014,0.012},
01009 //
01010 // multiplicity 4 (22 channels)
01011 //
01012 //  p n pi+ pi- (p n pi+ pi-)
01013  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01014     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.12, 0.38, 1.1,  3.5,
01015     5.9,  5.9,  5.1,  4.2,  3.7,  3.0,  2.6,  2.1,  1.8,  1.4 },
01016 
01017 //  p p pi- pi0 (n n pi+ pi0)
01018  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01019     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.03, 0.1,  0.24, 0.55,
01020     1.2,  1.5,  1.45, 1.25, 1.0,  0.9,  0.8,  0.7,  0.6,  0.53 },
01021 
01022 //  p n pi0 pi0 (p n pi0 pi0)
01023  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01024     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.07, 0.24, 0.66, 2.1,
01025     3.6,  3.6,  3.1,  2.5,  2.2,  1.8,  1.5,  1.2,  1.1,  0.84 },
01026 
01027 //  n n pi+ pi0 (p p pi- pi0)
01028  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01029     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.03, 0.1,  0.24, 0.55,
01030     1.2,  1.5,  1.45, 1.25, 1.0,  0.9,  0.8,  0.7,  0.6,  0.53 },
01031 
01032 //  p L K+ pi- (n L K0 pi+)
01033  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01034     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01035     0.0,  0.002,0.006,0.018,0.037,0.036,0.033,0.030,0.028,0.023 },
01036 
01037 //  p S0 K+ pi- (n S0 K0 pi+)
01038  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01039     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01040     0.0,  0.0,  0.002,0.007,0.021,0.025,0.022,0.020,0.018,0.017 },
01041 
01042 //  p L K0 pi0 (n L K+ pi0)
01043  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01044     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01045     0.0,  0.001,0.002,0.005,0.009,0.009,0.008,0.007,0.007,0.006 },
01046 
01047 //  p S0 K0 pi0 (n S0 K+ pi0)
01048  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01049     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01050     0.0,  0.0,  0.001,0.002,0.005,0.006,0.005,0.005,0.004,0.004 },
01051 
01052 //  p S+ K0 pi- (n S- K+ pi+)
01053  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01054     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01055     0.0,  0.0015,0.004,0.013,0.02,0.016,0.013,0.01, 0.009,0.007 },
01056 
01057 //  p S- K+ pi0 (n S+ K0 pi0)
01058  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01059     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01060     0.0,  0.003,0.008,0.026,0.042,0.042,0.035,0.029,0.023,0.018 },
01061 
01062 //  p S- K0 pi+ (n S+ K+ pi-)
01063  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01064     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01065     0.0,  0.003,0.008,0.026,0.042,0.042,0.035,0.029,0.023,0.018 },
01066 
01067 //  n L K+ pi0 (p L K0 pi0)
01068  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01069     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01070     0.0,  0.002,0.006,0.018,0.037,0.036,0.033,0.030,0.028,0.023 },
01071 
01072 //  n L K0 pi+ (p L K+ pi-)
01073  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01074     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01075     0.0,  0.001,0.002,0.005,0.009,0.009,0.008,0.007,0.007,0.006 },
01076 
01077 //  n S+ K+ pi- (p S- K0 pi+)
01078  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01079     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01080     0.0,  0.0,  0.002,0.007,0.021,0.025,0.022,0.020,0.018,0.017 },
01081 
01082 //  n S+ K0 pi0 (p S- K+ pi0)
01083  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01084     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01085     0.0,  0.0,  0.001,0.002,0.005,0.006,0.005,0.005,0.004,0.004 },
01086 
01087 //  n S0 K+ pi0 (p S0 K0 pi0)
01088  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01089     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01090     0.0,  0.003,0.008,0.026,0.042,0.042,0.035,0.029,0.023,0.018 },
01091 
01092 //  n S0 K0 pi+ (p S0 K+ pi-)
01093  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01094     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01095     0.0,  0.003,0.008,0.026,0.042,0.042,0.035,0.029,0.023,0.018 },
01096 
01097 //  n S- K+ pi+ (p S+ K0 pi-)
01098  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01099     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01100     0.0, 0.0015,0.004,0.013,0.02, 0.016,0.013,0.01, 0.009,0.007 },
01101 
01102 //  p n K+ K- (p n K0 K0bar)
01103  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01104     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01105     0.0,  0.0,  0.001,0.005,0.015,0.019,0.015,0.012,0.009,0.007 },
01106 
01107 //  p n K0 K0bar (p n K+ K-)
01108  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01109     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01110     0.0,  0.0,  0.001,0.005,0.015,0.019,0.015,0.012,0.009,0.007 },
01111 
01112 //  p p K0 K- (n n K+ K0bar)
01113  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01114     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01115     0.0,  0.0,  0.001,0.005,0.015,0.019,0.015,0.012,0.009,0.007 },
01116 
01117 //  n n K+ K0bar (p p K0 K-)
01118  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01119     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01120     0.0,  0.0,  0.001,0.005,0.015,0.019,0.015,0.012,0.009,0.007 },
01121 //
01122 // multiplicity 5 (38 channels)
01123 //
01124 //  p n pi+ pi- pi0 (p n pi+ pi- pi0)
01125  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01126     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.01, 0.04,
01127     0.3,  0.82, 1.35, 1.8,  1.8,  1.65, 1.5,  1.28, 1.12, 0.98 },
01128 
01129 //  p n pi0 pi0 pi0 (p n pi0 pi0 pi0)
01130  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01131     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.004,0.022,
01132     0.15, 0.41, 0.68, 0.9,  0.9,  0.82, 0.75, 0.64, 0.55, 0.49 },
01133 
01134 //  p p pi+ pi- pi- (n n pi+ pi+ pi-)
01135  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01136     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.09, 0.2,
01137     0.52, 1.2,  1.8,  2.0,  1.7,  1.5,  1.35, 1.2,  1.05, 0.9 },
01138 
01139 //  p p pi- pi0 pi0 (n n pi+ pi0 pi0)
01140  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01141     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.04, 0.1,
01142     0.26, 0.6,  0.9,  0.98, 0.86, 0.75, 0.68, 0.6,  0.52, 0.45 },
01143 
01144 //  n n pi+ pi+ pi- (p p pi+ pi- pi-)
01145  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01146     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.01, 0.04,
01147     0.3,  0.82, 1.35, 1.8,  1.8,  1.65, 1.5,  1.28, 1.12, 0.98 },
01148 
01149 //  n n pi+ pi0 pi0 (p p pi- pi0 pi0)
01150  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01151     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.003,0.02,
01152     0.15, 0.41, 0.68, 0.9,  0.9,  0.82, 0.75, 0.64, 0.56, 0.49 },
01153 
01154 //  p L K+ pi- pi0 (n L K0 pi+ pi0)
01155  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01156     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01157     0.0,  0.0,  0.002,0.009,0.023,0.025,0.022,0.018,0.015,0.013 },
01158 
01159 //  p L K0 pi+ pi- (n L K+ pi+ pi-)
01160  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01161     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01162     0.0,  0.0,  0.002,0.009,0.023,0.025,0.022,0.018,0.015,0.013 },
01163 
01164 //  p L K0 pi0 pi0 (n L K+ pi0 pi0)
01165  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01166     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01167     0.0,  0.0,  0.001,0.005,0.011,0.012,0.011,0.09, 0.07, 0.07 },
01168 
01169 //  p S0 K0 pi+ pi- (n S0 K+ pi+ pi-)
01170  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01171     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01172     0.0,  0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
01173 
01174 //  p S0 K0 pi0 pi0 (n S0 K+ pi0 pi0)
01175  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01176     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01177     0.0,  0.0,  0.001,0.003,0.007,0.008,0.007,0.006,0.005,0.004 },
01178 
01179 //  p S0 K+ pi- pi0 (n S0 K0 pi+ pi0)
01180  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01181     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01182     0.0,  0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
01183 
01184 //  p S+ K+ pi- pi- (n S- K0 pi+ pi+)
01185  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01186     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01187     0.0,  0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
01188 
01189 //  p S+ K0 pi- pi0 (n S- K+ pi+ pi0)
01190  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01191     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01192     0.0,  0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
01193 
01194 //  p S- K0 pi+ pi0 (n S+ K+ pi- pi0)
01195  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01196     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01197     0.0,  0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
01198 
01199 //  p S- K+ pi+ pi- (n S+ K0 pi+ pi-)
01200  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01201     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01202     0.0,  0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
01203 
01204 //  p S- K+ pi0 pi0 (n S+ K0 pi0 pi0)
01205  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01206     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01207     0.0,  0.0,  0.002,0.003,0.007,0.008,0.007,0.006,0.005,0.004 },
01208 
01209 //  n L K+ pi+ pi- (p L K0 pi+ pi-)
01210  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01211     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01212     0.0,  0.0,  0.002,0.009,0.023,0.025,0.022,0.018,0.015,0.013 },
01213 
01214 //  n L K+ pi0 pi0 (p L K0 pi0 pi0)
01215  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01216     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01217     0.0,  0.0,  0.001,0.004,0.011,0.012,0.011,0.009,0.007,0.006 },
01218 
01219 //  n L K0 pi+ pi0 (p L K+ pi- pi0)
01220  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01221     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01222     0.0,  0.0,  0.002,0.009,0.023,0.025,0.022,0.018,0.015,0.013 },
01223 
01224 //  n S0 K+ pi+ pi- (p S0 K0 pi+ pi-)
01225  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01226     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01227     0.0,  0.0,  0.001,0.005,0.011,0.012,0.011,0.009,0.007,0.006 },
01228 
01229 //  n S0 K+ pi0 pi0 (p S0 K0 pi0 pi0)
01230  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01231     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01232     0.0,  0.0,  0.0,  0.002,0.005,0.006,0.005,0.005,0.004,0.003 },
01233 
01234 //  n S0 K0 pi+ pi0 (p S0 K+ pi- pi0)
01235  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01236     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01237     0.0,  0.0,  0.001,0.005,0.011,0.012,0.011,0.009,0.007,0.006 },
01238 
01239 //  n S+ K0 pi+ pi- (p S- K+ pi+ pi-)
01240  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01241     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01242     0.0,  0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
01243 
01244 //  n S+ K0 pi0 pi0 (p S- K+ pi0 pi0)
01245  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01246     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01247     0.0,  0.0,  0.001,0.003,0.007,0.008,0.007,0.006,0.005,0.004 },
01248 
01249 //  n S+ K+ pi- pi0 (p S- K0 pi+ pi0)
01250  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01251     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01252     0.0,  0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
01253 
01254 //  n S- K+ pi+ pi0 (p S+ K0 pi- pi0)
01255  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01256     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01257     0.0,  0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
01258 
01259 //  n S- K0 pi+ pi+ (p S+ K+ pi- pi-)
01260  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01261     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01262     0.0,  0.001,0.003,0.006,0.014,0.016,0.015,0.012,0.010,0.008 },
01263 
01264 //  p n K+ K- pi0 (p n K0 K0bar pi0)
01265  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01266     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01267     0.0,  0.0,  0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
01268 
01269 //  p n K0 K0bar pi0 (p n K+ K- pi0)
01270  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01271     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01272     0.0,  0.0,  0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
01273 
01274 //  p n K0 K- pi+ (p n K+ K0bar pi-)
01275  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01276     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01277     0.0,  0.0,  0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
01278 
01279 //  p n K+ K0bar pi- (p n K0 K- pi+)
01280  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01281     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01282     0.0,  0.0,  0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
01283 
01284 //  p p K0 K0bar pi- (n n K+ K- pi+)
01285  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01286     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01287     0.0,  0.0,  0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
01288 
01289 //  p p K+ K- pi- (n n K0 K0bar pi+)
01290  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01291     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01292     0.0,  0.0,  0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
01293 
01294 //  p p K0 K- pi0 (n n K+ K0bar pi0) 
01295  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01296     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01297     0.0,  0.0,  0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
01298 
01299 //  n n K+ K- pi+ (p p K0 K0bar pi-)
01300  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01301     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01302     0.0,  0.0,  0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
01303 
01304 //  n n K0 K0bar pi+ (p p K+ K- pi-)
01305  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01306     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01307     0.0,  0.0,  0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
01308 
01309 //  n n K+ K0bar pi0 (p p K0 K- pi0)
01310  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01311     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01312     0.0,  0.0,  0.002,0.004,0.010,0.012,0.011,0.010,0.009,0.007 },
01313 //
01314 // multiplicity 6 (7 channels)
01315 //
01316 //  p n pi+ pi+ pi- pi- (p n pi+ pi+ pi- pi-)
01317  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01318     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01319     0.06, 0.1,  0.18, 0.38, 0.49, 0.46, 0.43, 0.40, 0.38, 0.36 },
01320 
01321 //  p n pi+ pi- pi0 pi0 (p n pi+ pi- pi0 pi0)
01322  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01323     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01324     0.03, 0.05, 0.09, 0.19, 0.25, 0.23, 0.22, 0.2,  0.19, 0.18 },
01325 
01326 //  p n pi0 pi0 pi0 pi0 (p n pi0 pi0 pi0 pi0)
01327  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01328     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01329     0.01, 0.02, 0.05, 0.1,  0.13, 0.12, 0.11, 0.1,  0.1,  0.09 },
01330 
01331 //  p p pi+ pi- pi- pi0 (n n pi+ pi+ pi- pi0)
01332  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01333     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01334     0.06, 0.1,  0.18, 0.38, 0.49, 0.46, 0.43, 0.40, 0.38, 0.36 },
01335 
01336 //  p p pi- pi0 pi0 pi0 (n n pi+ pi0 pi0 pi0)
01337  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01338     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01339     0.03, 0.05, 0.09, 0.19, 0.25, 0.23, 0.22, 0.2,  0.19, 0.18 },
01340 
01341 //  n n pi+ pi+ pi- pi0 (p p pi+ pi- pi- pi0)
01342  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01343     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01344     0.06, 0.1,  0.18, 0.38, 0.49, 0.46, 0.43, 0.40, 0.38, 0.36 },
01345 
01346 //  n n pi+ pi0 pi0 pi0 (p p pi- pi0 pi0 pi0)
01347  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01348     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01349     0.03, 0.05, 0.09, 0.19, 0.25, 0.23, 0.22, 0.2,  0.19, 0.18 },
01350 //
01351 // multiplicity 7 (9 channels)
01352 //
01353 //  p n pi+ pi+ pi- pi- pi0 (p n pi+ pi+ pi- pi- pi0)
01354  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01355     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01356     0.0,  0.0,  0.06, 0.17, 0.5,  0.7,  0.7,  0.69, 0.66, 0.62 },
01357 
01358 //  p n pi+ pi- pi0 pi0 pi0 (p n pi+ pi- pi0 pi0 pi0)
01359  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01360     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01361     0.0,  0.0,  0.03, 0.08, 0.25, 0.35, 0.35, 0.35, 0.33, 0.31 },
01362 
01363 //  p n pi0 pi0 pi0 pi0 pi0 (p n pi0 pi0 pi0 pi0 pi0)
01364  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01365     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01366     0.0,  0.0,  0.015,0.04, 0.12, 0.17, 0.18, 0.17, 0.16, 0.15 },
01367 
01368 //  p p pi+ pi+ pi- pi- pi- (n n pi+ pi+ pi+ pi- pi-)
01369  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01370     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01371     0.0,  0.0,  0.06, 0.19, 0.31, 0.41, 0.44, 0.47, 0.45, 0.45 },
01372 
01373 //  p p pi+ pi- pi- pi0 pi0 (n n pi+ pi+ pi- pi0 pi0)
01374  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01375     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01376     0.0,  0.0,  0.03, 0.1,  0.15, 0.2,  0.22, 0.23, 0.22, 0.22 },
01377 
01378 //  p p pi- pi0 pi0 pi0 pi0 (n n pi+ pi0 pi0 pi0 pi0)
01379  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01380     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01381     0.0,  0.0,  0.02, 0.05, 0.07, 0.1,  0.11, 0.12, 0.11, 0.11 },
01382 
01383 //  n n pi+ pi+ pi+ pi- pi- (p p pi+ pi+ pi- pi- pi-)
01384  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01385     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01386     0.0,  0.0,  0.06, 0.17, 0.5,  0.7,  0.7,  0.69, 0.66, 0.62 },
01387 
01388 //  n n pi+ pi+ pi- pi0 pi0 (p p pi+ pi- pi- pi0 pi0)
01389  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01390     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01391     0.0,  0.0,  0.03, 0.08, 0.25, 0.35, 0.35, 0.34, 0.33, 0.31 },
01392 
01393 //  n n pi+ pi0 pi0 pi0 pi0 (p p pi- pi0 pi0 pi0 pi0)
01394  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01395     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01396     0.0,  0.0,  0.02, 0.05, 0.07, 0.1,  0.11, 0.12, 0.11, 0.11 },
01397 //
01398 // multiplicity 8 (10 channels)
01399 //
01400 //  p n pi+ pi+ pi+ pi- pi- pi- (p n pi+ pi+ pi+ pi- pi- pi-)
01401  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01402     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01403     0.0,  0.003,0.01, 0.024,0.075,0.18, 0.27, 0.30, 0.27, 0.24 },
01404 
01405 //  p n pi+ pi+ pi- pi- pi0 pi0 (p n pi+ pi+ pi- pi- pi0 pi0)
01406  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01407     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01408     0.0,  0.003,0.01, 0.024,0.075,0.18, 0.27, 0.30, 0.27, 0.24 },
01409 
01410 //  p n pi+ pi- pi0 pi0 pi0 pi0 (p n pi+ pi- pi0 pi0 pi0 pi0)
01411  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01412     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01413     0.0,  0.003,0.006,0.015,0.045,0.12, 0.15, 0.18, 0.15, 0.15 },
01414 
01415 //  p n pi0 pi0 pi0 pi0 pi0 pi0 (p n pi0 pi0 pi0 pi0 pi0 pi0)
01416  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01417     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01418     0.0,  0.0,  0.003,0.01, 0.024,0.06, 0.09, 0.12, 0.09, 0.09 },
01419 
01420 //  p p pi+ pi+ pi- pi- pi- pi0 (n n pi+ pi+ pi+ pi- pi- pi0)
01421  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01422     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01423     0.0,  0.003,0.01, 0.024,0.075,0.18, 0.27, 0.30, 0.27, 0.24 },
01424 
01425 //  p p pi+ pi- pi- pi0 pi0 pi0 (n n pi+ pi+ pi- pi0 pi0 pi0)
01426  { 0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01427    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01428    0.0,  0.003,0.006,0.015,0.045,0.12, 0.15, 0.18, 0.15, 0.15 },
01429 
01430 //  p p pi- pi0 pi0 pi0 pi0 pi0 (n n pi+ pi0 pi0 pi0 pi0 pi0)
01431  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01432     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01433     0.0,  0.0,  0.003,0.01, 0.027,0.06, 0.09, 0.12, 0.09, 0.09 },
01434 
01435 //  n n pi+ pi+ pi+ pi- pi- pi0 (p p pi+ pi+ pi- pi- pi- pi0)
01436  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01437     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01438     0.0,  0.003,0.01, 0.024,0.075,0.18, 0.27, 0.30, 0.27, 0.24 },
01439 
01440 //  n n pi+ pi+ pi- pi0 pi0 pi0 (p p pi+ pi- pi- pi0 pi0 pi0)
01441  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01442     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01443     0.0,  0.003,0.006,0.015,0.045,0.12, 0.15, 0.18, 0.15, 0.15 },
01444 
01445 //  n n pi+ pi0 pi0 pi0 pi0 pi0 (p p pi- pi0 pi0 pi0 pi0 pi0)
01446  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01447     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01448     0.0,  0.0,  0.003,0.009,0.027,0.06, 0.09, 0.12, 0.09, 0.09 },
01449 //
01450 // multiplicity 9 (12 channels)
01451 //
01452 //  p n pi+ pi+ pi+ pi- pi- pi- pi0 (p n pi+ pi+ pi+ pi- pi- pi- pi0)
01453  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01454     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01455     0.0,  0.002,0.008,0.025,0.074,0.11, 0.14, 0.15, 0.15, 0.15 },
01456 
01457 //  p n pi+ pi+ pi- pi- pi0 pi0 pi0 (p n pi+ pi+ pi- pi- pi0 pi0 pi0)
01458  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01459     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01460     0.0,  0.001,0.005,0.015,0.045,0.07, 0.084,0.09, 0.09, 0.09 },
01461 
01462 //  p n pi+ pi- pi0 pi0 pi0 pi0 pi0 (p n pi+ pi- pi0 pi0 pi0 pi0 pi0)
01463  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01464     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01465     0.0,  0.001,0.003,0.009,0.027,0.042,0.05, 0.054,0.054,0.054 },
01466 
01467 //  p n pi0 pi0 pi0 pi0 pi0 pi0 pi0 (p n pi0 pi0 pi0 pi0 pi0 pi0 pi0)
01468  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01469     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01470     0.0,  0.0,  0.0,  0.005,0.016,0.025,0.03, 0.032,0.032,0.032 },
01471 
01472 //  p p pi+ pi+ pi+ pi- pi- pi- pi- (n n pi+ pi+ pi+ pi+ pi- pi- pi-)
01473  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01474     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01475     0.0,  0.002,0.006,0.02, 0.06, 0.15, 0.19, 0.22, 0.22, 0.22 },
01476 
01477 //  p p pi+ pi+ pi- pi- pi- pi0 pi0 (n n pi+ pi+ pi+ pi- pi- pi0 pi0)
01478  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01479     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01480     0.0,  0.002,0.006,0.02, 0.06, 0.15, 0.19, 0.22, 0.22, 0.22 },
01481 
01482 //  p p pi+ pi- pi- pi0 pi0 pi0 pi0 (n n pi+ pi+ pi- pi0 pi0 pi0 pi0)
01483  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01484     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01485     0.0,  0.001,0.004,0.012,0.036,0.09, 0.12, 0.13, 0.13, 0.13 },
01486 
01487 //  p p pi- pi0 pi0 pi0 pi0 pi0 pi0 (n n pi+ pi0 pi0 pi0 pi0 pi0 pi0)
01488  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01489     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01490     0.0,  0.001,0.002,0.007,0.022,0.054,0.072,0.078,0.078,0.078 },
01491 
01492 //  n n pi+ pi+ pi+ pi+ pi- pi- pi- (p p pi+ pi+ pi+ pi- pi- pi- pi-)
01493  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01494     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01495     0.0,  0.002,0.006,0.02, 0.06, 0.15, 0.19, 0.22, 0.22, 0.22 },
01496 
01497 //  n n pi+ pi+ pi+ pi- pi- pi0 pi0 (p p pi+ pi+ pi- pi- pi- pi0 pi0)
01498  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01499     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01500     0.0,  0.002,0.006,0.02, 0.06, 0.15, 0.19, 0.22, 0.22, 0.22 },
01501 
01502 //  n n pi+ pi+ pi- pi0 pi0 pi0 pi0 (p p pi+ pi- pi- pi0 pi0 pi0 pi0)
01503  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01504     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01505     0.0,  0.001,0.004,0.012,0.036,0.09, 0.12, 0.13, 0.13, 0.13 },
01506 
01507 //  n n pi+ pi0 pi0 pi0 pi0 pi0 pi0 (p p pi- pi0 pi0 pi0 pi0 pi0 pi0)
01508  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01509     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01510     0.0,  0.001,0.002,0.007,0.022,0.054,0.072,0.078,0.078,0.078 }};

Generated on Mon May 27 17:49:45 2013 for Geant4 by  doxygen 1.4.7