#include <G4INCLNucleus.hh>

Inheritance diagram for G4INCL::Nucleus:

Collaboration diagram for G4INCL::Nucleus:

Classes
struct	ConservationBalance
	Struct for conservation laws. More...

Public Member Functions
	Nucleus (G4int mass, G4int charge, Config const *const conf, const G4double universeRadius=-1.)

virtual	~Nucleus ()

	Nucleus (const Nucleus &rhs)
	Dummy copy constructor to silence Coverity warning. More...

Nucleus &	operator= (const Nucleus &rhs)
	Dummy assignment operator to silence Coverity warning. More...

void	initializeParticles ()

void	insertParticle (Particle *p)
	Insert a new particle (e.g. a projectile) in the nucleus. More...

void	applyFinalState (FinalState *)

G4int	getInitialA () const

G4int	getInitialZ () const

void	propagateParticles (G4double step)

G4int	getNumberOfEnteringProtons () const

G4int	getNumberOfEnteringNeutrons () const

G4double	computeSeparationEnergyBalance () const
	Outgoing - incoming separation energies. More...

G4bool	decayOutgoingDeltas ()
	Force the decay of outgoing deltas. More...

G4bool	decayInsideDeltas ()
	Force the decay of deltas inside the nucleus. More...

G4bool	decayOutgoingClusters ()
	Force the decay of unstable outgoing clusters. More...

G4bool	decayMe ()
	Force the phase-space decay of the Nucleus. More...

void	emitInsidePions ()
	Force emission of all pions inside the nucleus. More...

void	computeRecoilKinematics ()
	Compute the recoil momentum and spin of the nucleus. More...

ThreeVector	computeCenterOfMass () const
	Compute the current center-of-mass position. More...

G4double	computeTotalEnergy () const
	Compute the current total energy. More...

G4double	computeExcitationEnergy () const
	Compute the current excitation energy. More...

void	setIncomingAngularMomentum (const ThreeVector &j)
	Set the incoming angular-momentum vector. More...

const ThreeVector &	getIncomingAngularMomentum () const
	Get the incoming angular-momentum vector. More...

void	setIncomingMomentum (const ThreeVector &p)
	Set the incoming momentum vector. More...

const ThreeVector &	getIncomingMomentum () const
	Get the incoming momentum vector. More...

void	setInitialEnergy (const G4double e)
	Set the initial energy. More...

G4double	getInitialEnergy () const
	Get the initial energy. More...

G4double	getExcitationEnergy () const
	Get the excitation energy of the nucleus. More...

G4bool	containsDeltas ()
	Returns true if the nucleus contains any deltas. More...

std::string	print ()

Store *	getStore () const

void	setStore (Store *s)

G4double	getInitialInternalEnergy () const

G4bool	isEventTransparent () const
	Is the event transparent? More...

G4bool	hasRemnant () const
	Does the nucleus give a cascade remnant? More...

void	fillEventInfo (EventInfo *eventInfo)

G4bool	getTryCompoundNucleus ()

G4double	getTransmissionBarrier (Particle const *const p)
	Get the transmission barrier. More...

ConservationBalance	getConservationBalance (EventInfo const &theEventInfo, const G4bool afterRecoil) const
	Compute charge, mass, energy and momentum balance. More...

void	useFusionKinematics ()
	Adjust the kinematics for complete-fusion events. More...

G4double	getSurfaceRadius (Particle const *const particle) const
	Get the maximum allowed radius for a given particle. More...

G4double	getUniverseRadius () const
	Getter for theUniverseRadius. More...

void	setUniverseRadius (const G4double universeRadius)
	Setter for theUniverseRadius. More...

G4bool	isNucleusNucleusCollision () const
	Is it a nucleus-nucleus collision? More...

void	setNucleusNucleusCollision ()
	Set a nucleus-nucleus collision. More...

void	setParticleNucleusCollision ()
	Set a particle-nucleus collision. More...

void	setProjectileRemnant (ProjectileRemnant *const c)
	Set the projectile remnant. More...

ProjectileRemnant *	getProjectileRemnant () const
	Get the projectile remnant. More...

void	deleteProjectileRemnant ()
	Delete the projectile remnant. More...

void	finalizeProjectileRemnant (const G4double emissionTime)
	Finalise the projectile remnant. More...

void	updatePotentialEnergy (Particle *p) const
	Update the particle potential energy. More...

void	setDensity (NuclearDensity const *const d)
	Setter for theDensity. More...

NuclearDensity const *	getDensity () const
	Getter for theDensity. More...

NuclearPotential::INuclearPotential const *	getPotential () const
	Getter for thePotential. More...

Public Member Functions inherited from G4INCL::Cluster
	Cluster (const G4int Z, const G4int A, const G4bool createParticleSampler=true)
	Standard Cluster constructor. More...

template<class Iterator >
	Cluster (Iterator begin, Iterator end)

virtual	~Cluster ()

	Cluster (const Cluster &rhs)
	Copy constructor. More...

Cluster &	operator= (const Cluster &rhs)
	Assignment operator. More...

void	swap (Cluster &rhs)
	Helper method for the assignment operator. More...

ParticleSpecies	getSpecies () const

void	deleteParticles ()

void	clearParticles ()

void	setZ (const G4int Z)
	Set the charge number of the cluster. More...

void	setA (const G4int A)
	Set the mass number of the cluster. More...

G4double	getExcitationEnergy () const
	Get the excitation energy of the cluster. More...

void	setExcitationEnergy (const G4double e)
	Set the excitation energy of the cluster. More...

virtual G4double	getTableMass () const
	Get the real particle mass. More...

ParticleList const &	getParticles () const

void	removeParticle (Particle *const p)
	Remove a particle from the cluster components. More...

void	addParticle (Particle *const p)

void	updateClusterParameters ()
	Set total cluster mass, energy, size, etc. from the particles. More...

void	addParticles (ParticleList const &pL)
	Add a list of particles to the cluster. More...

ParticleList	getParticleList () const
	Returns the list of particles that make up the cluster. More...

std::string	print () const

void	internalBoostToCM ()
	Boost to the CM of the component particles. More...

void	putParticlesOffShell ()
	Put the cluster components off shell. More...

void	setPosition (const ThreeVector &position)
	Set the position of the cluster. More...

void	boost (const ThreeVector &aBoostVector)
	Boost the cluster with the indicated velocity. More...

void	freezeInternalMotion ()
	Freeze the internal motion of the particles. More...

virtual void	rotatePosition (const G4double angle, const ThreeVector &axis)
	Rotate position of all the particles. More...

virtual void	rotateMomentum (const G4double angle, const ThreeVector &axis)
	Rotate momentum of all the particles. More...

virtual void	makeProjectileSpectator ()
	Make all the components projectile spectators, too. More...

virtual void	makeTargetSpectator ()
	Make all the components target spectators, too. More...

virtual void	makeParticipant ()
	Make all the components participants, too. More...

ThreeVector const &	getSpin () const
	Get the spin of the nucleus. More...

void	setSpin (const ThreeVector &j)
	Set the spin of the nucleus. More...

G4INCL::ThreeVector	getAngularMomentum () const
	Get the total angular momentum (orbital + spin) More...

Public Member Functions inherited from G4INCL::Particle
	Particle ()

	Particle (ParticleType t, G4double energy, ThreeVector const &momentum, ThreeVector const &position)

	Particle (ParticleType t, ThreeVector const &momentum, ThreeVector const &position)

virtual	~Particle ()

	Particle (const Particle &rhs)
	Copy constructor. More...

Particle &	operator= (const Particle &rhs)
	Assignment operator. More...

G4INCL::ParticleType	getType () const

void	setType (ParticleType t)

G4bool	isNucleon () const

ParticipantType	getParticipantType () const

void	setParticipantType (ParticipantType const p)

G4bool	isParticipant () const

G4bool	isTargetSpectator () const

G4bool	isProjectileSpectator () const

G4bool	isPion () const
	Is this a pion? More...

G4bool	isResonance () const
	Is it a resonance? More...

G4bool	isDelta () const
	Is it a Delta? More...

G4int	getA () const
	Returns the baryon number. More...

G4int	getZ () const
	Returns the charge number. More...

G4double	getBeta () const

ThreeVector	boostVector () const

void	boost (const ThreeVector &aBoostVector)

void	lorentzContract (const ThreeVector &aBoostVector, const ThreeVector &refPos)
	Lorentz-contract the particle position around some center. More...

G4double	getMass () const
	Get the cached particle mass. More...

G4double	getINCLMass () const
	Get the INCL particle mass. More...

G4double	getRealMass () const
	Get the real particle mass. More...

void	setRealMass ()
	Set the mass of the Particle to its real mass. More...

void	setTableMass ()
	Set the mass of the Particle to its table mass. More...

void	setINCLMass ()
	Set the mass of the Particle to its table mass. More...

G4double	getEmissionQValueCorrection (const G4int AParent, const G4int ZParent) const
	Computes correction on the emission Q-value. More...

G4double	getTransferQValueCorrection (const G4int AFrom, const G4int ZFrom, const G4int ATo, const G4int ZTo) const
	Computes correction on the transfer Q-value. More...

G4double	getInvariantMass () const
	Get the the particle invariant mass. More...

G4double	getKineticEnergy () const
	Get the particle kinetic energy. More...

G4double	getPotentialEnergy () const
	Get the particle potential energy. More...

void	setPotentialEnergy (G4double v)
	Set the particle potential energy. More...

G4double	getEnergy () const

void	setMass (G4double mass)

void	setEnergy (G4double energy)

const G4INCL::ThreeVector &	getMomentum () const

virtual void	setMomentum (const G4INCL::ThreeVector &momentum)

const G4INCL::ThreeVector &	getPosition () const

G4double	getHelicity ()

void	setHelicity (G4double h)

void	propagate (G4double step)

G4int	getNumberOfCollisions () const
	Return the number of collisions undergone by the particle. More...

void	setNumberOfCollisions (G4int n)
	Set the number of collisions undergone by the particle. More...

void	incrementNumberOfCollisions ()
	Increment the number of collisions undergone by the particle. More...

G4int	getNumberOfDecays () const
	Return the number of decays undergone by the particle. More...

void	setNumberOfDecays (G4int n)
	Set the number of decays undergone by the particle. More...

void	incrementNumberOfDecays ()
	Increment the number of decays undergone by the particle. More...

void	setOutOfWell ()
	Mark the particle as out of its potential well. More...

G4bool	isOutOfWell () const
	Check if the particle is out of its potential well. More...

void	setEmissionTime (G4double t)

G4double	getEmissionTime ()

ThreeVector	getTransversePosition () const
	Transverse component of the position w.r.t. the momentum. More...

ThreeVector	getLongitudinalPosition () const
	Longitudinal component of the position w.r.t. the momentum. More...

const ThreeVector &	adjustMomentumFromEnergy ()
	Rescale the momentum to match the total energy. More...

G4double	adjustEnergyFromMomentum ()
	Recompute the energy to match the momentum. More...

G4bool	isCluster () const

void	setFrozenMomentum (const ThreeVector &momentum)
	Set the frozen particle momentum. More...

void	setFrozenEnergy (const G4double energy)
	Set the frozen particle momentum. More...

ThreeVector	getFrozenMomentum () const
	Get the frozen particle momentum. More...

G4double	getFrozenEnergy () const
	Get the frozen particle momentum. More...

ThreeVector	getPropagationVelocity () const
	Get the propagation velocity of the particle. More...

void	freezePropagation ()
	Freeze particle propagation. More...

void	thawPropagation ()
	Unfreeze particle propagation. More...

virtual void	rotatePositionAndMomentum (const G4double angle, const ThreeVector &axis)
	Rotate the particle position and momentum. More...

std::string	print () const

std::string	dump () const

long	getID () const

ParticleList const *	getParticles () const

G4double	getReflectionMomentum () const
	Return the reflection momentum. More...

void	setUncorrelatedMomentum (const G4double p)
	Set the uncorrelated momentum. More...

void	rpCorrelate ()
	Make the particle follow a strict r-p correlation. More...

void	rpDecorrelate ()
	Make the particle not follow a strict r-p correlation. More...

G4double	getCosRPAngle () const
	Get the cosine of the angle between position and momentum. More...

Private Member Functions
void	computeOneNucleonRecoilKinematics ()
	Compute the recoil kinematics for a 1-nucleon remnant. More...

	INCL_DECLARE_ALLOCATION_POOL (Nucleus)

Private Attributes
G4int	theInitialZ

G4int	theInitialA

G4int	theNpInitial
	The number of entering protons. More...

G4int	theNnInitial
	The number of entering neutrons. More...

G4double	initialInternalEnergy

ThreeVector	incomingAngularMomentum

ThreeVector	incomingMomentum

ThreeVector	initialCenterOfMass

G4bool	remnant

G4double	initialEnergy

Store *	theStore

G4bool	tryCN

G4int	projectileZ
	The charge number of the projectile. More...

G4int	projectileA
	The mass number of the projectile. More...

G4double	theUniverseRadius
	The radius of the universe. More...

G4bool	isNucleusNucleus
	true if running a nucleus-nucleus collision More...

ProjectileRemnant *	theProjectileRemnant
	Pointer to the quasi-projectile. More...

NuclearDensity const *	theDensity
	Pointer to the NuclearDensity object. More...

NuclearPotential::INuclearPotential const *	thePotential
	Pointer to the NuclearPotential object. More...

Additional Inherited Members
Protected Member Functions inherited from G4INCL::Cluster
	INCL_DECLARE_ALLOCATION_POOL (Cluster)

Protected Member Functions inherited from G4INCL::Particle
void	swap (Particle &rhs)
	Helper method for the assignment operator. More...

Protected Attributes inherited from G4INCL::Cluster
ParticleList	particles

G4double	theExcitationEnergy

ThreeVector	theSpin

ParticleSampler *	theParticleSampler

Protected Attributes inherited from G4INCL::Particle
G4int	theZ

G4int	theA

ParticipantType	theParticipantType

G4INCL::ParticleType	theType

G4double	theEnergy

G4double *	thePropagationEnergy

G4double	theFrozenEnergy

G4INCL::ThreeVector	theMomentum

G4INCL::ThreeVector *	thePropagationMomentum

G4INCL::ThreeVector	theFrozenMomentum

G4INCL::ThreeVector	thePosition

G4int	nCollisions

G4int	nDecays

G4double	thePotentialEnergy

long	ID

G4bool	rpCorrelated

G4double	uncorrelatedMomentum

Detailed Description

Definition at line 65 of file G4INCLNucleus.hh.

Constructor & Destructor Documentation

◆ Nucleus() [1/2]

G4INCL::Nucleus::Nucleus	(	G4int	mass,
		G4int	charge,
		Config const *const	conf,
		const G4double	universeRadius = `-1.`
	)

Definition at line 65 of file G4INCLNucleus.cc.

     : Cluster(charge,mass,true),
      theInitialZ(charge), theInitialA(mass),
      theNpInitial(0), theNnInitial(0),
      initialInternalEnergy(0.),
      incomingAngularMomentum(0.,0.,0.), incomingMomentum(0.,0.,0.),
      initialCenterOfMass(0.,0.,0.),
      remnant(true),
      initialEnergy(0.),
      tryCN(false),
      theUniverseRadius(universeRadius),
      isNucleusNucleus(false),
      theProjectileRemnant(NULL),
      theDensity(NULL),
      thePotential(NULL)
   {
     PotentialType potentialType;
     G4bool pionPotential;
     if(conf) {
       potentialType = conf->getPotentialType();
       pionPotential = conf->getPionPotential();
     } else { // By default we don't use energy dependent
       // potential. This is convenient for some tests.
       potentialType = IsospinPotential;
       pionPotential = true;
     }
 
     thePotential = NuclearPotential::createPotential(potentialType, theA, theZ, pionPotential);
 
     ParticleTable::setProtonSeparationEnergy(thePotential->getSeparationEnergy(Proton));
     ParticleTable::setNeutronSeparationEnergy(thePotential->getSeparationEnergy(Neutron));
 
     theDensity = NuclearDensityFactory::createDensity(theA, theZ);
 
     theParticleSampler->setPotential(thePotential);
     theParticleSampler->setDensity(theDensity);
 
     if(theUniverseRadius<0)
       theUniverseRadius = theDensity->getMaximumRadius();
     theStore = new Store(conf);
   }

Here is the call graph for this function:

◆ ~Nucleus()

G4INCL::Nucleus::~Nucleus ( )

virtual

Definition at line 107 of file G4INCLNucleus.cc.

                     {
     delete theStore;
     deleteProjectileRemnant();
     /* We don't delete the potential and the density here any more -- Factories
      * are caching them
     delete thePotential;
     delete theDensity;*/
   }

Here is the call graph for this function:

◆ Nucleus() [2/2]

G4INCL::Nucleus::Nucleus ( const Nucleus & rhs )

Dummy copy constructor to silence Coverity warning.

Member Function Documentation

◆ applyFinalState()

void G4INCL::Nucleus::applyFinalState ( FinalState * finalstate )

Apply reaction final state information to the nucleus.

Definition at line 131 of file G4INCLNucleus.cc.

                                                       {
     if(!finalstate) // do nothing if no final state was returned
       return;
 
     G4double totalEnergy = 0.0;
 
     FinalStateValidity const validity = finalstate->getValidity();
     if(validity == ValidFS) {
 
       ParticleList const &created = finalstate->getCreatedParticles();
       for(ParticleIter iter=created.begin(), e=created.end(); iter!=e; ++iter) {
         theStore->add((*iter));
         if(!(*iter)->isOutOfWell()) {
           totalEnergy += (*iter)->getEnergy() - (*iter)->getPotentialEnergy();
         }
       }
 
       ParticleList const &deleted = finalstate->getDestroyedParticles();
       for(ParticleIter iter=deleted.begin(), e=deleted.end(); iter!=e; ++iter) {
         theStore->particleHasBeenDestroyed(*iter);
       }
 
       ParticleList const &modified = finalstate->getModifiedParticles();
       for(ParticleIter iter=modified.begin(), e=modified.end(); iter!=e; ++iter) {
         theStore->particleHasBeenUpdated(*iter);
         totalEnergy += (*iter)->getEnergy() - (*iter)->getPotentialEnergy();
       }
 
       ParticleList const &out = finalstate->getOutgoingParticles();
       for(ParticleIter iter=out.begin(), e=out.end(); iter!=e; ++iter) {
         if((*iter)->isCluster()) {
       Cluster *clusterOut = dynamic_cast<Cluster*>((*iter));
 // assert(clusterOut);
 #ifdef INCLXX_IN_GEANT4_MODE
           if(!clusterOut)
             continue;
 #endif
           ParticleList const &components = clusterOut->getParticles();
           for(ParticleIter in=components.begin(), end=components.end(); in!=end; ++in)
             theStore->particleHasBeenEjected(*in);
         } else {
           theStore->particleHasBeenEjected(*iter);
         }
         totalEnergy += (*iter)->getEnergy();      // No potential here because the particle is gone
         theA -= (*iter)->getA();
         theZ -= (*iter)->getZ();
         theStore->addToOutgoing(*iter);
         (*iter)->setEmissionTime(theStore->getBook().getCurrentTime());
       }
 
       ParticleList const &entering = finalstate->getEnteringParticles();
       for(ParticleIter iter=entering.begin(), e=entering.end(); iter!=e; ++iter) {
         insertParticle(*iter);
         totalEnergy += (*iter)->getEnergy() - (*iter)->getPotentialEnergy();
       }
 
       // actually perform the removal of the scheduled avatars
       theStore->removeScheduledAvatars();
     } else if(validity == ParticleBelowFermiFS || validity == ParticleBelowZeroFS) {
       INCL_DEBUG("A Particle is entering below the Fermi sea:" << '\n' << finalstate->print() << '\n');
       tryCN = true;
       ParticleList const &entering = finalstate->getEnteringParticles();
       for(ParticleIter iter=entering.begin(), e=entering.end(); iter!=e; ++iter) {
         insertParticle(*iter);
       }
     }
 
     if(validity==ValidFS &&
         std::abs(totalEnergy - finalstate->getTotalEnergyBeforeInteraction()) > 0.1) {
       INCL_ERROR("Energy nonconservation! Energy at the beginning of the event = "
           << finalstate->getTotalEnergyBeforeInteraction()
           <<" and after interaction = "
           << totalEnergy << '\n'
           << finalstate->print());
     }
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ computeCenterOfMass()

ThreeVector G4INCL::Nucleus::computeCenterOfMass ( ) const

Compute the current center-of-mass position.

Returns: the center-of-mass position vector [fm].

Definition at line 259 of file G4INCLNucleus.cc.

                                                  {
     ThreeVector cm(0.,0.,0.);
     G4double totalMass = 0.0;
     ParticleList const &inside = theStore->getParticles();
     for(ParticleIter p=inside.begin(), e=inside.end(); p!=e; ++p) {
       const G4double mass = (*p)->getMass();
       cm += (*p)->getPosition() * mass;
       totalMass += mass;
     }
     cm /= totalMass;
     return cm;
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ computeExcitationEnergy()

G4double G4INCL::Nucleus::computeExcitationEnergy ( ) const

Compute the current excitation energy.

Returns: the excitation energy [MeV]

Definition at line 272 of file G4INCLNucleus.cc.

                                                   {
     const G4double totalEnergy = computeTotalEnergy();
     const G4double separationEnergies = computeSeparationEnergyBalance();
 
     return totalEnergy - initialInternalEnergy - separationEnergies;
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ computeOneNucleonRecoilKinematics()

void G4INCL::Nucleus::computeOneNucleonRecoilKinematics ( )

private

Compute the recoil kinematics for a 1-nucleon remnant.

Puts the remnant nucleon on mass shell and tries to enforce approximate energy conservation by modifying the masses of the outgoing particles.

Definition at line 494 of file G4INCLNucleus.cc.

                                                   {
     // We should be here only if the nucleus contains only one nucleon
 // assert(theStore->getParticles().size()==1);
 
     // No excitation energy!
     theExcitationEnergy = 0.0;
 
     // Move the nucleon to the outgoing list
     Particle *remN = theStore->getParticles().front();
     theA -= remN->getA();
     theZ -= remN->getZ();
     theStore->particleHasBeenEjected(remN);
     theStore->addToOutgoing(remN);
     remN->setEmissionTime(theStore->getBook().getCurrentTime());
 
     // Treat the special case of a remaining delta
     if(remN->isDelta()) {
       IAvatar *decay = new DecayAvatar(remN, 0.0, NULL);
       FinalState *fs = decay->getFinalState();
       // Eject the pion
       ParticleList const &created = fs->getCreatedParticles();
       for(ParticleIter j=created.begin(), e=created.end(); j!=e; ++j)
         theStore->addToOutgoing(*j);
       delete fs;
       delete decay;
     }
 
     // Do different things depending on how many outgoing particles we have
     ParticleList const &outgoing = theStore->getOutgoingParticles();
     if(outgoing.size() == 2) {
 
       INCL_DEBUG("Two particles in the outgoing channel, applying exact two-body kinematics" << '\n');
 
       // Can apply exact 2-body kinematics here. Keep the CM emission angle of
       // the first particle.
       Particle *p1 = outgoing.front(), *p2 = outgoing.back();
       const ThreeVector aBoostVector = incomingMomentum / initialEnergy;
       // Boost to the initial CM
       p1->boost(aBoostVector);
       const G4double sqrts = std::sqrt(initialEnergy*initialEnergy - incomingMomentum.mag2());
       const G4double pcm = KinematicsUtils::momentumInCM(sqrts, p1->getMass(), p2->getMass());
       const G4double scale = pcm/(p1->getMomentum().mag());
       // Reset the momenta
       p1->setMomentum(p1->getMomentum()*scale);
       p2->setMomentum(-p1->getMomentum());
       p1->adjustEnergyFromMomentum();
       p2->adjustEnergyFromMomentum();
       // Unboost
       p1->boost(-aBoostVector);
       p2->boost(-aBoostVector);
 
     } else {
 
       INCL_DEBUG("Trying to adjust final-state momenta to achieve energy and momentum conservation" << '\n');
 
       const G4int maxIterations=8;
       G4double totalEnergy, energyScale;
       G4double val=1.E+100, oldVal=1.E+100, oldOldVal=1.E+100, oldOldOldVal;
       ThreeVector totalMomentum, deltaP;
       std::vector<ThreeVector> minMomenta;  // use it to store the particle momenta that minimize the merit function
 
       // Reserve the vector size
       minMomenta.reserve(outgoing.size());
 
       // Compute the initial total momentum
       totalMomentum.setX(0.0);
       totalMomentum.setY(0.0);
       totalMomentum.setZ(0.0);
       for(ParticleIter i=outgoing.begin(), e=outgoing.end(); i!=e; ++i)
         totalMomentum += (*i)->getMomentum();
 
       // Compute the initial total energy
       totalEnergy = 0.0;
       for(ParticleIter i=outgoing.begin(), e=outgoing.end(); i!=e; ++i)
         totalEnergy += (*i)->getEnergy();
 
       // Iterative algorithm starts here:
       for(G4int iterations=0; iterations < maxIterations; ++iterations) {
 
         // Save the old merit-function values
         oldOldOldVal = oldOldVal;
         oldOldVal = oldVal;
         oldVal = val;
 
         if(iterations%2 == 0) {
           INCL_DEBUG("Momentum step" << '\n');
           // Momentum step: modify all the particle momenta
           deltaP = incomingMomentum - totalMomentum;
           G4double pOldTot = 0.0;
           for(ParticleIter i=outgoing.begin(), e=outgoing.end(); i!=e; ++i)
             pOldTot += (*i)->getMomentum().mag();
           for(ParticleIter i=outgoing.begin(), e=outgoing.end(); i!=e; ++i) {
             const ThreeVector mom = (*i)->getMomentum();
             (*i)->setMomentum(mom + deltaP*mom.mag()/pOldTot);
             (*i)->adjustEnergyFromMomentum();
           }
         } else {
           INCL_DEBUG("Energy step" << '\n');
           // Energy step: modify all the particle momenta
           energyScale = initialEnergy/totalEnergy;
           for(ParticleIter i=outgoing.begin(), e=outgoing.end(); i!=e; ++i) {
             const ThreeVector mom = (*i)->getMomentum();
             G4double pScale = ((*i)->getEnergy()*energyScale - std::pow((*i)->getMass(),2))/mom.mag2();
             if(pScale>0) {
               (*i)->setEnergy((*i)->getEnergy()*energyScale);
               (*i)->adjustMomentumFromEnergy();
             }
           }
         }
 
         // Compute the current total momentum and energy
         totalMomentum.setX(0.0);
         totalMomentum.setY(0.0);
         totalMomentum.setZ(0.0);
         totalEnergy = 0.0;
         for(ParticleIter i=outgoing.begin(), e=outgoing.end(); i!=e; ++i) {
           totalMomentum += (*i)->getMomentum();
           totalEnergy += (*i)->getEnergy();
         }
 
         // Merit factor
         val = std::pow(totalEnergy - initialEnergy,2) +
           0.25*(totalMomentum - incomingMomentum).mag2();
         INCL_DEBUG("Merit function: val=" << val << ", oldVal=" << oldVal << ", oldOldVal=" << oldOldVal << ", oldOldOldVal=" << oldOldOldVal << '\n');
 
         // Store the minimum
         if(val < oldVal) {
           INCL_DEBUG("New minimum found, storing the particle momenta" << '\n');
           minMomenta.clear();
           for(ParticleIter i=outgoing.begin(), e=outgoing.end(); i!=e; ++i)
             minMomenta.push_back((*i)->getMomentum());
         }
 
         // Stop the algorithm if the search diverges
         if(val > oldOldVal && oldVal > oldOldOldVal) {
           INCL_DEBUG("Search is diverging, breaking out of the iteration loop: val=" << val << ", oldVal=" << oldVal << ", oldOldVal=" << oldOldVal << ", oldOldOldVal=" << oldOldOldVal << '\n');
           break;
         }
       }
 
       // We should have made at least one successful iteration here
 // assert(minMomenta.size()==outgoing.size());
 
       // Apply the optimal momenta
       INCL_DEBUG("Applying the solution" << '\n');
       std::vector<ThreeVector>::const_iterator v = minMomenta.begin();
       for(ParticleIter i=outgoing.begin(), e=outgoing.end(); i!=e; ++i, ++v) {
         (*i)->setMomentum(*v);
         (*i)->adjustEnergyFromMomentum();
         INCL_DATABLOCK((*i)->print());
       }
 
     }
 
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ computeRecoilKinematics()

void G4INCL::Nucleus::computeRecoilKinematics ( )

Compute the recoil momentum and spin of the nucleus.

Definition at line 226 of file G4INCLNucleus.cc.

                                         {
     // If the remnant consists of only one nucleon, we need to apply a special
     // procedure to put it on mass shell.
     if(theA==1) {
       emitInsidePions();
       computeOneNucleonRecoilKinematics();
       remnant=false;
       return;
     }
 
     // Compute the recoil momentum and angular momentum
     theMomentum = incomingMomentum;
     theSpin = incomingAngularMomentum;
 
     ParticleList const &outgoing = theStore->getOutgoingParticles();
     for(ParticleIter p=outgoing.begin(), e=outgoing.end(); p!=e; ++p) {
       theMomentum -= (*p)->getMomentum();
       theSpin -= (*p)->getAngularMomentum();
     }
     if(theProjectileRemnant) {
       theMomentum -= theProjectileRemnant->getMomentum();
       theSpin -= theProjectileRemnant->getAngularMomentum();
     }
 
     // Subtract orbital angular momentum
     thePosition = computeCenterOfMass();
     theSpin -= (thePosition-initialCenterOfMass).vector(theMomentum);
 
     setMass(ParticleTable::getTableMass(theA,theZ) + theExcitationEnergy);
     adjustEnergyFromMomentum();
     remnant=true;
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ computeSeparationEnergyBalance()

G4double G4INCL::Nucleus::computeSeparationEnergyBalance ( ) const

inline

Outgoing - incoming separation energies.

Used by CDPP.

Definition at line 116 of file G4INCLNucleus.hh.

                                                     {
       G4double S = 0.0;
       ParticleList const &outgoing = theStore->getOutgoingParticles();
       for(ParticleIter i=outgoing.begin(), e=outgoing.end(); i!=e; ++i) {
         const ParticleType t = (*i)->getType();
         switch(t) {
           case Proton:
           case Neutron:
           case DeltaPlusPlus:
           case DeltaPlus:
           case DeltaZero:
           case DeltaMinus:
             S += thePotential->getSeparationEnergy(*i);
             break;
           case Composite:
             S += (*i)->getZ() * thePotential->getSeparationEnergy(Proton)
               + ((*i)->getA() - (*i)->getZ()) * thePotential->getSeparationEnergy(Neutron);
             break;
           case PiPlus:
             S += thePotential->getSeparationEnergy(Proton) - thePotential->getSeparationEnergy(Neutron);
             break;
           case PiMinus:
             S += thePotential->getSeparationEnergy(Neutron) - thePotential->getSeparationEnergy(Proton);
             break;
           default:
             break;
         }
       }
 
       S -= theNpInitial * thePotential->getSeparationEnergy(Proton);
       S -= theNnInitial * thePotential->getSeparationEnergy(Neutron);
       return S;
     }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ computeTotalEnergy()

G4double G4INCL::Nucleus::computeTotalEnergy ( ) const

Compute the current total energy.

Returns: the total energy [MeV]

Definition at line 212 of file G4INCLNucleus.cc.

                                              {
     G4double totalEnergy = 0.0;
     ParticleList const &inside = theStore->getParticles();
     for(ParticleIter p=inside.begin(), e=inside.end(); p!=e; ++p) {
       if((*p)->isNucleon()) // Ugly: we should calculate everything using total energies!
         totalEnergy += (*p)->getKineticEnergy() - (*p)->getPotentialEnergy();
       else if((*p)->isResonance())
         totalEnergy += (*p)->getEnergy() - (*p)->getPotentialEnergy() - ParticleTable::effectiveNucleonMass;
       else
         totalEnergy += (*p)->getEnergy() - (*p)->getPotentialEnergy();
     }
     return totalEnergy;
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ containsDeltas()

G4bool G4INCL::Nucleus::containsDeltas ( )

inline

Returns true if the nucleus contains any deltas.

Definition at line 231 of file G4INCLNucleus.hh.

                                    {
       ParticleList const &inside = theStore->getParticles();
       for(ParticleIter i=inside.begin(), e=inside.end(); i!=e; ++i)
         if((*i)->isDelta()) return true;
       return false;
     }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ decayInsideDeltas()

G4bool G4INCL::Nucleus::decayInsideDeltas ( )

Force the decay of deltas inside the nucleus.

Returns: true if any delta was forced to decay.

Definition at line 351 of file G4INCLNucleus.cc.

                                     {
     /* If there is a pion potential, do nothing (deltas will be counted as
      * excitation energy).
      * If, however, the remnant is unphysical (Z<0 or Z>A), force the deltas to
      * decay and get rid of all the pions. In case you're wondering, you can
      * end up with Z<0 or Z>A if the remnant contains more pi- than protons or
      * more pi+ than neutrons, respectively.
      */
     const G4bool unphysicalRemnant = (theZ<0 || theZ>theA);
     if(thePotential->hasPionPotential() && !unphysicalRemnant)
       return false;
 
     // Build a list of deltas (avoid modifying the list you are iterating on).
     ParticleList const &inside = theStore->getParticles();
     ParticleList deltas;
     for(ParticleIter i=inside.begin(), e=inside.end(); i!=e; ++i)
       if((*i)->isDelta()) deltas.push_back((*i));
 
     // Loop over the deltas, make them decay
     for(ParticleIter i=deltas.begin(), e=deltas.end(); i!=e; ++i) {
       INCL_DEBUG("Decay inside delta particle:" << '\n'
           << (*i)->print() << '\n');
       // Create a forced-decay avatar. Note the last boolean parameter. Note
       // also that if the remnant is unphysical we more or less explicitly give
       // up energy conservation and CDPP by passing a NULL pointer for the
       // nucleus.
       IAvatar *decay;
       if(unphysicalRemnant) {
         INCL_WARN("Forcing delta decay inside an unphysical remnant (A=" << theA
              << ", Z=" << theZ << "). Might lead to energy-violation warnings."
              << '\n');
         decay = new DecayAvatar((*i), 0.0, NULL, true);
       } else
         decay = new DecayAvatar((*i), 0.0, this, true);
       FinalState *fs = decay->getFinalState();
 
       // The pion can be ejected only if we managed to satisfy energy
       // conservation and if pion emission does not lead to negative excitation
       // energies.
       if(fs->getValidity()==ValidFS) {
         // Apply the final state to the nucleus
         applyFinalState(fs);
       }
       delete fs;
       delete decay;
     }
 
     // If the remnant is unphysical, emit all the pions
     if(unphysicalRemnant) {
       INCL_DEBUG("Remnant is unphysical: Z=" << theZ << ", A=" << theA << ", emitting all the pions" << '\n');
       emitInsidePions();
     }
 
     return true;
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ decayMe()

G4bool G4INCL::Nucleus::decayMe ( )

Force the phase-space decay of the Nucleus.

Only applied if Z==0 or Z==A.

Returns: true if the nucleus was forced to decay.

Definition at line 430 of file G4INCLNucleus.cc.

                           {
     // Do the phase-space decay only if Z=0 or Z=A
     if(theA<=1 || (theZ!=0 && theA!=theZ))
       return false;
 
     ParticleList decayProducts = ClusterDecay::decay(this);
     for(ParticleIter j=decayProducts.begin(), e=decayProducts.end(); j!=e; ++j)
       theStore->addToOutgoing(*j);
 
     return true;
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ decayOutgoingClusters()

G4bool G4INCL::Nucleus::decayOutgoingClusters ( )

Force the decay of unstable outgoing clusters.

Returns: true if any cluster was forced to decay.

Definition at line 407 of file G4INCLNucleus.cc.

                                         {
     ParticleList const &out = theStore->getOutgoingParticles();
     ParticleList clusters;
     for(ParticleIter i=out.begin(), e=out.end(); i!=e; ++i) {
       if((*i)->isCluster()) clusters.push_back((*i));
     }
     if(clusters.empty()) return false;
 
     for(ParticleIter i=clusters.begin(), e=clusters.end(); i!=e; ++i) {
       Cluster *cluster = dynamic_cast<Cluster*>(*i); // Can't avoid using a cast here
 // assert(cluster);
 #ifdef INCLXX_IN_GEANT4_MODE
       if(!cluster)
         continue;
 #endif
       cluster->deleteParticles(); // Don't need them
       ParticleList decayProducts = ClusterDecay::decay(cluster);
       for(ParticleIter j=decayProducts.begin(), end=decayProducts.end(); j!=end; ++j)
         theStore->addToOutgoing(*j);
     }
     return true;
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ decayOutgoingDeltas()

G4bool G4INCL::Nucleus::decayOutgoingDeltas ( )

Force the decay of outgoing deltas.

Returns: true if any delta was forced to decay.

Definition at line 299 of file G4INCLNucleus.cc.

                                       {
     ParticleList const &out = theStore->getOutgoingParticles();
     ParticleList deltas;
     for(ParticleIter i=out.begin(), e=out.end(); i!=e; ++i) {
       if((*i)->isDelta()) deltas.push_back((*i));
     }
     if(deltas.empty()) return false;
 
     for(ParticleIter i=deltas.begin(), e=deltas.end(); i!=e; ++i) {
       INCL_DEBUG("Decay outgoing delta particle:" << '\n'
           << (*i)->print() << '\n');
       const ThreeVector beta = -(*i)->boostVector();
       const G4double deltaMass = (*i)->getMass();
 
       // Set the delta momentum to zero and sample the decay in the CM frame.
       // This makes life simpler if we are using real particle masses.
       (*i)->setMomentum(ThreeVector());
       (*i)->setEnergy((*i)->getMass());
 
       // Use a DecayAvatar
       IAvatar *decay = new DecayAvatar((*i), 0.0, NULL);
       FinalState *fs = decay->getFinalState();
       Particle * const pion = fs->getCreatedParticles().front();
       Particle * const nucleon = fs->getModifiedParticles().front();
 
       // Adjust the decay momentum if we are using the real masses
       const G4double decayMomentum = KinematicsUtils::momentumInCM(deltaMass,
           nucleon->getTableMass(),
           pion->getTableMass());
       ThreeVector newMomentum = pion->getMomentum();
       newMomentum *= decayMomentum / newMomentum.mag();
 
       pion->setTableMass();
       pion->setMomentum(newMomentum);
       pion->adjustEnergyFromMomentum();
       pion->setEmissionTime(nucleon->getEmissionTime());
       pion->boost(beta);
 
       nucleon->setTableMass();
       nucleon->setMomentum(-newMomentum);
       nucleon->adjustEnergyFromMomentum();
       nucleon->boost(beta);
 
       theStore->addToOutgoing(pion);
 
       delete fs;
       delete decay;
     }
 
     return true;
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ deleteProjectileRemnant()

void G4INCL::Nucleus::deleteProjectileRemnant ( )

inline

Delete the projectile remnant.

Definition at line 337 of file G4INCLNucleus.hh.

                                    {
       delete theProjectileRemnant;
       theProjectileRemnant = NULL;
     }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ emitInsidePions()

void G4INCL::Nucleus::emitInsidePions ( )

Force emission of all pions inside the nucleus.

Definition at line 442 of file G4INCLNucleus.cc.

                                 {
     /* Forcing emissions of all pions in the nucleus. This probably violates
      * energy conservation (although the computation of the recoil kinematics
      * might sweep this under the carpet).
      */
     INCL_WARN("Forcing emissions of all pions in the nucleus." << '\n');
 
     // Emit the pions with this kinetic energy
     const G4double tinyPionEnergy = 0.1; // MeV
 
     // Push out the emitted pions
     ParticleList const &inside = theStore->getParticles();
     ParticleList toEject;
     for(ParticleIter i=inside.begin(), e=inside.end(); i!=e; ++i) {
       if((*i)->isPion()) {
         Particle * const thePion = *i;
         INCL_DEBUG("Forcing emission of the following particle: "
                    << thePion->print() << '\n');
         thePion->setEmissionTime(theStore->getBook().getCurrentTime());
         // Correction for real masses
         const G4double theQValueCorrection = thePion->getEmissionQValueCorrection(theA,theZ);
         const G4double kineticEnergyOutside = thePion->getKineticEnergy() - thePion->getPotentialEnergy() + theQValueCorrection;
         thePion->setTableMass();
         if(kineticEnergyOutside > 0.0)
           thePion->setEnergy(thePion->getMass()+kineticEnergyOutside);
         else
           thePion->setEnergy(thePion->getMass()+tinyPionEnergy);
         thePion->adjustMomentumFromEnergy();
         thePion->setPotentialEnergy(0.);
         theZ -= thePion->getZ();
         toEject.push_back(thePion);
       }
     }
     for(ParticleIter i=toEject.begin(), e=toEject.end(); i!=e; ++i) {
       theStore->particleHasBeenEjected(*i);
       theStore->addToOutgoing(*i);
     }
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ fillEventInfo()

void G4INCL::Nucleus::fillEventInfo ( EventInfo * eventInfo )

Fill the event info which contains INCL output data

Definition at line 650 of file G4INCLNucleus.cc.

                                                   {
     eventInfo->nParticles = 0;
     G4bool isNucleonAbsorption = false;
 
     G4bool isPionAbsorption = false;
     // It is possible to have pion absorption event only if the
     // projectile is pion.
     if(eventInfo->projectileType == PiPlus ||
        eventInfo->projectileType == PiMinus ||
        eventInfo->projectileType == PiZero) {
       isPionAbsorption = true;
     }
 
     // Forced CN
     eventInfo->forcedCompoundNucleus = tryCN;
 
     // Outgoing particles
     ParticleList const &outgoingParticles = getStore()->getOutgoingParticles();
 
     // Check if we have a nucleon absorption event: nucleon projectile
     // and no ejected particles.
     if(outgoingParticles.size() == 0 &&
        (eventInfo->projectileType == Proton ||
     eventInfo->projectileType == Neutron)) {
       isNucleonAbsorption = true;
     }
 
     // Reset the remnant counter
     eventInfo->nRemnants = 0;
     eventInfo->history.clear();
 
     for(ParticleIter i=outgoingParticles.begin(), e=outgoingParticles.end(); i!=e; ++i ) {
       // We have a pion absorption event only if the projectile is
       // pion and there are no ejected pions.
       if(isPionAbsorption) {
         if((*i)->isPion()) {
           isPionAbsorption = false;
         }
       }
 
       eventInfo->A[eventInfo->nParticles] = (*i)->getA();
       eventInfo->Z[eventInfo->nParticles] = (*i)->getZ();
       eventInfo->emissionTime[eventInfo->nParticles] = (*i)->getEmissionTime();
       eventInfo->EKin[eventInfo->nParticles] = (*i)->getKineticEnergy();
       ThreeVector mom = (*i)->getMomentum();
       eventInfo->px[eventInfo->nParticles] = mom.getX();
       eventInfo->py[eventInfo->nParticles] = mom.getY();
       eventInfo->pz[eventInfo->nParticles] = mom.getZ();
       eventInfo->theta[eventInfo->nParticles] = Math::toDegrees(mom.theta());
       eventInfo->phi[eventInfo->nParticles] = Math::toDegrees(mom.phi());
       eventInfo->origin[eventInfo->nParticles] = -1;
       eventInfo->history.push_back("");
       eventInfo->nParticles++;
     }
     eventInfo->nucleonAbsorption = isNucleonAbsorption;
     eventInfo->pionAbsorption = isPionAbsorption;
     eventInfo->nCascadeParticles = eventInfo->nParticles;
 
     // Projectile-like remnant characteristics
     if(theProjectileRemnant && theProjectileRemnant->getA()>0) {
       eventInfo->ARem[eventInfo->nRemnants] = theProjectileRemnant->getA();
       eventInfo->ZRem[eventInfo->nRemnants] = theProjectileRemnant->getZ();
       G4double eStar = theProjectileRemnant->getExcitationEnergy();
       if(std::abs(eStar)<1E-10)
         eStar = 0.0; // blame rounding and set the excitation energy to zero
       eventInfo->EStarRem[eventInfo->nRemnants] = eStar;
       if(eventInfo->EStarRem[eventInfo->nRemnants]<0.) {
     INCL_WARN("Negative excitation energy in projectile-like remnant! EStarRem = " << eventInfo->EStarRem[eventInfo->nRemnants] << '\n');
       }
       const ThreeVector &spin = theProjectileRemnant->getSpin();
       if(eventInfo->ARem[eventInfo->nRemnants]%2==0) { // even-A nucleus
     eventInfo->JRem[eventInfo->nRemnants] = (G4int) (spin.mag()/PhysicalConstants::hc + 0.5);
       } else { // odd-A nucleus
     eventInfo->JRem[eventInfo->nRemnants] = ((G4int) (spin.mag()/PhysicalConstants::hc)) + 0.5;
       }
       eventInfo->EKinRem[eventInfo->nRemnants] = theProjectileRemnant->getKineticEnergy();
       const ThreeVector &mom = theProjectileRemnant->getMomentum();
       eventInfo->pxRem[eventInfo->nRemnants] = mom.getX();
       eventInfo->pyRem[eventInfo->nRemnants] = mom.getY();
       eventInfo->pzRem[eventInfo->nRemnants] = mom.getZ();
       eventInfo->jxRem[eventInfo->nRemnants] = spin.getX() / PhysicalConstants::hc;
       eventInfo->jyRem[eventInfo->nRemnants] = spin.getY() / PhysicalConstants::hc;
       eventInfo->jzRem[eventInfo->nRemnants] = spin.getZ() / PhysicalConstants::hc;
       eventInfo->thetaRem[eventInfo->nRemnants] = Math::toDegrees(mom.theta());
       eventInfo->phiRem[eventInfo->nRemnants] = Math::toDegrees(mom.phi());
       eventInfo->nRemnants++;
     }
 
     // Target-like remnant characteristics
     if(hasRemnant()) {
       eventInfo->ARem[eventInfo->nRemnants] = getA();
       eventInfo->ZRem[eventInfo->nRemnants] = getZ();
       eventInfo->EStarRem[eventInfo->nRemnants] = getExcitationEnergy();
       if(eventInfo->EStarRem[eventInfo->nRemnants]<0.) {
     INCL_WARN("Negative excitation energy in target-like remnant! EStarRem = " << eventInfo->EStarRem[eventInfo->nRemnants] << '\n');
       }
       const ThreeVector &spin = getSpin();
       if(eventInfo->ARem[eventInfo->nRemnants]%2==0) { // even-A nucleus
     eventInfo->JRem[eventInfo->nRemnants] = (G4int) (spin.mag()/PhysicalConstants::hc + 0.5);
       } else { // odd-A nucleus
     eventInfo->JRem[eventInfo->nRemnants] = ((G4int) (spin.mag()/PhysicalConstants::hc)) + 0.5;
       }
       eventInfo->EKinRem[eventInfo->nRemnants] = getKineticEnergy();
       const ThreeVector &mom = getMomentum();
       eventInfo->pxRem[eventInfo->nRemnants] = mom.getX();
       eventInfo->pyRem[eventInfo->nRemnants] = mom.getY();
       eventInfo->pzRem[eventInfo->nRemnants] = mom.getZ();
       eventInfo->jxRem[eventInfo->nRemnants] = spin.getX() / PhysicalConstants::hc;
       eventInfo->jyRem[eventInfo->nRemnants] = spin.getY() / PhysicalConstants::hc;
       eventInfo->jzRem[eventInfo->nRemnants] = spin.getZ() / PhysicalConstants::hc;
       eventInfo->thetaRem[eventInfo->nRemnants] = Math::toDegrees(mom.theta());
       eventInfo->phiRem[eventInfo->nRemnants] = Math::toDegrees(mom.phi());
       eventInfo->nRemnants++;
     }
 
     // Global counters, flags, etc.
     Book const &theBook = theStore->getBook();
     eventInfo->nCollisions = theBook.getAcceptedCollisions();
     eventInfo->nBlockedCollisions = theBook.getBlockedCollisions();
     eventInfo->nDecays = theBook.getAcceptedDecays();
     eventInfo->nBlockedDecays = theBook.getBlockedDecays();
     eventInfo->firstCollisionTime = theBook.getFirstCollisionTime();
     eventInfo->firstCollisionXSec = theBook.getFirstCollisionXSec();
     eventInfo->firstCollisionSpectatorPosition = theBook.getFirstCollisionSpectatorPosition();
     eventInfo->firstCollisionSpectatorMomentum = theBook.getFirstCollisionSpectatorMomentum();
     eventInfo->firstCollisionIsElastic = theBook.getFirstCollisionIsElastic();
     eventInfo->nReflectionAvatars = theBook.getAvatars(SurfaceAvatarType);
     eventInfo->nCollisionAvatars = theBook.getAvatars(CollisionAvatarType);
     eventInfo->nDecayAvatars = theBook.getAvatars(DecayAvatarType);
     eventInfo->nEnergyViolationInteraction = theBook.getEnergyViolationInteraction();
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ finalizeProjectileRemnant()

void G4INCL::Nucleus::finalizeProjectileRemnant ( const G4double emissionTime )

Finalise the projectile remnant.

Complete the treatment of the projectile remnant. If it contains nucleons, assign its excitation energy and spin. Move stuff to the outgoing list, if appropriate.

Parameters

emissionTime the emission time of the projectile remnant

Definition at line 834 of file G4INCLNucleus.cc.

                                                                        {
     // Deal with the projectile remnant
     const G4int prA = theProjectileRemnant->getA();
     if(prA>=1) {
       // Set the mass
       const G4double aMass = theProjectileRemnant->getInvariantMass();
       theProjectileRemnant->setMass(aMass);
 
       // Compute the excitation energy from the invariant mass
       const G4double anExcitationEnergy = aMass
         - ParticleTable::getTableMass(prA, theProjectileRemnant->getZ());
 
       // Set the excitation energy
       theProjectileRemnant->setExcitationEnergy(anExcitationEnergy);
 
       // No spin!
       theProjectileRemnant->setSpin(ThreeVector());
 
       // Set the emission time
       theProjectileRemnant->setEmissionTime(anEmissionTime);
     }
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ getConservationBalance()

Nucleus::ConservationBalance G4INCL::Nucleus::getConservationBalance	(	EventInfo const &	theEventInfo,
		const G4bool	afterRecoil
	)		const

Compute charge, mass, energy and momentum balance.

Definition at line 782 of file G4INCLNucleus.cc.

                                                                                                                           {
     ConservationBalance theBalance;
     // Initialise balance variables with the incoming values
     theBalance.Z = theEventInfo.Zp + theEventInfo.Zt;
     theBalance.A = theEventInfo.Ap + theEventInfo.At;
 
     theBalance.energy = getInitialEnergy();
     theBalance.momentum = getIncomingMomentum();
 
     // Process outgoing particles
     ParticleList const &outgoingParticles = theStore->getOutgoingParticles();
     for(ParticleIter i=outgoingParticles.begin(), e=outgoingParticles.end(); i!=e; ++i ) {
       theBalance.Z -= (*i)->getZ();
       theBalance.A -= (*i)->getA();
       // For outgoing clusters, the total energy automatically includes the
       // excitation energy
       theBalance.energy -= (*i)->getEnergy(); // Note that outgoing particles should have the real mass
       theBalance.momentum -= (*i)->getMomentum();
     }
 
     // Projectile-like remnant contribution, if present
     if(theProjectileRemnant && theProjectileRemnant->getA()>0) {
       theBalance.Z -= theProjectileRemnant->getZ();
       theBalance.A -= theProjectileRemnant->getA();
       theBalance.energy -= ParticleTable::getTableMass(theProjectileRemnant->getA(),theProjectileRemnant->getZ()) +
         theProjectileRemnant->getExcitationEnergy();
       theBalance.energy -= theProjectileRemnant->getKineticEnergy();
       theBalance.momentum -= theProjectileRemnant->getMomentum();
     }
 
     // Target-like remnant contribution, if present
     if(hasRemnant()) {
       theBalance.Z -= getZ();
       theBalance.A -= getA();
       theBalance.energy -= ParticleTable::getTableMass(getA(),getZ()) +
         getExcitationEnergy();
       if(afterRecoil)
         theBalance.energy -= getKineticEnergy();
       theBalance.momentum -= getMomentum();
     }
 
     return theBalance;
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ getDensity()

NuclearDensity const* G4INCL::Nucleus::getDensity ( ) const

inline

Getter for theDensity.

Definition at line 365 of file G4INCLNucleus.hh.

365 { return theDensity; };

G4INCL::Nucleus::theDensity

NuclearDensity const * theDensity

Pointer to the NuclearDensity object.

Definition: G4INCLNucleus.hh:414

Here is the caller graph for this function:

◆ getExcitationEnergy()

G4double G4INCL::Nucleus::getExcitationEnergy ( ) const

inline

Get the excitation energy of the nucleus.

Method computeRecoilKinematics() should be called first.

Definition at line 228 of file G4INCLNucleus.hh.

228 { return theExcitationEnergy; }

G4INCL::Cluster::theExcitationEnergy

G4double theExcitationEnergy

Definition: G4INCLCluster.hh:444

Here is the caller graph for this function:

◆ getIncomingAngularMomentum()

const ThreeVector& G4INCL::Nucleus::getIncomingAngularMomentum ( ) const

inline

Get the incoming angular-momentum vector.

Definition at line 206 of file G4INCLNucleus.hh.

206 { return incomingAngularMomentum; }

G4INCL::Nucleus::incomingAngularMomentum

ThreeVector incomingAngularMomentum

Definition: G4INCLNucleus.hh:385

Here is the caller graph for this function:

◆ getIncomingMomentum()

const ThreeVector& G4INCL::Nucleus::getIncomingMomentum ( ) const

inline

Get the incoming momentum vector.

Definition at line 214 of file G4INCLNucleus.hh.

                                                    {
       return incomingMomentum;
     }

Here is the caller graph for this function:

◆ getInitialA()

G4int G4INCL::Nucleus::getInitialA ( ) const

inline

Definition at line 99 of file G4INCLNucleus.hh.

99 { return theInitialA; };

G4INCL::Nucleus::theInitialA

G4int theInitialA

Definition: G4INCLNucleus.hh:379

Here is the caller graph for this function:

◆ getInitialEnergy()

G4double G4INCL::Nucleus::getInitialEnergy ( ) const

inline

Get the initial energy.

Definition at line 222 of file G4INCLNucleus.hh.

222 { return initialEnergy; }

G4INCL::Nucleus::initialEnergy

G4double initialEnergy

Definition: G4INCLNucleus.hh:389

Here is the caller graph for this function:

◆ getInitialInternalEnergy()

G4double G4INCL::Nucleus::getInitialInternalEnergy ( ) const

inline

Definition at line 249 of file G4INCLNucleus.hh.

249 { return initialInternalEnergy; };

G4INCL::Nucleus::initialInternalEnergy

G4double initialInternalEnergy

Definition: G4INCLNucleus.hh:384

Here is the call graph for this function:

Here is the caller graph for this function:

◆ getInitialZ()

G4int G4INCL::Nucleus::getInitialZ ( ) const

inline

Definition at line 100 of file G4INCLNucleus.hh.

100 { return theInitialZ; };

G4INCL::Nucleus::theInitialZ

G4int theInitialZ

Definition: G4INCLNucleus.hh:379

Here is the call graph for this function:

Here is the caller graph for this function:

◆ getNumberOfEnteringNeutrons()

G4int G4INCL::Nucleus::getNumberOfEnteringNeutrons ( ) const

inline

Definition at line 110 of file G4INCLNucleus.hh.

110 { return theNnInitial; };

G4INCL::Nucleus::theNnInitial

G4int theNnInitial

The number of entering neutrons.

Definition: G4INCLNucleus.hh:383

◆ getNumberOfEnteringProtons()

G4int G4INCL::Nucleus::getNumberOfEnteringProtons ( ) const

inline

Definition at line 109 of file G4INCLNucleus.hh.

109 { return theNpInitial; };

G4INCL::Nucleus::theNpInitial

G4int theNpInitial

The number of entering protons.

Definition: G4INCLNucleus.hh:381

◆ getPotential()

NuclearPotential::INuclearPotential const* G4INCL::Nucleus::getPotential ( ) const

inline

Getter for thePotential.

Definition at line 368 of file G4INCLNucleus.hh.

368 { return thePotential; };

G4INCL::Nucleus::thePotential

NuclearPotential::INuclearPotential const * thePotential

Pointer to the NuclearPotential object.

Definition: G4INCLNucleus.hh:417

Here is the call graph for this function:

Here is the caller graph for this function:

◆ getProjectileRemnant()

ProjectileRemnant* G4INCL::Nucleus::getProjectileRemnant ( ) const

inline

Get the projectile remnant.

Definition at line 334 of file G4INCLNucleus.hh.

334 { return theProjectileRemnant; }

G4INCL::Nucleus::theProjectileRemnant

ProjectileRemnant * theProjectileRemnant

Pointer to the quasi-projectile.

Definition: G4INCLNucleus.hh:411

Here is the caller graph for this function:

◆ getStore()

Store* G4INCL::Nucleus::getStore ( ) const

inline

Definition at line 243 of file G4INCLNucleus.hh.

243 {return theStore; };

G4INCL::Nucleus::theStore

Store * theStore

Definition: G4INCLNucleus.hh:390

Here is the caller graph for this function:

◆ getSurfaceRadius()

G4double G4INCL::Nucleus::getSurfaceRadius ( Particle const *const particle ) const

inline

Get the maximum allowed radius for a given particle.

Calls the NuclearDensity::getMaxRFromP() method for nucleons and deltas, and the NuclearDensity::getTrasmissionRadius() method for pions.

Parameters

particle pointer to a particle

Returns: surface radius

Definition at line 298 of file G4INCLNucleus.hh.

                                                                      {
       if(particle->isPion())
         // Temporarily set RPION = RMAX
         return getUniverseRadius();
         //return 0.5*(theDensity->getTransmissionRadius(particle)+getUniverseRadius());
       else {
         const G4double pr = particle->getReflectionMomentum()/thePotential->getFermiMomentum(particle);
         if(pr>=1.)
           return getUniverseRadius();
         else
           return theDensity->getMaxRFromP(particle->getType(), pr);
       }
     }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ getTransmissionBarrier()

G4double G4INCL::Nucleus::getTransmissionBarrier ( Particle const *const p )

inline

Get the transmission barrier.

Definition at line 271 of file G4INCLNucleus.hh.

                                                               {
       const G4double theTransmissionRadius = theDensity->getTransmissionRadius(p);
       const G4double theParticleZ = p->getZ();
       return PhysicalConstants::eSquared*(theZ-theParticleZ)*theParticleZ/theTransmissionRadius;
     }

Here is the call graph for this function:

◆ getTryCompoundNucleus()

G4bool G4INCL::Nucleus::getTryCompoundNucleus ( )

inline

Definition at line 268 of file G4INCLNucleus.hh.

268 { return tryCN; }

G4INCL::Nucleus::tryCN

G4bool tryCN

Definition: G4INCLNucleus.hh:391

Here is the caller graph for this function:

◆ getUniverseRadius()

G4double G4INCL::Nucleus::getUniverseRadius ( ) const

inline

Getter for theUniverseRadius.

Definition at line 313 of file G4INCLNucleus.hh.

313 { return theUniverseRadius; }

G4INCL::Nucleus::theUniverseRadius

G4double theUniverseRadius

The radius of the universe.

Definition: G4INCLNucleus.hh:399

Here is the caller graph for this function:

◆ hasRemnant()

G4bool G4INCL::Nucleus::hasRemnant ( ) const

inline

Does the nucleus give a cascade remnant?

To be called after computeRecoilKinematics().

Definition at line 261 of file G4INCLNucleus.hh.

261 { return remnant; }

G4INCL::Nucleus::remnant

G4bool remnant

Definition: G4INCLNucleus.hh:387

Here is the call graph for this function:

Here is the caller graph for this function:

◆ INCL_DECLARE_ALLOCATION_POOL()

G4INCL::Nucleus::INCL_DECLARE_ALLOCATION_POOL ( Nucleus )

private

◆ initializeParticles()

void G4INCL::Nucleus::initializeParticles ( )

virtual

Call the Cluster method to generate the initial distribution of particles. At the beginning all particles are assigned as spectators.

Reimplemented from G4INCL::Cluster.

Definition at line 116 of file G4INCLNucleus.cc.

                                     {
     // Reset the variables connected with the projectile remnant
     delete theProjectileRemnant;
     theProjectileRemnant = NULL;
 
     Cluster::initializeParticles();
     for(ParticleIter i=particles.begin(), e=particles.end(); i!=e; ++i) {
       updatePotentialEnergy(*i);
     }
     theStore->add(particles);
     particles.clear();
     initialInternalEnergy = computeTotalEnergy();
     initialCenterOfMass = thePosition;
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ insertParticle()

void G4INCL::Nucleus::insertParticle ( Particle * p )

inline

Insert a new particle (e.g. a projectile) in the nucleus.

Definition at line 83 of file G4INCLNucleus.hh.

                                      {
       theZ += p->getZ();
       theA += p->getA();
       theStore->particleHasEntered(p);
       if(p->isNucleon()) {
         theNpInitial += Math::heaviside(ParticleTable::getIsospin(p->getType()));
         theNnInitial += Math::heaviside(-ParticleTable::getIsospin(p->getType()));
       }
       if(!p->isTargetSpectator()) theStore->getBook().incrementCascading();
     };

Here is the call graph for this function:

Here is the caller graph for this function:

◆ isEventTransparent()

G4bool G4INCL::Nucleus::isEventTransparent ( ) const

Is the event transparent?

To be called at the end of the cascade.

Definition at line 481 of file G4INCLNucleus.cc.

                                            {
 
     Book const &theBook = theStore->getBook();
     const G4int nEventCollisions = theBook.getAcceptedCollisions();
     const G4int nEventDecays = theBook.getAcceptedDecays();
     const G4int nEventClusters = theBook.getEmittedClusters();
     if(nEventCollisions==0 && nEventDecays==0 && nEventClusters==0)
       return true;
 
     return false;
 
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ isNucleusNucleusCollision()

G4bool G4INCL::Nucleus::isNucleusNucleusCollision ( ) const

inline

Is it a nucleus-nucleus collision?

Definition at line 319 of file G4INCLNucleus.hh.

319 { return isNucleusNucleus; }

G4INCL::Nucleus::isNucleusNucleus

G4bool isNucleusNucleus

true if running a nucleus-nucleus collision

Definition: G4INCLNucleus.hh:405

Here is the caller graph for this function:

◆ operator=()

Nucleus& G4INCL::Nucleus::operator= ( const Nucleus & rhs )

Dummy assignment operator to silence Coverity warning.

◆ print()

std::string G4INCL::Nucleus::print ( )

Print the nucleus info

Definition at line 279 of file G4INCLNucleus.cc.

   {
     std::stringstream ss;
     ss << "Particles in the nucleus:" << '\n'
       << "Inside:" << '\n';
     G4int counter = 1;
     ParticleList const &inside = theStore->getParticles();
     for(ParticleIter p=inside.begin(), e=inside.end(); p!=e; ++p) {
       ss << "index = " << counter << '\n'
         << (*p)->print();
       counter++;
     }
     ss <<"Outgoing:" << '\n';
     ParticleList const &outgoing = theStore->getOutgoingParticles();
     for(ParticleIter p=outgoing.begin(), e=outgoing.end(); p!=e; ++p)
       ss << (*p)->print();
 
     return ss.str();
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ propagateParticles()

void G4INCL::Nucleus::propagateParticles ( G4double step )

Propagate the particles one time step.

Parameters

step	length of the time step

Definition at line 208 of file G4INCLNucleus.cc.

                                              {
     INCL_WARN("Useless Nucleus::propagateParticles -method called." << '\n');
   }

Here is the caller graph for this function:

◆ setDensity()

void G4INCL::Nucleus::setDensity ( NuclearDensity const *const d )

inline

Setter for theDensity.

Definition at line 358 of file G4INCLNucleus.hh.

                                                     {
       theDensity=d;
       if(theParticleSampler)
         theParticleSampler->setDensity(theDensity);
     };

Here is the call graph for this function:

◆ setIncomingAngularMomentum()

void G4INCL::Nucleus::setIncomingAngularMomentum ( const ThreeVector & j )

inline

Set the incoming angular-momentum vector.

Definition at line 201 of file G4INCLNucleus.hh.

                                                           {
       incomingAngularMomentum = j;
     }

Here is the caller graph for this function:

◆ setIncomingMomentum()

void G4INCL::Nucleus::setIncomingMomentum ( const ThreeVector & p )

inline

Set the incoming momentum vector.

Definition at line 209 of file G4INCLNucleus.hh.

                                                    {
       incomingMomentum = p;
     }

Here is the caller graph for this function:

◆ setInitialEnergy()

void G4INCL::Nucleus::setInitialEnergy ( const G4double e )

inline

Set the initial energy.

Definition at line 219 of file G4INCLNucleus.hh.

219 { initialEnergy = e; }

e

Float_t e

Definition: extended/medical/dna/range/plot.C:37

G4INCL::Nucleus::initialEnergy

G4double initialEnergy

Definition: G4INCLNucleus.hh:389

Here is the caller graph for this function:

◆ setNucleusNucleusCollision()

void G4INCL::Nucleus::setNucleusNucleusCollision ( )

inline

Set a nucleus-nucleus collision.

Definition at line 322 of file G4INCLNucleus.hh.

322 { isNucleusNucleus=true; }

G4INCL::Nucleus::isNucleusNucleus

G4bool isNucleusNucleus

true if running a nucleus-nucleus collision

Definition: G4INCLNucleus.hh:405

Here is the caller graph for this function:

◆ setParticleNucleusCollision()

void G4INCL::Nucleus::setParticleNucleusCollision ( )

inline

Set a particle-nucleus collision.

Definition at line 325 of file G4INCLNucleus.hh.

325 { isNucleusNucleus=false; }

G4INCL::Nucleus::isNucleusNucleus

G4bool isNucleusNucleus

true if running a nucleus-nucleus collision

Definition: G4INCLNucleus.hh:405

Here is the caller graph for this function:

◆ setProjectileRemnant()

void G4INCL::Nucleus::setProjectileRemnant ( ProjectileRemnant *const c )

inline

Set the projectile remnant.

Definition at line 328 of file G4INCLNucleus.hh.

                                                            {
       delete theProjectileRemnant;
       theProjectileRemnant = c;
     }

Here is the caller graph for this function:

◆ setStore()

void G4INCL::Nucleus::setStore ( Store * s )

inline

Definition at line 244 of file G4INCLNucleus.hh.

                             {
       delete theStore;
       theStore = s;
     };

◆ setUniverseRadius()

void G4INCL::Nucleus::setUniverseRadius ( const G4double universeRadius )

inline

Setter for theUniverseRadius.

Definition at line 316 of file G4INCLNucleus.hh.

316 { theUniverseRadius=universeRadius; }

G4INCL::Nucleus::theUniverseRadius

G4double theUniverseRadius

The radius of the universe.

Definition: G4INCLNucleus.hh:399

◆ updatePotentialEnergy()

void G4INCL::Nucleus::updatePotentialEnergy ( Particle * p ) const

inline

Update the particle potential energy.

Definition at line 353 of file G4INCLNucleus.hh.

                                                          {
       p->setPotentialEnergy(thePotential->computePotentialEnergy(p));
     }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ useFusionKinematics()

void G4INCL::Nucleus::useFusionKinematics ( )

Adjust the kinematics for complete-fusion events.

Definition at line 826 of file G4INCLNucleus.cc.

                                     {
     setEnergy(initialEnergy);
     setMomentum(incomingMomentum);
     setSpin(incomingAngularMomentum);
     theExcitationEnergy = std::sqrt(theEnergy*theEnergy-theMomentum.mag2()) - getTableMass();
     setMass(getTableMass() + theExcitationEnergy);
   }

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

◆ incomingAngularMomentum

ThreeVector G4INCL::Nucleus::incomingAngularMomentum

private

Definition at line 385 of file G4INCLNucleus.hh.

◆ incomingMomentum

ThreeVector G4INCL::Nucleus::incomingMomentum

private

Definition at line 385 of file G4INCLNucleus.hh.

◆ initialCenterOfMass

ThreeVector G4INCL::Nucleus::initialCenterOfMass

private

Definition at line 386 of file G4INCLNucleus.hh.

◆ initialEnergy

G4double G4INCL::Nucleus::initialEnergy

private

Definition at line 389 of file G4INCLNucleus.hh.

◆ initialInternalEnergy

G4double G4INCL::Nucleus::initialInternalEnergy

private

Definition at line 384 of file G4INCLNucleus.hh.

◆ isNucleusNucleus

G4bool G4INCL::Nucleus::isNucleusNucleus

private

true if running a nucleus-nucleus collision

Tells INCL whether to make a projectile-like pre-fragment or not.

Definition at line 405 of file G4INCLNucleus.hh.

◆ projectileA

G4int G4INCL::Nucleus::projectileA

private

The mass number of the projectile.

Definition at line 396 of file G4INCLNucleus.hh.

◆ projectileZ

G4int G4INCL::Nucleus::projectileZ

private

The charge number of the projectile.

Definition at line 394 of file G4INCLNucleus.hh.

◆ remnant

G4bool G4INCL::Nucleus::remnant

private

Definition at line 387 of file G4INCLNucleus.hh.

◆ theDensity

NuclearDensity const* G4INCL::Nucleus::theDensity

private

Pointer to the NuclearDensity object.

Definition at line 414 of file G4INCLNucleus.hh.

◆ theInitialA

G4int G4INCL::Nucleus::theInitialA

private

Definition at line 379 of file G4INCLNucleus.hh.

◆ theInitialZ

G4int G4INCL::Nucleus::theInitialZ

private

Definition at line 379 of file G4INCLNucleus.hh.

◆ theNnInitial

G4int G4INCL::Nucleus::theNnInitial

private

The number of entering neutrons.

Definition at line 383 of file G4INCLNucleus.hh.

◆ theNpInitial

G4int G4INCL::Nucleus::theNpInitial

private

The number of entering protons.

Definition at line 381 of file G4INCLNucleus.hh.

◆ thePotential

NuclearPotential::INuclearPotential const* G4INCL::Nucleus::thePotential

private

Pointer to the NuclearPotential object.

Definition at line 417 of file G4INCLNucleus.hh.

◆ theProjectileRemnant

ProjectileRemnant* G4INCL::Nucleus::theProjectileRemnant

private

Pointer to the quasi-projectile.

Owned by the Nucleus object.

Definition at line 411 of file G4INCLNucleus.hh.

◆ theStore

Store* G4INCL::Nucleus::theStore

private

Definition at line 390 of file G4INCLNucleus.hh.

◆ theUniverseRadius

G4double G4INCL::Nucleus::theUniverseRadius

private

The radius of the universe.

Definition at line 399 of file G4INCLNucleus.hh.

◆ tryCN

G4bool G4INCL::Nucleus::tryCN

private

Definition at line 391 of file G4INCLNucleus.hh.

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

Classes

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ Nucleus() [1/2]

◆ ~Nucleus()

◆ Nucleus() [2/2]

Member Function Documentation

◆ applyFinalState()

◆ computeCenterOfMass()

◆ computeExcitationEnergy()

◆ computeOneNucleonRecoilKinematics()

◆ computeRecoilKinematics()

◆ computeSeparationEnergyBalance()

◆ computeTotalEnergy()

◆ containsDeltas()

◆ decayInsideDeltas()

◆ decayMe()

◆ decayOutgoingClusters()

◆ decayOutgoingDeltas()

◆ deleteProjectileRemnant()

◆ emitInsidePions()

◆ fillEventInfo()

◆ finalizeProjectileRemnant()

◆ getConservationBalance()

◆ getDensity()

◆ getExcitationEnergy()

◆ getIncomingAngularMomentum()

◆ getIncomingMomentum()

◆ getInitialA()

◆ getInitialEnergy()

◆ getInitialInternalEnergy()

◆ getInitialZ()

◆ getNumberOfEnteringNeutrons()

◆ getNumberOfEnteringProtons()

◆ getPotential()

◆ getProjectileRemnant()

◆ getStore()

◆ getSurfaceRadius()

◆ getTransmissionBarrier()

◆ getTryCompoundNucleus()

◆ getUniverseRadius()

◆ hasRemnant()

◆ INCL_DECLARE_ALLOCATION_POOL()

◆ initializeParticles()

◆ insertParticle()

◆ isEventTransparent()

◆ isNucleusNucleusCollision()

◆ operator=()

◆ print()

◆ propagateParticles()

◆ setDensity()

◆ setIncomingAngularMomentum()

◆ setIncomingMomentum()

◆ setInitialEnergy()

◆ setNucleusNucleusCollision()

◆ setParticleNucleusCollision()

◆ setProjectileRemnant()

◆ setStore()

◆ setUniverseRadius()

◆ updatePotentialEnergy()

◆ useFusionKinematics()

Member Data Documentation

◆ incomingAngularMomentum

◆ incomingMomentum

◆ initialCenterOfMass

◆ initialEnergy

◆ initialInternalEnergy

◆ isNucleusNucleus

◆ projectileA

◆ projectileZ

◆ remnant

◆ theDensity

◆ theInitialA

◆ theInitialZ

◆ theNnInitial

◆ theNpInitial