70 fpAuxiliaryPointVector(0),
76 fPreStepPointGlobalTime(0),
77 fPostStepPointGlobalTime(0),
78 fPreStepPointWeight(1.),
79 fPostStepPointWeight(1.)
84 fpAuxiliaryPointVector(0),
86 fRemainingEnergy(aTrack->GetKineticEnergy()),
90 fPreStepPointGlobalTime(aTrack->GetGlobalTime()),
91 fPostStepPointGlobalTime(aTrack->GetGlobalTime()),
92 fpPreStepPointVolume(aTrack->GetTouchableHandle()),
93 fpPostStepPointVolume(aTrack->GetNextTouchableHandle()),
94 fPreStepPointWeight(aTrack->GetWeight()),
95 fPostStepPointWeight(aTrack->GetWeight())
100 fpAuxiliaryPointVector(aStep->GetPointerToVectorOfAuxiliaryPoints()),
101 fTotEDep(aStep->GetTotalEnergyDeposit())
150 const std::map<G4String,G4AttDef>*
154 std::map<G4String,G4AttDef>* store
164 (*store)[ID] =
G4AttDef(ID,
"Auxiliary Point Position",
165 "Physics",
"G4BestUnit",
"G4ThreeVector");
167 (*store)[ID] =
G4AttDef(ID,
"Total Energy Deposit",
168 "Physics",
"G4BestUnit",
"G4double");
170 (*store)[ID] =
G4AttDef(ID,
"Remaining Energy",
171 "Physics",
"G4BestUnit",
"G4double");
173 (*store)[ID] =
G4AttDef(ID,
"Process Defined Step",
174 "Physics",
"",
"G4String");
176 (*store)[ID] =
G4AttDef(ID,
"Process Type Defined Step",
177 "Physics",
"",
"G4String");
179 (*store)[ID] =
G4AttDef(ID,
"Pre-step-point status",
180 "Physics",
"",
"G4String");
182 (*store)[ID] =
G4AttDef(ID,
"Post-step-point status",
183 "Physics",
"",
"G4String");
185 (*store)[ID] =
G4AttDef(ID,
"Pre-step-point global time",
186 "Physics",
"G4BestUnit",
"G4double");
188 (*store)[ID] =
G4AttDef(ID,
"Post-step-point global time",
189 "Physics",
"G4BestUnit",
"G4double");
191 (*store)[ID] =
G4AttDef(ID,
"Pre-step Volume Path",
192 "Physics",
"",
"G4String");
194 (*store)[ID] =
G4AttDef(ID,
"Post-step Volume Path",
195 "Physics",
"",
"G4String");
197 (*store)[ID] =
G4AttDef(ID,
"Pre-step-point weight",
198 "Physics",
"",
"G4double");
200 (*store)[ID] =
G4AttDef(ID,
"Post-step-point weight",
201 "Physics",
"",
"G4double");
217 case fUndefined: status =
"fUndefined";
break;
218 default: status =
"Not recognised";
break;
225 std::ostringstream oss;
226 G4int depth = th->GetHistoryDepth();
227 for (
G4int i = depth; i >= 0; --i) {
228 oss << th->GetVolume(i)->GetName()
229 <<
':' << th->GetCopyNumber(i);
230 if (i != 0) oss <<
'/';
241 std::vector<G4ThreeVector>::iterator iAux;
260 values->push_back(
G4AttValue(
"PDS",
"None",
""));
261 values->push_back(
G4AttValue(
"PTDS",
"None",
""));
279 values->push_back(
G4AttValue(
"PreVPath",
"None",
""));
285 values->push_back(
G4AttValue(
"PostVPath",
"None",
""));
289 std::ostringstream oss;
296 std::ostringstream oss;