Geant4  10.03.p02
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
PurgMagPhysicsList.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 // Code developed by:
27 // S.Larsson
28 //
29 // ********************************
30 // * *
31 // * PurgMagPhysicsList.cc *
32 // * *
33 // ********************************
34 //
35 //
36 // $Id: PurgMagPhysicsList.cc 100654 2016-10-31 10:14:05Z gcosmo $
37 //
38 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
39 
40 #include "PurgMagPhysicsList.hh"
41 
42 #include "G4SystemOfUnits.hh"
43 #include "G4ParticleDefinition.hh"
44 #include "G4ParticleWithCuts.hh"
45 #include "G4ProcessManager.hh"
46 #include "G4ParticleTypes.hh"
47 #include "G4ParticleTable.hh"
48 #include "G4Material.hh"
49 #include "G4UnitsTable.hh"
50 #include "G4ios.hh"
51 
52 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
53 
55 {
57  cutForGamma = defaultCutValue;
58  cutForElectron = defaultCutValue;
59  cutForPositron = defaultCutValue;
60  cutForProton = defaultCutValue;
61 
62  SetVerboseLevel(1);
63 }
64 
65 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
66 
68 {}
69 
70 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
71 
73 {
74  // In this method, static member functions should be called
75  // for all particles which you want to use.
76  // This ensures that objects of these particle types will be
77  // created in the program.
78 
82 }
83 
84 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
85 
87 {
88 
89  // gamma
91 
92  // optical photon
94 }
95  //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
96 
98 {
99  // leptons
102 }
103 
104 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
105 
107 {
108  // barions
111 }
112 
113 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
114 
116 {
118  ConstructEM();
120 }
121 
122 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
123 
124 #include "G4eMultipleScattering.hh"
125 #include "G4hMultipleScattering.hh"
126 // Bosons
127 #include "G4PhotoElectricEffect.hh"
128 #include "G4ComptonScattering.hh"
129 #include "G4GammaConversion.hh"
130 // Leptons
131 #include "G4eIonisation.hh"
132 #include "G4eBremsstrahlung.hh"
133 #include "G4eplusAnnihilation.hh"
134 
135 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
136 
138 {
140  particleIterator->reset();
141  while( (*particleIterator)() )
142  {
143  G4ParticleDefinition* particle = particleIterator->value();
144  G4ProcessManager* pmanager = particle->GetProcessManager();
145  G4String particleName = particle->GetParticleName();
146 
147 
148  if (particleName == "gamma") {
149  //gamma
150  pmanager->AddDiscreteProcess(new G4GammaConversion);
153 
154  } else if (particleName == "e-") {
155  //electron
156  pmanager->AddProcess(new G4eBremsstrahlung, -1,-1,3);
157  pmanager->AddProcess(new G4eIonisation, -1, 2,2);
158  pmanager->AddProcess(new G4eMultipleScattering, -1, 1,1);
159 
160  } else if (particleName == "e+") {
161  //positron
162  pmanager->AddProcess(new G4eBremsstrahlung, -1,-1,3);
163  pmanager->AddProcess(new G4eIonisation, -1, 2,2);
164  pmanager->AddProcess(new G4eMultipleScattering, -1, 1,1);
165  pmanager->AddProcess(new G4eplusAnnihilation, 0,-1,4);
166 
167  } else if (particleName == "proton") {
168  //proton
169  pmanager->AddProcess(new G4hMultipleScattering, -1,1,1);
170 
171  } else if (particleName == "anti_proton") {
172  //antiproton
173  pmanager->AddProcess(new G4hMultipleScattering, -1,1,1);
174  }
175  }
176 }
177 
178 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
179 
181 { }
182 
183 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
184 
186 {
187  if (verboseLevel >0){
188  G4cout << "PurgMagPhysicsList::SetCuts:";
189  G4cout << "CutLength : " << G4BestUnit(defaultCutValue,"Length") << G4endl;
190  }
191 
192  // set cut values for gamma at first and for e- second and next for e+,
193  // because some processes for e+/e- need cut values for gamma
194  SetCutValue(cutForGamma, "gamma");
195  SetCutValue(cutForElectron, "e-");
196  SetCutValue(cutForPositron, "e+");
197 
198  // set cut values for proton and anti_proton before all other hadrons
199  // because some processes for hadrons need cut values for proton/anti_proton
200  SetCutValue(cutForProton, "proton");
201  SetCutValue(cutForProton, "anti_proton");
202 
203  // SetCutValueForOthers(defaultCutValue);
204 
206 }
207 
208 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
209 
211 {
212  if (verboseLevel >0){
213  G4cout << "PurgMagPhysicsList::SetCuts:";
214  G4cout << "Gamma cut in energy: " << lowcut*MeV << " (MeV)" << G4endl;
215  }
216 
217  // G4Gamma::SetEnergyRange(lowcut,1e5);
218  SetGELowLimit(lowcut);
219 }
220 
221 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
222 
224 {
225  if (verboseLevel >0){
226 
227  G4cout << "PurgMagPhysicsList::SetCuts:";
228  G4cout << "Electron cut in energy: " << lowcut*MeV << " (MeV)" << G4endl;
229  }
230 
231  // G4Electron::SetEnergyRange(lowcut,1e5);
232  SetGELowLimit(lowcut);
233 }
234 
235 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
236 
238 {
239  if (verboseLevel >0){
240 
241  G4cout << "PurgMagPhysicsList::SetCuts:";
242  G4cout << "Positron cut in energy: " << lowcut*MeV << " (MeV)" << G4endl;
243  }
244 
245  G4cerr << "PurgMagPhysicsList::SetPositronLowLimit: Not currently able to set Positron LowLimit." << G4endl;
246  G4Exception("PurgMagPhysicsList::SetPositronLowLimit()","PurMag001",
247  FatalException,"Positron Low Limit: not implemented in PurgMagPhysicsList");
248  //
249  // G4Positron::SetEnergyRange(lowcut,1e5);
250 }
251 
252 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
253 
255 {
256  if (verboseLevel >0){
257 
258  G4cout << "PurgMagPhysicsList::SetCuts:";
259  G4cout << "Proton cut in energy: " << lowcut*MeV << " (MeV)" << G4endl;
260  }
261 
262  G4cerr << "PurgMagPhysicsList::SetProtonLowLimit: Not currently able to set Proton LowLimit." << G4endl;
263  G4Exception("PurgMagPhysicsList::SetProtonLowLimit()","PurMag002",
264  FatalException,"Proton Low Limit: not implemented in PurgMagPhysicsList");
265  //
266  // G4Proton::SetEnergyRange(lowcut,1e5);
267  // G4AntiProton::SetEnergyRange(lowcut,1e5);
268 }
269 
270 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
271 
273 {
274  if (verboseLevel >0){
275  G4cout << "PurgMagPhysicsList::SetGEPLowLimit:";
276  G4cout << "Gamma and Electron cut in energy: " << lowcut*MeV << " (MeV)" << G4endl;
277  }
278 
279  // G4Gamma::SetEnergyRange(lowcut,1e5);
280  // G4Electron::SetEnergyRange(lowcut,1e5);
281  // G4Positron::SetEnergyRange(lowcut,1e5);
282  this->SetGELowLimit(lowcut);
283 
284  G4cerr << " SetGEPLowLimit : Uncertain whether setting Positron low limit " << G4endl;
285 }
286 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
287 
288 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
290 {
291  if (verboseLevel >0){
292  G4cout << "PurgMagPhysicsList::SetGELowLimit:";
293  G4cout << "Gamma and Electron cut in energy: " << lowcut*MeV << " (MeV)" << G4endl;
294  }
295 
297 }
299 {
300  ResetCuts();
301  cutForGamma = val;
302 }
303 
304 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
305 
307 {
308  // ResetCuts();
309  cutForElectron = val;
310 }
311 
312 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
313 
315 {
316  // ResetCuts();
317  cutForPositron = val;
318 }
319 
320 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
321 
323 {
324  //ResetCuts();
325  cutForProton = val;
326 }
327 
328 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
329 
330 
331 
332 
333 
334 
static G4Electron * ElectronDefinition()
Definition: G4Electron.cc:89
void SetGEPLowLimit(G4double)
void SetPositronCut(G4double)
void SetCutValue(G4double aCut, const G4String &pname)
void SetEnergyRange(G4double lowedge, G4double highedge)
static G4Proton * ProtonDefinition()
Definition: G4Proton.cc:88
G4int AddDiscreteProcess(G4VProcess *aProcess, G4int ord=ordDefault)
#define G4BestUnit(a, b)
#define G4_USE_G4BESTUNIT_FOR_VERBOSE 1
void SetElectronCut(G4double)
void SetGELowLimit(G4double)
void SetGammaCut(G4double)
static G4AntiProton * AntiProtonDefinition()
Definition: G4AntiProton.cc:88
const G4String & GetParticleName() const
void SetGammaLowLimit(G4double)
void DumpCutValuesTable(G4int flag=1)
G4GLOB_DLL std::ostream G4cout
void SetProtonCut(G4double)
void SetVerboseLevel(G4int value)
G4int AddProcess(G4VProcess *aProcess, G4int ordAtRestDoIt=ordInActive, G4int ordAlongSteptDoIt=ordInActive, G4int ordPostStepDoIt=ordInActive)
void SetElectronLowLimit(G4double)
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
Definition: G4Exception.cc:41
static G4Positron * PositronDefinition()
Definition: G4Positron.cc:89
static G4ProductionCutsTable * GetProductionCutsTable()
void SetPositronLowLimit(G4double)
G4ProcessManager * GetProcessManager() const
std::vector< G4InuclElementaryParticle >::iterator particleIterator
Definition: G4BigBanger.cc:65
void ResetCuts()
obsolete methods
G4ParticleTable::G4PTblDicIterator * GetParticleIterator() const
#define G4endl
Definition: G4ios.hh:61
static constexpr double MeV
Definition: G4SIunits.hh:214
void SetProtonLowLimit(G4double)
static G4OpticalPhoton * OpticalPhotonDefinition()
double G4double
Definition: G4Types.hh:76
static constexpr double micrometer
Definition: G4SIunits.hh:100
G4GLOB_DLL std::ostream G4cerr
static G4Gamma * GammaDefinition()
Definition: G4Gamma.cc:81