Geant4  10.01
G4XnpElasticLowE.cc
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 //
27 // -------------------------------------------------------------------
28 // GEANT4 Class file
29 //
30 // For information related to this code contact:
31 //
32 // File name: G4XnpElasticLowE
33 //
34 // Author:
35 //
36 // Creation date: 15 April 1999
37 //
38 // Modifications:
39 //
40 // Neutron-Proton elastic cross section
41 // Linear interpolation from data table
42 //
43 // -------------------------------------------------------------------
44 
45 #include "globals.hh"
46 #include "G4ios.hh"
47 #include "G4SystemOfUnits.hh"
48 #include "G4XnpElasticLowE.hh"
49 #include "G4KineticTrack.hh"
50 #include "G4ParticleDefinition.hh"
51 #include "G4PhysicsVector.hh"
52 #include "G4PhysicsLnVector.hh"
53 #include "G4Proton.hh"
54 #include "G4Neutron.hh"
55 
58 
59 // Low energy limit of the cross-section table (in GeV)
60 // Units are assigned while filling the PhysicsVector
61 const G4double G4XnpElasticLowE::_eMinTable = 1.8964808;
63 
64 // Cross-sections in mb
65 // Units are assigned while filling the PhysicsVector
68 {
69  1500.0, //GF was 0.0
70  248.20, 93.38, 55.26, 44.50, 41.33, 38.48, 37.20, 35.98,
71  35.02, 34.47, 32.48, 30.76, 29.46, 28.53, 27.84, 27.20,
72  26.53, 25.95, 25.59, 25.46, 25.00, 24.49, 24.08, 23.86,
73  23.17, 22.70, 21.88, 21.48, 20.22, 19.75, 18.97, 18.39,
74  17.98, 17.63, 17.21, 16.72, 16.68, 16.58, 16.42, 16.22,
75  15.98, 15.71, 15.42, 15.14, 14.87, 14.65, 14.44, 14.26,
76  14.10, 13.95, 13.80, 13.64, 13.47, 13.29, 13.09, 12.89,
77  12.68, 12.47, 12.27, 12.06, 11.84, 11.76, 11.69, 11.60,
78  11.50, 11.41, 11.29, 11.17, 11.06, 10.93, 10.81, 10.68,
79  10.56, 10.44, 10.33, 10.21, 10.12, 10.03, 9.96, 9.89,
80  9.83, 9.80, 9.77, 9.75, 9.74, 9.74, 9.74, 9.76,
81  9.73, 9.70, 9.68, 9.65, 9.63, 9.60, 9.57, 9.55,
82  9.52, 9.49, 9.46, 9.43
83 };
84 
85 
87 {
88  // Cross-sections are available in the range (_eMin,_eMax)
89 
90  _eMin = _eMinTable * GeV;
91  _eMin = std::exp(std::log(_eMinTable)-_eStepLog)*GeV;
92  _eMax = std::exp(std::log(_eMinTable) + _tableSize * _eStepLog) * GeV;
93 
94  // Protections: validity limits must be compatible with available data
95 
96  if (_eMin < _lowLimit)
97  throw G4HadronicException(__FILE__, __LINE__, "G4XnpElasticLowE::G4XnpElasticLowE - Low energy limit not valid");
98 
99  if (_highLimit > _eMax)
100  throw G4HadronicException(__FILE__, __LINE__, "G4XnpElasticLowE::G4XnpElasticLowE - High energy limit not valid");
101 
103  G4int i;
104  for (i=0; i<_tableSize; i++)
105  {
106  G4double value = _sigmaTable[i] * millibarn;
107  _sigma->PutValue(i,value);
108  }
109 }
110 
111 
113 {
114  delete _sigma;
115  _sigma = 0;
116 }
117 
118 
120 {
121  return (this == (G4XnpElasticLowE *) &right);
122 }
123 
124 
126 {
127  return (this != (G4XnpElasticLowE *) &right);
128 }
129 
130 
132 {
133  G4double sigma = 0.;
134  G4double sqrtS = (trk1.Get4Momentum() + trk2.Get4Momentum()).mag();
135  G4bool dummy = false;
136 
139 
140  const G4ParticleDefinition* def1 = trk1.GetDefinition();
141  const G4ParticleDefinition* def2 = trk2.GetDefinition();
142  if ( (def1 == proton && def2 == neutron) ||
143  (def1 == neutron && def2 == proton) )
144  {
145  if (sqrtS >= _eMin && sqrtS <= _eMax)
146  {
147  sigma = _sigma->GetValue(sqrtS,dummy);
148  } else if ( sqrtS < _eMin )
149  {
150  sigma = _sigma->GetValue(_eMin,dummy);
151  }
152  }
153 
154  return sigma;
155 }
156 
158 {
159  // Dump the cross-section table
160  G4cout << Name() << "Cross-section table: " << G4endl;
161  G4bool dummy = false;
162  G4int i;
163 
164  for (i=0; i<_tableSize; i++)
165  {
167  G4double sigma = _sigma->GetValue(e,dummy) / millibarn;
168  G4cout << i << ") e = " << e << " GeV ---- Cross section = " << sigma << " mb " << G4endl;
169  }
170 
172 }
173 
174 
176 {
177  G4String name("npElasticLowE");
178  return name;
179 }
180 
181 
183 {
184  G4bool answer = InLimits(e,_lowLimit,_highLimit);
185 
186  return answer;
187 }
188 
189 
G4double GetValue(G4double theEnergy, G4bool &isOutRange) const
virtual void Print() const
virtual G4bool IsValid(G4double e) const
G4String name
Definition: TRTMaterials.hh:40
static G4Proton * ProtonDefinition()
Definition: G4Proton.cc:88
G4bool operator==(const G4XnpElasticLowE &right) const
G4PhysicsVector * _sigma
static const G4double _sigmaTable[101]
static const G4int _tableSize
static const G4double _eMinTable
G4double GetLowEdgeEnergy(size_t binNumber) const
static const G4double _highLimit
int G4int
Definition: G4Types.hh:78
G4GLOB_DLL std::ostream G4cout
virtual G4double CrossSection(const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
bool G4bool
Definition: G4Types.hh:79
virtual G4String Name() const
void PutValue(size_t index, G4double theValue)
static const double GeV
Definition: G4SIunits.hh:196
virtual ~G4XnpElasticLowE()
static const G4double _eStepLog
G4bool operator!=(const G4XnpElasticLowE &right) const
virtual void Print() const
static const double millibarn
Definition: G4SIunits.hh:96
#define G4endl
Definition: G4ios.hh:61
double G4double
Definition: G4Types.hh:76
const G4LorentzVector & Get4Momentum() const
const G4ParticleDefinition * GetDefinition() const
static G4Neutron * NeutronDefinition()
Definition: G4Neutron.cc:99
G4bool InLimits(G4double e, G4double eLow, G4double eHigh) const
static const G4double _lowLimit