Geant4  10.00.p03
G4MicroElecInelastic.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 // G4MicroElecInelastic.cc, 2011/08/29 A.Valentin, M. Raine
28 //
29 // Based on the following publications
30 //
31 // - Inelastic cross-sections of low energy electrons in silicon
32 // for the simulation of heavy ion tracks with theGeant4-DNA toolkit,
33 // NSS Conf. Record 2010, pp. 80-85.
34 // - Geant4 physics processes for microdosimetry simulation:
35 // very low energy electromagnetic models for electrons in Si,
36 // NIM B, vol. 288, pp. 66 - 73, 2012.
37 // - Geant4 physics processes for microdosimetry simulation:
38 // very low energy electromagnetic models for protons and
39 // heavy ions in Si, NIM B, vol. 287, pp. 124 - 129, 2012.
40 //
41 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
42 
43 
44 #include "G4MicroElecInelastic.hh"
45 #include "G4SystemOfUnits.hh"
46 
47 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
48 
49 using namespace std;
50 
52  G4ProcessType type):G4VEmProcess (processName, type),
53  isInitialised(false)
54 {
56 }
57 
58 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
59 
61 {}
62 
63 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
64 
66 {
67  return (&p == G4Electron::Electron() ||
68  &p == G4Proton::Proton() ||
69  (p.GetPDGCharge() != 0.0 && !p.IsShortLived() && p.GetParticleType() == "nucleus"));
70 }
71 
72 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
73 
75 {
76  if(!isInitialised)
77  {
78  isInitialised = true;
79  SetBuildTableFlag(false);
81 
82  if(name == "e-")
83  {
85  EmModel()->SetLowEnergyLimit(16.7*eV);
87 
88  AddEmModel(1, EmModel());
89  }
90 
91  else if(name == "proton")
92  {
96 
97  AddEmModel(1, EmModel());
98  }
99 
100  else if(name == "GenericIon")
101  {
103  EmModel()->SetLowEnergyLimit(50.*keV);
104  EmModel()->SetHighEnergyLimit(10000.*GeV);
105 
106  AddEmModel(1, EmModel());
107  }
108  }
109 }
110 
111 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
112 
114 {
115  // V.I. printout of models is perfored by model manager
116  // if this extra printout is needed it should be
117  // protected by verbosity level
118  /*
119  if (EmModel(2))
120  {
121  G4cout
122  << " Total cross sections computed from "
123  << EmModel(1)->GetName()
124  << " and "
125  << EmModel(2)->GetName()
126  << " models"
127  << G4endl;
128  }
129  else
130  {
131  G4cout
132  << " Total cross sections computed from "
133  << EmModel()->GetName()
134  << G4endl;
135  }
136  */
137 }
138 
139 //....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
virtual void InitialiseProcess(const G4ParticleDefinition *)
const G4String & GetParticleName() const
void SetHighEnergyLimit(G4double)
Definition: G4VEmModel.hh:683
void SetEmModel(G4VEmModel *, G4int index=1)
G4MicroElecInelastic(const G4String &processName="MicroElecIonisation", G4ProcessType type=fElectromagnetic)
bool G4bool
Definition: G4Types.hh:79
static G4Proton * Proton()
Definition: G4Proton.cc:93
void SetProcessSubType(G4int)
Definition: G4VProcess.hh:432
static const double GeV
Definition: G4SIunits.hh:196
const G4String & GetParticleType() const
static const double eV
Definition: G4SIunits.hh:194
void AddEmModel(G4int, G4VEmModel *, const G4Region *region=0)
virtual G4bool IsApplicable(const G4ParticleDefinition &)
static G4Electron * Electron()
Definition: G4Electron.cc:94
static const double keV
Definition: G4SIunits.hh:195
void SetLowEnergyLimit(G4double)
Definition: G4VEmModel.hh:690
G4double GetPDGCharge() const
G4ProcessType