| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 // 26 //
27 // Author: Mathieu Karamitros (kara (AT) cenbg 27 // Author: Mathieu Karamitros (kara (AT) cenbg . in2p3 . fr)
28 // 28 //
29 // History: 29 // History:
30 // ----------- 30 // -----------
31 // 10 Oct 2011 M.Karamitros created 31 // 10 Oct 2011 M.Karamitros created
32 // 32 //
33 // ------------------------------------------- 33 // -------------------------------------------------------------------
34 34
35 #include "G4ITReactionChange.hh" 35 #include "G4ITReactionChange.hh"
36 36
37 G4ITReactionChange::G4ITReactionChange() 37 G4ITReactionChange::G4ITReactionChange()
38 = default; 38 = default;
39 39
40 G4ITReactionChange::~G4ITReactionChange() 40 G4ITReactionChange::~G4ITReactionChange()
41 { 41 {
42 //dtor 42 //dtor
43 delete fSecondaries; 43 delete fSecondaries;
44 fSecondaries = nullptr; 44 fSecondaries = nullptr;
45 } 45 }
46 46
47 // Should not be used 47 // Should not be used
48 G4ITReactionChange::G4ITReactionChange(const G 48 G4ITReactionChange::G4ITReactionChange(const G4ITReactionChange& /*other*/)
49 { 49 {
50 } 50 }
51 51
52 // should not be used 52 // should not be used
53 G4ITReactionChange& G4ITReactionChange::operat 53 G4ITReactionChange& G4ITReactionChange::operator=(const G4ITReactionChange& rhs)
54 { 54 {
55 if (this == &rhs) return *this; // handle se 55 if (this == &rhs) return *this; // handle self assignment
56 //assignment operator 56 //assignment operator
57 return *this; 57 return *this;
58 } 58 }
59 59
60 void G4ITReactionChange::Initialize(const G4Tr 60 void G4ITReactionChange::Initialize(const G4Track& trackA,
61 const G4Tr 61 const G4Track& trackB,
62 G4VParticl 62 G4VParticleChange* particleChangeA,
63 G4VParticl 63 G4VParticleChange* particleChangeB)
64 { 64 {
65 fParticleChange.clear(); 65 fParticleChange.clear();
66 fParticleChange[&trackA] = particleChangeA; 66 fParticleChange[&trackA] = particleChangeA;
67 fParticleChange[&trackB] = particleChangeB; 67 fParticleChange[&trackB] = particleChangeB;
68 68
69 if ((particleChangeA != nullptr) || (particl 69 if ((particleChangeA != nullptr) || (particleChangeB != nullptr))
70 { 70 {
71 G4bool test = (particleChangeA != nullptr) 71 G4bool test = (particleChangeA != nullptr) && (particleChangeB != nullptr);
72 72
73 if (!test) 73 if (!test)
74 { 74 {
75 G4ExceptionDescription exceptionDescript 75 G4ExceptionDescription exceptionDescription;
76 exceptionDescription << "If you give for 76 exceptionDescription << "If you give for one track a particleChange, ";
77 exceptionDescription 77 exceptionDescription
78 << "G4ITReactionChange is expecting 78 << "G4ITReactionChange is expecting that you give for both ";
79 exceptionDescription << "reacting tracks 79 exceptionDescription << "reacting tracks a particleChange.";
80 G4Exception("G4ITReactionChange::Initial 80 G4Exception("G4ITReactionChange::Initialize", "ITReactionChange001",
81 FatalErrorInArgument, except 81 FatalErrorInArgument, exceptionDescription);
82 } 82 }
83 83
84 fParticleChangeIsSet = true; 84 fParticleChangeIsSet = true;
85 85
86 fParticleChange[&trackA]->Initialize(track 86 fParticleChange[&trackA]->Initialize(trackA);
87 fParticleChange[&trackB]->Initialize(track 87 fParticleChange[&trackB]->Initialize(trackB);
88 ; 88 ;
89 89
90 } 90 }
91 91
92 fSecondaries = nullptr; 92 fSecondaries = nullptr;
93 fNumberOfSecondaries = 0; 93 fNumberOfSecondaries = 0;
94 fKillParents = false; 94 fKillParents = false;
95 } 95 }
96 96
97 void G4ITReactionChange::AddSecondary(G4Track* 97 void G4ITReactionChange::AddSecondary(G4Track* aTrack)
98 { 98 {
99 if (fSecondaries == nullptr) fSecondaries = 99 if (fSecondaries == nullptr) fSecondaries = new std::vector<G4Track*>();
100 fSecondaries->push_back(aTrack); 100 fSecondaries->push_back(aTrack);
101 fNumberOfSecondaries++; 101 fNumberOfSecondaries++;
102 } 102 }
103 103
104 void G4ITReactionChange::UpdateStepInfo(G4Step 104 void G4ITReactionChange::UpdateStepInfo(G4Step* stepA, G4Step* stepB)
105 { 105 {
106 fParticleChange[stepA->GetTrack()]->UpdateSt 106 fParticleChange[stepA->GetTrack()]->UpdateStepForPostStep(stepA);
107 fParticleChange[stepB->GetTrack()]->UpdateSt 107 fParticleChange[stepB->GetTrack()]->UpdateStepForPostStep(stepB);
108 } 108 }
109 109
110 G4VParticleChange* G4ITReactionChange::GetPart 110 G4VParticleChange* G4ITReactionChange::GetParticleChange(const G4Track* track)
111 { 111 {
112 auto it = fParticleChange 112 auto it = fParticleChange
113 .find(track); 113 .find(track);
114 114
115 if (it == fParticleChange.end()) return null 115 if (it == fParticleChange.end()) return nullptr;
116 return it->second; 116 return it->second;
117 } 117 }
118 118
119 const G4Track* G4ITReactionChange::GetTrackA() 119 const G4Track* G4ITReactionChange::GetTrackA()
120 { 120 {
121 auto it = fParticleChange.begin(); 121 auto it = fParticleChange.begin();
122 if (it != fParticleChange.end()) 122 if (it != fParticleChange.end())
123 { 123 {
124 return it->first; 124 return it->first;
125 } 125 }
126 126
127 G4ExceptionDescription exceptionDescription; 127 G4ExceptionDescription exceptionDescription;
128 exceptionDescription 128 exceptionDescription
129 << "No track A found ! Have you initiali 129 << "No track A found ! Have you initialized the ReactionChange ?";
130 G4Exception("G4ITReactionChange::GetTrackA", 130 G4Exception("G4ITReactionChange::GetTrackA", "ITReactionChange001",
131 FatalErrorInArgument, exceptionD 131 FatalErrorInArgument, exceptionDescription);
132 return nullptr; 132 return nullptr;
133 } 133 }
134 134
135 const G4Track* G4ITReactionChange::GetTrackB() 135 const G4Track* G4ITReactionChange::GetTrackB()
136 { 136 {
137 auto it = fParticleChange.begin(); 137 auto it = fParticleChange.begin();
138 auto next = it++; 138 auto next = it++;
139 if (next == fParticleChange.end()) 139 if (next == fParticleChange.end())
140 { 140 {
141 G4ExceptionDescription exceptionDescriptio 141 G4ExceptionDescription exceptionDescription;
142 exceptionDescription 142 exceptionDescription
143 << "No track B found ! Have you initia 143 << "No track B found ! Have you initialized the ReactionChange ?";
144 G4Exception("G4ITReactionChange::GetTrackB 144 G4Exception("G4ITReactionChange::GetTrackB", "ITReactionChange002",
145 FatalErrorInArgument, exceptio 145 FatalErrorInArgument, exceptionDescription);
146 } 146 }
147 147
148 return it->first; 148 return it->first;
149 } 149 }
150 150