| 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 // $Id: G4ITReactionChange.cc 64057 2012-10-30 15:04:49Z gcosmo $
26 // 27 //
27 // Author: Mathieu Karamitros (kara (AT) cenbg 28 // Author: Mathieu Karamitros (kara (AT) cenbg . in2p3 . fr)
28 // 29 //
29 // History: 30 // History:
30 // ----------- 31 // -----------
31 // 10 Oct 2011 M.Karamitros created 32 // 10 Oct 2011 M.Karamitros created
32 // 33 //
33 // ------------------------------------------- 34 // -------------------------------------------------------------------
34 35
35 #include "G4ITReactionChange.hh" 36 #include "G4ITReactionChange.hh"
36 37
37 G4ITReactionChange::G4ITReactionChange() << 38 G4ITReactionChange::G4ITReactionChange():
38 = default; << 39 fSecondaries(0),
>> 40 fNumberOfSecondaries(0),
>> 41 fKillParents(false),
>> 42 fParticleChangeIsSet(false)
>> 43 {
>> 44 //ctor
>> 45 }
39 46
40 G4ITReactionChange::~G4ITReactionChange() 47 G4ITReactionChange::~G4ITReactionChange()
41 { 48 {
42 //dtor << 49 //dtor
43 delete fSecondaries; << 50 delete fSecondaries;
44 fSecondaries = nullptr; << 51 fSecondaries = 0 ;
45 } 52 }
46 53
47 // Should not be used 54 // Should not be used
48 G4ITReactionChange::G4ITReactionChange(const G << 55 G4ITReactionChange::G4ITReactionChange(const G4ITReactionChange& /*other*/) :
>> 56 fSecondaries(0),
>> 57 fNumberOfSecondaries(0),
>> 58 fKillParents(false),
>> 59 fParticleChangeIsSet(false)
49 { 60 {
>> 61 //copy ctor
50 } 62 }
51 63
52 // should not be used 64 // should not be used
53 G4ITReactionChange& G4ITReactionChange::operat 65 G4ITReactionChange& G4ITReactionChange::operator=(const G4ITReactionChange& rhs)
54 { 66 {
55 if (this == &rhs) return *this; // handle se << 67 if (this == &rhs) return *this; // handle self assignment
56 //assignment operator << 68 //assignment operator
57 return *this; << 69 return *this;
58 } 70 }
59 71
60 void G4ITReactionChange::Initialize(const G4Tr 72 void G4ITReactionChange::Initialize(const G4Track& trackA,
61 const G4Tr << 73 const G4Track& trackB,
62 G4VParticl << 74 G4VParticleChange* particleChangeA,
63 G4VParticl << 75 G4VParticleChange* particleChangeB
64 { << 76 )
65 fParticleChange.clear(); << 77 {
66 fParticleChange[&trackA] = particleChangeA; << 78 fParticleChange.clear();
67 fParticleChange[&trackB] = particleChangeB; << 79 fParticleChange[&trackA] = particleChangeA;
68 << 80 fParticleChange[&trackB] = particleChangeB;
69 if ((particleChangeA != nullptr) || (particl <<
70 { <<
71 G4bool test = (particleChangeA != nullptr) <<
72 81
73 if (!test) << 82 if(particleChangeA || particleChangeB)
74 { 83 {
75 G4ExceptionDescription exceptionDescript << 84 G4bool test = particleChangeA && particleChangeB ;
76 exceptionDescription << "If you give for << 85
77 exceptionDescription << 86 if(test == false)
78 << "G4ITReactionChange is expecting << 87 {
79 exceptionDescription << "reacting tracks << 88 G4ExceptionDescription exceptionDescription ;
80 G4Exception("G4ITReactionChange::Initial << 89 exceptionDescription << "If you give for one track a particleChange, ";
81 FatalErrorInArgument, except << 90 exceptionDescription << "G4ITReactionChange is expecting that you give for both ";
82 } << 91 exceptionDescription << "reacting tracks a particleChange.";
>> 92 G4Exception("G4ITReactionChange::Initialize","ITReactionChange001",
>> 93 FatalErrorInArgument,exceptionDescription);
>> 94 }
83 95
84 fParticleChangeIsSet = true; << 96 fParticleChangeIsSet = true;
85 97
86 fParticleChange[&trackA]->Initialize(track << 98 fParticleChange[&trackA]->Initialize(trackA);
87 fParticleChange[&trackB]->Initialize(track << 99 fParticleChange[&trackB]->Initialize(trackB);;
88 ; <<
89 100
90 } << 101 }
91 102
92 fSecondaries = nullptr; << 103 fSecondaries = 0;
93 fNumberOfSecondaries = 0; << 104 fNumberOfSecondaries = 0;
94 fKillParents = false; << 105 fKillParents = false;
95 } 106 }
96 107
97 void G4ITReactionChange::AddSecondary(G4Track* 108 void G4ITReactionChange::AddSecondary(G4Track* aTrack)
98 { 109 {
99 if (fSecondaries == nullptr) fSecondaries = << 110 if(fSecondaries==0) fSecondaries = new G4TrackFastVector();
100 fSecondaries->push_back(aTrack); << 111
101 fNumberOfSecondaries++; << 112 // add a secondary after size check
>> 113 if (G4TrackFastVectorSize > fNumberOfSecondaries)
>> 114 {
>> 115 fSecondaries->SetElement(fNumberOfSecondaries, aTrack);
>> 116 fNumberOfSecondaries++;
>> 117 }
>> 118 else
>> 119 {
>> 120 G4cerr << "G4ITReactionChange::AddSecondary() Warning ";
>> 121 G4cerr << "fSecondaries is full !! " << G4endl;
>> 122 G4cerr << " The object will not be added in fSecondaries" << G4endl;
>> 123 }
102 } 124 }
103 125
104 void G4ITReactionChange::UpdateStepInfo(G4Step 126 void G4ITReactionChange::UpdateStepInfo(G4Step* stepA, G4Step* stepB)
105 { 127 {
106 fParticleChange[stepA->GetTrack()]->UpdateSt << 128 fParticleChange[stepA->GetTrack()]->UpdateStepForPostStep(stepA);
107 fParticleChange[stepB->GetTrack()]->UpdateSt << 129 fParticleChange[stepB->GetTrack()]->UpdateStepForPostStep(stepB);
108 } 130 }
109 131
110 G4VParticleChange* G4ITReactionChange::GetPart 132 G4VParticleChange* G4ITReactionChange::GetParticleChange(const G4Track* track)
111 { 133 {
112 auto it = fParticleChange << 134 std::map<const G4Track*, G4VParticleChange*>::iterator it = fParticleChange.find(track);
113 .find(track); <<
114 135
115 if (it == fParticleChange.end()) return null << 136 if(it == fParticleChange.end()) return 0;
116 return it->second; << 137 else return it ->second;
117 } 138 }
118 139
119 const G4Track* G4ITReactionChange::GetTrackA() 140 const G4Track* G4ITReactionChange::GetTrackA()
120 { 141 {
121 auto it = fParticleChange.begin(); << 142 std::map<const G4Track*, G4VParticleChange*>::iterator it = fParticleChange.begin();
122 if (it != fParticleChange.end()) << 143 if(it != fParticleChange.end())
123 { << 144 {
124 return it->first; << 145 return it->first;
125 } << 146 }
126 147
127 G4ExceptionDescription exceptionDescription; << 148 G4ExceptionDescription exceptionDescription ;
128 exceptionDescription << 149 exceptionDescription << "No track A found ! Have you initialized the ReactionChange ?";
129 << "No track A found ! Have you initiali << 150 G4Exception("G4ITReactionChange::GetTrackA","ITReactionChange001",
130 G4Exception("G4ITReactionChange::GetTrackA", << 151 FatalErrorInArgument,exceptionDescription);
131 FatalErrorInArgument, exceptionD << 152 return 0;
132 return nullptr; <<
133 } 153 }
134 154
135 const G4Track* G4ITReactionChange::GetTrackB() 155 const G4Track* G4ITReactionChange::GetTrackB()
136 { 156 {
137 auto it = fParticleChange.begin(); << 157 std::map<const G4Track*, G4VParticleChange*>::iterator it = fParticleChange.begin();
138 auto next = it++; << 158 std::map<const G4Track*, G4VParticleChange*>::iterator next = it++;
139 if (next == fParticleChange.end()) << 159 if(next == fParticleChange.end())
140 { << 160 {
141 G4ExceptionDescription exceptionDescriptio << 161 G4ExceptionDescription exceptionDescription ;
142 exceptionDescription << 162 exceptionDescription << "No track B found ! Have you initialized the ReactionChange ?";
143 << "No track B found ! Have you initia << 163 G4Exception("G4ITReactionChange::GetTrackB","ITReactionChange002",
144 G4Exception("G4ITReactionChange::GetTrackB << 164 FatalErrorInArgument,exceptionDescription);
145 FatalErrorInArgument, exceptio << 165 }
146 } <<
147 166
148 return it->first; << 167 return it->first;
149 } 168 }
150 169