| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer << 3 // * DISCLAIMER *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th << 5 // * The following disclaimer summarizes all the specific disclaimers *
6 // * the Geant4 Collaboration. It is provided << 6 // * of contributors to this software. The specific disclaimers,which *
7 // * conditions of the Geant4 Software License << 7 // * govern, are listed with their locations in: *
8 // * LICENSE and available at http://cern.ch/ << 8 // * http://cern.ch/geant4/license *
9 // * include a list of copyright holders. <<
10 // * 9 // * *
11 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file << 14 // * use. *
16 // * for the full disclaimer and the limitatio <<
17 // * 15 // * *
18 // * This code implementation is the result << 16 // * This code implementation is the intellectual property of the *
19 // * technical work of the GEANT4 collaboratio << 17 // * GEANT4 collaboration. *
20 // * By using, copying, modifying or distri << 18 // * By copying, distributing or modifying the Program (or any work *
21 // * any work based on the software) you ag << 19 // * based on the Program) you indicate your acceptance of this *
22 // * use in resulting scientific publicati << 20 // * statement, and all its terms. *
23 // * acceptance of all terms of the Geant4 Sof <<
24 // ******************************************* 21 // ********************************************************************
25 // 22 //
26 23
27 #include "globals.hh" 24 #include "globals.hh"
28 #include "G4ios.hh" 25 #include "G4ios.hh"
29 #include "G4SystemOfUnits.hh" <<
30 #include "G4XNNTotalLowE.hh" 26 #include "G4XNNTotalLowE.hh"
31 #include "G4KineticTrack.hh" 27 #include "G4KineticTrack.hh"
32 #include "G4ParticleDefinition.hh" 28 #include "G4ParticleDefinition.hh"
>> 29 #include "G4ParticleDefinition.hh"
33 #include "G4Proton.hh" 30 #include "G4Proton.hh"
34 #include "G4Neutron.hh" 31 #include "G4Neutron.hh"
35 32
36 const G4int G4XNNTotalLowE::tableSize = 29; 33 const G4int G4XNNTotalLowE::tableSize = 29;
37 34
38 const G4double G4XNNTotalLowE::ss[29]= 35 const G4double G4XNNTotalLowE::ss[29]=
39 { 36 {
40 1877.05, 1879.58, 1882.69, 1884.8, 1886.95, 1 37 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 }; 38 };
42 39
43 const G4double G4XNNTotalLowE::ppTot[29] = 40 const G4double G4XNNTotalLowE::ppTot[29] =
44 { 41 {
45 2000, 600, 250, 180, 138, 92, 75, 57, 46, 35 42 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 }; 43 };
47 44
48 const G4double G4XNNTotalLowE::npTot[29] = 45 const G4double G4XNNTotalLowE::npTot[29] =
49 { 46 {
50 4250, 1380, 770, 585, 465, 300, 232, 175, 140, 47 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 }; 48 };
52 49
53 G4XNNTotalLowE::G4XNNTotalLowE() 50 G4XNNTotalLowE::G4XNNTotalLowE()
54 { 51 {
55 52
56 G4LowEXsection* pp = new G4LowEXsection(); 53 G4LowEXsection* pp = new G4LowEXsection();
57 G4LowEXsection* np = new G4LowEXsection(); 54 G4LowEXsection* np = new G4LowEXsection();
58 55
59 G4int i; 56 G4int i;
60 for (i=0; i<tableSize; i++) 57 for (i=0; i<tableSize; i++)
61 { 58 {
62 std::pair<double,double> it; << 59 G4Pair<double,double> it;
63 it.first=ss[i]; 60 it.first=ss[i];
64 it.second=ppTot[i]; pp->push_back(it); 61 it.second=ppTot[i]; pp->push_back(it);
65 it.second=npTot[i]; np->push_back(it); 62 it.second=npTot[i]; np->push_back(it);
66 } 63 }
67 theCrossSections[G4Proton::ProtonDefinition( << 64 theCrossSections[G4Proton::ProtonDefinition()->GetParticleName()] = pp;
68 theCrossSections[G4Neutron::NeutronDefinitio << 65 theCrossSections[G4Neutron::NeutronDefinition()->GetParticleName()] = np;
69 } 66 }
70 67
71 68
72 G4XNNTotalLowE::~G4XNNTotalLowE() 69 G4XNNTotalLowE::~G4XNNTotalLowE()
73 { 70 {
74 << 71 delete theCrossSections[G4Proton::ProtonDefinition()->GetParticleName()];
75 delete theCrossSections[G4Proton::ProtonDefi << 72 delete theCrossSections[G4Neutron::NeutronDefinition()->GetParticleName()];
76 delete theCrossSections[G4Neutron::NeutronDe <<
77 } 73 }
78 74
79 G4double G4XNNTotalLowE::CrossSection(const G4 75 G4double G4XNNTotalLowE::CrossSection(const G4KineticTrack& trk1, const G4KineticTrack& trk2) const
80 { 76 {
81 G4double result = 0.; 77 G4double result = 0.;
82 G4double sqrtS = (trk1.Get4Momentum() + trk2 78 G4double sqrtS = (trk1.Get4Momentum() + trk2.Get4Momentum()).mag();
83 const G4ParticleDefinition * key = FindKeyPa << 79 G4String key = FindKeyParticle(trk1,trk2);
84 80
85 if (theCrossSections.find(key)!= theCrossSec 81 if (theCrossSections.find(key)!= theCrossSections.end())
86 { 82 {
87 LowEMap::const_iterator iter; 83 LowEMap::const_iterator iter;
88 for (iter = theCrossSections.begin(); iter 84 for (iter = theCrossSections.begin(); iter != theCrossSections.end(); ++iter)
89 { 85 {
90 if ((*iter).first == key) 86 if ((*iter).first == key)
91 { 87 {
92 result = (*iter).second->CrossSection(sqrtS) 88 result = (*iter).second->CrossSection(sqrtS);
93 } 89 }
94 } 90 }
95 } 91 }
96 else 92 else
97 { 93 {
98 throw G4HadronicException(__FILE__, __LINE 94 throw G4HadronicException(__FILE__, __LINE__, "G4XNNTotalLowE: particle key out of range");
99 } 95 }
100 96
101 return result; 97 return result;
102 } 98 }
103 99
104 G4String G4XNNTotalLowE::Name() const 100 G4String G4XNNTotalLowE::Name() const
105 { 101 {
106 G4String name("NNTotalLowE"); 102 G4String name("NNTotalLowE");
107 return name; 103 return name;
108 } 104 }
109 105
110 106
111 G4bool G4XNNTotalLowE::IsValid(G4double e) con 107 G4bool G4XNNTotalLowE::IsValid(G4double e) const
112 { 108 {
113 G4bool result = e>0&&e<3*GeV; 109 G4bool result = e>0&&e<3*GeV;
114 return result; 110 return result;
115 } 111 }
116 112
117 113
118 114