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

Last change on this file since 968 was 961, checked in by garnier, 17 years ago

update processes
  • Property svn:executable set to *
File size: 7.0 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.4 2009/01/20 07:50:28 sincerti Exp $
27// GEANT4 tag $Name: geant4-09-02-ref-02 $
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
39//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
40
41G4FinalStateChargeDecrease::~G4FinalStateChargeDecrease()
42{}
43
44//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
45
46const G4FinalStateProduct& G4FinalStateChargeDecrease::GenerateFinalState(const G4Track& track, const G4Step& /* step */)
47{
48 product.Clear();
49
50 G4double inK = track.GetDynamicParticle()->GetKineticEnergy();
51
52 G4ParticleDefinition* definition = track.GetDefinition();
53 G4int finalStateIndex = cross.RandomSelect(inK,definition);
54
55 G4int n = NumberOfFinalStates(definition, finalStateIndex);
56 G4double waterBindingEnergy = WaterBindingEnergyConstant(definition, finalStateIndex);
57 G4double outgoingParticleBindingEnergy = OutgoingParticleBindingEnergyConstant(definition, finalStateIndex);
58
59 G4double outK = 0.;
60 if (definition==G4Proton::Proton())
61 outK = inK - n*(inK*electron_mass_c2/proton_mass_c2) - waterBindingEnergy + outgoingParticleBindingEnergy;
62 else
63 outK = inK - n*(inK*electron_mass_c2/(3728*MeV)) - waterBindingEnergy + outgoingParticleBindingEnergy;
64
65 if (outK<0)
66 {
67 G4String message;
68 message="ChargeDecreaseDingfelder::GenerateFinalState - Final kinetic energy is below 0! Process ";
69 G4Exception(message);
70 }
71
72 //SI - Added protection against total energy deposit
73 product.DoNotDepositEnergy();
74 //
75 product.KillPrimaryParticle();
76
77 product.AddEnergyDeposit(waterBindingEnergy);
78
79 G4DynamicParticle* aSecondary = new G4DynamicParticle(OutgoingParticleDefinition(definition, finalStateIndex),
80 track.GetDynamicParticle()->GetMomentumDirection(),
81 outK);
82
83 product.AddSecondary(aSecondary);
84
85 return product;
86}
87
88//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
89
90G4int G4FinalStateChargeDecrease::NumberOfFinalStates(G4ParticleDefinition* particleDefinition,
91 G4int finalStateIndex )
92
93{
94 if (particleDefinition == G4Proton::Proton()) return 1;
95
96 G4DNAGenericIonsManager*instance;
97 instance = G4DNAGenericIonsManager::Instance();
98
99 if (particleDefinition == instance->GetIon("alpha++") )
100 {
101 if (finalStateIndex==0) return 1;
102 return 2;
103 }
104
105 if (particleDefinition == instance->GetIon("alpha+") ) return 1;
106
107 return 0;
108}
109
110//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
111
112G4ParticleDefinition* G4FinalStateChargeDecrease::OutgoingParticleDefinition (G4ParticleDefinition* particleDefinition,
113 G4int finalStateIndex)
114{
115 G4DNAGenericIonsManager * instance(G4DNAGenericIonsManager::Instance());
116
117 if (particleDefinition == G4Proton::Proton()) return instance->GetIon("hydrogen");
118
119 if (particleDefinition == instance->GetIon("alpha++") )
120 {
121 if (finalStateIndex == 0) return instance->GetIon("alpha+");
122 return instance->GetIon("helium");
123 }
124
125 if (particleDefinition == instance->GetIon("alpha+") ) return instance->GetIon("helium");
126
127 return 0;
128}
129
130//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
131
132G4double G4FinalStateChargeDecrease::WaterBindingEnergyConstant(G4ParticleDefinition* particleDefinition,
133 G4int finalStateIndex )
134{
135 // Ionization energy of first water shell
136 // Rad. Phys. Chem. 59 p.267 by Dingf. et al.
137 // W + 10.79 eV -> W+ + e-
138
139 G4DNAGenericIonsManager * instance(G4DNAGenericIonsManager::Instance());
140
141 if(particleDefinition == G4Proton::Proton()) return 10.79*eV;
142
143 if (particleDefinition == instance->GetIon("alpha++") )
144 {
145 // Binding energy for W+ -> W++ + e- 10.79 eV
146 // Binding energy for W -> W+ + e- 10.79 eV
147 //
148 // Ionization energy of first water shell
149 // Rad. Phys. Chem. 59 p.267 by Dingf. et al.
150
151 if (finalStateIndex == 0) return 10.79*eV;
152
153 return 10.79*2*eV;
154 }
155
156 if (particleDefinition == instance->GetIon("alpha+") )
157 {
158 // Binding energy for W+ -> W++ + e- 10.79 eV
159 // Binding energy for W -> W+ + e- 10.79 eV
160 //
161 // Ionization energy of first water shell
162 // Rad. Phys. Chem. 59 p.267 by Dingf. et al.
163
164 return 10.79*eV;
165 }
166
167 return 0;
168}
169
170//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
171
172G4double G4FinalStateChargeDecrease::OutgoingParticleBindingEnergyConstant(G4ParticleDefinition* particleDefinition,
173 G4int finalStateIndex)
174{
175 G4DNAGenericIonsManager * instance(G4DNAGenericIonsManager::Instance());
176
177 if(particleDefinition == G4Proton::Proton()) return 13.6*eV;
178
179 if (particleDefinition == instance->GetIon("alpha++") )
180 {
181 // Binding energy for He+ -> He++ + e- 54.509 eV
182 // Binding energy for He -> He+ + e- 24.587 eV
183
184 if (finalStateIndex==0) return 54.509*eV;
185
186 return (54.509 + 24.587)*eV;
187 }
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 return 24.587*eV;
195 }
196
197 return 0;
198}
199
Note: See TracBrowser for help on using the repository browser.