Geant4  10.00.p03
G4MuElecInelastic.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 // G4MuElecInelastic.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 "G4MuElecInelastic.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  G4cout << G4endl;
58  G4cout << "*******************************************************************************" << G4endl;
59  G4cout << "*******************************************************************************" << G4endl;
60  G4cout << " The name of the class G4MuElecInelastic is changed to G4MicroElecInelastic. " << G4endl;
61  G4cout << " The obsolete class will be REMOVED with the next release of Geant4. " << G4endl;
62  G4cout << "*******************************************************************************" << G4endl;
63  G4cout << "*******************************************************************************" << G4endl;
64  G4cout << G4endl;
65 }
66 
67 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
68 
70 {}
71 
72 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
73 
75 {
76  return (&p == G4Electron::Electron() ||
77  &p == G4Proton::Proton() ||
78  (p.GetPDGCharge() != 0.0 && !p.IsShortLived() && p.GetParticleType() == "nucleus"));
79 }
80 
81 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
82 
84 {
85  if(!isInitialised)
86  {
87  isInitialised = true;
88  SetBuildTableFlag(false);
90 
91  if(name == "e-")
92  {
94  EmModel()->SetLowEnergyLimit(16.7*eV);
96 
97  AddEmModel(1, EmModel());
98  }
99 
100  else if(name == "proton")
101  {
103  EmModel()->SetLowEnergyLimit(50.*keV);
105 
106  AddEmModel(1, EmModel());
107  }
108 
109  else if(name == "GenericIon")
110  {
112  EmModel()->SetLowEnergyLimit(50.*keV);
113  EmModel()->SetHighEnergyLimit(10000.*GeV);
114 
115  AddEmModel(1, EmModel());
116  }
117  }
118 }
119 
120 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
121 
123 {
124  // V.I. printout of models is perfored by model manager
125  // if this extra printout is needed it should be
126  // protected by verbosity level
127  /*
128  if (EmModel(2))
129  {
130  G4cout
131  << " Total cross sections computed from "
132  << EmModel(1)->GetName()
133  << " and "
134  << EmModel(2)->GetName()
135  << " models"
136  << G4endl;
137  }
138  else
139  {
140  G4cout
141  << " Total cross sections computed from "
142  << EmModel()->GetName()
143  << G4endl;
144  }
145  */
146 }
147 
148 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
static const double MeV
Definition: G4SIunits.hh:193
virtual void PrintInfo()
virtual void InitialiseProcess(const G4ParticleDefinition *)
void SetBuildTableFlag(G4bool val)
G4VEmModel * EmModel(G4int index=1) const
G4String name
Definition: TRTMaterials.hh:40
const G4String & GetParticleName() const
virtual G4bool IsApplicable(const G4ParticleDefinition &)
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 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)
static G4Electron * Electron()
Definition: G4Electron.cc:94
#define G4endl
Definition: G4ios.hh:61
static const double keV
Definition: G4SIunits.hh:195
G4MuElecInelastic(const G4String &processName="MuElecIonisation", G4ProcessType type=fElectromagnetic)
virtual ~G4MuElecInelastic()
void SetLowEnergyLimit(G4double)
Definition: G4VEmModel.hh:690
G4double GetPDGCharge() const
G4ProcessType