Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/hadronic/models/im_r_matrix/src/G4XDeltaNstarTable.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/G4XDeltaNstarTable.cc (Version 11.3.0) and /processes/hadronic/models/im_r_matrix/src/G4XDeltaNstarTable.cc (Version 8.1)


  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 //      GEANT4 Class file                          28 //      GEANT4 Class file
 29 //                                                 29 //
 30 //                                                 30 //
 31 //      File name:     G4XDeltaNstarTable          31 //      File name:     G4XDeltaNstarTable
 32 //                                                 32 //
 33 //      Author:      Maria Grazia Pia (MariaGr     33 //      Author:      Maria Grazia Pia (MariaGrazia.Pia@genova.infn.it)  
 34 //                                                 34 // 
 35 //      Creation date: 4 June 1999                 35 //      Creation date: 4 June 1999
 36 //                                                 36 //
 37 //      Modifications:                             37 //      Modifications: 
 38 //                                                 38 //
 39 // Hadron Kinetic Model                            39 // Hadron Kinetic Model
 40 // p p -> Delta N* cross section tables            40 // p p -> Delta N* cross section tables
 41 //                                                 41 //
 42 // -------------------------------------------     42 // -------------------------------------------------------------------
 43                                                    43 
 44 #include "globals.hh"                              44 #include "globals.hh"
 45 #include "G4ios.hh"                                45 #include "G4ios.hh"
 46 #include "G4SystemOfUnits.hh"                  << 
 47 #include "G4XDeltaNstarTable.hh"                   46 #include "G4XDeltaNstarTable.hh"
 48 #include "G4PhysicsFreeVector.hh"                  47 #include "G4PhysicsFreeVector.hh"
 49                                                    48 
                                                   >>  49 
 50 const G4int G4XDeltaNstarTable::sizeDeltaNstar     50 const G4int G4XDeltaNstarTable::sizeDeltaNstar = 121;
 51                                                    51 
                                                   >>  52 
 52 // Energies (GeV) corresponding to the cross s     53 // Energies (GeV) corresponding to the cross section table
 53 // Units are assigned when filling the Physics <<  54 // Units are assigned while filling the PhysicsVector
 54                                                    55 
 55 const G4double G4XDeltaNstarTable::energyTable     56 const G4double G4XDeltaNstarTable::energyTable[121] =
 56 {                                                  57 {
 57   0.0,                                             58   0.0, 
 58   2.014, 2.014, 2.016, 2.018, 2.022, 2.026, 2.     59   2.014, 2.014, 2.016, 2.018, 2.022, 2.026, 2.031, 2.037, 2.044, 2.052,  
 59   2.061, 2.071, 2.082, 2.094, 2.107, 2.121, 2.     60   2.061, 2.071, 2.082, 2.094, 2.107, 2.121, 2.135, 2.151, 2.168, 2.185,  
 60   2.204, 2.223, 2.244, 2.265, 2.287, 2.311, 2.     61   2.204, 2.223, 2.244, 2.265, 2.287, 2.311, 2.335, 2.360, 2.386, 2.413,  
 61   2.441, 2.470, 2.500, 2.531, 2.562, 2.595, 2.     62   2.441, 2.470, 2.500, 2.531, 2.562, 2.595, 2.629, 2.664, 2.699, 2.736,
 62   2.773, 2.812, 2.851, 2.891, 2.933, 2.975, 3.     63   2.773, 2.812, 2.851, 2.891, 2.933, 2.975, 3.018, 3.062, 3.107, 3.153,  
 63   3.200, 3.248, 3.297, 3.347, 3.397, 3.449, 3.     64   3.200, 3.248, 3.297, 3.347, 3.397, 3.449, 3.502, 3.555, 3.610, 3.666,  
 64   3.722, 3.779, 3.838, 3.897, 3.957, 4.018, 4.     65   3.722, 3.779, 3.838, 3.897, 3.957, 4.018, 4.081, 4.144, 4.208, 4.273,  
 65   4.339, 4.406, 4.473, 4.542, 4.612, 4.683, 4.     66   4.339, 4.406, 4.473, 4.542, 4.612, 4.683, 4.754, 4.827, 4.900, 4.975,
 66   5.000, 6.134, 7.269, 8.403, 9.538, 10.672, 1     67   5.000, 6.134, 7.269, 8.403, 9.538, 10.672, 11.807, 12.941, 14.076, 15.210, 
 67  16.345, 17.479, 18.613, 19.748, 20.882, 22.01     68  16.345, 17.479, 18.613, 19.748, 20.882, 22.017, 23.151, 24.286, 25.420, 26.555, 
 68  27.689, 28.824, 29.958, 31.092, 32.227, 33.36     69  27.689, 28.824, 29.958, 31.092, 32.227, 33.361, 34.496, 35.630, 36.765, 37.899, 
 69  39.034, 40.168, 41.303, 42.437, 43.571, 44.70     70  39.034, 40.168, 41.303, 42.437, 43.571, 44.706, 45.840, 46.975, 48.109, 49.244
 70 };                                                 71 };
 71                                                    72 
 72 // Cross-sections in mb, from S.A. Bass et al.     73 // Cross-sections in mb, from S.A. Bass et al., Prog.Part.Nucl.Phys.41:225-370,1998 
 73 // Units are assigned when filling the Physics <<  74 // Units are assigned while filling the PhysicsVector
 74                                                    75 
 75 const G4double G4XDeltaNstarTable::sigmaDN1440     76 const G4double G4XDeltaNstarTable::sigmaDN1440[121] = 
 76 {                                                  77 {
 77   0.0,                                             78   0.0, 
 78   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,
 79   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,
 80   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.     81   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
 81   0.000, 0.000, 0.000, 0.001, 0.003, 0.006, 0.     82   0.000, 0.000, 0.000, 0.001, 0.003, 0.006, 0.011, 0.021,
 82   0.041, 0.000, 0.000, 0.002, 0.011, 0.047, 0.     83   0.041, 0.000, 0.000, 0.002, 0.011, 0.047, 0.131, 0.257,
 83   0.408, 0.568, 0.729, 0.886,  1.036,  1.178,      84   0.408, 0.568, 0.729, 0.886,  1.036,  1.178,  1.309,  1.431,
 84   1.542,  1.644,  1.736,  1.819,  1.894,  1.96     85   1.542,  1.644,  1.736,  1.819,  1.894,  1.960,  2.017,  2.068,
 85   2.111,  2.148,  2.178,  2.203,  2.222,  2.23     86   2.111,  2.148,  2.178,  2.203,  2.222,  2.237,  2.247,  2.253,
 86   2.254,  2.252,  2.247,  2.239,  2.228,  2.21     87   2.254,  2.252,  2.247,  2.239,  2.228,  2.215,  2.199,  2.181,
 87   2.161,  2.140,  2.117,  2.093,  2.068,  2.04     88   2.161,  2.140,  2.117,  2.093,  2.068,  2.042,  2.014,  1.986,
 88   1.977,  1.558,  1.224, 0.976, 0.797, 0.654,      89   1.977,  1.558,  1.224, 0.976, 0.797, 0.654, 0.548, 0.465,
 89   0.399, 0.346, 0.303, 0.268, 0.238, 0.213,  0     90   0.399, 0.346, 0.303, 0.268, 0.238, 0.213,  0.192,  0.173,
 90   0.158,  0.144,  0.132,  0.121,  0.112,  0.10     91   0.158,  0.144,  0.132,  0.121,  0.112,  0.103,  0.096,  0.089,
 91   0.083,  0.078,  0.073,  0.069,  0.065,  0.06     92   0.083,  0.078,  0.073,  0.069,  0.065,  0.061,  0.057,  0.054,
 92   0.051,  0.049,  0.046,  0.044,  0.042,  0.04     93   0.051,  0.049,  0.046,  0.044,  0.042,  0.040,  0.038,  0.036 
 93 };                                                 94 };
 94                                                    95 
 95 const G4double G4XDeltaNstarTable::sigmaDN1520     96 const G4double G4XDeltaNstarTable::sigmaDN1520[121] = 
 96 {                                                  97 {
 97   0.0,                                             98   0.0,
 98   0.000,  0.000,  0.000,  0.000,  0.000,  0.00     99   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
 99   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    100   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
100   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    101   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
101   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    102   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.001,  0.001,
102   0.003,  0.006,  0.013,  0.027,  0.000,  0.00    103   0.003,  0.006,  0.013,  0.027,  0.000,  0.001,  0.010,  0.048,
103   0.152,  0.318,  0.513,  0.713,  0.908,  1.09    104   0.152,  0.318,  0.513,  0.713,  0.908,  1.091,  1.262,  1.417,
104   1.559,  1.686,  1.801,  1.903,  1.992,  2.07    105   1.559,  1.686,  1.801,  1.903,  1.992,  2.071,  2.140,  2.198,
105   2.248,  2.289,  2.323,  2.349,  2.369,  2.38    106   2.248,  2.289,  2.323,  2.349,  2.369,  2.383,  2.392,  2.395,
106   2.394,  2.389,  2.380,  2.368,  2.352,  2.33    107   2.394,  2.389,  2.380,  2.368,  2.352,  2.334,  2.313,  2.290,
107   2.266,  2.239,  2.211,  2.182,  2.151,  2.11    108   2.266,  2.239,  2.211,  2.182,  2.151,  2.119,  2.087,  2.054,
108   2.043,  1.571,  1.211,  0.953,  0.769,  0.62    109   2.043,  1.571,  1.211,  0.953,  0.769,  0.627,  0.521,  0.440,
109   0.376,  0.325,  0.284,  0.250,  0.222,  0.19    110   0.376,  0.325,  0.284,  0.250,  0.222,  0.198,  0.178,  0.160,
110   0.146,  0.133,  0.121,  0.112,  0.103,  0.09    111   0.146,  0.133,  0.121,  0.112,  0.103,  0.095,  0.088,  0.082,
111   0.076,  0.071,  0.067,  0.063,  0.059,  0.05    112   0.076,  0.071,  0.067,  0.063,  0.059,  0.056,  0.052,  0.050,
112   0.047,  0.045,  0.042,  0.040,  0.038,  0.03    113   0.047,  0.045,  0.042,  0.040,  0.038,  0.036,  0.035,  0.033 
113 };                                                114 };
114                                                   115 
115 const G4double G4XDeltaNstarTable::sigmaDN1535    116 const G4double G4XDeltaNstarTable::sigmaDN1535[121] = 
116 {                                                 117 {
117   0.0,                                            118   0.0,
118   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    119   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
119   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    120   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
120   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    121   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
121   0.000,  0.000,  0.001,  0.001,  0.002,  0.00    122   0.000,  0.000,  0.001,  0.001,  0.002,  0.004,  0.006,  0.010,
122   0.014,  0.021,  0.030,  0.044,  0.000,  0.00    123   0.014,  0.021,  0.030,  0.044,  0.000,  0.000,  0.002,  0.010,
123   0.039,  0.097,  0.174,  0.257,  0.339,  0.41    124   0.039,  0.097,  0.174,  0.257,  0.339,  0.417,  0.490,  0.556,
124   0.615,  0.667,  0.714,  0.754,  0.790,  0.82    125   0.615,  0.667,  0.714,  0.754,  0.790,  0.820,  0.846,  0.867,
125   0.885,  0.899,  0.910,  0.918,  0.924,  0.92    126   0.885,  0.899,  0.910,  0.918,  0.924,  0.927,  0.928,  0.927,
126   0.924,  0.920,  0.914,  0.907,  0.899,  0.89    127   0.924,  0.920,  0.914,  0.907,  0.899,  0.890,  0.880,  0.870,
127   0.858,  0.847,  0.834,  0.822,  0.809,  0.79    128   0.858,  0.847,  0.834,  0.822,  0.809,  0.795,  0.782,  0.768,
128   0.764,  0.576,  0.439,  0.343,  0.275,  0.22    129   0.764,  0.576,  0.439,  0.343,  0.275,  0.223,  0.185,  0.156,
129   0.133,  0.115,  0.100,  0.088,  0.078,  0.07    130   0.133,  0.115,  0.100,  0.088,  0.078,  0.070,  0.062,  0.056,
130   0.051,  0.047,  0.043,  0.039,  0.036,  0.03    131   0.051,  0.047,  0.043,  0.039,  0.036,  0.033,  0.031,  0.029,
131   0.027,  0.025,  0.023,  0.022,  0.021,  0.01    132   0.027,  0.025,  0.023,  0.022,  0.021,  0.019,  0.018,  0.017,
132   0.016,  0.016,  0.015,  0.014,  0.013,  0.01    133   0.016,  0.016,  0.015,  0.014,  0.013,  0.013,  0.012,  0.012 
133 };                                                134 };
134                                                   135 
135 const G4double G4XDeltaNstarTable::sigmaDN1650    136 const G4double G4XDeltaNstarTable::sigmaDN1650[121] = 
136 {                                                 137 {
137   0.0,                                            138   0.0,
138   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    139   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
139   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    140   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
140   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    141   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
141   0.000,  0.000,  0.000,  0.000,  0.001,  0.00    142   0.000,  0.000,  0.000,  0.000,  0.001,  0.001,  0.002,  0.003,
142   0.004,  0.006,  0.008,  0.011,  0.014,  0.01    143   0.004,  0.006,  0.008,  0.011,  0.014,  0.019,  0.026,  0.000,
143   0.000,  0.001,  0.004,  0.018,  0.047,  0.08    144   0.000,  0.001,  0.004,  0.018,  0.047,  0.083,  0.122,  0.159,
144   0.194,  0.227,  0.256,  0.282,  0.305,  0.32    145   0.194,  0.227,  0.256,  0.282,  0.305,  0.326,  0.343,  0.359,
145   0.372,  0.383,  0.393,  0.400,  0.406,  0.41    146   0.372,  0.383,  0.393,  0.400,  0.406,  0.411,  0.415,  0.417,
146   0.418,  0.418,  0.418,  0.417,  0.415,  0.41    147   0.418,  0.418,  0.418,  0.417,  0.415,  0.412,  0.409,  0.405,
147   0.401,  0.397,  0.392,  0.387,  0.382,  0.37    148   0.401,  0.397,  0.392,  0.387,  0.382,  0.376,  0.371,  0.365,
148   0.363,  0.278,  0.214,  0.168,  0.135,  0.11    149   0.363,  0.278,  0.214,  0.168,  0.135,  0.110,  0.091,  0.077,
149   0.066,  0.057,  0.049,  0.043,  0.038,  0.03    150   0.066,  0.057,  0.049,  0.043,  0.038,  0.034,  0.031,  0.028,
150   0.025,  0.023,  0.021,  0.019,  0.018,  0.01    151   0.025,  0.023,  0.021,  0.019,  0.018,  0.016,  0.015,  0.014,
151   0.013,  0.012,  0.012,  0.011,  0.010,  0.01    152   0.013,  0.012,  0.012,  0.011,  0.010,  0.010,  0.009,  0.009,
152   0.008,  0.008,  0.007,  0.007,  0.007,  0.00    153   0.008,  0.008,  0.007,  0.007,  0.007,  0.006,  0.006,  0.006 
153 };                                                154 };
154                                                   155 
155 const G4double G4XDeltaNstarTable::sigmaDN1675    156 const G4double G4XDeltaNstarTable::sigmaDN1675[121] = 
156 {                                                 157 {
157   0.0,                                            158   0.0,
158   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    159   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
159   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    160   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
160   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    161   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
161   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    162   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
162   0.000,  0.000,  0.001,  0.001,  0.002,  0.00    163   0.000,  0.000,  0.001,  0.001,  0.002,  0.004,  0.007,  0.013,
163   0.000,  0.000,  0.003,  0.021,  0.072,  0.15    164   0.000,  0.000,  0.003,  0.021,  0.072,  0.155,  0.251,  0.350,
164   0.446,  0.538,  0.622,  0.699,  0.771,  0.83    165   0.446,  0.538,  0.622,  0.699,  0.771,  0.835,  0.893,  0.945,
165   0.992,  1.033,  1.068,  1.099,  1.126,  1.14    166   0.992,  1.033,  1.068,  1.099,  1.126,  1.148,  1.167,  1.182,
166   1.194,  1.203,  1.209,  1.212,  1.214,  1.21    167   1.194,  1.203,  1.209,  1.212,  1.214,  1.213,  1.210,  1.205,
167   1.199,  1.191,  1.182,  1.172,  1.161,  1.14    168   1.199,  1.191,  1.182,  1.172,  1.161,  1.148,  1.135,  1.122,
168   1.117,  0.891,  0.700,  0.556,  0.454,  0.37    169   1.117,  0.891,  0.700,  0.556,  0.454,  0.371,  0.310,  0.262,
169   0.225,  0.195,  0.170,  0.150,  0.133,  0.11    170   0.225,  0.195,  0.170,  0.150,  0.133,  0.119,  0.107,  0.097,
170   0.088,  0.080,  0.073,  0.067,  0.062,  0.05    171   0.088,  0.080,  0.073,  0.067,  0.062,  0.057,  0.053,  0.050,
171   0.046,  0.043,  0.040,  0.038,  0.036,  0.03    172   0.046,  0.043,  0.040,  0.038,  0.036,  0.034,  0.032,  0.030,
172   0.028,  0.027,  0.026,  0.024,  0.023,  0.02    173   0.028,  0.027,  0.026,  0.024,  0.023,  0.022,  0.021,  0.020 
173 };                                                174 };
174                                                   175 
175 const G4double G4XDeltaNstarTable::sigmaDN1680    176 const G4double G4XDeltaNstarTable::sigmaDN1680[121] = 
176 {                                                 177 {
177   0.0,                                            178   0.0,
178   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    179   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
179   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    180   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
180   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    181   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
181   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    182   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
182   0.000,  0.000,  0.000,  0.001,  0.001,  0.00    183   0.000,  0.000,  0.000,  0.001,  0.001,  0.002,  0.004,  0.009,
183   0.000,  0.000,  0.003,  0.019,  0.070,  0.15    184   0.000,  0.000,  0.003,  0.019,  0.070,  0.157,  0.258,  0.361,
184   0.461,  0.555,  0.641,  0.720,  0.792,  0.85    185   0.461,  0.555,  0.641,  0.720,  0.792,  0.857,  0.915,  0.966,
185   1.012,  1.051,  1.086,  1.116,  1.141,  1.16    186   1.012,  1.051,  1.086,  1.116,  1.141,  1.162,  1.179,  1.192,
186   1.202,  1.209,  1.214,  1.216,  1.215,  1.21    187   1.202,  1.209,  1.214,  1.216,  1.215,  1.212,  1.208,  1.202,
187   1.194,  1.185,  1.175,  1.163,  1.151,  1.13    188   1.194,  1.185,  1.175,  1.163,  1.151,  1.137,  1.123,  1.108,
188   1.103,  0.870,  0.678,  0.536,  0.437,  0.35    189   1.103,  0.870,  0.678,  0.536,  0.437,  0.355,  0.296,  0.250,
189   0.214,  0.185,  0.162,  0.142,  0.126,  0.11    190   0.214,  0.185,  0.162,  0.142,  0.126,  0.113,  0.101,  0.092,
190   0.083,  0.076,  0.069,  0.064,  0.059,  0.05    191   0.083,  0.076,  0.069,  0.064,  0.059,  0.054,  0.050,  0.047,
191   0.044,  0.041,  0.038,  0.036,  0.034,  0.03    192   0.044,  0.041,  0.038,  0.036,  0.034,  0.032,  0.030,  0.028,
192   0.027,  0.025,  0.024,  0.023,  0.022,  0.02    193   0.027,  0.025,  0.024,  0.023,  0.022,  0.021,  0.020,  0.019
193 };                                                194 };
194                                                   195 
195 const G4double G4XDeltaNstarTable::sigmaDN1700    196 const G4double G4XDeltaNstarTable::sigmaDN1700[121] = 
196 {                                                 197 {
197   0.0,                                            198   0.0,
198   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    199   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
199   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    200   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
200   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    201   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
201   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    202   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
202   0.000,  0.000,  0.000,  0.001,  0.001,  0.00    203   0.000,  0.000,  0.000,  0.001,  0.001,  0.002,  0.004,  0.007,
203   0.013,  0.000,  0.001,  0.006,  0.028,  0.07    204   0.013,  0.000,  0.001,  0.006,  0.028,  0.078,  0.142,  0.209,
204   0.273,  0.332,  0.385,  0.432,  0.475,  0.51    205   0.273,  0.332,  0.385,  0.432,  0.475,  0.512,  0.545,  0.573,
205   0.598,  0.620,  0.638,  0.653,  0.665,  0.67    206   0.598,  0.620,  0.638,  0.653,  0.665,  0.675,  0.683,  0.688,
206   0.692,  0.694,  0.695,  0.694,  0.692,  0.68    207   0.692,  0.694,  0.695,  0.694,  0.692,  0.689,  0.685,  0.680,
207   0.674,  0.668,  0.660,  0.653,  0.645,  0.63    208   0.674,  0.668,  0.660,  0.653,  0.645,  0.636,  0.627,  0.618,
208   0.615,  0.477,  0.368,  0.290,  0.234,  0.19    209   0.615,  0.477,  0.368,  0.290,  0.234,  0.190,  0.158,  0.134,
209   0.114,  0.099,  0.086,  0.076,  0.067,  0.06    210   0.114,  0.099,  0.086,  0.076,  0.067,  0.060,  0.054,  0.049,
210   0.044,  0.040,  0.037,  0.034,  0.031,  0.02    211   0.044,  0.040,  0.037,  0.034,  0.031,  0.029,  0.027,  0.025,
211   0.023,  0.022,  0.020,  0.019,  0.018,  0.01    212   0.023,  0.022,  0.020,  0.019,  0.018,  0.017,  0.016,  0.015,
212   0.014,  0.013,  0.013,  0.012,  0.012,  0.01    213   0.014,  0.013,  0.013,  0.012,  0.012,  0.011,  0.010,  0.010 
213 };                                                214 };
214                                                   215 
215 const G4double G4XDeltaNstarTable::sigmaDN1710    216 const G4double G4XDeltaNstarTable::sigmaDN1710[121] = 
216 {                                                 217 {
217   0.0,                                            218   0.0,
218   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    219   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
219   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    220   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
220   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    221   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
221   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    222   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
222   0.000,  0.000,  0.000,  0.000,  0.001,  0.00    223   0.000,  0.000,  0.000,  0.000,  0.001,  0.001,  0.001,  0.003,
223   0.005,  0.000,  0.000,  0.002,  0.009,  0.02    224   0.005,  0.000,  0.000,  0.002,  0.009,  0.028,  0.056,  0.086,
224   0.116,  0.144,  0.170,  0.194,  0.215,  0.23    225   0.116,  0.144,  0.170,  0.194,  0.215,  0.234,  0.251,  0.266,
225   0.279,  0.290,  0.300,  0.308,  0.315,  0.32    226   0.279,  0.290,  0.300,  0.308,  0.315,  0.321,  0.326,  0.329,
226   0.332,  0.334,  0.335,  0.336,  0.335,  0.33    227   0.332,  0.334,  0.335,  0.336,  0.335,  0.334,  0.333,  0.331,
227   0.329,  0.327,  0.324,  0.320,  0.317,  0.31    228   0.329,  0.327,  0.324,  0.320,  0.317,  0.313,  0.309,  0.305,
228   0.304,  0.239,  0.187,  0.161,  0.121,  0.09    229   0.304,  0.239,  0.187,  0.161,  0.121,  0.098,  0.082,  0.069,
229   0.059,  0.051,  0.045,  0.040,  0.035,  0.03    230   0.059,  0.051,  0.045,  0.040,  0.035,  0.031,  0.028,  0.025,
230   0.023,  0.021,  0.019,  0.018,  0.016,  0.01    231   0.023,  0.021,  0.019,  0.018,  0.016,  0.015,  0.014,  0.013,
231   0.012,  0.011,  0.011,  0.010,  0.009,  0.00    232   0.012,  0.011,  0.011,  0.010,  0.009,  0.009,  0.008,  0.008,
232   0.007,  0.007,  0.007,  0.006,  0.006,  0.00    233   0.007,  0.007,  0.007,  0.006,  0.006,  0.006,  0.006,  0.005 
233 };                                                234 };
234                                                   235 
235 const G4double G4XDeltaNstarTable::sigmaDN1720    236 const G4double G4XDeltaNstarTable::sigmaDN1720[121] = 
236 {                                                 237 {
237   0.0,                                            238   0.0,
238   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    239   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
239   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    240   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
240   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    241   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
241   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    242   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
242   0.000,  0.000,  0.001,  0.001,  0.002,  0.00    243   0.000,  0.000,  0.001,  0.001,  0.002,  0.003,  0.005,  0.008,
243   0.014,  0.000,  0.000,  0.001,  0.009,  0.03    244   0.014,  0.000,  0.000,  0.001,  0.009,  0.033,  0.075,  0.124,
244   0.175,  0.225,  0.271,  0.315,  0.355,  0.39    245   0.175,  0.225,  0.271,  0.315,  0.355,  0.391,  0.423,  0.452,
245   0.478,  0.501,  0.521,  0.538,  0.553,  0.56    246   0.478,  0.501,  0.521,  0.538,  0.553,  0.566,  0.576,  0.585,
246   0.592,  0.597,  0.601,  0.603,  0.604,  0.60    247   0.592,  0.597,  0.601,  0.603,  0.604,  0.604,  0.603,  0.601,
247   0.598,  0.595,  0.590,  0.585,  0.580,  0.57    248   0.598,  0.595,  0.590,  0.585,  0.580,  0.574,  0.568,  0.561,
248   0.559,  0.446,  0.351,  0.279,  0.228,  0.20    249   0.559,  0.446,  0.351,  0.279,  0.228,  0.202,  0.156,  0.132,
249   0.113,  0.098,  0.086,  0.076,  0.067,  0.06    250   0.113,  0.098,  0.086,  0.076,  0.067,  0.060,  0.054,  0.049,
250   0.044,  0.040,  0.037,  0.034,  0.031,  0.02    251   0.044,  0.040,  0.037,  0.034,  0.031,  0.029,  0.027,  0.025,
251   0.023,  0.022,  0.020,  0.019,  0.018,  0.01    252   0.023,  0.022,  0.020,  0.019,  0.018,  0.017,  0.016,  0.015,
252   0.014,  0.014,  0.013,  0.012,  0.012,  0.01    253   0.014,  0.014,  0.013,  0.012,  0.012,  0.011,  0.011,  0.010 
253 };                                                254 };
254                                                   255 
255 const G4double G4XDeltaNstarTable::sigmaDN1900    256 const G4double G4XDeltaNstarTable::sigmaDN1900[121] = 
256 {                                                 257 {
257   0.0,                                            258   0.0,
258   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    259   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
259   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    260   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
260   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    261   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
261   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    262   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
262   0.000,  0.000,  0.000,  0.001,  0.001,  0.00    263   0.000,  0.000,  0.000,  0.001,  0.001,  0.002,  0.002,  0.003,
263   0.005,  0.007,  0.009,  0.014,  0.021,  0.00    264   0.005,  0.007,  0.009,  0.014,  0.021,  0.000,  0.000,  0.001,
264   0.003,  0.009,  0.017,  0.028,  0.038,  0.04    265   0.003,  0.009,  0.017,  0.028,  0.038,  0.049,  0.060,  0.070,
265   0.080,  0.090,  0.099,  0.107,  0.115,  0.12    266   0.080,  0.090,  0.099,  0.107,  0.115,  0.122,  0.128,  0.134,
266   0.140,  0.145,  0.149,  0.153,  0.156,  0.15    267   0.140,  0.145,  0.149,  0.153,  0.156,  0.159,  0.162,  0.164,
267   0.166,  0.168,  0.169,  0.170,  0.170,  0.17    268   0.166,  0.168,  0.169,  0.170,  0.170,  0.171,  0.171,  0.171,
268   0.171,  0.155,  0.131,  0.110,  0.095,  0.07    269   0.171,  0.155,  0.131,  0.110,  0.095,  0.079,  0.067,  0.058,
269   0.051,  0.044,  0.039,  0.035,  0.031,  0.02    270   0.051,  0.044,  0.039,  0.035,  0.031,  0.028,  0.026,  0.023,
270   0.021,  0.019,  0.018,  0.016,  0.016,  0.01    271   0.021,  0.019,  0.018,  0.016,  0.016,  0.015,  0.013,  0.012,
271   0.011,  0.011,  0.010,  0.009,  0.009,  0.00    272   0.011,  0.011,  0.010,  0.009,  0.009,  0.008,  0.008,  0.008,
272   0.007,  0.007,  0.006,  0.006,  0.006,  0.00    273   0.007,  0.007,  0.006,  0.006,  0.006,  0.006,  0.005,  0.005 
273 };                                                274 };
274                                                   275 
275 const G4double G4XDeltaNstarTable::sigmaDN1990    276 const G4double G4XDeltaNstarTable::sigmaDN1990[121] = 
276 {                                                 277 {
277   0.0,                                            278   0.0,
278   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    279   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
279   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    280   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
280   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    281   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
281   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    282   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
282   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    283   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
283   0.000,  0.000,  0.001,  0.001,  0.002,  0.00    284   0.000,  0.000,  0.001,  0.001,  0.002,  0.004,  0.008,  0.016,
284   0.000,  0.000,  0.001,  0.005,  0.013,  0.02    285   0.000,  0.000,  0.001,  0.005,  0.013,  0.024,  0.035,  0.046,
285   0.057,  0.068,  0.078,  0.088,  0.097,  0.10    286   0.057,  0.068,  0.078,  0.088,  0.097,  0.105,  0.113,  0.120,
286   0.126,  0.132,  0.138,  0.143,  0.147,  0.15    287   0.126,  0.132,  0.138,  0.143,  0.147,  0.151,  0.155,  0.158,
287   0.161,  0.164,  0.166,  0.168,  0.169,  0.17    288   0.161,  0.164,  0.166,  0.168,  0.169,  0.171,  0.172,  0.172,
288   0.173,  0.164,  0.143,  0.122,  0.105,  0.08    289   0.173,  0.164,  0.143,  0.122,  0.105,  0.089,  0.076,  0.066,
289   0.058,  0.053,  0.045,  0.040,  0.036,  0.03    290   0.058,  0.053,  0.045,  0.040,  0.036,  0.032,  0.029,  0.027,
290   0.024,  0.022,  0.021,  0.019,  0.018,  0.01    291   0.024,  0.022,  0.021,  0.019,  0.018,  0.016,  0.015,  0.014,
291   0.013,  0.012,  0.012,  0.011,  0.010,  0.01    292   0.013,  0.012,  0.012,  0.011,  0.010,  0.010,  0.009,  0.009,
292   0.008,  0.008,  0.007,  0.007,  0.007,  0.00    293   0.008,  0.008,  0.007,  0.007,  0.007,  0.006,  0.006,  0.006 
293 };                                                294 };
294                                                   295 
295 const G4double G4XDeltaNstarTable::sigmaDN2090    296 const G4double G4XDeltaNstarTable::sigmaDN2090[121] = 
296 {                                                 297 {
297   0.0,                                            298   0.0,
298   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    299   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
299   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    300   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
300   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    301   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
301   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    302   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
302   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    303   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
303   0.000,  0.000,  0.001,  0.001,  0.002,  0.00    304   0.000,  0.000,  0.001,  0.001,  0.002,  0.003,  0.004,  0.006,
304   0.010,  0.000,  0.000,  0.001,  0.005,  0.01    305   0.010,  0.000,  0.000,  0.001,  0.005,  0.013,  0.022,  0.033,
305   0.043,  0.053,  0.063,  0.071,  0.079,  0.08    306   0.043,  0.053,  0.063,  0.071,  0.079,  0.086,  0.093,  0.099,
306   0.104,  0.108,  0.112,  0.116,  0.119,  0.12    307   0.104,  0.108,  0.112,  0.116,  0.119,  0.121,  0.123,  0.125,
307   0.126,  0.127,  0.128,  0.129,  0.129,  0.12    308   0.126,  0.127,  0.128,  0.129,  0.129,  0.129,  0.129,  0.128,
308   0.128,  0.110,  0.090,  0.073,  0.061,  0.04    309   0.128,  0.110,  0.090,  0.073,  0.061,  0.049,  0.041,  0.035,
309   0.030,  0.026,  0.023,  0.020,  0.018,  0.01    310   0.030,  0.026,  0.023,  0.020,  0.018,  0.016,  0.015,  0.013,
310   0.012,  0.011,  0.010,  0.009,  0.008,  0.00    311   0.012,  0.011,  0.010,  0.009,  0.008,  0.008,  0.007,  0.007,
311   0.006,  0.006,  0.006,  0.005,  0.005,  0.00    312   0.006,  0.006,  0.006,  0.005,  0.005,  0.005,  0.004,  0.004,
312   0.004,  0.004,  0.004,  0.003,  0.003,  0.00    313   0.004,  0.004,  0.004,  0.003,  0.003,  0.003,  0.003,  0.003 
313 };                                                314 };
314                                                   315 
315 const G4double G4XDeltaNstarTable::sigmaDN2190    316 const G4double G4XDeltaNstarTable::sigmaDN2190[121] = 
316 {                                                 317 {
317   0.0,                                            318   0.0,
318   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    319   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
319   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    320   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
320   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    321   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
321   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    322   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
322   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    323   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
323   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    324   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.001,  0.001,
324   0.002,  0.003,  0.006,  0.010,  0.000,  0.00    325   0.002,  0.003,  0.006,  0.010,  0.000,  0.000,  0.000,  0.003,
325   0.007,  0.013,  0.019,  0.026,  0.032,  0.03    326   0.007,  0.013,  0.019,  0.026,  0.032,  0.038,  0.044,  0.050,
326   0.055,  0.060,  0.065,  0.069,  0.073,  0.07    327   0.055,  0.060,  0.065,  0.069,  0.073,  0.077,  0.080,  0.083,
327   0.086,  0.088,  0.090,  0.092,  0.094,  0.09    328   0.086,  0.088,  0.090,  0.092,  0.094,  0.095,  0.096,  0.097,
328   0.098,  0.097,  0.086,  0.073,  0.062,  0.05    329   0.098,  0.097,  0.086,  0.073,  0.062,  0.053,  0.046,  0.039,
329   0.034,  0.030,  0.027,  0.024,  0.021,  0.01    330   0.034,  0.030,  0.027,  0.024,  0.021,  0.019,  0.017,  0.016,
330   0.015,  0.013,  0.012,  0.011,  0.010,  0.01    331   0.015,  0.013,  0.012,  0.011,  0.010,  0.010,  0.009,  0.008,
331   0.008,  0.007,  0.007,  0.006,  0.006,  0.00    332   0.008,  0.007,  0.007,  0.006,  0.006,  0.006,  0.005,  0.005,
332   0.005,  0.005,  0.004,  0.004,  0.004,  0.00    333   0.005,  0.005,  0.004,  0.004,  0.004,  0.004,  0.004,  0.003 
333 };                                                334 };
334                                                   335 
335 const G4double G4XDeltaNstarTable::sigmaDN2220    336 const G4double G4XDeltaNstarTable::sigmaDN2220[121] = 
336 {                                                 337 {
337   0.0,                                            338   0.0,
338   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    339   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
339   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    340   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
340   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    341   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
341   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    342   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
342   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    343   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
343   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    344   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
344   0.001,  0.001,  0.003,  0.006,  0.011,  0.00    345   0.001,  0.001,  0.003,  0.006,  0.011,  0.000,  0.000,  0.001,
345   0.005,  0.010,  0.017,  0.024,  0.030,  0.03    346   0.005,  0.010,  0.017,  0.024,  0.030,  0.037,  0.044,  0.050,
346   0.055,  0.061,  0.066,  0.070,  0.074,  0.07    347   0.055,  0.061,  0.066,  0.070,  0.074,  0.078,  0.082,  0.085,
347   0.088,  0.091,  0.093,  0.095,  0.097,  0.09    348   0.088,  0.091,  0.093,  0.095,  0.097,  0.099,  0.100,  0.101,
348   0.102,  0.103,  0.092,  0.079,  0.067,  0.05    349   0.102,  0.103,  0.092,  0.079,  0.067,  0.058,  0.050,  0.043,
349   0.038,  0.033,  0.029,  0.027,  0.023,  0.02    350   0.038,  0.033,  0.029,  0.027,  0.023,  0.021,  0.019,  0.017,
350   0.016,  0.014,  0.013,  0.012,  0.011,  0.01    351   0.016,  0.014,  0.013,  0.012,  0.011,  0.011,  0.010,  0.009,
351   0.009,  0.008,  0.007,  0.007,  0.007,  0.00    352   0.009,  0.008,  0.007,  0.007,  0.007,  0.006,  0.006,  0.006,
352   0.005,  0.005,  0.005,  0.005,  0.004,  0.00    353   0.005,  0.005,  0.005,  0.005,  0.004,  0.004,  0.004,  0.004, 
353 };                                                354 };
354                                                   355 
355 const G4double G4XDeltaNstarTable::sigmaDN2250    356 const G4double G4XDeltaNstarTable::sigmaDN2250[121] = 
356 {                                                 357 {
357   0.0,                                            358   0.0,
358   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    359   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
359   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    360   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
360   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    361   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
361   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    362   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
362   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    363   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
363   0.000,  0.000,  0.000,  0.000,  0.000,  0.00    364   0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,
364   0.001,  0.001,  0.002,  0.003,  0.006,  0.00    365   0.001,  0.001,  0.002,  0.003,  0.006,  0.000,  0.000,  0.000,
365   0.003,  0.007,  0.014,  0.021,  0.028,  0.03    366   0.003,  0.007,  0.014,  0.021,  0.028,  0.035,  0.042,  0.049,
366   0.055,  0.060,  0.066,  0.071,  0.075,  0.07    367   0.055,  0.060,  0.066,  0.071,  0.075,  0.079,  0.083,  0.087,
367   0.090,  0.093,  0.095,  0.098,  0.100,  0.10    368   0.090,  0.093,  0.095,  0.098,  0.100,  0.101,  0.103,  0.104,
368   0.105,  0.105,  0.093,  0.079,  0.067,  0.05    369   0.105,  0.105,  0.093,  0.079,  0.067,  0.057,  0.049,  0.042,
369   0.037,  0.032,  0.029,  0.025,  0.023,  0.02    370   0.037,  0.032,  0.029,  0.025,  0.023,  0.020,  0.018,  0.017,
370   0.015,  0.014,  0.013,  0.012,  0.011,  0.01    371   0.015,  0.014,  0.013,  0.012,  0.011,  0.010,  0.009,  0.009,
371   0.008,  0.008,  0.007,  0.007,  0.006,  0.00    372   0.008,  0.008,  0.007,  0.007,  0.006,  0.006,  0.006,  0.005,
372   0.005,  0.005,  0.005,  0.004,  0.004,  0.00    373   0.005,  0.005,  0.005,  0.004,  0.004,  0.004,  0.004,  0.004 
373 };                                                374 };
374                                                   375 
375                                                   376 
376 G4XDeltaNstarTable::G4XDeltaNstarTable()          377 G4XDeltaNstarTable::G4XDeltaNstarTable() 
377 {                                                 378 {
378   xMap["N(1440)0"] = (G4double*) sigmaDN1440;     379   xMap["N(1440)0"] = (G4double*) sigmaDN1440;
379   xMap["N(1440)+"] = (G4double*) sigmaDN1440;     380   xMap["N(1440)+"] = (G4double*) sigmaDN1440;
380                                                   381   
381   xMap["N(1520)0"] = (G4double*) sigmaDN1520;     382   xMap["N(1520)0"] = (G4double*) sigmaDN1520;
382   xMap["N(1520)+"] = (G4double*) sigmaDN1520;     383   xMap["N(1520)+"] = (G4double*) sigmaDN1520;
383                                                   384   
384   xMap["N(1535)0"] = (G4double*) sigmaDN1535;     385   xMap["N(1535)0"] = (G4double*) sigmaDN1535;
385   xMap["N(1535)+"] = (G4double*) sigmaDN1535;     386   xMap["N(1535)+"] = (G4double*) sigmaDN1535;
386                                                   387   
387   xMap["N(1650)0"] = (G4double*) sigmaDN1650;     388   xMap["N(1650)0"] = (G4double*) sigmaDN1650;
388   xMap["N(1650)+"] = (G4double*) sigmaDN1650;     389   xMap["N(1650)+"] = (G4double*) sigmaDN1650;
389                                                   390   
390   xMap["N(1675)0"] = (G4double*) sigmaDN1675;     391   xMap["N(1675)0"] = (G4double*) sigmaDN1675;
391   xMap["N(1675)+"] = (G4double*) sigmaDN1675;     392   xMap["N(1675)+"] = (G4double*) sigmaDN1675;
392                                                   393   
393   xMap["N(1680)0"] = (G4double*) sigmaDN1680;     394   xMap["N(1680)0"] = (G4double*) sigmaDN1680;
394   xMap["N(1680)+"] = (G4double*) sigmaDN1680;     395   xMap["N(1680)+"] = (G4double*) sigmaDN1680;
395                                                   396   
396   xMap["N(1700)0"] = (G4double*) sigmaDN1700;     397   xMap["N(1700)0"] = (G4double*) sigmaDN1700;
397   xMap["N(1700)+"] = (G4double*) sigmaDN1700;     398   xMap["N(1700)+"] = (G4double*) sigmaDN1700;
398                                                   399   
399   xMap["N(1710)0"] = (G4double*) sigmaDN1710;     400   xMap["N(1710)0"] = (G4double*) sigmaDN1710;
400   xMap["N(1710)+"] = (G4double*) sigmaDN1710;     401   xMap["N(1710)+"] = (G4double*) sigmaDN1710;
401                                                   402   
402   xMap["N(1720)0"] = (G4double*) sigmaDN1720;     403   xMap["N(1720)0"] = (G4double*) sigmaDN1720;
403   xMap["N(1720)+"] = (G4double*) sigmaDN1720;     404   xMap["N(1720)+"] = (G4double*) sigmaDN1720;
404                                                   405   
405   xMap["N(1900)0"] = (G4double*) sigmaDN1900;     406   xMap["N(1900)0"] = (G4double*) sigmaDN1900;
406   xMap["N(1900)+"] = (G4double*) sigmaDN1900;     407   xMap["N(1900)+"] = (G4double*) sigmaDN1900;
407                                                   408     
408   xMap["N(1990)0"] = (G4double*) sigmaDN1990;     409   xMap["N(1990)0"] = (G4double*) sigmaDN1990;
409   xMap["N(1990)+"] = (G4double*) sigmaDN1990;     410   xMap["N(1990)+"] = (G4double*) sigmaDN1990;
410                                                   411   
411   xMap["N(2090)0"] = (G4double*) sigmaDN2090;     412   xMap["N(2090)0"] = (G4double*) sigmaDN2090;
412   xMap["N(2090)+"] = (G4double*) sigmaDN2090;     413   xMap["N(2090)+"] = (G4double*) sigmaDN2090;
413                                                   414   
414   xMap["N(2190)0"] = (G4double*) sigmaDN2190;     415   xMap["N(2190)0"] = (G4double*) sigmaDN2190;
415   xMap["N(2190)+"] = (G4double*) sigmaDN2190;     416   xMap["N(2190)+"] = (G4double*) sigmaDN2190;
416                                                   417   
417   xMap["N(2220)0"] = (G4double*) sigmaDN2220;     418   xMap["N(2220)0"] = (G4double*) sigmaDN2220;
418   xMap["N(2220)+"] = (G4double*) sigmaDN2220;     419   xMap["N(2220)+"] = (G4double*) sigmaDN2220;
419                                                   420   
420   xMap["N(2250)0"] = (G4double*) sigmaDN2250;     421   xMap["N(2250)0"] = (G4double*) sigmaDN2250;
421   xMap["N(2250)+"] = (G4double*) sigmaDN2250;     422   xMap["N(2250)+"] = (G4double*) sigmaDN2250; 
422 }                                                 423 }
423                                                   424 
424                                                   425 
425 G4XDeltaNstarTable::~G4XDeltaNstarTable()         426 G4XDeltaNstarTable::~G4XDeltaNstarTable()
426 { }                                               427 { }
427                                                   428 
428                                                   429 
429 const G4PhysicsVector* G4XDeltaNstarTable::Cro    430 const G4PhysicsVector* G4XDeltaNstarTable::CrossSectionTable(const G4String& particleName) const
430 {                                                 431 {
431    // NOTE: the returned pointer is owned by t    432    // NOTE: the returned pointer is owned by the client
432                                                   433 
433   if (xMap.find(particleName) != xMap.end())      434   if (xMap.find(particleName) != xMap.end())
434     {                                             435     {
435       // Cross section table for the requested    436       // Cross section table for the requested particle available in the Map
436       G4PhysicsFreeVector* sigmaVector = new G    437       G4PhysicsFreeVector* sigmaVector = new G4PhysicsFreeVector(sizeDeltaNstar);
437       std::map <G4String, G4double*, std::less    438       std::map <G4String, G4double*, std::less<G4String> >::const_iterator iter;
438       G4double* sigmaPointer = 0;                 439       G4double* sigmaPointer = 0;
439       for (iter = xMap.begin(); iter != xMap.e    440       for (iter = xMap.begin(); iter != xMap.end(); ++iter)
440   {                                               441   {
441     G4String str = (*iter).first;                 442     G4String str = (*iter).first;
442           if (str == particleName)                443           if (str == particleName)
443       {                                           444       {
444         sigmaPointer = (*iter).second;            445         sigmaPointer = (*iter).second; 
445       }                                           446       }
446   }                                               447   }
447                                                   448 
448       G4int i;                                    449       G4int i;
449       for (i=0; i<sizeDeltaNstar; i++)            450       for (i=0; i<sizeDeltaNstar; i++)
450   {                                               451   {
451     G4double value = *(sigmaPointer + i) * 0.5    452     G4double value = *(sigmaPointer + i) * 0.5* millibarn;
452     G4double energy = energyTable[i] * GeV;       453     G4double energy = energyTable[i] * GeV;
453     sigmaVector->PutValue(i,energy,value);        454     sigmaVector->PutValue(i,energy,value);
454   }                                               455   }     
455       return sigmaVector;                         456       return sigmaVector;
456     }                                             457     }
457   else                                            458   else
458     // No cross section table for the requeste    459     // No cross section table for the requested particle is available in the Map
459     return 0;                                     460     return 0;
460 }                                                 461 }
461                                                   462 
462                                                   463 
463                                                   464 
464                                                   465 
465                                                   466 
466                                                   467 
467                                                   468 
468                                                   469 
469                                                   470 
470                                                   471 
471                                                   472 
472                                                   473 
473                                                   474