| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 //
27 //....oooOO0OOooo........oooOO0OOooo........oo
28 //....oooOO0OOooo........oooOO0OOooo........oo
29
30 #include "PhysicsList.hh"
31 #include "PhysicsListMessenger.hh"
32
33 #include "PhysListEmStandard.hh"
34
35 #include "G4EmStandardPhysics.hh"
36 #include "G4EmStandardPhysics_option1.hh"
37 #include "G4EmStandardPhysics_option2.hh"
38 #include "G4EmStandardPhysics_option3.hh"
39
40 #include "G4LossTableManager.hh"
41 #include "G4SystemOfUnits.hh"
42 #include "G4UnitsTable.hh"
43
44 //....oooOO0OOooo........oooOO0OOooo........oo
45
46 PhysicsList::PhysicsList() : G4VModularPhysics
47 {
48 pMessenger = new PhysicsListMessenger(this);
49 SetVerboseLevel(1);
50
51 // EM physics
52 emName = G4String("local");
53 emPhysicsList = new PhysListEmStandard(emNam
54
55 SetDefaultCutValue(1*mm);
56 G4LossTableManager::Instance();
57 }
58
59 //....oooOO0OOooo........oooOO0OOooo........oo
60
61 PhysicsList::~PhysicsList()
62 {
63 delete pMessenger;
64 }
65
66 //....oooOO0OOooo........oooOO0OOooo........oo
67
68 #include "G4BosonConstructor.hh"
69 #include "G4LeptonConstructor.hh"
70 #include "G4MesonConstructor.hh"
71 #include "G4BosonConstructor.hh"
72 #include "G4BaryonConstructor.hh"
73 #include "G4IonConstructor.hh"
74 #include "G4ShortLivedConstructor.hh"
75
76 void PhysicsList::ConstructParticle()
77 {
78 G4BosonConstructor pBosonConstructor;
79 pBosonConstructor.ConstructParticle();
80
81 G4LeptonConstructor pLeptonConstructor;
82 pLeptonConstructor.ConstructParticle();
83
84 G4MesonConstructor pMesonConstructor;
85 pMesonConstructor.ConstructParticle();
86
87 G4BaryonConstructor pBaryonConstructor;
88 pBaryonConstructor.ConstructParticle();
89
90 G4IonConstructor pIonConstructor;
91 pIonConstructor.ConstructParticle();
92
93 G4ShortLivedConstructor pShortLivedConstru
94 pShortLivedConstructor.ConstructParticle()
95 }
96
97 //....oooOO0OOooo........oooOO0OOooo........oo
98
99 #include "G4PhysicsListHelper.hh"
100 #include "G4Decay.hh"
101
102 void PhysicsList::ConstructProcess()
103 {
104 AddTransportation();
105
106 // electromagnetic Physics List
107 //
108 emPhysicsList->ConstructProcess();
109
110 // Add Decay Process
111 //
112 G4Decay* fDecayProcess = new G4Decay();
113 G4PhysicsListHelper* ph = G4PhysicsListHelpe
114
115 auto particleIterator=GetParticleIterator();
116 particleIterator->reset();
117 while( (*particleIterator)() ){
118 G4ParticleDefinition* particle = particleI
119 if (fDecayProcess->IsApplicable(*particle)
120 ph ->RegisterProcess(fDecayProcess, part
121 }
122 }
123
124 //....oooOO0OOooo........oooOO0OOooo........oo
125
126 void PhysicsList::AddPhysicsList(const G4Strin
127 {
128 if (verboseLevel>1) {
129 G4cout << "PhysicsList::AddPhysicsList: <"
130 }
131
132 if (name == emName) return;
133
134 if (name == "local") {
135
136 emName = name;
137 delete emPhysicsList;
138 emPhysicsList = new PhysListEmStandard(nam
139
140 } else if (name == "emstandard_opt0") {
141
142 emName = name;
143 delete emPhysicsList;
144 emPhysicsList = new G4EmStandardPhysics();
145
146 } else if (name == "emstandard_opt1") {
147
148 emName = name;
149 delete emPhysicsList;
150 emPhysicsList = new G4EmStandardPhysics_op
151
152 } else if (name == "emstandard_opt2") {
153
154 emName = name;
155 delete emPhysicsList;
156 emPhysicsList = new G4EmStandardPhysics_op
157
158 } else if (name == "emstandard_opt3") {
159
160 emName = name;
161 delete emPhysicsList;
162 emPhysicsList = new G4EmStandardPhysics_op
163
164 } else {
165
166 G4cout << "PhysicsList::AddPhysicsList: <"
167 << " is not defined"
168 << G4endl;
169 }
170 }
171
172 //....oooOO0OOooo........oooOO0OOooo........oo
173