| 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 // G4Step class implementation
27 //
28 // Authors:
29 // Katsuya Amako (e-mail: Katsuya.Amako@kek.jp)
30 // Takashi Sasaki (e-mail: Takashi.Sasaki@kek.jp)
31 // Revisions:
32 // Hisaya Kurashige, 1998-2007
33 // --------------------------------------------------------------------
34
35 #include "G4Step.hh"
36
37 // --------------------------------------------------------------------
38 G4Step::G4Step()
39 {
40 fpPreStepPoint = new G4StepPoint();
41 fpPostStepPoint = new G4StepPoint();
42
43 secondaryInCurrentStep = new std::vector<const G4Track*>;
44 }
45
46 // --------------------------------------------------------------------
47 G4Step::~G4Step()
48 {
49 delete fpPreStepPoint;
50 delete fpPostStepPoint;
51
52 delete secondaryInCurrentStep;
53 delete fSecondary;
54 }
55
56 // --------------------------------------------------------------------
57 G4Step::G4Step(const G4Step& right)
58 : fTotalEnergyDeposit(right.fTotalEnergyDeposit)
59 , fNonIonizingEnergyDeposit(right.fNonIonizingEnergyDeposit)
60 , fStepLength(right.fStepLength)
61 , fpTrack(right.fpTrack)
62 , fpSteppingControlFlag(right.fpSteppingControlFlag)
63 , fFirstStepInVolume(right.fFirstStepInVolume)
64 , fLastStepInVolume(right.fLastStepInVolume)
65 , nSecondaryByLastStep(right.nSecondaryByLastStep)
66 , secondaryInCurrentStep(right.secondaryInCurrentStep)
67 , fpVectorOfAuxiliaryPointsPointer(right.fpVectorOfAuxiliaryPointsPointer)
68 {
69 if(right.fpPreStepPoint != nullptr)
70 {
71 fpPreStepPoint = new G4StepPoint(*(right.fpPreStepPoint));
72 }
73 else
74 {
75 fpPreStepPoint = new G4StepPoint();
76 }
77 if(right.fpPostStepPoint != nullptr)
78 {
79 fpPostStepPoint = new G4StepPoint(*(right.fpPostStepPoint));
80 }
81 else
82 {
83 fpPostStepPoint = new G4StepPoint();
84 }
85
86 if(right.fSecondary != nullptr)
87 {
88 fSecondary = new G4TrackVector(*(right.fSecondary));
89 }
90 else
91 {
92 fSecondary = new G4TrackVector();
93 }
94
95 // secondaryInCurrentStep is cleared
96 secondaryInCurrentStep = new std::vector<const G4Track*>;
97 }
98
99 // --------------------------------------------------------------------
100 G4Step& G4Step::operator=(const G4Step& right)
101 {
102 if(this != &right)
103 {
104 fTotalEnergyDeposit = right.fTotalEnergyDeposit;
105 fNonIonizingEnergyDeposit = right.fNonIonizingEnergyDeposit;
106 fStepLength = right.fStepLength;
107 fpTrack = right.fpTrack;
108 fpSteppingControlFlag = right.fpSteppingControlFlag;
109 fFirstStepInVolume = right.fFirstStepInVolume;
110 fLastStepInVolume = right.fLastStepInVolume;
111 nSecondaryByLastStep = right.nSecondaryByLastStep;
112 secondaryInCurrentStep = right.secondaryInCurrentStep;
113 fpVectorOfAuxiliaryPointsPointer = right.fpVectorOfAuxiliaryPointsPointer;
114
115 delete fpPreStepPoint;
116
117 if(right.fpPreStepPoint != nullptr)
118 {
119 fpPreStepPoint = new G4StepPoint(*(right.fpPreStepPoint));
120 }
121 else
122 {
123 fpPreStepPoint = new G4StepPoint();
124 }
125
126 delete fpPostStepPoint;
127
128 if(right.fpPostStepPoint != nullptr)
129 {
130 fpPostStepPoint = new G4StepPoint(*(right.fpPostStepPoint));
131 }
132 else
133 {
134 fpPostStepPoint = new G4StepPoint();
135 }
136
137 if(fSecondary != nullptr)
138 {
139 fSecondary->clear();
140 delete fSecondary;
141 }
142 if(right.fSecondary != nullptr)
143 {
144 fSecondary = new G4TrackVector(*(right.fSecondary));
145 }
146 else
147 {
148 fSecondary = new G4TrackVector();
149 }
150
151 // secondaryInCurrentStep is not copied
152 if(secondaryInCurrentStep != nullptr)
153 {
154 secondaryInCurrentStep->clear();
155 delete secondaryInCurrentStep;
156 }
157 secondaryInCurrentStep = new std::vector<const G4Track*>;
158 }
159 return *this;
160 }
161
162 // --------------------------------------------------------------------
163 const std::vector<const G4Track*>* G4Step::GetSecondaryInCurrentStep() const
164 {
165 secondaryInCurrentStep->clear();
166 std::size_t nSecondary = fSecondary->size();
167 for(std::size_t i = nSecondaryByLastStep; i < nSecondary; ++i)
168 {
169 secondaryInCurrentStep->push_back((*fSecondary)[i]);
170 }
171 return secondaryInCurrentStep;
172 }
173