G4INCL::StandardPropagationModel Class Reference

#include <G4INCLStandardPropagationModel.hh>

Inheritance diagram for G4INCL::StandardPropagationModel:

Public Member Functions
	StandardPropagationModel (LocalEnergyType localEnergyType, LocalEnergyType localEnergyDeltaType)
virtual	~StandardPropagationModel ()
G4double	getCurrentTime ()
void	setNucleus (G4INCL::Nucleus *nucleus)
G4INCL::Nucleus *	getNucleus ()
G4double	shoot (ParticleSpecies const projectileSpecies, const G4double kineticEnergy, const G4double impactParameter, const G4double phi)
G4double	shootParticle (ParticleType const t, const G4double kineticEnergy, const G4double impactParameter, const G4double phi)
G4double	shootComposite (ParticleSpecies const s, const G4double kineticEnergy, const G4double impactParameter, const G4double phi)
void	setStoppingTime (G4double)
G4double	getStoppingTime ()
void	registerAvatar (G4INCL::IAvatar *anAvatar)
IAvatar *	generateBinaryCollisionAvatar (Particle *const p1, Particle *const p2) const
	Generate a two-particle avatar.
G4double	getReflectionTime (G4INCL::Particle const *const aParticle)
	Get the reflection time.
G4double	getTime (G4INCL::Particle const *const particleA, G4INCL::Particle const *const particleB, G4double *minDistOfApproach) const
void	generateUpdatedCollisions (const ParticleList &updatedParticles, const ParticleList &particles)
	Generate and register collisions between a list of updated particles and all the other particles.
void	generateCollisions (const ParticleList &particles, const ParticleList &except)
	Generate and register collisions among particles in a list, except between those in another list.
void	generateDecays (const ParticleList &particles)
	Generate decays for particles that can decay.
void	updateAvatars (const ParticleList &particles)
void	generateAllAvatars (G4bool excludeUpdated=false)
	(Re)Generate all possible avatars.
G4INCL::IAvatar *	propagate ()
Public Member Functions inherited from G4INCL::IPropagationModel
	IPropagationModel ()
virtual	~IPropagationModel ()
virtual void	setNucleus (G4INCL::Nucleus *nucleus)=0
virtual G4INCL::Nucleus *	getNucleus ()=0
virtual G4double	shoot (ParticleSpecies const projectileSpecies, const G4double kineticEnergy, const G4double impactParameter, const G4double phi)=0
virtual G4double	getCurrentTime ()=0
virtual void	setStoppingTime (G4double)=0
virtual G4double	getStoppingTime ()=0
virtual G4INCL::IAvatar *	propagate ()=0

Additional Inherited Members
virtual G4double	shootParticle (ParticleType const t, const G4double kineticEnergy, const G4double impactParameter, const G4double phi)=0
virtual G4double	shootComposite (ParticleSpecies const s, const G4double kineticEnergy, const G4double impactParameter, const G4double phi)=0

Detailed Description

Standard INCL4 particle propagation and avatar prediction

This class implements the standard INCL4 avatar prediction and particle propagation logic. The main idea is to predict all collisions between particles and their reflections from the potential wall. After this we select the avatar with the smallest time, propagate all particles to their positions at that time and return the avatar to the INCL kernel

See also

G4INCL::Kernel.

The particle trajectories in this propagation model are straight lines and all particles are assumed to move with constant velocity.

Definition at line 70 of file G4INCLStandardPropagationModel.hh.

Constructor & Destructor Documentation

StandardPropagationModel()

G4INCL::StandardPropagationModel::StandardPropagationModel	(	LocalEnergyType	localEnergyType,
		LocalEnergyType	localEnergyDeltaType
	)

Definition at line 64 of file G4INCLStandardPropagationModel.cc.

      :theNucleus(0), maximumTime(70.0), currentTime(0.0), firstAvatar(true),
      theLocalEnergyType(localEnergyType),
      theLocalEnergyDeltaType(localEnergyDeltaType)
    {
    }

~StandardPropagationModel()

G4INCL::StandardPropagationModel::~StandardPropagationModel ( )

virtual

Definition at line 71 of file G4INCLStandardPropagationModel.cc.

    {
      delete theNucleus;
    }

Member Function Documentation

generateAllAvatars()

void G4INCL::StandardPropagationModel::generateAllAvatars

(

G4bool

excludeUpdated = false

)

(Re)Generate all possible avatars.

Parameters

excludeUpdated

exclude collisions between updated particles.

Definition at line 428 of file G4INCLStandardPropagationModel.cc.

                                                                           {
      ParticleList particles = theNucleus->getStore()->getParticles();
      if(particles.empty()) { ERROR("No particles inside the nucleus!" << std::endl); }
      for(ParticleIter i = particles.begin(); i != particles.end(); ++i) {
        G4double time = this->getReflectionTime(*i);
        if(time <= maximumTime) registerAvatar(new SurfaceAvatar(*i, time, theNucleus));
      }
      ParticleList except;
      if(excludeUpdated)
        except = theNucleus->getUpdatedParticles();
      generateCollisions(particles,except);
      generateDecays(particles);
    }

Referenced by propagate(), shootComposite(), and shootParticle().

generateBinaryCollisionAvatar()

IAvatar * G4INCL::StandardPropagationModel::generateBinaryCollisionAvatar	(	Particle *const	p1,
		Particle *const	p2
	)		const

Generate a two-particle avatar.

Generate a two-particle avatar, if all the appropriate conditions are met.

Definition at line 265 of file G4INCLStandardPropagationModel.cc.

                                                                                                                   {
      // Is either particle a participant?
      if(!p1->isParticipant() && !p2->isParticipant() && p1->getParticipantType()==p2->getParticipantType()) return NULL;
 
      // Is it a pi-resonance collision (we don't treat them)?
      if((p1->isResonance() && p2->isPion()) || (p1->isPion() && p2->isResonance()))
        return NULL;
 
      // Will the avatar take place between now and the end of the cascade?
      G4double minDistOfApproachSquared = 0.0;
      G4double t = getTime(p1, p2, &minDistOfApproachSquared);
      if(t>maximumTime || t<currentTime) return NULL;
 
      // Local energy. Jump through some hoops to calculate the cross section
      // at the collision point, and clean up after yourself afterwards.
      G4bool hasLocalEnergy;
      if(p1->isPion() || p2->isPion())
        hasLocalEnergy = ((theLocalEnergyDeltaType == FirstCollisionLocalEnergy &&
              theNucleus->getStore()->getBook()->getAcceptedCollisions()==0) ||
            theLocalEnergyDeltaType == AlwaysLocalEnergy);
      else
        hasLocalEnergy = ((theLocalEnergyType == FirstCollisionLocalEnergy &&
              theNucleus->getStore()->getBook()->getAcceptedCollisions()==0) ||
            theLocalEnergyType == AlwaysLocalEnergy);
      const G4bool p1HasLocalEnergy = (hasLocalEnergy && !p1->isPion());
      const G4bool p2HasLocalEnergy = (hasLocalEnergy && !p2->isPion());
 
      Particle backupParticle1 = *p1;
      if(p1HasLocalEnergy) {
        p1->propagate(t - currentTime);
        if(p1->getPosition().mag() > theNucleus->getSurfaceRadius(p1)) {
          *p1 = backupParticle1;
          return NULL;
        }
        KinematicsUtils::transformToLocalEnergyFrame(theNucleus, p1);
      }
      Particle backupParticle2 = *p2;
      if(p2HasLocalEnergy) {
        p2->propagate(t - currentTime);
        if(p2->getPosition().mag() > theNucleus->getSurfaceRadius(p2)) {
          *p2 = backupParticle2;
          if(p1HasLocalEnergy) {
            *p1 = backupParticle1;
          }
          return NULL;
        }
        KinematicsUtils::transformToLocalEnergyFrame(theNucleus, p2);
      }
 
      // Compute the total cross section
      const G4double totalCrossSection = CrossSections::total(p1, p2);
      const G4double squareTotalEnergyInCM = KinematicsUtils::squareTotalEnergyInCM(p1,p2);
 
      // Restore particles to their state before the local-energy tweak
      if(p1HasLocalEnergy) {
        *p1 = backupParticle1;
      }
      if(p2HasLocalEnergy) {
        *p2 = backupParticle2;
      }
 
      // Is the CM energy > cutNN? (no cutNN on the first collision)
      if(theNucleus->getStore()->getBook()->getAcceptedCollisions()>0
          && p1->isNucleon() && p2->isNucleon()
          && squareTotalEnergyInCM < BinaryCollisionAvatar::cutNNSquared) return NULL;
 
      // Do the particles come close enough to each other?
      if(Math::tenPi*minDistOfApproachSquared > totalCrossSection) return NULL;
 
      // Bomb out if the two collision partners are the same particle
// assert(p1->getID() != p2->getID());
 
      // Return a new avatar, then!
      return new G4INCL::BinaryCollisionAvatar(t, totalCrossSection, theNucleus, p1, p2);
    }

Referenced by generateCollisions(), and generateUpdatedCollisions().

generateCollisions()

void G4INCL::StandardPropagationModel::generateCollisions	(	const ParticleList &	particles,
		const ParticleList &	except
	)

Generate and register collisions among particles in a list, except between those in another list.

This method generates all possible collisions among the particles. Each collision is generated only once. The collision is NOT generated if BOTH collision partners belong to the except list.

You should pass an empty list as the except parameter if you want to generate all possible collisions among particles.

Parameters

particles	list of particles
except	list of excluded particles

Definition at line 397 of file G4INCLStandardPropagationModel.cc.

                                                                                                               {
 
      G4bool haveExcept;
      haveExcept=(except.size()!=0);
 
      // Loop over all the particles
      for(ParticleIter p1 = particles.begin(); p1 != particles.end(); ++p1)
      {
        // Loop over the rest of the particles
        ParticleIter p2 = p1;
        for(++p2; p2 != particles.end(); ++p2)
        {
          // Skip the collision if both particles must be excluded
          if(haveExcept && (*p1)->isInList(except) && (*p2)->isInList(except)) continue;
 
          registerAvatar(generateBinaryCollisionAvatar(*p1,*p2));
        }
      }
 
    }

Referenced by generateAllAvatars().

generateDecays()

void G4INCL::StandardPropagationModel::generateDecays

(

const ParticleList &

particles

)

Generate decays for particles that can decay.

The list of particles given as an argument is allowed to contain also stable particles.

Parameters

particles

list of particles to (possibly) generate decays for

Definition at line 442 of file G4INCLStandardPropagationModel.cc.

                                                                               {
      for(ParticleIter i = particles.begin(); i != particles.end(); ++i) {
    if((*i)->isDelta()) {
    G4double decayTime = DeltaDecayChannel::computeDecayTime((*i));
      G4double time = currentTime + decayTime;
      if(time <= maximumTime) {
        registerAvatar(new DecayAvatar((*i), time, theNucleus));
      }
    }
      }
    }

Referenced by generateAllAvatars(), and propagate().

generateUpdatedCollisions()

void G4INCL::StandardPropagationModel::generateUpdatedCollisions	(	const ParticleList &	updatedParticles,
		const ParticleList &	particles
	)

Generate and register collisions between a list of updated particles and all the other particles.

This method does not generate collisions among the particles in updatedParticles; in other words, it generates a collision between one of the updatedParticles and one of the particles ONLY IF the latter does not belong to updatedParticles.

If you intend to generate all possible collisions among particles in a list, use generateCollisions().

Parameters

updatedParticles	list of updated particles
particles	list of particles

Definition at line 378 of file G4INCLStandardPropagationModel.cc.

                                                                                                                                {
 
      // Loop over all the updated particles
      for(ParticleIter updated = updatedParticles.begin(); updated != updatedParticles.end(); ++updated)
      {
        // Loop over all the particles
        for(ParticleIter particle = particles.begin(); particle != particles.end(); ++particle)
        {
          /* Consider the generation of a collision avatar only if (*particle)
           * is not one of the updated particles.
           * The criterion makes sure that you don't generate avatars between
           * updated particles. */
          if((*particle)->isInList(updatedParticles)) continue;
 
          registerAvatar(generateBinaryCollisionAvatar(*particle,*updated));
        }
      }
    }

Referenced by updateAvatars().

getCurrentTime()

G4double G4INCL::StandardPropagationModel::getCurrentTime ( )

virtual

Returns the current global time of the system.

Implements G4INCL::IPropagationModel.

Definition at line 251 of file G4INCLStandardPropagationModel.cc.

                                                      {
      return currentTime;
    }

getNucleus()

G4INCL::Nucleus * G4INCL::StandardPropagationModel::getNucleus ( )

virtual

Get the nucleus.

Implements G4INCL::IPropagationModel.

Definition at line 76 of file G4INCLStandardPropagationModel.cc.

    {
      return theNucleus;
    }

getReflectionTime()

G4double G4INCL::StandardPropagationModel::getReflectionTime

(

G4INCL::Particle const *const

aParticle

)

Get the reflection time.

Returns the reflection time of a particle on the potential wall.

Parameters

aParticle

pointer to the particle

Definition at line 341 of file G4INCLStandardPropagationModel.cc.

                                                                                               {
      Intersection theIntersection(
          IntersectionFactory::getLaterTrajectoryIntersection(
            aParticle->getPosition(),
            aParticle->getPropagationVelocity(),
            theNucleus->getSurfaceRadius(aParticle)));
      G4double time;
      if(theIntersection.exists) {
        time = currentTime + theIntersection.time;
      } else {
        ERROR("Imaginary reflection time for particle: " << std::endl
          << aParticle->print());
        time = 10000.0;
      }
      return time;
    }

Referenced by generateAllAvatars(), and updateAvatars().

getStoppingTime()

G4double G4INCL::StandardPropagationModel::getStoppingTime ( )

virtual

Get the current stopping time.

Implements G4INCL::IPropagationModel.

Definition at line 242 of file G4INCLStandardPropagationModel.cc.

                                                       {
      return maximumTime;
    }

getTime()

G4double G4INCL::StandardPropagationModel::getTime	(	G4INCL::Particle const *const	particleA,
		G4INCL::Particle const *const	particleB,
		G4double *	minDistOfApproach
	)		const

Get the predicted time of the collision between two particles.

Definition at line 358 of file G4INCLStandardPropagationModel.cc.

    {
      G4double time;
      G4INCL::ThreeVector t13 = particleA->getPropagationVelocity();
      t13 -= particleB->getPropagationVelocity();
      G4INCL::ThreeVector distance = particleA->getPosition();
      distance -= particleB->getPosition();
      const G4double t7 = t13.dot(distance);
      const G4double dt = t13.mag2();
      if(dt <= 1.0e-10) {
        (*minDistOfApproach) = 100000.0;
        return currentTime + 100000.0;
      } else {
        time = -t7/dt;
      }
      (*minDistOfApproach) = distance.mag2() + time * t7;
      return currentTime + time;
    }

Referenced by generateBinaryCollisionAvatar().

propagate()

G4INCL::IAvatar * G4INCL::StandardPropagationModel::propagate ( )

virtual

Propagate all particles and return the first avatar.

Implements G4INCL::IPropagationModel.

Definition at line 454 of file G4INCLStandardPropagationModel.cc.

    {
      // We update only the information related to particles that were updated
      // by the previous avatar.
#ifdef INCL_REGENERATE_AVATARS
#warning "The INCL_REGENERATE_AVATARS code has not been tested in a while. Use it at your peril."
      if(theNucleus->getUpdatedParticles().size()!=0 || theNucleus->getCreatedParticles().size()!=0) {
        // Regenerates the entire avatar list, skipping collisions between
        // updated particles
        theNucleus->getStore()->clearAvatars();
        theNucleus->getStore()->initialiseParticleAvatarConnections();
        generateAllAvatars(true);
      }
#else
      // Deltas are created by transforming nucleon into a delta for
      // efficiency reasons
      Particle * const blockedDelta = theNucleus->getBlockedDelta();
      ParticleList updatedParticles = theNucleus->getUpdatedParticles();
      if(blockedDelta)
        updatedParticles.push_back(blockedDelta);
      generateDecays(updatedParticles);
 
      ParticleList needNewAvatars = theNucleus->getUpdatedParticles();
      ParticleList created = theNucleus->getCreatedParticles();
      needNewAvatars.splice(needNewAvatars.end(), created);
      updateAvatars(needNewAvatars);
#endif
 
      G4INCL::IAvatar *theAvatar = theNucleus->getStore()->findSmallestTime();
      if(theAvatar == 0) return 0; // Avatar list is empty
      //      theAvatar->dispose();
 
      if(theAvatar->getTime() < currentTime) {
        ERROR("Avatar time = " << theAvatar->getTime() << ", currentTime = " << currentTime << std::endl);
        return 0;
      } else if(theAvatar->getTime() > currentTime) {
        theNucleus->getStore()->timeStep(theAvatar->getTime() - currentTime);
 
        currentTime = theAvatar->getTime();
        theNucleus->getStore()->getBook()->setCurrentTime(currentTime);
      }
 
      return theAvatar;
    }

registerAvatar()

void G4INCL::StandardPropagationModel::registerAvatar

(

G4INCL::IAvatar *

anAvatar

)

Add an avatar to the storage.

Definition at line 260 of file G4INCLStandardPropagationModel.cc.

    {
      if(anAvatar) theNucleus->getStore()->add(anAvatar);
    }

Referenced by generateAllAvatars(), generateCollisions(), generateDecays(), generateUpdatedCollisions(), and updateAvatars().

setNucleus()

void G4INCL::StandardPropagationModel::setNucleus ( G4INCL::Nucleus *  nucleus )

virtual

Set the nucleus for this propagation model.

Implements G4INCL::IPropagationModel.

Definition at line 255 of file G4INCLStandardPropagationModel.cc.

    {
      theNucleus = nucleus;
    }

setStoppingTime()

void G4INCL::StandardPropagationModel::setStoppingTime ( G4double  time )

virtual

Set the stopping time of the simulation.

Implements G4INCL::IPropagationModel.

Definition at line 246 of file G4INCLStandardPropagationModel.cc.

                                                                {
// assert(time>0.0);
      maximumTime = time;
    }

shoot()

G4double G4INCL::StandardPropagationModel::shoot ( ParticleSpecies const  projectileSpecies, const G4double  kineticEnergy, const G4double  impactParameter, const G4double  phi  )

virtual

Implements G4INCL::IPropagationModel.

Definition at line 81 of file G4INCLStandardPropagationModel.cc.

                                                                                                                                                                      {
      if(projectileSpecies.theType==Composite)
        return shootComposite(projectileSpecies, kineticEnergy, impactParameter, phi);
      else
        return shootParticle(projectileSpecies.theType, kineticEnergy, impactParameter, phi);
    }

shootComposite()

G4double G4INCL::StandardPropagationModel::shootComposite ( ParticleSpecies const  s, const G4double  kineticEnergy, const G4double  impactParameter, const G4double  phi  )

virtual

Implements G4INCL::IPropagationModel.

Definition at line 162 of file G4INCLStandardPropagationModel.cc.

                                                                                                                                                                     {
      theNucleus->setNucleusNucleusCollision();
      currentTime = 0.0;
 
      // Create the ProjectileRemnant object
      ProjectileRemnant *pr = new ProjectileRemnant(species, kineticEnergy);
 
      // Same stopping time as for nucleon-nucleus
      maximumTime = 29.8 * std::pow(theNucleus->getA(), 0.16);
 
      // If the incoming cluster is slow, use a larger stopping time
      const G4double rms = ParticleTable::getNuclearRadius(pr->getA(), pr->getZ());
      const G4double rMax = theNucleus->getUniverseRadius();
      const G4double distance = 2.*rMax + 2.725*rms;
      const G4double projectileVelocity = pr->boostVector().mag();
      const G4double traversalTime = distance / projectileVelocity;
      if(maximumTime < traversalTime)
        maximumTime = traversalTime;
      DEBUG("Cascade stopping time is " << maximumTime << std::endl);
 
      // If Coulomb is activated, do not process events with impact
      // parameter larger than the maximum impact parameter, taking into
      // account Coulomb distortion.
      if(impactParameter>CoulombDistortion::maxImpactParameter(pr,theNucleus)) {
        pr->deleteParticles();
        DEBUG("impactParameter>CoulombDistortion::maxImpactParameter" << std::endl);
        delete pr;
        return -1.;
      }
 
      // Position the cluster at the right impact parameter
      ThreeVector position(impactParameter * std::cos(phi),
          impactParameter * std::sin(phi),
          0.);
      pr->setPosition(position);
 
      /* Store the internal kinematics of the projectile remnant.
       *
       * Note that this is at variance with the Fortran version, which stores
       * the initial kinematics of the particles *after* putting the spectators
       * on mass shell, but *before* removing the same energy from the
       * participants. Due to the different code flow, doing so in the C++
       * version leads to wrong excitation energies for the forced compound
       * nucleus.
       */
      pr->storeComponents();
 
      // Fill in the relevant kinematic variables
      theNucleus->setIncomingAngularMomentum(pr->getAngularMomentum());
      theNucleus->setIncomingMomentum(pr->getMomentum());
      theNucleus->setInitialEnergy(pr->getEnergy()
          + ParticleTable::getTableMass(theNucleus->getA(),theNucleus->getZ()));
 
      generateAllAvatars();
      firstAvatar = false;
 
      // Get the entry avatars from Coulomb
      IAvatarList theAvatarList
        = CoulombDistortion::bringToSurface(pr, theNucleus);
 
      if(theAvatarList.empty()) {
        DEBUG("No ParticleEntryAvatar found, transparent event" << std::endl);
        pr->deleteParticles();
        delete pr;
        return -1.;
      }
 
      // Tell the Nucleus about the ProjectileRemnant
      theNucleus->setProjectileRemnant(pr);
 
      // Set the number of projectile particles
      theNucleus->setProjectileChargeNumber(pr->getZ());
      theNucleus->setProjectileMassNumber(pr->getA());
 
      // Register the ParticleEntryAvatars
      theNucleus->getStore()->addParticleEntryAvatars(theAvatarList);
 
      return pr->getTransversePosition().mag();
    }

Referenced by shoot().

shootParticle()

G4double G4INCL::StandardPropagationModel::shootParticle ( ParticleType const  t, const G4double  kineticEnergy, const G4double  impactParameter, const G4double  phi  )

virtual

Implements G4INCL::IPropagationModel.

Definition at line 88 of file G4INCLStandardPropagationModel.cc.

                                                                                                                                                              {
      theNucleus->setParticleNucleusCollision();
      currentTime = 0.0;
 
      // Create the projectile particle
      const G4double projectileMass = ParticleTable::getTableParticleMass(type);
      G4double energy = kineticEnergy + projectileMass;
      G4double momentumZ = std::sqrt(energy*energy - projectileMass*projectileMass);
      ThreeVector momentum(0.0, 0.0, momentumZ);
      Particle *p= new G4INCL::Particle(type, energy, momentum, ThreeVector());
 
      G4double temfin = 0.0;
      if( p->isNucleon() )
        temfin = 29.8 * std::pow(theNucleus->getA(), 0.16);
      else {
        const G4double tlab = p->getEnergy() - p->getMass();
        temfin = 30.18 * std::pow(theNucleus->getA(), 0.17*(1.0 - 5.7E-5*tlab));
      }
 
      maximumTime = temfin;
 
      // If the incoming particle is slow, use a larger stopping time
      const G4double rMax = theNucleus->getUniverseRadius();
      const G4double distance = 2.*rMax;
      const G4double projectileVelocity = p->boostVector().mag();
      const G4double traversalTime = distance / projectileVelocity;
      if(maximumTime < traversalTime)
        maximumTime = traversalTime;
      DEBUG("Cascade stopping time is " << maximumTime << std::endl);
 
      // If Coulomb is activated, do not process events with impact
      // parameter larger than the maximum impact parameter, taking into
      // account Coulomb distortion.
      if(impactParameter>CoulombDistortion::maxImpactParameter(p->getSpecies(), kineticEnergy, theNucleus)) {
        DEBUG("impactParameter>CoulombDistortion::maxImpactParameter" << std::endl);
        delete p;
        return -1.;
      }
 
      ThreeVector position(impactParameter * std::cos(phi),
          impactParameter * std::sin(phi),
          0.);
      p->setPosition(position);
 
      // Fill in the relevant kinematic variables
      theNucleus->setIncomingAngularMomentum(p->getAngularMomentum());
      theNucleus->setIncomingMomentum(p->getMomentum());
      theNucleus->setInitialEnergy(p->getEnergy()
          + ParticleTable::getTableMass(theNucleus->getA(),theNucleus->getZ()));
 
      // Reset the particle kinematics to the INCL values
      p->setINCLMass();
      p->setEnergy(p->getMass() + kineticEnergy);
      p->adjustMomentumFromEnergy();
 
      p->makeProjectileSpectator();
      generateAllAvatars();
      firstAvatar = false;
 
      // Get the entry avatars from Coulomb and put them in the Store
      ParticleEntryAvatar *theEntryAvatar = CoulombDistortion::bringToSurface(p, theNucleus);
      if(theEntryAvatar) {
        theNucleus->getStore()->addParticleEntryAvatar(theEntryAvatar);
 
        theNucleus->setProjectileChargeNumber(p->getZ());
        theNucleus->setProjectileMassNumber(p->getA());
 
        return p->getTransversePosition().mag();
      } else {
        delete p;
        return -1.;
      }
    }

Referenced by shoot().

updateAvatars()

void G4INCL::StandardPropagationModel::updateAvatars

(

const ParticleList &

particles

)

Update all avatars related to a particle.

Definition at line 418 of file G4INCLStandardPropagationModel.cc.

                                                                              {
 
      for(ParticleIter iter = particles.begin(); iter != particles.end(); ++iter) {
        G4double time = this->getReflectionTime(*iter);
        if(time <= maximumTime) registerAvatar(new SurfaceAvatar(*iter, time, theNucleus));
      }
      ParticleList const &p = theNucleus->getStore()->getParticles();
      generateUpdatedCollisions(particles, p);             // Predict collisions with spectators and participants
    }

Referenced by propagate().

---

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

• G4INCLStandardPropagationModel.hh
• G4INCLStandardPropagationModel.cc