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 // Authors: S. Guatelli and M. G. Pia, INFN Genova, Italy 27 // Previous authors: G. Guerrieri and M. G. Pi << 28 // Authors (since 2007): S. Guatelli, Universi << 29 // << 30 // 27 // 31 // << 28 // Based on code developed by the undergraduate student G. Guerrieri 32 // ********************************** << 29 // Note: this is a preliminary beta-version of the code; an improved 33 // * * << 30 // version will be distributed in the next Geant4 public release, compliant 34 // * G4HumanPhantomPhysicsList.cc << 31 // with the design in a forthcoming publication, and subject to a 35 // * * << 32 // design and code review. 36 // ********************************** << 33 37 // << 38 // << 39 #include "G4HumanPhantomPhysicsList.hh" 34 #include "G4HumanPhantomPhysicsList.hh" 40 35 41 #include "G4SystemOfUnits.hh" << 42 #include "G4ParticleDefinition.hh" 36 #include "G4ParticleDefinition.hh" 43 #include "G4ProductionCutsTable.hh" 37 #include "G4ProductionCutsTable.hh" 44 #include "G4ProcessManager.hh" 38 #include "G4ProcessManager.hh" 45 #include "G4ParticleTypes.hh" 39 #include "G4ParticleTypes.hh" 46 #include "G4UnitsTable.hh" 40 #include "G4UnitsTable.hh" 47 #include "G4ios.hh" 41 #include "G4ios.hh" 48 #include "G4VModularPhysicsList.hh" << 49 #include "G4EmStandardPhysics_option4.hh" << 50 #include "G4DecayPhysics.hh" << 51 42 52 G4HumanPhantomPhysicsList::G4HumanPhantomPhysi << 43 #include "G4MultipleScattering.hh" >> 44 // gamma >> 45 #include "G4LowEnergyRayleigh.hh" >> 46 #include "G4LowEnergyPhotoElectric.hh" >> 47 #include "G4LowEnergyCompton.hh" >> 48 #include "G4LowEnergyGammaConversion.hh" >> 49 // e- >> 50 #include "G4LowEnergyIonisation.hh" >> 51 #include "G4LowEnergyBremsstrahlung.hh" >> 52 // e+ >> 53 #include "G4eIonisation.hh" >> 54 #include "G4eBremsstrahlung.hh" >> 55 #include "G4eplusAnnihilation.hh" >> 56 >> 57 G4HumanPhantomPhysicsList::G4HumanPhantomPhysicsList(): G4VUserPhysicsList() 53 { 58 { 54 SetVerboseLevel(1); 59 SetVerboseLevel(1); 55 emPhysicsList = new G4EmStandardPhysics_opti << 56 decPhysicsList = new G4DecayPhysics(); << 57 // Alternatively you can substitute this physi << 58 // with the LowEnergy Livermore or LowEnergy P << 59 // emPhysicsList = new G4EmLivermorePhysics(); << 60 // Low Energy based on Livermore Evaluated Dat << 61 // << 62 // Penelope physics << 63 //emPhysicsList = new G4EmPenelopePhysics(); << 64 } 60 } 65 61 66 G4HumanPhantomPhysicsList::~G4HumanPhantomPhys 62 G4HumanPhantomPhysicsList::~G4HumanPhantomPhysicsList() 67 { 63 { 68 delete decPhysicsList; << 69 delete emPhysicsList; << 70 } 64 } 71 65 72 void G4HumanPhantomPhysicsList::ConstructParti 66 void G4HumanPhantomPhysicsList::ConstructParticle() 73 { 67 { 74 decPhysicsList -> ConstructParticle(); << 68 // In this method, static member functions should be called >> 69 // for all particles which you want to use. >> 70 // This ensures that objects of these particle types will be >> 71 // created in the program. >> 72 G4Geantino::GeantinoDefinition(); >> 73 >> 74 ConstructBosons(); >> 75 ConstructLeptons(); >> 76 } >> 77 >> 78 void G4HumanPhantomPhysicsList::ConstructBosons() >> 79 { >> 80 // photons >> 81 G4Gamma::GammaDefinition(); >> 82 } >> 83 >> 84 void G4HumanPhantomPhysicsList::ConstructLeptons() >> 85 { >> 86 // leptons >> 87 G4Electron::ElectronDefinition(); >> 88 G4Positron::PositronDefinition(); 75 } 89 } 76 90 77 void G4HumanPhantomPhysicsList::ConstructProce 91 void G4HumanPhantomPhysicsList::ConstructProcess() 78 { 92 { 79 AddTransportation(); 93 AddTransportation(); 80 emPhysicsList -> ConstructProcess(); << 94 ConstructEM(); 81 decPhysicsList -> ConstructProcess(); << 95 } >> 96 >> 97 void G4HumanPhantomPhysicsList::ConstructEM() >> 98 { >> 99 theParticleIterator->reset(); >> 100 >> 101 while( (*theParticleIterator)() ){ >> 102 >> 103 G4ParticleDefinition* particle = theParticleIterator->value(); >> 104 G4ProcessManager* pmanager = particle->GetProcessManager(); >> 105 G4String particleName = particle->GetParticleName(); >> 106 >> 107 // Processes >> 108 >> 109 if (particleName == "gamma") { >> 110 // Photon >> 111 pmanager->AddDiscreteProcess(new G4LowEnergyRayleigh); >> 112 pmanager->AddDiscreteProcess(new G4LowEnergyPhotoElectric); >> 113 pmanager->AddDiscreteProcess(new G4LowEnergyCompton); >> 114 pmanager->AddDiscreteProcess(new G4LowEnergyGammaConversion); >> 115 >> 116 } else if (particleName == "e-") { >> 117 // Electron >> 118 G4LowEnergyIonisation* loweIon = new G4LowEnergyIonisation("LowEnergyIoni"); >> 119 >> 120 G4LowEnergyBremsstrahlung* loweBrem = new G4LowEnergyBremsstrahlung("LowEnBrem"); >> 121 // Select the Bremsstrahlung angular distribution model (Tsai/2BN/2BS) >> 122 loweBrem->SetAngularGenerator("tsai"); >> 123 >> 124 pmanager->AddProcess(new G4MultipleScattering, -1, 1,1); >> 125 pmanager->AddProcess(loweIon, -1, 2,2); >> 126 pmanager->AddProcess(loweBrem, -1,-1,3); >> 127 >> 128 } else if (particleName == "e+") { >> 129 // Positron >> 130 pmanager->AddProcess(new G4MultipleScattering, -1, 1,1); >> 131 pmanager->AddProcess(new G4eIonisation, -1, 2,2); >> 132 pmanager->AddProcess(new G4eBremsstrahlung, -1,-1,3); >> 133 pmanager->AddProcess(new G4eplusAnnihilation, 0,-1,4); >> 134 >> 135 } >> 136 } 82 } 137 } 83 138 84 void G4HumanPhantomPhysicsList::SetCuts() 139 void G4HumanPhantomPhysicsList::SetCuts() 85 { 140 { 86 // The production threshold is fixed to 0.1 141 // The production threshold is fixed to 0.1 mm for all the particles 87 // Secondary particles with a range bigger t 142 // Secondary particles with a range bigger than 0.1 mm 88 // are generated; otherwise their energy is 143 // are generated; otherwise their energy is considered deposited locally 89 144 90 defaultCutValue = 1. * mm; << 145 defaultCutValue = 0.1 * mm; 91 146 92 const G4double cutForGamma = defaultCutValue 147 const G4double cutForGamma = defaultCutValue; 93 const G4double cutForElectron = defaultCutVa 148 const G4double cutForElectron = defaultCutValue; 94 const G4double cutForPositron = defaultCutVa 149 const G4double cutForPositron = defaultCutValue; 95 150 96 SetCutValue(cutForGamma, "gamma"); 151 SetCutValue(cutForGamma, "gamma"); 97 SetCutValue(cutForElectron, "e-"); 152 SetCutValue(cutForElectron, "e-"); 98 SetCutValue(cutForPositron, "e+"); 153 SetCutValue(cutForPositron, "e+"); 99 154 100 // Set the secondary production cut lower th 155 // Set the secondary production cut lower than 990. eV 101 // Very important for high precision of lowe 156 // Very important for high precision of lowenergy processes at low energies 102 157 103 G4double lowLimit = 250. * eV; 158 G4double lowLimit = 250. * eV; 104 G4double highLimit = 100. * GeV; 159 G4double highLimit = 100. * GeV; 105 G4ProductionCutsTable::GetProductionCutsTabl 160 G4ProductionCutsTable::GetProductionCutsTable()->SetEnergyRange(lowLimit, highLimit); 106 161 107 if (verboseLevel>0) DumpCutValuesTable(); 162 if (verboseLevel>0) DumpCutValuesTable(); 108 } 163 } 109 164