Geant4 Cross Reference |
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 // Previous authors: G. Guerrieri and M. G. Pia, INFN Genova, Italy 28 // Authors (since 2007): S. Guatelli, University of Wollongong, Australia 29 // 30 // 31 // 32 // ********************************** 33 // * * 34 // * G4HumanPhantomPhysicsList.cc * 35 // * * 36 // ********************************** 37 // 38 // 39 #include "G4HumanPhantomPhysicsList.hh" 40 41 #include "G4SystemOfUnits.hh" 42 #include "G4ParticleDefinition.hh" 43 #include "G4ProductionCutsTable.hh" 44 #include "G4ProcessManager.hh" 45 #include "G4ParticleTypes.hh" 46 #include "G4UnitsTable.hh" 47 #include "G4ios.hh" 48 #include "G4VModularPhysicsList.hh" 49 #include "G4EmStandardPhysics_option4.hh" 50 #include "G4DecayPhysics.hh" 51 52 G4HumanPhantomPhysicsList::G4HumanPhantomPhysicsList(): G4VModularPhysicsList() 53 { 54 SetVerboseLevel(1); 55 emPhysicsList = new G4EmStandardPhysics_option4(1); 56 decPhysicsList = new G4DecayPhysics(); 57 // Alternatively you can substitute this physics list 58 // with the LowEnergy Livermore or LowEnergy Penelope: 59 // emPhysicsList = new G4EmLivermorePhysics(); 60 // Low Energy based on Livermore Evaluated Data Libraries 61 // 62 // Penelope physics 63 //emPhysicsList = new G4EmPenelopePhysics(); 64 } 65 66 G4HumanPhantomPhysicsList::~G4HumanPhantomPhysicsList() 67 { 68 delete decPhysicsList; 69 delete emPhysicsList; 70 } 71 72 void G4HumanPhantomPhysicsList::ConstructParticle() 73 { 74 decPhysicsList -> ConstructParticle(); 75 } 76 77 void G4HumanPhantomPhysicsList::ConstructProcess() 78 { 79 AddTransportation(); 80 emPhysicsList -> ConstructProcess(); 81 decPhysicsList -> ConstructProcess(); 82 } 83 84 void G4HumanPhantomPhysicsList::SetCuts() 85 { 86 // The production threshold is fixed to 0.1 mm for all the particles 87 // Secondary particles with a range bigger than 0.1 mm 88 // are generated; otherwise their energy is considered deposited locally 89 90 defaultCutValue = 1. * mm; 91 92 const G4double cutForGamma = defaultCutValue; 93 const G4double cutForElectron = defaultCutValue; 94 const G4double cutForPositron = defaultCutValue; 95 96 SetCutValue(cutForGamma, "gamma"); 97 SetCutValue(cutForElectron, "e-"); 98 SetCutValue(cutForPositron, "e+"); 99 100 // Set the secondary production cut lower than 990. eV 101 // Very important for high precision of lowenergy processes at low energies 102 103 G4double lowLimit = 250. * eV; 104 G4double highLimit = 100. * GeV; 105 G4ProductionCutsTable::GetProductionCutsTable()->SetEnergyRange(lowLimit, highLimit); 106 107 if (verboseLevel>0) DumpCutValuesTable(); 108 } 109