| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 // 26 //
27 //-------------------------------------------- 27 //---------------------------------------------------------------------------
28 // 28 //
29 // ClassName: G4PhysListFactory 29 // ClassName: G4PhysListFactory
30 // 30 //
31 // Author: 21 April 2008 V. Ivanchenko 31 // Author: 21 April 2008 V. Ivanchenko
32 // 32 //
33 // Modified: 33 // Modified:
34 // 34 //
35 // 2023.04.12 A.Ribon added _HPT variants (i.e 35 // 2023.04.12 A.Ribon added _HPT variants (i.e. HP + thermal scattering)
36 // 2014.08.05 K.L.Genser used provision for Ha 36 // 2014.08.05 K.L.Genser used provision for Hadronic Physics Variant M in
37 // Shielding for ShieldingM 37 // Shielding for ShieldingM
38 // 38 //
39 //-------------------------------------------- 39 //----------------------------------------------------------------------------
40 // 40 //
41 41
42 #include "G4PhysListFactory.hh" 42 #include "G4PhysListFactory.hh"
43 #include "FTFP_BERT.hh" 43 #include "FTFP_BERT.hh"
44 #include "FTFP_BERT_HP.hh" 44 #include "FTFP_BERT_HP.hh"
45 #include "FTFP_BERT_TRV.hh" 45 #include "FTFP_BERT_TRV.hh"
46 #include "FTFP_BERT_ATL.hh" 46 #include "FTFP_BERT_ATL.hh"
47 #include "FTFQGSP_BERT.hh" 47 #include "FTFQGSP_BERT.hh"
48 #include "FTFP_INCLXX.hh" 48 #include "FTFP_INCLXX.hh"
49 #include "FTFP_INCLXX_HP.hh" 49 #include "FTFP_INCLXX_HP.hh"
50 #include "FTF_BIC.hh" 50 #include "FTF_BIC.hh"
51 #include "LBE.hh" 51 #include "LBE.hh"
52 #include "QBBC.hh" 52 #include "QBBC.hh"
53 #include "QGSP_BERT.hh" 53 #include "QGSP_BERT.hh"
54 #include "QGSP_BERT_HP.hh" 54 #include "QGSP_BERT_HP.hh"
55 #include "QGSP_BIC.hh" 55 #include "QGSP_BIC.hh"
56 #include "QGSP_BIC_HP.hh" 56 #include "QGSP_BIC_HP.hh"
57 #include "QGSP_BIC_HPT.hh" 57 #include "QGSP_BIC_HPT.hh"
58 #include "QGSP_BIC_AllHP.hh" 58 #include "QGSP_BIC_AllHP.hh"
59 #include "QGSP_FTFP_BERT.hh" 59 #include "QGSP_FTFP_BERT.hh"
60 #include "QGS_BIC.hh" 60 #include "QGS_BIC.hh"
61 #include "QGSP_INCLXX.hh" 61 #include "QGSP_INCLXX.hh"
62 #include "QGSP_INCLXX_HP.hh" 62 #include "QGSP_INCLXX_HP.hh"
63 #include "Shielding.hh" 63 #include "Shielding.hh"
64 #include "ShieldingLEND.hh" 64 #include "ShieldingLEND.hh"
65 #include "NuBeam.hh" 65 #include "NuBeam.hh"
66 #include "G4ThermalNeutrons.hh" 66 #include "G4ThermalNeutrons.hh"
67 67
68 #include "G4EmStandardPhysics.hh" 68 #include "G4EmStandardPhysics.hh"
69 #include "G4EmStandardPhysics_option1.hh" 69 #include "G4EmStandardPhysics_option1.hh"
70 #include "G4EmStandardPhysics_option2.hh" 70 #include "G4EmStandardPhysics_option2.hh"
71 #include "G4EmStandardPhysics_option3.hh" 71 #include "G4EmStandardPhysics_option3.hh"
72 #include "G4EmStandardPhysics_option4.hh" 72 #include "G4EmStandardPhysics_option4.hh"
73 #include "G4EmStandardPhysicsGS.hh" 73 #include "G4EmStandardPhysicsGS.hh"
74 #include "G4EmStandardPhysicsSS.hh" 74 #include "G4EmStandardPhysicsSS.hh"
75 #include "G4EmStandardPhysicsWVI.hh" 75 #include "G4EmStandardPhysicsWVI.hh"
76 #include "G4EmLowEPPhysics.hh" 76 #include "G4EmLowEPPhysics.hh"
77 #include "G4EmLivermorePhysics.hh" 77 #include "G4EmLivermorePhysics.hh"
78 #include "G4EmPenelopePhysics.hh" 78 #include "G4EmPenelopePhysics.hh"
79 #include "G4PhysListFactoryMessenger.hh" 79 #include "G4PhysListFactoryMessenger.hh"
80 #include "G4UImessenger.hh" 80 #include "G4UImessenger.hh"
81 81
82 G4PhysListFactory::G4PhysListFactory(G4int ver 82 G4PhysListFactory::G4PhysListFactory(G4int ver)
83 : defName("FTFP_BERT"),verbose(ver),theMesse 83 : defName("FTFP_BERT"),verbose(ver),theMessenger(nullptr)
84 { 84 {
85 nlists_hadr = 36; << 85 nlists_hadr = 33;
86 G4String ss[36] = { << 86 G4String ss[33] = {
87 "FTFP_BERT","FTFP_BERT_TRV","FTFP_BERT_ATL 87 "FTFP_BERT","FTFP_BERT_TRV","FTFP_BERT_ATL","FTFP_BERT_HP","FTFQGSP_BERT",
88 "FTFP_INCLXX","FTFP_INCLXX_HP","FTF_BIC"," << 88 "FTFP_INCLXX","FTFP_INCLXX_HP","FTF_BIC", "LBE","QBBC",
89 "QGSP_BERT","QGSP_BERT_HP","QGSP_BIC","QGS 89 "QGSP_BERT","QGSP_BERT_HP","QGSP_BIC","QGSP_BIC_HP","QGSP_BIC_AllHP",
90 "QGSP_FTFP_BERT","QGSP_INCLXX","QGSP_INCLX 90 "QGSP_FTFP_BERT","QGSP_INCLXX","QGSP_INCLXX_HP","QGS_BIC",
91 "Shielding","ShieldingLEND","ShieldingLIQM << 91 "Shielding","ShieldingLEND","ShieldingM","NuBeam",
92 "Shielding_HP","ShieldingLIQMD_HP","Shield << 92 "Shielding_HP","ShieldingM_HP",
93 "FTFP_BERT_HPT","FTFP_INCLXX_HPT","QGSP_BE << 93 "FTFP_BERT_HPT", "FTFP_INCLXX_HPT", "QGSP_BERT_HPT", "QGSP_BIC_HPT",
94 "QGSP_BIC_AllHPT","QGSP_INCLXX_HPT","Shiel << 94 "QGSP_BIC_AllHPT", "QGSP_INCLXX_HPT", "Shielding_HPT", "ShieldingM_HPT" };
95 "ShieldingM_HPT"}; << 95 for(size_t i=0; i<nlists_hadr; ++i) {
96 for (std::size_t i=0; i<nlists_hadr; ++i) { <<
97 listnames_hadr.push_back(ss[i]); 96 listnames_hadr.push_back(ss[i]);
98 } 97 }
99 98
100 nlists_em = 12; 99 nlists_em = 12;
101 G4String s1[12] = {"","_EMV","_EMX","_EMY"," 100 G4String s1[12] = {"","_EMV","_EMX","_EMY","_EMZ","_LIV","_PEN",
102 "__GS","__SS","_EM0","_WVI","__LE"}; 101 "__GS","__SS","_EM0","_WVI","__LE"};
103 for (std::size_t i=0; i<nlists_em; ++i) { << 102 for(size_t i=0; i<nlists_em; ++i) {
104 listnames_em.push_back(s1[i]); 103 listnames_em.push_back(s1[i]);
105 } 104 }
106 } 105 }
107 106
108 G4PhysListFactory::~G4PhysListFactory() 107 G4PhysListFactory::~G4PhysListFactory()
109 { 108 {
110 delete theMessenger; 109 delete theMessenger;
111 } 110 }
112 111
113 G4VModularPhysicsList* 112 G4VModularPhysicsList*
114 G4PhysListFactory::ReferencePhysList() 113 G4PhysListFactory::ReferencePhysList()
115 { 114 {
116 // instantiate PhysList by environment varia 115 // instantiate PhysList by environment variable "PHYSLIST"
117 G4String name = ""; 116 G4String name = "";
118 char* path = std::getenv("PHYSLIST"); 117 char* path = std::getenv("PHYSLIST");
119 if (path) { 118 if (path) {
120 name = G4String(path); 119 name = G4String(path);
121 } else { 120 } else {
122 name = defName; 121 name = defName;
123 G4cout << "### G4PhysListFactory WARNING: 122 G4cout << "### G4PhysListFactory WARNING: "
124 << " environment variable PHYSLIST is not 123 << " environment variable PHYSLIST is not defined"
125 << G4endl 124 << G4endl
126 << " Default Physics Lists " << name 125 << " Default Physics Lists " << name
127 << " is instantiated" 126 << " is instantiated"
128 << G4endl; 127 << G4endl;
129 } 128 }
130 return GetReferencePhysList(name); 129 return GetReferencePhysList(name);
131 } 130 }
132 131
133 G4VModularPhysicsList* 132 G4VModularPhysicsList*
134 G4PhysListFactory::GetReferencePhysList(const 133 G4PhysListFactory::GetReferencePhysList(const G4String& name)
135 { 134 {
136 // analysis on the string 135 // analysis on the string
137 std::size_t n = name.size(); << 136 size_t n = name.size();
138 137
139 // last characters in the string 138 // last characters in the string
140 size_t em_opt = 0; 139 size_t em_opt = 0;
141 G4String em_name = ""; 140 G4String em_name = "";
142 141
143 // check EM options 142 // check EM options
144 if(n > 4) { 143 if(n > 4) {
145 em_name = name.substr(n - 4, 4); 144 em_name = name.substr(n - 4, 4);
146 for (std::size_t i=1; i<nlists_em; ++i) { << 145 for(size_t i=1; i<nlists_em; ++i) {
147 if(listnames_em[i] == em_name) { 146 if(listnames_em[i] == em_name) {
148 em_opt = i; 147 em_opt = i;
149 n -= 4; 148 n -= 4;
150 break; 149 break;
151 } 150 }
152 } 151 }
153 if(0 == em_opt) { em_name = ""; } 152 if(0 == em_opt) { em_name = ""; }
154 } 153 }
155 154
156 // hadronic pHysics List 155 // hadronic pHysics List
157 G4String had_name = name.substr(0, n); 156 G4String had_name = name.substr(0, n);
158 157
159 if(0 < verbose) { 158 if(0 < verbose) {
160 G4cout << "G4PhysListFactory::GetReference 159 G4cout << "G4PhysListFactory::GetReferencePhysList <" << had_name
161 << em_name << "> EMoption= " << em_opt < 160 << em_name << "> EMoption= " << em_opt << G4endl;
162 } 161 }
163 G4VModularPhysicsList* p = nullptr; 162 G4VModularPhysicsList* p = nullptr;
164 if(had_name == "FTFP_BERT") {p << 163 if(had_name == "FTFP_BERT") {p = new FTFP_BERT(verbose);}
165 else if(had_name == "FTFP_BERT_HP") {p << 164 else if(had_name == "FTFP_BERT_HP") {p = new FTFP_BERT_HP(verbose);}
166 else if(had_name == "FTFP_BERT_TRV") {p << 165 else if(had_name == "FTFP_BERT_TRV") {p = new FTFP_BERT_TRV(verbose);}
167 else if(had_name == "FTFP_BERT_ATL") {p << 166 else if(had_name == "FTFP_BERT_ATL") {p = new FTFP_BERT_ATL(verbose);}
168 else if(had_name == "FTFQGSP_BERT") {p << 167 else if(had_name == "FTFQGSP_BERT") {p = new FTFQGSP_BERT(verbose);}
169 else if(had_name == "FTFP_INCLXX") {p << 168 else if(had_name == "FTFP_INCLXX") {p = new FTFP_INCLXX(verbose);}
170 else if(had_name == "FTFP_INCLXX_HP") {p << 169 else if(had_name == "FTFP_INCLXX_HP") {p = new FTFP_INCLXX_HP(verbose);}
171 else if(had_name == "FTF_BIC") {p << 170 else if(had_name == "FTF_BIC") {p = new FTF_BIC(verbose);}
172 else if(had_name == "LBE") {p << 171 else if(had_name == "LBE") {p = new LBE();}
173 else if(had_name == "QBBC") {p << 172 else if(had_name == "QBBC") {p = new QBBC(verbose);}
174 else if(had_name == "QGSP_BERT") {p << 173 else if(had_name == "QGSP_BERT") {p = new QGSP_BERT(verbose);}
175 else if(had_name == "QGSP_BERT_HP") {p << 174 else if(had_name == "QGSP_BERT_HP") {p = new QGSP_BERT_HP(verbose);}
176 else if(had_name == "QGSP_BIC") {p << 175 else if(had_name == "QGSP_BIC") {p = new QGSP_BIC(verbose);}
177 else if(had_name == "QGSP_BIC_HP") {p << 176 else if(had_name == "QGSP_BIC_HP") {p = new QGSP_BIC_HP(verbose);}
178 else if(had_name == "QGSP_BIC_AllHP") {p << 177 else if(had_name == "QGSP_BIC_AllHP") {p = new QGSP_BIC_AllHP(verbose);}
179 else if(had_name == "QGSP_FTFP_BERT") {p << 178 else if(had_name == "QGSP_FTFP_BERT") {p = new QGSP_FTFP_BERT(verbose);}
180 else if(had_name == "QGSP_INCLXX") {p << 179 else if(had_name == "QGSP_INCLXX") {p = new QGSP_INCLXX(verbose);}
181 else if(had_name == "QGSP_INCLXX_HP") {p << 180 else if(had_name == "QGSP_INCLXX_HP") {p = new QGSP_INCLXX_HP(verbose);}
182 else if(had_name == "QGS_BIC") {p << 181 else if(had_name == "QGS_BIC") {p = new QGS_BIC(verbose);}
183 else if(had_name == "Shielding") {p << 182 else if(had_name == "Shielding") {p = new Shielding(verbose);}
184 else if(had_name == "ShieldingLEND") {p << 183 else if(had_name == "ShieldingLEND") {p = new ShieldingLEND(verbose);}
185 else if(had_name == "ShieldingLIQMD") {p << 184 else if(had_name == "ShieldingM") {p = new Shielding(verbose,"HP","M");}
186 else if(had_name == "ShieldingM") {p << 185 else if(had_name == "NuBeam") {p = new NuBeam(verbose);}
187 else if(had_name == "NuBeam") {p << 186 else if(had_name == "Shielding_HP") {p = new Shielding(verbose);}
188 else if(had_name == "Shielding_HP") {p << 187 else if(had_name == "ShieldingM_HP") {p = new Shielding(verbose,"HP","M");}
189 else if(had_name == "ShieldingLIQMD_HP") {p << 188 else if(had_name == "FTFP_BERT_HPT") {p = new FTFP_BERT_HP(verbose);
190 else if(had_name == "ShieldingM_HP") {p << 189 p->RegisterPhysics(new G4ThermalNeutrons);}
191 else if(had_name == "FTFP_BERT_HPT") {p << 190 else if(had_name == "FTFP_INCLXX_HPT"){p = new FTFP_INCLXX_HP(verbose);
192 p << 191 p->RegisterPhysics(new G4ThermalNeutrons);}
193 else if(had_name == "FTFP_INCLXX_HPT") {p << 192 else if(had_name == "QGSP_BERT_HPT") {p = new QGSP_BERT_HP(verbose);
194 p << 193 p->RegisterPhysics(new G4ThermalNeutrons);}
195 else if(had_name == "QGSP_BERT_HPT") {p << 194 else if(had_name == "QGSP_BIC_HPT") {p = new QGSP_BIC_HPT(verbose);}
196 p << 195 else if(had_name == "QGSP_BIC_AllHPT"){p = new QGSP_BIC_AllHP(verbose);
197 else if(had_name == "QGSP_BIC_HPT") {p << 196 p->RegisterPhysics(new G4ThermalNeutrons);}
198 else if(had_name == "QGSP_BIC_AllHPT") {p << 197 else if(had_name == "QGSP_INCLXX_HPT"){p = new QGSP_INCLXX_HP(verbose);
199 p << 198 p->RegisterPhysics(new G4ThermalNeutrons);}
200 else if(had_name == "QGSP_INCLXX_HPT") {p << 199 else if(had_name == "Shielding_HPT") {p = new Shielding(verbose);
201 p << 200 p->RegisterPhysics(new G4ThermalNeutrons);}
202 else if(had_name == "Shielding_HPT") {p << 201 else if(had_name == "ShieldingM_HPT") {p = new Shielding(verbose,"HP","M");
203 p << 202 p->RegisterPhysics(new G4ThermalNeutrons);}
204 else if(had_name == "ShieldingLIQMD_HPT") {p <<
205 p <<
206 else if(had_name == "ShieldingM_HPT") {p <<
207 p <<
208 else { 203 else {
209 p = new FTFP_BERT(verbose); 204 p = new FTFP_BERT(verbose);
210 G4ExceptionDescription ed; 205 G4ExceptionDescription ed;
211 ed << "PhysicsList " << had_name << " is n 206 ed << "PhysicsList " << had_name << " is not known;"
212 << " the default FTFP_BERT is created"; 207 << " the default FTFP_BERT is created";
213 G4Exception("G4PhysListFactory: ","pl0003" 208 G4Exception("G4PhysListFactory: ","pl0003",JustWarning,ed,"");
214 } 209 }
215 if(nullptr != p) { 210 if(nullptr != p) {
216 if(0 < em_opt && had_name != "LBE") { 211 if(0 < em_opt && had_name != "LBE") {
217 if(1 == em_opt) { 212 if(1 == em_opt) {
218 p->ReplacePhysics(new G4EmStandardPhysics_op 213 p->ReplacePhysics(new G4EmStandardPhysics_option1(verbose));
219 } else if(2 == em_opt) { 214 } else if(2 == em_opt) {
220 p->ReplacePhysics(new G4EmStandardPhysics_op 215 p->ReplacePhysics(new G4EmStandardPhysics_option2(verbose));
221 } else if(3 == em_opt) { 216 } else if(3 == em_opt) {
222 p->ReplacePhysics(new G4EmStandardPhysics_op 217 p->ReplacePhysics(new G4EmStandardPhysics_option3(verbose));
223 } else if(4 == em_opt) { 218 } else if(4 == em_opt) {
224 p->ReplacePhysics(new G4EmStandardPhysics_op 219 p->ReplacePhysics(new G4EmStandardPhysics_option4(verbose));
225 } else if(5 == em_opt) { 220 } else if(5 == em_opt) {
226 p->ReplacePhysics(new G4EmLivermorePhysics(v 221 p->ReplacePhysics(new G4EmLivermorePhysics(verbose));
227 } else if(6 == em_opt) { 222 } else if(6 == em_opt) {
228 p->ReplacePhysics(new G4EmPenelopePhysics(ve 223 p->ReplacePhysics(new G4EmPenelopePhysics(verbose));
229 } else if(7 == em_opt) { 224 } else if(7 == em_opt) {
230 p->ReplacePhysics(new G4EmStandardPhysicsGS( 225 p->ReplacePhysics(new G4EmStandardPhysicsGS(verbose));
231 } else if(8 == em_opt) { 226 } else if(8 == em_opt) {
232 p->ReplacePhysics(new G4EmStandardPhysicsSS( 227 p->ReplacePhysics(new G4EmStandardPhysicsSS(verbose));
233 } else if(9 == em_opt) { 228 } else if(9 == em_opt) {
234 p->ReplacePhysics(new G4EmStandardPhysics(ve 229 p->ReplacePhysics(new G4EmStandardPhysics(verbose));
235 } else if(10 == em_opt) { 230 } else if(10 == em_opt) {
236 p->ReplacePhysics(new G4EmStandardPhysicsWVI 231 p->ReplacePhysics(new G4EmStandardPhysicsWVI(verbose));
237 } else if(11 == em_opt) { 232 } else if(11 == em_opt) {
238 p->ReplacePhysics(new G4EmLowEPPhysics(verbo 233 p->ReplacePhysics(new G4EmLowEPPhysics(verbose));
239 } 234 }
240 } 235 }
241 theMessenger = new G4PhysListFactoryMessen 236 theMessenger = new G4PhysListFactoryMessenger(p);
242 } 237 }
243 if(0 < verbose) G4cout << G4endl; 238 if(0 < verbose) G4cout << G4endl;
244 return p; 239 return p;
245 } 240 }
246 241
247 G4bool G4PhysListFactory::IsReferencePhysList( 242 G4bool G4PhysListFactory::IsReferencePhysList(const G4String& name) const
248 { 243 {
249 G4bool res = false; 244 G4bool res = false;
250 std::size_t n = name.size(); << 245 size_t n = name.size();
251 if(n > 4) { 246 if(n > 4) {
252 G4String em_name = name.substr(n - 4, 4); 247 G4String em_name = name.substr(n - 4, 4);
253 for (std::size_t i=1; i<nlists_em; ++i) { << 248 for(size_t i=1; i<nlists_em; ++i) {
254 if(listnames_em[i] == em_name) { 249 if(listnames_em[i] == em_name) {
255 n -= 4; 250 n -= 4;
256 break; 251 break;
257 } 252 }
258 } 253 }
259 } 254 }
260 G4String had_name = name.substr(0, n); 255 G4String had_name = name.substr(0, n);
261 for (std::size_t i=0; i<nlists_hadr; ++i) { << 256 for(size_t i=0; i<nlists_hadr; ++i) {
262 if(had_name == listnames_hadr[i]) { 257 if(had_name == listnames_hadr[i]) {
263 res = true; 258 res = true;
264 break; 259 break;
265 } 260 }
266 } 261 }
267 return res; 262 return res;
268 } 263 }
269 264
270 const std::vector<G4String>& 265 const std::vector<G4String>&
271 G4PhysListFactory::AvailablePhysLists() const 266 G4PhysListFactory::AvailablePhysLists() const
272 { 267 {
273 return listnames_hadr; 268 return listnames_hadr;
274 } 269 }
275 270
276 const std::vector<G4String>& 271 const std::vector<G4String>&
277 G4PhysListFactory::AvailablePhysListsEM() cons 272 G4PhysListFactory::AvailablePhysListsEM() const
278 { 273 {
279 return listnames_em; 274 return listnames_em;
280 } 275 }
281 276
282 277