//
// ********************************************************************
// * License and Disclaimer                                           *
// *                                                                  *
// * The  Geant4 software  is  copyright of the Copyright Holders  of *
// * the Geant4 Collaboration.  It is provided  under  the terms  and *
// * conditions of the Geant4 Software License,  included in the file *
// * LICENSE and available at  http://cern.ch/geant4/license .  These *
// * include a list of copyright holders.                             *
// *                                                                  *
// * Neither the authors of this software system, nor their employing *
// * institutes,nor the agencies providing financial support for this *
// * work  make  any representation or  warranty, express or implied, *
// * regarding  this  software system or assume any liability for its *
// * use.  Please see the license in the file  LICENSE  and URL above *
// * for the full disclaimer and the limitation of liability.         *
// *                                                                  *
// * This  code  implementation is the result of  the  scientific and *
// * technical work of the GEANT4 collaboration.                      *
// * By using,  copying,  modifying or  distributing the software (or *
// * any work based  on the software)  you  agree  to acknowledge its *
// * use  in  resulting  scientific  publications,  and indicate your *
// * acceptance of all terms of the Geant4 Software license.          *
// ********************************************************************
//
//

#include "globals.hh"
#include "G4ios.hh"
#include "G4XNNstarTable.hh"
#include "G4PhysicsFreeVector.hh"


const G4int G4XNNstarTable::sizeNNstar = 121;


// Energies (GeV) corresponding to the cross section table
// Units are assigned while filling the PhysicsVector

const G4double G4XNNstarTable::energyTable[121] =
{
  0.0, 
  2.014,  2.014,  2.016,  2.018,  2.022,  2.026,  2.031,  2.037,  2.044,  2.052,   
  2.061,  2.071,  2.082,  2.094,  2.107,  2.121,  2.135,  2.151,  2.168,  2.185,   
  2.204,  2.223,  2.244,  2.265,  2.287,  2.311,  2.335,  2.360,  2.386,  2.413,   
  2.441,  2.470,  2.500,  2.531,  2.562,  2.595,  2.629,  2.664,  2.699,  2.736,
  2.773,  2.812,  2.851,  2.891,  2.933,  2.975,  3.018,  3.062,  3.107,  3.153,   
  3.200,  3.248,  3.297,  3.347,  3.397,  3.449,  3.502,  3.555,  3.610,  3.666,   
  3.722,  3.779,  3.838,  3.897,  3.957,  4.018,  4.081,  4.144,  4.208,  4.273,   
  4.339,  4.406,  4.473,  4.542,  4.612,  4.683,  4.754,  4.827,  4.900,  4.975,
  5.000,  6.134,  7.269,  8.403,  9.538, 10.672, 11.807, 12.941, 14.076, 15.210,  
 16.345, 17.479, 18.613, 19.748, 20.882, 22.017, 23.151, 24.286, 25.420, 26.555, 
 27.689, 28.824, 29.958, 31.092, 32.227, 33.361, 34.496, 35.630, 36.765, 37.899,  
 39.034, 40.168, 41.303, 42.437, 43.571, 44.706, 45.840, 46.975, 48.109, 49.244
};

// Cross-sections in mb, from S.A. Bass et al., Prog.Part.Nucl.Phys.41:225-370,1998 
// Units are assigned while filling the PhysicsVector

const G4double G4XNNstarTable::sigmaNN1440[121] = 
{
  0.0, 
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,  
  0.000, 0.000, 0.000, 0.001, 0.001, 0.002, 0.003, 0.004, 0.006, 0.009,  
  0.014, 0.020, 0.031, 0.047, 0.072, 0.113, 0.178, 0.266, 0.370, 0.476,  
  0.575, 0.665, 0.744, 0.814, 0.874, 0.926, 0.971, 1.009, 1.040, 1.066,
  1.087, 1.103, 1.115, 1.124, 1.129, 1.132, 1.131, 1.129, 1.124, 1.117, 
  1.109, 1.099, 1.088, 1.075, 1.062, 1.048, 1.033, 1.017, 1.001, 0.984,  
  0.967, 0.950, 0.932, 0.914, 0.896, 0.879, 0.861, 0.843, 0.825, 0.807,  
  0.790, 0.773, 0.755, 0.738, 0.722, 0.705, 0.689, 0.673, 0.657, 0.641,
  0.636, 0.453, 0.336, 0.258, 0.204, 0.166, 0.137, 0.115, 0.098, 0.084,  
  0.073, 0.064, 0.057, 0.051, 0.046, 0.041, 0.037, 0.034, 0.031, 0.029,  
  0.026, 0.024, 0.022, 0.021, 0.019, 0.018, 0.017, 0.016, 0.015, 0.014, 
  0.013, 0.013, 0.012, 0.011, 0.011, 0.010, 0.010, 0.009, 0.009, 0.008
};

const G4double G4XNNstarTable::sigmaNN1520[121] = 
{
  0.0,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.001, 0.001, 
  0.001, 0.002, 0.003, 0.005, 0.008, 0.014, 0.024, 0.043, 0.079, 0.150, 
  0.280, 0.465, 0.662, 0.841, 0.995, 1.125, 1.235, 1.327, 1.410, 1.474,
  1.525, 1.567, 1.600, 1.625, 1.644, 1.651, 1.657, 1.659, 1.657, 1.651, 
  1.642, 1.629, 1.614, 1.597, 1.578, 1.558, 1.536, 1.512, 1.488, 1.463, 
  1.437, 1.411, 1.384, 1.357, 1.329, 1.302, 1.274, 1.247, 1.219, 1.192, 
  1.165, 1.138, 1.112, 1.086, 1.060, 1.035, 1.010, 0.985, 0.961, 0.938,
  0.930, 0.652, 0.479, 0.365, 0.287, 0.232, 0.191, 0.160, 0.135, 0.116, 
  0.101, 0.089, 0.078, 0.070, 0.062, 0.056, 0.051, 0.046, 0.042, 0.039, 
  0.036, 0.033, 0.031, 0.028, 0.026, 0.025, 0.023, 0.022, 0.020, 0.019, 
  0.018, 0.017, 0.016, 0.015, 0.015, 0.014, 0.013, 0.013, 0.012, 0.011
};

const G4double G4XNNstarTable::sigmaNN1535[121] = 
{
  0.0,
  0.000, 0.000, 0.001, 0.001, 0.001, 0.002, 0.002, 0.003,
  0.004, 0.005, 0.006, 0.008, 0.010, 0.012, 0.015, 0.019,
  0.024, 0.031, 0.039, 0.052, 0.069, 0.097, 0.145, 0.216,
  0.298, 0.378, 0.451, 0.513, 0.566, 0.610, 0.646, 0.675,
  0.699, 0.718, 0.732, 0.742, 0.749, 0.753, 0.754, 0.753,
  0.751, 0.746, 0.740, 0.733, 0.724, 0.715, 0.705, 0.694,
  0.683, 0.671, 0.659, 0.647, 0.634, 0.621, 0.608, 0.595,
  0.582, 0.569, 0.556, 0.543, 0.531, 0.518, 0.506, 0.493,
  0.481, 0.469, 0.458, 0.446, 0.435, 0.424, 0.413, 0.402,
  0.399, 0.276, 0.201, 0.153, 0.120, 0.096, 0.079, 0.066,
  0.056, 0.048, 0.042, 0.037, 0.032, 0.029, 0.026, 0.023,
  0.021, 0.019, 0.017, 0.016, 0.015, 0.014, 0.013, 0.012,
  0.011, 0.010, 0.009, 0.009, 0.008, 0.008, 0.007, 0.007,
  0.007, 0.006, 0.006, 0.006, 0.005, 0.005, 0.005, 0.005
};

const G4double G4XNNstarTable::sigmaNN1650[121] = 
{
  0.0,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.001, 0.001, 0.001, 0.001,
  0.002, 0.002, 0.002, 0.003, 0.004, 0.005, 0.006, 0.007,
  0.008, 0.010, 0.012, 0.015, 0.019, 0.023, 0.029, 0.038,
  0.051, 0.071, 0.104, 0.150, 0.201, 0.249, 0.290, 0.326,
  0.354, 0.378, 0.397, 0.412, 0.424, 0.434, 0.440, 0.445,
  0.448, 0.449, 0.449, 0.448, 0.445, 0.442, 0.438, 0.433,
  0.428, 0.422, 0.416, 0.409, 0.403, 0.395, 0.388, 0.381,
  0.373, 0.366, 0.358, 0.350, 0.343, 0.335, 0.327, 0.320,
  0.312, 0.305, 0.298, 0.291, 0.284, 0.277, 0.270, 0.263,
  0.261, 0.182, 0.133, 0.102, 0.080, 0.064, 0.053, 0.044,
  0.037, 0.032, 0.028, 0.024, 0.022, 0.019, 0.017, 0.015,
  0.014, 0.013, 0.012, 0.011, 0.010, 0.009, 0.008, 0.008,
  0.007, 0.007, 0.006, 0.006, 0.006, 0.005, 0.005, 0.005,
  0.004, 0.004, 0.004, 0.004, 0.004, 0.003, 0.003, 0.003
};

const G4double G4XNNstarTable::sigmaNN1675[121] = 
{
  0.0,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.001,
  0.001, 0.001, 0.002, 0.002, 0.003, 0.005, 0.009, 0.015,
  0.026, 0.048, 0.095, 0.189, 0.324, 0.463, 0.586, 0.691,
  0.780, 0.855, 0.919, 0.972, 1.016, 1.052, 1.081, 1.105,
  1.123, 1.136, 1.145, 1.151, 1.153, 1.153, 1.149, 1.144,
  1.136, 1.127, 1.116, 1.104, 1.090, 1.076, 1.061, 1.045,
  1.028, 1.011, 0.993, 0.975, 0.957, 0.939, 0.921, 0.902,
  0.884, 0.865, 0.847, 0.828, 0.810, 0.792, 0.775, 0.757,
  0.751, 0.538, 0.399, 0.307, 0.242, 0.196, 0.162, 0.136,
  0.115, 0.099, 0.086, 0.076, 0.067, 0.060, 0.053, 0.048,
  0.044, 0.040, 0.036, 0.033, 0.031, 0.028, 0.026, 0.024,
  0.023, 0.021, 0.020, 0.019, 0.018, 0.016, 0.016, 0.015,
  0.014, 0.013, 0.013, 0.012, 0.011, 0.011, 0.010, 0.010 
};

const G4double G4XNNstarTable::sigmaNN1680[121] = 
{
  0.0,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.001, 0.001, 0.001, 0.002, 0.003, 0.006, 0.010,
  0.018, 0.035, 0.073, 0.156, 0.294, 0.446, 0.580, 0.693,
  0.788, 0.867, 0.933, 0.988, 1.033, 1.070, 1.100, 1.124,
  1.142, 1.155, 1.163, 1.168, 1.170, 1.168, 1.164, 1.158,
  1.149, 1.139, 1.127, 1.114, 1.100, 1.085, 1.068, 1.052,
  1.034, 1.016, 0.998, 0.979, 0.960, 0.941, 0.922, 0.903,
  0.884, 0.865, 0.846, 0.827, 0.809, 0.790, 0.772, 0.754,
  0.748, 0.533, 0.394, 0.301, 0.238, 0.192, 0.158, 0.133,
  0.113, 0.097, 0.084, 0.074, 0.065, 0.058, 0.052, 0.047,
  0.042, 0.039, 0.035, 0.032, 0.030, 0.028, 0.026, 0.024,
  0.022, 0.021, 0.019, 0.018, 0.017, 0.016, 0.015, 0.014,
  0.013, 0.013, 0.012, 0.012, 0.011, 0.010, 0.010, 0.010
};

const G4double G4XNNstarTable::sigmaNN1700[121] = 
{
  0.0,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.001, 0.001, 0.001, 0.002, 0.003, 0.005, 0.008,
  0.013, 0.022, 0.038, 0.070, 0.140, 0.245, 0.347, 0.431,
  0.500, 0.556, 0.601, 0.637, 0.666, 0.689, 0.708, 0.721,
  0.731, 0.738, 0.742, 0.743, 0.743, 0.740, 0.736, 0.730,
  0.723, 0.716, 0.707, 0.697, 0.687, 0.676, 0.665, 0.654,
  0.642, 0.630, 0.617, 0.605, 0.593, 0.580, 0.568, 0.555,
  0.543, 0.531, 0.519, 0.507, 0.495, 0.483, 0.471, 0.460,
  0.456, 0.322, 0.236, 0.180, 0.142, 0.114, 0.094, 0.079,
  0.067, 0.057, 0.050, 0.044, 0.039, 0.034, 0.031, 0.028,
  0.025, 0.023, 0.021, 0.019, 0.018, 0.016, 0.015, 0.014,
  0.013, 0.012, 0.011, 0.011, 0.010, 0.009, 0.009, 0.008,
  0.008, 0.008, 0.007, 0.007, 0.007, 0.006, 0.006, 0.006
};

const G4double G4XNNstarTable::sigmaNN1710[121] = 
{
  0.0,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.001, 0.001, 0.001, 0.002, 0.003,
  0.005, 0.008, 0.013, 0.025, 0.052, 0.096, 0.144, 0.185,
  0.219, 0.246, 0.269, 0.288, 0.304, 0.316, 0.326, 0.334,
  0.340, 0.344, 0.347, 0.349, 0.349, 0.349, 0.348, 0.346,
  0.343, 0.340, 0.336, 0.332, 0.328, 0.323, 0.319, 0.313,
  0.308, 0.303, 0.297, 0.292, 0.286, 0.280, 0.275, 0.269,
  0.263, 0.257, 0.252, 0.246, 0.241, 0.235, 0.230, 0.224,
  0.223, 0.158, 0.117, 0.090, 0.071, 0.057, 0.047, 0.040,
  0.034, 0.029, 0.025, 0.022, 0.019, 0.017, 0.016, 0.014,
  0.013, 0.012, 0.011, 0.010, 0.009, 0.008, 0.008, 0.007,
  0.007, 0.006, 0.006, 0.005, 0.005, 0.005, 0.005, 0.004,
  0.004, 0.004, 0.004, 0.003, 0.003, 0.003, 0.003, 0.003
};

const G4double G4XNNstarTable::sigmaNN1720[121] = 
{
  0.0,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.001, 0.001, 0.002, 0.002, 0.003, 0.005, 0.007, 0.010,
  0.015, 0.023, 0.036, 0.061, 0.106, 0.174, 0.250, 0.321,
  0.382, 0.434, 0.477, 0.513, 0.543, 0.568, 0.589, 0.605,
  0.617, 0.627, 0.634, 0.639, 0.641, 0.642, 0.641, 0.639,
  0.635, 0.630, 0.625, 0.618, 0.611, 0.603, 0.595, 0.586,
  0.576, 0.567, 0.557, 0.547, 0.537, 0.527, 0.516, 0.506,
  0.496, 0.485, 0.475, 0.465, 0.455, 0.444, 0.435, 0.425,
  0.421, 0.302, 0.224, 0.172, 0.136, 0.117, 0.091, 0.076,
  0.065, 0.056, 0.049, 0.043, 0.038, 0.034, 0.030, 0.027,
  0.025, 0.022, 0.020, 0.019, 0.017, 0.016, 0.015, 0.014,
  0.013, 0.012, 0.011, 0.010, 0.010, 0.009, 0.009, 0.008,
  0.008, 0.007, 0.007, 0.007, 0.006, 0.006, 0.006, 0.006 
};

const G4double G4XNNstarTable::sigmaNN1900[121] = 
{
  0.0,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.001,
  0.001, 0.001, 0.001, 0.002, 0.002, 0.003, 0.003, 0.005,
  0.006, 0.008, 0.010, 0.014, 0.019, 0.026, 0.037, 0.054,
  0.074, 0.094, 0.114, 0.131, 0.147, 0.161, 0.173, 0.184,
  0.193, 0.201, 0.208, 0.213, 0.218, 0.221, 0.224, 0.226,
  0.228, 0.229, 0.229, 0.229, 0.229, 0.228, 0.227, 0.225,
  0.223, 0.221, 0.219, 0.217, 0.214, 0.212, 0.209, 0.206,
  0.203, 0.200, 0.197, 0.194, 0.190, 0.187, 0.184, 0.181,
  0.180, 0.137, 0.106, 0.083, 0.067, 0.056, 0.047, 0.039,
  0.034, 0.029, 0.026, 0.023, 0.020, 0.018, 0.016, 0.015,
  0.013, 0.012, 0.011, 0.010, 0.010, 0.009, 0.008, 0.008,
  0.007, 0.007, 0.006, 0.006, 0.006, 0.005, 0.005, 0.005,
  0.004, 0.004, 0.004, 0.004, 0.004, 0.003, 0.003, 0.003 
};

const G4double G4XNNstarTable::sigmaNN1990[121] = 
{
  0.0,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.001, 0.001, 0.002, 0.003, 0.005, 0.009,
  0.017, 0.030, 0.051, 0.076, 0.102, 0.127, 0.149, 0.168,
  0.185, 0.199, 0.211, 0.221, 0.230, 0.237, 0.243, 0.247,
  0.251, 0.254, 0.256, 0.258, 0.258, 0.259, 0.258, 0.258,
  0.257, 0.255, 0.254, 0.252, 0.249, 0.247, 0.245, 0.242,
  0.239, 0.236, 0.233, 0.230, 0.226, 0.223, 0.220, 0.216,
  0.215, 0.167, 0.131, 0.104, 0.085, 0.070, 0.059, 0.050,
  0.043, 0.038, 0.033, 0.029, 0.026, 0.023, 0.021, 0.019,
  0.017, 0.016, 0.014, 0.013, 0.012, 0.011, 0.010, 0.010,
  0.009, 0.009, 0.008, 0.007, 0.007, 0.007, 0.006, 0.006,
  0.006, 0.005, 0.005, 0.005, 0.005, 0.004, 0.004, 0.004
};

const G4double G4XNNstarTable::sigmaNN2090[121] = 
{
  0.0,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.001, 0.001, 0.001, 0.002, 0.002, 0.003, 0.005,
  0.007, 0.011, 0.016, 0.024, 0.036, 0.053, 0.071, 0.089,
  0.106, 0.120, 0.133, 0.142, 0.151, 0.158, 0.164, 0.169,
  0.172, 0.175, 0.178, 0.179, 0.180, 0.180, 0.180, 0.180,
  0.179, 0.178, 0.176, 0.175, 0.173, 0.171, 0.169, 0.166,
  0.164, 0.162, 0.159, 0.156, 0.154, 0.151, 0.148, 0.145,
  0.144, 0.107, 0.081, 0.063, 0.050, 0.041, 0.034, 0.028,
  0.024, 0.021, 0.018, 0.016, 0.014, 0.013, 0.011, 0.010,
  0.009, 0.008, 0.008, 0.007, 0.006, 0.006, 0.006, 0.005,
  0.005, 0.004, 0.004, 0.004, 0.004, 0.003, 0.003, 0.003,
  0.003, 0.003, 0.003, 0.003, 0.002, 0.002, 0.002, 0.002
};

const G4double G4XNNstarTable::sigmaNN2190[121] = 
{
  0.0,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.001,
  0.001, 0.002, 0.003, 0.005, 0.009, 0.015, 0.024, 0.036,
  0.050, 0.064, 0.078, 0.090, 0.100, 0.110, 0.118, 0.125,
  0.131, 0.136, 0.140, 0.143, 0.146, 0.149, 0.150, 0.152,
  0.153, 0.153, 0.153, 0.153, 0.153, 0.152, 0.152, 0.151,
  0.150, 0.148, 0.147, 0.145, 0.144, 0.142, 0.140, 0.139,
  0.138, 0.110, 0.087, 0.069, 0.056, 0.047, 0.039, 0.033,
  0.029, 0.025, 0.022, 0.019, 0.017, 0.015, 0.014, 0.013,
  0.012, 0.010, 0.009, 0.009, 0.008, 0.007, 0.007, 0.006,
  0.006, 0.006, 0.005, 0.005, 0.005, 0.004, 0.004, 0.004,
  0.004, 0.003, 0.003, 0.003, 0.003, 0.003, 0.003, 0.003
};

const G4double G4XNNstarTable::sigmaNN2220[121] = 
{
  0.0,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.001, 0.001, 0.001, 0.003, 0.005, 0.008, 0.015, 0.026,
  0.041, 0.058, 0.074, 0.089, 0.101, 0.113, 0.122, 0.131,
  0.138, 0.144, 0.149, 0.153, 0.156, 0.159, 0.161, 0.163,
  0.164, 0.165, 0.165, 0.165, 0.165, 0.165, 0.164, 0.163,
  0.162, 0.161, 0.160, 0.158, 0.157, 0.155, 0.153, 0.151,
  0.150, 0.121, 0.096, 0.077, 0.062, 0.052, 0.043, 0.037,
  0.032, 0.028, 0.024, 0.022, 0.019, 0.017, 0.015, 0.014,
  0.013, 0.011, 0.011, 0.010, 0.009, 0.008, 0.008, 0.007,
  0.007, 0.006, 0.006, 0.005, 0.005, 0.005, 0.005, 0.004,
  0.004, 0.004, 0.004, 0.004, 0.003, 0.003, 0.003, 0.003 
};

const G4double G4XNNstarTable::sigmaNN2250[121] = 
{
  0.0,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
  0.000, 0.001, 0.001, 0.002, 0.003, 0.005, 0.009, 0.016,
  0.027, 0.043, 0.060, 0.076, 0.090, 0.103, 0.114, 0.123,
  0.131, 0.138, 0.144, 0.149, 0.153, 0.156, 0.159, 0.161,
  0.163, 0.164, 0.165, 0.165, 0.165, 0.165, 0.165, 0.164,
  0.163, 0.162, 0.161, 0.159, 0.158, 0.156, 0.154, 0.152,
  0.152, 0.122, 0.096, 0.077, 0.062, 0.051, 0.043, 0.037,
  0.031, 0.027, 0.024, 0.022, 0.019, 0.017, 0.015, 0.014,
  0.012, 0.011, 0.010, 0.009, 0.009, 0.008, 0.008, 0.007,
  0.007, 0.006, 0.006, 0.005, 0.005, 0.005, 0.004, 0.004,
  0.004, 0.004, 0.004, 0.003, 0.003, 0.003, 0.003, 0.003
};


G4XNNstarTable::G4XNNstarTable() 
{
  xMap["N(1440)0"] = (G4double*) sigmaNN1440;
  xMap["N(1440)+"] = (G4double*) sigmaNN1440;
  
  xMap["N(1520)0"] = (G4double*) sigmaNN1520;
  xMap["N(1520)+"] = (G4double*) sigmaNN1520;
  
  xMap["N(1535)0"] = (G4double*) sigmaNN1535;
  xMap["N(1535)+"] = (G4double*) sigmaNN1535;
  
  xMap["N(1650)0"] = (G4double*) sigmaNN1650;
  xMap["N(1650)+"] = (G4double*) sigmaNN1650;
  
  xMap["N(1675)0"] = (G4double*) sigmaNN1675;
  xMap["N(1675)+"] = (G4double*) sigmaNN1675;
  
  xMap["N(1680)0"] = (G4double*) sigmaNN1680;
  xMap["N(1680)+"] = (G4double*) sigmaNN1680;
  
  xMap["N(1700)0"] = (G4double*) sigmaNN1700;
  xMap["N(1700)+"] = (G4double*) sigmaNN1700;
  
  xMap["N(1710)0"] = (G4double*) sigmaNN1710;
  xMap["N(1710)+"] = (G4double*) sigmaNN1710;
  
  xMap["N(1720)0"] = (G4double*) sigmaNN1720;
  xMap["N(1720)+"] = (G4double*) sigmaNN1720;
  
  xMap["N(1900)0"] = (G4double*) sigmaNN1900;
  xMap["N(1900)+"] = (G4double*) sigmaNN1900;
    
  xMap["N(1990)0"] = (G4double*) sigmaNN1990;
  xMap["N(1990)+"] = (G4double*) sigmaNN1990;
  
  xMap["N(2090)0"] = (G4double*) sigmaNN2090;
  xMap["N(2090)+"] = (G4double*) sigmaNN2090;
  
  xMap["N(2190)0"] = (G4double*) sigmaNN2190;
  xMap["N(2190)+"] = (G4double*) sigmaNN2190;
  
  xMap["N(2220)0"] = (G4double*) sigmaNN2220;
  xMap["N(2220)+"] = (G4double*) sigmaNN2220;
  
  xMap["N(2250)0"] = (G4double*) sigmaNN2250;
  xMap["N(2250)+"] = (G4double*) sigmaNN2250;  
}


G4XNNstarTable::~G4XNNstarTable()
{ }


const G4PhysicsVector* G4XNNstarTable::CrossSectionTable(const G4String& particleName) const 
{
  // NOTE: the returned pointer is owned by the client

  if (xMap.find(particleName) != xMap.end())
    {
      // Cross section table for the requested particle available in the Map
      G4PhysicsFreeVector* sigmaVector = new G4PhysicsFreeVector(sizeNNstar);
      std::map <G4String, G4double*, std::less<G4String> >::const_iterator iter;
      G4double* sigmaPointer = 0;
      for (iter = xMap.begin(); iter != xMap.end(); ++iter)
	{
	  G4String str = (*iter).first;
          if (str == particleName)
	    {
	      sigmaPointer = (*iter).second; 
	    }
	}   
      G4int i;
      for (i=0; i<sizeNNstar; i++)
	{
	  G4double value = *(sigmaPointer + i) * millibarn;
	  G4double energy = energyTable[i] * GeV;
	  sigmaVector->PutValue(i,energy,value);
	}	          
      return sigmaVector;
    }
  else
    // No cross section table for the requested particle is available in the Map
    return 0;
}