Geant4  10.02.p02
UltraPhysicsList.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 //
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 
59 
60 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
61 
63 
64 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
65 
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 
77 
78 }
79 
80 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
81 
83 {
84  // pseudo-particles
87 
88  // gamma
90 
91  // optical photon
93 
94 }
95 
96 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
97 
99 {
100  // leptons
109 }
110 
111 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
112 
114 {
115  // mesons
119 }
120 
121 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
122 
124 {
125 // barions
130 }
131 
132 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
133 
135 {
138  ConstructEM();
139  ConstructOp();
140 
141 }
142 
143 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
144 
145 #include "G4Decay.hh"
146 
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 "G4eMultipleScattering.hh"
167 #include "G4MuMultipleScattering.hh"
168 #include "G4hMultipleScattering.hh"
169 
170 #include "G4eIonisation.hh"
171 #include "G4eBremsstrahlung.hh"
172 #include "G4eplusAnnihilation.hh"
173 
174 #include "G4MuIonisation.hh"
175 #include "G4MuBremsstrahlung.hh"
176 #include "G4MuPairProduction.hh"
177 
178 #include "G4hIonisation.hh"
179 
180 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
181 
183 {
184  theParticleIterator->reset();
185  while( (*theParticleIterator)() ){
186  G4ParticleDefinition* particle = theParticleIterator->value();
187  G4ProcessManager* pmanager = particle->GetProcessManager();
188  G4String particleName = particle->GetParticleName();
189 
190  if (particleName == "gamma") {
191  // gamma
192  // Construct processes for gamma
193  pmanager->AddDiscreteProcess(new G4GammaConversion());
194  pmanager->AddDiscreteProcess(new G4ComptonScattering());
195  pmanager->AddDiscreteProcess(new G4PhotoElectricEffect());
196 
197  } else if (particleName == "e-") {
198  //electron
199  // Construct processes for electron
200  pmanager->AddProcess(new G4eMultipleScattering(),-1,1,1);
201  pmanager->AddProcess(new G4eIonisation(),-1,2,2);
202  pmanager->AddProcess(new G4eBremsstrahlung(),-1,-1,3);
203 
204  } else if (particleName == "e+") {
205  //positron
206  // Construct processes for positron
207  pmanager->AddProcess(new G4eMultipleScattering(),-1,1,1);
208  pmanager->AddProcess(new G4eIonisation(),-1,2,2);
209  pmanager->AddProcess(new G4eBremsstrahlung(),-1,-1,3);
210  pmanager->AddProcess(new G4eplusAnnihilation(),0,-1,4);
211 
212  } else if( particleName == "mu+" ||
213  particleName == "mu-" ) {
214  //muon
215  // Construct processes for muon
216  pmanager->AddProcess(new G4MuMultipleScattering(),-1,1,1);
217  pmanager->AddProcess(new G4MuIonisation(),-1,2,2);
218  pmanager->AddProcess(new G4MuBremsstrahlung(),-1,-1,3);
219  pmanager->AddProcess(new G4MuPairProduction(),-1,-1,4);
220 
221  } else {
222  if ((particle->GetPDGCharge() != 0.0) &&
223  (particle->GetParticleName() != "chargedgeantino")) {
224  // all others charged particles except geantino
225  pmanager->AddProcess(new G4hMultipleScattering(),-1,1,1);
226  pmanager->AddProcess(new G4hIonisation(),-1,2,2);
227  }
228  }
229  }
230 }
231 
232 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
233 
234 #include "G4Cerenkov.hh"
235 #include "G4Scintillation.hh"
236 #include "G4OpAbsorption.hh"
237 #include "G4OpRayleigh.hh"
238 #include "G4OpBoundaryProcess.hh"
239 
241 {
242  // this Cerenkov Process
243  G4Cerenkov* theCerenkovProcess = new G4Cerenkov("Cerenkov");
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  theAbsorptionProcess -> SetVerboseLevel(0);
254  theRayleighScatteringProcess -> SetVerboseLevel(0);
255  theBoundaryProcess -> SetVerboseLevel(0);
256 
257 
258 // Chose MaxNumPhotons that can be generated. Lets ignore this for now
259 // G4int MaxNumPhotons = 300;
260 // theCerenkovProcess->SetMaxNumPhotonsPerStep(MaxNumPhotons);
261  theCerenkovProcess->SetTrackSecondariesFirst(true);
262 
263  theParticleIterator->reset();
264  while( (*theParticleIterator)() ){
265  G4ParticleDefinition* particle = theParticleIterator->value();
266  G4ProcessManager* pmanager = particle->GetProcessManager();
267  G4String particleName = particle->GetParticleName();
268 
269  if (theCerenkovProcess->IsApplicable(*particle)) {
270  pmanager->AddProcess(theCerenkovProcess);
271  pmanager->SetProcessOrdering(theCerenkovProcess,idxPostStep);
272  }
273 
274 
275  if (particleName == "opticalphoton") {
276  G4cout << ">>>>>>>>>>>>>> AddDiscreteProcess to OpticalPhoton " << G4endl;
277  pmanager->AddDiscreteProcess(theAbsorptionProcess);
278  pmanager->AddDiscreteProcess(theRayleighScatteringProcess);
279  pmanager->AddDiscreteProcess(theBoundaryProcess);
280  }
281  }
282 }
283 
284 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
285 
287 {
288  if (verboseLevel >1){
289  G4cout << "UltraPhysicsList::SetCuts:";
290  }
291  // " G4VUserPhysicsList::SetCutsWithDefault" method sets
292  // the default cut value for all particle types
294 }
295 
296 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
static G4Geantino * GeantinoDefinition()
Definition: G4Geantino.cc:82
static G4Electron * ElectronDefinition()
Definition: G4Electron.cc:89
static G4MuonPlus * MuonPlusDefinition()
Definition: G4MuonPlus.cc:94
void SetTrackSecondariesFirst(const G4bool state)
Definition: G4Cerenkov.cc:146
virtual G4bool IsApplicable(const G4ParticleDefinition &)
Definition: G4Decay.cc:89
static G4Proton * ProtonDefinition()
Definition: G4Proton.cc:88
G4int AddDiscreteProcess(G4VProcess *aProcess, G4int ord=ordDefault)
static G4AntiProton * AntiProtonDefinition()
Definition: G4AntiProton.cc:88
G4ProcessManager * GetProcessManager() const
static G4AntiNeutron * AntiNeutronDefinition()
static G4PionZero * PionZeroDefinition()
Definition: G4PionZero.cc:103
const G4String & GetParticleName() const
G4GLOB_DLL std::ostream G4cout
static G4PionMinus * PionMinusDefinition()
Definition: G4PionMinus.cc:93
void SetVerboseLevel(G4int value)
G4int AddProcess(G4VProcess *aProcess, G4int ordAtRestDoIt=ordInActive, G4int ordAlongSteptDoIt=ordInActive, G4int ordPostStepDoIt=ordInActive)
static G4NeutrinoE * NeutrinoEDefinition()
Definition: G4NeutrinoE.cc:80
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 G4MuonMinus * MuonMinusDefinition()
Definition: G4MuonMinus.cc:95
static G4ChargedGeantino * ChargedGeantinoDefinition()
#define G4endl
Definition: G4ios.hh:61
static G4AntiNeutrinoE * AntiNeutrinoEDefinition()
static G4OpticalPhoton * OpticalPhotonDefinition()
G4double GetPDGCharge() const
static G4NeutrinoMu * NeutrinoMuDefinition()
Definition: G4NeutrinoMu.cc:80
static G4Neutron * NeutronDefinition()
Definition: G4Neutron.cc:99
#define theParticleIterator
static G4Gamma * GammaDefinition()
Definition: G4Gamma.cc:81
G4bool IsApplicable(const G4ParticleDefinition &aParticleType)
Definition: G4Cerenkov.cc:135