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 /// \file medical/electronScattering2/src/PhysicsList.cc 28 /// \brief Implementation of the PhysicsList class 29 30 #include "PhysicsList.hh" 31 32 #include "PhysicsListMessenger.hh" 33 #include "StepMax.hh" 34 35 #include "G4Decay.hh" 36 #include "G4EmBuilder.hh" 37 #include "G4EmLivermorePhysics.hh" 38 #include "G4EmLowEPPhysics.hh" 39 #include "G4EmPenelopePhysics.hh" 40 #include "G4EmStandardPhysics.hh" 41 #include "G4EmStandardPhysicsGS.hh" 42 #include "G4EmStandardPhysicsSS.hh" 43 #include "G4EmStandardPhysicsWVI.hh" 44 #include "G4EmStandardPhysics_option1.hh" 45 #include "G4EmStandardPhysics_option2.hh" 46 #include "G4EmStandardPhysics_option3.hh" 47 #include "G4EmStandardPhysics_option4.hh" 48 #include "G4LossTableManager.hh" 49 #include "G4ParticleDefinition.hh" 50 #include "G4ProcessManager.hh" 51 #include "G4SystemOfUnits.hh" 52 #include "G4UnitsTable.hh" 53 54 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 55 56 PhysicsList::PhysicsList() : G4VModularPhysicsList() 57 { 58 fMessenger = new PhysicsListMessenger(this); 59 60 // EM physics 61 fEmName = G4String("emstandard_opt4"); 62 fEmPhysicsList = new G4EmStandardPhysics_option4(1); 63 if (verboseLevel > -1) { 64 G4cout << "PhysicsList::Constructor with default list: <" << fEmName << ">" << G4endl; 65 } 66 67 G4LossTableManager::Instance(); 68 SetVerboseLevel(1); 69 } 70 71 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 72 73 PhysicsList::~PhysicsList() 74 { 75 delete fEmPhysicsList; 76 delete fMessenger; 77 } 78 79 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 80 81 void PhysicsList::ConstructParticle() 82 { 83 // minimal set of particles for EM physics 84 G4EmBuilder::ConstructMinimalEmSet(); 85 } 86 87 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 88 89 void PhysicsList::ConstructProcess() 90 { 91 AddTransportation(); 92 fEmPhysicsList->ConstructProcess(); 93 AddDecay(); 94 AddStepMax(); 95 } 96 97 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 98 99 void PhysicsList::AddDecay() 100 { 101 // Add Decay Process 102 103 G4Decay* fDecayProcess = new G4Decay(); 104 105 auto particleIterator = GetParticleIterator(); 106 particleIterator->reset(); 107 while ((*particleIterator)()) { 108 G4ParticleDefinition* particle = particleIterator->value(); 109 G4ProcessManager* pmanager = particle->GetProcessManager(); 110 111 if (fDecayProcess->IsApplicable(*particle) && !particle->IsShortLived()) { 112 pmanager->AddProcess(fDecayProcess); 113 114 // set ordering for PostStepDoIt and AtRestDoIt 115 pmanager->SetProcessOrdering(fDecayProcess, idxPostStep); 116 pmanager->SetProcessOrdering(fDecayProcess, idxAtRest); 117 } 118 } 119 } 120 121 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 122 123 void PhysicsList::AddStepMax() 124 { 125 // Step limitation seen as a process 126 StepMax* stepMaxProcess = new StepMax(); 127 128 auto particleIterator = GetParticleIterator(); 129 particleIterator->reset(); 130 while ((*particleIterator)()) { 131 G4ParticleDefinition* particle = particleIterator->value(); 132 G4ProcessManager* pmanager = particle->GetProcessManager(); 133 134 if (stepMaxProcess->IsApplicable(*particle)) { 135 pmanager->AddDiscreteProcess(stepMaxProcess); 136 } 137 } 138 } 139 140 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 141 142 void PhysicsList::AddPhysicsList(const G4String& name) 143 { 144 if (verboseLevel > -1) { 145 G4cout << "PhysicsList::AddPhysicsList: <" << name << ">" << G4endl; 146 } 147 148 if (name == fEmName) return; 149 150 if (name == "emstandard_opt0") { 151 fEmName = name; 152 delete fEmPhysicsList; 153 fEmPhysicsList = new G4EmStandardPhysics(1); 154 } 155 else if (name == "emstandard_opt1") { 156 fEmName = name; 157 delete fEmPhysicsList; 158 fEmPhysicsList = new G4EmStandardPhysics_option1(1); 159 } 160 else if (name == "emstandard_opt2") { 161 fEmName = name; 162 delete fEmPhysicsList; 163 fEmPhysicsList = new G4EmStandardPhysics_option2(1); 164 } 165 else if (name == "emstandard_opt3") { 166 fEmName = name; 167 delete fEmPhysicsList; 168 fEmPhysicsList = new G4EmStandardPhysics_option3(1); 169 } 170 else if (name == "emstandard_opt4") { 171 fEmName = name; 172 delete fEmPhysicsList; 173 fEmPhysicsList = new G4EmStandardPhysics_option4(1); 174 } 175 else if (name == "emlowenergy") { 176 fEmName = name; 177 delete fEmPhysicsList; 178 fEmPhysicsList = new G4EmLowEPPhysics(1); 179 } 180 else if (name == "emstandardSS") { 181 fEmName = name; 182 delete fEmPhysicsList; 183 fEmPhysicsList = new G4EmStandardPhysicsSS(1); 184 } 185 else if (name == "emstandardWVI") { 186 fEmName = name; 187 delete fEmPhysicsList; 188 fEmPhysicsList = new G4EmStandardPhysicsWVI(1); 189 } 190 else if (name == "emstandardGS") { 191 fEmName = name; 192 delete fEmPhysicsList; 193 fEmPhysicsList = new G4EmStandardPhysicsGS(1); 194 } 195 else if (name == "emlivermore") { 196 fEmName = name; 197 delete fEmPhysicsList; 198 fEmPhysicsList = new G4EmLivermorePhysics(1); 199 } 200 else if (name == "empenelope") { 201 fEmName = name; 202 delete fEmPhysicsList; 203 fEmPhysicsList = new G4EmPenelopePhysics(1); 204 } 205 else { 206 G4ExceptionDescription description; 207 description << " " 208 << "PhysicsList::AddPhysicsList: <" << name << "> is not defined"; 209 G4Exception("PhysicsList::AddPhysicsList", "electronScattering2_F001", FatalException, 210 description); 211 } 212 } 213 214 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 215