Geant4  10.03.p02
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4XNNTotalLowE.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 #include "globals.hh"
28 #include "G4ios.hh"
29 #include "G4SystemOfUnits.hh"
30 #include "G4XNNTotalLowE.hh"
31 #include "G4KineticTrack.hh"
32 #include "G4ParticleDefinition.hh"
33 #include "G4Proton.hh"
34 #include "G4Neutron.hh"
35 
36 const G4int G4XNNTotalLowE::tableSize = 29;
37 
38 const G4double G4XNNTotalLowE::ss[29]=
39 {
40  1877.05, 1879.58, 1882.69, 1884.8, 1886.95, 1891.89, 1895.27, 1899.9, 1904.67, 1913.36, 1921.18, 1933.56, 1949.82, 1990.11, 2025.14, 2059.03, 2101.89, 2166.47, 2201.01, 2236.36, 2289.27, 2377.17, 2426.86, 2500.18, 2602.91, 2733.62, 2925.49, 3002.71
41 };
42 
43 const G4double G4XNNTotalLowE::ppTot[29] =
44 {
45  2000, 600, 250, 180, 138, 92, 75, 57, 46, 35, 29.5, 25.5, 25, 24, 23.7, 25, 29, 39, 44, 47, 48, 48, 47.5, 47, 45.6, 45, 43.3, 42.9, 42.9
46 };
47 
48 const G4double G4XNNTotalLowE::npTot[29] =
49 {
50 4250, 1380, 770, 585, 465, 300, 232, 175, 140, 100, 80, 63, 50, 40, 35, 34, 34, 36.5, 37., 37.7, 38, 39, 39.8, 40.5, 40.7, 41, 41.2, 41.5, 41.5
51 };
52 
54 {
55 
57  G4LowEXsection* np = new G4LowEXsection();
58 
59  G4int i;
60  for (i=0; i<tableSize; i++)
61  {
62  std::pair<double,double> it;
63  it.first=ss[i];
64  it.second=ppTot[i]; pp->push_back(it);
65  it.second=npTot[i]; np->push_back(it);
66  }
67  theCrossSections[G4Proton::ProtonDefinition()] = pp;
68  theCrossSections[G4Neutron::NeutronDefinition()] = np;
69 }
70 
71 
73 {
74 
75  delete theCrossSections[G4Proton::ProtonDefinition()];
76  delete theCrossSections[G4Neutron::NeutronDefinition()];
77 }
78 
80 {
81  G4double result = 0.;
82  G4double sqrtS = (trk1.Get4Momentum() + trk2.Get4Momentum()).mag();
83  const G4ParticleDefinition * key = FindKeyParticle(trk1,trk2);
84 
85  if (theCrossSections.find(key)!= theCrossSections.end())
86  {
87  LowEMap::const_iterator iter;
88  for (iter = theCrossSections.begin(); iter != theCrossSections.end(); ++iter)
89  {
90  if ((*iter).first == key)
91  {
92  result = (*iter).second->CrossSection(sqrtS);
93  }
94  }
95  }
96  else
97  {
98  throw G4HadronicException(__FILE__, __LINE__, "G4XNNTotalLowE: particle key out of range");
99  }
100 
101  return result;
102 }
103 
105 {
106  G4String name("NNTotalLowE");
107  return name;
108 }
109 
110 
112 {
113  G4bool result = e>0&&e<3*GeV;
114  return result;
115 }
116 
117 
G4double G4ParticleHPJENDLHEData::G4double result
const XML_Char * name
Definition: expat.h:151
const G4ParticleDefinition * FindKeyParticle(const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
virtual ~G4XNNTotalLowE()
virtual G4bool IsValid(G4double e) const
static G4Proton * ProtonDefinition()
Definition: G4Proton.cc:88
int G4int
Definition: G4Types.hh:78
virtual G4String Name() const
bool G4bool
Definition: G4Types.hh:79
static constexpr double GeV
Definition: G4SIunits.hh:217
virtual G4double CrossSection(const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
double G4double
Definition: G4Types.hh:76
const G4LorentzVector & Get4Momentum() const
static G4Neutron * NeutronDefinition()
Definition: G4Neutron.cc:99