source: trunk/source/processes/electromagnetic/lowenergy/src/G4FinalStateChargeDecrease.cc@ 1229

Last change on this file since 1229 was 1228, checked in by garnier, 16 years ago

update geant4.9.3 tag
  • Property svn:executable set to *
File size: 7.7 KB
Line 
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// $Id: G4FinalStateChargeDecrease.cc,v 1.7 2009/06/11 15:47:08 mantero Exp $
27// GEANT4 tag $Name: geant4-09-03 $
28
29#include "G4FinalStateChargeDecrease.hh"
30
31//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
32
33G4FinalStateChargeDecrease::G4FinalStateChargeDecrease()
34{
35 lowEnergyLimit = 1 * keV;
36 highEnergyLimit = 10 * MeV;
37
38 G4cout << G4endl;
39 G4cout << "*******************************************************************************" << G4endl;
40 G4cout << "*******************************************************************************" << G4endl;
41 G4cout << " The class G4FinalStateChargeDecrease is NOT SUPPORTED ANYMORE. " << G4endl;
42 G4cout << " It will be REMOVED with the next major release of Geant4. " << G4endl;
43 G4cout << " Please consult: https://twiki.cern.ch/twiki/bin/view/Geant4/LoweProcesses" << G4endl;
44 G4cout << "*******************************************************************************" << G4endl;
45 G4cout << "*******************************************************************************" << G4endl;
46 G4cout << G4endl;
47}
48
49//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
50
51G4FinalStateChargeDecrease::~G4FinalStateChargeDecrease()
52{}
53
54//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
55
56const G4FinalStateProduct& G4FinalStateChargeDecrease::GenerateFinalState(const G4Track& track, const G4Step& /* step */)
57{
58 product.Clear();
59
60 G4double inK = track.GetDynamicParticle()->GetKineticEnergy();
61
62 G4ParticleDefinition* definition = track.GetDefinition();
63 G4int finalStateIndex = cross.RandomSelect(inK,definition);
64
65 G4int n = NumberOfFinalStates(definition, finalStateIndex);
66 G4double waterBindingEnergy = WaterBindingEnergyConstant(definition, finalStateIndex);
67 G4double outgoingParticleBindingEnergy = OutgoingParticleBindingEnergyConstant(definition, finalStateIndex);
68
69 G4double outK = 0.;
70 if (definition==G4Proton::Proton())
71 outK = inK - n*(inK*electron_mass_c2/proton_mass_c2) - waterBindingEnergy + outgoingParticleBindingEnergy;
72 else
73 outK = inK - n*(inK*electron_mass_c2/(3728*MeV)) - waterBindingEnergy + outgoingParticleBindingEnergy;
74
75 if (outK<0)
76 {
77 G4String message;
78 message="ChargeDecreaseDingfelder::GenerateFinalState - Final kinetic energy is below 0! Process ";
79 G4Exception(message);
80 }
81
82 //SI - Added protection against total energy deposit
83 product.DoNotDepositEnergy();
84 //
85 product.KillPrimaryParticle();
86
87 product.AddEnergyDeposit(waterBindingEnergy);
88
89 G4DynamicParticle* aSecondary = new G4DynamicParticle(OutgoingParticleDefinition(definition, finalStateIndex),
90 track.GetDynamicParticle()->GetMomentumDirection(),
91 outK);
92
93 product.AddSecondary(aSecondary);
94
95 return product;
96}
97
98//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
99
100G4int G4FinalStateChargeDecrease::NumberOfFinalStates(G4ParticleDefinition* particleDefinition,
101 G4int finalStateIndex )
102
103{
104 if (particleDefinition == G4Proton::Proton()) return 1;
105
106 G4DNAGenericIonsManager*instance;
107 instance = G4DNAGenericIonsManager::Instance();
108
109 if (particleDefinition == instance->GetIon("alpha++") )
110 {
111 if (finalStateIndex==0) return 1;
112 return 2;
113 }
114
115 if (particleDefinition == instance->GetIon("alpha+") ) return 1;
116
117 return 0;
118}
119
120//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
121
122G4ParticleDefinition* G4FinalStateChargeDecrease::OutgoingParticleDefinition (G4ParticleDefinition* particleDefinition,
123 G4int finalStateIndex)
124{
125 G4DNAGenericIonsManager * instance(G4DNAGenericIonsManager::Instance());
126
127 if (particleDefinition == G4Proton::Proton()) return instance->GetIon("hydrogen");
128
129 if (particleDefinition == instance->GetIon("alpha++") )
130 {
131 if (finalStateIndex == 0) return instance->GetIon("alpha+");
132 return instance->GetIon("helium");
133 }
134
135 if (particleDefinition == instance->GetIon("alpha+") ) return instance->GetIon("helium");
136
137 return 0;
138}
139
140//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
141
142G4double G4FinalStateChargeDecrease::WaterBindingEnergyConstant(G4ParticleDefinition* particleDefinition,
143 G4int finalStateIndex )
144{
145 // Ionization energy of first water shell
146 // Rad. Phys. Chem. 59 p.267 by Dingf. et al.
147 // W + 10.79 eV -> W+ + e-
148
149 G4DNAGenericIonsManager * instance(G4DNAGenericIonsManager::Instance());
150
151 if(particleDefinition == G4Proton::Proton()) return 10.79*eV;
152
153 if (particleDefinition == instance->GetIon("alpha++") )
154 {
155 // Binding energy for W+ -> W++ + e- 10.79 eV
156 // Binding energy for W -> W+ + e- 10.79 eV
157 //
158 // Ionization energy of first water shell
159 // Rad. Phys. Chem. 59 p.267 by Dingf. et al.
160
161 if (finalStateIndex == 0) return 10.79*eV;
162
163 return 10.79*2*eV;
164 }
165
166 if (particleDefinition == instance->GetIon("alpha+") )
167 {
168 // Binding energy for W+ -> W++ + e- 10.79 eV
169 // Binding energy for W -> W+ + e- 10.79 eV
170 //
171 // Ionization energy of first water shell
172 // Rad. Phys. Chem. 59 p.267 by Dingf. et al.
173
174 return 10.79*eV;
175 }
176
177 return 0;
178}
179
180//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
181
182G4double G4FinalStateChargeDecrease::OutgoingParticleBindingEnergyConstant(G4ParticleDefinition* particleDefinition,
183 G4int finalStateIndex)
184{
185 G4DNAGenericIonsManager * instance(G4DNAGenericIonsManager::Instance());
186
187 if(particleDefinition == G4Proton::Proton()) return 13.6*eV;
188
189 if (particleDefinition == instance->GetIon("alpha++") )
190 {
191 // Binding energy for He+ -> He++ + e- 54.509 eV
192 // Binding energy for He -> He+ + e- 24.587 eV
193
194 if (finalStateIndex==0) return 54.509*eV;
195
196 return (54.509 + 24.587)*eV;
197 }
198
199 if (particleDefinition == instance->GetIon("alpha+") )
200 {
201 // Binding energy for He+ -> He++ + e- 54.509 eV
202 // Binding energy for He -> He+ + e- 24.587 eV
203
204 return 24.587*eV;
205 }
206
207 return 0;
208}
209
Note: See TracBrowser for help on using the repository browser.