Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/hadronic/models/im_r_matrix/src/G4XNNstarTable.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /processes/hadronic/models/im_r_matrix/src/G4XNNstarTable.cc (Version 11.3.0) and /processes/hadronic/models/im_r_matrix/src/G4XNNstarTable.cc (Version 10.2.p3)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                         3 // * License and Disclaimer                                           *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th      5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided      6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License      7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/      8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.           9 // * include a list of copyright holders.                             *
 10 // *                                               10 // *                                                                  *
 11 // * Neither the authors of this software syst     11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file      15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitatio     16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                               17 // *                                                                  *
 18 // * This  code  implementation is the result      18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboratio     19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distri     20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  ag     21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publicati     22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Sof     23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // *******************************************     24 // ********************************************************************
 25 //                                                 25 //
 26 //                                                 26 //
 27                                                    27 
 28 #include "globals.hh"                              28 #include "globals.hh"
 29 #include "G4ios.hh"                                29 #include "G4ios.hh"
 30 #include "G4SystemOfUnits.hh"                      30 #include "G4SystemOfUnits.hh"
 31 #include "G4XNNstarTable.hh"                       31 #include "G4XNNstarTable.hh"
 32 #include "G4PhysicsFreeVector.hh"                  32 #include "G4PhysicsFreeVector.hh"
 33                                                    33 
 34 const G4int G4XNNstarTable::sizeNNstar = 121;      34 const G4int G4XNNstarTable::sizeNNstar = 121;
 35                                                    35 
 36 // Energies (GeV) corresponding to the cross s     36 // Energies (GeV) corresponding to the cross section table
 37 // Units are assigned when filling the Physics     37 // Units are assigned when filling the PhysicsVector
 38                                                    38 
 39 const G4double G4XNNstarTable::energyTable[121     39 const G4double G4XNNstarTable::energyTable[121] =
 40 {                                                  40 {
 41   0.0,                                             41   0.0, 
 42   2.014,  2.014,  2.016,  2.018,  2.022,  2.02     42   2.014,  2.014,  2.016,  2.018,  2.022,  2.026,  2.031,  2.037,  2.044,  2.052,   
 43   2.061,  2.071,  2.082,  2.094,  2.107,  2.12     43   2.061,  2.071,  2.082,  2.094,  2.107,  2.121,  2.135,  2.151,  2.168,  2.185,   
 44   2.204,  2.223,  2.244,  2.265,  2.287,  2.31     44   2.204,  2.223,  2.244,  2.265,  2.287,  2.311,  2.335,  2.360,  2.386,  2.413,   
 45   2.441,  2.470,  2.500,  2.531,  2.562,  2.59     45   2.441,  2.470,  2.500,  2.531,  2.562,  2.595,  2.629,  2.664,  2.699,  2.736,
 46   2.773,  2.812,  2.851,  2.891,  2.933,  2.97     46   2.773,  2.812,  2.851,  2.891,  2.933,  2.975,  3.018,  3.062,  3.107,  3.153,   
 47   3.200,  3.248,  3.297,  3.347,  3.397,  3.44     47   3.200,  3.248,  3.297,  3.347,  3.397,  3.449,  3.502,  3.555,  3.610,  3.666,   
 48   3.722,  3.779,  3.838,  3.897,  3.957,  4.01     48   3.722,  3.779,  3.838,  3.897,  3.957,  4.018,  4.081,  4.144,  4.208,  4.273,   
 49   4.339,  4.406,  4.473,  4.542,  4.612,  4.68     49   4.339,  4.406,  4.473,  4.542,  4.612,  4.683,  4.754,  4.827,  4.900,  4.975,
 50   5.000,  6.134,  7.269,  8.403,  9.538, 10.67     50   5.000,  6.134,  7.269,  8.403,  9.538, 10.672, 11.807, 12.941, 14.076, 15.210,  
 51  16.345, 17.479, 18.613, 19.748, 20.882, 22.01     51  16.345, 17.479, 18.613, 19.748, 20.882, 22.017, 23.151, 24.286, 25.420, 26.555, 
 52  27.689, 28.824, 29.958, 31.092, 32.227, 33.36     52  27.689, 28.824, 29.958, 31.092, 32.227, 33.361, 34.496, 35.630, 36.765, 37.899,  
 53  39.034, 40.168, 41.303, 42.437, 43.571, 44.70     53  39.034, 40.168, 41.303, 42.437, 43.571, 44.706, 45.840, 46.975, 48.109, 49.244
 54 };                                                 54 };
 55                                                    55 
 56 // Cross-sections in mb, from S.A. Bass et al.     56 // Cross-sections in mb, from S.A. Bass et al., Prog.Part.Nucl.Phys.41:225-370,1998 
 57 // Units are assigned when filling the Physics     57 // Units are assigned when filling the PhysicsVector
 58                                                    58 
 59 const G4double G4XNNstarTable::sigmaNN1440[121     59 const G4double G4XNNstarTable::sigmaNN1440[121] = 
 60 {                                                  60 {
 61   0.0,                                             61   0.0, 
 62   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.     62   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,  
 63   0.000, 0.000, 0.000, 0.001, 0.001, 0.002, 0.     63   0.000, 0.000, 0.000, 0.001, 0.001, 0.002, 0.003, 0.004, 0.006, 0.009,  
 64   0.014, 0.020, 0.031, 0.047, 0.072, 0.113, 0.     64   0.014, 0.020, 0.031, 0.047, 0.072, 0.113, 0.178, 0.266, 0.370, 0.476,  
 65   0.575, 0.665, 0.744, 0.814, 0.874, 0.926, 0.     65   0.575, 0.665, 0.744, 0.814, 0.874, 0.926, 0.971, 1.009, 1.040, 1.066,
 66   1.087, 1.103, 1.115, 1.124, 1.129, 1.132, 1.     66   1.087, 1.103, 1.115, 1.124, 1.129, 1.132, 1.131, 1.129, 1.124, 1.117, 
 67   1.109, 1.099, 1.088, 1.075, 1.062, 1.048, 1.     67   1.109, 1.099, 1.088, 1.075, 1.062, 1.048, 1.033, 1.017, 1.001, 0.984,  
 68   0.967, 0.950, 0.932, 0.914, 0.896, 0.879, 0.     68   0.967, 0.950, 0.932, 0.914, 0.896, 0.879, 0.861, 0.843, 0.825, 0.807,  
 69   0.790, 0.773, 0.755, 0.738, 0.722, 0.705, 0.     69   0.790, 0.773, 0.755, 0.738, 0.722, 0.705, 0.689, 0.673, 0.657, 0.641,
 70   0.636, 0.453, 0.336, 0.258, 0.204, 0.166, 0.     70   0.636, 0.453, 0.336, 0.258, 0.204, 0.166, 0.137, 0.115, 0.098, 0.084,  
 71   0.073, 0.064, 0.057, 0.051, 0.046, 0.041, 0.     71   0.073, 0.064, 0.057, 0.051, 0.046, 0.041, 0.037, 0.034, 0.031, 0.029,  
 72   0.026, 0.024, 0.022, 0.021, 0.019, 0.018, 0.     72   0.026, 0.024, 0.022, 0.021, 0.019, 0.018, 0.017, 0.016, 0.015, 0.014, 
 73   0.013, 0.013, 0.012, 0.011, 0.011, 0.010, 0.     73   0.013, 0.013, 0.012, 0.011, 0.011, 0.010, 0.010, 0.009, 0.009, 0.008
 74 };                                                 74 };
 75                                                    75 
 76 const G4double G4XNNstarTable::sigmaNN1520[121     76 const G4double G4XNNstarTable::sigmaNN1520[121] = 
 77 {                                                  77 {
 78   0.0,                                             78   0.0,
 79   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.     79   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 
 80   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.     80   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.001, 0.001, 
 81   0.001, 0.002, 0.003, 0.005, 0.008, 0.014, 0.     81   0.001, 0.002, 0.003, 0.005, 0.008, 0.014, 0.024, 0.043, 0.079, 0.150, 
 82   0.280, 0.465, 0.662, 0.841, 0.995, 1.125, 1.     82   0.280, 0.465, 0.662, 0.841, 0.995, 1.125, 1.235, 1.327, 1.410, 1.474,
 83   1.525, 1.567, 1.600, 1.625, 1.644, 1.651, 1.     83   1.525, 1.567, 1.600, 1.625, 1.644, 1.651, 1.657, 1.659, 1.657, 1.651, 
 84   1.642, 1.629, 1.614, 1.597, 1.578, 1.558, 1.     84   1.642, 1.629, 1.614, 1.597, 1.578, 1.558, 1.536, 1.512, 1.488, 1.463, 
 85   1.437, 1.411, 1.384, 1.357, 1.329, 1.302, 1.     85   1.437, 1.411, 1.384, 1.357, 1.329, 1.302, 1.274, 1.247, 1.219, 1.192, 
 86   1.165, 1.138, 1.112, 1.086, 1.060, 1.035, 1.     86   1.165, 1.138, 1.112, 1.086, 1.060, 1.035, 1.010, 0.985, 0.961, 0.938,
 87   0.930, 0.652, 0.479, 0.365, 0.287, 0.232, 0.     87   0.930, 0.652, 0.479, 0.365, 0.287, 0.232, 0.191, 0.160, 0.135, 0.116, 
 88   0.101, 0.089, 0.078, 0.070, 0.062, 0.056, 0.     88   0.101, 0.089, 0.078, 0.070, 0.062, 0.056, 0.051, 0.046, 0.042, 0.039, 
 89   0.036, 0.033, 0.031, 0.028, 0.026, 0.025, 0.     89   0.036, 0.033, 0.031, 0.028, 0.026, 0.025, 0.023, 0.022, 0.020, 0.019, 
 90   0.018, 0.017, 0.016, 0.015, 0.015, 0.014, 0.     90   0.018, 0.017, 0.016, 0.015, 0.015, 0.014, 0.013, 0.013, 0.012, 0.011
 91 };                                                 91 };
 92                                                    92 
 93 const G4double G4XNNstarTable::sigmaNN1535[121     93 const G4double G4XNNstarTable::sigmaNN1535[121] = 
 94 {                                                  94 {
 95   0.0,                                             95   0.0,
 96   0.000, 0.000, 0.001, 0.001, 0.001, 0.002, 0.     96   0.000, 0.000, 0.001, 0.001, 0.001, 0.002, 0.002, 0.003,
 97   0.004, 0.005, 0.006, 0.008, 0.010, 0.012, 0.     97   0.004, 0.005, 0.006, 0.008, 0.010, 0.012, 0.015, 0.019,
 98   0.024, 0.031, 0.039, 0.052, 0.069, 0.097, 0.     98   0.024, 0.031, 0.039, 0.052, 0.069, 0.097, 0.145, 0.216,
 99   0.298, 0.378, 0.451, 0.513, 0.566, 0.610, 0.     99   0.298, 0.378, 0.451, 0.513, 0.566, 0.610, 0.646, 0.675,
100   0.699, 0.718, 0.732, 0.742, 0.749, 0.753, 0.    100   0.699, 0.718, 0.732, 0.742, 0.749, 0.753, 0.754, 0.753,
101   0.751, 0.746, 0.740, 0.733, 0.724, 0.715, 0.    101   0.751, 0.746, 0.740, 0.733, 0.724, 0.715, 0.705, 0.694,
102   0.683, 0.671, 0.659, 0.647, 0.634, 0.621, 0.    102   0.683, 0.671, 0.659, 0.647, 0.634, 0.621, 0.608, 0.595,
103   0.582, 0.569, 0.556, 0.543, 0.531, 0.518, 0.    103   0.582, 0.569, 0.556, 0.543, 0.531, 0.518, 0.506, 0.493,
104   0.481, 0.469, 0.458, 0.446, 0.435, 0.424, 0.    104   0.481, 0.469, 0.458, 0.446, 0.435, 0.424, 0.413, 0.402,
105   0.399, 0.276, 0.201, 0.153, 0.120, 0.096, 0.    105   0.399, 0.276, 0.201, 0.153, 0.120, 0.096, 0.079, 0.066,
106   0.056, 0.048, 0.042, 0.037, 0.032, 0.029, 0.    106   0.056, 0.048, 0.042, 0.037, 0.032, 0.029, 0.026, 0.023,
107   0.021, 0.019, 0.017, 0.016, 0.015, 0.014, 0.    107   0.021, 0.019, 0.017, 0.016, 0.015, 0.014, 0.013, 0.012,
108   0.011, 0.010, 0.009, 0.009, 0.008, 0.008, 0.    108   0.011, 0.010, 0.009, 0.009, 0.008, 0.008, 0.007, 0.007,
109   0.007, 0.006, 0.006, 0.006, 0.005, 0.005, 0.    109   0.007, 0.006, 0.006, 0.006, 0.005, 0.005, 0.005, 0.005
110 };                                                110 };
111                                                   111 
112 const G4double G4XNNstarTable::sigmaNN1650[121    112 const G4double G4XNNstarTable::sigmaNN1650[121] = 
113 {                                                 113 {
114   0.0,                                            114   0.0,
115   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    115   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
116   0.000, 0.000, 0.000, 0.000, 0.001, 0.001, 0.    116   0.000, 0.000, 0.000, 0.000, 0.001, 0.001, 0.001, 0.001,
117   0.002, 0.002, 0.002, 0.003, 0.004, 0.005, 0.    117   0.002, 0.002, 0.002, 0.003, 0.004, 0.005, 0.006, 0.007,
118   0.008, 0.010, 0.012, 0.015, 0.019, 0.023, 0.    118   0.008, 0.010, 0.012, 0.015, 0.019, 0.023, 0.029, 0.038,
119   0.051, 0.071, 0.104, 0.150, 0.201, 0.249, 0.    119   0.051, 0.071, 0.104, 0.150, 0.201, 0.249, 0.290, 0.326,
120   0.354, 0.378, 0.397, 0.412, 0.424, 0.434, 0.    120   0.354, 0.378, 0.397, 0.412, 0.424, 0.434, 0.440, 0.445,
121   0.448, 0.449, 0.449, 0.448, 0.445, 0.442, 0.    121   0.448, 0.449, 0.449, 0.448, 0.445, 0.442, 0.438, 0.433,
122   0.428, 0.422, 0.416, 0.409, 0.403, 0.395, 0.    122   0.428, 0.422, 0.416, 0.409, 0.403, 0.395, 0.388, 0.381,
123   0.373, 0.366, 0.358, 0.350, 0.343, 0.335, 0.    123   0.373, 0.366, 0.358, 0.350, 0.343, 0.335, 0.327, 0.320,
124   0.312, 0.305, 0.298, 0.291, 0.284, 0.277, 0.    124   0.312, 0.305, 0.298, 0.291, 0.284, 0.277, 0.270, 0.263,
125   0.261, 0.182, 0.133, 0.102, 0.080, 0.064, 0.    125   0.261, 0.182, 0.133, 0.102, 0.080, 0.064, 0.053, 0.044,
126   0.037, 0.032, 0.028, 0.024, 0.022, 0.019, 0.    126   0.037, 0.032, 0.028, 0.024, 0.022, 0.019, 0.017, 0.015,
127   0.014, 0.013, 0.012, 0.011, 0.010, 0.009, 0.    127   0.014, 0.013, 0.012, 0.011, 0.010, 0.009, 0.008, 0.008,
128   0.007, 0.007, 0.006, 0.006, 0.006, 0.005, 0.    128   0.007, 0.007, 0.006, 0.006, 0.006, 0.005, 0.005, 0.005,
129   0.004, 0.004, 0.004, 0.004, 0.004, 0.003, 0.    129   0.004, 0.004, 0.004, 0.004, 0.004, 0.003, 0.003, 0.003
130 };                                                130 };
131                                                   131 
132 const G4double G4XNNstarTable::sigmaNN1675[121    132 const G4double G4XNNstarTable::sigmaNN1675[121] = 
133 {                                                 133 {
134   0.0,                                            134   0.0,
135   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    135   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
136   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    136   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
137   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    137   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.001,
138   0.001, 0.001, 0.002, 0.002, 0.003, 0.005, 0.    138   0.001, 0.001, 0.002, 0.002, 0.003, 0.005, 0.009, 0.015,
139   0.026, 0.048, 0.095, 0.189, 0.324, 0.463, 0.    139   0.026, 0.048, 0.095, 0.189, 0.324, 0.463, 0.586, 0.691,
140   0.780, 0.855, 0.919, 0.972, 1.016, 1.052, 1.    140   0.780, 0.855, 0.919, 0.972, 1.016, 1.052, 1.081, 1.105,
141   1.123, 1.136, 1.145, 1.151, 1.153, 1.153, 1.    141   1.123, 1.136, 1.145, 1.151, 1.153, 1.153, 1.149, 1.144,
142   1.136, 1.127, 1.116, 1.104, 1.090, 1.076, 1.    142   1.136, 1.127, 1.116, 1.104, 1.090, 1.076, 1.061, 1.045,
143   1.028, 1.011, 0.993, 0.975, 0.957, 0.939, 0.    143   1.028, 1.011, 0.993, 0.975, 0.957, 0.939, 0.921, 0.902,
144   0.884, 0.865, 0.847, 0.828, 0.810, 0.792, 0.    144   0.884, 0.865, 0.847, 0.828, 0.810, 0.792, 0.775, 0.757,
145   0.751, 0.538, 0.399, 0.307, 0.242, 0.196, 0.    145   0.751, 0.538, 0.399, 0.307, 0.242, 0.196, 0.162, 0.136,
146   0.115, 0.099, 0.086, 0.076, 0.067, 0.060, 0.    146   0.115, 0.099, 0.086, 0.076, 0.067, 0.060, 0.053, 0.048,
147   0.044, 0.040, 0.036, 0.033, 0.031, 0.028, 0.    147   0.044, 0.040, 0.036, 0.033, 0.031, 0.028, 0.026, 0.024,
148   0.023, 0.021, 0.020, 0.019, 0.018, 0.016, 0.    148   0.023, 0.021, 0.020, 0.019, 0.018, 0.016, 0.016, 0.015,
149   0.014, 0.013, 0.013, 0.012, 0.011, 0.011, 0.    149   0.014, 0.013, 0.013, 0.012, 0.011, 0.011, 0.010, 0.010 
150 };                                                150 };
151                                                   151 
152 const G4double G4XNNstarTable::sigmaNN1680[121    152 const G4double G4XNNstarTable::sigmaNN1680[121] = 
153 {                                                 153 {
154   0.0,                                            154   0.0,
155   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    155   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
156   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    156   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
157   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    157   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
158   0.000, 0.001, 0.001, 0.001, 0.002, 0.003, 0.    158   0.000, 0.001, 0.001, 0.001, 0.002, 0.003, 0.006, 0.010,
159   0.018, 0.035, 0.073, 0.156, 0.294, 0.446, 0.    159   0.018, 0.035, 0.073, 0.156, 0.294, 0.446, 0.580, 0.693,
160   0.788, 0.867, 0.933, 0.988, 1.033, 1.070, 1.    160   0.788, 0.867, 0.933, 0.988, 1.033, 1.070, 1.100, 1.124,
161   1.142, 1.155, 1.163, 1.168, 1.170, 1.168, 1.    161   1.142, 1.155, 1.163, 1.168, 1.170, 1.168, 1.164, 1.158,
162   1.149, 1.139, 1.127, 1.114, 1.100, 1.085, 1.    162   1.149, 1.139, 1.127, 1.114, 1.100, 1.085, 1.068, 1.052,
163   1.034, 1.016, 0.998, 0.979, 0.960, 0.941, 0.    163   1.034, 1.016, 0.998, 0.979, 0.960, 0.941, 0.922, 0.903,
164   0.884, 0.865, 0.846, 0.827, 0.809, 0.790, 0.    164   0.884, 0.865, 0.846, 0.827, 0.809, 0.790, 0.772, 0.754,
165   0.748, 0.533, 0.394, 0.301, 0.238, 0.192, 0.    165   0.748, 0.533, 0.394, 0.301, 0.238, 0.192, 0.158, 0.133,
166   0.113, 0.097, 0.084, 0.074, 0.065, 0.058, 0.    166   0.113, 0.097, 0.084, 0.074, 0.065, 0.058, 0.052, 0.047,
167   0.042, 0.039, 0.035, 0.032, 0.030, 0.028, 0.    167   0.042, 0.039, 0.035, 0.032, 0.030, 0.028, 0.026, 0.024,
168   0.022, 0.021, 0.019, 0.018, 0.017, 0.016, 0.    168   0.022, 0.021, 0.019, 0.018, 0.017, 0.016, 0.015, 0.014,
169   0.013, 0.013, 0.012, 0.012, 0.011, 0.010, 0.    169   0.013, 0.013, 0.012, 0.012, 0.011, 0.010, 0.010, 0.010
170 };                                                170 };
171                                                   171 
172 const G4double G4XNNstarTable::sigmaNN1700[121    172 const G4double G4XNNstarTable::sigmaNN1700[121] = 
173 {                                                 173 {
174   0.0,                                            174   0.0,
175   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    175   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
176   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    176   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
177   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    177   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
178   0.000, 0.001, 0.001, 0.001, 0.002, 0.003, 0.    178   0.000, 0.001, 0.001, 0.001, 0.002, 0.003, 0.005, 0.008,
179   0.013, 0.022, 0.038, 0.070, 0.140, 0.245, 0.    179   0.013, 0.022, 0.038, 0.070, 0.140, 0.245, 0.347, 0.431,
180   0.500, 0.556, 0.601, 0.637, 0.666, 0.689, 0.    180   0.500, 0.556, 0.601, 0.637, 0.666, 0.689, 0.708, 0.721,
181   0.731, 0.738, 0.742, 0.743, 0.743, 0.740, 0.    181   0.731, 0.738, 0.742, 0.743, 0.743, 0.740, 0.736, 0.730,
182   0.723, 0.716, 0.707, 0.697, 0.687, 0.676, 0.    182   0.723, 0.716, 0.707, 0.697, 0.687, 0.676, 0.665, 0.654,
183   0.642, 0.630, 0.617, 0.605, 0.593, 0.580, 0.    183   0.642, 0.630, 0.617, 0.605, 0.593, 0.580, 0.568, 0.555,
184   0.543, 0.531, 0.519, 0.507, 0.495, 0.483, 0.    184   0.543, 0.531, 0.519, 0.507, 0.495, 0.483, 0.471, 0.460,
185   0.456, 0.322, 0.236, 0.180, 0.142, 0.114, 0.    185   0.456, 0.322, 0.236, 0.180, 0.142, 0.114, 0.094, 0.079,
186   0.067, 0.057, 0.050, 0.044, 0.039, 0.034, 0.    186   0.067, 0.057, 0.050, 0.044, 0.039, 0.034, 0.031, 0.028,
187   0.025, 0.023, 0.021, 0.019, 0.018, 0.016, 0.    187   0.025, 0.023, 0.021, 0.019, 0.018, 0.016, 0.015, 0.014,
188   0.013, 0.012, 0.011, 0.011, 0.010, 0.009, 0.    188   0.013, 0.012, 0.011, 0.011, 0.010, 0.009, 0.009, 0.008,
189   0.008, 0.008, 0.007, 0.007, 0.007, 0.006, 0.    189   0.008, 0.008, 0.007, 0.007, 0.007, 0.006, 0.006, 0.006
190 };                                                190 };
191                                                   191 
192 const G4double G4XNNstarTable::sigmaNN1710[121    192 const G4double G4XNNstarTable::sigmaNN1710[121] = 
193 {                                                 193 {
194   0.0,                                            194   0.0,
195   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    195   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
196   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    196   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
197   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    197   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
198   0.000, 0.000, 0.000, 0.001, 0.001, 0.001, 0.    198   0.000, 0.000, 0.000, 0.001, 0.001, 0.001, 0.002, 0.003,
199   0.005, 0.008, 0.013, 0.025, 0.052, 0.096, 0.    199   0.005, 0.008, 0.013, 0.025, 0.052, 0.096, 0.144, 0.185,
200   0.219, 0.246, 0.269, 0.288, 0.304, 0.316, 0.    200   0.219, 0.246, 0.269, 0.288, 0.304, 0.316, 0.326, 0.334,
201   0.340, 0.344, 0.347, 0.349, 0.349, 0.349, 0.    201   0.340, 0.344, 0.347, 0.349, 0.349, 0.349, 0.348, 0.346,
202   0.343, 0.340, 0.336, 0.332, 0.328, 0.323, 0.    202   0.343, 0.340, 0.336, 0.332, 0.328, 0.323, 0.319, 0.313,
203   0.308, 0.303, 0.297, 0.292, 0.286, 0.280, 0.    203   0.308, 0.303, 0.297, 0.292, 0.286, 0.280, 0.275, 0.269,
204   0.263, 0.257, 0.252, 0.246, 0.241, 0.235, 0.    204   0.263, 0.257, 0.252, 0.246, 0.241, 0.235, 0.230, 0.224,
205   0.223, 0.158, 0.117, 0.090, 0.071, 0.057, 0.    205   0.223, 0.158, 0.117, 0.090, 0.071, 0.057, 0.047, 0.040,
206   0.034, 0.029, 0.025, 0.022, 0.019, 0.017, 0.    206   0.034, 0.029, 0.025, 0.022, 0.019, 0.017, 0.016, 0.014,
207   0.013, 0.012, 0.011, 0.010, 0.009, 0.008, 0.    207   0.013, 0.012, 0.011, 0.010, 0.009, 0.008, 0.008, 0.007,
208   0.007, 0.006, 0.006, 0.005, 0.005, 0.005, 0.    208   0.007, 0.006, 0.006, 0.005, 0.005, 0.005, 0.005, 0.004,
209   0.004, 0.004, 0.004, 0.003, 0.003, 0.003, 0.    209   0.004, 0.004, 0.004, 0.003, 0.003, 0.003, 0.003, 0.003
210 };                                                210 };
211                                                   211 
212 const G4double G4XNNstarTable::sigmaNN1720[121    212 const G4double G4XNNstarTable::sigmaNN1720[121] = 
213 {                                                 213 {
214   0.0,                                            214   0.0,
215   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    215   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
216   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    216   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
217   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    217   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
218   0.001, 0.001, 0.002, 0.002, 0.003, 0.005, 0.    218   0.001, 0.001, 0.002, 0.002, 0.003, 0.005, 0.007, 0.010,
219   0.015, 0.023, 0.036, 0.061, 0.106, 0.174, 0.    219   0.015, 0.023, 0.036, 0.061, 0.106, 0.174, 0.250, 0.321,
220   0.382, 0.434, 0.477, 0.513, 0.543, 0.568, 0.    220   0.382, 0.434, 0.477, 0.513, 0.543, 0.568, 0.589, 0.605,
221   0.617, 0.627, 0.634, 0.639, 0.641, 0.642, 0.    221   0.617, 0.627, 0.634, 0.639, 0.641, 0.642, 0.641, 0.639,
222   0.635, 0.630, 0.625, 0.618, 0.611, 0.603, 0.    222   0.635, 0.630, 0.625, 0.618, 0.611, 0.603, 0.595, 0.586,
223   0.576, 0.567, 0.557, 0.547, 0.537, 0.527, 0.    223   0.576, 0.567, 0.557, 0.547, 0.537, 0.527, 0.516, 0.506,
224   0.496, 0.485, 0.475, 0.465, 0.455, 0.444, 0.    224   0.496, 0.485, 0.475, 0.465, 0.455, 0.444, 0.435, 0.425,
225   0.421, 0.302, 0.224, 0.172, 0.136, 0.117, 0.    225   0.421, 0.302, 0.224, 0.172, 0.136, 0.117, 0.091, 0.076,
226   0.065, 0.056, 0.049, 0.043, 0.038, 0.034, 0.    226   0.065, 0.056, 0.049, 0.043, 0.038, 0.034, 0.030, 0.027,
227   0.025, 0.022, 0.020, 0.019, 0.017, 0.016, 0.    227   0.025, 0.022, 0.020, 0.019, 0.017, 0.016, 0.015, 0.014,
228   0.013, 0.012, 0.011, 0.010, 0.010, 0.009, 0.    228   0.013, 0.012, 0.011, 0.010, 0.010, 0.009, 0.009, 0.008,
229   0.008, 0.007, 0.007, 0.007, 0.006, 0.006, 0.    229   0.008, 0.007, 0.007, 0.007, 0.006, 0.006, 0.006, 0.006 
230 };                                                230 };
231                                                   231 
232 const G4double G4XNNstarTable::sigmaNN1900[121    232 const G4double G4XNNstarTable::sigmaNN1900[121] = 
233 {                                                 233 {
234   0.0,                                            234   0.0,
235   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    235   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
236   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    236   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
237   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    237   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.001,
238   0.001, 0.001, 0.001, 0.002, 0.002, 0.003, 0.    238   0.001, 0.001, 0.001, 0.002, 0.002, 0.003, 0.003, 0.005,
239   0.006, 0.008, 0.010, 0.014, 0.019, 0.026, 0.    239   0.006, 0.008, 0.010, 0.014, 0.019, 0.026, 0.037, 0.054,
240   0.074, 0.094, 0.114, 0.131, 0.147, 0.161, 0.    240   0.074, 0.094, 0.114, 0.131, 0.147, 0.161, 0.173, 0.184,
241   0.193, 0.201, 0.208, 0.213, 0.218, 0.221, 0.    241   0.193, 0.201, 0.208, 0.213, 0.218, 0.221, 0.224, 0.226,
242   0.228, 0.229, 0.229, 0.229, 0.229, 0.228, 0.    242   0.228, 0.229, 0.229, 0.229, 0.229, 0.228, 0.227, 0.225,
243   0.223, 0.221, 0.219, 0.217, 0.214, 0.212, 0.    243   0.223, 0.221, 0.219, 0.217, 0.214, 0.212, 0.209, 0.206,
244   0.203, 0.200, 0.197, 0.194, 0.190, 0.187, 0.    244   0.203, 0.200, 0.197, 0.194, 0.190, 0.187, 0.184, 0.181,
245   0.180, 0.137, 0.106, 0.083, 0.067, 0.056, 0.    245   0.180, 0.137, 0.106, 0.083, 0.067, 0.056, 0.047, 0.039,
246   0.034, 0.029, 0.026, 0.023, 0.020, 0.018, 0.    246   0.034, 0.029, 0.026, 0.023, 0.020, 0.018, 0.016, 0.015,
247   0.013, 0.012, 0.011, 0.010, 0.010, 0.009, 0.    247   0.013, 0.012, 0.011, 0.010, 0.010, 0.009, 0.008, 0.008,
248   0.007, 0.007, 0.006, 0.006, 0.006, 0.005, 0.    248   0.007, 0.007, 0.006, 0.006, 0.006, 0.005, 0.005, 0.005,
249   0.004, 0.004, 0.004, 0.004, 0.004, 0.003, 0.    249   0.004, 0.004, 0.004, 0.004, 0.004, 0.003, 0.003, 0.003 
250 };                                                250 };
251                                                   251 
252 const G4double G4XNNstarTable::sigmaNN1990[121    252 const G4double G4XNNstarTable::sigmaNN1990[121] = 
253 {                                                 253 {
254   0.0,                                            254   0.0,
255   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    255   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
256   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    256   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
257   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    257   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
258   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    258   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
259   0.000, 0.000, 0.001, 0.001, 0.002, 0.003, 0.    259   0.000, 0.000, 0.001, 0.001, 0.002, 0.003, 0.005, 0.009,
260   0.017, 0.030, 0.051, 0.076, 0.102, 0.127, 0.    260   0.017, 0.030, 0.051, 0.076, 0.102, 0.127, 0.149, 0.168,
261   0.185, 0.199, 0.211, 0.221, 0.230, 0.237, 0.    261   0.185, 0.199, 0.211, 0.221, 0.230, 0.237, 0.243, 0.247,
262   0.251, 0.254, 0.256, 0.258, 0.258, 0.259, 0.    262   0.251, 0.254, 0.256, 0.258, 0.258, 0.259, 0.258, 0.258,
263   0.257, 0.255, 0.254, 0.252, 0.249, 0.247, 0.    263   0.257, 0.255, 0.254, 0.252, 0.249, 0.247, 0.245, 0.242,
264   0.239, 0.236, 0.233, 0.230, 0.226, 0.223, 0.    264   0.239, 0.236, 0.233, 0.230, 0.226, 0.223, 0.220, 0.216,
265   0.215, 0.167, 0.131, 0.104, 0.085, 0.070, 0.    265   0.215, 0.167, 0.131, 0.104, 0.085, 0.070, 0.059, 0.050,
266   0.043, 0.038, 0.033, 0.029, 0.026, 0.023, 0.    266   0.043, 0.038, 0.033, 0.029, 0.026, 0.023, 0.021, 0.019,
267   0.017, 0.016, 0.014, 0.013, 0.012, 0.011, 0.    267   0.017, 0.016, 0.014, 0.013, 0.012, 0.011, 0.010, 0.010,
268   0.009, 0.009, 0.008, 0.007, 0.007, 0.007, 0.    268   0.009, 0.009, 0.008, 0.007, 0.007, 0.007, 0.006, 0.006,
269   0.006, 0.005, 0.005, 0.005, 0.005, 0.004, 0.    269   0.006, 0.005, 0.005, 0.005, 0.005, 0.004, 0.004, 0.004
270 };                                                270 };
271                                                   271 
272 const G4double G4XNNstarTable::sigmaNN2090[121    272 const G4double G4XNNstarTable::sigmaNN2090[121] = 
273 {                                                 273 {
274   0.0,                                            274   0.0,
275   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    275   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
276   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    276   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
277   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    277   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
278   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    278   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
279   0.000, 0.001, 0.001, 0.001, 0.002, 0.002, 0.    279   0.000, 0.001, 0.001, 0.001, 0.002, 0.002, 0.003, 0.005,
280   0.007, 0.011, 0.016, 0.024, 0.036, 0.053, 0.    280   0.007, 0.011, 0.016, 0.024, 0.036, 0.053, 0.071, 0.089,
281   0.106, 0.120, 0.133, 0.142, 0.151, 0.158, 0.    281   0.106, 0.120, 0.133, 0.142, 0.151, 0.158, 0.164, 0.169,
282   0.172, 0.175, 0.178, 0.179, 0.180, 0.180, 0.    282   0.172, 0.175, 0.178, 0.179, 0.180, 0.180, 0.180, 0.180,
283   0.179, 0.178, 0.176, 0.175, 0.173, 0.171, 0.    283   0.179, 0.178, 0.176, 0.175, 0.173, 0.171, 0.169, 0.166,
284   0.164, 0.162, 0.159, 0.156, 0.154, 0.151, 0.    284   0.164, 0.162, 0.159, 0.156, 0.154, 0.151, 0.148, 0.145,
285   0.144, 0.107, 0.081, 0.063, 0.050, 0.041, 0.    285   0.144, 0.107, 0.081, 0.063, 0.050, 0.041, 0.034, 0.028,
286   0.024, 0.021, 0.018, 0.016, 0.014, 0.013, 0.    286   0.024, 0.021, 0.018, 0.016, 0.014, 0.013, 0.011, 0.010,
287   0.009, 0.008, 0.008, 0.007, 0.006, 0.006, 0.    287   0.009, 0.008, 0.008, 0.007, 0.006, 0.006, 0.006, 0.005,
288   0.005, 0.004, 0.004, 0.004, 0.004, 0.003, 0.    288   0.005, 0.004, 0.004, 0.004, 0.004, 0.003, 0.003, 0.003,
289   0.003, 0.003, 0.003, 0.003, 0.002, 0.002, 0.    289   0.003, 0.003, 0.003, 0.003, 0.002, 0.002, 0.002, 0.002
290 };                                                290 };
291                                                   291 
292 const G4double G4XNNstarTable::sigmaNN2190[121    292 const G4double G4XNNstarTable::sigmaNN2190[121] = 
293 {                                                 293 {
294   0.0,                                            294   0.0,
295   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    295   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
296   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    296   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
297   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    297   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
298   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    298   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.001,
299   0.001, 0.002, 0.003, 0.005, 0.009, 0.015, 0.    299   0.001, 0.002, 0.003, 0.005, 0.009, 0.015, 0.024, 0.036,
300   0.050, 0.064, 0.078, 0.090, 0.100, 0.110, 0.    300   0.050, 0.064, 0.078, 0.090, 0.100, 0.110, 0.118, 0.125,
301   0.131, 0.136, 0.140, 0.143, 0.146, 0.149, 0.    301   0.131, 0.136, 0.140, 0.143, 0.146, 0.149, 0.150, 0.152,
302   0.153, 0.153, 0.153, 0.153, 0.153, 0.152, 0.    302   0.153, 0.153, 0.153, 0.153, 0.153, 0.152, 0.152, 0.151,
303   0.150, 0.148, 0.147, 0.145, 0.144, 0.142, 0.    303   0.150, 0.148, 0.147, 0.145, 0.144, 0.142, 0.140, 0.139,
304   0.138, 0.110, 0.087, 0.069, 0.056, 0.047, 0.    304   0.138, 0.110, 0.087, 0.069, 0.056, 0.047, 0.039, 0.033,
305   0.029, 0.025, 0.022, 0.019, 0.017, 0.015, 0.    305   0.029, 0.025, 0.022, 0.019, 0.017, 0.015, 0.014, 0.013,
306   0.012, 0.010, 0.009, 0.009, 0.008, 0.007, 0.    306   0.012, 0.010, 0.009, 0.009, 0.008, 0.007, 0.007, 0.006,
307   0.006, 0.006, 0.005, 0.005, 0.005, 0.004, 0.    307   0.006, 0.006, 0.005, 0.005, 0.005, 0.004, 0.004, 0.004,
308   0.004, 0.003, 0.003, 0.003, 0.003, 0.003, 0.    308   0.004, 0.003, 0.003, 0.003, 0.003, 0.003, 0.003, 0.003
309 };                                                309 };
310                                                   310 
311 const G4double G4XNNstarTable::sigmaNN2220[121    311 const G4double G4XNNstarTable::sigmaNN2220[121] = 
312 {                                                 312 {
313   0.0,                                            313   0.0,
314   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    314   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
315   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    315   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
316   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    316   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
317   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    317   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
318   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    318   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
319   0.001, 0.001, 0.001, 0.003, 0.005, 0.008, 0.    319   0.001, 0.001, 0.001, 0.003, 0.005, 0.008, 0.015, 0.026,
320   0.041, 0.058, 0.074, 0.089, 0.101, 0.113, 0.    320   0.041, 0.058, 0.074, 0.089, 0.101, 0.113, 0.122, 0.131,
321   0.138, 0.144, 0.149, 0.153, 0.156, 0.159, 0.    321   0.138, 0.144, 0.149, 0.153, 0.156, 0.159, 0.161, 0.163,
322   0.164, 0.165, 0.165, 0.165, 0.165, 0.165, 0.    322   0.164, 0.165, 0.165, 0.165, 0.165, 0.165, 0.164, 0.163,
323   0.162, 0.161, 0.160, 0.158, 0.157, 0.155, 0.    323   0.162, 0.161, 0.160, 0.158, 0.157, 0.155, 0.153, 0.151,
324   0.150, 0.121, 0.096, 0.077, 0.062, 0.052, 0.    324   0.150, 0.121, 0.096, 0.077, 0.062, 0.052, 0.043, 0.037,
325   0.032, 0.028, 0.024, 0.022, 0.019, 0.017, 0.    325   0.032, 0.028, 0.024, 0.022, 0.019, 0.017, 0.015, 0.014,
326   0.013, 0.011, 0.011, 0.010, 0.009, 0.008, 0.    326   0.013, 0.011, 0.011, 0.010, 0.009, 0.008, 0.008, 0.007,
327   0.007, 0.006, 0.006, 0.005, 0.005, 0.005, 0.    327   0.007, 0.006, 0.006, 0.005, 0.005, 0.005, 0.005, 0.004,
328   0.004, 0.004, 0.004, 0.004, 0.003, 0.003, 0.    328   0.004, 0.004, 0.004, 0.004, 0.003, 0.003, 0.003, 0.003 
329 };                                                329 };
330                                                   330 
331 const G4double G4XNNstarTable::sigmaNN2250[121    331 const G4double G4XNNstarTable::sigmaNN2250[121] = 
332 {                                                 332 {
333   0.0,                                            333   0.0,
334   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    334   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
335   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    335   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
336   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    336   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
337   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    337   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
338   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.    338   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
339   0.000, 0.001, 0.001, 0.002, 0.003, 0.005, 0.    339   0.000, 0.001, 0.001, 0.002, 0.003, 0.005, 0.009, 0.016,
340   0.027, 0.043, 0.060, 0.076, 0.090, 0.103, 0.    340   0.027, 0.043, 0.060, 0.076, 0.090, 0.103, 0.114, 0.123,
341   0.131, 0.138, 0.144, 0.149, 0.153, 0.156, 0.    341   0.131, 0.138, 0.144, 0.149, 0.153, 0.156, 0.159, 0.161,
342   0.163, 0.164, 0.165, 0.165, 0.165, 0.165, 0.    342   0.163, 0.164, 0.165, 0.165, 0.165, 0.165, 0.165, 0.164,
343   0.163, 0.162, 0.161, 0.159, 0.158, 0.156, 0.    343   0.163, 0.162, 0.161, 0.159, 0.158, 0.156, 0.154, 0.152,
344   0.152, 0.122, 0.096, 0.077, 0.062, 0.051, 0.    344   0.152, 0.122, 0.096, 0.077, 0.062, 0.051, 0.043, 0.037,
345   0.031, 0.027, 0.024, 0.022, 0.019, 0.017, 0.    345   0.031, 0.027, 0.024, 0.022, 0.019, 0.017, 0.015, 0.014,
346   0.012, 0.011, 0.010, 0.009, 0.009, 0.008, 0.    346   0.012, 0.011, 0.010, 0.009, 0.009, 0.008, 0.008, 0.007,
347   0.007, 0.006, 0.006, 0.005, 0.005, 0.005, 0.    347   0.007, 0.006, 0.006, 0.005, 0.005, 0.005, 0.004, 0.004,
348   0.004, 0.004, 0.004, 0.003, 0.003, 0.003, 0.    348   0.004, 0.004, 0.004, 0.003, 0.003, 0.003, 0.003, 0.003
349 };                                                349 };
350                                                   350 
351                                                   351 
352 G4XNNstarTable::G4XNNstarTable()                  352 G4XNNstarTable::G4XNNstarTable() 
353 {                                                 353 {
354   xMap["N(1440)0"] = (G4double*) sigmaNN1440;     354   xMap["N(1440)0"] = (G4double*) sigmaNN1440;
355   xMap["N(1440)+"] = (G4double*) sigmaNN1440;     355   xMap["N(1440)+"] = (G4double*) sigmaNN1440;
356                                                   356   
357   xMap["N(1520)0"] = (G4double*) sigmaNN1520;     357   xMap["N(1520)0"] = (G4double*) sigmaNN1520;
358   xMap["N(1520)+"] = (G4double*) sigmaNN1520;     358   xMap["N(1520)+"] = (G4double*) sigmaNN1520;
359                                                   359   
360   xMap["N(1535)0"] = (G4double*) sigmaNN1535;     360   xMap["N(1535)0"] = (G4double*) sigmaNN1535;
361   xMap["N(1535)+"] = (G4double*) sigmaNN1535;     361   xMap["N(1535)+"] = (G4double*) sigmaNN1535;
362                                                   362   
363   xMap["N(1650)0"] = (G4double*) sigmaNN1650;     363   xMap["N(1650)0"] = (G4double*) sigmaNN1650;
364   xMap["N(1650)+"] = (G4double*) sigmaNN1650;     364   xMap["N(1650)+"] = (G4double*) sigmaNN1650;
365                                                   365   
366   xMap["N(1675)0"] = (G4double*) sigmaNN1675;     366   xMap["N(1675)0"] = (G4double*) sigmaNN1675;
367   xMap["N(1675)+"] = (G4double*) sigmaNN1675;     367   xMap["N(1675)+"] = (G4double*) sigmaNN1675;
368                                                   368   
369   xMap["N(1680)0"] = (G4double*) sigmaNN1680;     369   xMap["N(1680)0"] = (G4double*) sigmaNN1680;
370   xMap["N(1680)+"] = (G4double*) sigmaNN1680;     370   xMap["N(1680)+"] = (G4double*) sigmaNN1680;
371                                                   371   
372   xMap["N(1700)0"] = (G4double*) sigmaNN1700;     372   xMap["N(1700)0"] = (G4double*) sigmaNN1700;
373   xMap["N(1700)+"] = (G4double*) sigmaNN1700;     373   xMap["N(1700)+"] = (G4double*) sigmaNN1700;
374                                                   374   
375   xMap["N(1710)0"] = (G4double*) sigmaNN1710;     375   xMap["N(1710)0"] = (G4double*) sigmaNN1710;
376   xMap["N(1710)+"] = (G4double*) sigmaNN1710;     376   xMap["N(1710)+"] = (G4double*) sigmaNN1710;
377                                                   377   
378   xMap["N(1720)0"] = (G4double*) sigmaNN1720;     378   xMap["N(1720)0"] = (G4double*) sigmaNN1720;
379   xMap["N(1720)+"] = (G4double*) sigmaNN1720;     379   xMap["N(1720)+"] = (G4double*) sigmaNN1720;
380                                                   380   
381   xMap["N(1900)0"] = (G4double*) sigmaNN1900;     381   xMap["N(1900)0"] = (G4double*) sigmaNN1900;
382   xMap["N(1900)+"] = (G4double*) sigmaNN1900;     382   xMap["N(1900)+"] = (G4double*) sigmaNN1900;
383                                                   383     
384   xMap["N(1990)0"] = (G4double*) sigmaNN1990;     384   xMap["N(1990)0"] = (G4double*) sigmaNN1990;
385   xMap["N(1990)+"] = (G4double*) sigmaNN1990;     385   xMap["N(1990)+"] = (G4double*) sigmaNN1990;
386                                                   386   
387   xMap["N(2090)0"] = (G4double*) sigmaNN2090;     387   xMap["N(2090)0"] = (G4double*) sigmaNN2090;
388   xMap["N(2090)+"] = (G4double*) sigmaNN2090;     388   xMap["N(2090)+"] = (G4double*) sigmaNN2090;
389                                                   389   
390   xMap["N(2190)0"] = (G4double*) sigmaNN2190;     390   xMap["N(2190)0"] = (G4double*) sigmaNN2190;
391   xMap["N(2190)+"] = (G4double*) sigmaNN2190;     391   xMap["N(2190)+"] = (G4double*) sigmaNN2190;
392                                                   392   
393   xMap["N(2220)0"] = (G4double*) sigmaNN2220;     393   xMap["N(2220)0"] = (G4double*) sigmaNN2220;
394   xMap["N(2220)+"] = (G4double*) sigmaNN2220;     394   xMap["N(2220)+"] = (G4double*) sigmaNN2220;
395                                                   395   
396   xMap["N(2250)0"] = (G4double*) sigmaNN2250;     396   xMap["N(2250)0"] = (G4double*) sigmaNN2250;
397   xMap["N(2250)+"] = (G4double*) sigmaNN2250;     397   xMap["N(2250)+"] = (G4double*) sigmaNN2250;  
398 }                                                 398 }
399                                                   399 
400                                                   400 
401 G4XNNstarTable::~G4XNNstarTable()                 401 G4XNNstarTable::~G4XNNstarTable()
402 { }                                               402 { }
403                                                   403 
404                                                   404 
405 const G4PhysicsVector* G4XNNstarTable::CrossSe    405 const G4PhysicsVector* G4XNNstarTable::CrossSectionTable(const G4String& particleName) const 
406 {                                                 406 {
407   // NOTE: the returned pointer is owned by th    407   // NOTE: the returned pointer is owned by the client
408                                                   408 
409   if (xMap.find(particleName) != xMap.end())      409   if (xMap.find(particleName) != xMap.end())
410     {                                             410     {
411       // Cross section table for the requested    411       // Cross section table for the requested particle available in the Map
412       G4PhysicsFreeVector* sigmaVector = new G    412       G4PhysicsFreeVector* sigmaVector = new G4PhysicsFreeVector(sizeNNstar);
413       std::map <G4String, G4double*, std::less    413       std::map <G4String, G4double*, std::less<G4String> >::const_iterator iter;
414       G4double* sigmaPointer = 0;                 414       G4double* sigmaPointer = 0;
415       for (iter = xMap.begin(); iter != xMap.e    415       for (iter = xMap.begin(); iter != xMap.end(); ++iter)
416   {                                               416   {
417     G4String str = (*iter).first;                 417     G4String str = (*iter).first;
418           if (str == particleName)                418           if (str == particleName)
419       {                                           419       {
420         sigmaPointer = (*iter).second;            420         sigmaPointer = (*iter).second; 
421       }                                           421       }
422   }                                               422   }   
423       G4int i;                                    423       G4int i;
424       for (i=0; i<sizeNNstar; i++)                424       for (i=0; i<sizeNNstar; i++)
425   {                                               425   {
426     G4double value = *(sigmaPointer + i) * mil    426     G4double value = *(sigmaPointer + i) * millibarn;
427     G4double energy = energyTable[i] * GeV;       427     G4double energy = energyTable[i] * GeV;
428     sigmaVector->PutValue(i,energy,value);        428     sigmaVector->PutValue(i,energy,value);
429   }                                               429   }           
430       return sigmaVector;                         430       return sigmaVector;
431     }                                             431     }
432   else                                            432   else
433     // No cross section table for the requeste    433     // No cross section table for the requested particle is available in the Map
434     return 0;                                     434     return 0;
435 }                                                 435 }
436                                                   436 
437                                                   437 
438                                                   438 
439                                                   439