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 64 of file G4INCLStore.hh.

Constructor & Destructor Documentation

Store()

G4INCL::Store::Store

(

Config const *const

config

)

Store constructor

Definition at line 46 of file G4INCLStore.cc.

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

~Store()

G4INCL::Store::~Store

(

)

Store destructor

Definition at line 55 of file G4INCLStore.cc.

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

Member Function Documentation

add() [1/2]

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 91 of file G4INCLStore.cc.

                            {
    // 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());
    }
 
  }

add() [2/2]

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 62 of file G4INCLStore.cc.

                             {
    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;
    }
  }

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

addIncomingParticle()

void G4INCL::Store::addIncomingParticle

(

Particle *const

p

)

Add a particle to the incoming list.

Parameters

p	particle to add

Definition at line 112 of file G4INCLStore.cc.

                                                    {
    incoming.push_back(p);
  }

Referenced by addParticleEntryAvatar().

addParticleEntryAvatar()

void G4INCL::Store::addParticleEntryAvatar

(

IAvatar *

a

)

Add one ParticleEntry avatar.

Definition at line 73 of file G4INCLStore.cc.

                                               {
    // 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());
    }
  }

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

addParticleEntryAvatars()

void G4INCL::Store::addParticleEntryAvatars

(

IAvatarList const &

al

)

Add one ParticleEntry avatar.

Definition at line 86 of file G4INCLStore.cc.

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

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

addToOutgoing() [1/2]

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 171 of file G4INCLStore.hh.

{ outgoing.push_back(p); }

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().

addToOutgoing() [2/2]

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 177 of file G4INCLStore.hh.

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

avatarComparisonPredicate()

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

inlinestatic

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 350 of file G4INCLStore.hh.

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

Referenced by findSmallestTime().

clear()

void G4INCL::Store::clear

(

)

Clear all avatars and particles from the store.

Warning! This actually deletes the objects as well!

Definition at line 342 of file G4INCLStore.cc.

                    {
    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
 
  }

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

clearAvatars()

void G4INCL::Store::clearAvatars

(

)

Clear avatars only.

Definition at line 318 of file G4INCLStore.cc.

                           {
    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();
 
  }

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

clearIncoming()

void G4INCL::Store::clearIncoming ( )

inline

Clear the incoming list.

Definition at line 130 of file G4INCLStore.hh.

                                {
      incoming.clear();
    }

Referenced by deleteIncoming().

clearOutgoing()

void G4INCL::Store::clearOutgoing

(

)

Clear all outgoing particles from the store.

Warning! This actually deletes the objects as well!

Definition at line 365 of file G4INCLStore.cc.

                            {
    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();
  }

Referenced by clear().

containsCollisions()

G4bool G4INCL::Store::containsCollisions

(

)

const

Definition at line 494 of file G4INCLStore.cc.

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

countCascading()

G4int G4INCL::Store::countCascading ( )

inline

Definition at line 245 of file G4INCLStore.hh.

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

deleteIncoming()

void G4INCL::Store::deleteIncoming ( )

inline

Clear the incoming list and delete the particles.

Definition at line 135 of file G4INCLStore.hh.

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

extractDynamicalSpectators()

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 213 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;
    }

findSmallestTime()

IAvatar * G4INCL::Store::findSmallestTime

(

)

Find the avatar that has the smallest time.

Definition at line 212 of file G4INCLStore.cc.

                                   {
    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;
  }

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

getAvatars()

IAvatarList const & G4INCL::Store::getAvatars ( ) const

inline

Return the list of avatars

Definition at line 111 of file G4INCLStore.hh.

                                          {
      return avatarList;
    }

getBook()

Book * G4INCL::Store::getBook ( )

inline

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

Definition at line 243 of file G4INCLStore.hh.

{return theBook; };

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::BinaryCollisionAvatar::postInteraction(), G4INCL::DecayAvatar::postInteraction(), G4INCL::InteractionAvatar::postInteraction(), G4INCL::SurfaceAvatar::postInteraction(), G4INCL::StandardPropagationModel::propagate(), and G4INCL::InteractionAvatar::shouldUseLocalEnergy().

getConfig()

Config const * G4INCL::Store::getConfig ( )

inline

Get the config object

Definition at line 257 of file G4INCLStore.hh.

{ return theConfig; };

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

getIncomingParticles()

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 198 of file G4INCLStore.hh.

{ return incoming; }

getLoadedA()

G4int G4INCL::Store::getLoadedA ( )

inline

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

Definition at line 309 of file G4INCLStore.hh.

{ return loadedA; };

getLoadedStoppingTime()

G4double G4INCL::Store::getLoadedStoppingTime ( )

inline

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

Definition at line 321 of file G4INCLStore.hh.

{ return loadedStoppingTime; };

getLoadedZ()

G4int G4INCL::Store::getLoadedZ ( )

inline

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

Definition at line 315 of file G4INCLStore.hh.

{ return loadedZ; };

getOutgoingParticles()

ParticleList const & G4INCL::Store::getOutgoingParticles ( ) const

inline

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

Definition at line 204 of file G4INCLStore.hh.

{ return outgoing; }

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().

getParticipants()

ParticleList G4INCL::Store::getParticipants

(

)

Get list of participants (active nucleons).

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

Definition at line 294 of file G4INCLStore.cc.

                                      {
    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;
  }

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

getParticles()

ParticleList const & G4INCL::Store::getParticles ( ) const

inline

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

Definition at line 237 of file G4INCLStore.hh.

{ return inside; }

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().

getSpectators()

ParticleList G4INCL::Store::getSpectators

(

)

Get list of spectators (active nucleons).

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

Definition at line 306 of file G4INCLStore.cc.

                                    {
    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;
  }

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

initialiseParticleAvatarConnections()

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 335 of file G4INCLStore.cc.

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

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

loadParticles()

void G4INCL::Store::loadParticles

(

std::string

filename

)

Load particle configuration from ASCII file (see avatarPredictionTest).

Definition at line 381 of file G4INCLStore.cc.

                                              {
    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);
        abort();
      }
 
      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();
  }

particleHasBeenDestroyed()

void G4INCL::Store::particleHasBeenDestroyed

(

long

ID

)

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

Definition at line 280 of file G4INCLStore.cc.

                                              {
    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);
  }

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

particleHasBeenEjected()

void G4INCL::Store::particleHasBeenEjected

(

long

ID

)

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 271 of file G4INCLStore.cc.

                                            {
    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);
  }

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

particleHasBeenUpdated()

void G4INCL::Store::particleHasBeenUpdated

(

long

particleID

)

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 164 of file G4INCLStore.cc.

                                                    {
    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;
    }
  }

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

particleHasEntered()

void G4INCL::Store::particleHasEntered

(

Particle *const

particle

)

Move a particle from incoming to inside.

Parameters

particle

pointer to a particle

Definition at line 289 of file G4INCLStore.cc.

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

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

printAvatars()

std::string G4INCL::Store::printAvatars

(

)

Print the list of avatars

Definition at line 485 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();
  }

printParticleConfiguration()

std::string G4INCL::Store::printParticleConfiguration

(

)

Print the nucleon configuration of the nucleus.

Definition at line 428 of file G4INCLStore.cc.

                                              {
    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();
  }

Referenced by writeParticles().

removeFromIncoming()

void G4INCL::Store::removeFromIncoming ( Particle *const  p )

inline

Add a particle to the incoming list.

Parameters

p	particle to add

Definition at line 127 of file G4INCLStore.hh.

{ incoming.remove(p); }

Referenced by particleHasEntered().

timeStep()

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 263 of file G4INCLStore.cc.

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

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

writeParticles()

void G4INCL::Store::writeParticles

(

std::string

filename

)

Print the nucleon configuration of the nucleus.

Definition at line 479 of file G4INCLStore.cc.

                                               {
    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