00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046
00047
00048
00049
00050
00051
00052
00053
00054
00055
00056 #include "G4IonYangFluctuationModel.hh"
00057
00058 #include "globals.hh"
00059 #include "G4PhysicalConstants.hh"
00060 #include "G4SystemOfUnits.hh"
00061 #include "G4DynamicParticle.hh"
00062 #include "G4ParticleDefinition.hh"
00063 #include "G4Material.hh"
00064
00065
00066
00067 G4IonYangFluctuationModel::G4IonYangFluctuationModel(const G4String& name)
00068 : G4VLowEnergyModel(name)
00069 {;}
00070
00071
00072
00073 G4IonYangFluctuationModel::~G4IonYangFluctuationModel()
00074 {;}
00075
00076
00077
00078 G4double G4IonYangFluctuationModel::TheValue(const G4DynamicParticle* particle,
00079 const G4Material* material)
00080 {
00081 G4double energy = particle->GetKineticEnergy() ;
00082 G4double mass = particle->GetMass() ;
00083 G4double charge = (particle->GetCharge())/eplus ;
00084
00085 G4double q = YangFluctuationModel(material,energy,mass,charge) ;
00086
00087 return q ;
00088 }
00089
00090
00091
00092 G4double G4IonYangFluctuationModel::TheValue(const G4ParticleDefinition* aParticle,
00093 const G4Material* material,
00094 G4double kineticEnergy)
00095 {
00096 G4double mass = aParticle->GetPDGMass() ;
00097 G4double charge = (aParticle->GetPDGCharge())/eplus ;
00098
00099 G4double q = YangFluctuationModel(material,kineticEnergy,mass,charge);
00100
00101 return q ;
00102 }
00103
00104
00105
00106 G4double G4IonYangFluctuationModel::HighEnergyLimit(const G4ParticleDefinition*,
00107 const G4Material*) const
00108 {
00109 return 1.0*TeV ;
00110 }
00111
00112
00113
00114 G4double G4IonYangFluctuationModel::LowEnergyLimit(const G4ParticleDefinition* ,
00115 const G4Material* ) const
00116 {
00117 return 0.0 ;
00118 }
00119
00120
00121
00122 G4double G4IonYangFluctuationModel::HighEnergyLimit(const G4ParticleDefinition* ) const
00123 {
00124 return 1.0*TeV ;
00125 }
00126
00127
00128
00129 G4double G4IonYangFluctuationModel::LowEnergyLimit(
00130 const G4ParticleDefinition* ) const
00131 {
00132 return 0.0 ;
00133 }
00134
00135
00136
00137 G4bool G4IonYangFluctuationModel::IsInCharge(const G4DynamicParticle*,
00138 const G4Material* ) const
00139 {
00140 return true ;
00141 }
00142
00143
00144
00145 G4bool G4IonYangFluctuationModel::IsInCharge(const G4ParticleDefinition* ,
00146 const G4Material* ) const
00147 {
00148 return true ;
00149 }
00150
00151
00152
00153 G4double G4IonYangFluctuationModel::YangFluctuationModel(const G4Material* material,
00154 G4double kineticEnergy,
00155 G4double particleMass,
00156 G4double charge) const
00157 {
00158
00159
00160
00161
00162 G4double energy = kineticEnergy *amu_c2/(particleMass*MeV) ;
00163
00164 G4int i = 0 ;
00165 G4double factor = 1.0 ;
00166
00167
00168
00169
00170
00171
00172 static G4double b[5][4] = {
00173 {0.1014, 0.3700, 0.9642, 3.987},
00174 {0.1955, 0.6941, 2.522, 1.040},
00175 {0.05058, 0.08975, 0.1419, 10.80},
00176 {0.05009, 0.08660, 0.2751, 3.787},
00177 {0.01273, 0.03458, 0.3951, 3.812}
00178 } ;
00179
00180
00181 if(1.5 > charge) {
00182 if( kStateGas != material->GetState() ) i = 1 ;
00183
00184
00185 } else {
00186 G4double zeff = (material->GetElectronDensity())/
00187 (material->GetTotNbOfAtomsPerVolume()) ;
00188 factor = charge * std::pow(charge/zeff, 0.3333) ;
00189
00190 if( kStateGas == material->GetState() ) {
00191 energy /= (charge * std::sqrt(charge)) ;
00192
00193 if(1 == (material->GetNumberOfElements())) {
00194 i = 2 ;
00195 } else {
00196 i = 3 ;
00197 }
00198
00199 } else {
00200 energy /= (charge * std::sqrt(charge*zeff)) ;
00201 i = 4 ;
00202 }
00203 }
00204
00205 G4double x = b[i][2] * (1.0 - std::exp( - energy * b[i][3] )) ;
00206
00207 G4double q = factor * x * b[i][0] /
00208 ((energy - b[i][1])*(energy - b[i][1]) + x*x) ;
00209
00210 return q ;
00211 }
00212
00213
00214
00215
00216
00217