Geant4  10.02.p01
PhysListEmPenelope.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: PhysListEmPenelope.cc 68585 2013-04-01 23:35:07Z adotti $
31 //
32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
33 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
34 
35 #include "PhysListEmPenelope.hh"
36 #include "G4ParticleDefinition.hh"
37 #include "G4ProcessManager.hh"
38 
39 // gamma
40 
41 #include "G4PhotoElectricEffect.hh"
43 
44 #include "G4ComptonScattering.hh"
46 
47 #include "G4GammaConversion.hh"
49 
50 #include "G4RayleighScattering.hh"
52 
53 // e-
54 
55 #include "G4eIonisation.hh"
58 
59 #include "G4eBremsstrahlung.hh"
61 
62 // e+
63 
64 #include "G4eplusAnnihilation.hh"
66 
67 // mu
68 
69 #include "G4MuIonisation.hh"
70 #include "G4MuBremsstrahlung.hh"
71 #include "G4MuPairProduction.hh"
72 
73 // hadrons, ions
74 
75 #include "G4hIonisation.hh"
76 #include "G4ionIonisation.hh"
77 
78 #include "G4SystemOfUnits.hh"
79 
80 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
81 
83  : G4VPhysicsConstructor(name)
84 { }
85 
86 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
87 
89 { }
90 
91 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
92 
94 {
95  // Add standard EM Processes
96 
97  aParticleIterator->reset();
98  while( (*aParticleIterator)() ){
99  G4ParticleDefinition* particle = aParticleIterator->value();
100  G4ProcessManager* pmanager = particle->GetProcessManager();
101  G4String particleName = particle->GetParticleName();
102 
103  //Applicability range for Penelope models
104  //for higher energies, the Standard models are used
105  G4double highEnergyLimit = 1*GeV;
106 
107  if (particleName == "gamma") {
108  // gamma
109 
112  photModel = new G4PenelopePhotoElectricModel();
113  photModel->SetHighEnergyLimit(highEnergyLimit);
114  phot->AddEmModel(0, photModel);
115  pmanager->AddDiscreteProcess(phot);
116 
119  comptModel = new G4PenelopeComptonModel();
120  comptModel->SetHighEnergyLimit(highEnergyLimit);
121  compt->AddEmModel(0, comptModel);
122  pmanager->AddDiscreteProcess(compt);
123 
124  G4GammaConversion* conv = new G4GammaConversion();
126  convModel = new G4PenelopeGammaConversionModel();
127  convModel->SetHighEnergyLimit(highEnergyLimit);
128  conv->AddEmModel(0, convModel);
129  pmanager->AddDiscreteProcess(conv);
130 
133  raylModel = new G4PenelopeRayleighModel();
134  raylModel->SetHighEnergyLimit(highEnergyLimit);
135  rayl->AddEmModel(0, raylModel);
136  pmanager->AddDiscreteProcess(rayl);
137 
138  } else if (particleName == "e-") {
139  //electron
140 
141  G4eIonisation* eIoni = new G4eIonisation();
143  eIoniModel = new G4PenelopeIonisationModel();
144  eIoniModel->SetHighEnergyLimit(highEnergyLimit);
145  eIoni->AddEmModel(0, eIoniModel, new G4UniversalFluctuation() );
146  pmanager->AddProcess(eIoni, -1,-1, 1);
147 
148  G4eBremsstrahlung* eBrem = new G4eBremsstrahlung();
150  eBremModel = new G4PenelopeBremsstrahlungModel();
151  eBremModel->SetHighEnergyLimit(highEnergyLimit);
152  eBrem->AddEmModel(0, eBremModel);
153  pmanager->AddProcess(eBrem, -1,-1, 2);
154 
155  } else if (particleName == "e+") {
156  //positron
157  G4eIonisation* eIoni = new G4eIonisation();
159  eIoniModel = new G4PenelopeIonisationModel();
160  eIoniModel->SetHighEnergyLimit(highEnergyLimit);
161  eIoni->AddEmModel(0, eIoniModel, new G4UniversalFluctuation() );
162  pmanager->AddProcess(eIoni, -1,-1, 1);
163 
164  G4eBremsstrahlung* eBrem = new G4eBremsstrahlung();
166  eBremModel = new G4PenelopeBremsstrahlungModel();
167  eBremModel->SetHighEnergyLimit(highEnergyLimit);
168  eBrem->AddEmModel(0, eBremModel);
169  pmanager->AddProcess(eBrem, -1,-1, 2);
170 
173  eAnniModel = new G4PenelopeAnnihilationModel();
174  eAnniModel->SetHighEnergyLimit(highEnergyLimit);
175  eAnni->AddEmModel(0, eAnniModel);
176  pmanager->AddProcess(eAnni, 0,-1, 3);
177 
178  } else if( particleName == "mu+" ||
179  particleName == "mu-" ) {
180  //muon
181  pmanager->AddProcess(new G4MuIonisation, -1,-1, 1);
182  pmanager->AddProcess(new G4MuBremsstrahlung, -1,-1, 2);
183  pmanager->AddProcess(new G4MuPairProduction, -1,-1, 3);
184 
185  } else if( particleName == "alpha" || particleName == "GenericIon" ) {
186  pmanager->AddProcess(new G4ionIonisation, -1,-1, 1);
187 
188  } else if ((!particle->IsShortLived()) &&
189  (particle->GetPDGCharge() != 0.0) &&
190  (particle->GetParticleName() != "chargedgeantino")) {
191  //all others charged particles except geantino
192  pmanager->AddProcess(new G4hIonisation, -1,-1, 1);
193  }
194  }
195 }
196 
197 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
198 
G4String name
Definition: TRTMaterials.hh:40
G4int AddDiscreteProcess(G4VProcess *aProcess, G4int ord=ordDefault)
PhysListEmPenelope(const G4String &name="penelope")
G4ProcessManager * GetProcessManager() const
const G4String & GetParticleName() const
void SetHighEnergyLimit(G4double)
Definition: G4VEmModel.hh:725
#define aParticleIterator
G4int AddProcess(G4VProcess *aProcess, G4int ordAtRestDoIt=ordInActive, G4int ordAlongSteptDoIt=ordInActive, G4int ordPostStepDoIt=ordInActive)
static const double GeV
Definition: G4SIunits.hh:214
void AddEmModel(G4int, G4VEmModel *, G4VEmFluctuationModel *fluc=0, const G4Region *region=0)
virtual void ConstructProcess()
void AddEmModel(G4int, G4VEmModel *, const G4Region *region=0)
double G4double
Definition: G4Types.hh:76
G4double GetPDGCharge() const