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)
	Generate a two-particle avatar. More...
G4double	getReflectionTime (G4INCL::Particle const *const aParticle)
	Get the reflection time. More...
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. More...
void	generateCollisions (const ParticleList &particles)
	Generate and register collisions among particles in a list, except between those in another list. More...
void	generateCollisions (const ParticleList &particles, const ParticleList &except)
	Generate and register collisions among particles in a list, except between those in another list. More...
void	generateDecays (const ParticleList &particles)
	Generate decays for particles that can decay. More...
void	updateAvatars (const ParticleList &particles)
void	generateAllAvatars ()
	(Re)Generate all possible avatars. More...
G4INCL::IAvatar *	propagate ()
Public Member Functions inherited from G4INCL::IPropagationModel
	IPropagationModel ()
virtual	~IPropagationModel ()

Additional Inherited Members

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 68 of file G4INCLStandardPropagationModel.hh.

Constructor & Destructor Documentation

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

Definition at line 62 of file G4INCLStandardPropagationModel.cc.

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

G4INCL::StandardPropagationModel::~StandardPropagationModel ( )

virtual

Definition at line 69 of file G4INCLStandardPropagationModel.cc.

    {
      delete theNucleus;
    }

Member Function Documentation

void G4INCL::StandardPropagationModel::generateAllAvatars

(

)

(Re)Generate all possible avatars.

Definition at line 433 of file G4INCLStandardPropagationModel.cc.

References generateCollisions(), generateDecays(), G4INCL::Store::getParticles(), getReflectionTime(), G4INCL::Nucleus::getStore(), and registerAvatar().

Referenced by shootComposite(), and shootParticle().

                                                      {
      ParticleList const &particles = theNucleus->getStore()->getParticles();
// assert(!particles.empty());
      for(ParticleIter i=particles.begin(), e=particles.end(); i!=e; ++i) {
        G4double time = this->getReflectionTime(*i);
        if(time <= maximumTime) registerAvatar(new SurfaceAvatar(*i, time, theNucleus));
      }
      generateCollisions(particles);
      generateDecays(particles);
    }

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

Generate a two-particle avatar.

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

Definition at line 261 of file G4INCLStandardPropagationModel.cc.

References G4INCL::AlwaysLocalEnergy, G4INCL::FirstCollisionLocalEnergy, G4INCL::Book::getAcceptedCollisions(), G4INCL::Store::getBook(), G4INCL::BinaryCollisionAvatar::getCutNNSquared(), G4INCL::Particle::getParticipantType(), G4INCL::Particle::getPosition(), G4INCL::Nucleus::getStore(), G4INCL::Nucleus::getSurfaceRadius(), getTime(), G4INCL::Particle::isNucleon(), G4INCL::Particle::isParticipant(), G4INCL::Particle::isPion(), G4INCL::Particle::isResonance(), G4INCL::ThreeVector::mag(), G4INCL::Particle::propagate(), G4INCL::KinematicsUtils::squareTotalEnergyInCM(), G4INCL::Math::tenPi, G4INCL::CrossSections::total(), and G4INCL::KinematicsUtils::transformToLocalEnergyFrame().

Referenced by generateCollisions(), and generateUpdatedCollisions().

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

      if(p1HasLocalEnergy) {
        backupParticle1 = *p1;
        p1->propagate(t - currentTime);
        if(p1->getPosition().mag() > theNucleus->getSurfaceRadius(p1)) {
          *p1 = backupParticle1;
          return NULL;
        }
        KinematicsUtils::transformToLocalEnergyFrame(theNucleus, p1);
      }
      if(p2HasLocalEnergy) {
        backupParticle2 = *p2;
        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::getCutNNSquared()) 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);
    }

void G4INCL::StandardPropagationModel::generateCollisions

(

const ParticleList &

particles

)

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.

Parameters

particles

list of particles

Definition at line 393 of file G4INCLStandardPropagationModel.cc.

References generateBinaryCollisionAvatar(), and registerAvatar().

Referenced by generateAllAvatars().

                                                                                   {
      // Loop over all the particles
      for(ParticleIter p1=particles.begin(), e=particles.end(); p1!=e; ++p1) {
        // Loop over the rest of the particles
        for(ParticleIter p2 = p1 + 1; p2 != particles.end(); ++p2) {
          registerAvatar(generateBinaryCollisionAvatar(*p1,*p2));
        }
      }
    }

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 403 of file G4INCLStandardPropagationModel.cc.

References generateBinaryCollisionAvatar(), and registerAvatar().

                                                                                                               {

      const G4bool haveExcept = !except.empty();

      // Loop over all the particles
      for(ParticleIter p1=particles.begin(), e=particles.end(); p1!=e; ++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));
        }
      }

    }

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 458 of file G4INCLStandardPropagationModel.cc.

References G4INCL::DeltaDecayChannel::computeDecayTime(), and registerAvatar().

Referenced by generateAllAvatars(), and propagate().

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

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 374 of file G4INCLStandardPropagationModel.cc.

References generateBinaryCollisionAvatar(), and registerAvatar().

Referenced by updateAvatars().

                                                                                                                                {

      // Loop over all the updated particles
      for(ParticleIter updated=updatedParticles.begin(), e=updatedParticles.end(); updated!=e; ++updated)
      {
        // Loop over all the particles
        for(ParticleIter particle=particles.begin(), end=particles.end(); particle!=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));
        }
      }
    }

G4double G4INCL::StandardPropagationModel::getCurrentTime ( )

virtual

Returns the current global time of the system.

Implements G4INCL::IPropagationModel.

Definition at line 247 of file G4INCLStandardPropagationModel.cc.

                                                      {
      return currentTime;
    }

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

virtual

Get the nucleus.

Implements G4INCL::IPropagationModel.

Definition at line 74 of file G4INCLStandardPropagationModel.cc.

    {
      return theNucleus;
    }

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 337 of file G4INCLStandardPropagationModel.cc.

References G4INCL::Intersection::exists, G4INCL::IntersectionFactory::getLaterTrajectoryIntersection(), G4INCL::Particle::getPosition(), G4INCL::Particle::getPropagationVelocity(), G4INCL::Nucleus::getSurfaceRadius(), INCL_ERROR, G4INCL::Particle::print(), and G4INCL::Intersection::time.

Referenced by generateAllAvatars(), and updateAvatars().

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

G4double G4INCL::StandardPropagationModel::getStoppingTime ( )

virtual

Get the current stopping time.

Implements G4INCL::IPropagationModel.

Definition at line 238 of file G4INCLStandardPropagationModel.cc.

                                                       {
      return maximumTime;
    }

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 354 of file G4INCLStandardPropagationModel.cc.

References G4INCL::ThreeVector::dot(), G4INCL::Particle::getPosition(), G4INCL::Particle::getPropagationVelocity(), and G4INCL::ThreeVector::mag2().

Referenced by generateBinaryCollisionAvatar().

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

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

virtual

Propagate all particles and return the first avatar.

Implements G4INCL::IPropagationModel.

Definition at line 470 of file G4INCLStandardPropagationModel.cc.

References G4INCL::Store::clearAvatars(), G4INCL::Store::findSmallestTime(), generateDecays(), G4INCL::Nucleus::getBlockedDelta(), G4INCL::Store::getBook(), G4INCL::Nucleus::getCreatedParticles(), G4INCL::Nucleus::getStore(), G4INCL::IAvatar::getTime(), G4INCL::Nucleus::getUpdatedParticles(), INCL_ERROR, G4INCL::Store::initialiseParticleAvatarConnections(), G4INCL::Book::setCurrentTime(), G4INCL::Store::timeStep(), and updateAvatars().

    {
      // 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();
        generateAllAvatarsExceptUpdated();
      }
#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 const &created = theNucleus->getCreatedParticles();
      needNewAvatars.insert(needNewAvatars.end(), created.begin(), created.end());
      updateAvatars(needNewAvatars);
#endif

      G4INCL::IAvatar *theAvatar = theNucleus->getStore()->findSmallestTime();
      if(theAvatar == 0) return 0; // Avatar list is empty
      //      theAvatar->dispose();

      if(theAvatar->getTime() < currentTime) {
        INCL_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;
    }

void G4INCL::StandardPropagationModel::registerAvatar

(

G4INCL::IAvatar *

anAvatar

)

Add an avatar to the storage.

Definition at line 256 of file G4INCLStandardPropagationModel.cc.

References G4INCL::Store::add(), and G4INCL::Nucleus::getStore().

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

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

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

virtual

Set the nucleus for this propagation model.

Implements G4INCL::IPropagationModel.

Definition at line 251 of file G4INCLStandardPropagationModel.cc.

    {
      theNucleus = nucleus;
    }

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

virtual

Set the stopping time of the simulation.

Implements G4INCL::IPropagationModel.

Definition at line 242 of file G4INCLStandardPropagationModel.cc.

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

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

virtual

Implements G4INCL::IPropagationModel.

Definition at line 79 of file G4INCLStandardPropagationModel.cc.

References G4INCL::Composite, shootComposite(), shootParticle(), and G4INCL::ParticleSpecies::theType.

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

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

virtual

Implements G4INCL::IPropagationModel.

Definition at line 160 of file G4INCLStandardPropagationModel.cc.

References G4INCL::Store::addParticleEntryAvatars(), G4INCL::Particle::boostVector(), G4INCL::CoulombDistortion::bringToSurface(), generateAllAvatars(), G4INCL::Particle::getA(), G4INCL::Cluster::getAngularMomentum(), G4INCL::Particle::getEnergy(), G4INCL::ParticleTable::getLargestNuclearRadius(), G4INCL::Particle::getMomentum(), G4INCL::Nucleus::getStore(), G4INCL::ParticleTable::getTableMass, G4INCL::Particle::getTransversePosition(), G4INCL::Nucleus::getUniverseRadius(), G4INCL::Particle::getZ(), INCL_DEBUG, G4INCL::ThreeVector::mag(), G4INCL::CoulombDistortion::maxImpactParameter(), position, gammaraytel::pr, G4INCL::Nucleus::setIncomingAngularMomentum(), G4INCL::Nucleus::setIncomingMomentum(), G4INCL::Nucleus::setInitialEnergy(), G4INCL::Nucleus::setNucleusNucleusCollision(), G4INCL::Cluster::setPosition(), G4INCL::Nucleus::setProjectileChargeNumber(), G4INCL::Nucleus::setProjectileMassNumber(), G4INCL::Nucleus::setProjectileRemnant(), and G4INCL::ProjectileRemnant::storeComponents().

Referenced by shoot().

                                                                                                                                                                     {
      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::getLargestNuclearRadius(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;
      INCL_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)) {
        INCL_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);

      // 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()) {
        INCL_DEBUG("No ParticleEntryAvatar found, transparent event" << std::endl);
        delete pr;
        return -1.;
      }

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

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

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

virtual

Implements G4INCL::IPropagationModel.

Definition at line 86 of file G4INCLStandardPropagationModel.cc.

References G4INCL::Store::addParticleEntryAvatar(), G4INCL::Particle::adjustMomentumFromEnergy(), G4INCL::Particle::boostVector(), G4INCL::CoulombDistortion::bringToSurface(), energy(), generateAllAvatars(), G4INCL::Particle::getA(), G4INCL::Particle::getAngularMomentum(), G4INCL::Particle::getEnergy(), G4INCL::Particle::getMass(), G4INCL::Particle::getMomentum(), G4INCL::Particle::getSpecies(), G4INCL::Nucleus::getStore(), G4INCL::ParticleTable::getTableMass, G4INCL::ParticleTable::getTableParticleMass, G4INCL::Particle::getTransversePosition(), G4INCL::Nucleus::getUniverseRadius(), G4INCL::Particle::getZ(), INCL_DEBUG, G4INCL::Particle::isNucleon(), G4INCL::ThreeVector::mag(), G4INCL::Particle::makeProjectileSpectator(), G4INCL::CoulombDistortion::maxImpactParameter(), position, G4INCL::Particle::setEnergy(), G4INCL::Particle::setINCLMass(), G4INCL::Nucleus::setIncomingAngularMomentum(), G4INCL::Nucleus::setIncomingMomentum(), G4INCL::Nucleus::setInitialEnergy(), G4INCL::Nucleus::setParticleNucleusCollision(), G4INCL::Particle::setPosition(), G4INCL::Nucleus::setProjectileChargeNumber(), and G4INCL::Nucleus::setProjectileMassNumber().

Referenced by shoot().

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

void G4INCL::StandardPropagationModel::updateAvatars

(

const ParticleList &

particles

)

Update all avatars related to a particle.

Definition at line 423 of file G4INCLStandardPropagationModel.cc.

References generateUpdatedCollisions(), G4INCL::Store::getParticles(), getReflectionTime(), G4INCL::Nucleus::getStore(), and registerAvatar().

Referenced by propagate().

                                                                              {

      for(ParticleIter iter=particles.begin(), e=particles.end(); iter!=e; ++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
    }

---

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

• G4INCLStandardPropagationModel.hh
• G4INCLStandardPropagationModel.cc