Geant4  10.00.p02
G4DNAIonisation.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 // $Id: G4DNAIonisation.cc 70171 2013-05-24 13:34:18Z gcosmo $
27 
28 #include "G4DNAIonisation.hh"
29 #include "G4SystemOfUnits.hh"
30 
31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
32 
33 using namespace std;
34 
36  G4ProcessType type):G4VEmProcess (processName, type),
37  isInitialised(false)
38 {
40 }
41 
42 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
43 
45 {}
46 
47 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
48 
50 {
53 
54  return
55  (
56  &p == G4Electron::Electron()
57  || &p == G4Proton::Proton()
58  || &p == instance->GetIon("hydrogen")
59  || &p == instance->GetIon("alpha++")
60  || &p == instance->GetIon("alpha+")
61  || &p == instance->GetIon("helium")
62  || &p == instance->GetIon("carbon")
63  || &p == instance->GetIon("nitrogen")
64  || &p == instance->GetIon("oxygen")
65  || &p == instance->GetIon("iron")
66  );
67 }
68 
69 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
70 
72 {
73  if(!isInitialised)
74  {
75  isInitialised = true;
76  SetBuildTableFlag(false);
77 
79 
80  if(name == "e-")
81  {
85 
86  AddEmModel(1, EmModel());
87  }
88 
89  if(name == "proton")
90  {
94 
96  EmModel(2)->SetLowEnergyLimit(500*keV);
98 
99  AddEmModel(1, EmModel(1));
100  AddEmModel(2, EmModel(2));
101  }
102 
103  if(name == "hydrogen")
104  {
108 
109  AddEmModel(1, EmModel());
110  }
111 
112  if(name == "alpha" || name == "alpha+" || name == "helium" )
113  {
117 
118  AddEmModel(1, EmModel());
119  }
120 
121  // Extension to HZE proposed by Z. Francis
122 
123  if(name == "carbon" || name == "nitrogen" || name == "oxygen" || name == "iron")
124  {
128 
129  AddEmModel(1, EmModel());
130  }
131 
132  }
133 }
134 
135 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
136 
138 {
139  if (EmModel(2))
140  {
141  G4cout
142  << " Total cross sections computed from "
143  << EmModel(1)->GetName()
144  << " and "
145  << EmModel(2)->GetName()
146  << " models"
147  << G4endl;
148  }
149  else
150  {
151  G4cout
152  << " Total cross sections computed from "
153  << EmModel()->GetName()
154  << G4endl;
155  }
156 }
157 
158 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
static const double MeV
Definition: G4SIunits.hh:193
void SetBuildTableFlag(G4bool val)
G4VEmModel * EmModel(G4int index=1) const
G4String name
Definition: TRTMaterials.hh:40
const G4String & GetParticleName() const
void SetHighEnergyLimit(G4double)
Definition: G4VEmModel.hh:683
void SetEmModel(G4VEmModel *, G4int index=1)
G4GLOB_DLL std::ostream G4cout
bool G4bool
Definition: G4Types.hh:79
static G4Proton * Proton()
Definition: G4Proton.cc:93
void SetProcessSubType(G4int)
Definition: G4VProcess.hh:432
static G4DNAGenericIonsManager * Instance(void)
G4int GetAtomicMass() const
static const double eV
Definition: G4SIunits.hh:194
virtual G4bool IsApplicable(const G4ParticleDefinition &)
void AddEmModel(G4int, G4VEmModel *, const G4Region *region=0)
virtual void PrintInfo()
static G4Electron * Electron()
Definition: G4Electron.cc:94
#define G4endl
Definition: G4ios.hh:61
const G4String & GetName() const
Definition: G4VEmModel.hh:753
static MCTruthManager * instance
static const double keV
Definition: G4SIunits.hh:195
virtual void InitialiseProcess(const G4ParticleDefinition *)
void SetLowEnergyLimit(G4double)
Definition: G4VEmModel.hh:690
virtual ~G4DNAIonisation()
G4DNAIonisation(const G4String &processName="DNAIonisation", G4ProcessType type=fElectromagnetic)
G4ParticleDefinition * GetIon(const G4String &name)
G4ProcessType