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" << 46 47 #include "G4Proton.hh" << 48 #include "G4GenericIon.hh" 47 #include "G4GenericIon.hh" 49 #include "G4Alpha.hh" << 50 #include "G4LowEnergyEmProcessSubType.hh" << 51 48 52 //....oooOO0OOooo........oooOO0OOooo........oo 49 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 53 50 54 using namespace std; 51 using namespace std; 55 52 56 G4MicroElecInelastic::G4MicroElecInelastic(con 53 G4MicroElecInelastic::G4MicroElecInelastic(const G4String& processName, 57 G4P << 54 G4ProcessType type):G4VEmProcess (processName, type), 58 : G4VEmProcess (processName, type), isInitia << 55 isInitialised(false) 59 { 56 { 60 SetProcessSubType(fLowEnergyIonisation); << 57 SetProcessSubType(53); 61 } 58 } 62 59 63 //....oooOO0OOooo........oooOO0OOooo........oo 60 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 64 61 65 G4MicroElecInelastic::~G4MicroElecInelastic() 62 G4MicroElecInelastic::~G4MicroElecInelastic() 66 {} 63 {} 67 64 68 //....oooOO0OOooo........oooOO0OOooo........oo 65 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 69 66 70 G4bool G4MicroElecInelastic::IsApplicable(cons 67 G4bool G4MicroElecInelastic::IsApplicable(const G4ParticleDefinition& p) 71 { 68 { 72 return (&p == G4Electron::Electron() || 69 return (&p == G4Electron::Electron() || 73 &p == G4Proton::Proton() || 70 &p == G4Proton::Proton() || 74 &p == G4Alpha::Alpha() || << 75 &p == G4GenericIon::GenericIonDefini 71 &p == G4GenericIon::GenericIonDefinition()); 76 } 72 } 77 73 78 //....oooOO0OOooo........oooOO0OOooo........oo 74 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 79 75 80 void G4MicroElecInelastic::InitialiseProcess(c 76 void G4MicroElecInelastic::InitialiseProcess(const G4ParticleDefinition* p) 81 { 77 { 82 if(!isInitialised) 78 if(!isInitialised) 83 { 79 { 84 isInitialised = true; 80 isInitialised = true; 85 SetBuildTableFlag(false); 81 SetBuildTableFlag(false); 86 G4String name = p->GetParticleName(); 82 G4String name = p->GetParticleName(); 87 83 88 if(!EmModel()) SetEmModel(new G4DummyModel << 84 if(name == "e-") 89 AddEmModel(2, EmModel()); << 85 { >> 86 if(!EmModel()) SetEmModel(new G4MicroElecInelasticModel); >> 87 EmModel()->SetLowEnergyLimit(16.7*eV); >> 88 EmModel()->SetHighEnergyLimit(100*MeV); >> 89 >> 90 AddEmModel(1, EmModel()); >> 91 } >> 92 >> 93 else if(name == "proton") >> 94 { >> 95 if(!EmModel()) SetEmModel(new G4MicroElecInelasticModel); >> 96 EmModel()->SetLowEnergyLimit(50.*keV); >> 97 EmModel()->SetHighEnergyLimit(10*GeV); >> 98 >> 99 AddEmModel(1, EmModel()); >> 100 } >> 101 >> 102 else if(name == "GenericIon") >> 103 { >> 104 if(!EmModel()) SetEmModel(new G4MicroElecInelasticModel); >> 105 EmModel()->SetLowEnergyLimit(50.*keV); >> 106 EmModel()->SetHighEnergyLimit(p->GetAtomicMass()*10.*GeV); >> 107 >> 108 AddEmModel(1, EmModel()); >> 109 } 90 } 110 } 91 } 111 } >> 112 >> 113 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... >> 114 >> 115 void G4MicroElecInelastic::PrintInfo() >> 116 { >> 117 // V.I. printout of models is perfored by model manager >> 118 // if this extra printout is needed it should be >> 119 // protected by verbosity level >> 120 /* >> 121 if (EmModel(2)) >> 122 { >> 123 G4cout >> 124 << " Total cross sections computed from " >> 125 << EmModel(1)->GetName() >> 126 << " and " >> 127 << EmModel(2)->GetName() >> 128 << " models" >> 129 << G4endl; >> 130 } >> 131 else >> 132 { >> 133 G4cout >> 134 << " Total cross sections computed from " >> 135 << EmModel()->GetName() >> 136 << G4endl; >> 137 } >> 138 */ >> 139 } 92 140 93 //....oooOO0OOooo........oooOO0OOooo........oo 141 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 94 142