G4INCL::Store Class Reference

#include <G4INCLStore.hh>

Public Member Functions
	Store (Config const *const config)
	~Store ()
void	add (Particle *p)
void	addParticleEntryAvatar (IAvatar *a)
	Add one ParticleEntry avatar.
void	addParticleEntryAvatars (IAvatarList const &al)
	Add one ParticleEntry avatar.
void	add (IAvatar *a)
IAvatarList const &	getAvatars () const
void	addIncomingParticle (Particle *const p)
void	removeFromIncoming (Particle *const p)
void	clearIncoming ()
	Clear the incoming list.
void	deleteIncoming ()
	Clear the incoming list and delete the particles.
void	particleHasBeenUpdated (long)
IAvatar *	findSmallestTime ()
void	timeStep (G4double step)
void	particleHasBeenEjected (long)
void	addToOutgoing (Particle *p)
	add the particle to the outgoing particle list.
void	addToOutgoing (ParticleList const &pl)
	Add a list of particles to the outgoing particle list.
void	particleHasBeenDestroyed (long)
void	particleHasEntered (Particle *const particle)
	Move a particle from incoming to inside.
ParticleList const &	getIncomingParticles () const
ParticleList const &	getOutgoingParticles () const
ParticleList	extractDynamicalSpectators ()
	Returns a list of dynamical spectators.
ParticleList const &	getParticles () const
Book *	getBook ()
G4int	countCascading ()
Config const *	getConfig ()
ParticleList	getParticipants ()
ParticleList	getSpectators ()
void	clear ()
void	clearOutgoing ()
void	clearAvatars ()
void	initialiseParticleAvatarConnections ()
	Initialise the particleAvatarConnections map.
void	loadParticles (std::string filename)
G4int	getLoadedA ()
G4int	getLoadedZ ()
G4double	getLoadedStoppingTime ()
std::string	printParticleConfiguration ()
void	writeParticles (std::string filename)
std::string	printAvatars ()
G4bool	containsCollisions () const
Static Public Member Functions
static G4bool	avatarComparisonPredicate (IAvatar *lhs, IAvatar *rhs)
	Comparison predicate for avatars.

---

Detailed Description

The purpose of the Store object is to act as a "particle manager" that keeps track ofall the particles in our simulation. It also tracks the avatars and their connections to particles.

Definition at line 62 of file G4INCLStore.hh.

---

Constructor & Destructor Documentation

G4INCL::Store::Store

(

Config const *const

config

)

Store constructor

Definition at line 44 of file G4INCLStore.cc.

                                          :
    theBook(new Book),
    loadedA(0),
    loadedZ(0),
    loadedStoppingTime(0.),
    theConfig(config)
  {
  }

G4INCL::Store::~Store

(

)

Store destructor

Definition at line 53 of file G4INCLStore.cc.

References clear(), and G4INCL::Book::reset().

                {
    theBook->reset();
    delete theBook;
    theBook = 0;
    clear();
  }

---

Member Function Documentation

void G4INCL::Store::add

(

IAvatar *

a

)

Add one avatar to the store

Avatars know about the particles they are associated with. Adding an avatar consists of the following steps: 1. Add the new avatar to the avatar list 2. Add any related new particles to the store by calling add(Particle*) (this should not happen, by the time we are adding avatars all particles should have already been added) 3. Connect the particles involved to the avatar in the map: particleAvatarConnections :: ParticleID -> [AvatarID] 4. Add the new avatar to the map: avatarParticleConnections :: AvatarID -> [ParticleID]

Definition at line 89 of file G4INCLStore.cc.

References add(), ERROR, G4INCL::IAvatar::getID(), and G4INCL::IAvatar::getParticles().

                            {
    // Add the avatar to the avatar map
    avatars[a->getID()]=a;
    avatarList.push_back(a);

    ParticleList pList = a->getParticles();
    for(ParticleIter i = pList.begin(); i != pList.end(); ++i) {
      // If one of the particles participating in this avatar haven't
      // been registered with the store we can do it now. On the other
      // hand, if this happens, it's probably a symptom of a bug
      // somewhere...
      if(particles.find((*i)->getID()) == particles.end()) {
       ERROR("Avatar was added before related particles. This is probably a bug." << std::endl);
        add((*i));
      }
      // Connect each particle to the avatar
      connectParticleAndAvatar((*i)->getID(), a->getID());
    }

  }

void G4INCL::Store::add

(

Particle *

p

)

Add one particle to the store.

Particle objects don't know anything about avatars so this method will only do two things: 1. add the particle to the particle map ParticleID -> Particle* 2. add an empty entry for this particle into map AvatarID -> [ParticleID]

Definition at line 60 of file G4INCLStore.cc.

References G4INCL::Particle::getID().

Referenced by add(), G4INCL::Nucleus::applyFinalState(), G4INCL::Nucleus::initializeParticles(), loadParticles(), particleHasEntered(), and G4INCL::StandardPropagationModel::registerAvatar().

                             {
    const long ID = p->getID();
    particles[ID]=p;
    inside.push_back(p);

    if(particleAvatarConnections.find(ID)==particleAvatarConnections.end()) {
      std::vector<long> *aIDs = new std::vector<long>;
      particleAvatarConnections[ID] = aIDs;
    }
  }

void G4INCL::Store::addIncomingParticle

(

Particle *const

p

)

Add a particle to the incoming list.

Parameters:

p

particle to add

Definition at line 110 of file G4INCLStore.cc.

Referenced by addParticleEntryAvatar().

                                                    {
    incoming.push_back(p);
  }

void G4INCL::Store::addParticleEntryAvatar

(

IAvatar *

a

)

Add one ParticleEntry avatar.

Definition at line 71 of file G4INCLStore.cc.

References addIncomingParticle(), G4INCL::IAvatar::getID(), and G4INCL::IAvatar::getParticles().

Referenced by addParticleEntryAvatars(), and G4INCL::StandardPropagationModel::shootParticle().

                                               {
    // Add the avatar to the avatar map
    avatars[a->getID()]=a;
    avatarList.push_back(a);

    ParticleList pList = a->getParticles();
    for(ParticleIter i = pList.begin(); i != pList.end(); ++i) {
      addIncomingParticle((*i));
      // Connect each particle to the avatar
      connectParticleAndAvatar((*i)->getID(), a->getID());
    }
  }

void G4INCL::Store::addParticleEntryAvatars

(

IAvatarList const &

al

)

Add one ParticleEntry avatar.

Definition at line 84 of file G4INCLStore.cc.

References addParticleEntryAvatar().

Referenced by G4INCL::StandardPropagationModel::shootComposite().

                                                           {
    for(IAvatarIter a=al.begin(); a!=al.end(); ++a)
      addParticleEntryAvatar(*a);
  }

void G4INCL::Store::addToOutgoing

(

ParticleList const &

pl

)

[inline]

Add a list of particles to the outgoing particle list.

Parameters:

pl

list of particles to be added

Definition at line 175 of file G4INCLStore.hh.

References addToOutgoing().

                                               {
      for(ParticleIter p=pl.begin(); p!=pl.end(); ++p)
        addToOutgoing(*p);
    }

void G4INCL::Store::addToOutgoing

(

Particle *

p

)

[inline]

add the particle to the outgoing particle list.

Parameters:

p

pointer to the particle to be added

Definition at line 169 of file G4INCLStore.hh.

Referenced by addToOutgoing(), G4INCL::Nucleus::applyFinalState(), G4INCL::Nucleus::decayMe(), G4INCL::Nucleus::decayOutgoingClusters(), G4INCL::Nucleus::decayOutgoingDeltas(), G4INCL::Nucleus::emitInsidePions(), G4INCL::Nucleus::finalizeProjectileRemnant(), and G4INCL::Nucleus::moveProjectileRemnantComponentsToOutgoing().

{ outgoing.push_back(p); }

static G4bool G4INCL::Store::avatarComparisonPredicate	(	IAvatar *	lhs,
		IAvatar *	rhs
	)			[inline, static]

Comparison predicate for avatars.

avatarComparisonPredicate is used by the std::sort or std::min_element functions to compare the avatar objects according to their time.

Parameters:

	lhs	pointer to the first avatar
	rhs	pointer to the second avatar

Returns:

true iff lhs' time is smaller than rhs'.

Definition at line 348 of file G4INCLStore.hh.

References G4INCL::IAvatar::getTime().

Referenced by findSmallestTime().

                                                                        {
      return (lhs->getTime() < rhs->getTime());
    }

void G4INCL::Store::clear

(

)

Clear all avatars and particles from the store.

Warning! This actually deletes the objects as well!

Definition at line 340 of file G4INCLStore.cc.

References clearAvatars(), clearOutgoing(), and WARN.

Referenced by loadParticles(), and ~Store().

                    {
    clearAvatars();
    inside.clear();

    for(std::map<long, Particle*>::iterator iter = particles.begin();
        iter != particles.end(); ++iter) {
      delete (*iter).second;
    }
    particles.clear();

    clearOutgoing();

    if( incoming.size() != 0 ) {
      WARN("Incoming list is not empty when Store::clear() is called" << std::endl);
    }
    incoming.clear();

#ifdef INCL_AVATAR_SEARCH_INCLSort
    avatarIterList.clear();
#endif

  }

void G4INCL::Store::clearAvatars

(

)

Clear avatars only.

Definition at line 316 of file G4INCLStore.cc.

Referenced by clear(), and G4INCL::StandardPropagationModel::propagate().

                           {
    for(std::map<long, IAvatar*>::iterator iter = avatars.begin();
        iter != avatars.end(); ++iter) {
      delete (*iter).second;
    }

    for(std::map<long, std::vector<long>*>::iterator iter = particleAvatarConnections.begin();
        iter != particleAvatarConnections.end(); ++iter) {
      delete (*iter).second;
    }

    particleAvatarConnections.clear();
    avatars.clear();
    avatarList.clear();

  }

void G4INCL::Store::clearIncoming

(

)

[inline]

Clear the incoming list.

Definition at line 128 of file G4INCLStore.hh.

Referenced by deleteIncoming().

                                {
      incoming.clear();
    }

void G4INCL::Store::clearOutgoing

(

)

Clear all outgoing particles from the store.

Warning! This actually deletes the objects as well!

Definition at line 363 of file G4INCLStore.cc.

Referenced by clear().

                            {
    for(ParticleIter iter = outgoing.begin();   iter != outgoing.end(); ++iter) {
      if((*iter)->isCluster()) {
        Cluster *c = dynamic_cast<Cluster *>(*iter);
// assert(c);
#ifdef INCLXX_IN_GEANT4_MODE
        if(!c)
          continue;
#endif
        c->deleteParticles();
      }
      delete (*iter);
    }
    outgoing.clear();
  }

G4bool G4INCL::Store::containsCollisions

(

)

const

Definition at line 492 of file G4INCLStore.cc.

References G4INCL::CollisionAvatarType.

                                         {
    IAvatarIter i;
    for(i = avatarList.begin(); i != avatarList.end(); ++i)
      if((*i)->getType()==CollisionAvatarType) return true;
    return false;
  }

G4int G4INCL::Store::countCascading

(

)

[inline]

Definition at line 243 of file G4INCLStore.hh.

References CLHEP::detail::n.

                           {
      G4int n=0;
      for(ParticleIter i=inside.begin(); i!=inside.end(); ++i) {
        if(!(*i)->isTargetSpectator())
          ++n;
      }
      return n;
    }

void G4INCL::Store::deleteIncoming

(

)

[inline]

Clear the incoming list and delete the particles.

Definition at line 133 of file G4INCLStore.hh.

References clearIncoming().

                                 {
      for(ParticleIter iter = incoming.begin(); iter != incoming.end(); ++iter) {
        delete (*iter);
      }
      clearIncoming();
    }

ParticleList G4INCL::Store::extractDynamicalSpectators

(

)

[inline]

Returns a list of dynamical spectators.

Looks in the outgoing list for particles without collisions and decays, removes them from outgoing and returns them in a list.

Returns:

the (possibly empty) list of dynamical spectators

Definition at line 211 of file G4INCLStore.hh.

                                              {
      ParticleList spectators;
      std::list<std::list<Particle*>::iterator> toBeErased;
      for(std::list<Particle*>::iterator p=outgoing.begin(); p!=outgoing.end(); ++p) {
        if((*p)->isProjectileSpectator()) {
// assert((*p)->isNucleon());
          spectators.push_back(*p); // add them to the list we will return
          toBeErased.push_back(p); // we will remove them from outgoing later
        }
      }

      // Now erase them from outgoing
      for(std::list<std::list<Particle*>::iterator>::iterator i=toBeErased.begin();
          i!=toBeErased.end(); ++i) {
        outgoing.erase(*i);
      }

      return spectators;
    }

IAvatar * G4INCL::Store::findSmallestTime

(

)

Find the avatar that has the smallest time.

Definition at line 210 of file G4INCLStore.cc.

References avatarComparisonPredicate(), and G4INCL::IAvatar::getID().

Referenced by G4INCL::StandardPropagationModel::propagate().

                                   {
    if(avatarList.empty()) return NULL;

#ifdef INCL_AVATAR_SEARCH_FullSort

    /* Full sort algorithm.
     *
     * Simple, but guaranteed to work.
     */
    avatarList.sort(Store::avatarComparisonPredicate);
    IAvatar *avatar = avatarList.front();

#elif defined(INCL_AVATAR_SEARCH_INCLSort)

    /* Partial sort algorithm used by INCL4.6.
     *
     * It nevers sorts the whole avatar list, but rather starts from the last
     * best avatar. It requires the avatarList to be updated by appending new
     * avatars at the end.
     */

    IAvatarIter best;
    if(avatarIterList.empty())
      best = avatarList.begin();
    else
      best = avatarIterList.back();
    G4double bestTime = (*best)->getTime();
    IAvatarIter a = best;

    for(++a; a!=avatarList.end(); ++a)
      if((*a)->getTime() < bestTime) {
        best = a;
        bestTime = (*best)->getTime();
        avatarIterList.push_back(best);
      }
    IAvatar *avatar = *best;

#elif defined(INCL_AVATAR_SEARCH_MinElement)

    /* Algorithm provided by the C++ stdlib. */
    IAvatar *avatar = *(std::min_element(avatarList.begin(), avatarList.end(),
          Store::avatarComparisonPredicate));

#else
#error Unrecognized INCL_AVATAR_SEARCH. Allowed values are: FullSort, INCLSort, MinElement.
#endif

    removeAvatarByID(avatar->getID());
    return avatar;
  }

IAvatarList const& G4INCL::Store::getAvatars

(

)

const [inline]

Return the list of avatars

Definition at line 109 of file G4INCLStore.hh.

                                          {
      return avatarList;
    }

Book* G4INCL::Store::getBook

(

)

[inline]

Return the pointer to the Book object which keeps track of various counters.

Definition at line 241 of file G4INCLStore.hh.

Referenced by G4INCL::Nucleus::applyFinalState(), G4INCL::Nucleus::decayOutgoingDeltas(), G4INCL::Nucleus::emitInsidePions(), G4INCL::Nucleus::fillEventInfo(), G4INCL::StandardPropagationModel::generateBinaryCollisionAvatar(), G4INCL::BinaryCollisionAvatar::getChannel(), G4INCL::INCL::initializeTarget(), G4INCL::Nucleus::insertParticle(), G4INCL::PauliStrictStandard::isBlocked(), G4INCL::SurfaceAvatar::postInteraction(), G4INCL::InteractionAvatar::postInteraction(), G4INCL::DecayAvatar::postInteraction(), G4INCL::BinaryCollisionAvatar::postInteraction(), G4INCL::StandardPropagationModel::propagate(), and G4INCL::InteractionAvatar::shouldUseLocalEnergy().

{return theBook; };

Config const* G4INCL::Store::getConfig

(

)

[inline]

Get the config object

Definition at line 255 of file G4INCLStore.hh.

Referenced by G4INCL::ClusteringModelIntercomparison::getCluster(), G4INCL::InteractionAvatar::postInteraction(), and G4INCL::InteractionAvatar::shouldUseLocalEnergy().

{ return theConfig; };

ParticleList const& G4INCL::Store::getIncomingParticles

(

)

const [inline]

Return the list of incoming particles (i.e. particles that have yet to enter the cascade).

Definition at line 196 of file G4INCLStore.hh.

{ return incoming; }

G4int G4INCL::Store::getLoadedA

(

)

[inline]

Get the value of the nucleus mass number that we read from file with loadParticles.

Definition at line 307 of file G4INCLStore.hh.

{ return loadedA; };

G4double G4INCL::Store::getLoadedStoppingTime

(

)

[inline]

Get the value of the stopping time that we read from file with loadParticles.

Definition at line 319 of file G4INCLStore.hh.

{ return loadedStoppingTime; };

G4int G4INCL::Store::getLoadedZ

(

)

[inline]

Get the value of the nucleus charge number that we read from file with loadParticles.

Definition at line 313 of file G4INCLStore.hh.

{ return loadedZ; };

ParticleList const& G4INCL::Store::getOutgoingParticles

(

)

const [inline]

Return the list of outgoing particles (i.e. particles that have left the cascade).

Definition at line 202 of file G4INCLStore.hh.

Referenced by G4INCL::Nucleus::computeRecoilKinematics(), G4INCL::Nucleus::computeSeparationEnergyBalance(), G4INCL::Nucleus::decayOutgoingClusters(), G4INCL::Nucleus::decayOutgoingDeltas(), G4INCL::Nucleus::fillEventInfo(), G4INCL::Nucleus::getConservationBalance(), G4INCL::Nucleus::isEventTransparent(), and G4INCL::Nucleus::print().

{ return outgoing; }

ParticleList G4INCL::Store::getParticipants

(

)

Get list of participants (active nucleons).

Warning: This (slow) method may be deprecated in the near future...

Definition at line 292 of file G4INCLStore.cc.

References WARN.

Referenced by G4INCL::Nucleus::dump(), and G4INCL::Nucleus::print().

                                      {
    WARN("Store::getParticipants is probably slow..." << std::endl);
    ParticleList result;
    for(std::map<long, Particle*>::iterator iter = particles.begin();
        iter != particles.end(); ++iter) {
      if((*iter).second->isParticipant()) {
        result.push_back((*iter).second);
      }
    }
    return result;
  }

ParticleList const& G4INCL::Store::getParticles

(

)

const [inline]

Return the list of "active" particles (i.e. particles that can participate in collisions).

Definition at line 235 of file G4INCLStore.hh.

Referenced by G4INCL::Nucleus::computeCenterOfMass(), G4INCL::Nucleus::computeTotalEnergy(), G4INCL::Nucleus::containsDeltas(), G4INCL::Nucleus::decayInsideDeltas(), G4INCL::Nucleus::emitInsidePions(), G4INCL::StandardPropagationModel::generateAllAvatars(), G4INCL::PauliStandard::getBlockingProbability(), G4INCL::ClusteringModelIntercomparison::getCluster(), G4INCL::CDPP::isBlocked(), and G4INCL::StandardPropagationModel::updateAvatars().

{ return inside; }

ParticleList G4INCL::Store::getSpectators

(

)

Get list of spectators (active nucleons).

Warning: This (slow) method may be deprecated in the near future...

Definition at line 304 of file G4INCLStore.cc.

References WARN.

Referenced by G4INCL::Nucleus::print().

                                    {
    WARN("Store::getSpectators is probably slow..." << std::endl);
    ParticleList result;
    for(std::map<long, Particle*>::iterator iter = particles.begin();
        iter != particles.end(); ++iter) {
      if(!(*iter).second->isParticipant()) {
        result.push_back((*iter).second);
      }
    }
    return result;
  }

void G4INCL::Store::initialiseParticleAvatarConnections

(

)

Initialise the particleAvatarConnections map.

Generate an empty avatar-ID vector for each particle in the inside list and fill in the relevant particle-avatar map entry.

Definition at line 333 of file G4INCLStore.cc.

Referenced by G4INCL::StandardPropagationModel::propagate().

                                                  {
    for(ParticleIter ip=inside.begin(); ip!=inside.end(); ++ip) {
      std::vector<long> *aIDs = new std::vector<long>;
      particleAvatarConnections[(*ip)->getID()] = aIDs;
    }
  }

void G4INCL::Store::loadParticles

(

std::string

filename

)

Load particle configuration from ASCII file (see avatarPredictionTest).

Definition at line 379 of file G4INCLStore.cc.

References add(), clear(), FATAL, G4INCL::Book::incrementCascading(), G4INCL::Particle::makeParticipant(), G4INCL::Neutron, G4INCL::Proton, G4INCL::Particle::setPotentialEnergy(), and G4INCL::UnknownParticle.

                                              {
    clear();
    G4int projectileA, projectileZ, A, Z;
    G4double stoppingTime, cutNN;
    G4int ID, type, isParticipant;
    G4double x, y, z;
    G4double px, py, pz, E, v;

    std::ifstream in(filename.c_str());
    in >> projectileA >> projectileZ >> A >> Z >> stoppingTime >> cutNN;
    loadedA = A;
    loadedZ = Z;
    loadedStoppingTime = stoppingTime;

    G4int readA = 0;
    G4int readZ = 0;
    while(1) {
      in >> ID >> type >> isParticipant >> x >> y >> z >> px >> py >> pz >> E >> v;
      if(!in.good()) break;
      ParticleType t;
      if(type == 1) {
        t = Proton;
        readZ++;
        readA++;
      }
      else if(type == -1) {
        t = Neutron;
        readA++;
      }
      else {
        FATAL("Unrecognized particle type while loading particles; type=" << type << std::endl);
        t = UnknownParticle;
      }

      Particle *p = new Particle(t, E, ThreeVector(px, py, pz),
                                 ThreeVector(x, y, z));
      p->setPotentialEnergy(v);
      if(isParticipant == 1) {
        p->makeParticipant();
        theBook->incrementCascading();
      }
      add(p);
    }
    
    in.close();
  }

void G4INCL::Store::particleHasBeenDestroyed

(

long

)

Remove the particle from the system. This also removes all avatars related to this particle.

Definition at line 278 of file G4INCLStore.cc.

References particleHasBeenUpdated().

Referenced by G4INCL::Nucleus::applyFinalState().

                                              {
    particleHasBeenUpdated(ID);
    // Have to destroy the particle here, the Store will forget about it
    Particle * const toDelete = particles.find(ID)->second;
    inside.remove(toDelete);
    delete toDelete;
    particles.erase(ID);
  }

void G4INCL::Store::particleHasBeenEjected

(

long

)

Mark the particle as ejected. This removes it from the list of inside particles and removes all avatars related to this particle.

Definition at line 269 of file G4INCLStore.cc.

References particleHasBeenUpdated().

Referenced by G4INCL::Nucleus::applyFinalState(), and G4INCL::Nucleus::emitInsidePions().

                                            {
    particleHasBeenUpdated(ID);
    // The particle will be destroyed when destroying the Store
    inside.remove(particles.find(ID)->second);
    particles.erase(ID);
    delete particleAvatarConnections.find(ID)->second;
    particleAvatarConnections.erase(ID);
  }

void G4INCL::Store::particleHasBeenUpdated

(

long

)

Notify the Store that a particle has been updated. This triggers the removal of obsolete avatars and their disconnection from the particle.

Definition at line 162 of file G4INCLStore.cc.

Referenced by G4INCL::Nucleus::applyFinalState(), particleHasBeenDestroyed(), and particleHasBeenEjected().

                                                    {
    std::vector<long> temp_aIDs;
    std::vector<long> *aIDs = particleAvatarConnections.find(particleID)->second;
    for(std::vector<long>::iterator i = aIDs->begin();
        i != aIDs->end(); ++i) {
      temp_aIDs.push_back((*i));
    }

    for(std::vector<long>::iterator i = temp_aIDs.begin();
        i != temp_aIDs.end(); ++i) {
      IAvatar *tmpAvatar = avatars.find(*i)->second;
      removeAvatarByID((*i));
      delete tmpAvatar;
    }
  }

void G4INCL::Store::particleHasEntered

(

Particle *const

particle

)

Move a particle from incoming to inside.

Parameters:

particle

pointer to a particle

Definition at line 287 of file G4INCLStore.cc.

References add(), and removeFromIncoming().

Referenced by G4INCL::Nucleus::insertParticle().

                                                          {
    removeFromIncoming(particle);
    add(particle);
  }

std::string G4INCL::Store::printAvatars

(

)

Print the list of avatars

Definition at line 483 of file G4INCLStore.cc.

                                {
    std::stringstream ss;
    IAvatarIter i;
    for(i = avatarList.begin(); i != avatarList.end(); ++i) {
      ss << (*i)->toString() << std::endl;
    }
    return ss.str();
  }

std::string G4INCL::Store::printParticleConfiguration

(

)

Print the nucleon configuration of the nucleus.

Definition at line 426 of file G4INCLStore.cc.

References G4INCL::Neutron, and G4INCL::Proton.

Referenced by writeParticles().

                                              {
    std::stringstream ss;
    G4int A = 0, Z = 0;
    for(ParticleIter i = inside.begin(); i != inside.end(); ++i) {
      if((*i)->getType() == Proton) {
        A++;
        Z++;
      }
      if((*i)->getType() == Neutron) {
        A++;
      }
    }
    // Note: Projectile A and Z are set to 0 (we don't really know
    // anything about them at this point).
    ss << "0 0 " << A << " " << Z << " "
              << "100.0" << " "
              << "0.0" << std::endl;

    for(ParticleIter i = inside.begin(); i != inside.end(); ++i) {
      G4int ID = (*i)->getID();
      G4int type = 0;
      if((*i)->getType() == Proton) {
        type = 1;
      }
      if((*i)->getType() == Neutron) {
        type = -1;
      }

      G4int isParticipant = 0;
      if((*i)->isParticipant()) {
        isParticipant = 1;
      }

      G4double x = (*i)->getPosition().getX();
      G4double y = (*i)->getPosition().getY();
      G4double z = (*i)->getPosition().getZ();
      G4double E = (*i)->getEnergy();
      G4double px = (*i)->getMomentum().getX();
      G4double py = (*i)->getMomentum().getY();
      G4double pz = (*i)->getMomentum().getZ();
      G4double V = (*i)->getPotentialEnergy();

      ss << ID << " " << type << " " << isParticipant << " "
                << x << " " << y << " " << z << " "
                << px << " " << py << " " << pz << " "
         << E << " " << V << std::endl;
    }

    return ss.str();
  }

void G4INCL::Store::removeFromIncoming

(

Particle *const

p

)

[inline]

Add a particle to the incoming list.

Parameters:

p

particle to add

Definition at line 125 of file G4INCLStore.hh.

Referenced by particleHasEntered().

{ incoming.remove(p); }

void G4INCL::Store::timeStep

(

G4double

step

)

Make one time step: propagate particles and subtract the length of the step from the avatar times.

Definition at line 261 of file G4INCLStore.cc.

Referenced by G4INCL::StandardPropagationModel::propagate().

                                    {
    for(std::map<long, Particle*>::iterator particleIter
          = particles.begin();
        particleIter != particles.end(); ++particleIter) {
      (*particleIter).second->propagate(step);
    }
  }

void G4INCL::Store::writeParticles

(

std::string

filename

)

Print the nucleon configuration of the nucleus.

Definition at line 477 of file G4INCLStore.cc.

References printParticleConfiguration().

                                               {
    std::ofstream out(filename.c_str());
    out << printParticleConfiguration();
    out.close();
  }

---

The documentation for this class was generated from the following files:

• G4INCLStore.hh
• G4INCLStore.cc

---