| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26
27 #include "PeriodicBoundaryPhysics.hh"
28
29 #include "PeriodicBoundaryProcess.hh"
30
31 #include "G4PhysicsConstructorFactory.hh"
32 #include "G4ProcessManager.hh"
33 #include "globals.hh"
34 //....oooOO0OOooo........oooOO0OOooo........oo
35
36 G4_DECLARE_PHYSCONSTR_FACTORY(PeriodicBoundary
37 //....oooOO0OOooo........oooOO0OOooo........oo
38
39 PeriodicBoundaryPhysics::PeriodicBoundaryPhysi
40
41 : G4VPhysicsConstructor(name), fPeriodicX(pe
42 {
43 verboseLevel = 0;
44 }
45
46 //....oooOO0OOooo........oooOO0OOooo........oo
47
48 void PeriodicBoundaryPhysics::ConstructParticl
49
50 //....oooOO0OOooo........oooOO0OOooo........oo
51
52 void PeriodicBoundaryPhysics::ConstructProcess
53 {
54 if (verboseLevel > 0) G4cout << "Constructin
55
56 auto* pbc =
57 new PeriodicBoundaryProcess("Cyclic", fNot
58
59 if (verboseLevel > 0) {
60 pbc->SetVerboseLevel(verboseLevel);
61 }
62
63 auto aParticleIterator = GetParticleIterator
64
65 aParticleIterator->reset();
66
67 G4ProcessManager* processManager = nullptr;
68
69 while ((*aParticleIterator)()) {
70 G4ParticleDefinition* particle = aParticle
71
72 G4String particleName = particle->GetParti
73
74 processManager = particle->GetProcessManag
75
76 if (!processManager) {
77 ThrowException(particleName);
78 return;
79 }
80
81 AddDiscreteProcess(pbc, *particle, process
82 }
83 }
84
85 //....oooOO0OOooo........oooOO0OOooo........oo
86
87 void PeriodicBoundaryPhysics::ThrowException(c
88 {
89 std::ostringstream o;
90 o << "Particle " << particleName << "without
91 G4Exception("G4PeriodicBoundaryPhysics::Cons
92 }
93
94 //....oooOO0OOooo........oooOO0OOooo........oo
95
96 void PeriodicBoundaryPhysics::AddDiscreteProce
97
98
99 {
100 if (periodicBoundaryProcess->IsApplicable(pa
101 processManager->AddDiscreteProcess(periodi
102 }
103 }
104
105 //....oooOO0OOooo........oooOO0OOooo........oo
106