Geant4 9.6.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4IonYangFluctuationModel.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// -------------------------------------------------------------------
28//
29// GEANT4 Class file
30//
31//
32// File name: G4IonYangFluctuationModel
33//
34// Author: V.Ivanchenko ([email protected])
35//
36// Creation date: 18 August 2000
37//
38// Modifications:
39// 18/08/2000 V.Ivanchenko First implementation
40// 04/09/2000 V.Ivanchenko Rename fluctuations
41// 03/10/2000 V.Ivanchenko CodeWizard clean up
42// 10/05/2001 V.Ivanchenko Clean up againist Linux compilation with -Wall
43//
44// -------------------------------------------------------------------
45// Class Description:
46//
47// The aproximation of additional ion energy loss fluctuations
48// Q.Yang et al., NIM B61(1991)149-155.
49//
50// Class Description: End
51//
52// -------------------------------------------------------------------
53//
54//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
55
57
58#include "globals.hh"
60#include "G4SystemOfUnits.hh"
61#include "G4DynamicParticle.hh"
63#include "G4Material.hh"
64
65//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
66
68 : G4VLowEnergyModel(name)
69{;}
70
71//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
72
74{;}
75
76//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
77
79 const G4Material* material)
80{
81 G4double energy = particle->GetKineticEnergy() ;
82 G4double mass = particle->GetMass() ;
83 G4double charge = (particle->GetCharge())/eplus ;
84
85 G4double q = YangFluctuationModel(material,energy,mass,charge) ;
86
87 return q ;
88}
89
90//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
91
93 const G4Material* material,
94 G4double kineticEnergy)
95{
96 G4double mass = aParticle->GetPDGMass() ;
97 G4double charge = (aParticle->GetPDGCharge())/eplus ;
98
99 G4double q = YangFluctuationModel(material,kineticEnergy,mass,charge);
100
101 return q ;
102}
103
104//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
105
107 const G4Material*) const
108{
109 return 1.0*TeV ;
110}
111
112//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
113
115 const G4Material* ) const
116{
117 return 0.0 ;
118}
119
120//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
121
123{
124 return 1.0*TeV ;
125}
126
127//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
128
130 const G4ParticleDefinition* ) const
131{
132 return 0.0 ;
133}
134
135//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
136
138 const G4Material* ) const
139{
140 return true ;
141}
142
143//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
144
146 const G4Material* ) const
147{
148 return true ;
149}
150
151//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
152
153G4double G4IonYangFluctuationModel::YangFluctuationModel(const G4Material* material,
154 G4double kineticEnergy,
155 G4double particleMass,
156 G4double charge) const
157{
158 // The aproximation of energy loss fluctuations
159 // Q.Yang et al., NIM B61(1991)149-155.
160
161 // Reduced energy in MeV/AMU
162 G4double energy = kineticEnergy *amu_c2/(particleMass*MeV) ;
163
164 G4int i = 0 ;
165 G4double factor = 1.0 ;
166
167 // The index of set of parameters i = 0 for protons(hadrons) in gases
168 // 1 for protons(hadrons) in solids
169 // 2 for ions in atomic gases
170 // 3 for ions in molecular gases
171 // 4 for ions in solids
172 static G4double b[5][4] = {
173 {0.1014, 0.3700, 0.9642, 3.987},
174 {0.1955, 0.6941, 2.522, 1.040},
175 {0.05058, 0.08975, 0.1419, 10.80},
176 {0.05009, 0.08660, 0.2751, 3.787},
177 {0.01273, 0.03458, 0.3951, 3.812}
178 } ;
179
180 // protons (hadrons)
181 if(1.5 > charge) {
182 if( kStateGas != material->GetState() ) i = 1 ;
183
184 // ions
185 } else {
186 G4double zeff = (material->GetElectronDensity())/
187 (material->GetTotNbOfAtomsPerVolume()) ;
188 factor = charge * std::pow(charge/zeff, 0.3333) ;
189
190 if( kStateGas == material->GetState() ) {
191 energy /= (charge * std::sqrt(charge)) ;
192
193 if(1 == (material->GetNumberOfElements())) {
194 i = 2 ;
195 } else {
196 i = 3 ;
197 }
198
199 } else {
200 energy /= (charge * std::sqrt(charge*zeff)) ;
201 i = 4 ;
202 }
203 }
204
205 G4double x = b[i][2] * (1.0 - std::exp( - energy * b[i][3] )) ;
206
207 G4double q = factor * x * b[i][0] /
208 ((energy - b[i][1])*(energy - b[i][1]) + x*x) ;
209
210 return q ;
211}
212
213
214
215
216
217
@ kStateGas
Definition: G4Material.hh:114
double G4double
Definition: G4Types.hh:64
int G4int
Definition: G4Types.hh:66
bool G4bool
Definition: G4Types.hh:67
G4double GetMass() const
G4double GetCharge() const
G4double GetKineticEnergy() const
G4double TheValue(const G4DynamicParticle *particle, const G4Material *material)
G4bool IsInCharge(const G4DynamicParticle *particle, const G4Material *material) const
G4IonYangFluctuationModel(const G4String &name)
G4double HighEnergyLimit(const G4ParticleDefinition *aParticle, const G4Material *material) const
G4double LowEnergyLimit(const G4ParticleDefinition *aParticle, const G4Material *material) const
G4double GetTotNbOfAtomsPerVolume() const
Definition: G4Material.hh:208
G4State GetState() const
Definition: G4Material.hh:180
size_t GetNumberOfElements() const
Definition: G4Material.hh:185
G4double GetElectronDensity() const
Definition: G4Material.hh:216
G4double GetPDGCharge() const