34 #define INCLXX_IN_GEANT4_MODE 1
48 loadedStoppingTime(0.),
105 for(
PAIter i=iterPair.first, last=iterPair.second; i!=last; ++i) {
111 INCL_WARN(
"Loop in Store::disconnectAvatarFromParticle fell through." << std::endl
112 <<
"This indicates an inconsistent state of the particleAvatarConnections map." << std::endl);
118 for(
ParticleIter particleIter = particlesRelatedToAvatar.begin(),
e = particlesRelatedToAvatar.end();
119 particleIter !=
e; ++particleIter) {
129 for(
PAIter i=iterPair.first, last=iterPair.second; i!=last; ++i) {
145 #ifdef INCL_AVATAR_SEARCH_FullSort
154 #elif defined(INCL_AVATAR_SEARCH_MinElement)
161 #error Unrecognized INCL_AVATAR_SEARCH. Allowed values are: FullSort, MinElement.
170 particleIter != particleEnd; ++particleIter) {
171 (*particleIter)->propagate(step);
210 INCL_WARN(
"Incoming list is not empty when Store::clear() is called" <<
'\n');
225 if((*iter)->isCluster()) {
228 #ifdef INCLXX_IN_GEANT4_MODE
241 G4int projectileA, projectileZ,
A,
Z;
243 G4int ID, type, isParticipant;
247 std::ifstream
in(filename.c_str());
248 in >> projectileA >> projectileZ >> A >> Z >> stoppingTime >> cutNN;
256 in >> ID >> type >> isParticipant >> x >> y >> z >> px >> py >> pz >> E >>
v;
257 if(!
in.good())
break;
264 else if(type == -1) {
269 INCL_FATAL(
"Unrecognized particle type while loading particles; type=" << type <<
'\n');
276 if(isParticipant == 1) {
287 std::stringstream ss;
290 if((*i)->getType() ==
Proton) {
294 if((*i)->getType() ==
Neutron) {
300 ss <<
"0 0 " << A <<
" " <<
Z <<
" "
305 G4int ID = (*i)->getID();
307 if((*i)->getType() ==
Proton) {
310 if((*i)->getType() ==
Neutron) {
314 G4int isParticipant = 0;
315 if((*i)->isParticipant()) {
323 G4double px = (*i)->getMomentum().getX();
324 G4double py = (*i)->getMomentum().getY();
325 G4double pz = (*i)->getMomentum().getZ();
326 G4double V = (*i)->getPotentialEnergy();
328 ss << ID <<
" " << type <<
" " << isParticipant <<
" "
329 << x <<
" " << y <<
" " << z <<
" "
330 << px <<
" " << py <<
" " << pz <<
" "
331 << E <<
" " << V <<
'\n';
338 std::ofstream out(filename.c_str());
344 std::stringstream ss;
346 ss << (*i)->toString() <<
'\n';