Geant4_10
G4PhysListFactory.cc
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 // $Id: G4PhysListFactory.cc 68750 2013-04-05 10:19:04Z gcosmo $
27 //
28 //---------------------------------------------------------------------------
29 //
30 // ClassName: G4PhysListFactory
31 //
32 // Author: 21 April 2008 V. Ivanchenko
33 //
34 // Modified:
35 //
36 //----------------------------------------------------------------------------
37 //
38 
39 #include "G4PhysListFactory.hh"
40 #include "FTFP_BERT.hh"
41 #include "FTFP_BERT_HP.hh"
42 #include "FTFP_BERT_TRV.hh"
43 #include "FTFP_INCLXX.hh"
44 #include "FTFP_INCLXX_HP.hh"
45 #include "FTF_BIC.hh"
46 #include "LBE.hh"
47 #include "QBBC.hh"
48 #include "QGSP_BERT.hh"
49 #include "QGSP_BERT_HP.hh"
50 #include "QGSP_BIC.hh"
51 #include "QGSP_BIC_HP.hh"
52 #include "QGSP_FTFP_BERT.hh"
53 #include "QGS_BIC.hh"
54 #include "QGSP_INCLXX.hh"
55 #include "QGSP_INCLXX_HP.hh"
56 #include "Shielding.hh"
57 
58 #include "G4EmStandardPhysics.hh"
63 #include "G4EmLivermorePhysics.hh"
64 #include "G4EmPenelopePhysics.hh"
65 
67  : defName("FTFP_BERT"),verbose(1)
68 {
69  nlists_hadr = 18;
70  G4String ss[18] = {
71  "FTFP_BERT","FTFP_BERT_TRV","FTFP_BERT_HP","FTFP_INCLXX",
72  "FTFP_INCLXX_HP","FTF_BIC", "LBE","QBBC",
73  "QGSP_BERT","QGSP_BERT_HP","QGSP_BIC","QGSP_BIC_HP",
74  "QGSP_FTFP_BERT","QGSP_INCLXX","QGSP_INCLXX_HP","QGS_BIC",
75  "Shielding","ShieldingLEND"};
76  for(size_t i=0; i<nlists_hadr; ++i) {
77  listnames_hadr.push_back(ss[i]);
78  }
79 
80  nlists_em = 7;
81  G4String s1[7] = {"","_EMV","_EMX","_EMY","_EMZ","_LIV","_PEN"};
82  for(size_t i=0; i<nlists_em; ++i) {
83  listnames_em.push_back(s1[i]);
84  }
85 }
86 
88 {}
89 
92 {
93  // instantiate PhysList by environment variable "PHYSLIST"
94  G4String name = "";
95  char* path = getenv("PHYSLIST");
96  if (path) {
97  name = G4String(path);
98  } else {
99  name = defName;
100  G4cout << "### G4PhysListFactory WARNING: "
101  << " environment variable PHYSLIST is not defined"
102  << G4endl
103  << " Default Physics Lists " << name
104  << " is instantiated"
105  << G4endl;
106  }
107  return GetReferencePhysList(name);
108 }
109 
112 {
113  // analysis on the string
114  size_t n = name.size();
115 
116  // last characters in the string
117  size_t em_opt = 0;
118  G4String em_name = "";
119 
120  // check EM options
121  if(n > 4) {
122  em_name = name.substr(n - 4, 4);
123  for(size_t i=1; i<nlists_em; ++i) {
124  if(listnames_em[i] == em_name) {
125  em_opt = i;
126  n -= 4;
127  break;
128  }
129  }
130  if(0 == em_opt) { em_name = ""; }
131  }
132 
133  // hadronic pHysics List
134  G4String had_name = name.substr(0, n);
135 
136  if(0 < verbose) {
137  G4cout << "G4PhysListFactory::GetReferencePhysList <" << had_name
138  << em_name << "> EMoption= " << em_opt << G4endl;
139  }
141  if(had_name == "FTFP_BERT") {p = new FTFP_BERT(verbose);}
142  else if(had_name == "FTFP_BERT_HP") {p = new FTFP_BERT_HP(verbose);}
143  else if(had_name == "FTFP_BERT_TRV") {p = new FTFP_BERT_TRV(verbose);}
144  else if(had_name == "FTFP_INCLXX") {p = new FTFP_INCLXX(verbose);}
145  else if(had_name == "FTFP_INCLXX_HP") {p = new FTFP_INCLXX_HP(verbose);}
146  else if(had_name == "FTF_BIC") {p = new FTF_BIC(verbose);}
147  else if(had_name == "LBE") {p = new LBE();}
148  else if(had_name == "QBBC") {p = new QBBC(verbose);}
149  else if(had_name == "QGSP_BERT") {p = new QGSP_BERT(verbose);}
150  else if(had_name == "QGSP_BERT_HP") {p = new QGSP_BERT_HP(verbose);}
151  else if(had_name == "QGSP_BIC") {p = new QGSP_BIC(verbose);}
152  else if(had_name == "QGSP_BIC_HP") {p = new QGSP_BIC_HP(verbose);}
153  else if(had_name == "QGSP_FTFP_BERT") {p = new QGSP_FTFP_BERT(verbose);}
154  else if(had_name == "QGSP_INCLXX") {p = new QGSP_INCLXX(verbose);}
155  else if(had_name == "QGSP_INCLXX_HP") {p = new QGSP_INCLXX_HP(verbose);}
156  else if(had_name == "QGS_BIC") {p = new QGS_BIC(verbose);}
157  else if(had_name == "Shielding") {p = new Shielding(verbose);}
158  else if(had_name == "ShieldingLEND") {p = new Shielding(verbose,"LEND");}
159  else {
160  G4cout << "### G4PhysListFactory WARNING: "
161  << "PhysicsList " << had_name << " is not known"
162  << G4endl;
163  }
164  if(p) {
165  G4cout << "<<< Reference Physics List " << had_name
166  << em_name << " is built" << G4endl;
167  G4int ver = p->GetVerboseLevel();
168  p->SetVerboseLevel(0);
169  if(0 < em_opt) {
170  if(1 == em_opt) {
172  } else if(2 == em_opt) {
174  } else if(3 == em_opt) {
176  } else if(4 == em_opt) {
178  } else if(5 == em_opt) {
179  p->ReplacePhysics(new G4EmLivermorePhysics(verbose));
180  } else if(6 == em_opt) {
181  p->ReplacePhysics(new G4EmPenelopePhysics(verbose));
182  }
183  }
184  p->SetVerboseLevel(ver);
185  }
186  G4cout << G4endl;
187  return p;
188 }
189 
191 {
192  G4bool res = false;
193  size_t n = name.size();
194  if(n > 4) {
195  G4String em_name = name.substr(n - 4, 4);
196  for(size_t i=1; i<nlists_em; ++i) {
197  if(listnames_em[i] == em_name) {
198  n -= 4;
199  break;
200  }
201  }
202  }
203  G4String had_name = name.substr(0, n);
204  for(size_t i=0; i<nlists_hadr; ++i) {
205  if(had_name == listnames_hadr[i]) {
206  res = true;
207  break;
208  }
209  }
210  return res;
211 }
212 
213 const std::vector<G4String>&
215 {
216  return listnames_hadr;
217 }
218 
219 const std::vector<G4String>&
221 {
222  return listnames_em;
223 }
224 
TQGSP_FTFP_BERT< G4VModularPhysicsList > QGSP_FTFP_BERT
TQGSP_BERT< G4VModularPhysicsList > QGSP_BERT
Definition: QGSP_BERT.hh:63
TQGSP_BIC_HP< G4VModularPhysicsList > QGSP_BIC_HP
Definition: QGSP_BIC_HP.hh:63
TFTFP_BERT_HP< G4VModularPhysicsList > FTFP_BERT_HP
Definition: FTFP_BERT_HP.hh:63
TShielding< G4VModularPhysicsList > Shielding
Definition: Shielding.hh:65
const char * p
Definition: xmltok.h:285
TLBE< G4VModularPhysicsList > LBE
Definition: LBE.hh:123
const XML_Char * name
Definition: expat.h:151
int G4int
Definition: G4Types.hh:78
Char_t n[5]
G4GLOB_DLL std::ostream G4cout
G4VModularPhysicsList * GetReferencePhysList(const G4String &)
G4bool IsReferencePhysList(const G4String &)
bool G4bool
Definition: G4Types.hh:79
TQGS_BIC< G4VModularPhysicsList > QGS_BIC
Definition: QGS_BIC.hh:63
const std::vector< G4String > & AvailablePhysListsEM() const
TFTFP_BERT_TRV< G4VModularPhysicsList > FTFP_BERT_TRV
void SetVerboseLevel(G4int value)
void ReplacePhysics(G4VPhysicsConstructor *)
TFTFP_BERT< G4VModularPhysicsList > FTFP_BERT
Definition: FTFP_BERT.hh:63
TINCLXXPhysicsListHelper< G4VModularPhysicsList, false, true > FTFP_INCLXX
Definition: FTFP_INCLXX.hh:41
TINCLXXPhysicsListHelper< G4VModularPhysicsList, false, false > QGSP_INCLXX
TINCLXXPhysicsListHelper< G4VModularPhysicsList, true, false > QGSP_INCLXX_HP
#define G4endl
Definition: G4ios.hh:61
G4VModularPhysicsList * ReferencePhysList()
TQGSP_BIC< G4VModularPhysicsList > QGSP_BIC
Definition: QGSP_BIC.hh:62
TFTF_BIC< G4VModularPhysicsList > FTF_BIC
Definition: FTF_BIC.hh:62
TQGSP_BERT_HP< G4VModularPhysicsList > QGSP_BERT_HP
Definition: QGSP_BERT_HP.hh:62
Definition: QBBC.hh:44
TINCLXXPhysicsListHelper< G4VModularPhysicsList, true, true > FTFP_INCLXX_HP
const std::vector< G4String > & AvailablePhysLists() const