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 /// \file electromagnetic/TestEm18/src/Physics 26 /// \file electromagnetic/TestEm18/src/PhysicsList.cc 27 /// \brief Implementation of the PhysicsList c 27 /// \brief Implementation of the PhysicsList class 28 // 28 // >> 29 // $Id: PhysicsList.cc 100275 2016-10-17 08:29:19Z gcosmo $ 29 // 30 // 30 //....oooOO0OOooo........oooOO0OOooo........oo 31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 31 //....oooOO0OOooo........oooOO0OOooo........oo 32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 32 33 33 #include "PhysicsList.hh" 34 #include "PhysicsList.hh" >> 35 #include "PhysicsListMessenger.hh" 34 36 >> 37 #include "PhysListEmStandard.hh" 35 #include "PhysListEmLivermore.hh" 38 #include "PhysListEmLivermore.hh" 36 #include "PhysListEmPenelope.hh" 39 #include "PhysListEmPenelope.hh" 37 #include "PhysListEmStandard.hh" << 38 #include "PhysicsListMessenger.hh" << 39 40 40 #include "G4LossTableManager.hh" 41 #include "G4LossTableManager.hh" 41 #include "G4SystemOfUnits.hh" << 42 #include "G4UnitsTable.hh" 42 #include "G4UnitsTable.hh" >> 43 #include "G4SystemOfUnits.hh" 43 44 44 //....oooOO0OOooo........oooOO0OOooo........oo 45 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 45 46 46 PhysicsList::PhysicsList() << 47 PhysicsList::PhysicsList() : G4VModularPhysicsList(), >> 48 fMessenger(0), >> 49 fEmPhysicsList(0) 47 { 50 { 48 G4LossTableManager::Instance(); 51 G4LossTableManager::Instance(); 49 fMessenger = new PhysicsListMessenger(this); << 52 fMessenger = new PhysicsListMessenger(this); 50 << 53 51 // EM physics 54 // EM physics 52 fEmName = G4String("standard"); 55 fEmName = G4String("standard"); 53 fEmPhysicsList = new PhysListEmStandard(fEmN 56 fEmPhysicsList = new PhysListEmStandard(fEmName); 54 << 57 55 SetDefaultCutValue(1. * mm); << 58 SetDefaultCutValue(1.*mm); 56 << 59 57 SetVerboseLevel(1); 60 SetVerboseLevel(1); 58 } 61 } 59 62 60 //....oooOO0OOooo........oooOO0OOooo........oo 63 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 61 64 62 PhysicsList::~PhysicsList() 65 PhysicsList::~PhysicsList() 63 { 66 { 64 delete fEmPhysicsList; 67 delete fEmPhysicsList; 65 delete fMessenger; << 68 delete fMessenger; 66 } 69 } 67 70 68 //....oooOO0OOooo........oooOO0OOooo........oo 71 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 69 72 70 void PhysicsList::AddPhysicsList(const G4Strin 73 void PhysicsList::AddPhysicsList(const G4String& name) 71 { 74 { 72 if (verboseLevel > 1) { << 75 if (verboseLevel>1) { 73 G4cout << "PhysicsList::AddPhysicsList: <" 76 G4cout << "PhysicsList::AddPhysicsList: <" << name << ">" << G4endl; 74 } 77 } 75 78 76 if (name == fEmName) return; 79 if (name == fEmName) return; 77 80 78 if (name == "standard") { 81 if (name == "standard") { >> 82 79 fEmName = name; 83 fEmName = name; 80 delete fEmPhysicsList; 84 delete fEmPhysicsList; 81 fEmPhysicsList = new PhysListEmStandard(na 85 fEmPhysicsList = new PhysListEmStandard(name); 82 } << 86 83 else if (name == "livermore") { << 87 } else if (name == "livermore") { >> 88 84 fEmName = name; 89 fEmName = name; 85 delete fEmPhysicsList; 90 delete fEmPhysicsList; 86 fEmPhysicsList = new PhysListEmLivermore(n 91 fEmPhysicsList = new PhysListEmLivermore(name); 87 } << 92 88 else if (name == "penelope") { << 93 } else if (name == "penelope") { >> 94 89 fEmName = name; 95 fEmName = name; 90 delete fEmPhysicsList; 96 delete fEmPhysicsList; 91 fEmPhysicsList = new PhysListEmPenelope(na 97 fEmPhysicsList = new PhysListEmPenelope(name); 92 } << 98 93 else { << 99 } else { >> 100 94 G4cout << "PhysicsList::AddPhysicsList: <" 101 G4cout << "PhysicsList::AddPhysicsList: <" << name << ">" 95 << " is not defined" << G4endl; << 102 << " is not defined" >> 103 << G4endl; 96 } 104 } 97 } 105 } 98 106 99 //....oooOO0OOooo........oooOO0OOooo........oo 107 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 100 108 101 // Bosons 109 // Bosons 102 #include "G4ChargedGeantino.hh" 110 #include "G4ChargedGeantino.hh" 103 #include "G4Gamma.hh" << 104 #include "G4Geantino.hh" 111 #include "G4Geantino.hh" >> 112 #include "G4Gamma.hh" 105 113 106 // leptons 114 // leptons 107 #include "G4Electron.hh" 115 #include "G4Electron.hh" 108 #include "G4MuonMinus.hh" << 109 #include "G4MuonPlus.hh" << 110 #include "G4Positron.hh" 116 #include "G4Positron.hh" 111 117 >> 118 #include "G4MuonPlus.hh" >> 119 #include "G4MuonMinus.hh" >> 120 112 // Mesons 121 // Mesons 113 #include "G4KaonMinus.hh" << 114 #include "G4KaonPlus.hh" << 115 #include "G4PionMinus.hh" << 116 #include "G4PionPlus.hh" 122 #include "G4PionPlus.hh" >> 123 #include "G4PionMinus.hh" >> 124 >> 125 #include "G4KaonPlus.hh" >> 126 #include "G4KaonMinus.hh" 117 127 118 // Baryons 128 // Baryons 119 #include "G4AntiNeutron.hh" << 129 #include "G4Proton.hh" 120 #include "G4AntiProton.hh" 130 #include "G4AntiProton.hh" 121 #include "G4Neutron.hh" 131 #include "G4Neutron.hh" 122 #include "G4Proton.hh" << 132 #include "G4AntiNeutron.hh" 123 133 124 // Nuclei 134 // Nuclei 125 #include "G4Alpha.hh" << 126 #include "G4Deuteron.hh" 135 #include "G4Deuteron.hh" 127 #include "G4GenericIon.hh" << 128 #include "G4Triton.hh" 136 #include "G4Triton.hh" >> 137 #include "G4Alpha.hh" >> 138 #include "G4GenericIon.hh" 129 139 130 void PhysicsList::ConstructParticle() 140 void PhysicsList::ConstructParticle() 131 { 141 { 132 // pseudo-particles << 142 // pseudo-particles 133 G4Geantino::GeantinoDefinition(); 143 G4Geantino::GeantinoDefinition(); 134 G4ChargedGeantino::ChargedGeantinoDefinition 144 G4ChargedGeantino::ChargedGeantinoDefinition(); 135 << 145 136 // gamma << 146 // gamma 137 G4Gamma::GammaDefinition(); 147 G4Gamma::GammaDefinition(); 138 148 139 // leptons << 149 // leptons 140 G4Electron::ElectronDefinition(); 150 G4Electron::ElectronDefinition(); 141 G4Positron::PositronDefinition(); 151 G4Positron::PositronDefinition(); 142 G4MuonPlus::MuonPlusDefinition(); 152 G4MuonPlus::MuonPlusDefinition(); 143 G4MuonMinus::MuonMinusDefinition(); 153 G4MuonMinus::MuonMinusDefinition(); 144 << 154 145 // mesons << 155 // mesons 146 G4PionPlus::PionPlusDefinition(); 156 G4PionPlus::PionPlusDefinition(); 147 G4PionMinus::PionMinusDefinition(); 157 G4PionMinus::PionMinusDefinition(); 148 G4KaonPlus::KaonPlusDefinition(); 158 G4KaonPlus::KaonPlusDefinition(); 149 G4KaonMinus::KaonMinusDefinition(); 159 G4KaonMinus::KaonMinusDefinition(); 150 << 160 151 // baryons << 161 // baryons 152 G4Proton::ProtonDefinition(); 162 G4Proton::ProtonDefinition(); 153 G4AntiProton::AntiProtonDefinition(); 163 G4AntiProton::AntiProtonDefinition(); 154 G4Neutron::NeutronDefinition(); 164 G4Neutron::NeutronDefinition(); 155 G4AntiNeutron::AntiNeutronDefinition(); 165 G4AntiNeutron::AntiNeutronDefinition(); 156 << 166 157 // ions << 167 // ions 158 G4Deuteron::DeuteronDefinition(); 168 G4Deuteron::DeuteronDefinition(); 159 G4Triton::TritonDefinition(); 169 G4Triton::TritonDefinition(); 160 G4Alpha::AlphaDefinition(); 170 G4Alpha::AlphaDefinition(); 161 G4GenericIon::GenericIonDefinition(); 171 G4GenericIon::GenericIonDefinition(); 162 } 172 } 163 173 164 //....oooOO0OOooo........oooOO0OOooo........oo 174 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 165 175 166 void PhysicsList::ConstructProcess() 176 void PhysicsList::ConstructProcess() 167 { 177 { 168 AddTransportation(); 178 AddTransportation(); 169 fEmPhysicsList->ConstructProcess(); 179 fEmPhysicsList->ConstructProcess(); 170 AddStepMax(); 180 AddStepMax(); 171 } 181 } 172 182 173 //....oooOO0OOooo........oooOO0OOooo........oo 183 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 174 184 175 #include "StepMax.hh" << 176 << 177 #include "G4ProcessManager.hh" 185 #include "G4ProcessManager.hh" >> 186 #include "StepMax.hh" 178 187 179 void PhysicsList::AddStepMax() 188 void PhysicsList::AddStepMax() 180 { 189 { 181 // Step limitation seen as a process 190 // Step limitation seen as a process 182 StepMax* stepMaxProcess = new StepMax(); 191 StepMax* stepMaxProcess = new StepMax(); 183 192 184 auto particleIterator = GetParticleIterator( << 193 auto particleIterator=GetParticleIterator(); 185 particleIterator->reset(); 194 particleIterator->reset(); 186 while ((*particleIterator)()) { << 195 while ((*particleIterator)()){ 187 G4ParticleDefinition* particle = particleI << 196 G4ParticleDefinition* particle = particleIterator->value(); 188 G4ProcessManager* pmanager = particle->Get << 197 G4ProcessManager* pmanager = particle->GetProcessManager(); 189 << 198 190 if (stepMaxProcess->IsApplicable(*particle << 199 if (stepMaxProcess->IsApplicable(*particle) && !particle->IsShortLived()) 191 pmanager->AddDiscreteProcess(stepMaxProc << 200 { 192 } << 201 pmanager ->AddDiscreteProcess(stepMaxProcess); >> 202 } 193 } 203 } 194 } 204 } 195 205 196 //....oooOO0OOooo........oooOO0OOooo........oo 206 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 197 207