source: trunk/examples/advanced/air_shower/src/UltraPhysicsList.cc@ 1194

Last change on this file since 1194 was 807, checked in by garnier, 17 years ago

update
File size: 10.6 KB
Line 
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//
27// --------------------------------------------------------------
28// GEANT 4 - ULTRA experiment example
29// --------------------------------------------------------------
30//
31// Code developed by:
32// B. Tome, M.C. Espirito-Santo, A. Trindade, P. Rodrigues
33//
34// ****************************************************
35// * UltraPhysicsList.cc
36// ****************************************************
37//
38// Ultra Physics List class; Standard and Low Energy EM processes are defined for
39// the relevant particles. Optical processes are declared.
40//
41#include "G4ios.hh"
42//#include "iomanip.h"
43#include "globals.hh"
44
45#include "UltraPhysicsList.hh"
46
47#include "G4ParticleDefinition.hh"
48#include "G4ParticleTypes.hh"
49#include "G4ParticleWithCuts.hh"
50#include "G4ParticleTable.hh"
51#include "G4Material.hh"
52#include "G4MaterialTable.hh"
53#include "G4ProcessManager.hh"
54#include "G4ProcessVector.hh"
55
56//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
57
58UltraPhysicsList::UltraPhysicsList() : G4VUserPhysicsList() {;}
59
60//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
61
62UltraPhysicsList::~UltraPhysicsList() {;}
63
64//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
65
66void UltraPhysicsList::ConstructParticle()
67{
68 // In this method, static member functions should be called
69 // for all particles which you want to use.
70 // This ensures that objects of these particle types will be
71 // created in the program.
72
73 ConstructBosons();
74 ConstructLeptons();
75 ConstructMesons();
76 ConstructBaryons();
77
78}
79
80//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
81
82void UltraPhysicsList::ConstructBosons()
83{
84 // pseudo-particles
85 G4Geantino::GeantinoDefinition();
86 G4ChargedGeantino::ChargedGeantinoDefinition();
87
88 // gamma
89 G4Gamma::GammaDefinition();
90
91 // optical photon
92 G4OpticalPhoton::OpticalPhotonDefinition();
93
94}
95
96//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
97
98void UltraPhysicsList::ConstructLeptons()
99{
100 // leptons
101 G4Electron::ElectronDefinition();
102 G4Positron::PositronDefinition();
103 G4NeutrinoE::NeutrinoEDefinition();
104 G4AntiNeutrinoE::AntiNeutrinoEDefinition();
105 G4MuonPlus::MuonPlusDefinition();
106 G4MuonMinus::MuonMinusDefinition();
107 G4NeutrinoMu::NeutrinoMuDefinition();
108 G4AntiNeutrinoMu::AntiNeutrinoMuDefinition();
109}
110
111//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
112
113void UltraPhysicsList::ConstructMesons()
114{
115 // mesons
116 G4PionPlus::PionPlusDefinition();
117 G4PionMinus::PionMinusDefinition();
118 G4PionZero::PionZeroDefinition();
119}
120
121//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
122
123void UltraPhysicsList::ConstructBaryons()
124{
125// barions
126 G4Proton::ProtonDefinition();
127 G4AntiProton::AntiProtonDefinition();
128 G4Neutron::NeutronDefinition();
129 G4AntiNeutron::AntiNeutronDefinition();
130}
131
132//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
133
134void UltraPhysicsList::ConstructProcess()
135{
136 AddTransportation();
137 ConstructGeneral();
138 ConstructEM();
139 ConstructOp();
140
141}
142
143//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
144
145#include "G4Decay.hh"
146
147void UltraPhysicsList::ConstructGeneral()
148{
149 G4Decay* theDecayProcess = new G4Decay();
150 theParticleIterator->reset();
151 while( (*theParticleIterator)() ){
152 G4ParticleDefinition* particle = theParticleIterator->value();
153 G4ProcessManager* pmanager = particle->GetProcessManager();
154 if (theDecayProcess->IsApplicable(*particle)) {
155 pmanager->AddDiscreteProcess(theDecayProcess);
156 }
157 }
158}
159
160//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
161
162#include "G4ComptonScattering.hh"
163#include "G4GammaConversion.hh"
164#include "G4PhotoElectricEffect.hh"
165
166#include "G4MultipleScattering.hh"
167
168#include "G4eIonisation.hh"
169#include "G4eBremsstrahlung.hh"
170#include "G4eplusAnnihilation.hh"
171
172#include "G4MuIonisation.hh"
173#include "G4MuBremsstrahlung.hh"
174#include "G4MuPairProduction.hh"
175
176#include "G4hIonisation.hh"
177
178//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
179
180void UltraPhysicsList::ConstructEM()
181{
182 theParticleIterator->reset();
183 while( (*theParticleIterator)() ){
184 G4ParticleDefinition* particle = theParticleIterator->value();
185 G4ProcessManager* pmanager = particle->GetProcessManager();
186 G4String particleName = particle->GetParticleName();
187
188 if (particleName == "gamma") {
189 // gamma
190 // Construct processes for gamma
191 pmanager->AddDiscreteProcess(new G4GammaConversion());
192 pmanager->AddDiscreteProcess(new G4ComptonScattering());
193 pmanager->AddDiscreteProcess(new G4PhotoElectricEffect());
194
195 } else if (particleName == "e-") {
196 //electron
197 // Construct processes for electron
198 pmanager->AddProcess(new G4MultipleScattering(),-1,1,1);
199 pmanager->AddProcess(new G4eIonisation(),-1,2,2);
200 pmanager->AddProcess(new G4eBremsstrahlung(),-1,-1,3);
201
202 } else if (particleName == "e+") {
203 //positron
204 // Construct processes for positron
205 pmanager->AddProcess(new G4MultipleScattering(),-1,1,1);
206 pmanager->AddProcess(new G4eIonisation(),-1,2,2);
207 pmanager->AddProcess(new G4eBremsstrahlung(),-1,-1,3);
208 pmanager->AddProcess(new G4eplusAnnihilation(),0,-1,4);
209
210 } else if( particleName == "mu+" ||
211 particleName == "mu-" ) {
212 //muon
213 // Construct processes for muon
214 pmanager->AddProcess(new G4MultipleScattering(),-1,1,1);
215 pmanager->AddProcess(new G4MuIonisation(),-1,2,2);
216 pmanager->AddProcess(new G4MuBremsstrahlung(),-1,-1,3);
217 pmanager->AddProcess(new G4MuPairProduction(),-1,-1,4);
218
219 } else {
220 if ((particle->GetPDGCharge() != 0.0) &&
221 (particle->GetParticleName() != "chargedgeantino")) {
222 // all others charged particles except geantino
223 pmanager->AddProcess(new G4MultipleScattering(),-1,1,1);
224 pmanager->AddProcess(new G4hIonisation(),-1,2,2);
225 }
226 }
227 }
228}
229
230//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
231
232#include "G4Cerenkov.hh"
233#include "G4Scintillation.hh"
234#include "G4OpAbsorption.hh"
235#include "G4OpRayleigh.hh"
236#include "G4OpBoundaryProcess.hh"
237
238void UltraPhysicsList::ConstructOp()
239{
240 // this Cerenkov Process
241 G4Cerenkov* theCerenkovProcess = new G4Cerenkov("Cerenkov");
242 // this is Scintillation process
243 G4Scintillation* theScintillationProcess = new G4Scintillation("Scintillation");
244 // this absorption process inside optical media
245 G4OpAbsorption* theAbsorptionProcess = new G4OpAbsorption();
246 // Rayleigh scattering for optical photons (aerogel radiators)
247 G4OpRayleigh* theRayleighScatteringProcess = new G4OpRayleigh();
248 // Boundary process definition Class
249 G4OpBoundaryProcess* theBoundaryProcess = new G4OpBoundaryProcess();
250
251 // Chose level 0 (no verbose)
252 theCerenkovProcess -> SetVerboseLevel(0);
253 theScintillationProcess -> SetVerboseLevel(0);
254 theAbsorptionProcess -> SetVerboseLevel(0);
255 theRayleighScatteringProcess -> SetVerboseLevel(0);
256 theBoundaryProcess -> SetVerboseLevel(0);
257
258
259// Chose MaxNumPhotons that can be generated. Lets ignore this for now
260// G4int MaxNumPhotons = 300;
261// theCerenkovProcess->SetMaxNumPhotonsPerStep(MaxNumPhotons);
262 theCerenkovProcess->SetTrackSecondariesFirst(true);
263
264 theScintillationProcess->SetTrackSecondariesFirst(true);
265 theScintillationProcess->SetScintillationYieldFactor(1.);
266 theScintillationProcess->SetScintillationExcitationRatio(0.0);
267
268 // Boundary model (UNIFIED OR GLISUR (OLD GEANT3)) For now only GEANT3
269 G4OpticalSurfaceModel themodel = unified;
270 theBoundaryProcess->SetModel(themodel);
271
272 theParticleIterator->reset();
273 while( (*theParticleIterator)() ){
274 G4ParticleDefinition* particle = theParticleIterator->value();
275 G4ProcessManager* pmanager = particle->GetProcessManager();
276 G4String particleName = particle->GetParticleName();
277
278 if (theCerenkovProcess->IsApplicable(*particle)) {
279 // pmanager->AddProcess(theCerenkovProcess);
280 // pmanager->SetProcessOrdering(theCerenkovProcess,idxPostStep);
281 }
282
283 if (theScintillationProcess->IsApplicable(*particle)) {
284 pmanager->AddProcess(theScintillationProcess);
285 pmanager->SetProcessOrderingToLast(theScintillationProcess, idxAtRest);
286 pmanager->SetProcessOrderingToLast(theScintillationProcess, idxPostStep);
287 }
288
289 if (particleName == "opticalphoton") {
290 G4cout << ">>>>>>>>>>>>>> AddDiscreteProcess to OpticalPhoton " << G4endl;
291 pmanager->AddDiscreteProcess(theAbsorptionProcess);
292 pmanager->AddDiscreteProcess(theRayleighScatteringProcess);
293 pmanager->AddDiscreteProcess(theBoundaryProcess);
294 }
295 }
296}
297
298//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
299
300void UltraPhysicsList::SetCuts()
301{
302 if (verboseLevel >1){
303 G4cout << "UltraPhysicsList::SetCuts:";
304 }
305 // " G4VUserPhysicsList::SetCutsWithDefault" method sets
306 // the default cut value for all particle types
307 SetCutsWithDefault();
308}
309
310//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
Note: See TracBrowser for help on using the repository browser.