#include <G4INCLStandardPropagationModel.hh>

Inheritance diagram for G4INCL::StandardPropagationModel:

Public Member Functions
void	generateAllAvatars ()
	(Re)Generate all possible avatars. More...

IAvatar *	generateBinaryCollisionAvatar (Particle const p1, Particle const p2)
	Generate a two-particle avatar. 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	generateUpdatedCollisions (const ParticleList &updatedParticles, const ParticleList &particles)
	Generate and register collisions between a list of updated particles and all the other particles. More...

G4double	getCurrentTime ()

G4INCL::Nucleus *	getNucleus ()

G4double	getReflectionTime (G4INCL::Particle const *const aParticle)
	Get the reflection time. More...

G4double	getStoppingTime ()

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

G4INCL::IAvatar *	propagate (FinalState const *const fs)

void	registerAvatar (G4INCL::IAvatar *anAvatar)

void	setNucleus (G4INCL::Nucleus *nucleus)

void	setStoppingTime (G4double)

G4double	shoot (ParticleSpecies const &projectileSpecies, const G4double kineticEnergy, const G4double impactParameter, const G4double phi)

G4double	shootComposite (ParticleSpecies const &s, const G4double kineticEnergy, const G4double impactParameter, const G4double phi)

G4double	shootParticle (ParticleType const t, const G4double kineticEnergy, const G4double impactParameter, const G4double phi)

	StandardPropagationModel (LocalEnergyType localEnergyType, LocalEnergyType localEnergyDeltaType, const G4double hTime=0.0)

void	updateAvatars (const ParticleList &particles)

virtual	~StandardPropagationModel ()

Private Attributes
Particle	backupParticle1

Particle	backupParticle2

G4double	currentTime

G4bool	firstAvatar

G4double	hadronizationTime

G4double	maximumTime

LocalEnergyType	theLocalEnergyDeltaType

LocalEnergyType	theLocalEnergyType

G4INCL::Nucleus *	theNucleus

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

Constructor & Destructor Documentation

◆ StandardPropagationModel()

G4INCL::StandardPropagationModel::StandardPropagationModel	(	LocalEnergyType	localEnergyType,
		LocalEnergyType	localEnergyDeltaType,
		const G4double	hTime = `0.0`
	)

Definition at line 65 of file G4INCLStandardPropagationModel.cc.

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

◆ ~StandardPropagationModel()

G4INCL::StandardPropagationModel::~StandardPropagationModel ( )

virtual

Definition at line 74 of file G4INCLStandardPropagationModel.cc.

    {
      delete theNucleus;
    }

References theNucleus.

Member Function Documentation

◆ generateAllAvatars()

void G4INCL::StandardPropagationModel::generateAllAvatars ( )

(Re)Generate all possible avatars.

Definition at line 437 of file G4INCLStandardPropagationModel.cc.

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

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

Referenced by shootComposite(), and shootParticle().

◆ generateBinaryCollisionAvatar()

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 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+hadronizationTime) 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->isMeson());
      const G4bool p2HasLocalEnergy = (hasLocalEnergy && !p2->isMeson());
 
      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);
    }

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

Referenced by generateCollisions(), and generateUpdatedCollisions().

◆ generateCollisions() [1/2]

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

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

References generateBinaryCollisionAvatar(), and registerAvatar().

Referenced by generateAllAvatars().

◆ generateCollisions() [2/2]

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

                                                                                                               {
 
      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 && except.contains(*p1) && except.contains(*p2)) continue;
 
          registerAvatar(generateBinaryCollisionAvatar(*p1,*p2));
        }
      }
 
    }

References G4INCL::UnorderedVector< T >::contains(), generateBinaryCollisionAvatar(), and registerAvatar().

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

                                                                               {
      for(ParticleIter i=particles.begin(), e=particles.end(); i!=e; ++i) {
           if((*i)->isDelta()) {
             G4double decayTime = DeltaDecayChannel::computeDecayTime((*i)); // time in fm/c
             G4double time = currentTime + decayTime;
             if(time <= maximumTime) {
               registerAvatar(new DecayAvatar((*i), time, theNucleus));
             }
           }
           else if((*i)->getType() == SigmaZero) {
             G4double decayTime = SigmaZeroDecayChannel::computeDecayTime((*i)); // time in fm/c
             G4double time = currentTime + decayTime;
             if(time <= maximumTime) {
               registerAvatar(new DecayAvatar((*i), time, theNucleus));
             }
           }
        if((*i)->isOmega()) {
          G4double decayTimeOmega = PionResonanceDecayChannel::computeDecayTime((*i));
          G4double timeOmega = currentTime + decayTimeOmega;
          if(timeOmega <= maximumTime) {
            registerAvatar(new DecayAvatar((*i), timeOmega, theNucleus));
          }
        }
      }
    }

References G4INCL::DeltaDecayChannel::computeDecayTime(), G4INCL::PionResonanceDecayChannel::computeDecayTime(), G4INCL::SigmaZeroDecayChannel::computeDecayTime(), currentTime, maximumTime, registerAvatar(), G4INCL::SigmaZero, and 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(), 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(updatedParticles.contains(*particle)) continue;
 
          registerAvatar(generateBinaryCollisionAvatar(*particle,*updated));
        }
      }
    }

References G4INCL::UnorderedVector< T >::contains(), generateBinaryCollisionAvatar(), and registerAvatar().

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

References currentTime.

◆ getNucleus()

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

virtual

Get the nucleus.

Implements G4INCL::IPropagationModel.

Definition at line 79 of file G4INCLStandardPropagationModel.cc.

    {
      return theNucleus;
    }

References 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 {
        INCL_ERROR("Imaginary reflection time for particle: " << '\n'
          << aParticle->print());
        time = 10000.0;
      }
      return time;
    }

References currentTime, G4INCL::Intersection::exists, G4INCL::IntersectionFactory::getLaterTrajectoryIntersection(), G4INCL::Particle::getPosition(), G4INCL::Particle::getPropagationVelocity(), G4INCL::Nucleus::getSurfaceRadius(), INCL_ERROR, G4INCL::Particle::print(), theNucleus, and G4INCL::Intersection::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;
    }

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

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

Referenced by generateBinaryCollisionAvatar().

◆ propagate()

G4INCL::IAvatar * G4INCL::StandardPropagationModel::propagate ( FinalState const *const fs )

virtual

Propagate all particles and return the first avatar.

Implements G4INCL::IPropagationModel.

Definition at line 491 of file G4INCLStandardPropagationModel.cc.

    {
      if(fs) {
        // 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(!fs->getModifiedParticles().empty() || !fs->getEnteringParticles().empty() || !fs->getCreatedParticles().empty()) {
          // Regenerates the entire avatar list, skipping collisions between
          // updated particles
          theNucleus->getStore()->clearAvatars();
          theNucleus->getStore()->initialiseParticleAvatarConnections();
          generateAllAvatarsExceptUpdated(fs);
        }
#else
        ParticleList const &updatedParticles = fs->getModifiedParticles();
        if(fs->getValidity()==PauliBlockedFS) {
          // This final state might represents the outcome of a Pauli-blocked delta
          // decay
// assert(updatedParticles.empty() || (updatedParticles.size()==1 && updatedParticles.front()->isResonance()));
// assert(fs->getEnteringParticles().empty() && fs->getCreatedParticles().empty() && fs->getOutgoingParticles().empty() && fs->getDestroyedParticles().empty());
          generateDecays(updatedParticles);
        } else {
          ParticleList const &entering = fs->getEnteringParticles();
          generateDecays(updatedParticles);
          generateDecays(entering);
 
          ParticleList const &created = fs->getCreatedParticles();
          if(created.empty() && entering.empty())
            updateAvatars(updatedParticles);
          else {
            ParticleList updatedParticlesCopy = updatedParticles;
            updatedParticlesCopy.insert(updatedParticlesCopy.end(), entering.begin(), entering.end());
            updatedParticlesCopy.insert(updatedParticlesCopy.end(), created.begin(), created.end());
            updateAvatars(updatedParticlesCopy);
          }
        }
#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 << '\n');
        return 0;
      } else if(theAvatar->getTime() > currentTime) {
        theNucleus->getStore()->timeStep(theAvatar->getTime() - currentTime);
 
        currentTime = theAvatar->getTime();
        theNucleus->getStore()->getBook().setCurrentTime(currentTime);
      }
 
      return theAvatar;
    }

References G4INCL::Store::clearAvatars(), currentTime, G4INCL::Store::findSmallestTime(), generateDecays(), G4INCL::Store::getBook(), G4INCL::FinalState::getCreatedParticles(), G4INCL::FinalState::getEnteringParticles(), G4INCL::FinalState::getModifiedParticles(), G4INCL::Nucleus::getStore(), G4INCL::IAvatar::getTime(), G4INCL::FinalState::getValidity(), INCL_ERROR, G4INCL::PauliBlockedFS, G4INCL::Book::setCurrentTime(), theNucleus, G4INCL::Store::timeStep(), and updateAvatars().

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

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

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

References theNucleus.

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

References maximumTime.

◆ shoot()

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

virtual

Implements G4INCL::IPropagationModel.

Definition at line 84 of file G4INCLStandardPropagationModel.cc.

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

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

◆ shootComposite()

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

virtual

Implements G4INCL::IPropagationModel.

Definition at line 168 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::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 << '\n');
 
      // 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" << '\n');
        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(),theNucleus->getS()));
 
      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" << '\n');
        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);
 
      // Register the ParticleEntryAvatars
      theNucleus->getStore()->addParticleEntryAvatars(theAvatarList);
 
      return pr->getTransversePosition().mag();
    }

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

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 91 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;
      G4double TLab;
      if( p->isMeson()) {
        temfin = 30.18 * std::pow(theNucleus->getA(), 0.17);
        TLab = p->getKineticEnergy();
      } else {
        temfin = 29.8 * std::pow(theNucleus->getA(), 0.16);
        TLab = p->getKineticEnergy()/p->getA();
      }
 
      // energy-dependent stopping time above 2 AGeV
      if(TLab>2000.)
        temfin *= (5.8E4-TLab)/5.6E4;
 
      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 << '\n');
 
      // 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" << '\n');
        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(),theNucleus->getS()));
 
      // 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);
 
        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 427 of file G4INCLStandardPropagationModel.cc.

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

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

Referenced by propagate().

Field Documentation

◆ backupParticle1

Particle G4INCL::StandardPropagationModel::backupParticle1

private

Definition at line 200 of file G4INCLStandardPropagationModel.hh.

Referenced by generateBinaryCollisionAvatar().

◆ backupParticle2

Particle G4INCL::StandardPropagationModel::backupParticle2

private

Definition at line 200 of file G4INCLStandardPropagationModel.hh.

Referenced by generateBinaryCollisionAvatar().

◆ currentTime

G4double G4INCL::StandardPropagationModel::currentTime

private

Definition at line 196 of file G4INCLStandardPropagationModel.hh.

Referenced by generateBinaryCollisionAvatar(), generateDecays(), getCurrentTime(), getReflectionTime(), getTime(), propagate(), shootComposite(), and shootParticle().

◆ firstAvatar

G4bool G4INCL::StandardPropagationModel::firstAvatar

private

Definition at line 198 of file G4INCLStandardPropagationModel.hh.

Referenced by shootComposite(), and shootParticle().

◆ hadronizationTime

G4double G4INCL::StandardPropagationModel::hadronizationTime

private

Definition at line 197 of file G4INCLStandardPropagationModel.hh.

Referenced by generateBinaryCollisionAvatar().

◆ maximumTime

G4double G4INCL::StandardPropagationModel::maximumTime

private

Definition at line 195 of file G4INCLStandardPropagationModel.hh.

Referenced by generateAllAvatars(), generateBinaryCollisionAvatar(), generateDecays(), getStoppingTime(), setStoppingTime(), shootComposite(), shootParticle(), and updateAvatars().

◆ theLocalEnergyDeltaType

LocalEnergyType G4INCL::StandardPropagationModel::theLocalEnergyDeltaType

private

Definition at line 199 of file G4INCLStandardPropagationModel.hh.

Referenced by generateBinaryCollisionAvatar().

◆ theLocalEnergyType

LocalEnergyType G4INCL::StandardPropagationModel::theLocalEnergyType

private

Definition at line 199 of file G4INCLStandardPropagationModel.hh.

Referenced by generateBinaryCollisionAvatar().

◆ theNucleus

G4INCL::Nucleus* G4INCL::StandardPropagationModel::theNucleus

private

Definition at line 194 of file G4INCLStandardPropagationModel.hh.

Referenced by generateAllAvatars(), generateBinaryCollisionAvatar(), generateDecays(), getNucleus(), getReflectionTime(), propagate(), registerAvatar(), setNucleus(), shootComposite(), shootParticle(), updateAvatars(), and ~StandardPropagationModel().

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

source/processes/hadronic/models/inclxx/incl_physics/include/G4INCLStandardPropagationModel.hh
source/processes/hadronic/models/inclxx/incl_physics/src/G4INCLStandardPropagationModel.cc

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ StandardPropagationModel()

◆ ~StandardPropagationModel()

Member Function Documentation

◆ generateAllAvatars()

◆ generateBinaryCollisionAvatar()

◆ generateCollisions() [1/2]

◆ generateCollisions() [2/2]

◆ generateDecays()

◆ generateUpdatedCollisions()

◆ getCurrentTime()

◆ getNucleus()

◆ getReflectionTime()

◆ getStoppingTime()

◆ getTime()

◆ propagate()

◆ registerAvatar()

◆ setNucleus()

◆ setStoppingTime()

◆ shoot()

◆ shootComposite()

◆ shootParticle()

◆ updateAvatars()

Field Documentation

◆ backupParticle1

◆ backupParticle2

◆ currentTime

◆ firstAvatar

◆ hadronizationTime

◆ maximumTime

◆ theLocalEnergyDeltaType

◆ theLocalEnergyType

◆ theNucleus