G4RPGPionInelastic.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 "G4RPGPionInelastic.hh"
00030 #include "G4SystemOfUnits.hh"
00031 #include "Randomize.hh"
00032 
00033 G4RPGPionInelastic::G4RPGPionInelastic(const G4String& modelName)
00034  :G4RPGInelastic(modelName)
00035 {
00036   SetMinEnergy( 0.0 );
00037   SetMaxEnergy( 30.*GeV );
00038 
00039   // Initialize t32_dSigma_dMult, t12_dSigma_dMult,
00040   //   pi-nucleon inelastic cross sections for a given multiplicity 
00041   //   for |T_z| = 3/2 and 1/2, respectively 
00042 
00043   G4int i, k, j;
00044   G4int start, stop;
00045 
00046   for (j = 0; j < 8; j++) {
00047     start = pipPindex[j][0];
00048     stop = pipPindex[j][1] + 1;
00049     for (k = 0; k < 30; k++) {
00050       t32_dSigma_dMult[j][k] = 0.0;
00051       for (i = start; i < stop; i++) t32_dSigma_dMult[j][k] += pipPCrossSections[i][k];
00052     }
00053 
00054     start = pimPindex[j][0];
00055     stop = pimPindex[j][1] + 1;
00056     for (k = 0; k < 30; k++) {
00057       t12_dSigma_dMult[j][k] = 0.0;
00058       for (i = start; i < stop; i++) t12_dSigma_dMult[j][k] += pimPCrossSections[i][k];
00059     }
00060   }
00061 
00062   // Initialize total cross section array
00063   for (k = 0; k < 30; k++) {
00064     pipPtot[k] = 0.0;
00065     pimPtot[k] = 0.0;
00066     for (j = 0; j < 8; j++) {
00067       pipPtot[k] += t32_dSigma_dMult[j][k];
00068       pimPtot[k] += t12_dSigma_dMult[j][k];
00069     }
00070   }
00071 
00072   //  printCrossSections();
00073 }
00074 
00075 
00076 /*
00077 void G4RPGPionInelastic::printCrossSections() const
00078 {
00079   G4cout << " pi+ p 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 << pipPtot[t] << "  " ;
00084     G4cout << G4endl;
00085   }
00086   G4cout << G4endl;
00087 
00088   G4cout << " pi- p 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 << pimPtot[t] << "  " ;
00093     G4cout << G4endl;
00094   }
00095 }
00096 */
00097 
00098 
00099 G4int G4RPGPionInelastic::GetMultiplicityT12(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 = t12_dSigma_dMult[j][k]
00110          + fraction*(t12_dSigma_dMult[j][k+1] - t12_dSigma_dMult[j][k]);
00111       sigma.push_back(multint);
00112   }
00113 
00114   return sampleFlat(sigma) + 2;
00115 }
00116 
00117 
00118 G4int G4RPGPionInelastic::GetMultiplicityT32(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 = t32_dSigma_dMult[j][k]
00129          + fraction*(t32_dSigma_dMult[j][k+1] - t32_dSigma_dMult[j][k]);
00130       sigma.push_back(multint);
00131   }
00132 
00133   return sampleFlat(sigma) + 2;
00134 }
00135 
00136 
00137 std::vector<G4int> 
00138 G4RPGPionInelastic::GetFSPartTypesForT32(G4int mult, G4double KE, G4int tzindex) 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 = pipPindex[mult-2][0];
00149   G4int stop = pipPindex[mult-2][1];
00150 
00151   for(i = start; i < stop; i++) {
00152       sigint = pipPCrossSections[i][k]
00153           + fraction*(pipPCrossSections[i][k+1] - pipPCrossSections[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(T32_2bfs[tzindex][channel][i]);
00163   } else if (mult == 3) {
00164     for(i = 0; i < mult; i++) kinds.push_back(T32_3bfs[tzindex][channel][i]);
00165   } else if (mult == 4) {
00166     for(i = 0; i < mult; i++) kinds.push_back(T32_4bfs[tzindex][channel][i]);
00167   } else if (mult == 5) {
00168     for(i = 0; i < mult; i++) kinds.push_back(T32_5bfs[tzindex][channel][i]);
00169   } else if (mult == 6) {
00170     for(i = 0; i < mult; i++) kinds.push_back(T32_6bfs[tzindex][channel][i]);
00171   } else if (mult == 7) {
00172     for(i = 0; i < mult; i++) kinds.push_back(T32_7bfs[tzindex][channel][i]);
00173   } else if (mult == 8) {
00174     for(i = 0; i < mult; i++) kinds.push_back(T32_8bfs[tzindex][channel][i]);
00175   } else if (mult == 9) {
00176     for(i = 0; i < mult; i++) kinds.push_back(T32_9bfs[tzindex][channel][i]);
00177   } else {
00178     G4cout << " Illegal multiplicity " << G4endl;
00179   }
00180 
00181   return kinds;
00182 }
00183 
00184 
00185 std::vector<G4int> 
00186 G4RPGPionInelastic::GetFSPartTypesForT12(G4int mult, G4double KE, G4int tzindex) const
00187 {
00188   G4int i;
00189   G4double sigint(0.);
00190   std::vector<G4double> sigma;
00191 
00192   std::pair<G4int, G4double> epair = interpolateEnergy(KE);
00193   G4int k = epair.first;
00194   G4double fraction = epair.second;
00195 
00196   G4int start = pimPindex[mult-2][0];
00197   G4int stop = pimPindex[mult-2][1];
00198 
00199   for(i = start; i < stop; i++) {
00200       sigint = pimPCrossSections[i][k]
00201           + fraction*(pimPCrossSections[i][k+1] - pimPCrossSections[i][k]);
00202       sigma.push_back(sigint);
00203   }
00204 
00205   G4int channel = sampleFlat(sigma);
00206 
00207   std::vector<G4int> kinds;
00208 
00209   if (mult == 2) {
00210     for(i = 0; i < mult; i++) kinds.push_back(T12_2bfs[tzindex][channel][i]);
00211   } else if (mult == 3) {
00212     for(i = 0; i < mult; i++) kinds.push_back(T12_3bfs[tzindex][channel][i]);
00213   } else if (mult == 4) {
00214     for(i = 0; i < mult; i++) kinds.push_back(T12_4bfs[tzindex][channel][i]);
00215   } else if (mult == 5) {
00216     for(i = 0; i < mult; i++) kinds.push_back(T12_5bfs[tzindex][channel][i]);
00217   } else if (mult == 6) {
00218     for(i = 0; i < mult; i++) kinds.push_back(T12_6bfs[tzindex][channel][i]);
00219   } else if (mult == 7) {
00220     for(i = 0; i < mult; i++) kinds.push_back(T12_7bfs[tzindex][channel][i]);
00221   } else if (mult == 8) {
00222     for(i = 0; i < mult; i++) kinds.push_back(T12_8bfs[tzindex][channel][i]);
00223   } else if (mult == 9) {
00224     for(i = 0; i < mult; i++) kinds.push_back(T12_9bfs[tzindex][channel][i]);
00225   } else {
00226     G4cout << " Illegal multiplicity " << G4endl;
00227   }
00228 
00229   return kinds;
00230 }
00231 
00232 
00234 //                                                                           //
00235 //   pi+ p and pi- n (|Tz| = 3/2) cross sections                             //
00236 //   and final state particle types                                          //
00237 //                                                                           //
00239 
00240 // Total pi+ cross section as a function of kinetic energy
00241 G4double G4RPGPionInelastic::pipPtot[30];
00242 
00243 // pi+ multiplicities as a function of kinetic energy
00244 G4double G4RPGPionInelastic::t32_dSigma_dMult[8][30];
00245 
00246 const G4int G4RPGPionInelastic::pipPindex[8][2] =
00247  {{0, 1}, {2, 8}, {9,23}, {24,47}, {48,52}, {53,58}, {59,65}, {66,73}};  
00248 
00249 // Outgoing particle types of a given multiplicity
00250 // T32_nbfs = final state types for pi+ p and pi- n
00251 
00252 const G4int G4RPGPionInelastic::T32_2bfs[2][2][2] =
00253   {{{pro,pip}, {sp,kp}},
00254 
00255    {{neu,pim}, {sm,k0}}};
00256 
00257 const G4int G4RPGPionInelastic::T32_3bfs[2][7][3] =
00258   {{{pro,pip,pi0}, {neu,pip,pip}, {sp,kp,pi0}, {sp,k0,pip}, 
00259     {s0,kp,pip},   {lam,kp,pip},  {pro,kp,k0b}},
00260 
00261    {{neu,pim,pi0}, {pro,pim,pim}, {sm,k0,pi0}, {sm,kp,pim},
00262     {s0,k0,pim},   {lam,k0,pim},  {neu,k0,km}}};
00263 
00264 const G4int G4RPGPionInelastic::T32_4bfs[2][15][4] =
00265   {{{pro,pip,pip,pim},{pro,pip,pi0,pi0},{neu,pip,pip,pi0},
00266     {sp,kp,pip,pim},  {sp,kp,pi0,pi0},  {sp,k0,pip,pi0},
00267     {s0,k0,pip,pip},  {s0,kp,pip,pi0},  {lam,kp,pip,pi0},
00268     {lam,k0,pip,pip}, {sm,kp,pip,pip},  {pro,pip,kp,km},
00269     {pro,pip,k0,k0b}, {pro,pi0,kp,k0b}, {neu,pip,kp,k0b}},
00270 
00271    {{neu,pip,pim,pim},{neu,pim,pi0,pi0},{pro,pim,pim,pi0},
00272     {sm,k0,pip,pim},  {sm,k0,pi0,pi0},  {sm,kp,pim,pi0},
00273     {s0,kp,pim,pim},  {s0,k0,pim,pi0},  {lam,k0,pim,pi0},
00274     {lam,kp,pim,pim}, {sp,k0,pim,pim},  {neu,pim,k0,k0b},
00275     {neu,pim,kp,km},  {neu,pi0,k0,km},  {pro,pim,k0,km}}};
00276 
00277 const G4int G4RPGPionInelastic::T32_5bfs[2][24][5] =
00278   {{{pro,pip,pip,pim,pi0}, {pro,pip,pi0,pi0,pi0}, {neu,pip,pip,pip,pim},
00279     {neu,pip,pip,pi0,pi0}, {sp,kp,pip,pim,pi0},   {sp,kp,pi0,pi0,pi0},
00280     {sp,k0,pip,pip,pim},   {sp,k0,pip,pi0,pi0},   {lam,k0,pip,pip,pi0},
00281     {lam,kp,pip,pip,pim},  {lam,kp,pip,pi0,pi0},  {s0,kp,pip,pip,pim},
00282     {s0,kp,pip,pi0,pi0},   {s0,k0,pip,pip,pi0},   {sm,kp,pip,pip,pi0},
00283     {sm,k0,pip,pip,pip},   {pro,pip,pim,kp,k0b},  {pro,pip,pip,k0,km},
00284     {pro,pip,pi0,kp,km},   {pro,pip,pi0,k0,k0b},  {pro,pi0,pi0,kp,k0b},
00285     {neu,pip,pip,kp,km},   {neu,pip,pip,k0,k0b},  {neu,pip,pi0,kp,k0b}},
00286 
00287    {{neu,pip,pim,pim,pi0}, {neu,pim,pi0,pi0,pi0}, {pro,pip,pim,pim,pim},
00288     {pro,pim,pim,pi0,pi0}, {sm,k0,pip,pim,pi0},   {sm,k0,pi0,pi0,pi0},
00289     {sm,kp,pip,pim,pim},   {sm,kp,pim,pi0,pi0},   {lam,kp,pim,pim,pi0},
00290     {lam,k0,pip,pim,pim},  {lam,k0,pim,pi0,pi0},  {s0,k0,pip,pim,pim},
00291     {s0,k0,pim,pi0,pi0},   {s0,kp,pim,pim,pi0},   {sp,k0,pim,pim,pi0},
00292     {sp,kp,pim,pim,pim},   {neu,pip,pim,k0,km},   {neu,pim,pim,kp,k0b},
00293     {neu,pim,pi0,k0,k0b},  {neu,pim,pi0,kp,km},   {neu,pi0,pi0,k0,km},
00294     {pro,pim,pim,k0,k0b},  {pro,pim,pim,kp,km},   {pro,pim,pi0,k0,km}}};
00295 
00296 const G4int G4RPGPionInelastic::T32_6bfs[2][5][6] =
00297 {{{pro,pip,pip,pip,pim,pim}, {pro,pip,pip,pim,pi0,pi0},
00298   {pro,pip,pi0,pi0,pi0,pi0}, {neu,pip,pip,pi0,pi0,pi0},
00299   {neu,pip,pip,pip,pim,pi0}},
00300 
00301  {{neu,pip,pip,pim,pim,pim}, {neu,pip,pim,pim,pi0,pi0},
00302   {neu,pim,pi0,pi0,pi0,pi0}, {pro,pim,pim,pi0,pi0,pi0},
00303   {pro,pip,pim,pim,pim,pi0}}};
00304 
00305 const G4int G4RPGPionInelastic::T32_7bfs[2][6][7] =
00306 {{{pro,pip,pip,pip,pim,pim,pi0}, {pro,pip,pip,pim,pi0,pi0,pi0},
00307   {pro,pip,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pip,pip,pip,pim,pim},
00308   {neu,pip,pip,pip,pim,pi0,pi0}, {neu,pip,pip,pi0,pi0,pi0,pi0}},
00309 
00310  {{neu,pip,pip,pim,pim,pim,pi0}, {neu,pip,pim,pim,pi0,pi0,pi0},
00311   {neu,pim,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pip,pim,pim,pim,pim},
00312   {pro,pip,pim,pim,pim,pi0,pi0}, {pro,pim,pim,pi0,pi0,pi0,pi0}}};
00313 
00314 const G4int G4RPGPionInelastic::T32_8bfs[2][7][8] =
00315 {{{pro,pip,pip,pip,pip,pim,pim,pim}, {pro,pip,pip,pip,pim,pim,pi0,pi0},
00316   {pro,pip,pip,pim,pi0,pi0,pi0,pi0}, {pro,pip,pi0,pi0,pi0,pi0,pi0,pi0},
00317   {neu,pip,pip,pip,pip,pim,pim,pi0}, {neu,pip,pip,pip,pim,pi0,pi0,pi0},
00318   {neu,pip,pip,pi0,pi0,pi0,pi0,pi0}},
00319 
00320  {{neu,pip,pip,pip,pim,pim,pim,pim}, {neu,pip,pip,pim,pim,pim,pi0,pi0},
00321   {neu,pip,pim,pim,pi0,pi0,pi0,pi0}, {neu,pim,pi0,pi0,pi0,pi0,pi0,pi0},
00322   {pro,pip,pip,pim,pim,pim,pim,pi0}, {pro,pip,pim,pim,pim,pi0,pi0,pi0},
00323   {pro,pim,pim,pi0,pi0,pi0,pi0,pi0}}};
00324 
00325 const G4int G4RPGPionInelastic::T32_9bfs[2][8][9] =
00326 {{{pro,pip,pip,pip,pip,pim,pim,pim,pi0}, {pro,pip,pip,pip,pim,pim,pi0,pi0,pi0},
00327   {pro,pip,pip,pim,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
00328   {neu,pip,pip,pip,pip,pip,pim,pim,pim}, {neu,pip,pip,pip,pip,pim,pim,pi0,pi0},
00329   {neu,pip,pip,pip,pim,pi0,pi0,pi0,pi0}, {neu,pip,pip,pi0,pi0,pi0,pi0,pi0,pi0}},
00330 
00331  {{neu,pip,pip,pip,pim,pim,pim,pim,pi0}, {neu,pip,pip,pim,pim,pim,pi0,pi0,pi0},
00332   {neu,pip,pim,pim,pi0,pi0,pi0,pi0,pi0}, {neu,pim,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
00333   {pro,pip,pip,pip,pim,pim,pim,pim,pim}, {pro,pip,pip,pim,pim,pim,pim,pi0,pi0},
00334   {pro,pip,pim,pim,pim,pi0,pi0,pi0,pi0}, {pro,pim,pim,pi0,pi0,pi0,pi0,pi0,pi0}}};
00335 
00336 //
00337 // Cross sections (in mb) for pi+ p -> 2-9 body final states
00338 //
00339 // first index:    0-1: channels for mult = 2
00340 //                 2-8: channels for mult = 3
00341 //                9-23: channels for mult = 4
00342 //               24-47: channels for mult = 5
00343 //               48-52: channels for mult = 6
00344 //               53-58: channels for mult = 7
00345 //               59-65: channels for mult = 8
00346 //               66-73: channels for mult = 9
00347 //
00348 // second index: kinetic energy
00349 //
00350 
00351 const G4float G4RPGPionInelastic::pipPCrossSections[74][30] = {
00352 //
00353 // multiplicity 2 (2 channels)
00354 //
00355 //  p pi+ (n pi-)
00356  { 0.00, 1.20, 2.50, 3.80, 5.00, 7.00, 9.00, 15.0, 30.0, 64.0,
00357   130.0,190.0,130.0, 55.7, 27.2, 14.0, 8.50, 13.0, 18.0, 11.0,
00358    8.50, 7.00, 6.20, 5.60, 5.00, 4.50, 4.20, 4.00, 3.80, 3.60 },
00359 
00360 //  S+ K+ (S- K0)
00361  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00362     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.16, 0.60, 0.32,
00363    0.19,  0.1, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01 },
00364 //
00365 // multiplicity 3 (7 channels)
00366 //
00367 //  p pi+ pi0 (n pi- pi0)
00368  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00369     0.0,  0.0,  0.0,  0.2,  0.6,  2.4,  8.8, 10.0, 12.0,  6.2,
00370    4.00, 2.40, 1.69, 1.10, 0.73, 0.49, 0.41, 0.31, 0.24, 0.15 },
00371 
00372 //  n pi+ pi+ (p pi- pi-)
00373  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00374     0.0,  0.0,  0.0, 0.10, 0.20, 0.60, 1.50, 2.30, 3.60, 3.00,
00375    2.30, 1.70, 1.30, 0.95, 0.69, 0.46, 0.38, 0.27, 0.20, 0.15 },
00376 
00377 //  S+ K+ pi0 (S- K0 pi0)
00378  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00379     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.005,0.12,
00380    0.14, 0.09, 0.07, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01 },
00381 
00382 //  S+ K0 pi+
00383  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00384     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.005,0.12,
00385    0.14, 0.09, 0.07, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01 },
00386 
00387 //  S0 K+ pi+
00388  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00389     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.005,0.12,
00390    0.14, 0.09, 0.07, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01 },
00391 
00392 //  L K+ pi+
00393  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00394     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.005,0.12,
00395    0.14, 0.09, 0.07, 0.06, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01 },
00396 
00397 //  p K+ K0bar
00398  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00399     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.04,
00400    0.06, 0.05, 0.04, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01 },
00401 //
00402 // multiplicity 4 (15 channels)
00403 //
00404 //  p pi+ pi+ pi-
00405  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00406     0.0,  0.0,  0.0,  0.0,  0.0, 0.06, 0.20, 0.78, 2.20, 3.20,
00407    3.50, 3.10, 2.70, 2.30, 2.00, 1.50, 1.40, 1.20, 1.00, 0.90 },
00408 
00409 //  p pi+ pi0 pi0
00410  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00411     0.0,  0.0,  0.0,  0.0,  0.0, 0.04, 0.13, 0.52, 1.50, 2.20,
00412    2.40, 2.00, 1.80, 1.60, 1.30, 1.10, 1.00, 0.80, 0.70, 0.60 },
00413 
00414 //  n pi+ pi+ pi0
00415  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00416     0.0,  0.0,  0.0,  0.0,  0.0, 0.04, 0.13, 0.52, 1.50, 2.20,
00417    2.40, 2.00, 1.80, 1.60, 1.30, 1.10, 1.00, 0.80, 0.70, 0.60 },
00418 
00419 //  S+ K+ pi+ pi-
00420  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00421     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00422    0.003,0.04, 0.12, 0.08, 0.05, 0.03, 0.02, 0.01,0.007,0.004 },
00423 
00424 //  S+ K+ pi0 pi0
00425  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00426     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00427    0.003,0.04, 0.12, 0.08, 0.05, 0.03, 0.02, 0.01,0.007,0.004 },
00428 
00429 //  S+ K0 pi+ pi0
00430  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00431     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.002,
00432    0.015,0.06, 0.06, 0.05, 0.04,0.032,0.028, 0.02,0.017,0.014 },
00433 
00434 //  S0 K0 pi+ pi+
00435  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00436     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00437    0.002,0.01, 0.02, 0.02,0.015,0.012,0.011,0.009,0.008,0.007 },
00438 
00439 //  S0 K+ pi+ pi0
00440  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00441     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00442    0.002,0.01, 0.02, 0.02,0.015,0.012,0.011,0.009,0.008,0.007 },
00443 
00444 //  L K+ pi+ pi0
00445  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00446     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00447    0.003,0.02, 0.04, 0.04, 0.03, 0.03, 0.02, 0.02,0.016,0.014 },
00448 
00449 //  L K0 pi+ pi+
00450  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00451     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00452    0.003,0.02, 0.04, 0.04, 0.03, 0.03, 0.02, 0.02,0.016,0.014 },
00453 
00454 //  S- K+ pi+ pi+
00455  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00456     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.002,
00457    0.02, 0.16, 0.18, 0.13, 0.09, 0.06,0.045, 0.03,0.025, 0.02 },
00458 
00459 //  p pi+ K+ K-
00460  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00461     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00462    0.02, 0.13, 0.11, 0.10, 0.08, 0.07, 0.06, 0.05,0.045, 0.04 },
00463 
00464 //  p pi+ K0 K0bar
00465  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00466     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00467   0.005, 0.05, 0.15, 0.14, 0.11, 0.09, 0.07, 0.05,0.045, 0.04 },
00468 
00469 //  p pi0 K+ K0bar
00470  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00471     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00472    0.01, 0.05,0.075, 0.08,0.075, 0.07, 0.06, 0.05,0.045, 0.04 },
00473 
00474 //  n pi+ K+ K0bar
00475  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00476     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00477    0.01, 0.05, 0.07,0.065,0.055,0.048, 0.04, 0.03,0.027,0.023 },
00478 //
00479 // multiplicity 5 (24 channels)
00480 //
00481 //  p pi+ pi+ pi- pi0
00482  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00483     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.001, 0.05, 0.30, 2.00,
00484    3.20, 3.70, 3.00, 2.50, 2.10, 1.60, 1.40, 1.10, 0.89, 0.70 },
00485 
00486 //  p pi+ pi0 pi0 pi0
00487  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00488     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.001, 0.03, 0.21, 1.40,
00489    2.20, 2.60, 2.10, 1.80, 1.50, 1.10, 1.00, 0.80, 0.62, 0.50 },
00490 
00491 //  n pi+ pi+ pi+ pi-
00492  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00493     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.007, 0.02, 0.05, 0.19,
00494    0.35, 0.65, 0.90, 0.87, 0.71, 0.55, 0.42, 0.31, 0.24, 0.18 },
00495 
00496 //  n pi+ pi+ pi0 pi0
00497  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00498     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.007, 0.02, 0.05, 0.19,
00499    0.35, 0.65, 0.90, 0.87, 0.71, 0.55, 0.42, 0.31, 0.24, 0.18 },
00500 
00501 //  S+ K+ pi+ pi- pi0
00502  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00503     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00504   0.003,0.013, 0.05, 0.09, 0.08, 0.06, 0.05,  0.04,0.036,0.03 },
00505 
00506 //  S+ K+ pi0 pi0 pi0
00507  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00508     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00509   0.001,0.006, 0.02, 0.04, 0.04, 0.03,0.025, 0.02, 0.02, 0.01 },
00510 
00511 //  S+ K0 pi+ pi+ pi-
00512  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00513     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00514   0.001,0.005,0.015,0.036,0.034,0.029,0.024, 0.02,0.017,0.014 },
00515 
00516 //  S+ K0 pi+ pi0 pi0
00517  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00518     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00519   0.001,0.005,0.015,0.036,0.034,0.029,0.024, 0.02,0.017,0.014 },
00520 
00521 //  L K0 pi+ pi+ pi0
00522  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00523     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00524   0.002,0.013, 0.04,0.052,0.059,0.053, 0.05,0.043,0.037, 0.03 },
00525 
00526 //  L K+ pi+ pi+ pi-
00527  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00528     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00529   0.001,0.005,0.018, 0.04, 0.05,0.041,0.038,0.032,0.028,0.024 },
00530 
00531 //  L K+ pi+ pi0 pi0
00532  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00533     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00534   0.001,0.005,0.018, 0.04, 0.05,0.041,0.038,0.032,0.028,0.024 },
00535 
00536 //  S0 K+ pi+ pi+ pi-
00537  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00538     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00539   0.001,0.002,0.005, 0.01,0.016,0.014,0.012, 0.01,0.009,0.008 },
00540 
00541 //  S0 K+ pi+ pi0 pi0
00542  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00543     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00544   0.001,0.002,0.005, 0.01,0.016,0.014,0.012, 0.01,0.009,0.008 },
00545 
00546 //  S0 K0 pi+ pi+ pi0
00547  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00548     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00549   0.001,0.002,0.005, 0.01,0.016,0.014,0.012, 0.01,0.009,0.008 },
00550 
00551 //  S- K+ pi+ pi+ pi0
00552  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00553     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00554   0.002,0.005,0.015,0.025, 0.02,0.017,0.015,0.013,0.011,0.009 },
00555 
00556 //  S- K0 pi+ pi+ pi+
00557  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00558     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00559   0.002,0.005,0.015,0.025, 0.02,0.017,0.015,0.013,0.011,0.009 },
00560 
00561 //  p pi+ pi- K+ K0bar
00562  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00563     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00564   0.001,0.005, 0.02,0.065, 0.08, 0.07, 0.06,0.054,0.048, 0.04 },
00565 
00566 //  p pi+ pi+ K0 K-
00567  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00568     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00569   0.001,0.005, 0.02,0.045, 0.05,0.047, 0.04,0.033, 0.03,0.026 },
00570 
00571 //  p pi+ pi0 K+ K-
00572  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00573     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00574   0.001, 0.02, 0.06, 0.06,0.054,0.048,0.042,0.038,0.035, 0.03 },
00575 
00576 //  p pi+ pi0 K0 K0bar
00577  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00578     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00579   0.001, 0.01, 0.04, 0.05, 0.05, 0.04, 0.04,0.035,0.032, 0.03 },
00580 
00581 //  p pi0 pi0 K+ K0bar
00582  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00583     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00584   0.001,0.005, 0.01, 0.03, 0.04,0.035, 0.03,0.025,0.025, 0.02 },
00585 
00586 //  n pi+ pi+ K+ K-
00587  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00588     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00589     0.0, 0.02, 0.06, 0.05,0.042,0.038,0.035, 0.03,0.027,0.022 },
00590 
00591 //  n pi+ pi+ K0 K0bar
00592  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00593     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00594     0.0, 0.01, 0.04, 0.03,0.028,0.024,0.023, 0.02,0.018,0.014 },
00595 
00596 //  n pi+ pi0 K+ K0bar
00597  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00598     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00599     0.0, 0.01, 0.04, 0.03,0.028,0.024,0.023, 0.02,0.018,0.014 },
00600 //
00601 // multiplicity 6 (5 channels)
00602 //
00603 //  p pi+ pi+ pi+ pi- pi-
00604  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00605     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.001, 0.02,
00606    0.08, 0.20, 0.31, 0.40, 0.42, 0.42, 0.40, 0.32, 0.29, 0.23 },
00607 
00608 //  p pi+ pi+ pi- pi0 pi0
00609  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00610     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.001, 0.02,
00611    0.08, 0.20, 0.31, 0.40, 0.42, 0.42, 0.40, 0.32, 0.29, 0.23 },
00612 
00613 //  p pi+ pi0 pi0 pi0 pi0
00614  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00615     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.001, 0.02,
00616    0.08, 0.20, 0.31, 0.40, 0.42, 0.42, 0.40, 0.32, 0.29, 0.23 },
00617 
00618 //  n pi+ pi+ pi0 pi0 pi0
00619  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00620     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.001, 0.02,
00621    0.08, 0.20, 0.31, 0.40, 0.42, 0.42, 0.40, 0.32, 0.29, 0.23 },
00622 
00623 //  n pi+ pi+ pi+ pi- pi0
00624  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00625     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.001, 0.02,
00626    0.08, 0.20, 0.31, 0.40, 0.42, 0.42, 0.40, 0.32, 0.29, 0.23 },
00627 //
00628 // multiplicity 7 (6 channels)
00629 //
00630 //  p pi+ pi+ pi+ pi- pi- pi0
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00633   0.018, 0.10, 0.36, 0.96, 0.96, 0.96, 0.90, 0.84, 0.78, 0.72 },
00634 
00635 //  p pi+ pi+ pi- pi0 pi0 pi0
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00638   0.018, 0.10, 0.36, 0.96, 0.96, 0.96, 0.90, 0.84, 0.78, 0.72 },
00639 
00640 //  p pi+ pi0 pi0 pi0 pi0 pi0
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00643   0.018, 0.10, 0.36, 0.96, 0.96, 0.96, 0.90, 0.84, 0.78, 0.72 },
00644 
00645 //  n pi+ pi+ pi+ pi+ pi- pi-
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00648    0.01, 0.04, 0.06, 0.12, 0.24, 0.30, 0.30, 0.26, 0.24, 0.22 },
00649 
00650 //  n pi+ pi+ pi+ pi- pi0 pi0
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00653    0.01, 0.04, 0.06, 0.12, 0.24, 0.30, 0.30, 0.26, 0.24, 0.22 },
00654 
00655 //  n pi+ pi+ pi0 pi0 pi0 pi0
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00658    0.01, 0.04, 0.06, 0.12, 0.24, 0.30, 0.30, 0.26, 0.24, 0.22 },
00659 //
00660 // multiplicity 8 (7 channels)
00661 //
00662 //  p pi+ pi+ pi+ pi+ pi- pi- pi-
00663  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00664     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00665   0.015, 0.03,0.045,0.075, 0.12, 0.16, 0.16, 0.15, 0.15, 0.14 },
00666 
00667 //  p pi+ pi+ pi+ pi- pi- pi0 pi0
00668  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00669     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00670   0.015, 0.03,0.045,0.075, 0.12, 0.16, 0.16, 0.15, 0.15, 0.14 },
00671 
00672 //  p pi+ pi+ pi- pi0 pi0 pi0 pi0
00673  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00674     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00675   0.015, 0.03,0.045,0.075, 0.12, 0.16, 0.16, 0.15, 0.15, 0.14 },
00676 
00677 //  p pi+ pi0 pi0 pi0 pi0 pi0 pi0
00678  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00679     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00680   0.015, 0.03,0.045,0.075, 0.12, 0.16, 0.16, 0.15, 0.15, 0.14 },
00681 
00682 //  n pi+ pi+ pi+ pi+ pi- pi- pi0
00683  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00684     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00685   0.036,0.072, 0.11, 0.18, 0.28, 0.40, 0.40, 0.36, 0.36, 0.32 },
00686 
00687 //  n pi+ pi+ pi+ pi- pi0 pi0 pi0
00688  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00689     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00690   0.036,0.072, 0.11, 0.18, 0.28, 0.40, 0.40, 0.36, 0.36, 0.32 },
00691 
00692 //  n pi+ pi+ pi0 pi0 pi0 pi0 pi0
00693  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00694     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00695   0.036,0.072, 0.11, 0.18, 0.28, 0.40, 0.40, 0.36, 0.36, 0.32 },
00696 //
00697 // multiplicity 9 (8 channels)
00698 //
00699 //  p pi+ pi+ pi+ pi+ pi- pi- pi- pi0
00700  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00701     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00702   0.012,0.024,0.036, 0.06, 0.11, 0.18, 0.26, 0.36, 0.36, 0.36 },
00703 
00704 //  p pi+ pi+ pi+ pi- pi- pi0 pi0 pi0
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00707   0.012,0.024,0.036, 0.06, 0.11, 0.18, 0.26, 0.36, 0.36, 0.36 },
00708 
00709 //  p pi+ pi+ pi- pi0 pi0 pi0 pi0 pi0
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00712   0.012,0.024,0.036, 0.06, 0.11, 0.18, 0.26, 0.36, 0.36, 0.36 },
00713 
00714 //  p pi+ pi0 pi0 pi0 pi0 pi0 pi0 pi0
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00717   0.012,0.024,0.036, 0.06, 0.11, 0.18, 0.26, 0.36, 0.36, 0.36 },
00718 
00719 //  n pi+ pi+ pi+ pi+ pi+ pi- pi- pi-
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00722     0.0,  0.0,0.015, 0.03,0.045,0.075, 0.10, 0.15, 0.15, 0.15 },
00723 
00724 //  n pi+ pi+ pi+ pi+ pi- pi- pi0 pi0
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00727     0.0,  0.0,0.015, 0.03,0.045,0.075, 0.10, 0.15, 0.15, 0.15 },
00728 
00729 //  n pi+ pi+ pi+ pi- pi0 pi0 pi0 pi0
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00732     0.0,  0.0,0.015, 0.03,0.045,0.075, 0.10, 0.15, 0.15, 0.15 },
00733 
00734 //  n pi+ pi+ pi0 pi0 pi0 pi0 pi0 pi0
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.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00737     0.0,  0.0,0.015, 0.03,0.045,0.075, 0.10, 0.15, 0.15, 0.15 } };
00738 
00740 //                                                                           //
00741 //   pi- p and pi+ n (|Tz| = 1/2) cross sections                             //
00742 //   and final state particle types                                          //
00743 //                                                                           //
00745 
00746 // Total pi- cross section as a function of kinetic energy
00747 G4double G4RPGPionInelastic::pimPtot[30];
00748 
00749 // pi- multiplicities as a function of kinetic energy
00750 G4double G4RPGPionInelastic::t12_dSigma_dMult[8][30];
00751 
00752 const G4int G4RPGPionInelastic::pimPindex[8][2] =
00753  {{0, 4}, {5,17}, {18,39}, {40,70}, {71,76}, {77,83}, {84,91}, {92,100}};  
00754 
00755 // Outgoing particle types of a given multiplicity
00756 // T12_nbfs = final state types for pi- p and pi+ n
00757 
00758 const G4int G4RPGPionInelastic::T12_2bfs[2][5][2] =
00759   {{{pro,pim}, {neu,pi0}, {lam,k0}, {s0,k0}, {sm,kp}},
00760 
00761    {{neu,pip}, {pro,pi0}, {lam,kp}, {s0,kp}, {sp,k0}}};
00762 
00763 const G4int G4RPGPionInelastic::T12_3bfs[2][13][3] =
00764   {{{pro,pim,pi0}, {neu,pip,pim}, {neu,pi0,pi0}, {lam,k0,pi0}, 
00765     {lam,kp,pim},  {sm,k0,pip},   {sm,kp,pi0},   {sp,k0,pim},
00766     {s0,kp,pim},   {s0,k0,pi0},   {pro,k0,km},   {neu,kp,km},
00767     {neu,k0,k0b}},
00768 
00769    {{neu,pip,pi0}, {pro,pip,pim}, {pro,pi0,pi0}, {lam,kp,pi0},
00770     {lam,k0,pip},  {sp,kp,pim},   {sp,k0,pi0},   {sm,kp,pip},
00771     {s0,k0,pip},   {s0,kp,pi0},   {neu,kp,k0b},  {pro,k0,k0b},
00772     {pro,kp,km}}};
00773 
00774 const G4int G4RPGPionInelastic::T12_4bfs[2][22][4] =
00775   {{{pro,pip,pim,pim}, {pro,pim,pi0,pi0}, {neu,pip,pim,pi0},
00776     {neu,pi0,pi0,pi0}, {lam,k0,pip,pim},  {lam,k0,pi0,pi0},
00777     {lam,kp,pim,pi0},  {s0,k0,pip,pim},   {s0,k0,pi0,pi0},
00778     {s0,kp,pim,pi0},   {sp,kp,pim,pim},   {sp,k0,pim,pi0},
00779     {sm,kp,pip,pim},   {sm,kp,pi0,pi0},   {sm,k0,pip,pi0},
00780     {pro,pim,kp,km},   {pro,pim,k0,k0b},  {pro,pi0,k0,km},
00781     {neu,pip,k0,km},   {neu,pi0,k0,k0b},  {neu,pi0,kp,km},
00782     {neu,pim,kp,k0b}},
00783 
00784    {{neu,pip,pip,pim},  {neu,pip,pi0,pi0}, {pro,pip,pim,pi0},
00785     {pro,pi0,pi0,pi0},  {lam,kp,pip,pim},  {lam,kp,pi0,pi0},
00786     {lam,k0,pip,pi0},   {s0,kp,pip,pim},   {s0,kp,pi0,pi0},
00787     {s0,k0,pip,pi0},    {sm,k0,pip,pip},   {sm,kp,pip,pi0},
00788     {sp,k0,pip,pim},    {sp,k0,pi0,pi0},   {sp,kp,pim,pi0},
00789     {neu,pip,k0,k0b},   {neu,pip,kp,km},   {neu,pi0,kp,k0b},
00790     {pro,pim,kp,k0b},   {pro,pi0,kp,km},   {pro,pi0,k0,k0b},
00791     {pro,pip,k0,km}}};
00792 
00793 const G4int G4RPGPionInelastic::T12_5bfs[2][31][5] =
00794   {{{pro,pip,pim,pim,pi0}, {pro,pim,pi0,pi0,pi0}, {neu,pip,pip,pim,pim},
00795     {neu,pip,pim,pi0,pi0}, {neu,pi0,pi0,pi0,pi0}, {lam,k0,pip,pim,pi0},
00796     {lam,kp,pim,pi0,pi0},  {lam,kp,pip,pim,pim},  {lam,k0,pi0,pi0,pi0},
00797     {s0,kp,pip,pim,pim},   {s0,kp,pim,pi0,pi0},   {s0,k0,pip,pim,pi0},
00798     {s0,k0,pi0,pi0,pi0},   {sp,k0,pip,pim,pim},   {sp,k0,pim,pi0,pi0},
00799     {sp,kp,pim,pim,pi0},   {sm,k0,pip,pip,pim},   {sm,k0,pip,pi0,pi0},
00800     {sm,kp,pip,pim,pi0},   {sm,kp,pi0,pi0,pi0},   {pro,pim,pi0,kp,km},
00801     {pro,pim,pi0,k0,k0b},  {pro,pip,pim,k0,km},   {pro,pi0,pi0,k0,km},
00802     {pro,pim,pim,kp,k0b},  {neu,pip,pim,kp,km},   {neu,pip,pim,k0,k0b},
00803     {neu,pip,pi0,k0,km},   {neu,pim,pi0,kp,k0b},  {neu,pi0,pi0,k0,k0b},
00804     {neu,pi0,pi0,kp,km}}, 
00805 
00806    {{neu,pip,pip,pim,pi0}, {neu,pip,pi0,pi0,pi0}, {pro,pip,pip,pim,pim},
00807     {pro,pip,pim,pi0,pi0}, {pro,pi0,pi0,pi0,pi0}, {lam,kp,pip,pim,pi0},
00808     {lam,k0,pip,pi0,pi0},  {lam,k0,pip,pip,pim},  {lam,kp,pi0,pi0,pi0},
00809     {s0,k0,pip,pip,pim},   {s0,k0,pip,pi0,pi0},   {s0,kp,pip,pim,pi0},
00810     {s0,kp,pi0,pi0,pi0},   {sm,kp,pip,pip,pim},   {sm,kp,pip,pi0,pi0},
00811     {sm,k0,pip,pip,pi0},   {sp,kp,pip,pim,pim},   {sp,kp,pim,pi0,pi0},
00812     {sp,k0,pip,pim,pi0},   {sp,k0,pi0,pi0,pi0},   {neu,pip,pi0,k0,k0b},
00813     {neu,pip,pi0,kp,km},   {neu,pip,pim,kp,k0b},  {neu,pi0,pi0,kp,k0b},
00814     {neu,pip,pip,k0,km},   {pro,pip,pim,k0,k0b},  {pro,pip,pim,kp,km},
00815     {pro,pim,pi0,kp,k0b},  {pro,pip,pi0,k0,km},   {pro,pi0,pi0,kp,km},
00816     {pro,pi0,pi0,k0,k0b}}};
00817 
00818 const G4int G4RPGPionInelastic::T12_6bfs[2][6][6] =
00819 {{{pro,pip,pip,pim,pim,pim}, {pro,pip,pim,pim,pi0,pi0},
00820   {pro,pim,pi0,pi0,pi0,pi0}, {neu,pip,pip,pim,pim,pi0},
00821   {neu,pip,pim,pi0,pi0,pi0}, {neu,pi0,pi0,pi0,pi0,pi0}},
00822 
00823  {{neu,pip,pip,pip,pim,pim}, {neu,pip,pip,pim,pi0,pi0},
00824   {neu,pip,pi0,pi0,pi0,pi0}, {pro,pip,pip,pim,pim,pi0},
00825   {pro,pip,pim,pi0,pi0,pi0}, {pro,pi0,pi0,pi0,pi0,pi0}}};
00826 
00827 const G4int G4RPGPionInelastic::T12_7bfs[2][7][7] =
00828 {{{pro,pip,pip,pim,pim,pim,pi0}, {pro,pip,pim,pim,pi0,pi0,pi0},
00829   {pro,pim,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pip,pip,pim,pim,pim},
00830   {neu,pip,pip,pim,pim,pi0,pi0}, {neu,pip,pim,pi0,pi0,pi0,pi0},
00831   {neu,pi0,pi0,pi0,pi0,pi0,pi0}},
00832 
00833  {{neu,pip,pip,pip,pim,pim,pi0}, {neu,pip,pip,pim,pi0,pi0,pi0},
00834   {neu,pip,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pip,pip,pim,pim,pim},
00835   {pro,pip,pip,pim,pim,pi0,pi0}, {pro,pip,pim,pi0,pi0,pi0,pi0},
00836   {pro,pi0,pi0,pi0,pi0,pi0,pi0}}};
00837 
00838 const G4int G4RPGPionInelastic::T12_8bfs[2][8][8] =
00839 {{{pro,pip,pip,pip,pim,pim,pim,pim}, {pro,pip,pip,pim,pim,pim,pi0,pi0},
00840   {pro,pip,pim,pim,pi0,pi0,pi0,pi0}, {pro,pim,pi0,pi0,pi0,pi0,pi0,pi0},
00841   {neu,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pim,pi0,pi0,pi0,pi0,pi0},
00842   {neu,pip,pip,pim,pim,pi0,pi0,pi0}, {neu,pip,pip,pip,pim,pim,pim,pi0}},
00843 
00844  {{neu,pip,pip,pip,pip,pim,pim,pim}, {neu,pip,pip,pip,pim,pim,pi0,pi0},
00845   {neu,pip,pip,pim,pi0,pi0,pi0,pi0}, {neu,pip,pi0,pi0,pi0,pi0,pi0,pi0},
00846   {pro,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pim,pi0,pi0,pi0,pi0,pi0},
00847   {pro,pip,pip,pim,pim,pi0,pi0,pi0}, {pro,pip,pip,pip,pim,pim,pim,pi0}}};
00848 
00849 const G4int G4RPGPionInelastic::T12_9bfs[2][9][9] =
00850 {{{pro,pip,pip,pip,pim,pim,pim,pim,pi0}, {pro,pip,pip,pim,pim,pim,pi0,pi0,pi0},
00851   {pro,pip,pim,pim,pi0,pi0,pi0,pi0,pi0}, {pro,pim,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
00852   {neu,pi0,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pim,pi0,pi0,pi0,pi0,pi0,pi0},
00853   {neu,pip,pip,pim,pim,pi0,pi0,pi0,pi0}, {neu,pip,pip,pip,pim,pim,pim,pi0,pi0},
00854   {neu,pip,pip,pip,pip,pim,pim,pim,pim}},
00855 
00856  {{neu,pip,pip,pip,pip,pim,pim,pim,pi0}, {neu,pip,pip,pip,pim,pim,pi0,pi0,pi0},
00857   {neu,pip,pip,pim,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
00858   {pro,pi0,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pim,pi0,pi0,pi0,pi0,pi0,pi0},
00859   {pro,pip,pip,pim,pim,pi0,pi0,pi0,pi0}, {pro,pip,pip,pip,pim,pim,pim,pi0,pi0},
00860   {pro,pip,pip,pip,pip,pim,pim,pim,pim}}};
00861 
00862 //
00863 // Cross sections (in mb) for pi- p -> 2-9 body final states
00864 //
00865 // first index:    0-4: channels for mult = 2
00866 //                5-17: channels for mult = 3
00867 //               18-39: channels for mult = 4
00868 //               40-70: channels for mult = 5
00869 //               71-76: channels for mult = 6
00870 //               77-83: channels for mult = 7
00871 //               84-91: channels for mult = 8
00872 //              92-100: channels for mult = 9
00873 //
00874 // second index: kinetic energy
00875 //
00876 const G4float G4RPGPionInelastic::pimPCrossSections[101][30] = {
00877 //
00878 // multiplicity 2 (5 channels)
00879 //
00880 //  pi- p (pi+ n)
00881 //
00882  {  0.0,  1.1,  1.2,  1.4,  1.5,  1.8,  2.0,  3.0,  3.4,  7.0,
00883    14.0, 24.0, 14.7, 10.5, 11.8, 20.0, 14.0, 25.0, 12.0,  9.5,
00884     8.0,  7.0,  6.0,  5.7,  5.0,  4.6,  4.3,  4.0,  3.8,  3.7 },
00885 
00886 //  n pi0  (p pi0)
00887  {  0.0,  2.4,  2.8,  3.3,  4.5,  5.7,  6.3,  9.0, 11.0, 17.0,
00888    30.0, 43.0, 30.0, 16.5, 11.0,  7.0,  4.3,  5.0,  2.0,  0.9,
00889     0.5, 0.24, 0.15,0.094,0.061,0.048,0.035,0.023,0.018,0.014 },
00890 
00891 //  L K0  (L K+)
00892  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00893     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.65, 0.28, 0.16,
00894    0.13,0.075, 0.05,0.032,0.022,0.015,0.011,0.008,0.006,0.004 },
00895 
00896 //  S0 K0  (S0 K+)
00897  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00898     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.28, 0.18, 0.11,
00899   0.091,0.055,0.037,0.025,0.018,0.012,0.008,0.005,0.004,0.003 },
00900 
00901 //  S- K+  (S+ K0)
00902  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00903     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.2, 0.24, 0.09,
00904    0.04,0.012,0.004,0.002,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0 },
00905 //
00906 // multiplicity 3 (13 channels)
00907 //
00908 //  p pi- pi0  (n pi+ pi0)
00909  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00910     0.0,  0.0, 0.05, 0.18, 0.86,  4.4,  5.2,  6.6,  5.4,  4.4,
00911     3.5,  2.5,  2.0,  1.4, 0.97, 0.68, 0.55, 0.36,  0.3, 0.22 },
00912 
00913 //  n pi+ pi-  (p pi+ pi-)
00914  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00915     0.0,  0.0, 0.09,  0.9,  3.3,  6.6,  8.9,  8.2,  8.8,  6.6,
00916     5.2,  3.8,  2.9,  1.9,  1.3,  0.9, 0.75, 0.38, 0.31, 0.22 },
00917 
00918 //  n pi0 pi0  (p pi0 pi0)
00919  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00920     0.0,  0.0, 0.22, 0.74,  1.8,  2.7,  3.0,  2.5,  1.3, 0.64,
00921    0.41, 0.24, 0.15, 0.10,0.065,0.042,  0.0,  0.0,  0.0,  0.0 },
00922 
00923 //  L K0 pi0  (L K+ pi0)
00924  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00925     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.049, 0.18,
00926    0.14, 0.09, 0.07, 0.05, 0.03, 0.02, 0.02, 0.01,0.007,0.004 },
00927 
00928 //  L K+ pi-  (L K0 pi+)
00929  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00930     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.035, 0.12,
00931    0.13, 0.11, 0.08, 0.06, 0.04, 0.02, 0.02, 0.01,0.007,0.004 },
00932 
00933 //  S- K0 pi+  (S+ K+ pi-)
00934  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00935     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.10,
00936    0.13, 0.06, 0.03, 0.02, 0.01,0.003,0.001,  0.0,  0.0,  0.0 },
00937 
00938 //  S- K+ pi0  (S+ K0 pi0)
00939  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00940     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.05,
00941    0.05, 0.03, 0.02, 0.01,0.005,0.002,0.001,  0.0,  0.0,  0.0 },
00942 
00943 //  S+ K0 pi-  (S- K+ pi+)
00944  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00945     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.04,
00946    0.06, 0.04, 0.02, 0.01, 0.01,0.005,0.003,0.001,  0.0,  0.0 },
00947 
00948 //  S0 K+ pi-  (S0 K0 pi+)
00949  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00950     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.04,
00951    0.07, 0.04, 0.03, 0.02, 0.02, 0.01,0.004,0.002,0.001,  0.0 },
00952 
00953 //  S0 K0 pi0  (S0 K+ pi0)
00954  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00955     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.02, 0.09,
00956    0.07, 0.05, 0.03, 0.02, 0.02, 0.01, 0.01,0.005,0.003,0.002 },
00957 
00958 //  p K0 K-  (n K+ K0bar)
00959  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00960     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.03,
00961    0.08, 0.07, 0.05, 0.04, 0.03, 0.02, 0.02, 0.01, 0.01, 0.01 },
00962 
00963 //  n K+ K-  (p K0 K0bar)
00964  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00965     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.04,
00966    0.11, 0.28, 0.12, 0.07, 0.04, 0.02, 0.01,  0.0,  0.0,  0.0 },
00967 
00968 //  n K0 K0bar  (p K+ K-)
00969  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
00970     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 0.06,
00971    0.10, 0.15, 0.18, 0.10, 0.05, 0.02, 0.01,  0.0,  0.0,  0.0 },
00972 //
00973 // multiplicity 4 (22 channels)
00974 //
00975 //  p pi+ pi- pi-  (n pi+ pi+ 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.01, 0.05,  0.2, 0.48, 0.88,  1.8,
00978     2.0,  2.2,  2.0,  1.9,  1.6,  1.4,  1.3,  1.1, 1.05,  1.0 },
00979 
00980 //  p pi- pi0 pi0  (n pi+ pi0 pi0)
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.01, 0.06, 0.18, 0.42,  1.0,  1.9,
00983     2.4,  2.4,  2.0,  1.8,  1.6,  1.4,  1.3,  1.2,  1.1,  1.0 },
00984 
00985 //  n pi+ pi- pi0  (p pi+ pi- pi0)
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.1,  0.3,  0.5,  1.0,  2.0,  2.0,
00988     1.0,  0.8,  0.6,  0.5,  0.4,  0.3, 0.25,  0.2, 0.15,  0.1 },
00989 
00990 //  n pi0 pi0 pi0  (p pi0 pi0 pi0)
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.1, 0.55,  1.0, 0.87,  0.7,  0.6,
00993     0.5,  0.4, 0.35, 0.29, 0.21, 0.14,  0.1, 0.06, 0.03,  0.0 },
00994 
00995 //  L K0 pi+ pi-  (L K+ pi+ pi-)
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.014,
00998    0.08, 0.12, 0.09, 0.08, 0.07, 0.05, 0.05, 0.04, 0.03, 0.03 },
00999 
01000 //  L K0 pi0 pi0  (L K+ pi0 pi0)
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.014,
01003    0.08, 0.12, 0.09, 0.08, 0.07, 0.05, 0.05, 0.04, 0.03, 0.03 },
01004 
01005 //  L K+ pi- pi0  (L K0 pi+ pi0)
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.009,
01008    0.06,  0.1, 0.09, 0.07, 0.06, 0.05, 0.04, 0.03, 0.02, 0.02 },
01009 
01010 //  S0 K0 pi+ pi-  (S0 K+ pi+ pi-)
01011  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01012     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.007,
01013    0.04, 0.06, 0.05, 0.04, 0.04, 0.03, 0.03, 0.02, 0.01, 0.01 },
01014 
01015 //  S0 K0 pi0 pi0  (S0 K+ pi0 pi0)
01016  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01017     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.007,
01018     0.04,0.06, 0.05, 0.04, 0.04, 0.03, 0.03, 0.02, 0.01, 0.01 },
01019 
01020 //  S0 K+ pi- pi0  (S0 K0 pi+ pi0)
01021  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01022     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.005,
01023    0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
01024 
01025 //  S+ K+ pi- pi-  (S- K0 pi+ pi+)
01026  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01027     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.005,
01028    0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
01029 
01030 //  S+ K0 pi- pi0  (S- K+ pi+ pi0)
01031  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01032     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.005,
01033    0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
01034 
01035 //  S- K+ pi+ pi-  (S+ K0 pi+ pi-)
01036  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01037     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.005,
01038    0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
01039 
01040 //  S- K+ pi0 pi0  (S+ K0 pi0 pi0)
01041  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01042     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.005,
01043    0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
01044 
01045 //  S- K0 pi+ pi0  (S+ K+ pi- pi0)
01046  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01047     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.005,
01048    0.03, 0.05, 0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.01, 0.01 },
01049 
01050 //  p pi- K+ K-  (n pi+ K0 K0bar)
01051  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01052     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01053    0.003,0.06,  0.1, 0.09, 0.08, 0.07, 0.07, 0.06, 0.06, 0.05 },
01054 
01055 //  p pi- K0 K0bar  (n pi+ K+ K-)
01056  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01057     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01058    0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
01059 
01060 //  p pi0 K0 K-  (n pi0 K+ K0bar)
01061  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01062     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01063    0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
01064 
01065 //  n pi+ K0 K-  (p pi- K+ K0bar)
01066  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01067     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01068    0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
01069 
01070 //  n pi0 K0 K0bar  (p pi0 K+ K-)
01071  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01072     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01073    0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
01074 
01075 //  n pi0 K+ K-  (p pi0 K0 K0bar)
01076  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01077     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01078    0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
01079 
01080 //  n pi- K+ K0bar  (p pi+ K0 K-)
01081  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01082     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01083    0.03, 0.05, 0.07, 0.06, 0.06, 0.06, 0.06, 0.05, 0.05, 0.05 },
01084 //
01085 // multiplicity 5 (31 channels)
01086 //
01087 //  p pi+ pi- pi- pi0  (n pi+ pi+ pi- 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.007,0.022, 0.10, 0.70,
01090     1.3,  1.9,  2.2,  2.0,  1.7,  1.4,  1.2, 0.90, 0.76, 0.62 },
01091 
01092 //  p pi- pi0 pi0 pi0  (n pi+ pi0 pi0 pi0)
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.007,0.022, 0.10, 0.70,
01095     1.3,  1.9,  2.2,  2.0,  1.7,  1.4,  1.2, 0.90, 0.76, 0.62 },
01096 
01097 //  n pi+ pi+ pi- pi-  (p pi+ pi+ pi- 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.014,0.029, 0.10, 0.30,
01100    0.56, 0.93,  1.2,  1.2,  1.2, 0.94, 0.74, 0.53, 0.40, 0.30 },
01101 
01102 //  n pi+ pi- pi0 pi0  (p pi+ pi- pi0 pi0)
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.006, 0.02, 0.09, 0.63,
01105     1.2,  1.7,  2.0,  1.8,  1.5,  1.3,  1.1, 0.80, 0.70, 0.60 },
01106 
01107 //  n pi0 pi0 pi0 pi0  (p pi0 pi0 pi0 pi0)
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.006, 0.02, 0.08, 0.60,
01110     1.0,  1.5,  2.0,  1.6,  1.4,  1.1,  1.0, 0.70, 0.60, 0.50 },
01111 
01112 //  L K0 pi+ pi- pi0  (L K+ pi+ pi- pi0)
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.007, 0.05, 0.07, 0.08, 0.08, 0.07,0.063, 0.06,0.053,0.048 },
01116 
01117 //  L K+ pi- pi0 pi0  (L K0 pi+ pi0 pi0)
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.003,0.025,0.035, 0.04, 0.04,0.035, 0.03, 0.03,0.025,0.024 },
01121 
01122 //  L K+ pi+ pi- pi-  (L K0 pi+ pi+ pi-)
01123  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01124     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01125   0.002, 0.01,0.036, 0.04, 0.04,0.033,0.031,0.029,0.025,0.022 },
01126 
01127 //  L K0 pi0 pi0 pi0  (L K+ pi0 pi0 pi0)
01128  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01129     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01130   0.001,0.005,0.018, 0.02, 0.02,0.017,0.016,0.015,0.012,0.011 },
01131 
01132 //  S0 K+ pi+ pi- pi-  (S0 K0 pi+ pi+ pi-)
01133  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01134     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01135     0.0,0.003,0.009, 0.01,0.009,0.008,0.007,0.006,0.006,0.005 },
01136 
01137 //  S0 K+ pi- pi0 pi0  (S0 K0 pi+ pi0 pi0)
01138  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01139     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01140     0.0,0.003,0.009, 0.01,0.009,0.008,0.007,0.006,0.006,0.005 },
01141 
01142 //  S0 K0 pi+ pi- pi0  (S0 K+ pi+ pi- pi0)
01143  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01144     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01145   0.003, 0.02, 0.03, 0.04, 0.04,0.035, 0.03, 0.03,0.025, 0.02 },
01146 
01147 //  S0 K0 pi0 pi0 pi0  (S0 K+ pi0 pi0 pi0)
01148  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01149     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01150   0.001, 0.01,0.015, 0.02, 0.02,0.017,0.015,0.015,0.012, 0.01 },
01151 
01152 //  S+ K0 pi+ pi- pi-  (S- K+ pi+ pi+ pi-)
01153  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01154     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01155   0.001,0.002,0.005,0.006,0.006,0.005,0.004,0.004,0.003,0.003 },
01156 
01157 //  S+ K0 pi- pi0 pi0  (S- K+ pi+ pi0 pi0)
01158  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01159     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01160   0.001,0.002,0.005,0.006,0.006,0.005,0.004,0.004,0.003,0.003 },
01161 
01162 //  S+ K+ pi- pi- pi0  (S- K0 pi+ pi+ pi0)
01163  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01164     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01165   0.001,0.002,0.006,0.009,0.008,0.007,0.006,0.005,0.005,0.004 },
01166 
01167 //  S- K0 pi+ pi+ pi-  (S+ K+ pi+ pi- pi-)
01168  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01169     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01170   0.002,0.005,0.015,0.014,0.012,0.009,0.008,0.006,0.005,0.004 },
01171 
01172 //  S- K0 pi+ pi0 pi0  (S+ K+ pi- pi0 pi0)
01173  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01174     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01175   0.001,0.002,0.008,0.007,0.006,0.005,0.004,0.003,0.003,0.002 },
01176 
01177 //  S- K+ pi+ pi- pi0  (S+ K0 pi+ pi- pi0)
01178  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01179     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01180   0.002,0.007,0.023,0.025,0.021,0.018,0.015,0.013,0.010,0.009 },
01181 
01182 //  S- K+ pi0 pi0 pi0  (S+ K0 pi0 pi0 pi0)
01183  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01184     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01185   0.001,0.004,0.011,0.013,0.011,0.009,0.008,0.007,0.006,0.005 },
01186 
01187 //  p pi- pi0 K+ K-  (n pi+ pi0 K0 K0bar)
01188  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01189     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01190     0.0, 0.01, 0.07, 0.07, 0.06,0.055, 0.05,0.047,0.042, 0.04 },
01191 
01192 //  p pi- pi0 K0 K0bar  (n pi+ pi0 K+ K-)
01193  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01194     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01195     0.0,0.003, 0.02, 0.02,0.018,0.017,0.015,0.012,0.011, 0.01 },
01196 
01197 //  p pi+ pi- K0 K-  (n pi+ pi- K+ K0bar)
01198  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01199     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01200     0.0,0.005, 0.02, 0.02,0.018,0.016,0.015,0.013,0.012,0.011 },
01201 
01202 //  p pi0 pi0 K0 K-  (n pi0 pi0 K+ K0bar)
01203  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01204     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01205     0.0,0.005, 0.02, 0.02,0.018,0.016,0.015,0.013,0.012,0.011 },
01206 
01207 //  p pi- pi- K+ K0bar  (n pi+ pi+ K0 K-)
01208  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01209     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01210     0.0,0.007, 0.04, 0.05, 0.04,0.035, 0.03, 0.03,0.027,0.025 },
01211 
01212 //  n pi+ pi- K+ K-  (p pi+ pi- K0 K0bar)
01213  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01214     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01215     0.0, 0.01, 0.07,0.055, 0.05,0.042, 0.04,0.035,0.032,0.029 },
01216 
01217 //  n pi+ pi- K0 K0bar  (p pi+ pi- K+ K-)
01218  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01219     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01220     0.0,0.005,0.035,0.027,0.025,0.021, 0.02,0.017,0.016,0.014 },
01221 
01222 //  n pi+ pi0 K0 K-  (p pi- pi0 K+ K0bar)
01223  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01224     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01225     0.0, 0.01, 0.07,0.055, 0.05,0.042, 0.04,0.035,0.032,0.029 },
01226 
01227 //  n pi- pi0 K+ K0bar  (p pi+ pi0 K0 K-)
01228  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01229     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01230     0.0, 0.01, 0.07,0.055, 0.05,0.042, 0.04,0.035,0.032,0.029 },
01231 
01232 //  n pi0 pi0 K0 K0bar  (p pi0 pi0 K+ K-)
01233  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01234     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01235     0.0,0.005,0.035,0.027,0.025,0.021, 0.02,0.017,0.016,0.014 },
01236 
01237 //  n pi0 pi0 K+ K-  (p pi0 pi0 K0 K0bar)
01238  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01239     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01240     0.0,0.005,0.035,0.027,0.025,0.021, 0.02,0.017,0.016,0.014 },
01241 //
01242 // multiplicity 6 (6 channels)
01243 //
01244 //  p pi+ pi+ pi- pi- pi-  (n pi+ pi+ pi+ pi- pi-)
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.01, 0.03, 0.04,
01247    0.06, 0.11, 0.16, 0.22, 0.31, 0.34,  0.3, 0.24, 0.19, 0.16 },
01248 
01249 //  p pi+ pi- pi- pi0 pi0  (n pi+ pi+ pi- pi0 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.01, 0.04,  0.1,
01252    0.14, 0.19, 0.24, 0.31, 0.37,  0.4, 0.38,  0.4, 0.33, 0.32 },
01253 
01254 //  p pi- pi0 pi0 pi0 pi0  (n pi+ pi0 pi0 pi0 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.01, 0.04,  0.1,
01257    0.14, 0.19, 0.24, 0.31, 0.37,  0.4, 0.38,  0.4, 0.33, 0.32 },
01258 
01259 //  n pi+ pi+ pi- pi- pi0  (p pi+ pi+ pi- pi- pi0)
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.01, 0.04,  0.1,
01262    0.14, 0.19, 0.24, 0.31, 0.37,  0.4, 0.38, 0.4,  0.33, 0.32 },
01263 
01264 //  n pi+ pi- pi0 pi0 pi0  (p pi+ pi- pi0 pi0 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.01, 0.04,  0.1,
01267    0.14, 0.19, 0.24, 0.31, 0.37,  0.4, 0.38,  0.4, 0.33, 0.32 },
01268 
01269 //  n pi0 pi0 pi0 pi0 pi0  (p pi0 pi0 pi0 pi0 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.01, 0.04,  0.1,
01272    0.14, 0.19, 0.24, 0.31, 0.37,  0.4, 0.38,  0.4, 0.33, 0.32 },
01273 //
01274 // multiplicity 7 (7 channels)
01275 //
01276 //  p pi+ pi+ pi- pi- pi- pi0
01277  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01278     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.014, 0.03,
01279    0.05, 0.11, 0.19, 0.37, 0.67, 0.75,  0.7, 0.59, 0.52, 0.47 },
01280 
01281 //  p pi+ pi- pi- pi0 pi0 pi0
01282  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01283     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.008, 0.02,
01284    0.03, 0.07, 0.11, 0.22,  0.4, 0.45,  0.4,  0.4,  0.3,  0.3 },
01285 
01286 //  p pi- pi0 pi0 pi0 pi0 pi0
01287  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01288     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.008, 0.02,
01289    0.03, 0.07, 0.11, 0.22,  0.4, 0.45,  0.4,  0.4,  0.3,  0.3 },
01290 
01291 //  n pi+ pi+ pi+ pi- pi- pi-
01292  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01293     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.002,0.006,
01294   0.012,0.029,0.062, 0.15, 0.29, 0.29, 0.24, 0.18, 0.16, 0.13 },
01295 
01296 //  n pi+ pi+ pi- pi- pi0 pi0
01297  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01298     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.008, 0.02,
01299    0.03, 0.07, 0.11, 0.22,  0.4, 0.45,  0.4,  0.4,  0.3,  0.3 },
01300 
01301 //  n pi+ pi- pi0 pi0 pi0 pi0
01302  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01303     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.008, 0.02,
01304    0.03, 0.07, 0.11, 0.22,  0.4, 0.45,  0.4,  0.4,  0.3,  0.3 },
01305 
01306 //  n pi0 pi0 pi0 pi0 pi0 pi0
01307  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01308     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,0.008, 0.02,
01309    0.03, 0.07, 0.11, 0.22,  0.4, 0.45,  0.4,  0.4,  0.3,  0.3 },
01310 //
01311 // multiplicity 8 (8 channels)
01312 //
01313 //  p pi+ pi+ pi+ pi- pi- pi- pi-  (n pi+ pi+ pi+ pi+ pi- pi- pi-)
01314  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01315     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01316   0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
01317 
01318 //  p pi+ pi+ pi- pi- pi- pi0 pi0  (n pi+ pi+ pi+ pi- pi- pi0 pi0)
01319  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01320     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01321   0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
01322 
01323 //  p pi+ pi- pi- pi0 pi0 pi0 pi0  (n pi+ pi+ pi- pi0 pi0 pi0 pi0)
01324  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01325     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01326   0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
01327 
01328 //  p pi- pi0 pi0 pi0 pi0 pi0 pi0  (n pi+ pi0 pi0 pi0 pi0 pi0 pi0)
01329  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01330     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01331   0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
01332 
01333 //  n pi0 pi0 pi0 pi0 pi0 pi0 pi0  (p pi0 pi0 pi0 pi0 pi0 pi0 pi0)
01334  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01335     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01336   0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
01337 
01338 //  n pi+ pi- pi0 pi0 pi0 pi0 pi0  (p pi+ pi- pi0 pi0 pi0 pi0 pi0)
01339  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01340     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01341   0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
01342 
01343 //  n pi+ pi+ pi- pi- pi0 pi0 pi0  (p pi+ pi+ pi- pi- pi0 pi0 pi0)
01344  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01345     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01346   0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
01347 
01348 //  n pi+ pi+ pi+ pi- pi- pi- pi0  (p pi+ pi+ pi+ pi- pi- pi- pi0)
01349  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01350     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01351   0.002,0.008, 0.02, 0.06, 0.09, 0.11, 0.11, 0.10, 0.10, 0.09 },
01352 //
01353 // multiplicity 9 (9 channels)
01354 //
01355 //  p pi+ pi+ pi+ pi- pi- pi- pi- pi0
01356  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01357     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01358     0.0,  0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
01359 
01360 //  p pi+ pi+ pi- pi- pi- pi0 pi0 pi0
01361  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01362     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01363     0.0,  0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
01364 
01365 //  p pi+ pi- pi- pi0 pi0 pi0 pi0 pi0
01366  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01367     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01368     0.0,  0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
01369 
01370 //  p pi- pi0 pi0 pi0 pi0 pi0 pi0 pi0
01371  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01372     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01373     0.0,  0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
01374 
01375 //  n pi0 pi0 pi0 pi0 pi0 pi0 pi0 pi0
01376  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01377     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01378     0.0,  0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
01379 
01380 //  n pi+ pi- pi0 pi0 pi0 pi0 pi0 pi0
01381  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01382     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01383     0.0,  0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
01384 
01385 //  n pi+ pi+ pi- pi- pi0 pi0 pi0 pi0
01386  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01387     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01388     0.0,  0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
01389 
01390 //  n pi+ pi+ pi+ pi- pi- pi- pi0 pi0
01391  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01392     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01393     0.0,  0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 },
01394 
01395 //  n pi+ pi+ pi+ pi+ pi- pi- pi- pi-
01396  {  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01397     0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,
01398     0.0,  0.0,0.012,0.036,0.084, 0.14, 0.18, 0.18, 0.18, 0.17 } };
01399 
01400  /* end of file */
01401  

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