File indexing completed on 2025-02-23 09:22:01
0001
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021
0022
0023
0024
0025
0026
0027
0028
0029
0030 #include "PhysicsList.hh"
0031
0032 #include "ChemistryList.hh"
0033 #include "ParallelWorldPhysics.hh"
0034
0035 #include "G4DecayPhysics.hh"
0036 #include "G4EmDNAPhysics.hh"
0037 #include "G4EmDNAPhysics_option2.hh"
0038 #include "G4EmDNAPhysics_option4.hh"
0039 #include "G4EmDNAPhysics_option6.hh"
0040 #include "G4EmParameters.hh"
0041 #include "G4RadioactiveDecayPhysics.hh"
0042 #include "G4SystemOfUnits.hh"
0043
0044
0045
0046 PhysicsList::PhysicsList(G4int phylist)
0047 {
0048 SetDefaultCutValue(1.0 * micrometer);
0049 SetVerboseLevel(1);
0050
0051 UseCoupledTransportation();
0052
0053 if (phylist == 0) {
0054 RegisterPhysics(new G4EmDNAPhysics());
0055 }
0056 else if (phylist == 2) {
0057 RegisterPhysics(new G4EmDNAPhysics_option2());
0058 }
0059 else if (phylist == 4) {
0060 RegisterPhysics(new G4EmDNAPhysics_option4());
0061 }
0062 else if (phylist == 6) {
0063 RegisterPhysics(new G4EmDNAPhysics_option6());
0064 }
0065 else {
0066 G4ExceptionDescription errmsg;
0067 errmsg << "Recommend to option 2, 4, 6 or default" << G4endl;
0068 G4Exception("PhysicsList::PhysicsList", "", FatalException, errmsg);
0069 }
0070 RegisterPhysics(new G4DecayPhysics());
0071 RegisterPhysics(new G4RadioactiveDecayPhysics());
0072
0073 G4bool useParallelPhysicsWorld = false;
0074 if (useParallelPhysicsWorld) {
0075 RegisterPhysics(new ParallelWorldPhysics("DNAWorld", true));
0076 }
0077
0078 RegisterPhysics(new ChemistryList());
0079
0080 G4ProductionCutsTable::GetProductionCutsTable()->SetEnergyRange(100 * eV, 1 * GeV);
0081 }
0082
0083