| 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 //--------------------------------------------
28 //
29 // ClassName: G4PhysListFactory
30 //
31 // Author: 21 April 2008 V. Ivanchenko
32 //
33 // Modified:
34 //
35 // 2023.04.12 A.Ribon added _HPT variants (i.e
36 // 2014.08.05 K.L.Genser used provision for Ha
37 // Shielding for ShieldingM
38 //
39 //--------------------------------------------
40 //
41
42 #include "G4PhysListFactory.hh"
43 #include "FTFP_BERT.hh"
44 #include "FTFP_BERT_HP.hh"
45 #include "FTFP_BERT_TRV.hh"
46 #include "FTFP_BERT_ATL.hh"
47 #include "FTFQGSP_BERT.hh"
48 #include "FTFP_INCLXX.hh"
49 #include "FTFP_INCLXX_HP.hh"
50 #include "FTF_BIC.hh"
51 #include "LBE.hh"
52 #include "QBBC.hh"
53 #include "QGSP_BERT.hh"
54 #include "QGSP_BERT_HP.hh"
55 #include "QGSP_BIC.hh"
56 #include "QGSP_BIC_HP.hh"
57 #include "QGSP_BIC_HPT.hh"
58 #include "QGSP_BIC_AllHP.hh"
59 #include "QGSP_FTFP_BERT.hh"
60 #include "QGS_BIC.hh"
61 #include "QGSP_INCLXX.hh"
62 #include "QGSP_INCLXX_HP.hh"
63 #include "Shielding.hh"
64 #include "ShieldingLEND.hh"
65 #include "NuBeam.hh"
66 #include "G4ThermalNeutrons.hh"
67
68 #include "G4EmStandardPhysics.hh"
69 #include "G4EmStandardPhysics_option1.hh"
70 #include "G4EmStandardPhysics_option2.hh"
71 #include "G4EmStandardPhysics_option3.hh"
72 #include "G4EmStandardPhysics_option4.hh"
73 #include "G4EmStandardPhysicsGS.hh"
74 #include "G4EmStandardPhysicsSS.hh"
75 #include "G4EmStandardPhysicsWVI.hh"
76 #include "G4EmLowEPPhysics.hh"
77 #include "G4EmLivermorePhysics.hh"
78 #include "G4EmPenelopePhysics.hh"
79 #include "G4PhysListFactoryMessenger.hh"
80 #include "G4UImessenger.hh"
81
82 G4PhysListFactory::G4PhysListFactory(G4int ver
83 : defName("FTFP_BERT"),verbose(ver),theMesse
84 {
85 nlists_hadr = 36;
86 G4String ss[36] = {
87 "FTFP_BERT","FTFP_BERT_TRV","FTFP_BERT_ATL
88 "FTFP_INCLXX","FTFP_INCLXX_HP","FTF_BIC","
89 "QGSP_BERT","QGSP_BERT_HP","QGSP_BIC","QGS
90 "QGSP_FTFP_BERT","QGSP_INCLXX","QGSP_INCLX
91 "Shielding","ShieldingLEND","ShieldingLIQM
92 "Shielding_HP","ShieldingLIQMD_HP","Shield
93 "FTFP_BERT_HPT","FTFP_INCLXX_HPT","QGSP_BE
94 "QGSP_BIC_AllHPT","QGSP_INCLXX_HPT","Shiel
95 "ShieldingM_HPT"};
96 for (std::size_t i=0; i<nlists_hadr; ++i) {
97 listnames_hadr.push_back(ss[i]);
98 }
99
100 nlists_em = 12;
101 G4String s1[12] = {"","_EMV","_EMX","_EMY","
102 "__GS","__SS","_EM0","_WVI","__LE"};
103 for (std::size_t i=0; i<nlists_em; ++i) {
104 listnames_em.push_back(s1[i]);
105 }
106 }
107
108 G4PhysListFactory::~G4PhysListFactory()
109 {
110 delete theMessenger;
111 }
112
113 G4VModularPhysicsList*
114 G4PhysListFactory::ReferencePhysList()
115 {
116 // instantiate PhysList by environment varia
117 G4String name = "";
118 char* path = std::getenv("PHYSLIST");
119 if (path) {
120 name = G4String(path);
121 } else {
122 name = defName;
123 G4cout << "### G4PhysListFactory WARNING:
124 << " environment variable PHYSLIST is not
125 << G4endl
126 << " Default Physics Lists " << name
127 << " is instantiated"
128 << G4endl;
129 }
130 return GetReferencePhysList(name);
131 }
132
133 G4VModularPhysicsList*
134 G4PhysListFactory::GetReferencePhysList(const
135 {
136 // analysis on the string
137 std::size_t n = name.size();
138
139 // last characters in the string
140 size_t em_opt = 0;
141 G4String em_name = "";
142
143 // check EM options
144 if(n > 4) {
145 em_name = name.substr(n - 4, 4);
146 for (std::size_t i=1; i<nlists_em; ++i) {
147 if(listnames_em[i] == em_name) {
148 em_opt = i;
149 n -= 4;
150 break;
151 }
152 }
153 if(0 == em_opt) { em_name = ""; }
154 }
155
156 // hadronic pHysics List
157 G4String had_name = name.substr(0, n);
158
159 if(0 < verbose) {
160 G4cout << "G4PhysListFactory::GetReference
161 << em_name << "> EMoption= " << em_opt <
162 }
163 G4VModularPhysicsList* p = nullptr;
164 if(had_name == "FTFP_BERT") {p
165 else if(had_name == "FTFP_BERT_HP") {p
166 else if(had_name == "FTFP_BERT_TRV") {p
167 else if(had_name == "FTFP_BERT_ATL") {p
168 else if(had_name == "FTFQGSP_BERT") {p
169 else if(had_name == "FTFP_INCLXX") {p
170 else if(had_name == "FTFP_INCLXX_HP") {p
171 else if(had_name == "FTF_BIC") {p
172 else if(had_name == "LBE") {p
173 else if(had_name == "QBBC") {p
174 else if(had_name == "QGSP_BERT") {p
175 else if(had_name == "QGSP_BERT_HP") {p
176 else if(had_name == "QGSP_BIC") {p
177 else if(had_name == "QGSP_BIC_HP") {p
178 else if(had_name == "QGSP_BIC_AllHP") {p
179 else if(had_name == "QGSP_FTFP_BERT") {p
180 else if(had_name == "QGSP_INCLXX") {p
181 else if(had_name == "QGSP_INCLXX_HP") {p
182 else if(had_name == "QGS_BIC") {p
183 else if(had_name == "Shielding") {p
184 else if(had_name == "ShieldingLEND") {p
185 else if(had_name == "ShieldingLIQMD") {p
186 else if(had_name == "ShieldingM") {p
187 else if(had_name == "NuBeam") {p
188 else if(had_name == "Shielding_HP") {p
189 else if(had_name == "ShieldingLIQMD_HP") {p
190 else if(had_name == "ShieldingM_HP") {p
191 else if(had_name == "FTFP_BERT_HPT") {p
192 p
193 else if(had_name == "FTFP_INCLXX_HPT") {p
194 p
195 else if(had_name == "QGSP_BERT_HPT") {p
196 p
197 else if(had_name == "QGSP_BIC_HPT") {p
198 else if(had_name == "QGSP_BIC_AllHPT") {p
199 p
200 else if(had_name == "QGSP_INCLXX_HPT") {p
201 p
202 else if(had_name == "Shielding_HPT") {p
203 p
204 else if(had_name == "ShieldingLIQMD_HPT") {p
205 p
206 else if(had_name == "ShieldingM_HPT") {p
207 p
208 else {
209 p = new FTFP_BERT(verbose);
210 G4ExceptionDescription ed;
211 ed << "PhysicsList " << had_name << " is n
212 << " the default FTFP_BERT is created";
213 G4Exception("G4PhysListFactory: ","pl0003"
214 }
215 if(nullptr != p) {
216 if(0 < em_opt && had_name != "LBE") {
217 if(1 == em_opt) {
218 p->ReplacePhysics(new G4EmStandardPhysics_op
219 } else if(2 == em_opt) {
220 p->ReplacePhysics(new G4EmStandardPhysics_op
221 } else if(3 == em_opt) {
222 p->ReplacePhysics(new G4EmStandardPhysics_op
223 } else if(4 == em_opt) {
224 p->ReplacePhysics(new G4EmStandardPhysics_op
225 } else if(5 == em_opt) {
226 p->ReplacePhysics(new G4EmLivermorePhysics(v
227 } else if(6 == em_opt) {
228 p->ReplacePhysics(new G4EmPenelopePhysics(ve
229 } else if(7 == em_opt) {
230 p->ReplacePhysics(new G4EmStandardPhysicsGS(
231 } else if(8 == em_opt) {
232 p->ReplacePhysics(new G4EmStandardPhysicsSS(
233 } else if(9 == em_opt) {
234 p->ReplacePhysics(new G4EmStandardPhysics(ve
235 } else if(10 == em_opt) {
236 p->ReplacePhysics(new G4EmStandardPhysicsWVI
237 } else if(11 == em_opt) {
238 p->ReplacePhysics(new G4EmLowEPPhysics(verbo
239 }
240 }
241 theMessenger = new G4PhysListFactoryMessen
242 }
243 if(0 < verbose) G4cout << G4endl;
244 return p;
245 }
246
247 G4bool G4PhysListFactory::IsReferencePhysList(
248 {
249 G4bool res = false;
250 std::size_t n = name.size();
251 if(n > 4) {
252 G4String em_name = name.substr(n - 4, 4);
253 for (std::size_t i=1; i<nlists_em; ++i) {
254 if(listnames_em[i] == em_name) {
255 n -= 4;
256 break;
257 }
258 }
259 }
260 G4String had_name = name.substr(0, n);
261 for (std::size_t i=0; i<nlists_hadr; ++i) {
262 if(had_name == listnames_hadr[i]) {
263 res = true;
264 break;
265 }
266 }
267 return res;
268 }
269
270 const std::vector<G4String>&
271 G4PhysListFactory::AvailablePhysLists() const
272 {
273 return listnames_hadr;
274 }
275
276 const std::vector<G4String>&
277 G4PhysListFactory::AvailablePhysListsEM() cons
278 {
279 return listnames_em;
280 }
281
282