| 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 // Authors: Francesco Longo, franzlongo1969@gm <<
27 // 26 //
28 // Code based on the hadrontherapy && radiopro << 27 // $Id: GammaRayTelPhysicsList.cc,v 1.7 2006/06/29 15:56:55 gunter Exp $
>> 28 // GEANT4 tag $Name: geant4-08-03-patch-01 $
>> 29 //
>> 30 //
29 31
30 #include "GammaRayTelPhysicsList.hh" 32 #include "GammaRayTelPhysicsList.hh"
31 #include "GammaRayTelPhysicsListMessenger.hh" <<
32 #include "G4PhysListFactory.hh" <<
33 #include "G4VPhysicsConstructor.hh" <<
34 <<
35 // Physics lists (contained inside the Geant4 <<
36 #include "G4EmLivermorePhysics.hh" <<
37 #include "G4EmLivermorePolarizedPhysics.hh" <<
38 #include "G4EmPenelopePhysics.hh" <<
39 #include "G4EmStandardPhysics_option3.hh" <<
40 #include "G4EmStandardPhysics_option4.hh" // t <<
41 <<
42 #include "G4Decay.hh" <<
43 #include "G4DecayPhysics.hh" <<
44 #include "G4HadronDElasticPhysics.hh" <<
45 #include "G4HadronElasticPhysics.hh" <<
46 #include "G4HadronElasticPhysicsHP.hh" <<
47 #include "G4HadronPhysicsQGSP_BIC_HP.hh" <<
48 #include "G4IonBinaryCascadePhysics.hh" <<
49 #include "G4IonFluctuations.hh" <<
50 #include "G4IonParametrisedLossModel.hh" <<
51 #include "G4LossTableManager.hh" <<
52 #include "G4ProcessManager.hh" <<
53 #include "G4RadioactiveDecayPhysics.hh" <<
54 #include "G4StepLimiter.hh" <<
55 #include "G4SystemOfUnits.hh" <<
56 #include "G4UnitsTable.hh" <<
57 <<
58 //....oooOO0OOooo........oooOO0OOooo........oo <<
59 <<
60 GammaRayTelPhysicsList::GammaRayTelPhysicsList <<
61 G4LossTableManager::Instance(); <<
62 <<
63 constexpr auto DEFAULT_CUT_VALUE{100 * mic <<
64 defaultCutValue = DEFAULT_CUT_VALUE; <<
65 <<
66 SetVerboseLevel(1); <<
67 <<
68 constexpr auto ENERGY_LOWER_BOUND{250 * eV <<
69 constexpr auto ENERGY_UPPER_BOUND{1 * GeV} <<
70 <<
71 G4ProductionCutsTable::GetProductionCutsTa <<
72 SetDefaultCutValue (defaultCutValue); <<
73 DumpCutValuesTable(); <<
74 <<
75 helIsRegisted = false; <<
76 bicIsRegisted = false; <<
77 biciIsRegisted = false; <<
78 locIonIonInelasticIsRegistered = false; <<
79 radioactiveDecayIsRegisted = false; <<
80 <<
81 pMessenger = new GammaRayTelPhysicsListMes <<
82 <<
83 SetVerboseLevel(1); <<
84 <<
85 // EM physics <<
86 emPhysicsList = new G4EmStandardPhysics_op <<
87 emName = G4String("emstandard_opt3"); <<
88 <<
89 // Decay physics and all particles <<
90 decPhysicsList = new G4DecayPhysics(); <<
91 } <<
92 <<
93 //....oooOO0OOooo........oooOO0OOooo........oo <<
94 <<
95 GammaRayTelPhysicsList::~GammaRayTelPhysicsLis <<
96 delete pMessenger; <<
97 delete emPhysicsList; <<
98 delete decPhysicsList; <<
99 <<
100 for (auto &hadronPhy : hadronPhys) { <<
101 delete hadronPhy; <<
102 } <<
103 } <<
104 33
105 //....oooOO0OOooo........oooOO0OOooo........oo << 34 #include "globals.hh"
106 << 35 #include "G4ParticleDefinition.hh"
107 void GammaRayTelPhysicsList::AddPackage(const << 36 #include "G4ParticleWithCuts.hh"
108 G4PhysListFactory factory; << 37 #include "G4ProcessManager.hh"
109 auto *phys = factory.GetReferencePhysList( << 38 #include "G4ProcessVector.hh"
110 << 39 #include "G4ParticleTypes.hh"
111 G4int i = 0; << 40 #include "G4ParticleTable.hh"
112 const G4VPhysicsConstructor *element = phy << 41
113 auto *tmp = const_cast<G4VPhysicsConstruct << 42 #include "G4Material.hh"
114 << 43 #include "G4MaterialTable.hh"
115 while (element != nullptr) { << 44 #include "G4ios.hh"
116 RegisterPhysics(tmp); << 45 #include <iomanip>
117 element = phys->GetPhysics(++i); << 46 #include "GammaRayTelParticles.hh"
118 tmp = const_cast<G4VPhysicsConstructor << 47 #include "GammaRayTelGeneralPhysics.hh"
119 } << 48 #include "GammaRayTelEMPhysics.hh"
120 } << 49 #include "GammaRayTelMuonPhysics.hh"
121 << 50 #include "GammaRayTelHadronPhysics.hh"
122 //....oooOO0OOooo........oooOO0OOooo........oo << 51 #include "GammaRayTelIonPhysics.hh"
123 << 52
124 void GammaRayTelPhysicsList::ConstructParticle << 53 GammaRayTelPhysicsList::GammaRayTelPhysicsList(): G4VModularPhysicsList()
125 decPhysicsList->ConstructParticle(); << 54 {
126 } << 55 // default cut value (1.0mm)
127 << 56 defaultCutValue = 1.0*mm;
128 //....oooOO0OOooo........oooOO0OOooo........oo << 57 SetVerboseLevel(1);
129 << 58
130 void GammaRayTelPhysicsList::ConstructProcess( << 59 RegisterPhysics( new GammaRayTelParticles("particles") );
131 // transportation << 60
132 // << 61 // General Physics
133 AddTransportation(); << 62 RegisterPhysics( new GammaRayTelGeneralPhysics("general") );
134 << 63
135 // electromagnetic physics list << 64 // EM Physics
136 // << 65 RegisterPhysics( new GammaRayTelEMPhysics("electromagnetic"));
137 emPhysicsList->ConstructProcess(); << 66
138 emConfigurator.AddModels(); << 67 // Muon Physics
139 << 68 RegisterPhysics( new GammaRayTelMuonPhysics("muon"));
140 // decay physics list << 69
141 // << 70 // Hadron Physics
142 decPhysicsList->ConstructProcess(); << 71 RegisterPhysics( new GammaRayTelHadronPhysics("hadron"));
143 << 72
144 // hadronic physics lists << 73 // Ion Physics
145 for (auto &hadronPhy : hadronPhys) { << 74 RegisterPhysics( new GammaRayTelIonPhysics("ion"));
146 hadronPhy->ConstructProcess(); << 75
147 } << 76
148 << 77 }
149 // step limitation (as a full process) << 78
150 // AddStepMax(); << 79 GammaRayTelPhysicsList::~GammaRayTelPhysicsList()
151 } << 80 {
152 << 81 }
153 //....oooOO0OOooo........oooOO0OOooo........oo << 82
154 << 83 #include "G4Region.hh"
155 void GammaRayTelPhysicsList::AddPhysicsList(co << 84 #include "G4RegionStore.hh"
156 if (verboseLevel > 1) { << 85 #include "G4ProductionCuts.hh"
157 G4cout << "PhysicsList::AddPhysicsList << 86
158 } << 87 void GammaRayTelPhysicsList::SetCuts()
159 if (name == emName) { << 88 {
160 return; << 89 if (verboseLevel >0){
161 } << 90 G4cout << "GammaRayTelPhysicsList::SetCuts: default cut length : "
162 << 91 << G4BestUnit(defaultCutValue,"Length") << G4endl;
163 //////////////////////////////// << 92 }
164 // ELECTROMAGNETIC MODELS // << 93
165 //////////////////////////////// << 94 // These values are used as the default production thresholds
166 << 95 // for the world volume.
167 if (name == "standard_opt3") { << 96 SetCutsWithDefault();
168 emName = name; << 97
169 delete emPhysicsList; << 98
170 emPhysicsList = new G4EmStandardPhysic << 99 // Production thresholds for detector regions
171 G4cout << "THE FOLLOWING ELECTROMAGNET << 100
172 } else if (name == "LowE_Livermore") { << 101 G4String regName[] = {"Calorimeter","Tracker"};
173 emName = name; << 102 G4double cutValue[] = {1*mm, 0.1*mm};
174 delete emPhysicsList; << 103
175 emPhysicsList = new G4EmLivermorePhysi << 104 for(G4int i=0;i<2;i++)
176 G4cout << "THE FOLLOWING ELECTROMAGNET << 105 {
177 } else if (name == "LowE_Penelope") { << 106 G4Region* reg = G4RegionStore::GetInstance()->GetRegion(regName[i]);
178 emName = name; << 107 G4ProductionCuts* cuts = new G4ProductionCuts;
179 delete emPhysicsList; << 108 cuts->SetProductionCut(cutValue[i]);
180 emPhysicsList = new G4EmPenelopePhysic << 109 reg->SetProductionCuts(cuts);
181 G4cout << "THE FOLLOWING ELECTROMAGNET << 110 }
182 } else if (name == "LowE_Polarized") { <<
183 emName = name; <<
184 delete emPhysicsList; <<
185 emPhysicsList = new G4EmLivermorePolar <<
186 G4cout << "THE FOLLOWING ELECTROMAGNET <<
187 } else if (name == "standard_opt4") { <<
188 emName = name; <<
189 delete emPhysicsList; <<
190 emPhysicsList = new G4EmStandardPhysic <<
191 G4cout << "THE FOLLOWING ELECTROMAGNET <<
192 <<
193 ///////////////////////// <<
194 // HADRONIC MODELS // <<
195 ///////////////////////// <<
196 <<
197 } else if (name == "elastic" && !helIsRegi <<
198 G4cout << "THE FOLLOWING HADRONIC ELAS <<
199 hadronPhys.push_back(new G4HadronElast <<
200 helIsRegisted = true; <<
201 } else if (name == "DElastic" && !helIsReg <<
202 hadronPhys.push_back(new G4HadronDElas <<
203 helIsRegisted = true; <<
204 } else if (name == "HPElastic" && !helIsRe <<
205 hadronPhys.push_back(new G4HadronElast <<
206 helIsRegisted = true; <<
207 } else if (name == "binary" && !bicIsRegis <<
208 hadronPhys.push_back(new G4HadronPhysi <<
209 bicIsRegisted = true; <<
210 G4cout << "THE FOLLOWING HADRONIC INEL <<
211 } else if (name == "binary_ion" && !biciIs <<
212 hadronPhys.push_back(new G4IonBinaryCa <<
213 biciIsRegisted = true; <<
214 G4cout << "THE FOLLOWING HADRONIC INEL <<
215 } else if (name == "radioactive_decay" && <<
216 hadronPhys.push_back(new G4Radioactive <<
217 radioactiveDecayIsRegisted = true; <<
218 G4cout << "THE FOLLOWING HADRONIC INEL <<
219 } else { <<
220 G4cout << "PhysicsList::AddPhysicsList <<
221 } <<
222 } <<
223 <<
224 //....oooOO0OOooo........oooOO0OOooo........oo <<
225 <<
226 void GammaRayTelPhysicsList::SetCutForGamma(G4 <<
227 SetParticleCuts(cut, G4Gamma::Gamma()); <<
228 } 111 }
229 112
230 //....oooOO0OOooo........oooOO0OOooo........oo <<
231 113
232 void GammaRayTelPhysicsList::SetCutForElectron <<
233 SetParticleCuts(cut, G4Electron::Electron( <<
234 } <<
235 114
236 //....oooOO0OOooo........oooOO0OOooo........oo <<
237 115
238 void GammaRayTelPhysicsList::SetCutForPositron <<
239 SetParticleCuts(cut, G4Positron::Positron( <<
240 } <<
241 116