Geant4_10
PhysicsList.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 //
28 //
29 //
30 // $Id: PhysicsList.cc 67268 2013-02-13 11:38:40Z ihrivnac $
31 //
32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
33 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
34 
35 #include "PhysicsList.hh"
36 #include "PhysicsListMessenger.hh"
37 
38 #include "PhysListEmStandard.hh"
39 #include "PhysListEmLivermore.hh"
40 #include "PhysListEmPenelope.hh"
41 
42 #include "G4LossTableManager.hh"
43 #include "G4UnitsTable.hh"
44 #include "G4SystemOfUnits.hh"
45 
46 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
47 
50  fEmPhysicsList(0),
51  fMessenger(0)
52 {
54 
55  fCurrentDefaultCut = 1.0*mm;
56  fCutForGamma = fCurrentDefaultCut;
57  fCutForElectron = fCurrentDefaultCut;
58  fCutForPositron = fCurrentDefaultCut;
59 
60  fMessenger = new PhysicsListMessenger(this);
61 
62  SetVerboseLevel(1);
63 
64  // EM physics
65  fEmName = G4String("standard");
66  fEmPhysicsList = new PhysListEmStandard(fEmName);
67 
68  //add new units for cross sections
69  //
70  new G4UnitDefinition( "mm2/g", "mm2/g","Surface/Mass", mm2/g);
71  new G4UnitDefinition( "um2/mg", "um2/mg","Surface/Mass", um*um/mg);
72 }
73 
74 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
75 
77 {
78  delete fMessenger;
79 }
80 
81 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
82 
83 // Bosons
84 #include "G4ChargedGeantino.hh"
85 #include "G4Geantino.hh"
86 #include "G4Gamma.hh"
87 #include "G4OpticalPhoton.hh"
88 
89 // leptons
90 #include "G4MuonPlus.hh"
91 #include "G4MuonMinus.hh"
92 #include "G4NeutrinoMu.hh"
93 #include "G4AntiNeutrinoMu.hh"
94 
95 #include "G4Electron.hh"
96 #include "G4Positron.hh"
97 #include "G4NeutrinoE.hh"
98 #include "G4AntiNeutrinoE.hh"
99 
100 // Mesons
101 #include "G4PionPlus.hh"
102 #include "G4PionMinus.hh"
103 #include "G4PionZero.hh"
104 #include "G4Eta.hh"
105 #include "G4EtaPrime.hh"
106 
107 #include "G4KaonPlus.hh"
108 #include "G4KaonMinus.hh"
109 #include "G4KaonZero.hh"
110 #include "G4AntiKaonZero.hh"
111 #include "G4KaonZeroLong.hh"
112 #include "G4KaonZeroShort.hh"
113 
114 // Baryons
115 #include "G4Proton.hh"
116 #include "G4AntiProton.hh"
117 #include "G4Neutron.hh"
118 #include "G4AntiNeutron.hh"
119 
120 // Nuclei
121 #include "G4Deuteron.hh"
122 #include "G4Triton.hh"
123 #include "G4Alpha.hh"
124 #include "G4GenericIon.hh"
125 
126 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
127 
129 {
130 // pseudo-particles
133 
134 // gamma
136 
137 // optical photon
139 
140 // leptons
145 
150 
151 // mesons
163 
164 // barions
169 
170 // ions
175 }
176 
177 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
178 
179 #include "G4EmProcessOptions.hh"
180 
182 {
183  // Transportation
184  //
186 
187  // Electromagnetic physics list
188  //
189  fEmPhysicsList->ConstructProcess();
190 
191  // Em options
192  //
193  G4EmProcessOptions emOptions;
194  emOptions.SetIntegral(false);
195 }
196 
197 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
198 
200 {
201  if (verboseLevel>0) {
202  G4cout << "PhysicsList::AddPhysicsList: <" << name << ">" << G4endl;
203  }
204 
205  if (name == fEmName) return;
206 
207  if (name == "standard") {
208  fEmName = name;
209  delete fEmPhysicsList;
210  fEmPhysicsList = new PhysListEmStandard(name);
211 
212  } else if (name == "livermore") {
213  fEmName = name;
214  delete fEmPhysicsList;
215  fEmPhysicsList = new PhysListEmLivermore(name);
216 
217  } else if (name == "penelope") {
218  fEmName = name;
219  delete fEmPhysicsList;
220  fEmPhysicsList = new PhysListEmPenelope(name);
221 
222  } else {
223  G4cout << "PhysicsList::AddPhysicsList: <" << name << ">"
224  << " is not defined"
225  << G4endl;
226  }
227 }
228 
229 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
230 
231 #include "G4Gamma.hh"
232 #include "G4Electron.hh"
233 #include "G4Positron.hh"
234 
236 {
237  // fixe lower limit for cut
239 
240  // set cut values for gamma at first and for e- second and next for e+,
241  // because some processes for e+/e- need cut values for gamma
242  SetCutValue(fCutForGamma, "gamma");
243  SetCutValue(fCutForElectron, "e-");
244  SetCutValue(fCutForPositron, "e+");
245 }
246 
247 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
248 
250 {
251  fCutForGamma = cut;
252  SetParticleCuts(fCutForGamma, G4Gamma::Gamma());
253 }
254 
255 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
256 
258 {
259  fCutForElectron = cut;
260  SetParticleCuts(fCutForElectron, G4Electron::Electron());
261 }
262 
263 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
264 
266 {
267  fCutForPositron = cut;
268  SetParticleCuts(fCutForPositron, G4Positron::Positron());
269 }
270 
271 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
static G4Geantino * GeantinoDefinition()
Definition: G4Geantino.cc:82
static G4Electron * ElectronDefinition()
Definition: G4Electron.cc:89
static G4Triton * TritonDefinition()
Definition: G4Triton.cc:90
static G4KaonPlus * KaonPlusDefinition()
Definition: G4KaonPlus.cc:108
static G4GenericIon * GenericIonDefinition()
Definition: G4GenericIon.cc:88
void ConstructParticle()
Definition: PhysicsList.cc:117
static G4MuonPlus * MuonPlusDefinition()
Definition: G4MuonPlus.cc:94
static G4LossTableManager * Instance()
void SetCutValue(G4double aCut, const G4String &pname)
void SetCutForGamma(G4double)
Definition: PhysicsList.cc:231
void SetEnergyRange(G4double lowedge, G4double highedge)
static G4KaonMinus * KaonMinusDefinition()
Definition: G4KaonMinus.cc:108
static G4Proton * ProtonDefinition()
Definition: G4Proton.cc:88
void SetCutForPositron(G4double)
Definition: PhysicsList.cc:247
const XML_Char * name
Definition: expat.h:151
static G4KaonZero * KaonZeroDefinition()
Definition: G4KaonZero.cc:99
static G4AntiKaonZero * AntiKaonZeroDefinition()
static G4KaonZeroShort * KaonZeroShortDefinition()
static G4AntiProton * AntiProtonDefinition()
Definition: G4AntiProton.cc:88
void SetCutForElectron(G4double)
Definition: PhysicsList.cc:239
static G4AntiNeutron * AntiNeutronDefinition()
static G4PionZero * PionZeroDefinition()
Definition: G4PionZero.cc:99
void AddPhysicsList(const G4String &name)
Definition: PhysicsList.cc:191
function g(Y1, Y2, PT2)
Definition: hijing1.383.f:5205
void SetParticleCuts(G4double cut, G4ParticleDefinition *particle, G4Region *region=0)
G4GLOB_DLL std::ostream G4cout
static G4PionMinus * PionMinusDefinition()
Definition: G4PionMinus.cc:93
static G4KaonZeroLong * KaonZeroLongDefinition()
static G4NeutrinoE * NeutrinoEDefinition()
Definition: G4NeutrinoE.cc:80
static G4AntiNeutrinoMu * AntiNeutrinoMuDefinition()
static G4PionPlus * PionPlusDefinition()
Definition: G4PionPlus.cc:93
void SetVerboseLevel(G4int value)
static G4Gamma * Gamma()
Definition: G4Gamma.cc:86
void SetCuts()
Definition: PhysicsList.cc:219
static G4Positron * PositronDefinition()
Definition: G4Positron.cc:89
static G4EtaPrime * EtaPrimeDefinition()
Definition: G4EtaPrime.cc:100
static G4ProductionCutsTable * GetProductionCutsTable()
static G4Positron * Positron()
Definition: G4Positron.cc:94
static G4MuonMinus * MuonMinusDefinition()
Definition: G4MuonMinus.cc:95
virtual void ConstructProcess()=0
void SetIntegral(G4bool val)
static G4ChargedGeantino * ChargedGeantinoDefinition()
static G4Electron * Electron()
Definition: G4Electron.cc:94
#define G4endl
Definition: G4ios.hh:61
static G4AntiNeutrinoE * AntiNeutrinoEDefinition()
static G4OpticalPhoton * OpticalPhotonDefinition()
double G4double
Definition: G4Types.hh:76
void ConstructProcess()
Definition: PhysicsList.cc:170
static G4NeutrinoMu * NeutrinoMuDefinition()
Definition: G4NeutrinoMu.cc:80
static G4Deuteron * DeuteronDefinition()
Definition: G4Deuteron.cc:89
static G4Alpha * AlphaDefinition()
Definition: G4Alpha.cc:84
static G4Neutron * NeutronDefinition()
Definition: G4Neutron.cc:99
static G4Eta * EtaDefinition()
Definition: G4Eta.cc:104
static G4Gamma * GammaDefinition()
Definition: G4Gamma.cc:81