Geant4: G4CascadeT1GamNChannel.cc Source 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 // GEANT4 tag: $Name: not supported by cvs2svn $
00028 //
00029 
00030 #include "G4CascadeGamPChannel.hh"
00031 #include "G4CascadeGamNChannel.hh"
00032 #include "G4InuclParticleNames.hh"
00033 using namespace G4InuclParticleNames;
00034 
00035 namespace {
00036   // gamma p : Outgoing particle types of a given multiplicity
00037   static const G4int gamp2bfs[6][2] =
00038   {{gam,pro}, {pro,pi0}, {neu,pip}, {lam,kpl}, {s0, kpl}, {sp,k0}};
00039 
00040   static const G4int gamp3bfs[6][3] =
00041   {{pro,pip,pim}, {pro,pi0,pi0}, {neu,pip,pi0},
00042    {pro,kpl,kmi},  {pro,k0,k0b},  {neu,kpl,k0b}};
00043 
00044   static const G4int gamp4bfs[4][4] =
00045   {{pro,pip,pim,pi0}, {pro,pi0,pi0,pi0}, {neu,pip,pip,pim}, {neu,pip,pi0,pi0}};
00046 
00047   static const G4int gamp5bfs[5][5] =
00048   {{pro,pip,pip,pim,pim}, {pro,pip,pim,pi0,pi0}, {pro,pi0,pi0,pi0,pi0},
00049    {neu,pip,pip,pim,pi0}, {neu,pip,pi0,pi0,pi0}};
00050 
00051   static const G4int gamp6bfs[6][6] =
00052   {{pro,pip,pip,pim,pim,pi0}, {pro,pip,pim,pi0,pi0,pi0}, {pro,pi0,pi0,pi0,pi0,pi0},
00053    {neu,pip,pip,pip,pim,pim}, {neu,pip,pip,pim,pi0,pi0}, {neu,pip,pi0,pi0,pi0,pi0}};
00054 
00055   static const G4int gamp7bfs[7][7] =
00056   {{pro,pip,pip,pip,pim,pim,pim}, {pro,pip,pip,pim,pim,pi0,pi0},
00057    {pro,pip,pim,pi0,pi0,pi0,pi0}, {pro,pi0,pi0,pi0,pi0,pi0,pi0},
00058    {neu,pip,pip,pim,pi0,pi0,pi0}, {neu,pip,pip,pip,pim,pim,pi0},
00059    {neu,pip,pi0,pi0,pi0,pi0,pi0}};
00060 
00061   static const G4int gamp8bfs[7][8] =
00062   {{pro,pip,pip,pim,pim,pi0,pi0,pi0}, {pro,pip,pip,pip,pim,pim,pim,pi0},
00063    {pro,pip,pim,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pip,pip,pim,pim,pi0,pi0},
00064    {neu,pip,pip,pim,pi0,pi0,pi0,pi0}, {neu,pip,pip,pip,pip,pim,pim,pim},
00065    {neu,pip,pi0,pi0,pi0,pi0,pi0,pi0}};
00066 
00067   static const G4int gamp9bfs[7][9] =
00068   {{pro,pip,pip,pip,pim,pim,pim,pi0,pi0}, {pro,pip,pip,pim,pim,pi0,pi0,pi0,pi0},
00069    {pro,pip,pip,pip,pip,pim,pim,pim,pim}, {pro,pip,pim,pi0,pi0,pi0,pi0,pi0,pi0},
00070    {neu,pip,pip,pip,pim,pim,pi0,pi0,pi0}, {neu,pip,pip,pip,pip,pim,pim,pim,pi0},
00071    {neu,pip,pip,pim,pi0,pi0,pi0,pi0,pi0}};
00072 }
00073 
00074 namespace {
00075   // gamma n : Outgoing particle types of a given multiplicity
00076   static const G4int gamn2bfs[6][2] =
00077   {{gam,neu}, {neu,pi0}, {pro,pim}, {lam,k0}, {s0, k0}, {sm,kpl}};
00078 
00079   static const G4int gamn3bfs[6][3] =
00080   {{neu,pip,pim}, {neu,pi0,pi0}, {pro,pim,pi0},
00081    {neu,k0,k0b},  {neu,kpl,kmi},   {pro,k0,kmi}};
00082 
00083   static const G4int gamn4bfs[4][4] =
00084   {{neu,pip,pim,pi0}, {neu,pi0,pi0,pi0}, {pro,pip,pim,pim}, {pro,pim,pi0,pi0}};
00085 
00086   static const G4int gamn5bfs[5][5] =
00087   {{neu,pip,pip,pim,pim}, {neu,pip,pim,pi0,pi0}, {neu,pi0,pi0,pi0,pi0},
00088    {pro,pip,pim,pim,pi0}, {pro,pim,pi0,pi0,pi0}};
00089 
00090   static const G4int gamn6bfs[6][6] =
00091   {{neu,pip,pip,pim,pim,pi0}, {neu,pip,pim,pi0,pi0,pi0}, {neu,pi0,pi0,pi0,pi0,pi0},
00092    {pro,pip,pip,pim,pim,pim}, {pro,pip,pim,pim,pi0,pi0}, {pro,pim,pi0,pi0,pi0,pi0}};
00093 
00094   static const G4int gamn7bfs[7][7] =
00095   {{neu,pip,pip,pip,pim,pim,pim}, {neu,pip,pip,pim,pim,pi0,pi0},
00096    {neu,pip,pim,pi0,pi0,pi0,pi0}, {neu,pi0,pi0,pi0,pi0,pi0,pi0},
00097    {pro,pip,pim,pim,pi0,pi0,pi0}, {pro,pip,pip,pim,pim,pim,pi0},
00098    {pro,pim,pi0,pi0,pi0,pi0,pi0}};
00099 
00100   static const G4int gamn8bfs[7][8] =
00101   {{neu,pip,pip,pim,pim,pi0,pi0,pi0}, {neu,pip,pip,pip,pim,pim,pim,pi0},
00102    {neu,pip,pim,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pip,pim,pim,pim,pi0,pi0},
00103    {pro,pip,pim,pim,pi0,pi0,pi0,pi0}, {pro,pip,pip,pip,pim,pim,pim,pim},
00104    {pro,pim,pi0,pi0,pi0,pi0,pi0,pi0}};
00105 
00106   static const G4int gamn9bfs[7][9] =
00107   {{neu,pip,pip,pip,pim,pim,pim,pi0,pi0}, {neu,pip,pip,pim,pim,pi0,pi0,pi0,pi0},
00108    {neu,pip,pip,pip,pip,pim,pim,pim,pim}, {neu,pip,pim,pi0,pi0,pi0,pi0,pi0,pi0},
00109    {pro,pip,pip,pim,pim,pim,pi0,pi0,pi0}, {pro,pip,pip,pip,pim,pim,pim,pim,pi0},
00110    {pro,pip,pim,pim,pi0,pi0,pi0,pi0,pi0}};
00111 }
00112 
00113 
00114 namespace {
00115   // Total gamma p cross section as a function of kinetic energy
00116   static const G4double gampTotXSec[30] =
00117    {0.0002, 0.0002, 0.0002, 0.0002, 0.0002, 0.0002, 0.0002, 0.0002, 0.0002,
00118     0.0002, 0.0002, 0.0723, 0.2327, 0.54,   0.24,   0.19,   0.28,   0.21,
00119     0.18,   0.146,  0.138,  0.13,   0.1277, 0.124,  0.122,  0.12,   0.1185,
00120     0.117,  0.115,  0.115};
00121 
00122   static const G4double gampCrossSections[48][30] = {
00123   //
00124   // multiplicity 2 (6 channels)
00125   //
00126   // gamma p elastic ( = gamma n elastic )
00127    {0.0002,  0.0002, 0.0002, 0.0002, 0.0002, 0.0002, 0.0002, 0.0002, 0.0002,
00128     0.0002,  0.0002, 0.0003, 0.0007, 0.0027, 0.0014, 0.0010, 0.0012, 0.0005,
00129     0.0003,  0.0002, 0.0002, 0.0002, 0.0002, 0.0002, 0.0001, 0.0001, 0.0001,
00130     0.0001,  0.0001, 0.0001},
00131  
00132   // p pi0 ( = n pi0 )
00133    {0.0,   0.0,   0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,   0.0,
00134     0.0,   0.01,  0.08,   0.3,    0.1143, 0.0226, 0.0122, 0.0106, 0.009, 0.006,
00135     0.004, 0.003, 0.0017, 0.0009, 0.0005, 0.0002, 0.0,    0.0,    0.0,   0.0}, 
00136 
00137   // n pi+ ( = p pi- )
00138    {0.0,   0.0,   0.0,    0.0,    0.0,    0.0,    0.0,   0.0,    0.0,   0.0,
00139     0.0,   0.062, 0.152,  0.2373, 0.1143, 0.093,  0.105, 0.0228, 0.011, 0.006,
00140     0.003, 0.003, 0.0017, 0.0009, 0.0005, 0.0001, 0.0,   0.0,    0.0,   0.0},
00141 
00142   // L K+ ( = L K0 )
00143    {0.0,    0.0,    0.0,    0.0,    0.0,   0.0,  0.0,  0.0,    0.0,   0.0,
00144     0.0,    0.0,    0.0,    0.0,    0.0,   0.0,  0.0,  0.0014, 0.002, 0.0007,
00145     0.0004, 0.0003, 0.0002, 0.0001, 0.0,   0.0,  0.0,  0.0,    0.0,   0.0}, 
00146 
00147   // S0 K+ ( = S0 K0 )
00148    {0.0,    0.0,    0.0,    0.0,    0.0,  0.0,  0.0,  0.0, 0.0,    0.0,
00149     0.0,    0.0,    0.0,    0.0,    0.0,  0.0,  0.0,  0.0, 0.0012, 0.0006,
00150     0.0003, 0.0002, 0.0001, 0.0001, 0.0,  0.0,  0.0,  0.0, 0.0,    0.0},
00151  
00152   // S+ K0  ( = S- K+ )
00153    {0.0,    0.0,    0.0,    0.0,    0.0,   0.0,  0.0,  0.0, 0.0,    0.0,
00154     0.0,    0.0,    0.0,    0.0,    0.0,   0.0,  0.0,  0.0, 0.0015, 0.0006,
00155     0.0003, 0.0002, 0.0001, 0.0001, 0.0,   0.0,  0.0,  0.0, 0.0,    0.0},
00156   //
00157   // multiplicity 3 (6 channels)
00158   //
00159   // p pi+ pi- ( = n pi+ pi- )
00160    {0.0,   0.0,   0.0,   0.0,   0.0,   0.0,   0.0,    0.0,   0.0,   0.0,
00161     0.0,   0.0,   0.0,   0.0,   0.004, 0.05,  0.077, 0.077,  0.061, 0.040,
00162     0.033, 0.026, 0.022, 0.018, 0.016, 0.015, 0.014, 0.0135, 0.013, 0.012},
00163 
00164   // p pi0 pi0 ( = n pi0 pi0 )
00165    {0.0,  0.0,    0.0,    0.0,    0.0,   0.0,    0.0,    0.0,    0.0,   0.0,
00166     0.0,  0.0,    0.0,    0.0,    0.002, 0.008,  0.0317, 0.0318, 0.026, 0.016,
00167     0.01, 0.0089, 0.0076, 0.0071, 0.006, 0.0065, 0.0066, 0.007,  0.006, 0.006},
00168 
00169   // n pi+ pi0 ( = p pi- pi0 ) 
00170    {0.0,   0.0,   0.0,    0.0,    0.0,    0.0,    0.0,    0.0,   0.0,   0.0,
00171     0.0,   0.0,   0.0,    0.0,    0.004,  0.0154, 0.0488, 0.057, 0.044, 0.029,
00172     0.017, 0.013, 0.0081, 0.0043, 0.0017, 0.0006, 0.0002, 0.0,   0.0,   0.0},
00173 
00174   // p K+ K- ( = n K0 K0b )
00175    {0.0,    0.0,   0.0,   0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00176     0.0,    0.0,   0.0,   0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00177     0.0004, 0.001, 0.001, 0.0008, 0.0007, 0.0006, 0.0005, 0.0003, 0.0002, 0.0001},   
00178 
00179   // p K0 K0b  (= n K+ K- )
00180    {0.0,    0.0,   0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00181     0.0,    0.0,   0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00182     0.0004, 0.001, 0.0009, 0.0006, 0.0007, 0.0006, 0.0005, 0.0003, 0.0002, 0.0001},  
00183 
00184   // n K+ K0b  (= p K0 K- )
00185    {0.0,    0.0,   0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00186     0.0,    0.0,   0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00187     0.0004, 0.001, 0.0009, 0.0006, 0.0007, 0.0006, 0.0005, 0.0003, 0.0002, 0.0001},
00188   //
00189   // multiplicity 4 (4 channels)
00190   //
00191   // p pi+ pi- pi0 ( = n pi+ pi- pi0 )
00192    {0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00193     0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0023, 0.0049, 0.014,  0.024,
00194     0.025,  0.0225, 0.0196, 0.0162, 0.0125, 0.0088, 0.0059, 0.003,  0.0013, 0.0004},
00195 
00196   // p pi0 pi0 pi0 ( = n pi0 pi0 pi0 )
00197    {0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00198     0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0004, 0.0008, 0.0016, 0.0025,
00199     0.004,  0.0024, 0.0023, 0.0021, 0.0016, 0.0014, 0.0009, 0.0005, 0.0002, 0.0001},
00200 
00201   // n pi+ pi+ pi- ( = p pi+ pi- pi- )
00202    {0.0,   0.0,   0.0,    0.0,    0.0,   0.0,    0.0,    0.0,    0.0,    0.0,
00203     0.0,   0.0,   0.0,    0.0,    0.0,   0.0,    0.0007, 0.0016, 0.004,  0.0088,
00204     0.011, 0.010, 0.0081, 0.0064, 0.0046, 0.003, 0.0018, 0.0008, 0.0003, 0.0001},
00205 
00206   // n pi+ pi0 pi0 (= p pi- pi0 pi0 )
00207    {0.0,    0.0,   0.0,    0.0,   0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00208     0.0,    0.0,   0.0,    0.0,   0.0,    0.0,    0.0007, 0.0016, 0.003,  0.006,
00209     0.0106, 0.007, 0.0056, 0.005, 0.0035, 0.0025, 0.0017, 0.0008, 0.0003, 0.0001},
00210   //
00211   // multiplicity 5 (5 channels)
00212   //
00213   // p pi+ pi+ pi- pi- ( = n pi+ pi+ pi- pi- )
00214    {0.0,    0.0,    0.0,    0.0,    0.0,   0.0,    0.0,    0.0,    0.0,    0.0,
00215     0.0,    0.0,    0.0,    0.0,    0.0,   0.0,    0.0,    0.0,    0.0004, 0.0012,
00216     0.0032, 0.0049, 0.0055, 0.0052, 0.005, 0.0048, 0.0046, 0.0042, 0.0038, 0.0033},
00217 
00218   // p pi+ pi- pi0 pi0 ( = n pi+ pi- pi0 pi0 )
00219    {0.0,    0.0,    0.0,    0.0,    0.0,    0.0,   0.0,    0.0,    0.0,    0.0,
00220     0.0,    0.0,    0.0,    0.0,    0.0,    0.0,   0.0,    0.0,    0.0005, 0.0016,
00221     0.0038, 0.0064, 0.0076, 0.0081, 0.0075, 0.008, 0.0086, 0.0084, 0.0076, 0.0066},   
00222 
00223   // p pi0 pi0 pi0 pi0 ( = n pi0 pi0 pi0 pi0 )
00224    {0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00225     0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0001,
00226     0.0005, 0.0005, 0.0006, 0.0007, 0.0007, 0.0007, 0.0008, 0.0007, 0.0006, 0.0005},   
00227 
00228   // n pi+ pi+ pi- pi0 ( = p pi+ pi- pi- pi0 )
00229    {0.0,   0.0,    0.0,    0.0,    0.0,   0.0,    0.0,    0.0,    0.0,    0.0,
00230     0.0,   0.0,    0.0,    0.0,    0.0,   0.0,    0.0,    0.0,    0.0005, 0.0016,
00231     0.006, 0.0064, 0.0076, 0.0081, 0.009, 0.0093, 0.0086, 0.0084, 0.0076, 0.0066},    
00232 
00233   // n pi+ pi0 pi0 pi0 ( = p pi- pi0 pi0 pi0 )
00234    {0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00235     0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0006,
00236     0.0013, 0.0022, 0.0026, 0.0027, 0.0025, 0.0027, 0.0028, 0.0028, 0.0026, 0.0022},
00237   //
00238   // multiplicity 6 (6 channels)
00239   //
00240   // p pi+ pi+ pi- pi- pi0 ( = n pi+ pi+ pi- pi- pi0 )
00241    {0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00242     0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0002,
00243     0.0014, 0.0038, 0.0062, 0.0083, 0.0090, 0.0086, 0.0082, 0.0075, 0.0068, 0.0059},
00244 
00245   // p pi+ pi- pi0 pi0 pi0 ( = n pi+ pi- pi0 pi0 pi0 )
00246    {0.0,    0.0,    0.0,    0.0,   0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00247     0.0,    0.0,    0.0,    0.0,   0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00248     0.0009, 0.0017, 0.0029, 0.005, 0.0045, 0.0049, 0.0051, 0.0050, 0.0046, 0.0040},
00249 
00250   // p pi0 pi0 pi0 pi0 pi0 ( = n pi0 pi0 pi0 pi0 pi0 )
00251    {0.0,  0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00252     0.0,  0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00253     0.0,  0.0001, 0.0002, 0.0002, 0.0002, 0.0003, 0.0003, 0.0002, 0.0002, 0.0002},
00254 
00255   // n pi+ pi+ pi+ pi- pi- ( = p pi+ pi+ pi- pi- pi- )
00256    {0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00257     0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0001,
00258     0.0003, 0.0009, 0.0026, 0.0032, 0.0035, 0.0033, 0.0031, 0.0028, 0.0025, 0.0022},
00259 
00260   // n pi+ pi+ pi- pi0 pi0 ( = p pi+ pi- pi- pi0 pi0 )
00261    {0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00262     0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0002,
00263     0.0005, 0.0017, 0.0054, 0.0073, 0.0079, 0.0085, 0.0088, 0.0084, 0.0075, 0.0066},
00264 
00265   // n pi+ pi0 pi0 pi0 pi0 ( = p pi- pi0 pi0 pi0 pi0 )
00266    {0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00267     0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00268     0.0001, 0.0003, 0.0009, 0.0012, 0.0013, 0.0013, 0.0014, 0.0014, 0.0012, 0.0011},
00269   //
00270   // multiplicity 7 (7 channels)
00271   //
00272   // p pi+ pi+ pi+ pi- pi- pi- ( = n pi+ pi+ pi+ pi- pi- pi- )
00273    {0.0,  0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00274     0.0,  0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00275     0.0,  0.0001, 0.0004, 0.0006, 0.0009, 0.0009, 0.0009, 0.0009, 0.0008, 0.0008},
00276 
00277   // p pi+ pi+ pi- pi- pi0 pi0 ( = n pi+ pi+ pi- pi- pi0 pi0 )
00278    {0.0,  0.0,    0.0,    0.0,    0.0,   0.0,   0.0,    0.0,   0.0,    0.0,
00279     0.0,  0.0,    0.0,    0.0,    0.0,   0.0,   0.0,    0.0,   0.0,    0.0,
00280     0.0,  0.0003, 0.0012, 0.0021, 0.004, 0.004, 0.0047, 0.005, 0.0045, 0.004},
00281 
00282   // p pi+ pi- pi0 pi0 pi0 pi0 ( = n pi+ pi- pi0 pi0 pi0 pi0 )
00283    {0.0,  0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00284     0.0,  0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00285     0.0,  0.0001, 0.0004, 0.0007, 0.0013, 0.0014, 0.0016, 0.0017, 0.0015, 0.0013},
00286 
00287   // p pi0 pi0 pi0 pi0 pi0 pi0 ( = n pi0 pi0 pi0 pi0 pi0 pi0 )
00288    {0.0, 0.0, 0.0, 0.0, 0.0, 0.0,    0.0,    0.0,    0.0, 0.0,
00289     0.0, 0.0, 0.0, 0.0, 0.0, 0.0,    0.0,    0.0,    0.0, 0.0,
00290     0.0, 0.0, 0.0, 0.0, 0.0, 0.0001, 0.0001, 0.0001, 0.0, 0.0},
00291 
00292   // n pi+ pi+ pi- pi0 pi0 pi0 ( = p pi+ pi- pi- pi0 pi0 pi0 )
00293    {0.0,  0.0,    0.0,    0.0,    0.0,   0.0,    0.0,    0.0,    0.0,    0.0,
00294     0.0,  0.0,    0.0,    0.0,    0.0,   0.0,    0.0,    0.0,    0.0,    0.0,
00295     0.0,  0.0003, 0.0008, 0.0014, 0.002, 0.0023, 0.0025, 0.0027, 0.0025, 0.0025},
00296 
00297   // n pi+ pi+ pi+ pi- pi- pi0 ( = p pi+ pi+ pi- pi- pi- pi0 )
00298    {0.0,  0.0,    0.0,    0.0,   0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00299     0.0,  0.0,    0.0,    0.0,   0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00300     0.0,  0.0003, 0.0008, 0.0014, 0.002, 0.0023, 0.0025, 0.0027, 0.0025, 0.0025},
00301 
00302   // n pi+ pi0 pi0 pi0 pi0 pi0 ( = p pi- pi0 pi0 pi0 pi0 pi0 )
00303    {0.0, 0.0,  0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00304     0.0, 0.0,  0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00305     0.0, 0.0,  0.0001, 0.0002, 0.0002, 0.0003, 0.0003, 0.0003, 0.0002, 0.0002},
00306   //
00307   // multiplicity 8 (7 channels)
00308   //
00309   // p pi+ pi+ pi- pi- pi0 pi0 pi0 ( = n pi+ pi+ pi- pi- pi0 pi0 pi0 )
00310    {0.0,  0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00311     0.0,  0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00312     0.0,  0.0001, 0.0003, 0.0006, 0.0015, 0.0026, 0.0033, 0.0034, 0.0034, 0.0033},
00313 
00314   // p pi+ pi+ pi+ pi- pi- pi- pi0 ( = n pi+ pi+ pi+ pi- pi- pi- pi0 )
00315    {0.0,  0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00316     0.0,  0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00317     0.0,  0.0001, 0.0002, 0.0005, 0.0011, 0.0017, 0.0023, 0.0022, 0.0022, 0.0022},
00318 
00319   // p pi+ pi- pi0 pi0 pi0 pi0 pi0 ( = n pi+ pi- pi0 pi0 pi0 pi0 pi0 )
00320    {0.0, 0.0,  0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00321     0.0, 0.0,  0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00322     0.0, 0.0,  0.0001, 0.0002, 0.0004, 0.0005, 0.0007, 0.0007, 0.0007, 0.0007},
00323 
00324   // p pi0 pi0 pi0 pi0 pi0 pi0 pi0 ( = n pi0 pi0 pi0 pi0 pi0 pi0 pi0 )
00325   // negligible 
00326 
00327   // n pi+ pi+ pi+ pi- pi- pi0 pi0 ( = p pi+ pi+ pi- pi- pi- pi0 pi0 )
00328    {0.0,  0.0,    0.0,    0.0,    0.0,   0.0,    0.0,    0.0,    0.0,    0.0,
00329     0.0,  0.0,    0.0,    0.0,    0.0,   0.0,    0.0,    0.0,    0.0,    0.0,
00330     0.0,  0.0001, 0.0003, 0.0006, 0.002, 0.0026, 0.0033, 0.0034, 0.0034, 0.0033},
00331 
00332   // n pi+ pi+ pi- pi0 pi0 pi0 pi0 ( = p pi+ pi- pi- pi0 pi0 pi0 pi0 )
00333    {0.0, 0.0, 0.0,    0.0,   0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00334     0.0, 0.0, 0.0,    0.0,   0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00335     0.0, 0.0, 0.0001, 0.0003, 0.001, 0.0013, 0.0017, 0.0017, 0.0017, 0.0016},
00336 
00337   // n pi+ pi+ pi+ pi+ pi- pi- pi- ( = p pi+ pi+ pi+ pi- pi- pi- pi- )
00338    {0.0, 0.0, 0.0,  0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00339     0.0, 0.0, 0.0,  0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00340     0.0, 0.0, 0.0,  0.0001, 0.0002, 0.0004, 0.0006, 0.0005, 0.0005, 0.0005},
00341 
00342   // n pi+ pi0 pi0 pi0 pi0 pi0 pi0 ( = p pi- pi0 pi0 pi0 pi0 pi0 pi0 )
00343    {0.0, 0.0, 0.0, 0.0,  0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00344     0.0, 0.0, 0.0, 0.0,  0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00345     0.0, 0.0, 0.0, 0.0,  0.0001, 0.0001, 0.0001, 0.0001, 0.0001, 0.0001},
00346   //
00347   // multiplicity 9 (7 channels)
00348   //
00349   // p pi+ pi+ pi+ pi- pi- pi- pi0 pi0 ( = n pi+ pi+ pi+ pi- pi- pi- pi0 pi0 )
00350    {0.0, 0.0,  0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00351     0.0, 0.0,  0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00352     0.0, 0.0,  0.0002, 0.0005, 0.0013, 0.0019, 0.0023, 0.0026, 0.0026, 0.0024},
00353 
00354   // p pi+ pi+ pi- pi- pi0 pi0 pi0 pi0 ( = n pi+ pi+ pi- pi- pi0 pi0 pi0 pi0 )
00355    {0.0, 0.0,  0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00356     0.0, 0.0,  0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00357     0.0, 0.0,  0.0002, 0.0003, 0.0011, 0.0014, 0.0018, 0.0020, 0.0020, 0.0018},
00358 
00359   // p pi+ pi+ pi+ pi+ pi- pi- pi- pi- ( = n pi+ pi+ pi+ pi+ pi- pi- pi- pi- )
00360    {0.0, 0.0, 0.0,  0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
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.0001, 0.0002, 0.0003, 0.0003, 0.0003, 0.0003, 0.0003},
00363 
00364   // p pi+ pi- pi0 pi0 pi0 pi0 pi0 pi0 ( = n pi+ pi- pi0 pi0 pi0 pi0 pi0 pi0 )
00365    {0.0, 0.0, 0.0,  0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00366     0.0, 0.0, 0.0,  0.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00367     0.0, 0.0, 0.0,  0.0001, 0.0001, 0.0002, 0.0003, 0.0002, 0.0002, 0.0002},
00368 
00369   // p pi0 pi0 pi0 pi0 pi0 pi0 pi0 pi0 ( = n pi0 pi0 pi0 pi0 pi0 pi0 pi0 pi0 )
00370   // negligible
00371 
00372   // n pi+ pi+ pi+ pi- pi- pi0 pi0 pi0 ( = p pi+ pi+ pi- pi- pi- pi0 pi0 pi0 )
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.0,    0.0,    0.0,    0.0,    0.0,    0.0,    0.0,
00375     0.0, 0.0,  0.0002, 0.0005, 0.0013, 0.0018, 0.0023, 0.0026, 0.0026, 0.0024},
00376 
00377   // n pi+ pi+ pi+ pi+ pi- pi- pi- pi0 ( = p pi+ pi+ pi+ pi- pi- pi- pi- 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.0,    0.0,
00380     0.0, 0.0,  0.0001, 0.0002, 0.0007, 0.0009, 0.0011, 0.0013, 0.0013, 0.0012},
00381 
00382   // n pi+ pi+ pi- pi0 pi0 pi0 pi0 pi0 ( = p pi+ pi- pi- pi0 pi0 pi0 pi0 pi0 )
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.0,    0.0,
00385     0.0, 0.0,  0.0001, 0.0002, 0.0004, 0.0006, 0.0008, 0.0009, 0.0009, 0.0008}};
00386 
00387   // n pi+ pi0 pi0 pi0 pi0 pi0 pi0 pi0 ( = p pi- pi0 pi0 pi0 pi0 pi0 pi0 pi0 )
00388   // negligible
00389 }
00390 
00391 // Initialize both |T Tz> = |1/2 +-1/2> channels using gamma-p cross section
00392 // tables
00393 
00394 G4CascadeGamPChannelData::data_t
00395 G4CascadeGamPChannelData::data(gamp2bfs, gamp3bfs, gamp4bfs, gamp5bfs,
00396                                gamp6bfs, gamp7bfs, gamp8bfs, gamp9bfs,
00397                                gampCrossSections, gampTotXSec, gam*pro,
00398                                "GammaProton");
00399 
00400 G4CascadeGamNChannelData::data_t
00401 G4CascadeGamNChannelData::data(gamn2bfs, gamn3bfs, gamn4bfs, gamn5bfs,
00402                                gamn6bfs, gamn7bfs, gamn8bfs, gamn9bfs,
00403                                gampCrossSections, gampTotXSec, gam*neu,
00404                                "GammaNeutron");
00405 
00406