| 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 * Based on 'G4pbc'.
28 * Copyright (c) 2020 Amentum Pty Ltd
29 * team@amentum.space
30 * The original open-source version of this code
31 * may be found at https://github.com/amentumspace/g4pbc
32 * Permission is hereby granted, free of charge, to any person obtaining a copy of this software and
33 * associated documentation files (the "Software"), to deal in the Software without restriction,
34 * including without limitation the rights to use, copy, modify, merge, publish, distribute,
35 * sublicense, and/or sell copies of the Software, and to permit persons to whom the Software
36 * is furnished to do so, subject to the following conditions:
37 * The above copyright notice and this permission notice shall be included in all copies
38 * or substantial portions of the Software.
39 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING BUT
40 * NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
41 * NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM,
42 * DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
43 * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
44 *
45 *
46 */
47 #include "ParticleChangeForPeriodic.hh"
48
49 #include "G4DynamicParticle.hh"
50 #include "G4Step.hh"
51 #include "G4Track.hh"
52
53 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
54
55 ParticleChangeForPeriodic::ParticleChangeForPeriodic() : G4VParticleChange() {}
56
57 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
58
59 G4Step* ParticleChangeForPeriodic::UpdateStepForPostStep(G4Step* pStep)
60 {
61 G4StepPoint* pPostStepPoint = pStep->GetPostStepPoint();
62 pPostStepPoint->SetMomentumDirection(fProposedMomentumDirection);
63 pPostStepPoint->SetPolarization(fProposedPolarization);
64 pPostStepPoint->SetPosition(fProposedPosition);
65
66 if (isParentWeightProposed) {
67 pPostStepPoint->SetWeight(theParentWeight);
68 }
69
70 pStep->AddTotalEnergyDeposit(theLocalEnergyDeposit);
71 pStep->AddNonIonizingEnergyDeposit(theNonIonizingEnergyDeposit);
72
73 return pStep;
74 }
75
76 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
77
78 void ParticleChangeForPeriodic::AddSecondary(G4DynamicParticle* aParticle)
79 {
80 auto* aTrack = new G4Track(aParticle, fTrack->GetGlobalTime(), fTrack->GetPosition());
81
82 aTrack->SetTouchableHandle(fTrack->GetTouchableHandle());
83
84 G4VParticleChange::AddSecondary(aTrack);
85 }
86
87 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
88
89 void ParticleChangeForPeriodic::DumpInfo() const
90 {
91 G4VParticleChange::DumpInfo();
92 G4int oldprc = G4cout.precision(3);
93
94 G4cout << " Momentum Direction: " << std::setw(20) << fProposedMomentumDirection << G4endl;
95 G4cout << " Polarization: " << std::setw(20) << fProposedPolarization << G4endl;
96 G4cout << " Position: " << std::setw(20) << fProposedPosition << G4endl;
97 G4cout.precision(oldprc);
98 }
99
100 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
101