| Geant4 Cross Reference |
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 #include "globals.hh" 27 #include "globals.hh"
28 #include "G4ios.hh" 28 #include "G4ios.hh"
29 #include "G4SystemOfUnits.hh" 29 #include "G4SystemOfUnits.hh"
30 #include "G4XNNTotalLowE.hh" 30 #include "G4XNNTotalLowE.hh"
31 #include "G4KineticTrack.hh" 31 #include "G4KineticTrack.hh"
32 #include "G4ParticleDefinition.hh" 32 #include "G4ParticleDefinition.hh"
33 #include "G4Proton.hh" 33 #include "G4Proton.hh"
34 #include "G4Neutron.hh" 34 #include "G4Neutron.hh"
35 35
36 const G4int G4XNNTotalLowE::tableSize = 29; 36 const G4int G4XNNTotalLowE::tableSize = 29;
37 37
38 const G4double G4XNNTotalLowE::ss[29]= 38 const G4double G4XNNTotalLowE::ss[29]=
39 { 39 {
40 1877.05, 1879.58, 1882.69, 1884.8, 1886.95, 1 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 }; 41 };
42 42
43 const G4double G4XNNTotalLowE::ppTot[29] = 43 const G4double G4XNNTotalLowE::ppTot[29] =
44 { 44 {
45 2000, 600, 250, 180, 138, 92, 75, 57, 46, 35 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 }; 46 };
47 47
48 const G4double G4XNNTotalLowE::npTot[29] = 48 const G4double G4XNNTotalLowE::npTot[29] =
49 { 49 {
50 4250, 1380, 770, 585, 465, 300, 232, 175, 140, 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 }; 51 };
52 52
53 G4XNNTotalLowE::G4XNNTotalLowE() 53 G4XNNTotalLowE::G4XNNTotalLowE()
54 { 54 {
55 55
56 G4LowEXsection* pp = new G4LowEXsection(); 56 G4LowEXsection* pp = new G4LowEXsection();
57 G4LowEXsection* np = new G4LowEXsection(); 57 G4LowEXsection* np = new G4LowEXsection();
58 58
59 G4int i; 59 G4int i;
60 for (i=0; i<tableSize; i++) 60 for (i=0; i<tableSize; i++)
61 { 61 {
62 std::pair<double,double> it; 62 std::pair<double,double> it;
63 it.first=ss[i]; 63 it.first=ss[i];
64 it.second=ppTot[i]; pp->push_back(it); 64 it.second=ppTot[i]; pp->push_back(it);
65 it.second=npTot[i]; np->push_back(it); 65 it.second=npTot[i]; np->push_back(it);
66 } 66 }
67 theCrossSections[G4Proton::ProtonDefinition( 67 theCrossSections[G4Proton::ProtonDefinition()] = pp;
68 theCrossSections[G4Neutron::NeutronDefinitio 68 theCrossSections[G4Neutron::NeutronDefinition()] = np;
69 } 69 }
70 70
71 71
72 G4XNNTotalLowE::~G4XNNTotalLowE() 72 G4XNNTotalLowE::~G4XNNTotalLowE()
73 { 73 {
74 74
75 delete theCrossSections[G4Proton::ProtonDefi 75 delete theCrossSections[G4Proton::ProtonDefinition()];
76 delete theCrossSections[G4Neutron::NeutronDe 76 delete theCrossSections[G4Neutron::NeutronDefinition()];
77 } 77 }
78 78
79 G4double G4XNNTotalLowE::CrossSection(const G4 79 G4double G4XNNTotalLowE::CrossSection(const G4KineticTrack& trk1, const G4KineticTrack& trk2) const
80 { 80 {
81 G4double result = 0.; 81 G4double result = 0.;
82 G4double sqrtS = (trk1.Get4Momentum() + trk2 82 G4double sqrtS = (trk1.Get4Momentum() + trk2.Get4Momentum()).mag();
83 const G4ParticleDefinition * key = FindKeyPa 83 const G4ParticleDefinition * key = FindKeyParticle(trk1,trk2);
84 84
85 if (theCrossSections.find(key)!= theCrossSec 85 if (theCrossSections.find(key)!= theCrossSections.end())
86 { 86 {
87 LowEMap::const_iterator iter; 87 LowEMap::const_iterator iter;
88 for (iter = theCrossSections.begin(); iter 88 for (iter = theCrossSections.begin(); iter != theCrossSections.end(); ++iter)
89 { 89 {
90 if ((*iter).first == key) 90 if ((*iter).first == key)
91 { 91 {
92 result = (*iter).second->CrossSection(sqrtS) 92 result = (*iter).second->CrossSection(sqrtS);
93 } 93 }
94 } 94 }
95 } 95 }
96 else 96 else
97 { 97 {
98 throw G4HadronicException(__FILE__, __LINE 98 throw G4HadronicException(__FILE__, __LINE__, "G4XNNTotalLowE: particle key out of range");
99 } 99 }
100 100
101 return result; 101 return result;
102 } 102 }
103 103
104 G4String G4XNNTotalLowE::Name() const 104 G4String G4XNNTotalLowE::Name() const
105 { 105 {
106 G4String name("NNTotalLowE"); 106 G4String name("NNTotalLowE");
107 return name; 107 return name;
108 } 108 }
109 109
110 110
111 G4bool G4XNNTotalLowE::IsValid(G4double e) con 111 G4bool G4XNNTotalLowE::IsValid(G4double e) const
112 { 112 {
113 G4bool result = e>0&&e<3*GeV; 113 G4bool result = e>0&&e<3*GeV;
114 return result; 114 return result;
115 } 115 }
116 116
117 117
118 118