Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
PhysListEmLivermore.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: PhysListEmLivermore.cc 100279 2016-10-17 08:34:54Z gcosmo $
31 //
32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
33 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
34 
35 #include "PhysListEmLivermore.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"
65 
66 // mu
67 
68 #include "G4MuIonisation.hh"
69 #include "G4MuBremsstrahlung.hh"
70 #include "G4MuPairProduction.hh"
71 
72 // hadrons, ions
73 
74 #include "G4hIonisation.hh"
75 #include "G4ionIonisation.hh"
76 
77 #include "G4SystemOfUnits.hh"
78 
79 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
80 
82  : G4VPhysicsConstructor(name)
83 { }
84 
85 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
86 
88 { }
89 
90 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
91 
93 {
94  // Add standard EM Processes
95 
97  particleIterator->reset();
98  while( (*particleIterator)() ){
99  G4ParticleDefinition* particle = particleIterator->value();
100  G4ProcessManager* pmanager = particle->GetProcessManager();
101  G4String particleName = particle->GetParticleName();
102 
103  //Applicability range for Livermore 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 G4LivermorePhotoElectricModel();
113  photModel->SetHighEnergyLimit(highEnergyLimit);
114  phot->AddEmModel(0, photModel);
115  pmanager->AddDiscreteProcess(phot);
116 
119  comptModel = new G4LivermoreComptonModel();
120  comptModel->SetHighEnergyLimit(highEnergyLimit);
121  compt->AddEmModel(0, comptModel);
122  pmanager->AddDiscreteProcess(compt);
123 
124  G4GammaConversion* conv = new G4GammaConversion();
126  convModel = new G4LivermoreGammaConversionModel();
127  convModel->SetHighEnergyLimit(highEnergyLimit);
128  conv->AddEmModel(0, convModel);
129  pmanager->AddDiscreteProcess(conv);
130 
133  raylModel = new G4LivermoreRayleighModel();
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 G4LivermoreIonisationModel();
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 G4LivermoreBremsstrahlungModel();
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  pmanager->AddProcess(new G4eIonisation, -1,-1, 1);
158  pmanager->AddProcess(new G4eBremsstrahlung, -1,-1, 2);
159  pmanager->AddProcess(new G4eplusAnnihilation, 0,-1, 3);
160 
161  } else if( particleName == "mu+" ||
162  particleName == "mu-" ) {
163  //muon
164  pmanager->AddProcess(new G4MuIonisation, -1,-1, 1);
165  pmanager->AddProcess(new G4MuBremsstrahlung, -1,-1, 2);
166  pmanager->AddProcess(new G4MuPairProduction, -1,-1, 3);
167 
168  } else if( particleName == "alpha" || particleName == "GenericIon" ) {
169  pmanager->AddProcess(new G4ionIonisation, -1,-1, 1);
170 
171  } else if ((!particle->IsShortLived()) &&
172  (particle->GetPDGCharge() != 0.0) &&
173  (particle->GetParticleName() != "chargedgeantino")) {
174  //all others charged particles except geantino
175  pmanager->AddProcess(new G4hIonisation, -1,-1, 1);
176  }
177  }
178 }
179 
180 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
181 
const XML_Char * name
Definition: expat.h:151
virtual void ConstructProcess()
G4int AddDiscreteProcess(G4VProcess *aProcess, G4int ord=ordDefault)
const G4String & GetParticleName() const
void SetHighEnergyLimit(G4double)
Definition: G4VEmModel.hh:732
void AddEmModel(G4int, G4VEmModel *, G4VEmFluctuationModel *fluc=0, const G4Region *region=nullptr)
G4ParticleTable::G4PTblDicIterator * GetParticleIterator() const
G4int AddProcess(G4VProcess *aProcess, G4int ordAtRestDoIt=ordInActive, G4int ordAlongSteptDoIt=ordInActive, G4int ordPostStepDoIt=ordInActive)
void AddEmModel(G4int, G4VEmModel *, const G4Region *region=nullptr)
G4ProcessManager * GetProcessManager() const
static constexpr double GeV
Definition: G4SIunits.hh:217
std::vector< G4InuclElementaryParticle >::iterator particleIterator
Definition: G4BigBanger.cc:65
double G4double
Definition: G4Types.hh:76
G4double GetPDGCharge() const
PhysListEmLivermore(const G4String &name="livermore")