Geant4_10
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 "G4ParticleDefinition.hh"
34 #include "G4Proton.hh"
35 #include "G4Neutron.hh"
36 
37 const G4int G4XNNTotalLowE::tableSize = 29;
38 
39 const G4double G4XNNTotalLowE::ss[29]=
40 {
41  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
42 };
43 
44 const G4double G4XNNTotalLowE::ppTot[29] =
45 {
46  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
47 };
48 
49 const G4double G4XNNTotalLowE::npTot[29] =
50 {
51 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
52 };
53 
55 {
56 
58  G4LowEXsection* np = new G4LowEXsection();
59 
60  G4int i;
61  for (i=0; i<tableSize; i++)
62  {
63  std::pair<double,double> it;
64  it.first=ss[i];
65  it.second=ppTot[i]; pp->push_back(it);
66  it.second=npTot[i]; np->push_back(it);
67  }
68  theCrossSections[G4Proton::ProtonDefinition()] = pp;
69  theCrossSections[G4Neutron::NeutronDefinition()] = np;
70 }
71 
72 
74 {
75 
76  delete theCrossSections[G4Proton::ProtonDefinition()];
77  delete theCrossSections[G4Neutron::NeutronDefinition()];
78 }
79 
81 {
82  G4double result = 0.;
83  G4double sqrtS = (trk1.Get4Momentum() + trk2.Get4Momentum()).mag();
84  G4ParticleDefinition * key = FindKeyParticle(trk1,trk2);
85 
86  if (theCrossSections.find(key)!= theCrossSections.end())
87  {
88  LowEMap::const_iterator iter;
89  for (iter = theCrossSections.begin(); iter != theCrossSections.end(); ++iter)
90  {
91  if ((*iter).first == key)
92  {
93  result = (*iter).second->CrossSection(sqrtS);
94  }
95  }
96  }
97  else
98  {
99  throw G4HadronicException(__FILE__, __LINE__, "G4XNNTotalLowE: particle key out of range");
100  }
101 
102  return result;
103 }
104 
106 {
107  G4String name("NNTotalLowE");
108  return name;
109 }
110 
111 
113 {
114  G4bool result = e>0&&e<3*GeV;
115  return result;
116 }
117 
118 
virtual ~G4XNNTotalLowE()
virtual G4bool IsValid(G4double e) const
static G4Proton * ProtonDefinition()
Definition: G4Proton.cc:88
G4double G4NeutronHPJENDLHEData::G4double result
const XML_Char * name
Definition: expat.h:151
int G4int
Definition: G4Types.hh:78
virtual G4String Name() const
bool G4bool
Definition: G4Types.hh:79
G4ParticleDefinition * FindKeyParticle(const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
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