Geant4_10
ExP01PhysicsList.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 // $Id: ExP01PhysicsList.cc 68025 2013-03-13 13:43:46Z gcosmo $
30 //
31 //
32 //
33 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
34 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
35 
36 #include "globals.hh"
37 #include "ExP01PhysicsList.hh"
38 
39 #include "G4ProcessManager.hh"
40 #include "G4ParticleTypes.hh"
41 #include "G4SystemOfUnits.hh"
42 
43 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
44 
46 {
47  defaultCutValue = 1.0*cm;
48  SetVerboseLevel(1);
49 }
50 
51 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
52 
54 {}
55 
56 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
57 
59 {
60  // In this method, static member functions should be called
61  // for all particles which you want to use.
62  // This ensures that objects of these particle types will be
63  // created in the program.
64 
69 }
70 
71 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
72 
74 {
75  // pseudo-particles
78 
79  // gamma
81 }
82 
83 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
84 
86 {
87  // leptons
88  // e+/-
91  // mu+/-
94  // nu_e
97  // nu_mu
100 }
101 
102 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
103 
105 {
106  // mesons
107  // light mesons
119 }
120 
121 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
122 
124 {
125  // barions
128 
131 }
132 
133 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
134 
136 {
138  ConstructEM();
140 }
141 
142 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
143 
144 #include "G4ComptonScattering.hh"
145 #include "G4GammaConversion.hh"
146 #include "G4PhotoElectricEffect.hh"
147 
148 #include "G4eMultipleScattering.hh"
149 #include "G4eIonisation.hh"
150 #include "G4eBremsstrahlung.hh"
151 #include "G4eplusAnnihilation.hh"
152 
153 #include "G4MuMultipleScattering.hh"
154 #include "G4MuIonisation.hh"
155 #include "G4MuBremsstrahlung.hh"
156 #include "G4MuPairProduction.hh"
157 
158 #include "G4hMultipleScattering.hh"
159 #include "G4hIonisation.hh"
160 
161 #include "G4StepLimiter.hh"
162 #include "G4UserSpecialCuts.hh"
163 
164 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
165 
167 {
168  theParticleIterator->reset();
169  while( (*theParticleIterator)() ){
170  G4ParticleDefinition* particle = theParticleIterator->value();
171  G4ProcessManager* pmanager = particle->GetProcessManager();
172  G4String particleName = particle->GetParticleName();
173 
174  if (particleName == "gamma") {
175  // gamma
178  pmanager->AddDiscreteProcess(new G4GammaConversion);
179 
180  } else if (particleName == "e-") {
181  //electron
182  pmanager->AddProcess(new G4eMultipleScattering,-1, 1,1);
183  pmanager->AddProcess(new G4eIonisation, -1, 2,2);
184  pmanager->AddProcess(new G4eBremsstrahlung, -1, 3,3);
185 
186  } else if (particleName == "e+") {
187  //positron
188  pmanager->AddProcess(new G4eMultipleScattering,-1, 1,1);
189  pmanager->AddProcess(new G4eIonisation, -1, 2,2);
190  pmanager->AddProcess(new G4eBremsstrahlung, -1, 3,3);
191  pmanager->AddProcess(new G4eplusAnnihilation, 0,-1,4);
192 
193  } else if( particleName == "mu+" ||
194  particleName == "mu-" ) {
195  //muon
196  pmanager->AddProcess(new G4MuMultipleScattering,-1, 1,1);
197  pmanager->AddProcess(new G4MuIonisation, -1, 2,2);
198  pmanager->AddProcess(new G4MuBremsstrahlung, -1, 3,3);
199  pmanager->AddProcess(new G4MuPairProduction, -1, 4,4);
200 
201  } else if ((!particle->IsShortLived()) &&
202  (particle->GetPDGCharge() != 0.0) &&
203  (particle->GetParticleName() != "chargedgeantino")) {
204  //all others charged particles except geantino
205  pmanager->AddProcess(new G4hMultipleScattering,-1, 1,1);
206  pmanager->AddProcess(new G4hIonisation, -1, 2,2);
207  //step limit
208  pmanager->AddProcess(new G4StepLimiter, -1,-1,3);
210  }
211  }
212 }
213 
214 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
215 
216 #include "G4Decay.hh"
218 {
219  // Add Decay Process
220  G4Decay* theDecayProcess = new G4Decay();
221  theParticleIterator->reset();
222  while( (*theParticleIterator)() ){
223  G4ParticleDefinition* particle = theParticleIterator->value();
224  G4ProcessManager* pmanager = particle->GetProcessManager();
225  if (theDecayProcess->IsApplicable(*particle)) {
226  pmanager ->AddProcess(theDecayProcess);
227  // set ordering for PostStepDoIt and AtRestDoIt
228  pmanager ->SetProcessOrdering(theDecayProcess, idxPostStep);
229  pmanager ->SetProcessOrdering(theDecayProcess, idxAtRest);
230  }
231  }
232 }
233 
234 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
235 
237 {
238  //G4VUserPhysicsList::SetCutsWithDefault method sets
239  //the default cut value for all particle types
240  //
242 
244 }
245 
246 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
247 
static G4Geantino * GeantinoDefinition()
Definition: G4Geantino.cc:82
static G4Electron * ElectronDefinition()
Definition: G4Electron.cc:89
static G4KaonPlus * KaonPlusDefinition()
Definition: G4KaonPlus.cc:108
static G4MuonPlus * MuonPlusDefinition()
Definition: G4MuonPlus.cc:94
virtual G4bool IsApplicable(const G4ParticleDefinition &)
Definition: G4Decay.cc:89
static G4KaonMinus * KaonMinusDefinition()
Definition: G4KaonMinus.cc:108
static G4Proton * ProtonDefinition()
Definition: G4Proton.cc:88
G4int AddDiscreteProcess(G4VProcess *aProcess, G4int ord=ordDefault)
static G4KaonZero * KaonZeroDefinition()
Definition: G4KaonZero.cc:99
static G4AntiKaonZero * AntiKaonZeroDefinition()
static G4KaonZeroShort * KaonZeroShortDefinition()
static G4AntiProton * AntiProtonDefinition()
Definition: G4AntiProton.cc:88
G4ProcessManager * GetProcessManager() const
static G4AntiNeutron * AntiNeutronDefinition()
static G4PionZero * PionZeroDefinition()
Definition: G4PionZero.cc:99
const G4String & GetParticleName() const
void DumpCutValuesTable(G4int flag=1)
static G4PionMinus * PionMinusDefinition()
Definition: G4PionMinus.cc:93
static G4KaonZeroLong * KaonZeroLongDefinition()
void SetVerboseLevel(G4int value)
G4int AddProcess(G4VProcess *aProcess, G4int ordAtRestDoIt=ordInActive, G4int ordAlongSteptDoIt=ordInActive, G4int ordPostStepDoIt=ordInActive)
static G4NeutrinoE * NeutrinoEDefinition()
Definition: G4NeutrinoE.cc:80
virtual void ConstructParticle()
static G4AntiNeutrinoMu * AntiNeutrinoMuDefinition()
static G4PionPlus * PionPlusDefinition()
Definition: G4PionPlus.cc:93
void SetProcessOrdering(G4VProcess *aProcess, G4ProcessVectorDoItIndex idDoIt, G4int ordDoIt=ordDefault)
static G4Positron * PositronDefinition()
Definition: G4Positron.cc:89
static G4EtaPrime * EtaPrimeDefinition()
Definition: G4EtaPrime.cc:100
Definition of the ExP01PhysicsList class.
virtual void ConstructProcess()
static G4MuonMinus * MuonMinusDefinition()
Definition: G4MuonMinus.cc:95
static G4ChargedGeantino * ChargedGeantinoDefinition()
virtual void SetCuts()
static G4AntiNeutrinoE * AntiNeutrinoEDefinition()
G4double GetPDGCharge() const
static G4NeutrinoMu * NeutrinoMuDefinition()
Definition: G4NeutrinoMu.cc:80
static G4Neutron * NeutronDefinition()
Definition: G4Neutron.cc:99
#define theParticleIterator
static G4Eta * EtaDefinition()
Definition: G4Eta.cc:104
static G4Gamma * GammaDefinition()
Definition: G4Gamma.cc:81