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 /// \file electromagnetic/TestEm9/src/PhysicsList.cc 27 /// \brief Implementation of the PhysicsList class 28 // 29 // 30 //--------------------------------------------------------------------------- 31 // 32 // ClassName: PhysicsList 33 // 34 // Author: V.Ivanchenko 14.10.2002 35 // 36 // Modified: 37 // 17.11.06 Use components from physics_lists subdirectory (V.Ivanchenko) 38 // 24.10.12 Migrated to new stopping and ion physics (A.Ribon) 39 // 40 //---------------------------------------------------------------------------- 41 // 42 43 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 44 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 45 46 #include "PhysicsList.hh" 47 48 #include "PhysListEmStandard.hh" 49 #include "PhysicsListMessenger.hh" 50 #include "StepMax.hh" 51 52 #include "G4DecayPhysics.hh" 53 #include "G4Electron.hh" 54 #include "G4EmExtraPhysics.hh" 55 #include "G4EmLivermorePhysics.hh" 56 #include "G4EmLowEPPhysics.hh" 57 #include "G4EmPenelopePhysics.hh" 58 #include "G4EmStandardPhysics.hh" 59 #include "G4EmStandardPhysicsGS.hh" 60 #include "G4EmStandardPhysicsSS.hh" 61 #include "G4EmStandardPhysicsWVI.hh" 62 #include "G4EmStandardPhysics_option1.hh" 63 #include "G4EmStandardPhysics_option2.hh" 64 #include "G4EmStandardPhysics_option3.hh" 65 #include "G4EmStandardPhysics_option4.hh" 66 #include "G4Gamma.hh" 67 #include "G4HadronElasticPhysics.hh" 68 #include "G4HadronInelasticQBBC.hh" 69 #include "G4IonPhysics.hh" 70 #include "G4LossTableManager.hh" 71 #include "G4ParticleTable.hh" 72 #include "G4ParticleTypes.hh" 73 #include "G4Positron.hh" 74 #include "G4ProcessManager.hh" 75 #include "G4Proton.hh" 76 #include "G4RegionStore.hh" 77 #include "G4StoppingPhysics.hh" 78 #include "G4SystemOfUnits.hh" 79 #include "G4UnitsTable.hh" 80 81 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 82 83 PhysicsList::PhysicsList() 84 : G4VModularPhysicsList(), 85 fEmPhysicsList(nullptr), 86 fDecayPhysicsList(nullptr), 87 fStepMaxProcess(nullptr), 88 fMessenger(nullptr) 89 { 90 G4LossTableManager::Instance(); 91 SetDefaultCutValue(1 * mm); 92 93 fMessenger = new PhysicsListMessenger(this); 94 fStepMaxProcess = new StepMax(); 95 96 // Initilise flags 97 98 SetVerboseLevel(1); 99 100 fHelIsRegisted = false; 101 fBicIsRegisted = false; 102 fGnucIsRegisted = false; 103 fStopIsRegisted = false; 104 105 // EM physics 106 fEmName = G4String("emstandard"); 107 fEmPhysicsList = new G4EmStandardPhysics(); 108 109 // Decay Physics is always defined 110 fDecayPhysicsList = new G4DecayPhysics(); 111 } 112 113 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 114 115 PhysicsList::~PhysicsList() 116 { 117 delete fMessenger; 118 delete fDecayPhysicsList; 119 delete fEmPhysicsList; 120 delete fStepMaxProcess; 121 for (size_t i = 0; i < fHadronPhys.size(); i++) { 122 delete fHadronPhys[i]; 123 } 124 } 125 126 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 127 128 void PhysicsList::ConstructParticle() 129 { 130 fDecayPhysicsList->ConstructParticle(); 131 } 132 133 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 134 135 void PhysicsList::ConstructProcess() 136 { 137 AddTransportation(); 138 fEmPhysicsList->ConstructProcess(); 139 fDecayPhysicsList->ConstructProcess(); 140 for (size_t i = 0; i < fHadronPhys.size(); ++i) { 141 fHadronPhys[i]->ConstructProcess(); 142 } 143 AddStepMax(); 144 } 145 146 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 147 148 void PhysicsList::AddPhysicsList(const G4String& name) 149 { 150 if (verboseLevel > 1) G4cout << "PhysicsList::AddPhysicsList: <" << name << ">" << G4endl; 151 152 if (name == fEmName) return; 153 154 if (name == "emstandard") { 155 fEmName = name; 156 delete fEmPhysicsList; 157 fEmPhysicsList = new G4EmStandardPhysics(verboseLevel); 158 if (verboseLevel > 0) G4cout << "PhysicsList::Set " << name << " EM physics" << G4endl; 159 } 160 else if (name == "emstandard_opt1") { 161 fEmName = name; 162 delete fEmPhysicsList; 163 fEmPhysicsList = new G4EmStandardPhysics_option1(verboseLevel); 164 if (verboseLevel > 0) G4cout << "PhysicsList::Set " << name << " EM physics" << G4endl; 165 } 166 else if (name == "emstandard_opt2") { 167 fEmName = name; 168 delete fEmPhysicsList; 169 fEmPhysicsList = new G4EmStandardPhysics_option2(verboseLevel); 170 if (verboseLevel > 0) G4cout << "PhysicsList::Set " << name << " EM physics" << G4endl; 171 } 172 else if (name == "emstandard_opt3") { 173 fEmName = name; 174 delete fEmPhysicsList; 175 fEmPhysicsList = new G4EmStandardPhysics_option3(verboseLevel); 176 if (verboseLevel > 0) G4cout << "PhysicsList::Set " << name << " EM physics" << G4endl; 177 } 178 else if (name == "emstandard_opt4") { 179 fEmName = name; 180 delete fEmPhysicsList; 181 fEmPhysicsList = new G4EmStandardPhysics_option4(verboseLevel); 182 if (verboseLevel > 0) G4cout << "PhysicsList::Set " << name << " EM physics" << G4endl; 183 } 184 else if (name == "local") { 185 fEmName = name; 186 delete fEmPhysicsList; 187 fEmPhysicsList = new PhysListEmStandard(verboseLevel); 188 if (verboseLevel > 0) G4cout << "PhysicsList::Set " << name << " EM physics" << G4endl; 189 } 190 else if (name == "emlivermore") { 191 fEmName = name; 192 delete fEmPhysicsList; 193 fEmPhysicsList = new G4EmLivermorePhysics(verboseLevel); 194 } 195 else if (name == "empenelope") { 196 fEmName = name; 197 delete fEmPhysicsList; 198 fEmPhysicsList = new G4EmPenelopePhysics(verboseLevel); 199 } 200 else if (name == "emlowenergy") { 201 fEmName = name; 202 delete fEmPhysicsList; 203 fEmPhysicsList = new G4EmLowEPPhysics(verboseLevel); 204 } 205 else if (name == "emstandardGS") { 206 fEmName = name; 207 delete fEmPhysicsList; 208 fEmPhysicsList = new G4EmStandardPhysicsGS(verboseLevel); 209 } 210 else if (name == "emstandardSS") { 211 fEmName = name; 212 delete fEmPhysicsList; 213 fEmPhysicsList = new G4EmStandardPhysicsSS(verboseLevel); 214 } 215 else if (name == "emstandardWVI") { 216 fEmName = name; 217 delete fEmPhysicsList; 218 fEmPhysicsList = new G4EmStandardPhysicsWVI(); 219 } 220 else if (name == "elastic" && !fHelIsRegisted) { 221 fHadronPhys.push_back(new G4HadronElasticPhysics()); 222 fHelIsRegisted = true; 223 if (verboseLevel > 0) G4cout << "PhysicsList::Add hadron elastic physics" << G4endl; 224 } 225 else if (name == "binary" && !fBicIsRegisted) { 226 fHadronPhys.push_back(new G4HadronInelasticQBBC()); 227 fHadronPhys.push_back(new G4IonPhysics()); 228 fBicIsRegisted = true; 229 if (verboseLevel > 0) 230 G4cout << "PhysicsList::Add hadron inelastic physics from <QBBC>" << G4endl; 231 } 232 else if (name == "gamma_nuc" && !fGnucIsRegisted) { 233 fHadronPhys.push_back(new G4EmExtraPhysics()); 234 fGnucIsRegisted = true; 235 if (verboseLevel > 0) G4cout << "PhysicsList::Add gamma- and electro-nuclear physics" << G4endl; 236 } 237 else if (name == "stopping" && !fStopIsRegisted) { 238 fHadronPhys.push_back(new G4StoppingPhysics()); 239 fStopIsRegisted = true; 240 if (verboseLevel > 0) G4cout << "PhysicsList::Add stopping physics" << G4endl; 241 } 242 else { 243 G4cout << "PhysicsList::AddPhysicsList: <" << name << ">" 244 << " is not defined" << G4endl; 245 } 246 } 247 248 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 249 250 void PhysicsList::AddStepMax() 251 { 252 // Step limitation seen as a process 253 254 auto particleIterator = GetParticleIterator(); 255 particleIterator->reset(); 256 while ((*particleIterator)()) { 257 G4ParticleDefinition* particle = particleIterator->value(); 258 G4ProcessManager* pmanager = particle->GetProcessManager(); 259 260 if (fStepMaxProcess->IsApplicable(*particle) && !particle->IsShortLived()) { 261 pmanager->AddDiscreteProcess(fStepMaxProcess); 262 } 263 } 264 } 265 266 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 267