| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 //
27 //
28 //
29 // Geant4 Hadron Charge Exchange Process -- so
30 //
31 // Created 21 April 2006 V.Ivanchenko
32 //
33 // Modified:
34 // 24-Apr-06 V.Ivanchenko add neutron scatteri
35 // 07-Jun-06 V.Ivanchenko fix problem of rotat
36 // 25-Jul-06 V.Ivanchenko add 19 MeV low energ
37 // 23-Jan-07 V.Ivanchenko add cross section in
38 // and do not use CHIPS
39 // 14-Sep-12 M.Kelsey -- Pass subType code to
40 // 06-Aug-15 A.Ribon migrating to G4Pow
41
42 #include "G4ChargeExchangeProcess.hh"
43 #include "globals.hh"
44 #include "G4SystemOfUnits.hh"
45 #include "G4CrossSectionDataStore.hh"
46 #include "G4CrossSectionElastic.hh"
47 #include "G4ComponentGGHadronNucleusXsc.hh"
48 #include "G4Element.hh"
49 #include "G4ElementVector.hh"
50 #include "G4IsotopeVector.hh"
51 #include "G4Neutron.hh"
52 #include "G4Proton.hh"
53 #include "G4PhysicsLinearVector.hh"
54
55 #include "G4Pow.hh"
56
57
58 G4ChargeExchangeProcess::G4ChargeExchangeProce
59 : G4HadronicProcess(procName,fChargeExchange
60 {
61 G4cout << "###=== The class G4ChargeExchange
62 G4cout << "###=== It will be removed at the
63 thEnergy = 20.*MeV;
64 pPDG = 0;
65 verboseLevel= 1;
66 AddDataSet( new G4CrossSectionElastic( new G
67 theProton = G4Proton::Proton();
68 theNeutron = G4Neutron::Neutron();
69 theAProton = G4AntiProton::AntiProton();
70 theANeutron = G4AntiNeutron::AntiNeutron();
71 thePiPlus = G4PionPlus::PionPlus();
72 thePiMinus = G4PionMinus::PionMinus();
73 thePiZero = G4PionZero::PionZero();
74 theKPlus = G4KaonPlus::KaonPlus();
75 theKMinus = G4KaonMinus::KaonMinus();
76 theK0S = G4KaonZeroShort::KaonZeroShort
77 theK0L = G4KaonZeroLong::KaonZeroLong()
78 theL = G4Lambda::Lambda();
79 theAntiL = G4AntiLambda::AntiLambda();
80 theSPlus = G4SigmaPlus::SigmaPlus();
81 theASPlus = G4AntiSigmaPlus::AntiSigmaPlus
82 theSMinus = G4SigmaMinus::SigmaMinus();
83 theASMinus = G4AntiSigmaMinus::AntiSigmaMin
84 theS0 = G4SigmaZero::SigmaZero();
85 theAS0 = G4AntiSigmaZero::AntiSigmaZero
86 theXiMinus = G4XiMinus::XiMinus();
87 theXi0 = G4XiZero::XiZero();
88 theAXiMinus = G4AntiXiMinus::AntiXiMinus();
89 theAXi0 = G4AntiXiZero::AntiXiZero();
90 theOmega = G4OmegaMinus::OmegaMinus();
91 theAOmega = G4AntiOmegaMinus::AntiOmegaMin
92 theD = G4Deuteron::Deuteron();
93 theT = G4Triton::Triton();
94 theA = G4Alpha::Alpha();
95 theHe3 = G4He3::He3();
96 }
97
98 G4ChargeExchangeProcess::~G4ChargeExchangeProc
99 {
100 if (factors) delete factors;
101 }
102
103 void G4ChargeExchangeProcess::
104 BuildPhysicsTable(const G4ParticleDefinition&
105 {
106 if(first) {
107 first = false;
108 theParticle = &aParticleType;
109 pPDG = theParticle->GetPDGEncoding();
110
111 store = G4HadronicProcess::GetCrossSection
112
113 const size_t n = 10;
114 if(theParticle == thePiPlus || theParticle
115 theParticle == theKPlus || theParticle
116 theParticle == theK0S || theParticle
117
118 G4double F[n] = {0.33,0.27,0.29,0.31,0.2
119 factors = new G4PhysicsLinearVector(0.0,
120 for(size_t i=0; i<n; i++) {factors->PutV
121
122 } else {
123
124 G4double F[n] = {0.50,0.45,0.40,0.35,0.3
125 factors = new G4PhysicsLinearVector(0.0,
126 for(size_t i=0; i<n; i++) {factors->PutV
127 }
128
129 if(verboseLevel>1)
130 G4cout << "G4ChargeExchangeProcess for "
131 << theParticle->GetParticleName()
132 << G4endl;
133 }
134 G4HadronicProcess::BuildPhysicsTable(aPartic
135 }
136
137 G4double G4ChargeExchangeProcess::GetElementCr
138 const G4Dyna
139 const G4Element* elm,
140 const G4Material* mat)
141 {
142 // gives the microscopic cross section in GE
143 G4double Z = elm->GetZ();
144 G4int iz = G4int(Z);
145 G4double x = 0.0;
146
147 // The process is effective only above the t
148 if(iz == 1 || dp->GetKineticEnergy() < thEne
149
150 if(verboseLevel>1)
151 G4cout << "G4ChargeExchangeProcess compute
152 << elm->GetName()
153 << G4endl;
154 x = store->GetCrossSection(dp, elm, mat);
155
156 if(verboseLevel>1)
157 G4cout << "G4ChargeExchangeProcess cross(m
158 << " E(MeV)= " << dp->GetKineticEn
159 << " " << theParticle->GetParticleName()
160 << " in Z= " << iz
161 << G4endl;
162 G4bool b;
163 G4double A = elm->GetN();
164 G4double ptot = dp->GetTotalMomentum();
165 x *= factors->GetValue(ptot, b)/G4Pow::GetIn
166 if(theParticle == thePiPlus || theParticle =
167 theParticle == theKPlus || theParticle =
168 { x *= (1.0 - Z/A); }
169
170 else if(theParticle == thePiMinus || thePart
171 theParticle == theKMinus || thePart
172 { x *= Z/A; }
173
174 if(theParticle->GetPDGMass() < GeV) {
175 if(ptot > 2.*GeV) x *= 4.0*GeV*GeV/(ptot*p
176 }
177
178 if(verboseLevel>1)
179 G4cout << "Corrected cross(mb)= " << x/mil
180
181 return x;
182 }
183
184 G4bool G4ChargeExchangeProcess::
185 IsApplicable(const G4ParticleDefinition& aPart
186 {
187 const G4ParticleDefinition* p = &aParticleTy
188 return (p == thePiPlus || p == thePiMinus ||
189 p == theProton || p == theNeutron ||
190 p == theAProton|| p == theANeutron||
191 p == theKPlus || p == theKMinus ||
192 p == theK0S || p == theK0L ||
193 p == theL);
194 }
195
196 void G4ChargeExchangeProcess::
197 DumpPhysicsTable(const G4ParticleDefinition& a
198 {
199 store->DumpPhysicsTable(aParticleType);
200 }
201