source: trunk/examples/extended/analysis/AnaEx01/src/AnaEx01PhysicsList.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: 12.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// $Id: AnaEx01PhysicsList.cc,v 1.5 2006/06/29 16:33:55 gunter Exp $
28// GEANT4 tag $Name: geant4-09-04-beta-01 $
29//
30//
31
32//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
33//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
34
35#include "AnaEx01PhysicsList.hh"
36
37#include "G4ParticleDefinition.hh"
38#include "G4ParticleWithCuts.hh"
39#include "G4ProcessManager.hh"
40#include "G4ProcessVector.hh"
41#include "G4ParticleTypes.hh"
42#include "G4ParticleTable.hh"
43#include "G4Material.hh"
44#include "G4ios.hh"             
45
46//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
47
48AnaEx01PhysicsList::AnaEx01PhysicsList():  G4VUserPhysicsList()
49{
50  currentDefaultCut = defaultCutValue = 1.0*mm;
51  cutForGamma       = defaultCutValue;
52  cutForElectron    = defaultCutValue;
53  cutForProton      = defaultCutValue;
54
55 SetVerboseLevel(1);
56}
57
58//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
59
60AnaEx01PhysicsList::~AnaEx01PhysicsList()
61{}
62
63//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
64
65void AnaEx01PhysicsList::ConstructParticle()
66{
67  // In this method, static member functions should be called
68  // for all particles which you want to use.
69  // This ensures that objects of these particle types will be
70  // created in the program.
71
72  ConstructBosons();
73  ConstructLeptons();
74  ConstructMesons();
75  ConstructBaryons();
76}
77
78//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
79
80void AnaEx01PhysicsList::ConstructBosons()
81{
82  // pseudo-particles
83  G4Geantino::GeantinoDefinition();
84  G4ChargedGeantino::ChargedGeantinoDefinition();
85
86  // gamma
87  G4Gamma::GammaDefinition();
88
89  // optical photon
90  G4OpticalPhoton::OpticalPhotonDefinition();
91}
92 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
93
94void AnaEx01PhysicsList::ConstructLeptons()
95{
96  // leptons
97  G4Electron::ElectronDefinition();
98  G4Positron::PositronDefinition();
99  G4MuonPlus::MuonPlusDefinition();
100  G4MuonMinus::MuonMinusDefinition();
101
102  G4NeutrinoE::NeutrinoEDefinition();
103  G4AntiNeutrinoE::AntiNeutrinoEDefinition();
104  G4NeutrinoMu::NeutrinoMuDefinition();
105  G4AntiNeutrinoMu::AntiNeutrinoMuDefinition();
106}
107
108//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
109
110void AnaEx01PhysicsList::ConstructMesons()
111{
112 //  mesons
113  G4PionPlus::PionPlusDefinition();
114  G4PionMinus::PionMinusDefinition();
115  G4PionZero::PionZeroDefinition();
116  G4Eta::EtaDefinition();
117  G4EtaPrime::EtaPrimeDefinition();
118  G4KaonPlus::KaonPlusDefinition();
119  G4KaonMinus::KaonMinusDefinition();
120  G4KaonZero::KaonZeroDefinition();
121  G4AntiKaonZero::AntiKaonZeroDefinition();
122  G4KaonZeroLong::KaonZeroLongDefinition();
123  G4KaonZeroShort::KaonZeroShortDefinition();
124}
125
126//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
127
128void AnaEx01PhysicsList::ConstructBaryons()
129{
130//  barions
131  G4Proton::ProtonDefinition();
132  G4AntiProton::AntiProtonDefinition();
133  G4Neutron::NeutronDefinition();
134  G4AntiNeutron::AntiNeutronDefinition();
135}
136
137
138//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
139
140void AnaEx01PhysicsList::ConstructProcess()
141{
142  AddTransportation();
143  ConstructEM();
144  ConstructGeneral();
145}
146
147//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
148
149#include "G4ComptonScattering.hh"
150#include "G4GammaConversion.hh"
151#include "G4PhotoElectricEffect.hh"
152
153#include "G4MultipleScattering.hh"
154
155#include "G4eIonisation.hh"
156#include "G4eBremsstrahlung.hh"
157#include "G4eplusAnnihilation.hh"
158
159#include "G4MuIonisation.hh"
160#include "G4MuBremsstrahlung.hh"
161#include "G4MuPairProduction.hh"
162
163#include "G4hIonisation.hh"
164
165//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
166
167void AnaEx01PhysicsList::ConstructEM()
168{
169  theParticleIterator->reset();
170  while( (*theParticleIterator)() ){
171    G4ParticleDefinition* particle = theParticleIterator->value();
172    G4ProcessManager* pmanager = particle->GetProcessManager();
173    G4String particleName = particle->GetParticleName();
174     
175    if (particleName == "gamma") {
176    // gamma
177      pmanager->AddDiscreteProcess(new G4GammaConversion());
178      pmanager->AddDiscreteProcess(new G4ComptonScattering());     
179      pmanager->AddDiscreteProcess(new G4PhotoElectricEffect());
180
181    } else if (particleName == "e-") {
182    //electron
183      G4VProcess* theeminusMultipleScattering = new G4MultipleScattering();
184      G4VProcess* theeminusIonisation         = new G4eIonisation();
185      G4VProcess* theeminusBremsstrahlung     = new G4eBremsstrahlung();
186      //
187      // add processes
188      pmanager->AddProcess(theeminusMultipleScattering);
189      pmanager->AddProcess(theeminusIonisation);
190      pmanager->AddProcess(theeminusBremsstrahlung);
191      //     
192      // set ordering for AlongStepDoIt
193      pmanager->SetProcessOrdering(theeminusMultipleScattering, idxAlongStep,1);
194      pmanager->SetProcessOrdering(theeminusIonisation,         idxAlongStep,2);
195      //
196      // set ordering for PostStepDoIt
197      pmanager->SetProcessOrdering(theeminusMultipleScattering, idxPostStep,1);
198      pmanager->SetProcessOrdering(theeminusIonisation,         idxPostStep,2);
199      pmanager->SetProcessOrdering(theeminusBremsstrahlung,     idxPostStep,3);
200
201    } else if (particleName == "e+") {
202    //positron
203      G4VProcess* theeplusMultipleScattering = new G4MultipleScattering();
204      G4VProcess* theeplusIonisation         = new G4eIonisation();
205      G4VProcess* theeplusBremsstrahlung     = new G4eBremsstrahlung();
206      G4VProcess* theeplusAnnihilation       = new G4eplusAnnihilation();
207      //
208      // add processes
209      pmanager->AddProcess(theeplusMultipleScattering);
210      pmanager->AddProcess(theeplusIonisation);
211      pmanager->AddProcess(theeplusBremsstrahlung);
212      pmanager->AddProcess(theeplusAnnihilation);
213      //
214      // set ordering for AtRestDoIt
215      pmanager->SetProcessOrderingToFirst(theeplusAnnihilation, idxAtRest);
216      //
217      // set ordering for AlongStepDoIt
218      pmanager->SetProcessOrdering(theeplusMultipleScattering, idxAlongStep,1);
219      pmanager->SetProcessOrdering(theeplusIonisation,         idxAlongStep,2);
220      //
221      // set ordering for PostStepDoIt
222      pmanager->SetProcessOrdering(theeplusMultipleScattering, idxPostStep,1);
223      pmanager->SetProcessOrdering(theeplusIonisation,         idxPostStep,2);
224      pmanager->SetProcessOrdering(theeplusBremsstrahlung,     idxPostStep,3);
225      pmanager->SetProcessOrdering(theeplusAnnihilation,       idxPostStep,4);
226 
227    } else if( particleName == "mu+" || 
228               particleName == "mu-"    ) {
229    //muon 
230      G4VProcess* aMultipleScattering = new G4MultipleScattering();
231      G4VProcess* aBremsstrahlung     = new G4MuBremsstrahlung();
232      G4VProcess* aPairProduction     = new G4MuPairProduction();
233      G4VProcess* anIonisation        = new G4MuIonisation();
234      //
235      // add processes
236      pmanager->AddProcess(anIonisation);
237      pmanager->AddProcess(aMultipleScattering);
238      pmanager->AddProcess(aBremsstrahlung);
239      pmanager->AddProcess(aPairProduction);
240      //
241      // set ordering for AlongStepDoIt
242      pmanager->SetProcessOrdering(aMultipleScattering, idxAlongStep,1);
243      pmanager->SetProcessOrdering(anIonisation,        idxAlongStep,2);
244      //
245      // set ordering for PostStepDoIt
246      pmanager->SetProcessOrdering(aMultipleScattering, idxPostStep,1);
247      pmanager->SetProcessOrdering(anIonisation,        idxPostStep,2);
248      pmanager->SetProcessOrdering(aBremsstrahlung,     idxPostStep,3);
249      pmanager->SetProcessOrdering(aPairProduction,     idxPostStep,4);
250
251     } else if ((!particle->IsShortLived()) &&
252               (particle->GetPDGCharge() != 0.0) && 
253               (particle->GetParticleName() != "chargedgeantino")) {
254     // all others charged particles except geantino     
255     G4VProcess* aMultipleScattering = new G4MultipleScattering();
256     G4VProcess* anIonisation        = new G4hIonisation();
257     //
258     // add processes
259     pmanager->AddProcess(anIonisation);
260     pmanager->AddProcess(aMultipleScattering);
261     //
262     // set ordering for AlongStepDoIt
263     pmanager->SetProcessOrdering(aMultipleScattering, idxAlongStep,1);
264     pmanager->SetProcessOrdering(anIonisation,        idxAlongStep,2);
265     //
266     // set ordering for PostStepDoIt
267     pmanager->SetProcessOrdering(aMultipleScattering, idxPostStep,1);
268     pmanager->SetProcessOrdering(anIonisation,        idxPostStep,2);
269    }
270  }
271}
272
273//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
274
275#include "G4Decay.hh"
276
277void AnaEx01PhysicsList::ConstructGeneral()
278{
279  // Add Decay Process
280   G4Decay* theDecayProcess = new G4Decay();
281  theParticleIterator->reset();
282  while( (*theParticleIterator)() ){
283    G4ParticleDefinition* particle = theParticleIterator->value();
284    G4ProcessManager* pmanager = particle->GetProcessManager();
285    if (theDecayProcess->IsApplicable(*particle)) { 
286      pmanager ->AddProcess(theDecayProcess);
287      // set ordering for PostStepDoIt and AtRestDoIt
288      pmanager ->SetProcessOrdering(theDecayProcess, idxPostStep);
289      pmanager ->SetProcessOrdering(theDecayProcess, idxAtRest);
290    }
291  }
292}
293
294//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
295
296void AnaEx01PhysicsList::SetCuts()
297{
298  // reactualise cutValues
299  if (currentDefaultCut != defaultCutValue)
300    {
301     if(cutForGamma    == currentDefaultCut) cutForGamma    = defaultCutValue;
302     if(cutForElectron == currentDefaultCut) cutForElectron = defaultCutValue;
303     if(cutForProton   == currentDefaultCut) cutForProton   = defaultCutValue;
304     currentDefaultCut = defaultCutValue;
305    }
306   
307  if (verboseLevel >0){
308    G4cout << "AnaEx01PhysicsList::SetCuts:";
309    G4cout << "CutLength : " << G4BestUnit(defaultCutValue,"Length") << G4endl;   
310  } 
311
312  // set cut values for gamma at first and for e- second and next for e+,
313  // because some processes for e+/e- need cut values for gamma
314  SetCutValue(cutForGamma, "gamma");
315  SetCutValue(cutForElectron, "e-");
316  SetCutValue(cutForElectron, "e+");
317 
318  // set cut values for proton and anti_proton before all other hadrons
319  // because some processes for hadrons need cut values for proton/anti_proton
320  SetCutValue(cutForProton, "proton");
321  SetCutValue(cutForProton, "anti_proton");
322 
323  SetCutValueForOthers(defaultCutValue);
324
325  if (verboseLevel>0) DumpCutValuesTable();
326}
327
328//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
329
330
331void AnaEx01PhysicsList::SetCutForGamma(G4double cut)
332{
333  ResetCuts();
334  cutForGamma = cut;
335}
336
337void AnaEx01PhysicsList::SetCutForElectron(G4double cut)
338{
339  ResetCuts();
340  cutForElectron = cut;
341}
342
343void AnaEx01PhysicsList::SetCutForProton(G4double cut)
344{
345  ResetCuts();
346  cutForProton = cut;
347}
348
349G4double AnaEx01PhysicsList::GetCutForGamma() const
350{
351  return cutForGamma;
352}
353
354G4double AnaEx01PhysicsList::GetCutForElectron() const
355{
356  return cutForElectron;
357}
358
359G4double AnaEx01PhysicsList::GetCutForProton() const
360{
361  return cutForGamma;
362}
363
364//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
365
366
367
368
369
370
371
372
373
374
375
Note: See TracBrowser for help on using the repository browser.