Geant4 Cross Reference |
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 #include "PeriodicBoundaryPhysics.hh" 28 29 #include "PeriodicBoundaryProcess.hh" 30 31 #include "G4PhysicsConstructorFactory.hh" 32 #include "G4ProcessManager.hh" 33 #include "globals.hh" 34 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 35 36 G4_DECLARE_PHYSCONSTR_FACTORY(PeriodicBoundaryPhysics); 37 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 38 39 PeriodicBoundaryPhysics::PeriodicBoundaryPhysics(const G4String& name, G4bool per_x, G4bool per_y, 40 G4bool per_z) 41 : G4VPhysicsConstructor(name), fPeriodicX(per_x), fPeriodicY(per_y), fPeriodicZ(per_z) 42 { 43 verboseLevel = 0; 44 } 45 46 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 47 48 void PeriodicBoundaryPhysics::ConstructParticle() {} 49 50 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 51 52 void PeriodicBoundaryPhysics::ConstructProcess() 53 { 54 if (verboseLevel > 0) G4cout << "Constructing cyclic boundary physics process" << G4endl; 55 56 auto* pbc = 57 new PeriodicBoundaryProcess("Cyclic", fNotDefined, fPeriodicX, fPeriodicY, fPeriodicZ); 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 = aParticleIterator->value(); 71 72 G4String particleName = particle->GetParticleName(); 73 74 processManager = particle->GetProcessManager(); 75 76 if (!processManager) { 77 ThrowException(particleName); 78 return; 79 } 80 81 AddDiscreteProcess(pbc, *particle, processManager); 82 } 83 } 84 85 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 86 87 void PeriodicBoundaryPhysics::ThrowException(const G4String& particleName) 88 { 89 std::ostringstream o; 90 o << "Particle " << particleName << "without a Process Manager"; 91 G4Exception("G4PeriodicBoundaryPhysics::ConstructProcess()", "", FatalException, o.str().c_str()); 92 } 93 94 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 95 96 void PeriodicBoundaryPhysics::AddDiscreteProcess(PeriodicBoundaryProcess* periodicBoundaryProcess, 97 G4ParticleDefinition& particle, 98 G4ProcessManager* processManager) 99 { 100 if (periodicBoundaryProcess->IsApplicable(particle)) { 101 processManager->AddDiscreteProcess(periodicBoundaryProcess); 102 } 103 } 104 105 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 106