| 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 // G4PhysicsModelCatalog class implementation
27 //
28 // Author: M.Asai (SLAC), 26 September 2013
29 //
30 // Revised in August 2021 by A.Ribon (CERN).
31 // -------------------------------------------
32
33 #include "G4PhysicsModelCatalog.hh"
34 #include "G4Threading.hh"
35
36 G4bool G4PhysicsModelCatalog::isInitialized =
37 std::vector< G4int >* G4PhysicsModelCatalog
38 std::vector< G4String >* G4PhysicsModelCatalog
39
40 // -------------------------------------------
41 void G4PhysicsModelCatalog::Initialize() {
42 if(isInitialized)
43 {
44 return;
45 }
46 if ( theVectorOfModelIDs == nullptr && the
47 static std::vector< G4int > aVectorOfInts;
48 theVectorOfModelIDs = &aVectorOfInts;
49 static std::vector< G4String > aVectorOfSt
50 theVectorOfModelNames = &aVectorOfStrings;
51
52 // NOTE:
53 // The goal is that, starting from Geant4
54 // identifiers (modelID, index, name) rema
55 // physics list, application, and version
56 // Therefore, after Geant4 11.0, you can o
57 // a new model is added, or when a pre-exi
58 // in this catalogue is included - at the
59 // (rather than inserting it in the middle
60 //
61 // For the modelID values, these are the c
62 // - Values below 1000 are excluded becaus
63 // sub-process ID values.
64 // - Whenever resonable, modelID values sh
65 // them, to allow for eventual, future v
66 // tunings of the same model - to keep m
67 // original model.
68 // - modelID values are between 10'000 and
69 // subdivision in 3 categories (identifi
70 // digit):
71 // * Values between 10'000 and 19'999 a
72 // * Values between 20'000 and 29'999 a
73 // * Values between 30'000 and 39'999 a
74 // (i.e. included neither in EM model
75 // Note that larger values of modelID ar
76 // handle nor less computing performant
77 // (we remind that, for plotting, the in
78 // its modelID, should be conveniently
79 // the modelID is recommended).
80
81 // =======================================
82 // ================= 1st EM MODELS : from
83 // =======================================
84
85 InsertModel( 10000, "model_EM" );
86
87 // e- production
88 InsertModel( 10010, "model_DeltaElectron"
89 InsertModel( 10011, "model_DeltaEBelowCut"
90 InsertModel( 10012, "model_PhotoElectron"
91 InsertModel( 10013, "model_ComptonElectron
92 InsertModel( 10014, "model_TripletElectron
93
94 // gamma production
95 InsertModel( 10020, "model_Bremsstrahlung"
96 InsertModel( 10021, "model_SplitBremsstrah
97 InsertModel( 10022, "model_ComptonGamma" )
98 InsertModel( 10023, "model_Annihilation" )
99 InsertModel( 10024, "model_TripletGamma" )
100 InsertModel( 10025, "model_GammaGammaEntan
101
102 // e+e- pair production
103 InsertModel( 10030, "model_EplusEminisPair
104
105 // atomic de-excitation
106 InsertModel( 10040, "model_Fluorescence" )
107 InsertModel( 10041, "model_gammaPIXE" );
108 InsertModel( 10050, "model_AugerElectron"
109 InsertModel( 10051, "model_ePIXE" );
110
111 // recoil ions
112 InsertModel( 10060, "model_IonRecoil" );
113
114 // DNA models
115 InsertModel( 11000, "model_DNA" );
116 InsertModel( 11001, "model_Ritchie1994eSol
117 InsertModel( 11002, "model_Terrisol1990eSo
118 InsertModel( 11003, "model_Meesungnoen2002
119 InsertModel( 11004, "model_Kreipl2009eSolv
120 InsertModel( 11005, "model_MeesungnoenSoli
121
122 // =======================================
123 // ============= 2nd HADRONIC MODELS : fro
124 // =======================================
125
126 // ---------------------------------------
127 // --- Gamma- , Lepto- , Neutrino-nuclear
128 // ---------------------------------------
129 // - 20'000 - 20'099 : Electromagnetic
130 // - 20'100 - 20'199 : Gamma-nuclear m
131 // - 20'200 - 20-299 : Electron/positr
132 // - 20'300 - 20'399 : Muon-nuclear mo
133 // - 20'400 - 20'499 : Tau-nuclear mod
134 // - 20'500 - 20'999 : Neutrino models
135 // ...
136
137 // --- EM dissociation models: 20'000 - 20
138
139 // Class G4EMDissociation
140 InsertModel( 20000, "model_projectileEMDis
141 InsertModel( 20001, "model_targetEMDissoci
142
143 // --- Gamma-nuclear models: 20'100 - 20'1
144
145 // Class G4LENDorBERTModel
146 InsertModel( 20100, "model_LENDorBERTModel
147
148 // Class G4LowEGammaNuclearModel
149 InsertModel( 20150, "model_GammaNPreco" );
150
151 // --- Charged-lepton - nuclear models: 20
152
153 // Class G4ElectroVDNuclearModel
154 InsertModel( 20200, "model_G4ElectroVDNucl
155
156 // Class G4MuonVDNuclearModel
157 InsertModel( 20300, "model_G4MuonVDNuclear
158
159 // --- Neutrino models: 20'500 - 20'999 --
160
161 // Class G4NeutrinoElectronCcModel
162 InsertModel( 20510, "model_nu-e-inelastic"
163
164 // Class G4NeutrinoNucleusModel
165 InsertModel( 20520, "model_neutrino-nucleu
166
167 // The following classes derives from G4Ne
168 // Class G4ANuElNucleusCcModel
169 InsertModel( 20530, "model_ANuElNuclCcMode
170 // Class G4ANuElNucleusNcModel
171 InsertModel( 20540, "model_ANuElNuclNcMode
172 // Class G4ANuMuNucleusCcModel
173 InsertModel( 20550, "model_ANuMuNuclCcMode
174 // Class G4ANuMuNucleusNcModel
175 InsertModel( 20560, "model_ANuMuNuclNcMode
176 // Class G4NuElNucleusCcModel
177 InsertModel( 20570, "model_NuElNuclCcModel
178 // Class G4NuElNucleusNcModel
179 InsertModel( 20580, "model_NuElNuclNcModel
180 // Class G4NuMuNucleusCcModel
181 InsertModel( 20590, "model_NuMuNuclCcModel
182 // Class G4NuMuNucleusNcModel
183 InsertModel( 20600, "model_NuMuNuclNcModel
184
185 // ---------------------------------------
186 // --- Elastic, Charge-Exchange, Quasi-Ela
187 // ---------------------------------------
188 // - 21'000 - 21'199 : Elastic
189 // - 21'200 - 21'299 : Charge-Exchange
190 // - 21'300 - 21'499 : Quasi-Elastic
191 // - 21'500 - 21'999 : specialized Dif
192 // ---------------------------------------
193
194 // --- Elastic models: 21'000 - 21'199 ---
195
196 // Class G4HadronElastic
197 InsertModel( 21000, "model_hElasticLHEP" )
198 // Class G4AntiNuclElastic
199 InsertModel( 21010, "model_AntiAElastic" )
200 // Class G4ChipsElasticModel
201 InsertModel( 21020, "model_hElasticCHIPS"
202 // Class G4DiffuseElastic
203 InsertModel( 21030, "model_DiffuseElastic"
204 // Class G4DiffuseElasticV2
205 InsertModel( 21040, "model_DiffuseElasticV
206 // Class G4NuclNuclDiffuseElastic
207 InsertModel( 21050, "model_NNDiffuseElasti
208 // Class G4ElasticHadrNucleusHE
209 InsertModel( 21060, "model_hElasticGlauber
210 // Class G4hhElastic
211 InsertModel( 21070, "model_HadrHadrElastic
212 // Class G4LowEHadronElastic
213 InsertModel( 21080, "model_hLowEElastic" )
214 // Class G4LEHadronProtonElastic
215 InsertModel( 21090, "model_G4LEHadronProto
216 // Class G4LEnp
217 InsertModel( 21100, "model_G4LEnp" );
218 // Class G4LEpp
219 InsertModel( 21110, "model_G4LEpp" );
220 // Class G4NeutronElectronElModel
221 InsertModel( 21120, "model_n-e-elastic" );
222 // Class G4NeutrinoElectronNcModel
223 InsertModel( 21130, "model_nu-e-elastic" )
224
225 // --- Charge exchange : 21'200 - 21'299 -
226
227 // Class: G4ChargeExchange
228 InsertModel( 21200, "model_ChargeExchange"
229
230 // --- Quasi-Elastic : 21'300 - 21'499 ---
231
232 // Class: G4QuasiElasticChannel (which use
233 InsertModel( 21300, "model_QuasiElastic" )
234
235 // --- Special diffraction generators : 21
236
237 // Class: G4LMsdGenerator
238 InsertModel( 21500, "model_LMsdGenerator"
239
240 // ---------------------------------------
241 // --- High energy-models (e.g. string mod
242 // ---------------------------------------
243 // - 22'000 - 22'099 : G4TheoFSGenerat
244 // - 22'100 - 22'199 : FTF
245 // - 22'200 - 22'299 : QGS
246 // - ...
247 // For gamma - nuclear : QGS strin
248 // For e-/e+/mu-/mu+ - nuclear : FTF strin
249 InsertModel( 22000, "model_TheoFSGenerator
250 // FTF string formation + Lund string frag
251 InsertModel( 22100, "model_FTFP" );
252 // FTF string formation + Lund string frag
253 InsertModel( 22150, "model_FTFB" );
254 // FTF string formation + QGS string fragm
255 InsertModel( 22175, "model_FTFQGSP" );
256 // QGS string formation + QGS string fragm
257 InsertModel( 22200, "model_QGSP" );
258 // QGS string formation + QGS string fragm
259 InsertModel( 22250, "model_QGSB" );
260
261 // ---------------------------------------
262 // --- Intermediate energy models : 23'000
263 // ---------------------------------------
264 // - 23'000 - 23'099 : BERT
265 // - 23'100 - 23'199 : BIC
266 // - 23'200 - 23'299 : INCL
267 // - 23'300 - 23'399 : QMD
268 // ...
269 // Class: G4CascadeInterface
270 InsertModel( 23000, "model_BertiniCascade"
271
272 // The names are similar, but not identica
273 // Class: G4BinaryCascade
274 InsertModel( 23100, "model_G4BinaryCascade
275 // Class: G4BinaryLightIonReaction
276 InsertModel( 23110, "model_G4BinaryLightIo
277
278 // Class: G4INCLXXInterface
279 InsertModel( 23200, "model_INCLXXCascade"
280
281 // Class: G4QMDReaction
282 InsertModel( 23300, "model_QMDModel" );
283
284 // ---------------------------------------
285 // --- Pre-equilibrium/De-excitation model
286 // ---------------------------------------
287 // - Pre-equilibrium : 24'000 - 24'099
288 // * 24'000 - 24'049 : precompound
289 // * 24'050 - 24'099 : internal BER
290 // - de-excitation : 24'100 - 24'999
291 // * 24'100 - 24'149 : Evaporation
292 // * 24'150 - 24'199 : Photon evapo
293 // * 24'200 - 24'299 : GEM evaporat
294 // * 24'300 - 24'349 : Fermi BreakU
295 // * 24'350 - 24'399 : Multifragmen
296 // * 24'400 - 24'449 : Ablation
297 // * 24'450 - 24'499 : Fission
298 // * 24'500 - 24'599 : ABLA
299 // * 24'600 - 24'699 : internal BER
300 // ...
301
302 // --- Pre-equilibrium: 24'000 - 24'099 --
303
304 // Class: G4PreCompoundModel
305 InsertModel( 24000, "model_PRECO" );
306
307 // Class: G4LowEIonFragmentation
308 InsertModel( 24010, "model_LowEIonPreco" )
309
310 // Class: G4NonEquilibriumEvaporator (i.e.
311 InsertModel( 24050, "model_G4NonEquilibriu
312
313 // --- De-excitation: 24'100 - 24'999 ---
314
315 // --- Evaporation: 24'100- 24'149 --
316
317 // Class: G4EvaporationChannel
318 InsertModel( 24100, "model_G4EvaporationCh
319
320 // Class: G4UnstableFragmentBreakUp
321 InsertModel( 24120, "model_G4UnstableFragm
322
323 // --- Photon-Evaporation: 24'150 - 24
324
325 // Class: G4PhotonEvaporation
326 InsertModel( 24150, "model_G4PhotonEvapora
327
328 // Class: G4NeutronRadCapture
329 InsertModel( 24160, "model_nRadCapture" );
330
331 // --- GEM evaporation : 24'200 - 24'2
332
333 // Class: G4GEMChannel
334 InsertModel( 24200, "model_G4GEMChannel" )
335
336 // Class: G4GEMChannelVI
337 InsertModel( 24210, "model_G4GEMChannelVI"
338
339 // --- Fermi BreakUp : 24'300 - 24'349
340
341 // Class: G4FermiBreakUpVI
342 InsertModel( 24300, "model_G4FermiBreakUpV
343
344 // --- Multifragmentation : 24'350 - 2
345
346 // Class: G4StatMF
347 InsertModel( 24350, "model_G4StatMF" );
348
349 // --- Ablation : 24'400 - 24'449 ---
350
351 // Class: G4WilsonAblationModel
352 InsertModel( 24400, "model_G4WilsonAblatio
353
354 // --- Fission : 24'450 - 24'499 ---
355
356 // Class: G4CompetitiveFission
357 InsertModel( 24450, "model_G4CompetitiveFi
358
359 // --- ABLA : 24'500 - 24'599 ---
360
361 // Class: G4AblaInterface
362 InsertModel( 24500, "model_ABLA" );
363
364 // --- internal BERT de-excitation : 2
365
366 // Class: G4EquilibriumEvaporator
367 InsertModel( 24600, "model_G4EquilibriumEv
368
369 // --- Other types of de-excitation :
370
371 // ...
372
373 // ---------------------------------------
374 // --- Low-energy data-driven : 25'000 - 2
375 // ---------------------------------------
376 // - 25'000 - 25'199 : ParticleHP
377 // - 25'200 - 25'299 : LEND
378 // ...
379 // - 25'500 - 25'999 : RadioactiveDeca
380
381 // --- ParticleHP : 25'000 - 25'199 ---
382
383 // Classes: G4ParticleHPCapture , G4Partic
384 InsertModel( 25000, "model_NeutronHPCaptur
385
386 // Classes: G4ParticleHPElastic , G4Partic
387 InsertModel( 25010, "model_NeutronHPElasti
388
389 // Classes: G4ParticleHPFission , G4Partic
390 InsertModel( 25020, "model_NeutronHPFissio
391
392 // Inelastic: the following classes inheri
393 // G4ParticleHPInelasticCompFS,
394
395 // Class G4ParticleHPNInelasticFS
396 InsertModel( 25030, "model_G4ParticleHPNIn
397 // Class model_G4ParticleHPNXInelasticFS
398 InsertModel( 25031, "model_G4ParticleHPNXI
399 // Class G4ParticleHP2NDInelasticFS
400 InsertModel( 25032, "model_G4ParticleHP2ND
401 // Class G4ParticleHP2NInelasticFS
402 InsertModel( 25033, "model_G4ParticleHP2NI
403 // Class G4ParticleHP3NInelasticFS
404 InsertModel( 25034, "model_G4ParticleHP3NI
405 // Class G4ParticleHPNAInelasticFS
406 InsertModel( 25035, "model_G4ParticleHPNAI
407 // Class G4ParticleHPN3AInelasticFS
408 InsertModel( 25036, "model_G4ParticleHPN3A
409 // Class G4ParticleHP2NAInelasticFS
410 InsertModel( 25037, "model_G4ParticleHP2NA
411 // Class G4ParticleHP3NAInelasticFS
412 InsertModel( 25038, "model_G4ParticleHP3NA
413 // Class G4ParticleHPNPInelasticFS
414 InsertModel( 25039, "model_G4ParticleHPNPI
415 // Class G4ParticleHPN2AInelasticFS
416 InsertModel( 25040, "model_G4ParticleHPN2A
417 // Clas G4ParticleHP2N2AInelasticFS
418 InsertModel( 25041, "model_G4ParticleHP2N2
419 // Class G4ParticleHPNDInelasticFS
420 InsertModel( 25042, "model_G4ParticleHPNDI
421 // Class G4ParticleHPNTInelasticFS
422 InsertModel( 25043, "model_G4ParticleHPNTI
423 // Class G4ParticleHPNHe3InelasticFS
424 InsertModel( 25044, "model_G4ParticleHPNHe
425 // Class G4ParticleHPND2AInelasticFS
426 InsertModel( 25045, "model_G4ParticleHPND2
427 // Class G4ParticleHPNT2AInelasticFS
428 InsertModel( 25046, "model_G4ParticleHPNT2
429 // Class G4ParticleHP4NInelasticFS
430 InsertModel( 25047, "model_G4ParticleHP4NI
431 // Class G4ParticleHP2NPInelasticFS
432 InsertModel( 25048, "model_G4ParticleHP2NP
433 // Class G4ParticleHP3NPInelasticFS
434 InsertModel( 25049, "model_G4ParticleHP3NP
435 // Class G4ParticleHPN2PInelasticFS
436 InsertModel( 25050, "model_G4ParticleHPN2P
437 // Class G4ParticleHPNPAInelasticFS
438 InsertModel( 25051, "model_G4ParticleHPNPA
439 // Class G4ParticleHPPInelasticFS
440 InsertModel( 25052, "model_G4ParticleHPPIn
441 // Class G4ParticleHPDInelasticFS
442 InsertModel( 25053, "model_G4ParticleHPDIn
443 // Class G4ParticleHPTInelasticFS
444 InsertModel( 25054, "model_G4ParticleHPTIn
445 // Class G4ParticleHPHe3InelasticFS
446 InsertModel( 25055, "model_G4ParticleHPHe3
447 // Class G4ParticleHPAInelasticFS
448 InsertModel( 25056, "model_G4ParticleHPAIn
449 // Class G4ParticleHP2AInelasticFS
450 InsertModel( 25057, "model_G4ParticleHP2AI
451 // Class G4ParticleHP3AInelasticFS
452 InsertModel( 25058, "model_G4ParticleHP3AI
453 // Class G4ParticleHP2PInelasticFS
454 InsertModel( 25059, "model_G4ParticleHP2PI
455 // Class G4ParticleHPPAInelasticFS
456 InsertModel( 25060, "model_G4ParticleHPPAI
457 // Class G4ParticleHPD2AInelasticFS
458 InsertModel( 25061, "model_G4ParticleHPD2A
459 // Class G4ParticleHPT2AInelasticFS
460 InsertModel( 25062, "model_G4ParticleHPT2A
461 // Class G4ParticleHPPDInelasticFS
462 InsertModel( 25063, "model_G4ParticleHPPDI
463 // Class G4ParticleHPPTInelasticFS
464 InsertModel( 25064, "model_G4ParticleHPPTI
465 // Class G4ParticleHPDAInelasticFS
466 InsertModel( 25065, "model_G4ParticleHPDAI
467
468 // --- LEND : 25'200 - 25'299 ---
469
470 // Class: G4LENDModel
471 InsertModel( 25200, "model_LENDModel" );
472
473 // The following classes inherit from G4LE
474
475 // Class: G4LENDCapture
476 InsertModel( 25210, "model_LENDCapture" );
477 // Class: G4LENDElastic
478 InsertModel( 25220, "model_LENDElastic" );
479 // Class: G4LENDFission
480 InsertModel( 25230, "model_LENDFission" );
481 // Class: G4LENDInelastic
482 InsertModel( 25240, "model_LENDInelastic"
483 // Class: G4LENDCombinedModel
484 InsertModel( 25250, "model_LENDCombinedMod
485 // Class: G4LENDGammaModel
486 InsertModel( 25260, "model_LENDGammaModel"
487
488 // --- Radioactive Decay : 25'500 - 25'999
489
490 // 25'510 + 10*G4RadioactiveDecayMode
491 InsertModel( 25510, "model_RDM_IT" );
492 InsertModel( 25520, "model_RDM_BetaMinus"
493 InsertModel( 25530, "model_RDM_BetaPlus" )
494 InsertModel( 25540, "model_RDM_KshellEC" )
495 InsertModel( 25550, "model_RDM_LshellEC" )
496 InsertModel( 25560, "model_RDM_MshellEC" )
497 InsertModel( 25570, "model_RDM_NshellEC" )
498 InsertModel( 25580, "model_RDM_Alpha" );
499 InsertModel( 25590, "model_RDM_Proton" );
500 InsertModel( 25600, "model_RDM_Neutron" );
501 InsertModel( 25610, "model_RDM_SpFission"
502 InsertModel( 25620, "model_RDM_BDProton" )
503 InsertModel( 25630, "model_RDM_BDNeutron"
504 InsertModel( 25640, "model_RDM_Beta2Minus"
505 InsertModel( 25650, "model_RDM_Beta2Plus"
506 InsertModel( 25660, "model_RDM_Proton2" );
507 InsertModel( 25670, "model_RDM_Neutron2" )
508 InsertModel( 25680, "model_RDM_Triton" );
509
510 InsertModel( 25810, "model_RDM_AtomicRelax
511
512 // ---------------------------------------
513 // --- Others HAD (everything not include
514 // ---------------------------------------
515 // - 26'000 - 26'099 : Stopping
516 // - 26'100 - 26'199 : Fission
517 // - 26'200 - 26'299 : Abration
518 // - 26'300 - 26'399 : Coalescence
519 // ...
520
521 // --- Stopping : 26'000 - 26'099 ---
522
523 // Below are classes that derives from G4H
524 // The names are the same as before Geant4
525 // Classes that use it: G4HadronStoppingPr
526
527 // Class: G4HadronicAbsorptionBertini
528 InsertModel( 26000, "model_hBertiniCapture
529 InsertModel( 26001, "model_hBertiniCapture
530 InsertModel( 26002, "model_hBertiniCapture
531 // Class: G4HadronicAbsorptionFritiof
532 InsertModel( 26010, "model_hFritiofCapture
533 InsertModel( 26011, "model_hFritiofCapture
534 InsertModel( 26012, "model_hFritiofCapture
535 // Class: G4HadronicAbsorptionFritiofWithB
536 InsertModel( 26020, "model_hFritiofWithBin
537 InsertModel( 26021, "model_hFritiofWithBin
538 InsertModel( 26022, "model_hFritiofWithBin
539 // Class: G4MuonMinusCapture
540 InsertModel( 26030, "model_muMinusCaptureA
541 InsertModel( 26031, "model_muMinusCaptureA
542 InsertModel( 26032, "model_muMinusCaptureA
543
544 // --- Fission : 26'100 - 26'199 ---
545
546 // Class G4LFission
547 InsertModel( 26100, "model_G4LFission" );
548
549 // LLNL fission (related classes: G4Fissio
550 InsertModel( 26110, "model_G4LLNLFission"
551
552 // --- Abration : 26'200 - 26'299 ---
553
554 // Class G4WilsonAbrasionModel
555 InsertModel( 26200, "model_G4WilsonAbrasio
556
557 // --- Coalescence : 26'300 - 26'399 ---
558
559 // Class G4CRCoalescence
560 InsertModel( 26300, "model_G4CRCoalescence
561
562 // =======================================
563 // === 3rd OTHER (i.e. non-EM and non-HAD)
564 // =======================================
565
566 // -------------------------------
567 // --- Biasing : 30'000-30'999 ---
568 // -------------------------------
569
570 // The name is the same as before Geant4 1
571 // Classes that use it: G4BOptrForceCollis
572 InsertModel( 30010, "model_GenBiasForceCol
573
574 // ----------------------------------
575 // --- Channeling : 31'000-31'999 ---
576 // ----------------------------------
577
578 // The name is the same as before Geant4 1
579 // Classes that use it: G4Channeling , G4C
580 InsertModel( 31010, "model_channeling" );
581
582 // --- Others ... ---
583
584 // =======================================
585 // ================== 4th MODELS ADDED AFT
586 // =======================================
587 // PLEASE ADD MODELS ONLY BELOW HERE, WITH
588 // IF YOU ARE NOT SURE, PLEASE CONTACT ONE
589 // GEANT4 PHYSICS WORKING GROUPS.
590
591 // Class: G4HadronicAbsorptionINCLXX
592 InsertModel( 26040, "model_hINCLXXCaptureA
593 InsertModel( 26041, "model_hINCLXXCaptureA
594 InsertModel( 26042, "model_hINCLXXCaptureA
595
596 // Class: G4LightIonQMDReaction
597 InsertModel( 23310, "model_LightIonQMDMode
598
599 // Class: G4NuDEXEvaporation (to be added
600 // InsertModel( 24700, "model_nuDEX_evapor
601
602 // Class: G4NuDEXNeutronCapture
603 InsertModel( 25300, "model_nuDEX_neutronCa
604
605 // ...
606
607 SanityCheck();
608 isInitialized = true;
609
610 // The following call is commented out bec
611 // the verbosity level, but this would imp
612 // category to other categories which is n
613 // Anyhow, the call G4PhysicsModelCatalog
614 // can be easily placed elsewhere (in eith
615 // Geant4 classes).
616 //PrintAllInformation();
617 }
618 }
619
620 // -------------------------------------------
621 void G4PhysicsModelCatalog::SanityCheck() {
622 if ( theVectorOfModelIDs->size() != theVecto
623 G4ExceptionDescription ed;
624 ed << "theVectorOfModelIDs' size=" << theV
625 << " is NOT the same as theVectorOfMode
626 G4Exception( "G4PhysicsModelCatalog::Sanit
627 FatalException, ed, "should be the same!"
628 } else {
629 G4bool isModelIDOutsideRange = false;
630 G4bool isModelIDRepeated = false;
631 G4bool isModelNameRepeated = false;
632 for ( int idx = 0; idx < Entries(); ++idx
633 G4int modelID = (*theVectorOfModelIDs)[
634 G4String modelName = (*theVectorOfModelN
635 if ( modelID < GetMinAllowedModelIDValue
636 isModelIDOutsideRange = true;
637 }
638 for ( int jdx = idx + 1; jdx < Entries()
639 if(modelID == (*theVectorOfModelIDs)[j
640 {
641 isModelIDRepeated = true;
642 }
643 if(modelName == (*theVectorOfModelName
644 {
645 isModelNameRepeated = true;
646 }
647 }
648 }
649 if ( isModelIDOutsideRange || isModelIDRep
650 G4ExceptionDescription ed;
651 if(isModelIDOutsideRange)
652 {
653 ed << "theVectorOfModelIDs has NOT all
654 << GetMinAllowedModelIDValue() << "
655 << GetMaxAllowedModelIDValue();
656 }
657 if(isModelIDRepeated)
658 {
659 ed << "theVectorOfModelIDs has NOT all
660 }
661 if(isModelNameRepeated)
662 {
663 ed << "theVectorOfModelNames has NOT a
664 }
665 G4Exception( "G4PhysicsModelCatalog::San
666 FatalException, ed, "cannot continue!"
667 }
668 }
669 return;
670 }
671
672 // -------------------------------------------
673
674 // -------------------------------------------
675 const G4String G4PhysicsModelCatalog::GetModel
676 G4String modelName = "Undefined";
677 if ( modelID >= GetMinAllowedModelIDValue()
678 for ( int idx = 0; idx < Entries(); ++idx
679 if ( (*theVectorOfModelIDs)[ idx ] == mo
680 modelName = (*theVectorOfModelNames)[ idx ];
681 break;
682 }
683 }
684 }
685 return modelName;
686 }
687
688 // -------------------------------------------
689 const G4String G4PhysicsModelCatalog::GetModel
690 return ( modelIndex >= 0 && modelIndex < E
691 ? (*theVectorOfModelNames)[ modelIndex
692 }
693
694 // -------------------------------------------
695 G4int G4PhysicsModelCatalog::GetModelID( const
696 return ( modelIndex >= 0 && modelIndex < E
697 ? (*theVectorOfModelIDs)[ modelIndex ]
698 }
699
700 // -------------------------------------------
701 G4int G4PhysicsModelCatalog::GetModelID( const
702 if(!isInitialized)
703 {
704 Initialize();
705 }
706 G4int modelID = -1;
707 for ( G4int idx = 0; idx < Entries(); ++idx
708 if ( (*theVectorOfModelNames)[ idx ] == mo
709 modelID = (*theVectorOfModelIDs)[ idx ];
710 break;
711 }
712 }
713 return modelID;
714 }
715
716 // -------------------------------------------
717 G4int G4PhysicsModelCatalog::GetModelIndex( co
718 G4int modelIndex = -1;
719 if ( modelID >= GetMinAllowedModelIDValue()
720 for ( G4int idx = 0; idx < Entries(); ++id
721 if ( (*theVectorOfModelIDs)[ idx ] == mo
722 modelIndex = idx;
723 break;
724 }
725 }
726 }
727 return modelIndex;
728 }
729
730 // -------------------------------------------
731 G4int G4PhysicsModelCatalog::GetModelIndex( co
732 G4int modelIndex = -1;
733 for ( G4int idx = 0; idx < Entries(); ++idx
734 if ( (*theVectorOfModelNames)[ idx ] == mo
735 modelIndex = idx;
736 break;
737 }
738 }
739 return modelIndex;
740 }
741
742 // -------------------------------------------
743 G4int G4PhysicsModelCatalog::Entries() {
744 // It is enough to check the size of one of
745 return ( theVectorOfModelIDs != nullptr ) ?
746 }
747
748 // -------------------------------------------
749 void G4PhysicsModelCatalog::PrintAllInformatio
750 G4cout << G4endl
751 << " ======================================
752 << " === G4PhysicsModelCatalog::PrintAllInf
753 << " ======================================
754 << " SIZE (i.e. number of models in the cat
755 for ( int idx = 0; idx < Entries(); ++idx )
756 G4int modelID = (*theVectorOfModelIDs)[ id
757 G4String modelName = (*theVectorOfModelNam
758 G4cout << "\t index=" << idx << "\t modelN
759 << "\t modelID=" << modelID << G4endl;
760 }
761 G4cout << " ================================
762 << " ======================================
763 << " ======================================
764 << G4endl;
765 }
766