| 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 // Author: Mathieu Karamitros (kara (AT) cenbg
28 //
29 // History:
30 // -----------
31 // 10 Oct 2011 M.Karamitros created
32 //
33 // -------------------------------------------
34
35 #include "G4ITReactionChange.hh"
36
37 G4ITReactionChange::G4ITReactionChange()
38 = default;
39
40 G4ITReactionChange::~G4ITReactionChange()
41 {
42 //dtor
43 delete fSecondaries;
44 fSecondaries = nullptr;
45 }
46
47 // Should not be used
48 G4ITReactionChange::G4ITReactionChange(const G
49 {
50 }
51
52 // should not be used
53 G4ITReactionChange& G4ITReactionChange::operat
54 {
55 if (this == &rhs) return *this; // handle se
56 //assignment operator
57 return *this;
58 }
59
60 void G4ITReactionChange::Initialize(const G4Tr
61 const G4Tr
62 G4VParticl
63 G4VParticl
64 {
65 fParticleChange.clear();
66 fParticleChange[&trackA] = particleChangeA;
67 fParticleChange[&trackB] = particleChangeB;
68
69 if ((particleChangeA != nullptr) || (particl
70 {
71 G4bool test = (particleChangeA != nullptr)
72
73 if (!test)
74 {
75 G4ExceptionDescription exceptionDescript
76 exceptionDescription << "If you give for
77 exceptionDescription
78 << "G4ITReactionChange is expecting
79 exceptionDescription << "reacting tracks
80 G4Exception("G4ITReactionChange::Initial
81 FatalErrorInArgument, except
82 }
83
84 fParticleChangeIsSet = true;
85
86 fParticleChange[&trackA]->Initialize(track
87 fParticleChange[&trackB]->Initialize(track
88 ;
89
90 }
91
92 fSecondaries = nullptr;
93 fNumberOfSecondaries = 0;
94 fKillParents = false;
95 }
96
97 void G4ITReactionChange::AddSecondary(G4Track*
98 {
99 if (fSecondaries == nullptr) fSecondaries =
100 fSecondaries->push_back(aTrack);
101 fNumberOfSecondaries++;
102 }
103
104 void G4ITReactionChange::UpdateStepInfo(G4Step
105 {
106 fParticleChange[stepA->GetTrack()]->UpdateSt
107 fParticleChange[stepB->GetTrack()]->UpdateSt
108 }
109
110 G4VParticleChange* G4ITReactionChange::GetPart
111 {
112 auto it = fParticleChange
113 .find(track);
114
115 if (it == fParticleChange.end()) return null
116 return it->second;
117 }
118
119 const G4Track* G4ITReactionChange::GetTrackA()
120 {
121 auto it = fParticleChange.begin();
122 if (it != fParticleChange.end())
123 {
124 return it->first;
125 }
126
127 G4ExceptionDescription exceptionDescription;
128 exceptionDescription
129 << "No track A found ! Have you initiali
130 G4Exception("G4ITReactionChange::GetTrackA",
131 FatalErrorInArgument, exceptionD
132 return nullptr;
133 }
134
135 const G4Track* G4ITReactionChange::GetTrackB()
136 {
137 auto it = fParticleChange.begin();
138 auto next = it++;
139 if (next == fParticleChange.end())
140 {
141 G4ExceptionDescription exceptionDescriptio
142 exceptionDescription
143 << "No track B found ! Have you initia
144 G4Exception("G4ITReactionChange::GetTrackB
145 FatalErrorInArgument, exceptio
146 }
147
148 return it->first;
149 }
150