ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
PhysicsList.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file PhysicsList.cc
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 // -------------------------------------------------------------------
27 // -------------------------------------------------------------------
28 
29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
30 
31 #include "PhysicsList.hh"
32 #include "G4SystemOfUnits.hh"
33 
34 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
35 
36 PhysicsList::PhysicsList(): G4VUserPhysicsList()
37 {
38  defaultCutValue = 1*micrometer;
39  cutForGamma = defaultCutValue;
40  cutForElectron = defaultCutValue;
41  cutForPositron = defaultCutValue;
42  cutForProton = defaultCutValue;
43 
44  SetVerboseLevel(1);
45 }
46 
47 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
48 
49 PhysicsList::~PhysicsList()
50 {}
51 
52 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
53 
54 void PhysicsList::ConstructParticle()
55 {
56  ConstructBosons();
57  ConstructLeptons();
58  ConstructBarions();
59 }
60 
61 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
62 
63 void PhysicsList::ConstructBosons()
64 {
65  // gamma
66  G4Gamma::GammaDefinition();
67 }
68  //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
69 
70 void PhysicsList::ConstructLeptons()
71 {
72  // leptons
73  G4Electron::ElectronDefinition();
74  G4Positron::PositronDefinition();
75 }
76 
77 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
78 
79 void PhysicsList::ConstructBarions()
80 {
81  // baryons
82  G4Proton::ProtonDefinition();
83  G4GenericIon::GenericIonDefinition();
84 }
85 
86 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
87 
88 void PhysicsList::ConstructProcess()
89 {
90  AddTransportation();
91  ConstructEM();
92  ConstructGeneral();
93 }
94 
95 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
96 
97 // Geant4-MicroElec MODELS
98 
99 #include "G4MicroElecElastic.hh"
100 #include "G4MicroElecElasticModel.hh"
101 
102 #include "G4MicroElecInelastic.hh"
103 #include "G4MicroElecInelasticModel.hh"
104 //
105 
106 #include "G4LossTableManager.hh"
107 #include "G4EmConfigurator.hh"
108 #include "G4VEmModel.hh"
109 #include "G4DummyModel.hh"
110 #include "G4eIonisation.hh"
111 #include "G4hIonisation.hh"
112 #include "G4ionIonisation.hh"
113 #include "G4eMultipleScattering.hh"
114 #include "G4hMultipleScattering.hh"
115 #include "G4BraggModel.hh"
116 #include "G4BraggIonModel.hh"
117 #include "G4BetheBlochModel.hh"
118 #include "G4UrbanMscModel.hh"
119 #include "G4MollerBhabhaModel.hh"
120 #include "G4IonFluctuations.hh"
121 #include "G4UniversalFluctuation.hh"
122 
123 #include "G4ElectronCapture.hh"
124 
125 #include "G4UAtomicDeexcitation.hh"
126 
127 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
128 
129 void PhysicsList::ConstructEM()
130 {
131 
132  auto particleIterator=GetParticleIterator();
133  particleIterator->reset();
134 
135  while( (*particleIterator)() )
136  {
137 
138  G4ParticleDefinition* particle = particleIterator->value();
139  G4ProcessManager* pmanager = particle->GetProcessManager();
140  G4String particleName = particle->GetParticleName();
141 
142  // *********************************
143  // 1) Processes for the World region
144  // *********************************
145 
146  if (particleName == "e-") {
147 
148  // STANDARD msc is active in the world
149  G4eMultipleScattering* msc = new G4eMultipleScattering();
150  msc->AddEmModel(1, new G4UrbanMscModel());
151  pmanager->AddProcess(msc, -1, 1, -1);
152 
153  // STANDARD ionisation is active in the world
154  G4eIonisation* eion = new G4eIonisation();
155  pmanager->AddProcess(eion, -1, 2, 2);
156 
157  // MicroElec elastic is not active in the world
158  G4MicroElecElastic* theMicroElecElasticProcess = new G4MicroElecElastic("e-_G4MicroElecElastic");
159  theMicroElecElasticProcess->SetEmModel(new G4DummyModel(),1);
160  pmanager->AddDiscreteProcess(theMicroElecElasticProcess);
161 
162  // MicroElec ionisation is not active in the world
163  G4MicroElecInelastic* microelecioni = new G4MicroElecInelastic("e-_G4MicroElecInelastic");
164  microelecioni->SetEmModel(new G4DummyModel(),1);
165  pmanager->AddDiscreteProcess(microelecioni);
166 
167  // THE FOLLOWING PROCESS WILL KILL ALL ELECTRONS BELOW A SELECTED ENERY THRESHOLD
168  // Capture of low-energy e-
169  G4ElectronCapture* ecap = new G4ElectronCapture("Target", 16.7*eV);
170  pmanager->AddDiscreteProcess(ecap);
171 
172  } else if ( particleName == "proton" ) {
173 
174  // STANDARD msc is active in the world
175  G4hMultipleScattering* msc = new G4hMultipleScattering();
176  msc->AddEmModel(1, new G4UrbanMscModel());
177  pmanager->AddProcess(msc, -1, 1, -1);
178 
179  // STANDARD ionisation is active in the world
180  G4hIonisation* hion = new G4hIonisation();
181  pmanager->AddProcess(hion, -1, 2, 2);
182 
183  // MicroElec ionisation is not active in the world
184  G4MicroElecInelastic* microelecioni = new G4MicroElecInelastic("p_G4MicroElecInelastic");
185  microelecioni->SetEmModel(new G4DummyModel(),1);
186  microelecioni->SetEmModel(new G4DummyModel(),2);
187  pmanager->AddDiscreteProcess(microelecioni);
188 
189  } else if (particleName == "GenericIon") {
190 
191  // STANDARD msc is active in the world
192  G4hMultipleScattering* msc = new G4hMultipleScattering();
193  msc->AddEmModel(1, new G4UrbanMscModel());
194  pmanager->AddProcess(new G4hMultipleScattering, -1, 1, -1);
195 
196  // STANDARD ionisation is active in the world
197  G4ionIonisation* hion = new G4ionIonisation();
198  pmanager->AddProcess(hion, -1, 2, 2);
199 
200  // MicroElec ionisation is not active in the world
201  G4MicroElecInelastic* microelecioni = new G4MicroElecInelastic("ion_G4MicroElecInelastic");
202  microelecioni->SetEmModel(new G4DummyModel(),1);
203  microelecioni->SetEmModel(new G4DummyModel(),2);
204  pmanager->AddDiscreteProcess(microelecioni);
205  }
206  }
207 
208  // **************************************
209  // 2) Define processes for Target region
210  // **************************************
211 
212  // STANDARD EM processes should be inactivated when corresponding MicroElec processes are used
213  // - STANDARD EM e- processes are inactivated below 100 MeV
214  // - STANDARD EM proton & ion processes are inactivated below standEnergyLimit
215  //
216  G4EmConfigurator* em_config = G4LossTableManager::Instance()->EmConfigurator();
217 
218  G4VEmModel* mod;
219  // *** e-
220 
221  // ---> STANDARD EM processes are inactivated below 100 MeV
222 
223  G4UrbanMscModel* msc = new G4UrbanMscModel();
224  msc->SetActivationLowEnergyLimit(100*MeV);
225  em_config->SetExtraEmModel("e-","msc",msc,"Target");
226 
227  mod = new G4MollerBhabhaModel();
228  mod->SetActivationLowEnergyLimit(100*MeV);
229  em_config->SetExtraEmModel("e-","eIoni",mod,"Target",0.0,10*TeV, new G4UniversalFluctuation());
230 
231  // ---> MicroElec processes activated
232 
233  mod = new G4MicroElecElasticModel();
234  em_config->SetExtraEmModel("e-","e-_G4MicroElecElastic",mod,"Target",0.0,100*MeV);
235 
236  mod = new G4MicroElecInelasticModel();
237  em_config->SetExtraEmModel("e-","e-_G4MicroElecInelastic",mod,"Target",16.7*eV,100*MeV);
238 
239  // *** proton
240 
241  // ---> STANDARD EM processes inactivated below standEnergyLimit
242 
243  // STANDARD msc is still active
244  // Inactivate following STANDARD processes
245 
246  mod = new G4BraggModel();
247  mod->SetActivationHighEnergyLimit(50*keV);
248  em_config->SetExtraEmModel("proton","hIoni",mod,"Target",0.0,2*MeV, new G4IonFluctuations());
249 
250  mod = new G4BetheBlochModel();
251  mod->SetActivationLowEnergyLimit(10*GeV);
252  em_config->SetExtraEmModel("proton","hIoni",mod,"Target",2*MeV,10*TeV, new G4IonFluctuations());
253 
254  // ---> MicroElec processes activated
255  mod = new G4MicroElecInelasticModel();
256  mod->SetActivationLowEnergyLimit(50*keV);
257  em_config->SetExtraEmModel("proton","p_G4MicroElecInelastic",mod,"Target",0.0,10*GeV);
258  em_config->SetExtraEmModel("proton","p_G4MicroElecInelastic",new G4DummyModel,
259  "Target",10*GeV,10*TeV);
260 
261  // *** ion
262 
263  // ---> STANDARD EM processes inactivated below standEnergyLimit
264 
265  // STANDARD msc is still active
266  // Inactivate following STANDARD processes
267 
268  mod = new G4BraggIonModel();
269  mod->SetActivationHighEnergyLimit(50*keV);
270  em_config->SetExtraEmModel("GenericIon","ionIoni",mod,"Target",0.0,2*MeV, new G4IonFluctuations());
271 
272  mod = new G4BetheBlochModel();
273  mod->SetActivationLowEnergyLimit(10*GeV);
274  em_config->SetExtraEmModel("GenericIon","ionIoni",mod,"Target",2*MeV,10*TeV, new G4IonFluctuations());
275 
276  // ---> MicroElec processes activated
277  mod = new G4MicroElecInelasticModel();
278  mod->SetActivationLowEnergyLimit(50*keV);
279  em_config->SetExtraEmModel("GenericIon","ion_G4MicroElecInelastic",mod,
280  "Target",0.0,10*GeV);
281  em_config->SetExtraEmModel("GenericIon","ion_G4MicroElecInelastic",new G4DummyModel,
282  "Target",10*GeV,10*TeV);
283 
284  // Deexcitation
285  //
286  G4VAtomDeexcitation* de = new G4UAtomicDeexcitation();
287  G4LossTableManager::Instance()->SetAtomDeexcitation(de);
288  de->SetFluo(true);
289  de->SetAuger(true);
290  de->SetPIXE(true);
291  de->InitialiseForNewRun();
292 
293 
294 }
295 
296 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
297 
298 void PhysicsList::ConstructGeneral()
299 { }
300 
301 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
302 
303 void PhysicsList::SetCuts()
304 {
305  if (verboseLevel >0)
306  {
307  G4cout << "PhysicsList::SetCuts:";
308  G4cout << "CutLength : " << G4BestUnit(defaultCutValue,"Length") << G4endl;
309  }
310 
311  // set cut values for gamma at first and for e- second and next for e+,
312  // because some processes for e+/e- need cut values for gamma
313  SetCutValue(cutForGamma, "gamma");
314  SetCutValue(cutForElectron, "e-");
315  SetCutValue(cutForPositron, "e+");
316  SetCutValue(cutForProton, "proton");
317 
318  if (verboseLevel>0) { DumpCutValuesTable(); }
319 }