source: trunk/examples/extended/electromagnetic/TestEm6/src/PhysicsList.cc

Last change on this file was 1337, checked in by garnier, 14 years ago

tag geant4.9.4 beta 1 + modifs locales
File size: 11.0 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// $Id: PhysicsList.cc,v 1.15 2010/03/31 09:26:09 vnivanch Exp $
27// GEANT4 tag $Name: geant4-09-04-beta-01 $
28//
29//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
30//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
31
32#include "PhysicsList.hh"
33#include "PhysicsListMessenger.hh"
34
35#include "G4ParticleDefinition.hh"
36#include "G4ProcessManager.hh"
37#include "G4ParticleTypes.hh"
38#include "G4ParticleTable.hh"
39
40#include "G4Proton.hh"
41#include "G4AntiProton.hh"
42#include "G4Neutron.hh"
43#include "G4AntiNeutron.hh"
44
45#include "G4GammaConversionToMuons.hh"
46
47#include "G4eMultipleScattering.hh"
48#include "G4MuMultipleScattering.hh"
49#include "G4hMultipleScattering.hh"
50
51#include "G4eIonisation.hh"
52#include "G4eBremsstrahlung.hh"
53#include "G4eplusAnnihilation.hh"
54#include "G4AnnihiToMuPair.hh"
55
56#include "G4MuIonisation.hh"
57#include "G4MuBremsstrahlung.hh"
58#include "G4MuPairProduction.hh"
59
60#include "G4hIonisation.hh"
61#include "G4hhIonisation.hh"
62#include "G4ionIonisation.hh"
63#include "G4eeToHadrons.hh"
64
65#include "G4Decay.hh"
66#include "G4EmProcessOptions.hh"
67
68#include "G4StepLimiter.hh"
69
70#include "G4ProcessTable.hh"
71
72#include "G4ComptonScattering.hh"
73#include "G4GammaConversion.hh"
74#include "G4PhotoElectricEffect.hh"
75
76
77//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
78
79PhysicsList::PhysicsList()
80: G4VUserPhysicsList()
81{
82  defaultCutValue = 1.*km;
83  pMes = new PhysicsListMessenger(this);
84  SetVerboseLevel(1);
85}
86
87//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
88
89PhysicsList::~PhysicsList()
90{
91  delete pMes;
92}
93
94//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
95
96void PhysicsList::ConstructParticle()
97{
98  // In this method, static member functions should be called
99  // for all particles which you want to use.
100  // This ensures that objects of these particle types will be
101  // created in the program.
102 
103  ConstructBosons();
104  ConstructLeptons();
105  ConstructHadrons();
106}
107
108//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
109
110void PhysicsList::ConstructBosons()
111{
112  // pseudo-particles
113  G4Geantino::GeantinoDefinition();
114  G4ChargedGeantino::ChargedGeantinoDefinition();
115
116  // gamma
117  G4Gamma::GammaDefinition();
118}
119
120//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
121
122void PhysicsList::ConstructLeptons()
123{
124  // leptons
125  G4Electron::ElectronDefinition();
126  G4Positron::PositronDefinition();
127  G4MuonPlus::MuonPlusDefinition();
128  G4MuonMinus::MuonMinusDefinition();
129
130  G4NeutrinoE::NeutrinoEDefinition();
131  G4AntiNeutrinoE::AntiNeutrinoEDefinition();
132  G4NeutrinoMu::NeutrinoMuDefinition();
133  G4AntiNeutrinoMu::AntiNeutrinoMuDefinition();
134}
135//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
136
137void PhysicsList::ConstructHadrons()
138{
139  // mesons
140  G4PionPlus::PionPlusDefinition();
141  G4PionMinus::PionMinusDefinition();
142  G4PionZero::PionZeroDefinition();
143  G4Eta::EtaDefinition();
144  G4EtaPrime::EtaPrimeDefinition();
145  G4KaonPlus::KaonPlusDefinition();
146  G4KaonMinus::KaonMinusDefinition();
147  G4KaonZero::KaonZeroDefinition();
148  G4AntiKaonZero::AntiKaonZeroDefinition();
149  G4KaonZeroLong::KaonZeroLongDefinition();
150  G4KaonZeroShort::KaonZeroShortDefinition();
151
152  // baryons
153  G4Proton::ProtonDefinition();
154  G4AntiProton::AntiProtonDefinition();
155  G4Neutron::NeutronDefinition();
156  G4AntiNeutron::AntiNeutronDefinition();
157}
158
159//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
160
161void PhysicsList::ConstructProcess()
162{
163  AddTransportation();
164  ConstructEM();
165  ConstructGeneral();
166}
167
168//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
169
170void PhysicsList::ConstructEM()
171{
172  theParticleIterator->reset();
173  while( (*theParticleIterator)() ){
174    G4ParticleDefinition* particle = theParticleIterator->value();
175    G4ProcessManager* pmanager = particle->GetProcessManager();
176    G4String particleName = particle->GetParticleName();
177
178    if (particleName == "gamma") {
179      // gamma    allow  only   gamma -> mu+mu-
180      pmanager->AddDiscreteProcess(new G4GammaConversionToMuons);
181      pmanager->AddDiscreteProcess(new G4PhotoElectricEffect);
182      pmanager->AddDiscreteProcess(new G4ComptonScattering);
183      pmanager->AddDiscreteProcess(new G4GammaConversion);
184         
185    } else if (particleName == "e-") {
186      //electron
187      pmanager->AddProcess(new G4eMultipleScattering,-1, 1,1);
188      pmanager->AddProcess(new G4eIonisation,       -1, 2,2);
189      pmanager->AddProcess(new G4eBremsstrahlung,   -1, 3,3);
190      pmanager->AddProcess(new G4StepLimiter,       -1,-1,4);     
191
192    } else if (particleName == "e+") {
193      //positron
194      // to make the process of e+e- annihilation more visible,
195      // do not enable the other standard processes:
196      //pmanager->AddProcess(new G4eMultipleScattering,-1, 1,1);
197      //pmanager->AddProcess(new G4eIonisation,       -1, 2,2);
198      //pmanager->AddProcess(new G4eBremsstrahlung,   -1, 3,3);
199      //pmanager->AddProcess(new G4eplusAnnihilation,  0,-1,4);
200
201      pmanager->AddDiscreteProcess(new G4AnnihiToMuPair);
202      pmanager->AddDiscreteProcess(new G4eeToHadrons);
203      pmanager->AddDiscreteProcess(new G4StepLimiter);
204       
205    } else if( particleName == "mu+" ||
206               particleName == "mu-"    ) {
207      //muon
208      pmanager->AddProcess(new G4MuMultipleScattering,-1, 1,1);
209      pmanager->AddProcess(new G4MuIonisation,      -1, 2,2);
210      pmanager->AddProcess(new G4MuBremsstrahlung,  -1, 3,3);
211      pmanager->AddProcess(new G4MuPairProduction,  -1, 4,4);
212      pmanager->AddProcess(new G4StepLimiter,       -1,-1,5);           
213     
214    } else if( particleName == "anti_proton") {
215      pmanager->AddProcess(new G4hMultipleScattering,-1, 1,1);
216      pmanager->AddProcess(new G4hhIonisation,      -1, 2,2);
217      pmanager->AddProcess(new G4StepLimiter,       -1,-1,3);       
218
219    } else if( particleName == "GenericIon") {
220      pmanager->AddProcess(new G4hMultipleScattering,-1, 1,1);
221      pmanager->AddProcess(new G4ionIonisation,     -1, 2,2);
222      pmanager->AddProcess(new G4StepLimiter,       -1,-1,3);       
223     
224    } else if( particle->GetPDGCharge() != 0.0 && !particle->IsShortLived()
225            && particleName != "chargedgeantino") {
226      pmanager->AddProcess(new G4hMultipleScattering,-1, 1,1);
227      pmanager->AddProcess(new G4hIonisation,       -1, 2,2);
228      pmanager->AddProcess(new G4StepLimiter,       -1,-1,3);       
229    }
230  }
231 
232  G4EmProcessOptions opt;
233  opt.SetVerbose(1);
234  opt.SetMinEnergy(100*eV);       
235  opt.SetMaxEnergy(1000*TeV);     
236  opt.SetDEDXBinning(13*7);     
237  opt.SetLambdaBinning(13*7);   
238}
239
240//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
241
242void PhysicsList::ConstructGeneral()
243{
244  // Add Decay Process
245  G4Decay* theDecayProcess = new G4Decay();
246  theParticleIterator->reset();
247  while ((*theParticleIterator)()){
248      G4ParticleDefinition* particle = theParticleIterator->value();
249      G4ProcessManager* pmanager = particle->GetProcessManager();
250      if (theDecayProcess->IsApplicable(*particle) && !particle->IsShortLived()) {
251        pmanager ->AddProcess(theDecayProcess);
252        // set ordering for PostStepDoIt and AtRestDoIt
253        pmanager ->SetProcessOrdering(theDecayProcess, idxPostStep);
254        pmanager ->SetProcessOrdering(theDecayProcess, idxAtRest);
255      }
256  }
257}
258
259//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
260
261void PhysicsList::SetCuts()
262{
263  if (verboseLevel >0){
264    G4cout << "PhysicsList::SetCuts:";
265    G4cout << "CutLength : " << G4BestUnit(defaultCutValue,"Length") << G4endl;
266  } 
267
268  // set cut values for gamma at first and for e- second and next for e+,
269  // because some processes for e+/e- need cut values for gamma
270  SetCutValue(defaultCutValue, "gamma");
271  SetCutValue(defaultCutValue, "e-");
272  SetCutValue(defaultCutValue, "e+");
273
274  if (verboseLevel>0) DumpCutValuesTable();
275}
276
277//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
278
279void PhysicsList::SetGammaToMuPairFac(G4double fac)
280{
281  G4ProcessTable* theProcessTable = G4ProcessTable::GetProcessTable();
282  G4GammaConversionToMuons* theGammaToMuPairProcess = (G4GammaConversionToMuons*) theProcessTable->FindProcess("GammaToMuPair","gamma");
283  if(theGammaToMuPairProcess) theGammaToMuPairProcess->SetCrossSecFactor(fac);
284  else G4cout << "Warning. No process GammaToMuPair found, SetGammaToMuPairFac was ignored" << G4endl;
285}
286
287//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
288
289void PhysicsList::SetAnnihiToMuPairFac(G4double fac)
290{
291  G4ProcessTable* theProcessTable = G4ProcessTable::GetProcessTable();
292  G4AnnihiToMuPair* theAnnihiToMuPairProcess = (G4AnnihiToMuPair*) theProcessTable->FindProcess("AnnihiToMuPair","e+");
293  if(theAnnihiToMuPairProcess) theAnnihiToMuPairProcess->SetCrossSecFactor(fac);
294  else G4cout << "Warning. No process AnnihiToMuPair found, SetAnnihiToMuPairFac was ignored" << G4endl;
295}
296
297//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
298
299void PhysicsList::SetAnnihiToHadronFac(G4double fac)
300{
301  G4ProcessTable* theProcessTable = G4ProcessTable::GetProcessTable();
302  G4eeToHadrons* eehadProcess = (G4eeToHadrons*) theProcessTable->FindProcess("ee2hadr","e+");
303  if(eehadProcess) eehadProcess->SetCrossSecFactor(fac);
304  else G4cout << "Warning. No process ee2hadr found, SetAnnihiToHadronFac was ignored" << G4endl;
305}
306
307//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
Note: See TracBrowser for help on using the repository browser.