Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer 3 // * License and Disclaimer * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. 9 // * include a list of copyright holders. * 10 // * 10 // * * 11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. * 17 // * 17 // * * 18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************* 24 // ******************************************************************** 25 // 25 // 26 // 26 // 27 // G4MicroElecInelastic.cc, 2011/08/29 A.Valen 27 // G4MicroElecInelastic.cc, 2011/08/29 A.Valentin, M. Raine 28 // 28 // 29 // Based on the following publications 29 // Based on the following publications 30 // 30 // 31 // - Inelastic cross-sections of low 31 // - Inelastic cross-sections of low energy electrons in silicon 32 // for the simulation of heavy ion tracks 32 // for the simulation of heavy ion tracks with theGeant4-DNA toolkit, 33 // NSS Conf. Record 2010, pp. 80-85. 33 // NSS Conf. Record 2010, pp. 80-85. 34 // - Geant4 physics processes for microdo 34 // - Geant4 physics processes for microdosimetry simulation: 35 // very low energy electromagnetic models 35 // very low energy electromagnetic models for electrons in Si, 36 // NIM B, vol. 288, pp. 66 - 73, 2012. 36 // NIM B, vol. 288, pp. 66 - 73, 2012. 37 // - Geant4 physics processes for microdo 37 // - Geant4 physics processes for microdosimetry simulation: 38 // very low energy electromagnetic models 38 // very low energy electromagnetic models for protons and 39 // heavy ions in Si, NIM B, vol. 287, pp. 39 // heavy ions in Si, NIM B, vol. 287, pp. 124 - 129, 2012. 40 // 40 // 41 //....oooOO0OOooo........oooOO0OOooo........oo 41 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 42 42 >> 43 43 #include "G4MicroElecInelastic.hh" 44 #include "G4MicroElecInelastic.hh" 44 #include "G4DummyModel.hh" << 45 #include "G4SystemOfUnits.hh" 45 #include "G4SystemOfUnits.hh" 46 #include "G4Electron.hh" << 47 #include "G4Proton.hh" << 48 #include "G4GenericIon.hh" << 49 #include "G4Alpha.hh" << 50 #include "G4LowEnergyEmProcessSubType.hh" << 51 46 52 //....oooOO0OOooo........oooOO0OOooo........oo 47 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 53 48 54 using namespace std; 49 using namespace std; 55 50 56 G4MicroElecInelastic::G4MicroElecInelastic(con 51 G4MicroElecInelastic::G4MicroElecInelastic(const G4String& processName, 57 G4P << 52 G4ProcessType type):G4VEmProcess (processName, type), 58 : G4VEmProcess (processName, type), isInitia << 53 isInitialised(false) 59 { 54 { 60 SetProcessSubType(fLowEnergyIonisation); << 55 SetProcessSubType(53); 61 } 56 } 62 57 63 //....oooOO0OOooo........oooOO0OOooo........oo 58 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 64 59 65 G4MicroElecInelastic::~G4MicroElecInelastic() 60 G4MicroElecInelastic::~G4MicroElecInelastic() 66 {} 61 {} 67 62 68 //....oooOO0OOooo........oooOO0OOooo........oo 63 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 69 64 70 G4bool G4MicroElecInelastic::IsApplicable(cons 65 G4bool G4MicroElecInelastic::IsApplicable(const G4ParticleDefinition& p) 71 { 66 { 72 return (&p == G4Electron::Electron() || 67 return (&p == G4Electron::Electron() || 73 &p == G4Proton::Proton() || << 68 &p == G4Proton::Proton() || 74 &p == G4Alpha::Alpha() || << 69 (p.GetPDGCharge() != 0.0 && !p.IsShortLived() && p.GetParticleType() == "nucleus")); 75 &p == G4GenericIon::GenericIonDefini << 76 } 70 } 77 71 78 //....oooOO0OOooo........oooOO0OOooo........oo 72 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 79 73 80 void G4MicroElecInelastic::InitialiseProcess(c 74 void G4MicroElecInelastic::InitialiseProcess(const G4ParticleDefinition* p) 81 { 75 { 82 if(!isInitialised) 76 if(!isInitialised) 83 { 77 { 84 isInitialised = true; 78 isInitialised = true; 85 SetBuildTableFlag(false); 79 SetBuildTableFlag(false); 86 G4String name = p->GetParticleName(); 80 G4String name = p->GetParticleName(); 87 81 88 if(!EmModel()) SetEmModel(new G4DummyModel << 82 if(name == "e-") 89 AddEmModel(2, EmModel()); << 83 { >> 84 if(!EmModel()) SetEmModel(new G4MicroElecInelasticModel); >> 85 EmModel()->SetLowEnergyLimit(16.7*eV); >> 86 EmModel()->SetHighEnergyLimit(100*MeV); >> 87 >> 88 AddEmModel(1, EmModel()); >> 89 } >> 90 >> 91 else if(name == "proton") >> 92 { >> 93 if(!EmModel()) SetEmModel(new G4MicroElecInelasticModel); >> 94 EmModel()->SetLowEnergyLimit(50.*keV); >> 95 EmModel()->SetHighEnergyLimit(10*GeV); >> 96 >> 97 AddEmModel(1, EmModel()); >> 98 } >> 99 >> 100 else if(name == "GenericIon") >> 101 { >> 102 if(!EmModel()) SetEmModel(new G4MicroElecInelasticModel); >> 103 EmModel()->SetLowEnergyLimit(50.*keV); >> 104 EmModel()->SetHighEnergyLimit(10000.*GeV); >> 105 >> 106 AddEmModel(1, EmModel()); >> 107 } 90 } 108 } 91 } 109 } >> 110 >> 111 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... >> 112 >> 113 void G4MicroElecInelastic::PrintInfo() >> 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 } 92 138 93 //....oooOO0OOooo........oooOO0OOooo........oo 139 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 94 140