#include <G4INCLNucleus.hh>

Inheritance diagram for G4INCL::Nucleus:

Data Structures
struct	ConservationBalance
	Struct for conservation laws. More...

Public Member Functions
void	addParticle (Particle *const p)

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

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

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

void	applyFinalState (FinalState *)

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

ThreeVector	boostVector () const

void	clearParticles ()

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

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

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

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

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

G4bool	containsAntiKaon ()
	Returns true if the nucleus contains any anti Kaons. More...

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

G4bool	containsEtas ()
	Returns true if the nucleus contains any etas. More...

G4bool	containsKaon ()
	Returns true if the nucleus contains any Kaons. More...

G4bool	containsLambda ()
	Returns true if the nucleus contains any Lambda. More...

G4bool	containsOmegas ()
	Returns true if the nucleus contains any omegas. More...

G4bool	containsSigma ()
	Returns true if the nucleus contains any Sigma. More...

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

G4bool	decayInsideStrangeParticles ()
	Force the transformation of strange particles into a Lambda;. More...

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

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

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

G4bool	decayOutgoingNeutralKaon ()
	Force the transformation of outgoing Neutral Kaon into propation eigenstate. More...

G4bool	decayOutgoingPionResonances (G4double timeThreshold)
	Force the decay of outgoing PionResonances (eta/omega). More...

G4bool	decayOutgoingSigmaZero (G4double timeThreshold)
	Force the decay of outgoing Neutral Sigma. More...

void	deleteParticles ()

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

std::string	dump () const

G4bool	emitInsideKaon ()
	Force emission of all Kaon inside the nucleus. More...

G4int	emitInsideLambda ()
	Force emission of all Lambda (desexitation code with strangeness not implanted yet) More...

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

void	emitInsideStrangeParticles ()
	Force emission of all strange particles inside the nucleus. More...

void	fillEventInfo (EventInfo *eventInfo)

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

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

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

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

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

G4double	getBeta () const

std::vector< G4int >	getBiasCollisionVector () const
	Get the vector list of biased vertices on the particle path. More...

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

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

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

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

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

G4double	getEmissionTime ()

G4double	getEnergy () const

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

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

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

G4double	getHelicity ()

long	getID () const

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

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

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

G4int	getInitialA () const

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

G4double	getInitialInternalEnergy () const

G4int	getInitialS () const

G4int	getInitialZ () const

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

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

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

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

const G4INCL::ThreeVector &	getMomentum () const

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

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

G4int	getNumberOfEnteringKaons () const

G4int	getNumberOfEnteringNeutrons () const

G4int	getNumberOfEnteringPions () const

G4int	getNumberOfEnteringProtons () const

G4int	getNumberOfKaon () const
	Number of Kaon inside de nucleus. More...

ParticipantType	getParticipantType () const

G4double	getParticleBias () const
	Get the particle bias. More...

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

ParticleList const &	getParticles () const

const G4INCL::ThreeVector &	getPosition () const

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

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

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

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

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

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

G4int	getS () const
	Returns the strangeness number. More...

ParticleSpecies	getSpecies () const
	Get the particle species. More...

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

Store *	getStore () const

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

virtual G4double	getTableMass () const
	Get the real particle mass. 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	getTransferQValueCorrection (const G4int AFrom, const G4int ZFrom, const G4int SFrom, const G4int ATo, const G4int ZTo, const G4int STo) const
	Computes correction on the transfer Q-value for hypernuclei. More...

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

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

G4bool	getTryCompoundNucleus ()

G4INCL::ParticleType	getType () const

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

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

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

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

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

void	initializeParticles ()

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

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

G4bool	isAntiKaon () const
	Is this an antiKaon? More...

G4bool	isBaryon () const
	Is this a Baryon? More...

G4bool	isCluster () const

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

G4bool	isEta () const
	Is this an eta? More...

G4bool	isEtaPrime () const
	Is this an etaprime? More...

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

G4bool	isHyperon () const
	Is this an Hyperon? More...

G4bool	isKaon () const
	Is this a Kaon? More...

G4bool	isLambda () const
	Is this a Lambda? More...

G4bool	isMeson () const
	Is this a Meson? More...

G4bool	isNucleon () const

G4bool	isNucleonorLambda () const
	Is this a Nucleon or a Lambda? More...

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

G4bool	isOmega () const
	Is this an omega? More...

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

G4bool	isParticipant () const

G4bool	isPhoton () const
	Is this a photon? More...

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

G4bool	isProjectileSpectator () const

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

G4bool	isSigma () const
	Is this a Sigma? More...

G4bool	isStrange () const
	Is this an Strange? More...

G4bool	isTargetSpectator () const

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

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

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

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

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

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

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

std::string	print ()

std::string	print () const

void	propagate (G4double step)

void	propagateParticles (G4double step)

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

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

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

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

virtual void	rotatePositionAndMomentum (const G4double angle, const ThreeVector &axis)
	Rotate the particle position and 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...

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

void	setBiasCollisionVector (std::vector< G4int > BiasCollisionVector)
	Set the vector list of biased vertices on the particle path. More...

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

void	setEmissionTime (G4double t)

void	setEnergy (G4double energy)

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

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

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

void	setHelicity (G4double h)

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

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

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

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

void	setMass (G4double mass)

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

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

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

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

void	setNumberOfKaon (const G4int NK)

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

void	setParticipantType (ParticipantType const p)

void	setParticleBias (G4double ParticleBias)
	Set the particle bias. More...

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

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

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

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

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

void	setS (const G4int S)
	Set the strangess number of the cluster. More...

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

void	setStore (Store *str)

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

void	setType (ParticleType t)

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

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

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

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

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

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

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

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

virtual	~Nucleus ()

Static Public Member Functions
static void	FillINCLBiasVector (G4double newBias)

static G4double	getBiasFromVector (std::vector< G4int > VectorBias)

static G4double	getTotalBias ()
	General bias vector function. More...

static std::vector< G4int >	MergeVectorBias (Particle const const p1, Particle const const p2)

static std::vector< G4int >	MergeVectorBias (std::vector< G4int > p1, Particle const *const p2)

static void	setINCLBiasVector (std::vector< G4double > NewVector)

Static Public Attributes
static std::vector< G4double >	INCLBiasVector
	Time ordered vector of all bias applied. More...

static G4ThreadLocal G4int	nextBiasedCollisionID = 0

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

Protected Attributes
long	ID

G4int	nCollisions

G4int	nDecays

ParticleList	particles

G4bool	rpCorrelated

G4int	theA

G4double	theEnergy

G4double	theExcitationEnergy

G4double	theFrozenEnergy

G4INCL::ThreeVector	theFrozenMomentum

G4INCL::ThreeVector	theMomentum

G4int	theNKaon
	The number of Kaons inside the nucleus (update during the cascade) More...

ParticipantType	theParticipantType

G4double	theParticleBias

ParticleSampler *	theParticleSampler

G4INCL::ThreeVector	thePosition

G4double	thePotentialEnergy

G4double *	thePropagationEnergy

G4INCL::ThreeVector *	thePropagationMomentum

G4int	theS

ThreeVector	theSpin

G4INCL::ParticleType	theType

G4int	theZ

G4double	uncorrelatedMomentum

Private Member Functions
G4double	computeDynamicalPotential ()
	Compute the dynamical cluster potential. More...

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

Private Attributes
G4double	emissionTime

ThreeVector	incomingAngularMomentum

ThreeVector	incomingMomentum

ThreeVector	initialCenterOfMass

G4double	initialEnergy

G4double	initialInternalEnergy

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

G4bool	outOfWell

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

G4int	projectileS
	The strangeness number of the projectile. More...

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

G4bool	remnant

std::vector< G4int >	theBiasCollisionVector
	Time ordered vector of all biased vertices on the particle path. More...

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

G4double	theHelicity

G4int	theInitialA

G4int	theInitialS

G4int	theInitialZ

G4double	theMass

G4int	theNkaonminusInitial

G4int	theNkaonplusInitial
	The number of entering kaons. More...

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

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

G4int	theNpionminusInitial

G4int	theNpionplusInitial
	The number of entering pions. More...

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

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

Store *	theStore

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

G4bool	tryCN

Static Private Attributes
static G4ThreadLocal long	nextID = 1

Detailed Description

Definition at line 65 of file G4INCLNucleus.hh.

Constructor & Destructor Documentation

◆ Nucleus() [1/2]

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

Definition at line 68 of file G4INCLNucleus.cc.

    : Cluster(charge,mass,strangess,true),
     theInitialZ(charge), theInitialA(mass), theInitialS(strangess),
     theNpInitial(0), theNnInitial(0), 
     theNpionplusInitial(0), theNpionminusInitial(0), 
     theNkaonplusInitial(0), theNkaonminusInitial(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, theS);
 
    theParticleSampler->setPotential(thePotential);
    theParticleSampler->setDensity(theDensity);
 
    if(theUniverseRadius<0)
      theUniverseRadius = theDensity->getMaximumRadius();
    theStore = new Store(conf);
  }

References G4INCL::NuclearDensityFactory::createDensity(), G4INCL::NuclearPotential::createPotential(), G4INCL::NuclearDensity::getMaximumRadius(), G4INCL::Config::getPionPotential(), G4INCL::Config::getPotentialType(), G4INCL::NuclearPotential::INuclearPotential::getSeparationEnergy(), G4INCL::IsospinPotential, G4INCL::Neutron, G4INCL::Proton, G4INCL::ParticleSampler::setDensity(), G4INCL::ParticleTable::setNeutronSeparationEnergy(), G4INCL::ParticleSampler::setPotential(), G4INCL::ParticleTable::setProtonSeparationEnergy(), G4INCL::Particle::theA, theDensity, G4INCL::Cluster::theParticleSampler, thePotential, G4INCL::Particle::theS, theStore, theUniverseRadius, and G4INCL::Particle::theZ.

◆ ~Nucleus()

G4INCL::Nucleus::~Nucleus ( )

virtual

Definition at line 112 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;*/
  }

References deleteProjectileRemnant(), and theStore.

◆ Nucleus() [2/2]

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

Dummy copy constructor to silence Coverity warning.

Member Function Documentation

◆ addParticle()

void G4INCL::Cluster::addParticle ( Particle *const p )

inlineinherited

Add one particle to the cluster. This updates the cluster mass, energy, size, etc.

Definition at line 178 of file G4INCLCluster.hh.

                                         {
      particles.push_back(p);
      theEnergy += p->getEnergy();
      thePotentialEnergy += p->getPotentialEnergy();
      theMomentum += p->getMomentum();
      thePosition += p->getPosition();
      theA += p->getA();
      theZ += p->getZ();
      theS += p->getS();
      nCollisions += p->getNumberOfCollisions();
    }

Referenced by G4INCL::Cluster::Cluster(), and G4INCL::ProjectileRemnant::reset().

◆ addParticles()

void G4INCL::Cluster::addParticles ( ParticleList const & pL )

inlineinherited

Add a list of particles to the cluster.

Definition at line 213 of file G4INCLCluster.hh.

                                              {
      particles = pL;
      updateClusterParameters();
    }

References G4INCL::Cluster::particles, and G4INCL::Cluster::updateClusterParameters().

◆ adjustEnergyFromMomentum()

G4double G4INCL::Particle::adjustEnergyFromMomentum ( )

inherited

Recompute the energy to match the momentum.

Definition at line 155 of file G4INCLParticle.cc.

                                              {
    theEnergy = std::sqrt(theMomentum.mag2() + theMass*theMass);
    return theEnergy;
  }

References G4INCL::ThreeVector::mag2(), G4INCL::Particle::theEnergy, G4INCL::Particle::theMass, and G4INCL::Particle::theMomentum.

◆ adjustMomentumFromEnergy()

const ThreeVector & G4INCL::Particle::adjustMomentumFromEnergy ( )

inherited

Rescale the momentum to match the total energy.

Definition at line 142 of file G4INCLParticle.cc.

                                                        {
    const G4double p2 = theMomentum.mag2();
    G4double newp2 = theEnergy*theEnergy - theMass*theMass;
    if( newp2<0.0 ) {
      INCL_ERROR("Particle has E^2 < m^2." << '\n' << print());
      newp2 = 0.0;
      theEnergy = theMass;
    }
 
    theMomentum *= std::sqrt(newp2/p2);
    return theMomentum;
  }

References INCL_ERROR, G4INCL::ThreeVector::mag2(), G4INCL::Particle::print(), G4INCL::Particle::theEnergy, G4INCL::Particle::theMass, and G4INCL::Particle::theMomentum.

Referenced by G4INCL::Cluster::Cluster(), emitInsideKaon(), emitInsideLambda(), emitInsidePions(), emitInsideStrangeParticles(), G4INCL::CrossSections::interactionDistanceKbarN(), G4INCL::CrossSections::interactionDistanceKN(), G4INCL::CrossSections::interactionDistanceNN(), G4INCL::CrossSections::interactionDistancePiN(), G4INCL::CrossSections::interactionDistanceYN(), G4INCL::ParticleEntryChannel::particleEnters(), G4INCL::TransmissionChannel::particleLeaves(), G4INCL::StandardPropagationModel::shootParticle(), and G4INCL::KinematicsUtils::transformToLocalEnergyFrame().

◆ applyFinalState()

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

Apply reaction final state information to the nucleus.

Definition at line 136 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();
        theS -= (*iter)->getS();
        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());
    }
  }

Referenced by G4INCL::INCL::cascade(), decayInsideDeltas(), decayInsideStrangeParticles(), and G4INCL::INCL::makeCompoundNucleus().

◆ boost()

void G4INCL::Cluster::boost ( const ThreeVector & aBoostVector )

inlineinherited

Boost the cluster with the indicated velocity.

The Cluster is boosted as a whole, just like any Particle object; moreover, the internal components (particles list) are also boosted, according to Alain Boudard's off-shell recipe.

Parameters

aBoostVector the velocity to boost to [c]

Definition at line 344 of file G4INCLCluster.hh.

                                                {
      Particle::boost(aBoostVector);
      for(ParticleIter p=particles.begin(), e=particles.end(); p!=e; ++p) {
        (*p)->boost(aBoostVector);
        // Apply Lorentz contraction to the particle position
        (*p)->lorentzContract(aBoostVector,thePosition);
        (*p)->rpCorrelate();
      }
 
      INCL_DEBUG("Cluster was boosted with (bx,by,bz)=("
          << aBoostVector.getX() << ", " << aBoostVector.getY() << ", " << aBoostVector.getZ() << "):"
          << '\n' << print());
 
    }

References G4INCL::Particle::boost(), G4INCL::ThreeVector::getX(), G4INCL::ThreeVector::getY(), G4INCL::ThreeVector::getZ(), INCL_DEBUG, G4INCL::Cluster::particles, G4INCL::Cluster::print(), and G4INCL::Particle::thePosition.

Referenced by G4INCL::ProjectileRemnant::ProjectileRemnant(), G4INCL::INCL::RecoilCMFunctor::RecoilCMFunctor(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::threeBodyDecay(), and G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::twoBodyDecay().

◆ boostVector()

ThreeVector G4INCL::Particle::boostVector ( ) const

inlineinherited

Returns a three vector we can give to the boost() -method.

In order to go to the particle rest frame you need to multiply the boost vector by -1.0.

Definition at line 412 of file G4INCLParticle.hh.

                                    {
      return theMomentum / theEnergy;
    }

References G4INCL::Particle::theEnergy, and G4INCL::Particle::theMomentum.

Referenced by G4INCL::PhaseSpaceKopylov::generate(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::phaseSpaceDecay(), G4INCL::StandardPropagationModel::shootComposite(), G4INCL::StandardPropagationModel::shootParticle(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::threeBodyDecay(), and G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::twoBodyDecay().

◆ clearParticles()

void G4INCL::Cluster::clearParticles ( )

inlineinherited

Definition at line 143 of file G4INCLCluster.hh.

143{ particles.clear(); }

References G4INCL::Cluster::particles.

Referenced by G4INCL::Cluster::deleteParticles().

◆ computeCenterOfMass()

ThreeVector G4INCL::Nucleus::computeCenterOfMass ( ) const

Compute the current center-of-mass position.

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

Definition at line 268 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;
  }

References cm, G4INCL::Store::getParticles(), and theStore.

Referenced by computeRecoilKinematics().

◆ computeDynamicalPotential()

G4double G4INCL::Cluster::computeDynamicalPotential ( )

inlineprivateinherited

Compute the dynamical cluster potential.

Alain Boudard's boost prescription for low-energy beams requires to define a "dynamical potential" that allows us to conserve momentum and energy when boosting the projectile cluster.

Definition at line 439 of file G4INCLCluster.hh.

                                         {
      G4double theDynamicalPotential = 0.0;
      for(ParticleIter p=particles.begin(), e=particles.end(); p!=e; ++p) {
        theDynamicalPotential += (*p)->getEnergy();
      }
      theDynamicalPotential -= getTableMass();
      theDynamicalPotential /= theA;
 
      return theDynamicalPotential;
    }

References G4INCL::Cluster::getTableMass(), G4INCL::Cluster::particles, and G4INCL::Particle::theA.

Referenced by G4INCL::Cluster::putParticlesOffShell().

◆ computeExcitationEnergy()

G4double G4INCL::Nucleus::computeExcitationEnergy ( ) const

Compute the current excitation energy.

Returns: the excitation energy [MeV]

Definition at line 281 of file G4INCLNucleus.cc.

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

References computeSeparationEnergyBalance(), computeTotalEnergy(), and initialInternalEnergy.

Referenced by G4INCL::INCL::postCascade().

◆ 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 853 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();
    theS -= remN->getS();
    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());
      }
 
    }
 
  }

Referenced by computeRecoilKinematics().

◆ computeRecoilKinematics()

void G4INCL::Nucleus::computeRecoilKinematics ( )

Compute the recoil momentum and spin of the nucleus.

Definition at line 235 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,theS) + theExcitationEnergy);
    adjustEnergyFromMomentum();
    remnant=true;
  }

References G4INCL::Particle::adjustEnergyFromMomentum(), computeCenterOfMass(), computeOneNucleonRecoilKinematics(), emitInsidePions(), G4INCL::Cluster::getAngularMomentum(), G4INCL::Particle::getMomentum(), G4INCL::Store::getOutgoingParticles(), G4INCL::ParticleTable::getTableMass, incomingAngularMomentum, incomingMomentum, initialCenterOfMass, remnant, G4INCL::Particle::setMass(), G4INCL::Particle::theA, G4INCL::Cluster::theExcitationEnergy, G4INCL::Particle::theMomentum, G4INCL::Particle::thePosition, theProjectileRemnant, G4INCL::Particle::theS, G4INCL::Cluster::theSpin, theStore, and G4INCL::Particle::theZ.

Referenced by G4INCL::INCL::postCascade().

◆ computeSeparationEnergyBalance()

G4double G4INCL::Nucleus::computeSeparationEnergyBalance ( ) const

inline

Outgoing - incoming separation energies.

Used by CDPP.

Definition at line 128 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:
          case Lambda:
          case PiPlus:
          case PiMinus:
          case KPlus:
          case KMinus:
          case KZero:
          case KZeroBar:
          case KShort:
          case KLong:
          case SigmaPlus:
          case SigmaZero:
          case SigmaMinus:
            S += thePotential->getSeparationEnergy(*i);
            break;
          case Composite:
            S += (*i)->getZ() * thePotential->getSeparationEnergy(Proton)
              + ((*i)->getA() + (*i)->getS() - (*i)->getZ()) * thePotential->getSeparationEnergy(Neutron) 
              - (*i)->getS() * thePotential->getSeparationEnergy(Lambda);
            break;
          default:
            break;
        }
      }
 
      S -= theNpInitial * thePotential->getSeparationEnergy(Proton);
      S -= theNnInitial * thePotential->getSeparationEnergy(Neutron);
      S -= theNpionplusInitial*thePotential->getSeparationEnergy(PiPlus);;
      S -= theNkaonplusInitial*thePotential->getSeparationEnergy(KPlus);
      S -= theNpionminusInitial*thePotential->getSeparationEnergy(PiMinus);
      S -= theNkaonminusInitial*thePotential->getSeparationEnergy(KMinus);
      return S;
    }

References G4INCL::Composite, G4INCL::DeltaMinus, G4INCL::DeltaPlus, G4INCL::DeltaPlusPlus, G4INCL::DeltaZero, G4INCL::Store::getOutgoingParticles(), G4INCL::NuclearPotential::INuclearPotential::getSeparationEnergy(), G4INCL::KLong, G4INCL::KMinus, G4INCL::KPlus, G4INCL::KShort, G4INCL::KZero, G4INCL::KZeroBar, G4INCL::Lambda, G4INCL::Neutron, G4INCL::PiMinus, G4INCL::PiPlus, G4INCL::Proton, S(), G4INCL::SigmaMinus, G4INCL::SigmaPlus, G4INCL::SigmaZero, theNkaonminusInitial, theNkaonplusInitial, theNnInitial, theNpInitial, theNpionminusInitial, theNpionplusInitial, thePotential, and theStore.

Referenced by computeExcitationEnergy(), and G4INCL::CDPP::isBlocked().

◆ computeTotalEnergy()

G4double G4INCL::Nucleus::computeTotalEnergy ( ) const

Compute the current total energy.

Returns: the total energy [MeV]

Definition at line 218 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 if((*p)->isHyperon())
        totalEnergy += (*p)->getEnergy() - (*p)->getPotentialEnergy() - ParticleTable::getRealMass((*p)->getType());
      else
        totalEnergy += (*p)->getEnergy() - (*p)->getPotentialEnergy();
    }
 
    return totalEnergy;
  }

References G4INCL::ParticleTable::effectiveNucleonMass, G4INCL::Store::getParticles(), G4INCL::ParticleTable::getRealMass(), and theStore.

Referenced by computeExcitationEnergy(), and initializeParticles().

◆ containsAntiKaon()

G4bool G4INCL::Nucleus::containsAntiKaon ( )

inline

Returns true if the nucleus contains any anti Kaons.

Definition at line 295 of file G4INCLNucleus.hh.

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

References G4INCL::Store::getParticles(), and theStore.

Referenced by G4INCL::INCL::postCascade().

◆ containsDeltas()

G4bool G4INCL::Nucleus::containsDeltas ( )

inline

Returns true if the nucleus contains any deltas.

Definition at line 287 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;
    }

References G4INCL::Store::getParticles(), and theStore.

Referenced by G4INCL::INCL::postCascade().

◆ containsEtas()

G4bool G4INCL::Nucleus::containsEtas ( )

inline

Returns true if the nucleus contains any etas.

Definition at line 327 of file G4INCLNucleus.hh.

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

References G4INCL::Store::getParticles(), and theStore.

◆ containsKaon()

G4bool G4INCL::Nucleus::containsKaon ( )

inline

Returns true if the nucleus contains any Kaons.

Definition at line 319 of file G4INCLNucleus.hh.

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

References G4INCL::Store::getParticles(), and theStore.

Referenced by G4INCL::INCL::postCascade().

◆ containsLambda()

G4bool G4INCL::Nucleus::containsLambda ( )

inline

Returns true if the nucleus contains any Lambda.

Definition at line 303 of file G4INCLNucleus.hh.

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

References G4INCL::Store::getParticles(), and theStore.

Referenced by G4INCL::INCL::postCascade().

◆ containsOmegas()

G4bool G4INCL::Nucleus::containsOmegas ( )

inline

Returns true if the nucleus contains any omegas.

Definition at line 335 of file G4INCLNucleus.hh.

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

References G4INCL::Store::getParticles(), and theStore.

◆ containsSigma()

G4bool G4INCL::Nucleus::containsSigma ( )

inline

Returns true if the nucleus contains any Sigma.

Definition at line 311 of file G4INCLNucleus.hh.

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

References G4INCL::Store::getParticles(), and theStore.

Referenced by G4INCL::INCL::postCascade().

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

References applyFinalState(), G4INCL::ClusterDecay::decay(), emitInsidePions(), G4INCL::Store::getParticles(), G4INCL::FinalState::getValidity(), G4INCL::NuclearPotential::INuclearPotential::hasPionPotential(), INCL_DEBUG, INCL_WARN, G4INCL::Particle::theA, thePotential, theStore, G4INCL::Particle::theZ, and G4INCL::ValidFS.

Referenced by G4INCL::INCL::postCascade().

◆ decayInsideStrangeParticles()

G4bool G4INCL::Nucleus::decayInsideStrangeParticles ( )

Force the transformation of strange particles into a Lambda;.

Returns: true if any strange particles was forced to absorb.

Definition at line 417 of file G4INCLNucleus.cc.

                                              {
     
    /* Transform each strange particles into a lambda
     * Every Kaon (KPlus and KZero) are emited
     */
    const G4bool unphysicalRemnant = (theZ<0 || theZ>theA);
    if(unphysicalRemnant){
        emitInsideStrangeParticles();
        INCL_WARN("Remnant is unphysical: Z=" << theZ << ", A=" << theA << ", too much strange particles? -> all emit" << '\n');
        return false;
    }
 
    /* Build a list of particles with a strangeness == -1 except Lambda,
     * and two other one for proton and neutron
     */
    ParticleList const &inside = theStore->getParticles();
    ParticleList stranges;
    ParticleList protons;
    ParticleList neutrons;
    for(ParticleIter i=inside.begin(), e=inside.end(); i!=e; ++i){
        if((*i)->isSigma() || (*i)->isAntiKaon()) stranges.push_back((*i));
        else if((*i)->isNucleon() && (*i)->getZ() == 1) protons.push_back((*i));
        else if((*i)->isNucleon() && (*i)->getZ() == 0) neutrons.push_back((*i));
    }
    
    if((stranges.size() > protons.size()) || (stranges.size() > neutrons.size())){
        INCL_WARN("Remnant is unphysical: Nproton=" << protons.size() << ", Nneutron=" << neutrons.size() << ", Strange particles : " << stranges.size() <<  '\n');
        emitInsideStrangeParticles();
        return false;
    }
    
    // Loop over the strange particles, make them absorbe
    ParticleIter protonIter = protons.begin();
    ParticleIter neutronIter = neutrons.begin();
    for(ParticleIter i=stranges.begin(), e=stranges.end(); i!=e; ++i) {
      INCL_DEBUG("Absorbe inside strange particles:" << '\n'
          << (*i)->print() << '\n');
      IAvatar *decay;
      if((*i)->getType() == SigmaMinus){
          decay = new DecayAvatar((*protonIter), (*i), 0.0, this, true);
          ++protonIter;
      }
      else if((*i)->getType() == SigmaPlus){
          decay = new DecayAvatar((*neutronIter), (*i), 0.0, this, true);
          ++neutronIter;
      }
      else if(Random::shoot()*(protons.size() + neutrons.size()) < protons.size()){
          decay = new DecayAvatar((*protonIter), (*i), 0.0, this, true);
          ++protonIter;
      }
      else {
          decay = new DecayAvatar((*neutronIter), (*i), 0.0, this, true);
          ++neutronIter;
      }
      FinalState *fs = decay->getFinalState();
      applyFinalState(fs);
      delete fs;
      delete decay;
    }
 
    return true;
  }

References applyFinalState(), G4INCL::ClusterDecay::decay(), emitInsideStrangeParticles(), G4INCL::Store::getParticles(), INCL_DEBUG, INCL_WARN, G4INCL::Random::shoot(), G4INCL::SigmaMinus, G4INCL::SigmaPlus, G4INCL::Particle::theA, theStore, and G4INCL::Particle::theZ.

Referenced by G4INCL::INCL::postCascade().

◆ decayMe()

G4bool G4INCL::Nucleus::decayMe ( )

Force the phase-space decay of the Nucleus.

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

Returns: true if the nucleus was forced to decay.

Definition at line 656 of file G4INCLNucleus.cc.

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

References G4INCL::Store::addToOutgoing(), G4INCL::ClusterDecay::decay(), G4INCL::Particle::getBiasCollisionVector(), G4INCL::Particle::theA, G4INCL::Particle::theS, theStore, and G4INCL::Particle::theZ.

Referenced by G4INCL::INCL::makeCompoundNucleus(), and G4INCL::INCL::postCascade().

◆ decayOutgoingClusters()

G4bool G4INCL::Nucleus::decayOutgoingClusters ( )

Force the decay of unstable outgoing clusters.

Returns: true if any cluster was forced to decay.

Definition at line 631 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){
        (*j)->setBiasCollisionVector(cluster->getBiasCollisionVector());
        theStore->addToOutgoing(*j);
      }
    }
    return true;
  }

References G4INCL::Store::addToOutgoing(), G4INCL::ClusterDecay::decay(), G4INCL::Cluster::deleteParticles(), G4INCL::Particle::getBiasCollisionVector(), G4INCL::Store::getOutgoingParticles(), and theStore.

Referenced by G4INCL::INCL::makeCompoundNucleus(), and G4INCL::INCL::postCascade().

◆ decayOutgoingDeltas()

G4bool G4INCL::Nucleus::decayOutgoingDeltas ( )

Force the decay of outgoing deltas.

Returns: true if any delta was forced to decay.

Definition at line 308 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);
      pion->setBiasCollisionVector(nucleon->getBiasCollisionVector());
 
      nucleon->setTableMass();
      nucleon->setMomentum(-newMomentum);
      nucleon->adjustEnergyFromMomentum();
      nucleon->boost(beta);
 
      theStore->addToOutgoing(pion);
 
      delete fs;
      delete decay;
    }
 
    return true;
  }

References G4INCL::Store::addToOutgoing(), anonymous_namespace{G4PionRadiativeDecayChannel.cc}::beta, G4INCL::ClusterDecay::decay(), G4INCL::FinalState::getCreatedParticles(), G4INCL::FinalState::getModifiedParticles(), G4INCL::Store::getOutgoingParticles(), INCL_DEBUG, G4INCL::ThreeVector::mag(), G4INCL::KinematicsUtils::momentumInCM(), G4InuclParticleNames::nucleon(), G4InuclParticleNames::pion(), and theStore.

Referenced by G4INCL::INCL::makeCompoundNucleus(), and G4INCL::INCL::postCascade().

◆ decayOutgoingNeutralKaon()

G4bool G4INCL::Nucleus::decayOutgoingNeutralKaon ( )

Force the transformation of outgoing Neutral Kaon into propation eigenstate.

Returns: true if any kaon was forced to decay.

Definition at line 608 of file G4INCLNucleus.cc.

                                           {
        ParticleList const &out = theStore->getOutgoingParticles();
        ParticleList neutralkaon;
        for(ParticleIter i=out.begin(), e=out.end(); i!=e; ++i) {
            if((*i)->getType() == KZero  || (*i)->getType() == KZeroBar)  neutralkaon.push_back((*i));
        }
        if(neutralkaon.empty()) return false;
        
        for(ParticleIter i=neutralkaon.begin(), e=neutralkaon.end(); i!=e; ++i) {
            INCL_DEBUG("Transform outgoing neutral kaon:" << '\n'
                       << (*i)->print() << '\n');
            
            // Use a DecayAvatar
            IAvatar *decay = new DecayAvatar((*i), 0.0, NULL);
            FinalState *fs = decay->getFinalState();
            
            delete fs;
            delete decay;
        }
        
        return true;
    }

References G4INCL::ClusterDecay::decay(), G4INCL::Store::getOutgoingParticles(), INCL_DEBUG, G4INCL::KZero, G4INCL::KZeroBar, and theStore.

Referenced by G4INCL::INCL::postCascade().

◆ decayOutgoingPionResonances()

G4bool G4INCL::Nucleus::decayOutgoingPionResonances ( G4double timeThreshold )

Force the decay of outgoing PionResonances (eta/omega).

Returns: true if any eta was forced to decay.

Definition at line 480 of file G4INCLNucleus.cc.

                                                                    {
        ParticleList const &out = theStore->getOutgoingParticles();
        ParticleList pionResonances;
        for(ParticleIter i=out.begin(), e=out.end(); i!=e; ++i) {
//            if((*i)->isEta() || (*i)->isOmega()) pionResonances.push_back((*i));
            if(((*i)->isEta() && timeThreshold > ParticleTable::getWidth(Eta)) || ((*i)->isOmega() && timeThreshold > ParticleTable::getWidth(Omega))) pionResonances.push_back((*i));
        }
        if(pionResonances.empty()) return false;
        
        for(ParticleIter i=pionResonances.begin(), e=pionResonances.end(); i!=e; ++i) {
            INCL_DEBUG("Decay outgoing pionResonances particle:" << '\n'
                       << (*i)->print() << '\n');
            const ThreeVector beta = -(*i)->boostVector();
            const G4double pionResonanceMass = (*i)->getMass();
            
            // Set the pionResonance 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 theModifiedParticle = fs->getModifiedParticles().front();
            ParticleList const &created = fs->getCreatedParticles();
            Particle * const theCreatedParticle1 = created.front();
                            
            if (created.size() == 1) {
                
                // Adjust the decay momentum if we are using the real masses
                const G4double decayMomentum = KinematicsUtils::momentumInCM(pionResonanceMass,theModifiedParticle->getTableMass(),theCreatedParticle1->getTableMass());
                ThreeVector newMomentum = theCreatedParticle1->getMomentum();
                newMomentum *= decayMomentum / newMomentum.mag();
                
                theCreatedParticle1->setTableMass();
                theCreatedParticle1->setMomentum(newMomentum);
                theCreatedParticle1->adjustEnergyFromMomentum();
                //theCreatedParticle1->setEmissionTime(nucleon->getEmissionTime());
                theCreatedParticle1->boost(beta);
                theCreatedParticle1->setBiasCollisionVector(theModifiedParticle->getBiasCollisionVector());
                
                theModifiedParticle->setTableMass();
                theModifiedParticle->setMomentum(-newMomentum);
                theModifiedParticle->adjustEnergyFromMomentum();
                theModifiedParticle->boost(beta);
                
                theStore->addToOutgoing(theCreatedParticle1);
            }
            else if (created.size() == 2) {
                Particle * const theCreatedParticle2 = created.back();
                
                theCreatedParticle1->boost(beta);
                theCreatedParticle1->setBiasCollisionVector(theModifiedParticle->getBiasCollisionVector());
                theCreatedParticle2->boost(beta);
                theCreatedParticle2->setBiasCollisionVector(theModifiedParticle->getBiasCollisionVector());
                theModifiedParticle->boost(beta);
                
                theStore->addToOutgoing(theCreatedParticle1);
                theStore->addToOutgoing(theCreatedParticle2);
            }
            else {
                INCL_ERROR("Wrong number (< 2) of created particles during the decay of a pion resonance");
            }
            delete fs;
            delete decay;
        }
        
        return true;
    }

References G4INCL::Store::addToOutgoing(), G4INCL::Particle::adjustEnergyFromMomentum(), anonymous_namespace{G4PionRadiativeDecayChannel.cc}::beta, G4INCL::Particle::boost(), G4INCL::ClusterDecay::decay(), G4INCL::Eta, G4INCL::Particle::getBiasCollisionVector(), G4INCL::FinalState::getCreatedParticles(), G4INCL::FinalState::getModifiedParticles(), G4INCL::Particle::getMomentum(), G4INCL::Store::getOutgoingParticles(), G4INCL::Particle::getTableMass(), G4INCL::ParticleTable::getWidth(), INCL_DEBUG, INCL_ERROR, G4INCL::ThreeVector::mag(), G4INCL::KinematicsUtils::momentumInCM(), G4INCL::Omega, G4INCL::Particle::setBiasCollisionVector(), G4INCL::Particle::setMomentum(), G4INCL::Particle::setTableMass(), and theStore.

Referenced by G4INCL::INCL::makeCompoundNucleus(), and G4INCL::INCL::postCascade().

◆ decayOutgoingSigmaZero()

G4bool G4INCL::Nucleus::decayOutgoingSigmaZero ( G4double timeThreshold )

Force the decay of outgoing Neutral Sigma.

Returns: true if any Sigma was forced to decay.

Definition at line 551 of file G4INCLNucleus.cc.

                                                               {
        ParticleList const &out = theStore->getOutgoingParticles();
        ParticleList neutralsigma;
        for(ParticleIter i=out.begin(), e=out.end(); i!=e; ++i) {
            if((*i)->getType() == SigmaZero && timeThreshold > ParticleTable::getWidth(SigmaZero))  neutralsigma.push_back((*i));
        }
        if(neutralsigma.empty()) return false;
        
        for(ParticleIter i=neutralsigma.begin(), e=neutralsigma.end(); i!=e; ++i) {
            INCL_DEBUG("Decay outgoing neutral sigma:" << '\n'
                       << (*i)->print() << '\n');
            const ThreeVector beta = -(*i)->boostVector();
            const G4double neutralsigmaMass = (*i)->getMass();
            
            // Set the neutral sigma 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 theModifiedParticle = fs->getModifiedParticles().front();
            ParticleList const &created = fs->getCreatedParticles();
            Particle * const theCreatedParticle = created.front();
                            
            if (created.size() == 1) {
                
                // Adjust the decay momentum if we are using the real masses
                const G4double decayMomentum = KinematicsUtils::momentumInCM(neutralsigmaMass,theModifiedParticle->getTableMass(),theCreatedParticle->getTableMass());
                ThreeVector newMomentum = theCreatedParticle->getMomentum();
                newMomentum *= decayMomentum / newMomentum.mag();
                
                theCreatedParticle->setTableMass();
                theCreatedParticle->setMomentum(newMomentum);
                theCreatedParticle->adjustEnergyFromMomentum();
                theCreatedParticle->boost(beta);
                theCreatedParticle->setBiasCollisionVector(theModifiedParticle->getBiasCollisionVector());
                
                theModifiedParticle->setTableMass();
                theModifiedParticle->setMomentum(-newMomentum);
                theModifiedParticle->adjustEnergyFromMomentum();
                theModifiedParticle->boost(beta);
                
                theStore->addToOutgoing(theCreatedParticle);
            }
            else {
                INCL_ERROR("Wrong number (!= 1) of created particles during the decay of a sigma zero");
            }
            delete fs;
            delete decay;
        }
        
        return true;
    }

References G4INCL::Store::addToOutgoing(), G4INCL::Particle::adjustEnergyFromMomentum(), anonymous_namespace{G4PionRadiativeDecayChannel.cc}::beta, G4INCL::Particle::boost(), G4INCL::ClusterDecay::decay(), G4INCL::Particle::getBiasCollisionVector(), G4INCL::FinalState::getCreatedParticles(), G4INCL::FinalState::getModifiedParticles(), G4INCL::Particle::getMomentum(), G4INCL::Store::getOutgoingParticles(), G4INCL::Particle::getTableMass(), G4INCL::ParticleTable::getWidth(), INCL_DEBUG, INCL_ERROR, G4INCL::ThreeVector::mag(), G4INCL::KinematicsUtils::momentumInCM(), G4INCL::Particle::setBiasCollisionVector(), G4INCL::Particle::setMomentum(), G4INCL::Particle::setTableMass(), G4INCL::SigmaZero, and theStore.

Referenced by G4INCL::INCL::postCascade().

◆ deleteParticles()

void G4INCL::Cluster::deleteParticles ( )

inlineinherited

Definition at line 136 of file G4INCLCluster.hh.

                           {
      for(ParticleIter p=particles.begin(), e=particles.end(); p!=e; ++p) {
        delete (*p);
      }
      clearParticles();
    }

References G4INCL::Cluster::clearParticles(), and G4INCL::Cluster::particles.

Referenced by G4INCL::Store::clearOutgoing(), decayOutgoingClusters(), G4INCL::ProjectileRemnant::reset(), and G4INCL::ProjectileRemnant::~ProjectileRemnant().

◆ deleteProjectileRemnant()

void G4INCL::Nucleus::deleteProjectileRemnant ( )

inline

Delete the projectile remnant.

Definition at line 442 of file G4INCLNucleus.hh.

                                   {
      delete theProjectileRemnant;
      theProjectileRemnant = NULL;
    }

References theProjectileRemnant.

Referenced by ~Nucleus().

◆ dump()

std::string G4INCL::Particle::dump ( ) const

inlineinherited

Definition at line 968 of file G4INCLParticle.hh.

                           {
      std::stringstream ss;
      ss << "(particle " << ID << " ";
      ss << ParticleTable::getName(theType);
      ss << '\n'
        << thePosition.dump()
        << '\n'
        << theMomentum.dump()
        << '\n'
        << theEnergy << ")" << '\n';
      return ss.str();
    };

References G4INCL::ThreeVector::dump(), G4INCL::ParticleTable::getName(), G4INCL::Particle::ID, G4INCL::Particle::theEnergy, G4INCL::Particle::theMomentum, G4INCL::Particle::thePosition, and G4INCL::Particle::theType.

Referenced by G4INCL::BinaryCollisionAvatar::dump(), G4INCL::DecayAvatar::dump(), G4INCL::ParticleEntryAvatar::dump(), and G4INCL::SurfaceAvatar::dump().

◆ emitInsideKaon()

G4bool G4INCL::Nucleus::emitInsideKaon ( )

Force emission of all Kaon inside the nucleus.

Definition at line 796 of file G4INCLNucleus.cc.

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

Referenced by G4INCL::INCL::makeCompoundNucleus(), and G4INCL::INCL::postCascade().

◆ emitInsideLambda()

G4int G4INCL::Nucleus::emitInsideLambda ( )

Force emission of all Lambda (desexitation code with strangeness not implanted yet)

Definition at line 753 of file G4INCLNucleus.cc.

                                  {
    /* Forcing emissions of all Lambda in the nucleus.
     * This probably violates energy conservation
     * (although the computation of the recoil kinematics
     * might sweep this under the carpet).
     */
    INCL_DEBUG("Forcing emissions of all Lambda in the nucleus." << '\n');
 
    // Emit the Lambda with this kinetic energy
    const G4double tinyEnergy = 0.1; // MeV
 
    // Push out the emitted Lambda
    ParticleList const &inside = theStore->getParticles();
    ParticleList toEject;
    for(ParticleIter i=inside.begin(), e=inside.end(); i!=e; ++i) {
      if((*i)->isLambda()) {
        Particle * const theLambda = *i;
        INCL_DEBUG("Forcing emission of the following particle: "
                   << theLambda->print() << '\n');
        theLambda->setEmissionTime(theStore->getBook().getCurrentTime());
        // Correction for real masses
        const G4double theQValueCorrection = theLambda->getEmissionQValueCorrection(theA,theZ,theS); // Does it work for strange particles? Should be check
        const G4double kineticEnergyOutside = theLambda->getKineticEnergy() - theLambda->getPotentialEnergy() + theQValueCorrection;
        theLambda->setTableMass();
        if(kineticEnergyOutside > 0.0)
          theLambda->setEnergy(theLambda->getMass()+kineticEnergyOutside);
        else
          theLambda->setEnergy(theLambda->getMass()+tinyEnergy);
        theLambda->adjustMomentumFromEnergy();
        theLambda->setPotentialEnergy(0.);
        theA -= theLambda->getA();
        theS -= theLambda->getS();
        toEject.push_back(theLambda);
      }
    }
    for(ParticleIter i=toEject.begin(), e=toEject.end(); i!=e; ++i) {
      theStore->particleHasBeenEjected(*i);
      theStore->addToOutgoing(*i);
      (*i)->setParticleBias(Particle::getTotalBias());
    }
    return toEject.size();
  }

Referenced by G4INCL::INCL::postCascade().

◆ emitInsidePions()

void G4INCL::Nucleus::emitInsidePions ( )

Force emission of all pions inside the nucleus.

Definition at line 670 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,theS);
        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);
      (*i)->setParticleBias(Particle::getTotalBias());
    }
  }

Referenced by computeRecoilKinematics(), and decayInsideDeltas().

◆ emitInsideStrangeParticles()

void G4INCL::Nucleus::emitInsideStrangeParticles ( )

Force emission of all strange particles inside the nucleus.

Definition at line 710 of file G4INCLNucleus.cc.

                                           {
    /* Forcing emissions of Sigmas and antiKaons.
     * This probably violates energy conservation
     * (although the computation of the recoil kinematics
     * might sweep this under the carpet).
     */
    INCL_DEBUG("Forcing emissions of all strange particles in the nucleus." << '\n');
 
    // Emit the strange particles with this kinetic energy
    const G4double tinyEnergy = 0.1; // MeV
 
    // Push out the emitted strange particles
    ParticleList const &inside = theStore->getParticles();
    ParticleList toEject;
    for(ParticleIter i=inside.begin(), e=inside.end(); i!=e; ++i) {
      if((*i)->isSigma() || (*i)->isAntiKaon()) {
        Particle * const theParticle = *i;
        INCL_DEBUG("Forcing emission of the following particle: "
                   << theParticle->print() << '\n');
        theParticle->setEmissionTime(theStore->getBook().getCurrentTime());
        // Correction for real masses
        const G4double theQValueCorrection = theParticle->getEmissionQValueCorrection(theA,theZ,theS); // Does it work for strange particles? should be check
        const G4double kineticEnergyOutside = theParticle->getKineticEnergy() - theParticle->getPotentialEnergy() + theQValueCorrection;
        theParticle->setTableMass();
        if(kineticEnergyOutside > 0.0)
          theParticle->setEnergy(theParticle->getMass()+kineticEnergyOutside);
        else
          theParticle->setEnergy(theParticle->getMass()+tinyEnergy);
        theParticle->adjustMomentumFromEnergy();
        theParticle->setPotentialEnergy(0.);
        theA -= theParticle->getA();
        theZ -= theParticle->getZ();
        theS -= theParticle->getS();
        toEject.push_back(theParticle);
      }
    }
    for(ParticleIter i=toEject.begin(), e=toEject.end(); i!=e; ++i) {
      theStore->particleHasBeenEjected(*i);
      theStore->addToOutgoing(*i);
      (*i)->setParticleBias(Particle::getTotalBias());
    }
  }

Referenced by decayInsideStrangeParticles(), and G4INCL::INCL::postCascade().

◆ fillEventInfo()

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

Fill the event info which contains INCL output data

Definition at line 1010 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->ParticleBias[eventInfo->nParticles] = (*i)->getParticleBias();
      
      eventInfo->A[eventInfo->nParticles] = (*i)->getA();
      eventInfo->Z[eventInfo->nParticles] = (*i)->getZ();
      eventInfo->S[eventInfo->nParticles] = (*i)->getS();
      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("");
      if ((*i)->getType() != Composite) {
        ParticleSpecies pt((*i)->getType());
        eventInfo->PDGCode[eventInfo->nParticles] = pt.getPDGCode();
      }
      else {
        ParticleSpecies pt((*i)->getA(), (*i)->getZ(), (*i)->getS());
        eventInfo->PDGCode[eventInfo->nParticles] = pt.getPDGCode();
      }
      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();
      eventInfo->SRem[eventInfo->nRemnants] = theProjectileRemnant->getS();
      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->SRem[eventInfo->nRemnants] = getS();
      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] << " eventNumber=" << eventInfo->eventNumber << '\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();
  }

Referenced by G4INCL::INCL::makeCompoundNucleus(), and G4INCL::INCL::postCascade().

◆ FillINCLBiasVector()

void G4INCL::Particle::FillINCLBiasVector ( G4double newBias )

staticinherited

Definition at line 202 of file G4INCLParticle.cc.

                                                   {
// assert(G4int(Particle::INCLBiasVector.size())==nextBiasedCollisionID);
    //assert(G4int(Particle::INCLBiasVector.Size())==nextBiasedCollisionID);
// assert(std::fabs(newBias - 1.) > 1E-6);
    Particle::INCLBiasVector.push_back(newBias);
    //Particle::INCLBiasVector.Push_back(newBias);
    Particle::nextBiasedCollisionID++;
  }

References G4INCL::Particle::INCLBiasVector, and G4INCL::Particle::nextBiasedCollisionID.

Referenced by G4INCL::InteractionAvatar::postInteraction().

◆ 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 1211 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(), theProjectileRemnant->getS());
 
      // Set the excitation energy
      theProjectileRemnant->setExcitationEnergy(anExcitationEnergy);
 
      // No spin!
      theProjectileRemnant->setSpin(ThreeVector());
 
      // Set the emission time
      theProjectileRemnant->setEmissionTime(anEmissionTime);
    }
  }

References G4INCL::Particle::getA(), G4INCL::Particle::getInvariantMass(), G4INCL::Particle::getS(), G4INCL::ParticleTable::getTableMass, G4INCL::Particle::getZ(), G4INCL::Particle::setEmissionTime(), G4INCL::Cluster::setExcitationEnergy(), G4INCL::Particle::setMass(), G4INCL::Cluster::setSpin(), and theProjectileRemnant.

Referenced by G4INCL::INCL::makeCompoundNucleus(), and G4INCL::INCL::makeProjectileRemnant().

◆ freezeInternalMotion()

void G4INCL::Cluster::freezeInternalMotion ( )

inlineinherited

Freeze the internal motion of the particles.

Each particle is assigned a frozen momentum four-vector determined by the collective cluster velocity. This is used for propagation, but not for dynamics. Normal propagation is restored by calling the Particle::thawPropagation() method, which should be done in InteractionAvatar::postInteraction.

Definition at line 367 of file G4INCLCluster.hh.

                                {
      const ThreeVector &normMomentum = theMomentum / getMass();
      for(ParticleIter p=particles.begin(), e=particles.end(); p!=e; ++p) {
        const G4double pMass = (*p)->getMass();
        const ThreeVector frozenMomentum = normMomentum * pMass;
        const G4double frozenEnergy = std::sqrt(frozenMomentum.mag2()+pMass*pMass);
        (*p)->setFrozenMomentum(frozenMomentum);
        (*p)->setFrozenEnergy(frozenEnergy);
        (*p)->freezePropagation();
      }
    }

References G4INCL::Particle::getMass(), G4INCL::ThreeVector::mag2(), G4INCL::Cluster::particles, and G4INCL::Particle::theMomentum.

Referenced by G4INCL::ProjectileRemnant::ProjectileRemnant().

◆ freezePropagation()

void G4INCL::Particle::freezePropagation ( )

inlineinherited

Freeze particle propagation.

Make the particle use theFrozenMomentum and theFrozenEnergy for propagation. The normal state can be restored by calling the thawPropagation() method.

Definition at line 908 of file G4INCLParticle.hh.

                             {
      thePropagationMomentum = &theFrozenMomentum;
      thePropagationEnergy = &theFrozenEnergy;
    }

References G4INCL::Particle::theFrozenEnergy, G4INCL::Particle::theFrozenMomentum, G4INCL::Particle::thePropagationEnergy, and G4INCL::Particle::thePropagationMomentum.

◆ getA()

G4int G4INCL::Particle::getA ( ) const

inlineinherited

Returns the baryon number.

Definition at line 393 of file G4INCLParticle.hh.

393{ return theA; }

References G4INCL::Particle::theA.

Referenced by G4INCL::ProjectileRemnant::addDynamicalSpectator(), G4INCL::Cluster::addParticle(), G4INCL::ClusteringModelIntercomparison::clusterCanEscape(), computeOneNucleonRecoilKinematics(), G4INCL::INCL::continueCascade(), G4INCL::ClusterDecay::decay(), emitInsideLambda(), emitInsideStrangeParticles(), fillEventInfo(), G4INCL::ParticleEntryChannel::fillFinalState(), G4INCL::TransmissionChannel::fillFinalState(), finalizeProjectileRemnant(), G4INCL::ClusteringModelIntercomparison::findClusterStartingFrom(), G4INCL::PauliStandard::getBlockingProbability(), G4INCL::ClusteringModelIntercomparison::getCluster(), getConservationBalance(), G4INCL::SurfaceAvatar::getTransmissionProbability(), G4INCL::NuclearDensity::getTransmissionRadius(), G4INCL::TransmissionChannel::initializeKineticEnergyOutside(), insertParticle(), G4INCL::ClusterDecay::isStable(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::phaseSpaceDecay(), G4INCL::INCL::postCascade(), G4INCL::INCL::preCascade(), G4INCL::INCL::RecoilCMFunctor::RecoilCMFunctor(), G4INCL::INCL::RecoilFunctor::RecoilFunctor(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::recursiveDecay(), G4INCL::ProjectileRemnant::removeParticle(), G4INCL::INCL::RecoilCMFunctor::scaleParticleCMMomenta(), G4INCL::INCL::RecoilFunctor::scaleParticleEnergies(), G4INCL::InteractionAvatar::ViolationEMomentumFunctor::scaleParticleMomenta(), G4INCL::StandardPropagationModel::shootComposite(), G4INCL::StandardPropagationModel::shootParticle(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::threeBodyDecay(), and G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::twoBodyDecay().

◆ getAngularMomentum()

G4INCL::ThreeVector G4INCL::Cluster::getAngularMomentum ( ) const

inlinevirtualinherited

Get the total angular momentum (orbital + spin)

Reimplemented from G4INCL::Particle.

Definition at line 428 of file G4INCLCluster.hh.

                                                 {
      return Particle::getAngularMomentum() + getSpin();
    }

References G4INCL::Particle::getAngularMomentum(), and G4INCL::Cluster::getSpin().

Referenced by computeRecoilKinematics(), G4INCL::INCL::makeCompoundNucleus(), and G4INCL::StandardPropagationModel::shootComposite().

◆ getBeta()

G4double G4INCL::Particle::getBeta ( ) const

inlineinherited

Definition at line 401 of file G4INCLParticle.hh.

                             {
      const G4double P = theMomentum.mag();
      return P/theEnergy;
    }

References G4INCL::ThreeVector::mag(), P, G4INCL::Particle::theEnergy, and G4INCL::Particle::theMomentum.

◆ getBiasCollisionVector()

std::vector< G4int > G4INCL::Particle::getBiasCollisionVector ( ) const

inlineinherited

Get the vector list of biased vertices on the particle path.

Definition at line 1038 of file G4INCLParticle.hh.

1038{ return theBiasCollisionVector; }

G4INCL::Particle::theBiasCollisionVector

std::vector< G4int > theBiasCollisionVector

Time ordered vector of all biased vertices on the particle path.

Definition: G4INCLParticle.hh:1108

References G4INCL::Particle::theBiasCollisionVector.

Referenced by G4INCL::ClusterDecay::decay(), decayMe(), decayOutgoingClusters(), decayOutgoingPionResonances(), decayOutgoingSigmaZero(), and G4INCL::Particle::MergeVectorBias().

◆ getBiasFromVector()

G4double G4INCL::Particle::getBiasFromVector ( std::vector< G4int > VectorBias )

staticinherited

Definition at line 211 of file G4INCLParticle.cc.

                                                                  {
    if(VectorBias.empty()) return 1.;
    
    G4double ParticleBias = 1.;
    
    for(G4int i=0; i<G4int(VectorBias.size()); i++){
        ParticleBias *= Particle::INCLBiasVector[G4int(VectorBias[i])];
    }
    
    return ParticleBias;
  }

References G4INCL::Particle::INCLBiasVector.

Referenced by G4INCL::BinaryCollisionAvatar::getChannel(), G4INCL::ParticleList::getParticleListBias(), and G4INCL::Particle::setBiasCollisionVector().

◆ getConservationBalance()

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

Compute charge, mass, energy and momentum balance.

Definition at line 1155 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.S = theEventInfo.Sp + theEventInfo.St;
 
    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();
      theBalance.S -= (*i)->getS();
      // 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.S -= theProjectileRemnant->getS();
      theBalance.energy -= ParticleTable::getTableMass(theProjectileRemnant->getA(),theProjectileRemnant->getZ(),theProjectileRemnant->getS()) +
        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.S -= getS();
      theBalance.energy -= ParticleTable::getTableMass(getA(),getZ(),getS()) +
        getExcitationEnergy();
      if(afterRecoil)
        theBalance.energy -= getKineticEnergy();
      theBalance.momentum -= getMomentum();
    }
 
    return theBalance;
  }

References G4INCL::Nucleus::ConservationBalance::A, G4INCL::EventInfo::Ap, G4INCL::EventInfo::At, G4INCL::Nucleus::ConservationBalance::energy, G4INCL::Particle::getA(), G4INCL::Cluster::getExcitationEnergy(), getExcitationEnergy(), getIncomingMomentum(), getInitialEnergy(), G4INCL::Particle::getKineticEnergy(), G4INCL::Particle::getMomentum(), G4INCL::Store::getOutgoingParticles(), G4INCL::Particle::getS(), G4INCL::ParticleTable::getTableMass, G4INCL::Particle::getZ(), hasRemnant(), G4INCL::Nucleus::ConservationBalance::momentum, G4INCL::Nucleus::ConservationBalance::S, G4INCL::EventInfo::Sp, G4INCL::EventInfo::St, theProjectileRemnant, theStore, G4INCL::Nucleus::ConservationBalance::Z, G4INCL::EventInfo::Zp, and G4INCL::EventInfo::Zt.

Referenced by G4INCL::INCL::RecoilFunctor::operator()(), and G4INCL::INCL::RecoilCMFunctor::operator()().

◆ getCosRPAngle()

G4double G4INCL::Particle::getCosRPAngle ( ) const

inlineinherited

Get the cosine of the angle between position and momentum.

Definition at line 1014 of file G4INCLParticle.hh.

                                   {
      const G4double norm = thePosition.mag2()*thePropagationMomentum->mag2();
      if(norm>0.)
        return thePosition.dot(*thePropagationMomentum) / std::sqrt(norm);
      else
        return 1.;
    }

References G4INCL::ThreeVector::dot(), G4INCL::ThreeVector::mag2(), G4INCL::Particle::thePosition, and G4INCL::Particle::thePropagationMomentum.

Referenced by G4INCL::SurfaceAvatar::initializeRefractionVariables(), and G4INCL::ParticleEntryChannel::particleEnters().

◆ getDensity()

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

inline

Getter for theDensity.

Definition at line 470 of file G4INCLNucleus.hh.

470{ return theDensity; };

References theDensity.

Referenced by G4INCL::CoulombNonRelativistic::distortOut(), G4INCL::ClusteringModelIntercomparison::findClusterStartingFrom(), and G4INCL::ClusteringModelIntercomparison::getCluster().

◆ getEmissionQValueCorrection() [1/2]

G4double G4INCL::Particle::getEmissionQValueCorrection	(	const G4int	AParent,
		const G4int	ZParent
	)		const

inlineinherited

Computes correction on the emission Q-value.

Computes the correction that must be applied to INCL particles in order to obtain the correct Q-value for particle emission from a given nucleus. For absorption, the correction is obviously equal to minus the value returned by this function.

Parameters

AParent	the mass number of the emitting nucleus
ZParent	the charge number of the emitting nucleus

Returns: the correction

Definition at line 602 of file G4INCLParticle.hh.

                                                                                         {
      const G4int SParent = 0;
      const G4int ADaughter = AParent - theA;
      const G4int ZDaughter = ZParent - theZ;
      const G4int SDaughter = 0;
 
      // Note the minus sign here
      G4double theQValue;
      if(isCluster())
        theQValue = -ParticleTable::getTableQValue(theA, theZ, theS, ADaughter, ZDaughter, SDaughter);
      else {
        const G4double massTableParent = ParticleTable::getTableMass(AParent,ZParent,SParent);
        const G4double massTableDaughter = ParticleTable::getTableMass(ADaughter,ZDaughter,SDaughter);
        const G4double massTableParticle = getTableMass();
        theQValue = massTableParent - massTableDaughter - massTableParticle;
      }
 
      const G4double massINCLParent = ParticleTable::getINCLMass(AParent,ZParent,SParent);
      const G4double massINCLDaughter = ParticleTable::getINCLMass(ADaughter,ZDaughter,SDaughter);
      const G4double massINCLParticle = getINCLMass();
 
      // The rhs corresponds to the INCL Q-value
      return theQValue - (massINCLParent-massINCLDaughter-massINCLParticle);
    }

References G4INCL::Particle::getINCLMass(), G4INCL::ParticleTable::getINCLMass(), G4INCL::Particle::getTableMass(), G4INCL::ParticleTable::getTableMass, G4INCL::ParticleTable::getTableQValue(), G4INCL::Particle::isCluster(), G4INCL::Particle::theA, G4INCL::Particle::theS, and G4INCL::Particle::theZ.

Referenced by emitInsideKaon(), emitInsideLambda(), emitInsidePions(), emitInsideStrangeParticles(), G4INCL::ParticleEntryChannel::fillFinalState(), G4INCL::TransmissionChannel::fillFinalState(), G4INCL::SurfaceAvatar::getTransmissionProbability(), and G4INCL::TransmissionChannel::initializeKineticEnergyOutside().

◆ getEmissionQValueCorrection() [2/2]

G4double G4INCL::Particle::getEmissionQValueCorrection	(	const G4int	AParent,
		const G4int	ZParent,
		const G4int	SParent
	)		const

inlineinherited

Computes correction on the emission Q-value for hypernuclei.

Computes the correction that must be applied to INCL particles in order to obtain the correct Q-value for particle emission from a given nucleus. For absorption, the correction is obviously equal to minus the value returned by this function.

Parameters

AParent	the mass number of the emitting nucleus
ZParent	the charge number of the emitting nucleus
SParent	the strangess number of the emitting nucleus

Returns: the correction

Definition at line 678 of file G4INCLParticle.hh.

                                                                                                              {
      const G4int ADaughter = AParent - theA;
      const G4int ZDaughter = ZParent - theZ;
      const G4int SDaughter = SParent - theS;
 
      // Note the minus sign here
      G4double theQValue;
      if(isCluster())
        theQValue = -ParticleTable::getTableQValue(theA, theZ, theS, ADaughter, ZDaughter, SDaughter);
      else {
        const G4double massTableParent = ParticleTable::getTableMass(AParent,ZParent,SParent);
        const G4double massTableDaughter = ParticleTable::getTableMass(ADaughter,ZDaughter,SDaughter);
        const G4double massTableParticle = getTableMass();
        theQValue = massTableParent - massTableDaughter - massTableParticle;
      }
 
      const G4double massINCLParent = ParticleTable::getINCLMass(AParent,ZParent,SParent);
      const G4double massINCLDaughter = ParticleTable::getINCLMass(ADaughter,ZDaughter,SDaughter);
      const G4double massINCLParticle = getINCLMass();
 
      // The rhs corresponds to the INCL Q-value
      return theQValue - (massINCLParent-massINCLDaughter-massINCLParticle);
    }

References G4INCL::Particle::getINCLMass(), G4INCL::ParticleTable::getINCLMass(), G4INCL::Particle::getTableMass(), G4INCL::ParticleTable::getTableMass, G4INCL::ParticleTable::getTableQValue(), G4INCL::Particle::isCluster(), G4INCL::Particle::theA, G4INCL::Particle::theS, and G4INCL::Particle::theZ.

◆ getEmissionTime()

G4double G4INCL::Particle::getEmissionTime ( )

inlineinherited

Definition at line 865 of file G4INCLParticle.hh.

865{ return emissionTime; };

G4INCL::Particle::emissionTime

G4double emissionTime

Definition: G4INCLParticle.hh:1104

References G4INCL::Particle::emissionTime.

Referenced by G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::twoBodyDecay().

◆ getEnergy()

G4double G4INCL::Particle::getEnergy ( ) const

inlineinherited

Get the energy of the particle in MeV.

Definition at line 770 of file G4INCLParticle.hh.

    {
      return theEnergy;
    };

References G4INCL::Particle::theEnergy.

Referenced by G4INCL::ProjectileRemnant::addDynamicalSpectator(), G4INCL::Cluster::addParticle(), G4INCL::DeltaDecayChannel::computeDecayTime(), G4INCL::PionResonanceDecayChannel::computeDecayTime(), G4INCL::EtaNElasticChannel::fillFinalState(), G4INCL::EtaNToPiNChannel::fillFinalState(), G4INCL::OmegaNElasticChannel::fillFinalState(), G4INCL::OmegaNToPiNChannel::fillFinalState(), G4INCL::ParticleEntryChannel::fillFinalState(), G4INCL::ReflectionChannel::fillFinalState(), G4INCL::TransmissionChannel::fillFinalState(), G4INCL::ClusteringModelIntercomparison::getCluster(), G4INCL::TransmissionChannel::initializeKineticEnergyOutside(), G4INCL::KinematicsUtils::makeBoostVector(), G4INCL::INCL::makeCompoundNucleus(), G4INCL::KinematicsUtils::momentumInCM(), G4INCL::ParticleEntryChannel::particleEnters(), G4INCL::InteractionAvatar::preInteraction(), G4INCL::InteractionAvatar::preInteractionBlocking(), G4INCL::ProjectileRemnant::removeParticle(), G4INCL::StandardPropagationModel::shootComposite(), G4INCL::StandardPropagationModel::shootParticle(), G4INCL::KinematicsUtils::squareTotalEnergyInCM(), G4INCL::KinematicsUtils::transformToLocalEnergyFrame(), and G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::twoBodyDecay().

◆ getExcitationEnergy()

G4double G4INCL::Nucleus::getExcitationEnergy ( ) const

inline

Get the excitation energy of the nucleus.

Method computeRecoilKinematics() should be called first.

Definition at line 284 of file G4INCLNucleus.hh.

284{ return theExcitationEnergy; }

References G4INCL::Cluster::theExcitationEnergy.

Referenced by fillEventInfo(), getConservationBalance(), G4INCL::INCL::postCascade(), G4INCL::INCL::RecoilCMFunctor::scaleParticleCMMomenta(), and G4INCL::INCL::RecoilFunctor::scaleParticleEnergies().

◆ getFrozenEnergy()

G4double G4INCL::Particle::getFrozenEnergy ( ) const

inlineinherited

Get the frozen particle momentum.

Definition at line 897 of file G4INCLParticle.hh.

897{ return theFrozenEnergy; }

References G4INCL::Particle::theFrozenEnergy.

◆ getFrozenMomentum()

ThreeVector G4INCL::Particle::getFrozenMomentum ( ) const

inlineinherited

Get the frozen particle momentum.

Definition at line 894 of file G4INCLParticle.hh.

894{ return theFrozenMomentum; }

References G4INCL::Particle::theFrozenMomentum.

◆ getHelicity()

G4double G4INCL::Particle::getHelicity ( )

inlineinherited

Definition at line 826 of file G4INCLParticle.hh.

826{ return theHelicity; };

G4INCL::Particle::theHelicity

G4double theHelicity

Definition: G4INCLParticle.hh:1103

References G4INCL::Particle::theHelicity.

Referenced by G4INCL::DeltaDecayChannel::sampleAngles().

◆ getID()

long G4INCL::Particle::getID ( ) const

inlineinherited

Definition at line 981 of file G4INCLParticle.hh.

981{ return ID; };

References G4INCL::Particle::ID.

Referenced by G4INCL::ParticleEntryChannel::fillFinalState(), G4INCL::ReflectionChannel::fillFinalState(), G4INCL::PauliStandard::getBlockingProbability(), G4INCL::SurfaceAvatar::getChannel(), G4INCL::ClusteringModelIntercomparison::getCluster(), G4INCL::ProjectileRemnant::getStoredMomentum(), G4INCL::ParticleEntryChannel::particleEnters(), and G4INCL::ProjectileRemnant::reset().

◆ getINCLMass()

G4double G4INCL::Particle::getINCLMass ( ) const

inlineinherited

Get the INCL particle mass.

Definition at line 454 of file G4INCLParticle.hh.

                                        {
      switch(theType) {
        case Proton:
        case Neutron:
        case PiPlus:
        case PiMinus:
        case PiZero:
        case Lambda:
        case SigmaPlus:
        case SigmaZero:
        case SigmaMinus:
        case KPlus:
        case KZero:
        case KZeroBar:
        case KShort:
        case KLong:
        case KMinus:
        case Eta:
        case Omega:
        case EtaPrime:
        case Photon:
          return ParticleTable::getINCLMass(theType);
          break;
 
        case DeltaPlusPlus:
        case DeltaPlus:
        case DeltaZero:
        case DeltaMinus:
          return theMass;
          break;
 
        case Composite:
          return ParticleTable::getINCLMass(theA,theZ,theS);
          break;
 
        default:
          INCL_ERROR("Particle::getINCLMass: Unknown particle type." << '\n');
          return 0.0;
          break;
      }
    }

References G4INCL::Composite, G4INCL::DeltaMinus, G4INCL::DeltaPlus, G4INCL::DeltaPlusPlus, G4INCL::DeltaZero, G4INCL::Eta, G4INCL::EtaPrime, G4INCL::ParticleTable::getINCLMass(), INCL_ERROR, G4INCL::KLong, G4INCL::KMinus, G4INCL::KPlus, G4INCL::KShort, G4INCL::KZero, G4INCL::KZeroBar, G4INCL::Lambda, G4INCL::Neutron, G4INCL::Omega, G4INCL::Photon, G4INCL::PiMinus, G4INCL::PiPlus, G4INCL::PiZero, G4INCL::Proton, G4INCL::SigmaMinus, G4INCL::SigmaPlus, G4INCL::SigmaZero, G4INCL::Particle::theA, G4INCL::Particle::theMass, G4INCL::Particle::theS, G4INCL::Particle::theType, and G4INCL::Particle::theZ.

Referenced by G4INCL::ParticleEntryChannel::fillFinalState(), G4INCL::Particle::getEmissionQValueCorrection(), and G4INCL::Particle::setINCLMass().

◆ getIncomingAngularMomentum()

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

inline

Get the incoming angular-momentum vector.

Definition at line 262 of file G4INCLNucleus.hh.

262{ return incomingAngularMomentum; }

References incomingAngularMomentum.

Referenced by G4INCL::INCL::makeCompoundNucleus().

◆ getIncomingMomentum()

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

inline

Get the incoming momentum vector.

Definition at line 270 of file G4INCLNucleus.hh.

                                                   {
      return incomingMomentum;
    }

References incomingMomentum.

Referenced by getConservationBalance(), G4INCL::INCL::makeCompoundNucleus(), G4INCL::INCL::RecoilCMFunctor::RecoilCMFunctor(), and G4INCL::INCL::RecoilFunctor::scaleParticleEnergies().

◆ getInitialA()

G4int G4INCL::Nucleus::getInitialA ( ) const

inline

Definition at line 108 of file G4INCLNucleus.hh.

108{ return theInitialA; };

References theInitialA.

◆ getInitialEnergy()

G4double G4INCL::Nucleus::getInitialEnergy ( ) const

inline

Get the initial energy.

Definition at line 278 of file G4INCLNucleus.hh.

278{ return initialEnergy; }

References initialEnergy.

Referenced by getConservationBalance(), and G4INCL::INCL::RecoilCMFunctor::RecoilCMFunctor().

◆ getInitialInternalEnergy()

G4double G4INCL::Nucleus::getInitialInternalEnergy ( ) const

inline

Definition at line 353 of file G4INCLNucleus.hh.

353{ return initialInternalEnergy; };

References initialInternalEnergy.

Referenced by G4INCL::CDPP::isBlocked().

◆ getInitialS()

G4int G4INCL::Nucleus::getInitialS ( ) const

inline

Definition at line 110 of file G4INCLNucleus.hh.

110{ return theInitialS; };

References theInitialS.

◆ getInitialZ()

G4int G4INCL::Nucleus::getInitialZ ( ) const

inline

Definition at line 109 of file G4INCLNucleus.hh.

109{ return theInitialZ; };

References theInitialZ.

◆ getInvariantMass()

G4double G4INCL::Particle::getInvariantMass ( ) const

inlineinherited

Get the the particle invariant mass.

Uses the relativistic invariant

$m = \sqrt{E^2 - {\vec p}^2}$

Definition at line 748 of file G4INCLParticle.hh.

                                      {
      const G4double mass = std::pow(theEnergy, 2) - theMomentum.dot(theMomentum);
      if(mass < 0.0) {
        INCL_ERROR("E*E - p*p is negative." << '\n');
        return 0.0;
      } else {
        return std::sqrt(mass);
      }
    };

References G4INCL::ThreeVector::dot(), INCL_ERROR, G4INCL::Particle::theEnergy, and G4INCL::Particle::theMomentum.

Referenced by finalizeProjectileRemnant(), and G4INCL::Particle::Particle().

◆ getKineticEnergy()

G4double G4INCL::Particle::getKineticEnergy ( ) const

inlineinherited

Get the particle kinetic energy.

Definition at line 759 of file G4INCLParticle.hh.

759{ return theEnergy - theMass; }

References G4INCL::Particle::theEnergy, and G4INCL::Particle::theMass.

Referenced by G4INCL::NuclearPotential::NuclearPotentialEnergyIsospin::computePotentialEnergy(), G4INCL::NuclearPotential::NuclearPotentialEnergyIsospinSmooth::computePotentialEnergy(), emitInsideKaon(), emitInsideLambda(), emitInsidePions(), emitInsideStrangeParticles(), fillEventInfo(), G4INCL::ParticleEntryChannel::fillFinalState(), G4INCL::SurfaceAvatar::getChannel(), getConservationBalance(), G4INCL::KinematicsUtils::getLocalEnergy(), G4INCL::SurfaceAvatar::getTransmissionProbability(), G4INCL::CoulombDistortion::maxImpactParameter(), G4INCL::CoulombNonRelativistic::minimumDistance(), G4INCL::ParticleEntryChannel::particleEnters(), G4INCL::CDPP::processOneParticle(), G4INCL::INCL::RecoilFunctor::RecoilFunctor(), and G4INCL::StandardPropagationModel::shootParticle().

◆ getLongitudinalPosition()

ThreeVector G4INCL::Particle::getLongitudinalPosition ( ) const

inlineinherited

Longitudinal component of the position w.r.t. the momentum.

Definition at line 873 of file G4INCLParticle.hh.

                                                {
      return *thePropagationMomentum * (thePosition.dot(*thePropagationMomentum)/thePropagationMomentum->mag2());
    }

References G4INCL::ThreeVector::dot(), G4INCL::ThreeVector::mag2(), G4INCL::Particle::thePosition, and G4INCL::Particle::thePropagationMomentum.

Referenced by G4INCL::CoulombNonRelativistic::coulombDeviation(), and G4INCL::Particle::getTransversePosition().

◆ getMass()

G4double G4INCL::Particle::getMass ( ) const

inlineinherited

Get the cached particle mass.

Definition at line 451 of file G4INCLParticle.hh.

451{ return theMass; }

References G4INCL::Particle::theMass.

◆ getMomentum()

const G4INCL::ThreeVector & G4INCL::Particle::getMomentum ( ) const

inlineinherited

Get the momentum vector.

Definition at line 794 of file G4INCLParticle.hh.

    {
      return theMomentum;
    };

References G4INCL::Particle::theMomentum.

◆ getNumberOfCollisions()

G4int G4INCL::Particle::getNumberOfCollisions ( ) const

inlineinherited

Return the number of collisions undergone by the particle.

Definition at line 834 of file G4INCLParticle.hh.

834{ return nCollisions; }

References G4INCL::Particle::nCollisions.

Referenced by G4INCL::Cluster::addParticle().

◆ getNumberOfDecays()

G4int G4INCL::Particle::getNumberOfDecays ( ) const

inlineinherited

Return the number of decays undergone by the particle.

Definition at line 843 of file G4INCLParticle.hh.

843{ return nDecays; }

G4INCL::Particle::nDecays

G4int nDecays

Definition: G4INCLParticle.hh:1091

References G4INCL::Particle::nDecays.

◆ getNumberOfEnteringKaons()

G4int G4INCL::Nucleus::getNumberOfEnteringKaons ( ) const

inline

Definition at line 122 of file G4INCLNucleus.hh.

122{ return theNkaonplusInitial+theNkaonminusInitial; };

References theNkaonminusInitial, and theNkaonplusInitial.

◆ getNumberOfEnteringNeutrons()

G4int G4INCL::Nucleus::getNumberOfEnteringNeutrons ( ) const

inline

Definition at line 120 of file G4INCLNucleus.hh.

120{ return theNnInitial; };

References theNnInitial.

◆ getNumberOfEnteringPions()

G4int G4INCL::Nucleus::getNumberOfEnteringPions ( ) const

inline

Definition at line 121 of file G4INCLNucleus.hh.

121{ return theNpionplusInitial+theNpionminusInitial; };

References theNpionminusInitial, and theNpionplusInitial.

◆ getNumberOfEnteringProtons()

G4int G4INCL::Nucleus::getNumberOfEnteringProtons ( ) const

inline

Definition at line 119 of file G4INCLNucleus.hh.

119{ return theNpInitial; };

References theNpInitial.

◆ getNumberOfKaon()

G4int G4INCL::Particle::getNumberOfKaon ( ) const

inlineinherited

Number of Kaon inside de nucleus.

Put in the Particle class in order to calculate the "correct" mass of composit particle.

Definition at line 1053 of file G4INCLParticle.hh.

1053{ return theNKaon; };

References G4INCL::Particle::theNKaon.

Referenced by G4INCL::ParticleEntryChannel::fillFinalState(), G4INCL::SurfaceAvatar::getChannel(), and G4INCL::BinaryCollisionAvatar::postInteraction().

◆ getParticipantType()

ParticipantType G4INCL::Particle::getParticipantType ( ) const

inlineinherited

Definition at line 310 of file G4INCLParticle.hh.

                                               {
      return theParticipantType;
    }

References G4INCL::Particle::theParticipantType.

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

◆ getParticleBias()

G4double G4INCL::Particle::getParticleBias ( ) const

inlineinherited

Get the particle bias.

Definition at line 1032 of file G4INCLParticle.hh.

1032{ return theParticleBias; };

G4INCL::Particle::theParticleBias

G4double theParticleBias

Definition: G4INCLParticle.hh:1098

References G4INCL::Particle::theParticleBias.

◆ getParticleList()

ParticleList G4INCL::Cluster::getParticleList ( ) const

inlineinherited

Returns the list of particles that make up the cluster.

Definition at line 219 of file G4INCLCluster.hh.

219{ return particles; }

References G4INCL::Cluster::particles.

◆ getParticles()

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

inlineinherited

Get the list of particles in the cluster.

Definition at line 169 of file G4INCLCluster.hh.

169{ return particles; }

References G4INCL::Cluster::particles.

Referenced by applyFinalState(), G4INCL::CoulombNone::bringToSurface(), G4INCL::INCL::makeCompoundNucleus(), G4INCL::INCL::makeProjectileRemnant(), G4INCL::INCL::postCascade(), and G4INCL::SurfaceAvatar::postInteraction().

◆ getPosition()

const G4INCL::ThreeVector & G4INCL::Particle::getPosition ( ) const

inlineinherited

Set the position vector.

Definition at line 816 of file G4INCLParticle.hh.

    {
      return thePosition;
    };

References G4INCL::Particle::thePosition.

◆ getPotential()

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

inline

Getter for thePotential.

Definition at line 473 of file G4INCLNucleus.hh.

473{ return thePotential; };

References thePotential.

Referenced by G4INCL::ParticleEntryChannel::fillFinalState(), G4INCL::PauliStandard::getBlockingProbability(), G4INCL::SurfaceAvatar::getChannel(), G4INCL::CDPP::isBlocked(), and G4INCL::ParticleEntryChannel::particleEnters().

◆ getPotentialEnergy()

G4double G4INCL::Particle::getPotentialEnergy ( ) const

inlineinherited

Get the particle potential energy.

Definition at line 762 of file G4INCLParticle.hh.

762{ return thePotentialEnergy; }

References G4INCL::Particle::thePotentialEnergy.

Referenced by G4INCL::Cluster::addParticle(), emitInsideKaon(), emitInsideLambda(), emitInsidePions(), emitInsideStrangeParticles(), G4INCL::ReflectionChannel::fillFinalState(), G4INCL::TransmissionChannel::fillFinalState(), G4INCL::ClusteringModelIntercomparison::getCluster(), G4INCL::KinematicsUtils::getLocalEnergy(), G4INCL::SurfaceAvatar::getTransmissionProbability(), G4INCL::TransmissionChannel::initializeKineticEnergyOutside(), G4INCL::InteractionAvatar::preInteractionBlocking(), and G4INCL::CDPP::processOneParticle().

◆ getProjectileRemnant()

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

inline

Get the projectile remnant.

Definition at line 439 of file G4INCLNucleus.hh.

439{ return theProjectileRemnant; }

References theProjectileRemnant.

Referenced by G4INCL::ParticleEntryChannel::fillFinalState(), G4INCL::INCL::makeCompoundNucleus(), G4INCL::INCL::makeProjectileRemnant(), and G4INCL::INCL::postCascade().

◆ getPropagationVelocity()

ThreeVector G4INCL::Particle::getPropagationVelocity ( ) const

inlineinherited

Get the propagation velocity of the particle.

Definition at line 900 of file G4INCLParticle.hh.

900{ return (*thePropagationMomentum)/(*thePropagationEnergy); }

References G4INCL::Particle::thePropagationMomentum.

Referenced by G4INCL::CoulombNone::bringToSurface(), G4INCL::StandardPropagationModel::getReflectionTime(), and G4INCL::StandardPropagationModel::getTime().

◆ getRealMass()

G4double G4INCL::Particle::getRealMass ( ) const

inlineinherited

Get the real particle mass.

Definition at line 540 of file G4INCLParticle.hh.

                                        {
      switch(theType) {
        case Proton:
        case Neutron:
        case PiPlus:
        case PiMinus:
        case PiZero:
        case Lambda:
        case SigmaPlus:
        case SigmaZero:
        case SigmaMinus:
        case KPlus:
        case KZero:
        case KZeroBar:
        case KShort:
        case KLong:
        case KMinus:
        case Eta:
        case Omega:
        case EtaPrime:
        case Photon:
          return ParticleTable::getRealMass(theType);
          break;
 
        case DeltaPlusPlus:
        case DeltaPlus:
        case DeltaZero:
        case DeltaMinus:
          return theMass;
          break;
 
        case Composite:
          return ParticleTable::getRealMass(theA,theZ,theS);
          break;
 
        default:
          INCL_ERROR("Particle::getRealMass: Unknown particle type." << '\n');
          return 0.0;
          break;
      }
    }

References G4INCL::Composite, G4INCL::DeltaMinus, G4INCL::DeltaPlus, G4INCL::DeltaPlusPlus, G4INCL::DeltaZero, G4INCL::Eta, G4INCL::EtaPrime, G4INCL::ParticleTable::getRealMass(), INCL_ERROR, G4INCL::KLong, G4INCL::KMinus, G4INCL::KPlus, G4INCL::KShort, G4INCL::KZero, G4INCL::KZeroBar, G4INCL::Lambda, G4INCL::Neutron, G4INCL::Omega, G4INCL::Photon, G4INCL::PiMinus, G4INCL::PiPlus, G4INCL::PiZero, G4INCL::Proton, G4INCL::SigmaMinus, G4INCL::SigmaPlus, G4INCL::SigmaZero, G4INCL::Particle::theA, G4INCL::Particle::theMass, G4INCL::Particle::theS, G4INCL::Particle::theType, and G4INCL::Particle::theZ.

Referenced by G4INCL::Cluster::getTableMass(), and G4INCL::Particle::setRealMass().

◆ getReflectionMomentum()

G4double G4INCL::Particle::getReflectionMomentum ( ) const

inlineinherited

Return the reflection momentum.

The reflection momentum is used by calls to getSurfaceRadius to compute the radius of the sphere where the nucleon moves. It is necessary to introduce fuzzy r-p correlations.

Definition at line 997 of file G4INCLParticle.hh.

                                           {
      if(rpCorrelated)
        return theMomentum.mag();
      else
        return uncorrelatedMomentum;
    }

References G4INCL::ThreeVector::mag(), G4INCL::Particle::rpCorrelated, G4INCL::Particle::theMomentum, and G4INCL::Particle::uncorrelatedMomentum.

Referenced by G4INCL::KinematicsUtils::getLocalEnergy(), and getSurfaceRadius().

◆ getS()

G4int G4INCL::Particle::getS ( ) const

inlineinherited

Returns the strangeness number.

Definition at line 399 of file G4INCLParticle.hh.

399{ return theS; }

References G4INCL::Particle::theS.

Referenced by G4INCL::ProjectileRemnant::addDynamicalSpectator(), G4INCL::Cluster::addParticle(), G4INCL::ClusteringModelIntercomparison::clusterCanEscape(), computeOneNucleonRecoilKinematics(), G4INCL::ClusterDecay::decay(), emitInsideKaon(), emitInsideLambda(), emitInsideStrangeParticles(), fillEventInfo(), G4INCL::ParticleEntryChannel::fillFinalState(), G4INCL::TransmissionChannel::fillFinalState(), finalizeProjectileRemnant(), getConservationBalance(), G4INCL::SurfaceAvatar::getTransmissionProbability(), G4INCL::TransmissionChannel::initializeKineticEnergyOutside(), insertParticle(), G4INCL::ClusterDecay::isStable(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::phaseSpaceDecay(), G4INCL::INCL::preCascade(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::recursiveDecay(), G4INCL::ProjectileRemnant::removeParticle(), G4INCL::INCL::RecoilCMFunctor::scaleParticleCMMomenta(), G4INCL::INCL::RecoilFunctor::scaleParticleEnergies(), G4INCL::StandardPropagationModel::shootComposite(), G4INCL::StandardPropagationModel::shootParticle(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::threeBodyDecay(), and G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::twoBodyDecay().

◆ getSpecies()

ParticleSpecies G4INCL::Cluster::getSpecies ( ) const

inlinevirtualinherited

Get the particle species.

Reimplemented from G4INCL::Particle.

Definition at line 132 of file G4INCLCluster.hh.

                                       {
      return ParticleSpecies(theA, theZ, theS);
    }

References G4INCL::Particle::theA, G4INCL::Particle::theS, and G4INCL::Particle::theZ.

◆ getSpin()

ThreeVector const & G4INCL::Cluster::getSpin ( ) const

inlineinherited

Get the spin of the nucleus.

Definition at line 422 of file G4INCLCluster.hh.

422{ return theSpin; }

References G4INCL::Cluster::theSpin.

Referenced by fillEventInfo(), and G4INCL::Cluster::getAngularMomentum().

◆ getStore()

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

inline

Definition at line 347 of file G4INCLNucleus.hh.

347{return theStore; };

References theStore.

Referenced by G4INCL::AvatarDumpAction::afterAvatarUserAction(), G4INCL::INCL::continueCascade(), G4INCL::InteractionAvatar::enforceEnergyConservation(), fillEventInfo(), G4INCL::StandardPropagationModel::generateAllAvatars(), G4INCL::StandardPropagationModel::generateBinaryCollisionAvatar(), G4INCL::PauliStandard::getBlockingProbability(), G4INCL::BinaryCollisionAvatar::getChannel(), G4INCL::SurfaceAvatar::getChannel(), G4INCL::ClusteringModelIntercomparison::getCluster(), G4INCL::SurfaceAvatar::getTransmissionProbability(), G4INCL::INCL::initializeTarget(), G4INCL::CDPP::isBlocked(), G4INCL::PauliStrictStandard::isBlocked(), G4INCL::INCL::makeCompoundNucleus(), G4INCL::INCL::makeProjectileRemnant(), G4INCL::INCL::postCascade(), G4INCL::BinaryCollisionAvatar::postInteraction(), G4INCL::InteractionAvatar::postInteraction(), G4INCL::SurfaceAvatar::postInteraction(), G4INCL::DecayAvatar::postInteraction(), G4INCL::StandardPropagationModel::propagate(), G4INCL::StandardPropagationModel::registerAvatar(), G4INCL::StandardPropagationModel::shootComposite(), G4INCL::StandardPropagationModel::shootParticle(), G4INCL::InteractionAvatar::shouldUseLocalEnergy(), and G4INCL::StandardPropagationModel::updateAvatars().

◆ 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 402 of file G4INCLNucleus.hh.

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

References G4INCL::NuclearPotential::INuclearPotential::getFermiMomentum(), G4INCL::NuclearDensity::getMaxRFromP(), G4INCL::Particle::getReflectionMomentum(), G4INCL::Particle::getType(), getUniverseRadius(), G4INCL::Particle::isLambda(), G4INCL::Particle::isNucleon(), G4INCL::Particle::isResonance(), theDensity, and thePotential.

Referenced by G4INCL::InteractionAvatar::bringParticleInside(), G4INCL::StandardPropagationModel::generateBinaryCollisionAvatar(), G4INCL::StandardPropagationModel::getReflectionTime(), and G4INCL::InteractionAvatar::postInteraction().

◆ getTableMass()

virtual G4double G4INCL::Cluster::getTableMass ( ) const

inlinevirtualinherited

Get the real particle mass.

Overloads the Particle method.

Reimplemented from G4INCL::Particle.

Definition at line 164 of file G4INCLCluster.hh.

164{ return getRealMass(); }

G4INCL::Particle::getRealMass

G4double getRealMass() const

Get the real particle mass.

Definition: G4INCLParticle.hh:540

References G4INCL::Particle::getRealMass().

Referenced by G4INCL::Cluster::computeDynamicalPotential(), and useFusionKinematics().

◆ getTotalBias()

G4double G4INCL::Particle::getTotalBias ( )

staticinherited

General bias vector function.

Definition at line 300 of file G4INCLParticle.cc.

                                  {
      G4double TotalBias = 1.;
      for(G4int i=0; i<G4int(INCLBiasVector.size());i++) TotalBias *= Particle::INCLBiasVector[i];
      return TotalBias;
  }

References G4INCL::Particle::INCLBiasVector.

Referenced by emitInsideKaon(), emitInsideLambda(), emitInsidePions(), emitInsideStrangeParticles(), G4INCL::INCL::postCascade(), and G4INCL::EventInfo::remnantToParticle().

◆ getTransferQValueCorrection() [1/2]

G4double G4INCL::Particle::getTransferQValueCorrection	(	const G4int	AFrom,
		const G4int	ZFrom,
		const G4int	ATo,
		const G4int	ZTo
	)		const

inlineinherited

Computes correction on the transfer Q-value.

Computes the correction that must be applied to INCL particles in order to obtain the correct Q-value for particle transfer from a given nucleus to another.

Assumes that the receving nucleus is INCL's target nucleus, with the INCL separation energy.

Parameters

AFrom	the mass number of the donating nucleus
ZFrom	the charge number of the donating nucleus
ATo	the mass number of the receiving nucleus
ZTo	the charge number of the receiving nucleus

Returns: the correction

Definition at line 642 of file G4INCLParticle.hh.

                                                                                                                       {
      const G4int SFrom = 0;
      const G4int STo = 0;
      const G4int AFromDaughter = AFrom - theA;
      const G4int ZFromDaughter = ZFrom - theZ;
      const G4int SFromDaughter = 0;
      const G4int AToDaughter = ATo + theA;
      const G4int ZToDaughter = ZTo + theZ;
      const G4int SToDaughter = 0;
      const G4double theQValue = ParticleTable::getTableQValue(AToDaughter,ZToDaughter,SToDaughter,AFromDaughter,ZFromDaughter,SFromDaughter,AFrom,ZFrom,SFrom);
 
      const G4double massINCLTo = ParticleTable::getINCLMass(ATo,ZTo,STo);
      const G4double massINCLToDaughter = ParticleTable::getINCLMass(AToDaughter,ZToDaughter,SToDaughter);
      /* Note that here we have to use the table mass in the INCL Q-value. We
       * cannot use theMass, because at this stage the particle is probably
       * still off-shell; and we cannot use getINCLMass(), because it leads to
       * violations of global energy conservation.
       */
      const G4double massINCLParticle = getTableMass();
 
      // The rhs corresponds to the INCL Q-value for particle absorption
      return theQValue - (massINCLToDaughter-massINCLTo-massINCLParticle);
    }

References G4INCL::ParticleTable::getINCLMass(), G4INCL::Particle::getTableMass(), G4INCL::ParticleTable::getTableQValue(), G4INCL::Particle::theA, and G4INCL::Particle::theZ.

◆ getTransferQValueCorrection() [2/2]

G4double G4INCL::Particle::getTransferQValueCorrection	(	const G4int	AFrom,
		const G4int	ZFrom,
		const G4int	SFrom,
		const G4int	ATo,
		const G4int	ZTo,
		const G4int	STo
	)		const

inlineinherited

Computes correction on the transfer Q-value for hypernuclei.

Computes the correction that must be applied to INCL particles in order to obtain the correct Q-value for particle transfer from a given nucleus to another.

Assumes that the receving nucleus is INCL's target nucleus, with the INCL separation energy.

Parameters

AFrom	the mass number of the donating nucleus
ZFrom	the charge number of the donating nucleus
SFrom	the strangess number of the donating nucleus
ATo	the mass number of the receiving nucleus
ZTo	the charge number of the receiving nucleus
STo	the strangess number of the receiving nucleus

Returns: the correction

Definition at line 719 of file G4INCLParticle.hh.

                                                                                                                                                            {
      const G4int AFromDaughter = AFrom - theA;
      const G4int ZFromDaughter = ZFrom - theZ;
      const G4int SFromDaughter = SFrom - theS;
      const G4int AToDaughter = ATo + theA;
      const G4int ZToDaughter = ZTo + theZ;
      const G4int SToDaughter = STo + theS;
      const G4double theQValue = ParticleTable::getTableQValue(AToDaughter,ZToDaughter,SFromDaughter,AFromDaughter,ZFromDaughter,SToDaughter,AFrom,ZFrom,SFrom);
 
      const G4double massINCLTo = ParticleTable::getINCLMass(ATo,ZTo,STo);
      const G4double massINCLToDaughter = ParticleTable::getINCLMass(AToDaughter,ZToDaughter,SToDaughter);
      /* Note that here we have to use the table mass in the INCL Q-value. We
       * cannot use theMass, because at this stage the particle is probably
       * still off-shell; and we cannot use getINCLMass(), because it leads to
       * violations of global energy conservation.
       */
      const G4double massINCLParticle = getTableMass();
 
      // The rhs corresponds to the INCL Q-value for particle absorption
      return theQValue - (massINCLToDaughter-massINCLTo-massINCLParticle);
    }

References G4INCL::ParticleTable::getINCLMass(), G4INCL::Particle::getTableMass(), G4INCL::ParticleTable::getTableQValue(), G4INCL::Particle::theA, G4INCL::Particle::theS, and G4INCL::Particle::theZ.

◆ getTransmissionBarrier()

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

inline

Get the transmission barrier.

Definition at line 375 of file G4INCLNucleus.hh.

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

References G4INCL::PhysicalConstants::eSquared, G4INCL::NuclearDensity::getTransmissionRadius(), G4INCL::Particle::getZ(), theDensity, and G4INCL::Particle::theZ.

Referenced by G4INCL::SurfaceAvatar::getTransmissionProbability(), and G4INCL::InteractionAvatar::postInteraction().

◆ getTransversePosition()

ThreeVector G4INCL::Particle::getTransversePosition ( ) const

inlineinherited

Transverse component of the position w.r.t. the momentum.

Definition at line 868 of file G4INCLParticle.hh.

                                              {
      return thePosition - getLongitudinalPosition();
    }

References G4INCL::Particle::getLongitudinalPosition(), and G4INCL::Particle::thePosition.

Referenced by G4INCL::CoulombNonRelativistic::coulombDeviation(), G4INCL::StandardPropagationModel::shootComposite(), and G4INCL::StandardPropagationModel::shootParticle().

◆ getTryCompoundNucleus()

G4bool G4INCL::Nucleus::getTryCompoundNucleus ( )

inline

Definition at line 372 of file G4INCLNucleus.hh.

372{ return tryCN; }

References tryCN.

Referenced by G4INCL::INCL::continueCascade(), G4INCL::INCL::postCascade(), and G4INCL::INCL::updateGlobalInfo().

◆ getType()

G4INCL::ParticleType G4INCL::Particle::getType ( ) const

inlineinherited

Get the particle type.

See also: G4INCL::ParticleType

Definition at line 178 of file G4INCLParticle.hh.

                                       {
      return theType;
    };

References G4INCL::Particle::theType.

◆ getUniverseRadius()

G4double G4INCL::Nucleus::getUniverseRadius ( ) const

inline

Getter for theUniverseRadius.

Definition at line 418 of file G4INCLNucleus.hh.

418{ return theUniverseRadius; }

References theUniverseRadius.

Referenced by G4INCL::PauliStandard::getBlockingProbability(), G4INCL::ClusteringModelIntercomparison::getCluster(), getSurfaceRadius(), G4INCL::StandardPropagationModel::shootComposite(), and G4INCL::StandardPropagationModel::shootParticle().

◆ getZ()

G4int G4INCL::Particle::getZ ( ) const

inlineinherited

Returns the charge number.

Definition at line 396 of file G4INCLParticle.hh.

396{ return theZ; }

References G4INCL::Particle::theZ.

Referenced by G4INCL::ProjectileRemnant::addDynamicalSpectator(), G4INCL::Cluster::addParticle(), G4INCL::CoulombNonRelativistic::bringToSurface(), computeOneNucleonRecoilKinematics(), G4INCL::ClusterDecay::decay(), G4INCL::CoulombNonRelativistic::distortOut(), emitInsideKaon(), emitInsidePions(), emitInsideStrangeParticles(), fillEventInfo(), G4INCL::ParticleEntryChannel::fillFinalState(), G4INCL::TransmissionChannel::fillFinalState(), finalizeProjectileRemnant(), G4INCL::ClusteringModelIntercomparison::findClusterStartingFrom(), G4INCL::PauliStandard::getBlockingProbability(), G4INCL::ClusteringModelIntercomparison::getCluster(), getConservationBalance(), getTransmissionBarrier(), G4INCL::SurfaceAvatar::getTransmissionProbability(), G4INCL::NuclearDensity::getTransmissionRadius(), G4INCL::TransmissionChannel::initializeKineticEnergyOutside(), insertParticle(), G4INCL::ClusterDecay::isStable(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::phaseSpaceDecay(), G4INCL::INCL::preCascade(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::recursiveDecay(), G4INCL::ProjectileRemnant::removeParticle(), G4INCL::INCL::RecoilCMFunctor::scaleParticleCMMomenta(), G4INCL::INCL::RecoilFunctor::scaleParticleEnergies(), G4INCL::StandardPropagationModel::shootComposite(), G4INCL::StandardPropagationModel::shootParticle(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::threeBodyDecay(), and G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::twoBodyDecay().

◆ hasRemnant()

G4bool G4INCL::Nucleus::hasRemnant ( ) const

inline

Does the nucleus give a cascade remnant?

To be called after computeRecoilKinematics().

Definition at line 365 of file G4INCLNucleus.hh.

365{ return remnant; }

References remnant.

Referenced by fillEventInfo(), getConservationBalance(), and G4INCL::INCL::postCascade().

◆ incrementNumberOfCollisions()

void G4INCL::Particle::incrementNumberOfCollisions ( )

inlineinherited

Increment the number of collisions undergone by the particle.

Definition at line 840 of file G4INCLParticle.hh.

840{ nCollisions++; }

References G4INCL::Particle::nCollisions.

◆ incrementNumberOfDecays()

void G4INCL::Particle::incrementNumberOfDecays ( )

inlineinherited

Increment the number of decays undergone by the particle.

Definition at line 849 of file G4INCLParticle.hh.

849{ nDecays++; }

References G4INCL::Particle::nDecays.

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

References G4INCL::Store::add(), computeTotalEnergy(), initialCenterOfMass, initialInternalEnergy, G4INCL::Cluster::initializeParticles(), G4INCL::Cluster::particles, G4INCL::Particle::thePosition, theProjectileRemnant, theStore, and updatePotentialEnergy().

Referenced by G4INCL::INCL::initializeTarget().

◆ 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();
      theS += p->getS();
      theStore->particleHasEntered(p);
      if(p->isNucleon()) {
        theNpInitial += Math::heaviside(ParticleTable::getIsospin(p->getType()));
        theNnInitial += Math::heaviside(-ParticleTable::getIsospin(p->getType()));
      }
      if(p->isPion()) {
        theNpionplusInitial += Math::heaviside(ParticleTable::getIsospin(p->getType()));
        theNpionminusInitial += Math::heaviside(-ParticleTable::getIsospin(p->getType()));
      }
      if(p->isKaon() || p->isAntiKaon()) {
        theNkaonplusInitial += Math::heaviside(ParticleTable::getIsospin(p->getType()));
        theNkaonminusInitial += Math::heaviside(-ParticleTable::getIsospin(p->getType()));
      }
      if(!p->isTargetSpectator()) theStore->getBook().incrementCascading();
    };

Referenced by applyFinalState().

◆ internalBoostToCM()

void G4INCL::Cluster::internalBoostToCM ( )

inlineinherited

Boost to the CM of the component particles.

The position of all particles in the particles list is shifted so that their centre of mass is in the origin and their total momentum is zero.

Definition at line 254 of file G4INCLCluster.hh.

                             {
 
      // First compute the current CM position and total momentum
      ThreeVector theCMPosition, theTotalMomentum;
//    G4double theTotalEnergy = 0.0;
      for(ParticleIter p=particles.begin(), e=particles.end(); p!=e; ++p) {
        theCMPosition += (*p)->getPosition();
        theTotalMomentum += (*p)->getMomentum();
//      theTotalEnergy += (*p)->getEnergy();
      }
      theCMPosition /= theA;
// assert((unsigned int)theA==particles.size());
 
      // Now determine the CM velocity of the particles
      // commented out because currently unused, see below
      // ThreeVector betaCM = theTotalMomentum / theTotalEnergy;
 
      // The new particle positions and momenta are scaled by a factor of
      // \f$\sqrt{A/(A-1)}\f$, so that the resulting density distributions in
      // the CM have the same variance as the one we started with.
      const G4double rescaling = std::sqrt(((G4double)theA)/((G4double)(theA-1)));
 
      // Loop again to boost and reposition
      for(ParticleIter p=particles.begin(), e=particles.end(); p!=e; ++p) {
        // \bug{We should do the following, but the Fortran version actually
        // does not!
        // (*p)->boost(betaCM);
        // Here is what the Fortran version does:}
        (*p)->setMomentum(((*p)->getMomentum()-theTotalMomentum/theA)*rescaling);
 
        // Set the CM position of the particles
        (*p)->setPosition(((*p)->getPosition()-theCMPosition)*rescaling);
      }
 
      // Set the global cluster kinematic variables
      thePosition.setX(0.0);
      thePosition.setY(0.0);
      thePosition.setZ(0.0);
      theMomentum.setX(0.0);
      theMomentum.setY(0.0);
      theMomentum.setZ(0.0);
      theEnergy = getMass();
 
      INCL_DEBUG("Cluster boosted to internal CM:" << '\n' << print());
 
    }

References G4INCL::Particle::getMass(), INCL_DEBUG, G4INCL::Cluster::particles, G4INCL::Cluster::print(), G4INCL::ThreeVector::setX(), G4INCL::ThreeVector::setY(), G4INCL::ThreeVector::setZ(), G4INCL::Particle::theA, G4INCL::Particle::theEnergy, G4INCL::Particle::theMomentum, and G4INCL::Particle::thePosition.

Referenced by G4INCL::ProjectileRemnant::ProjectileRemnant().

◆ isAntiKaon()

G4bool G4INCL::Particle::isAntiKaon ( ) const

inlineinherited

Is this an antiKaon?

Definition at line 372 of file G4INCLParticle.hh.

372{ return (theType == KZeroBar || theType == KMinus); }

References G4INCL::KMinus, G4INCL::KZeroBar, and G4INCL::Particle::theType.

◆ isBaryon()

G4bool G4INCL::Particle::isBaryon ( ) const

inlineinherited

Is this a Baryon?

Definition at line 387 of file G4INCLParticle.hh.

387{ return (isNucleon() || isResonance() || isHyperon()); }

G4INCL::Particle::isResonance

G4bool isResonance() const

Is it a resonance?

Definition: G4INCLParticle.hh:358

G4INCL::Particle::isHyperon

G4bool isHyperon() const

Is this an Hyperon?

Definition: G4INCLParticle.hh:381

G4INCL::Particle::isNucleon

G4bool isNucleon() const

Definition: G4INCLParticle.hh:303

References G4INCL::Particle::isHyperon(), G4INCL::Particle::isNucleon(), and G4INCL::Particle::isResonance().

Referenced by G4INCL::CDPP::processOneParticle().

◆ isCluster()

G4bool G4INCL::Particle::isCluster ( ) const

inlineinherited

Definition at line 883 of file G4INCLParticle.hh.

                             {
      return (theType == Composite);
    }

References G4INCL::Composite, and G4INCL::Particle::theType.

Referenced by G4INCL::Particle::getEmissionQValueCorrection(), and G4INCL::SurfaceAvatar::postInteraction().

◆ isDelta()

G4bool G4INCL::Particle::isDelta ( ) const

inlineinherited

Is it a Delta?

Definition at line 361 of file G4INCLParticle.hh.

                                  {
      return (theType==DeltaPlusPlus || theType==DeltaPlus ||
          theType==DeltaZero || theType==DeltaMinus); }

References G4INCL::DeltaMinus, G4INCL::DeltaPlus, G4INCL::DeltaPlusPlus, G4INCL::DeltaZero, and G4INCL::Particle::theType.

◆ isEta()

G4bool G4INCL::Particle::isEta ( ) const

inlineinherited

Is this an eta?

Definition at line 346 of file G4INCLParticle.hh.

346{ return (theType == Eta); }

References G4INCL::Eta, and G4INCL::Particle::theType.

Referenced by G4INCL::CrossSectionsMultiPionsAndResonances::elastic(), G4INCL::CrossSectionsStrangeness::elastic(), G4INCL::CrossSectionsMultiPionsAndResonances::etaNElastic(), G4INCL::CrossSectionsMultiPionsAndResonances::etaNToPiN(), G4INCL::CrossSectionsMultiPionsAndResonances::etaNToPiPiN(), G4INCL::BinaryCollisionAvatar::getChannel(), G4INCL::DecayAvatar::getChannel(), G4INCL::Particle::isMeson(), G4INCL::CrossSectionsMultiPionsAndResonances::total(), and G4INCL::CrossSectionsStrangeness::total().

◆ isEtaPrime()

G4bool G4INCL::Particle::isEtaPrime ( ) const

inlineinherited

Is this an etaprime?

Definition at line 352 of file G4INCLParticle.hh.

352{ return (theType == EtaPrime); }

References G4INCL::EtaPrime, and G4INCL::Particle::theType.

Referenced by G4INCL::Particle::isMeson(), G4INCL::CrossSectionsMultiPionsAndResonances::total(), and G4INCL::CrossSectionsStrangeness::total().

◆ isEventTransparent()

G4bool G4INCL::Nucleus::isEventTransparent ( ) const

Is the event transparent?

To be called at the end of the cascade.

Definition at line 840 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;
 
  }

References G4INCL::Book::getAcceptedCollisions(), G4INCL::Book::getAcceptedDecays(), G4INCL::Store::getBook(), G4INCL::Book::getEmittedClusters(), and theStore.

Referenced by G4INCL::INCL::postCascade().

◆ isHyperon()

G4bool G4INCL::Particle::isHyperon ( ) const

inlineinherited

Is this an Hyperon?

Definition at line 381 of file G4INCLParticle.hh.

381{ return (isLambda() || isSigma()); }

G4INCL::Particle::isLambda

G4bool isLambda() const

Is this a Lambda?

Definition: G4INCLParticle.hh:375

G4INCL::Particle::isSigma

G4bool isSigma() const

Is this a Sigma?

Definition: G4INCLParticle.hh:366

References G4INCL::Particle::isLambda(), and G4INCL::Particle::isSigma().

Referenced by G4INCL::CrossSectionsStrangeness::elastic(), G4INCL::Particle::isBaryon(), G4INCL::Particle::isStrange(), and G4INCL::CrossSectionsStrangeness::NYelastic().

◆ isKaon()

G4bool G4INCL::Particle::isKaon ( ) const

inlineinherited

Is this a Kaon?

Definition at line 369 of file G4INCLParticle.hh.

369{ return (theType == KPlus || theType == KZero); }

References G4INCL::KPlus, G4INCL::KZero, and G4INCL::Particle::theType.

Referenced by G4INCL::CrossSectionsStrangeness::elastic(), G4INCL::ParticleEntryChannel::fillFinalState(), G4INCL::BinaryCollisionAvatar::getChannel(), G4INCL::SurfaceAvatar::getChannel(), insertParticle(), G4INCL::Particle::isMeson(), G4INCL::Particle::isStrange(), G4INCL::CrossSectionsStrangeness::NKelastic(), G4INCL::CrossSectionsStrangeness::NKToNK(), G4INCL::CrossSectionsStrangeness::NKToNK2pi(), G4INCL::CrossSectionsStrangeness::NKToNKpi(), and G4INCL::CrossSectionsStrangeness::total().

◆ isLambda()

G4bool G4INCL::Particle::isLambda ( ) const

inlineinherited

Is this a Lambda?

Definition at line 375 of file G4INCLParticle.hh.

375{ return (theType == Lambda); }

References G4INCL::Lambda, and G4INCL::Particle::theType.

Referenced by G4INCL::NNToMissingStrangenessChannel::fillFinalState(), G4INCL::NpiToMissingStrangenessChannel::fillFinalState(), G4INCL::BinaryCollisionAvatar::getChannel(), getSurfaceRadius(), G4INCL::Particle::isHyperon(), G4INCL::Particle::isNucleonorLambda(), G4INCL::CrossSectionsStrangeness::NLToNS(), and G4INCL::CrossSectionsStrangeness::total().

◆ isMeson()

G4bool G4INCL::Particle::isMeson ( ) const

inlineinherited

Is this a Meson?

Definition at line 384 of file G4INCLParticle.hh.

384{ return (isPion() || isKaon() || isAntiKaon() || isEta() || isEtaPrime() || isOmega()); }

G4INCL::Particle::isEtaPrime

G4bool isEtaPrime() const

Is this an etaprime?

Definition: G4INCLParticle.hh:352

G4INCL::Particle::isOmega

G4bool isOmega() const

Is this an omega?

Definition: G4INCLParticle.hh:349

G4INCL::Particle::isEta

G4bool isEta() const

Is this an eta?

Definition: G4INCLParticle.hh:346

G4INCL::Particle::isPion

G4bool isPion() const

Is this a pion?

Definition: G4INCLParticle.hh:343

G4INCL::Particle::isAntiKaon

G4bool isAntiKaon() const

Is this an antiKaon?

Definition: G4INCLParticle.hh:372

G4INCL::Particle::isKaon

G4bool isKaon() const

Is this a Kaon?

Definition: G4INCLParticle.hh:369

References G4INCL::Particle::isAntiKaon(), G4INCL::Particle::isEta(), G4INCL::Particle::isEtaPrime(), G4INCL::Particle::isKaon(), G4INCL::Particle::isOmega(), and G4INCL::Particle::isPion().

Referenced by G4INCL::StandardPropagationModel::generateBinaryCollisionAvatar(), G4INCL::InteractionAvatar::preInteractionLocalEnergy(), G4INCL::CDPP::processOneParticle(), and G4INCL::StandardPropagationModel::shootParticle().

◆ isNucleon()

G4bool G4INCL::Particle::isNucleon ( ) const

inlineinherited

Is this a nucleon?

Definition at line 303 of file G4INCLParticle.hh.

                             {
      if(theType == G4INCL::Proton || theType == G4INCL::Neutron)
    return true;
      else
    return false;
    };

References G4INCL::Neutron, G4INCL::Proton, and G4INCL::Particle::theType.

◆ isNucleonorLambda()

G4bool G4INCL::Particle::isNucleonorLambda ( ) const

inlineinherited

Is this a Nucleon or a Lambda?

Definition at line 378 of file G4INCLParticle.hh.

378{ return (isNucleon() || isLambda()); }

References G4INCL::Particle::isLambda(), and G4INCL::Particle::isNucleon().

Referenced by G4INCL::ParticleEntryChannel::fillFinalState(), and G4INCL::SurfaceAvatar::getChannel().

◆ isNucleusNucleusCollision()

G4bool G4INCL::Nucleus::isNucleusNucleusCollision ( ) const

inline

Is it a nucleus-nucleus collision?

Definition at line 424 of file G4INCLNucleus.hh.

424{ return isNucleusNucleus; }

References isNucleusNucleus.

Referenced by G4INCL::ParticleEntryChannel::fillFinalState(), G4INCL::SurfaceAvatar::getChannel(), and G4INCL::INCL::makeCompoundNucleus().

◆ isOmega()

G4bool G4INCL::Particle::isOmega ( ) const

inlineinherited

Is this an omega?

Definition at line 349 of file G4INCLParticle.hh.

349{ return (theType == Omega); }

References G4INCL::Omega, and G4INCL::Particle::theType.

Referenced by G4INCL::CrossSectionsMultiPionsAndResonances::elastic(), G4INCL::BinaryCollisionAvatar::getChannel(), G4INCL::DecayAvatar::getChannel(), G4INCL::Particle::isMeson(), G4INCL::CrossSectionsMultiPionsAndResonances::omegaNElastic(), G4INCL::CrossSectionsMultiPionsAndResonances::omegaNInelastic(), G4INCL::CrossSectionsMultiPionsAndResonances::omegaNToPiN(), G4INCL::CrossSectionsMultiPionsAndResonances::total(), and G4INCL::CrossSectionsStrangeness::total().

◆ isOutOfWell()

G4bool G4INCL::Particle::isOutOfWell ( ) const

inlineinherited

Check if the particle is out of its potential well.

Definition at line 862 of file G4INCLParticle.hh.

862{ return outOfWell; }

G4INCL::Particle::outOfWell

G4bool outOfWell

Definition: G4INCLParticle.hh:1105

References G4INCL::Particle::outOfWell.

Referenced by G4INCL::NuclearPotential::INuclearPotential::computeKaonPotentialEnergy(), G4INCL::NuclearPotential::INuclearPotential::computePionPotentialEnergy(), and G4INCL::NuclearPotential::INuclearPotential::computePionResonancePotentialEnergy().

◆ isParticipant()

G4bool G4INCL::Particle::isParticipant ( ) const

inlineinherited

Definition at line 318 of file G4INCLParticle.hh.

                                 {
      return (theParticipantType==Participant);
    }

References G4INCL::Participant, and G4INCL::Particle::theParticipantType.

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

◆ isPhoton()

G4bool G4INCL::Particle::isPhoton ( ) const

inlineinherited

Is this a photon?

Definition at line 355 of file G4INCLParticle.hh.

355{ return (theType == Photon); }

References G4INCL::Photon, and G4INCL::Particle::theType.

◆ isPion()

G4bool G4INCL::Particle::isPion ( ) const

inlineinherited

Is this a pion?

Definition at line 343 of file G4INCLParticle.hh.

343{ return (theType == PiPlus || theType == PiZero || theType == PiMinus); }

References G4INCL::PiMinus, G4INCL::PiPlus, G4INCL::PiZero, and G4INCL::Particle::theType.

◆ isProjectileSpectator()

G4bool G4INCL::Particle::isProjectileSpectator ( ) const

inlineinherited

Definition at line 326 of file G4INCLParticle.hh.

                                         {
      return (theParticipantType==ProjectileSpectator);
    }

References G4INCL::ProjectileSpectator, and G4INCL::Particle::theParticipantType.

Referenced by G4INCL::SurfaceAvatar::getChannel().

◆ isResonance()

G4bool G4INCL::Particle::isResonance ( ) const

inlineinherited

Is it a resonance?

Definition at line 358 of file G4INCLParticle.hh.

358{ return isDelta(); }

G4INCL::Particle::isDelta

G4bool isDelta() const

Is it a Delta?

Definition: G4INCLParticle.hh:361

References G4INCL::Particle::isDelta().

Referenced by G4INCL::NDeltaToDeltaLKChannel::fillFinalState(), G4INCL::NDeltaToDeltaSKChannel::fillFinalState(), G4INCL::NDeltaToNSKChannel::fillFinalState(), G4INCL::StandardPropagationModel::generateBinaryCollisionAvatar(), G4INCL::SurfaceAvatar::getChannel(), getSurfaceRadius(), G4INCL::Particle::isBaryon(), G4INCL::Particle::Particle(), G4INCL::CDPP::processOneParticle(), and G4INCL::Particle::setType().

◆ isSigma()

G4bool G4INCL::Particle::isSigma ( ) const

inlineinherited

Is this a Sigma?

Definition at line 366 of file G4INCLParticle.hh.

366{ return (theType == SigmaPlus || theType == SigmaZero || theType == SigmaMinus); }

References G4INCL::SigmaMinus, G4INCL::SigmaPlus, G4INCL::SigmaZero, and G4INCL::Particle::theType.

Referenced by G4INCL::BinaryCollisionAvatar::getChannel(), G4INCL::DecayAvatar::getChannel(), G4INCL::Particle::isHyperon(), G4INCL::CrossSectionsStrangeness::NSToNL(), G4INCL::CrossSectionsStrangeness::NSToNS(), and G4INCL::CrossSectionsStrangeness::total().

◆ isStrange()

G4bool G4INCL::Particle::isStrange ( ) const

inlineinherited

Is this an Strange?

Definition at line 390 of file G4INCLParticle.hh.

390{ return (isKaon() || isAntiKaon() || isHyperon()); }

References G4INCL::Particle::isAntiKaon(), G4INCL::Particle::isHyperon(), and G4INCL::Particle::isKaon().

◆ isTargetSpectator()

G4bool G4INCL::Particle::isTargetSpectator ( ) const

inlineinherited

Definition at line 322 of file G4INCLParticle.hh.

                                     {
      return (theParticipantType==TargetSpectator);
    }

References G4INCL::TargetSpectator, and G4INCL::Particle::theParticipantType.

Referenced by G4INCL::SurfaceAvatar::getChannel(), insertParticle(), G4INCL::SurfaceAvatar::postInteraction(), and G4INCL::BinaryCollisionAvatar::preInteraction().

◆ lorentzContract()

void G4INCL::Particle::lorentzContract	(	const ThreeVector &	aBoostVector,
		const ThreeVector &	refPos
	)

inlineinherited

Lorentz-contract the particle position around some center.

Apply Lorentz contraction to the position component along the direction of the boost vector.

Parameters

aBoostVector	the boost vector (velocity) [c]
refPos	the reference position

Definition at line 440 of file G4INCLParticle.hh.

                                                                                     {
      const G4double beta2 = aBoostVector.mag2();
      const G4double gamma = 1.0 / std::sqrt(1.0 - beta2);
      const ThreeVector theRelativePosition = thePosition - refPos;
      const ThreeVector transversePosition = theRelativePosition - aBoostVector * (theRelativePosition.dot(aBoostVector) / aBoostVector.mag2());
      const ThreeVector longitudinalPosition = theRelativePosition - transversePosition;
 
      thePosition = refPos + transversePosition + longitudinalPosition / gamma;
    }

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

◆ makeParticipant()

virtual void G4INCL::Cluster::makeParticipant ( )

inlinevirtualinherited

Make all the components participants, too.

Reimplemented from G4INCL::Particle.

Definition at line 414 of file G4INCLCluster.hh.

                                   {
      Particle::makeParticipant();
      for(ParticleIter p=particles.begin(), e=particles.end(); p!=e; ++p) {
        (*p)->makeParticipant();
      }
    }

References G4INCL::Particle::makeParticipant(), and G4INCL::Cluster::particles.

◆ makeProjectileSpectator()

virtual void G4INCL::Cluster::makeProjectileSpectator ( )

inlinevirtualinherited

Make all the components projectile spectators, too.

Reimplemented from G4INCL::Particle.

Definition at line 398 of file G4INCLCluster.hh.

                                           {
      Particle::makeProjectileSpectator();
      for(ParticleIter p=particles.begin(), e=particles.end(); p!=e; ++p) {
        (*p)->makeProjectileSpectator();
      }
    }

References G4INCL::Particle::makeProjectileSpectator(), and G4INCL::Cluster::particles.

Referenced by G4INCL::ProjectileRemnant::ProjectileRemnant().

◆ makeTargetSpectator()

virtual void G4INCL::Cluster::makeTargetSpectator ( )

inlinevirtualinherited

Make all the components target spectators, too.

Reimplemented from G4INCL::Particle.

Definition at line 406 of file G4INCLCluster.hh.

                                       {
      Particle::makeTargetSpectator();
      for(ParticleIter p=particles.begin(), e=particles.end(); p!=e; ++p) {
        (*p)->makeTargetSpectator();
      }
    }

References G4INCL::Particle::makeTargetSpectator(), and G4INCL::Cluster::particles.

◆ MergeVectorBias() [1/2]

std::vector< G4int > G4INCL::Particle::MergeVectorBias	(	Particle const *const	p1,
		Particle const *const	p2
	)

staticinherited

Definition at line 223 of file G4INCLParticle.cc.

                                                                                                {
    std::vector<G4int> MergedVectorBias;
    std::vector<G4int> VectorBias1 = p1->getBiasCollisionVector();
    std::vector<G4int> VectorBias2 = p2->getBiasCollisionVector();
    G4int i = 0;
    G4int j = 0;
    if(VectorBias1.size()==0 && VectorBias2.size()==0) return MergedVectorBias;
    else if(VectorBias1.size()==0) return VectorBias2;
    else if(VectorBias2.size()==0) return VectorBias1;
        
    while(i < G4int(VectorBias1.size()) || j < G4int(VectorBias2.size())){
        if(VectorBias1[i]==VectorBias2[j]){
            MergedVectorBias.push_back(VectorBias1[i]);
            i++;
            j++;
            if(i == G4int(VectorBias1.size())){
                for(;j<G4int(VectorBias2.size());j++) MergedVectorBias.push_back(VectorBias2[j]);
            }
            else if(j == G4int(VectorBias2.size())){
                for(;i<G4int(VectorBias1.size());i++) MergedVectorBias.push_back(VectorBias1[i]);
            }
        } else if(VectorBias1[i]<VectorBias2[j]){
            MergedVectorBias.push_back(VectorBias1[i]);
            i++;
            if(i == G4int(VectorBias1.size())){
                for(;j<G4int(VectorBias2.size());j++) MergedVectorBias.push_back(VectorBias2[j]);
            }
        }
        else {
            MergedVectorBias.push_back(VectorBias2[j]);
            j++;
            if(j == G4int(VectorBias2.size())){
                for(;i<G4int(VectorBias1.size());i++) MergedVectorBias.push_back(VectorBias1[i]);
            }
        }
    }
    return MergedVectorBias;
  }

References G4INCL::Particle::getBiasCollisionVector().

Referenced by G4INCL::BinaryCollisionAvatar::getChannel(), G4INCL::ParticleList::getParticleListBias(), and G4INCL::ParticleList::getParticleListBiasVector().

◆ MergeVectorBias() [2/2]

std::vector< G4int > G4INCL::Particle::MergeVectorBias	(	std::vector< G4int >	p1,
		Particle const *const	p2
	)

staticinherited

Definition at line 262 of file G4INCLParticle.cc.

                                                                                          {
    std::vector<G4int> MergedVectorBias;
    std::vector<G4int> VectorBias = p2->getBiasCollisionVector();
    G4int i = 0;
    G4int j = 0;
    if(p1.size()==0 && VectorBias.size()==0) return MergedVectorBias;
    else if(p1.size()==0) return VectorBias;
    else if(VectorBias.size()==0) return p1;
        
    while(i < G4int(p1.size()) || j < G4int(VectorBias.size())){
        if(p1[i]==VectorBias[j]){
            MergedVectorBias.push_back(p1[i]);
            i++;
            j++;
            if(i == G4int(p1.size())){
                for(;j<G4int(VectorBias.size());j++) MergedVectorBias.push_back(VectorBias[j]);
            }
            else if(j == G4int(VectorBias.size())){
                for(;i<G4int(p1.size());i++) MergedVectorBias.push_back(p1[i]);
            }
        } else if(p1[i]<VectorBias[j]){
            MergedVectorBias.push_back(p1[i]);
            i++;
            if(i == G4int(p1.size())){
                for(;j<G4int(VectorBias.size());j++) MergedVectorBias.push_back(VectorBias[j]);
            }
        }
        else {
            MergedVectorBias.push_back(VectorBias[j]);
            j++;
            if(j == G4int(VectorBias.size())){
                for(;i<G4int(p1.size());i++) MergedVectorBias.push_back(p1[i]);
            }
        }
    }
    return MergedVectorBias;
  }

References G4INCL::Particle::getBiasCollisionVector().

◆ operator=()

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

Dummy assignment operator to silence Coverity warning.

◆ print() [1/2]

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

Print the nucleus info

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

References G4INCL::Store::getOutgoingParticles(), G4INCL::Store::getParticles(), and theStore.

◆ print() [2/2]

std::string G4INCL::Cluster::print ( ) const

inlineinherited

Definition at line 221 of file G4INCLCluster.hh.

                            {
      std::stringstream ss;
      ss << "Cluster (ID = " << ID << ") type = ";
      ss << ParticleTable::getName(theType);
      ss << '\n'
        << "   A = " << theA << '\n'
        << "   Z = " << theZ << '\n'
        << "   S = " << theS << '\n'
        << "   mass = " << getMass() << '\n'
        << "   energy = " << theEnergy << '\n'
        << "   momentum = "
        << theMomentum.print()
        << '\n'
        << "   position = "
        << thePosition.print()
        << '\n'
        << "Contains the following particles:"
        << '\n';
      for(ParticleIter i=particles.begin(), e=particles.end(); i!=e; ++i)
        ss << (*i)->print();
      ss << '\n';
      return ss.str();
    }

References G4INCL::Particle::getMass(), G4INCL::ParticleTable::getName(), G4INCL::Particle::ID, G4INCL::Cluster::particles, G4INCL::ThreeVector::print(), G4INCL::Particle::theA, G4INCL::Particle::theEnergy, G4INCL::Particle::theMomentum, G4INCL::Particle::thePosition, G4INCL::Particle::theS, G4INCL::Particle::theType, and G4INCL::Particle::theZ.

Referenced by G4INCL::Cluster::boost(), G4INCL::SurfaceAvatar::getChannel(), G4INCL::Cluster::initializeParticles(), G4INCL::Cluster::internalBoostToCM(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::phaseSpaceDecay(), G4INCL::Cluster::putParticlesOffShell(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::recursiveDecay(), G4INCL::ProjectileRemnant::removeParticle(), G4INCL::ProjectileRemnant::reset(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::threeBodyDecay(), and G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::twoBodyDecay().

◆ propagate()

void G4INCL::Particle::propagate ( G4double step )

inlineinherited

Definition at line 829 of file G4INCLParticle.hh.

                                  {
      thePosition += ((*thePropagationMomentum)*(step/(*thePropagationEnergy)));
    };

References G4INCL::Particle::thePosition.

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

◆ propagateParticles()

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

Propagate the particles one time step.

Parameters

step	length of the time step

Definition at line 214 of file G4INCLNucleus.cc.

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

References INCL_WARN.

◆ putParticlesOffShell()

void G4INCL::Cluster::putParticlesOffShell ( )

inlineinherited

Put the cluster components off shell.

The Cluster components are put off shell in such a way that their total energy equals the cluster mass.

Definition at line 306 of file G4INCLCluster.hh.

                                {
      // Compute the dynamical potential
      const G4double theDynamicalPotential = computeDynamicalPotential();
      INCL_DEBUG("The dynamical potential is " << theDynamicalPotential << " MeV" << '\n');
 
      for(ParticleIter p=particles.begin(), e=particles.end(); p!=e; ++p) {
        const G4double energy = (*p)->getEnergy() - theDynamicalPotential;
        const ThreeVector &momentum = (*p)->getMomentum();
        // Here particles are put off-shell so that we can satisfy the energy-
        // and momentum-conservation laws
        (*p)->setEnergy(energy);
        (*p)->setMass(std::sqrt(energy*energy - momentum.mag2()));
      }
      INCL_DEBUG("Cluster components are now off shell:" << '\n'
            << print());
    }

References G4INCL::Cluster::computeDynamicalPotential(), G4INCL::KinematicsUtils::energy(), INCL_DEBUG, G4INCL::ThreeVector::mag2(), G4INCL::Cluster::particles, and G4INCL::Cluster::print().

Referenced by G4INCL::ProjectileRemnant::ProjectileRemnant().

◆ removeParticle()

void G4INCL::Cluster::removeParticle ( Particle *const p )

inlineinherited

Remove a particle from the cluster components.

Definition at line 172 of file G4INCLCluster.hh.

172{ particles.remove(p); }

G4INCL::UnorderedVector::remove

void remove(const T &t)

Definition: G4INCLUnorderedVector.hh:88

References G4INCL::Cluster::particles, and G4INCL::UnorderedVector< T >::remove().

Referenced by G4INCL::ProjectileRemnant::removeParticle().

◆ rotateMomentum()

void G4INCL::Cluster::rotateMomentum	(	const G4double	angle,
		const ThreeVector &	axis
	)

virtualinherited

Rotate momentum of all the particles.

This includes the cluster components. Overloads Particle::rotateMomentum().

Parameters

angle	the rotation angle
axis	a unit vector representing the rotation axis

Reimplemented from G4INCL::Particle.

Definition at line 64 of file G4INCLCluster.cc.

                                                                            {
    Particle::rotateMomentum(angle, axis);
    for(ParticleIter p=particles.begin(), e=particles.end(); p!=e; ++p) {
      (*p)->rotateMomentum(angle, axis);
    }
  }

References angle, G4INCL::Cluster::particles, and G4INCL::Particle::rotateMomentum().

◆ rotatePosition()

void G4INCL::Cluster::rotatePosition	(	const G4double	angle,
		const ThreeVector &	axis
	)

virtualinherited

Rotate position of all the particles.

This includes the cluster components. Overloads Particle::rotateMomentum().

Parameters

angle	the rotation angle
axis	a unit vector representing the rotation axis

Reimplemented from G4INCL::Particle.

Definition at line 57 of file G4INCLCluster.cc.

                                                                            {
    Particle::rotatePosition(angle, axis);
    for(ParticleIter p=particles.begin(), e=particles.end(); p!=e; ++p) {
      (*p)->rotatePosition(angle, axis);
    }
  }

References angle, G4INCL::Cluster::particles, and G4INCL::Particle::rotatePosition().

◆ rotatePositionAndMomentum()

virtual void G4INCL::Particle::rotatePositionAndMomentum	(	const G4double	angle,
		const ThreeVector &	axis
	)

inlinevirtualinherited

Rotate the particle position and momentum.

Parameters

angle	the rotation angle
axis	a unit vector representing the rotation axis

Definition at line 929 of file G4INCLParticle.hh.

                                                                                          {
      rotatePosition(angle, axis);
      rotateMomentum(angle, axis);
    }

References angle, G4INCL::Particle::rotateMomentum(), and G4INCL::Particle::rotatePosition().

Referenced by G4INCL::CoulombNonRelativistic::coulombDeviation().

◆ rpCorrelate()

void G4INCL::Particle::rpCorrelate ( )

inlineinherited

Make the particle follow a strict r-p correlation.

Definition at line 1008 of file G4INCLParticle.hh.

1008{ rpCorrelated = true; }

References G4INCL::Particle::rpCorrelated.

Referenced by G4INCL::InteractionAvatar::bringParticleInside(), G4INCL::ParticleEntryAvatar::postInteraction(), and G4INCL::SurfaceAvatar::postInteraction().

◆ rpDecorrelate()

void G4INCL::Particle::rpDecorrelate ( )

inlineinherited

Make the particle not follow a strict r-p correlation.

Definition at line 1011 of file G4INCLParticle.hh.

1011{ rpCorrelated = false; }

References G4INCL::Particle::rpCorrelated.

◆ setA()

void G4INCL::Cluster::setA ( const G4int A )

inlineinherited

Set the mass number of the cluster.

Definition at line 149 of file G4INCLCluster.hh.

149{ theA = A; }

A

const G4double A[17]

Definition: G4WaterStopping.cc:53

References A, and G4INCL::Particle::theA.

Referenced by G4INCL::INCL::makeCompoundNucleus(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::phaseSpaceDecay(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::threeBodyDecay(), and G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::twoBodyDecay().

◆ setBiasCollisionVector()

void G4INCL::Particle::setBiasCollisionVector ( std::vector< G4int > BiasCollisionVector )

inlineinherited

Set the vector list of biased vertices on the particle path.

Definition at line 1041 of file G4INCLParticle.hh.

                                                                      {
      this->theBiasCollisionVector = BiasCollisionVector;
      this->setParticleBias(Particle::getBiasFromVector(BiasCollisionVector));
      }

References G4INCL::Particle::getBiasFromVector(), G4INCL::Particle::setParticleBias(), and G4INCL::Particle::theBiasCollisionVector.

Referenced by decayOutgoingPionResonances(), decayOutgoingSigmaZero(), and G4INCL::SurfaceAvatar::postInteraction().

◆ setDensity()

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

inline

Setter for theDensity.

Definition at line 463 of file G4INCLNucleus.hh.

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

References G4INCL::ParticleSampler::setDensity(), theDensity, and G4INCL::Cluster::theParticleSampler.

◆ setEmissionTime()

void G4INCL::Particle::setEmissionTime ( G4double t )

inlineinherited

Definition at line 864 of file G4INCLParticle.hh.

864{ emissionTime = t; }

References G4INCL::Particle::emissionTime.

Referenced by computeOneNucleonRecoilKinematics(), emitInsideKaon(), emitInsideLambda(), emitInsidePions(), emitInsideStrangeParticles(), finalizeProjectileRemnant(), and G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::twoBodyDecay().

◆ setEnergy()

void G4INCL::Particle::setEnergy ( G4double energy )

inlineinherited

Set the energy of the particle in MeV.

Definition at line 786 of file G4INCLParticle.hh.

    {
      this->theEnergy = energy;
    };

References G4INCL::KinematicsUtils::energy(), and G4INCL::Particle::theEnergy.

Referenced by emitInsideKaon(), emitInsideLambda(), emitInsidePions(), emitInsideStrangeParticles(), G4INCL::DeltaProductionChannel::fillFinalState(), G4INCL::EtaNElasticChannel::fillFinalState(), G4INCL::EtaNToPiNChannel::fillFinalState(), G4INCL::NSToNLChannel::fillFinalState(), G4INCL::NSToNSChannel::fillFinalState(), G4INCL::OmegaNElasticChannel::fillFinalState(), G4INCL::OmegaNToPiNChannel::fillFinalState(), G4INCL::CrossSections::interactionDistanceKbarN(), G4INCL::CrossSections::interactionDistanceKN(), G4INCL::CrossSections::interactionDistanceNN(), G4INCL::CrossSections::interactionDistancePiN(), G4INCL::CrossSections::interactionDistanceYN(), G4INCL::INCL::makeCompoundNucleus(), G4INCL::ParticleEntryChannel::particleEnters(), G4INCL::TransmissionChannel::particleLeaves(), G4INCL::INCL::RecoilCMFunctor::scaleParticleCMMomenta(), G4INCL::INCL::RecoilFunctor::scaleParticleEnergies(), G4INCL::StandardPropagationModel::shootParticle(), G4INCL::KinematicsUtils::transformToLocalEnergyFrame(), and useFusionKinematics().

◆ setExcitationEnergy()

void G4INCL::Cluster::setExcitationEnergy ( const G4double e )

inlineinherited

Set the excitation energy of the cluster.

Definition at line 158 of file G4INCLCluster.hh.

158{ theExcitationEnergy=e; }

References G4INCL::Cluster::theExcitationEnergy.

Referenced by finalizeProjectileRemnant(), G4INCL::INCL::makeCompoundNucleus(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::phaseSpaceDecay(), G4INCL::INCL::postCascade(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::recursiveDecay(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::threeBodyDecay(), and G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::twoBodyDecay().

◆ setFrozenEnergy()

void G4INCL::Particle::setFrozenEnergy ( const G4double energy )

inlineinherited

Set the frozen particle momentum.

Definition at line 891 of file G4INCLParticle.hh.

891{ theFrozenEnergy = energy; }

References G4INCL::KinematicsUtils::energy(), and G4INCL::Particle::theFrozenEnergy.

◆ setFrozenMomentum()

void G4INCL::Particle::setFrozenMomentum ( const ThreeVector & momentum )

inlineinherited

Set the frozen particle momentum.

Definition at line 888 of file G4INCLParticle.hh.

888{ theFrozenMomentum = momentum; }

References G4INCL::Particle::theFrozenMomentum.

◆ setHelicity()

void G4INCL::Particle::setHelicity ( G4double h )

inlineinherited

Definition at line 827 of file G4INCLParticle.hh.

827{ theHelicity = h; };

References G4INCL::Particle::theHelicity.

Referenced by G4INCL::DeltaDecayChannel::fillFinalState(), G4INCL::DeltaProductionChannel::fillFinalState(), G4INCL::NDeltaEtaProductionChannel::fillFinalState(), and G4INCL::NDeltaOmegaProductionChannel::fillFinalState().

◆ setINCLBiasVector()

void G4INCL::Particle::setINCLBiasVector ( std::vector< G4double > NewVector )

staticinherited

Definition at line 306 of file G4INCLParticle.cc.

                                                                {
      Particle::INCLBiasVector = NewVector;
  }

References G4INCL::Particle::INCLBiasVector.

◆ setINCLMass()

void G4INCL::Particle::setINCLMass ( )

inlineinherited

Set the mass of the Particle to its table mass.

Definition at line 589 of file G4INCLParticle.hh.

589{ setMass(getINCLMass()); }

References G4INCL::Particle::getINCLMass(), and G4INCL::Particle::setMass().

Referenced by G4INCL::Cluster::Cluster(), G4INCL::ParticleEntryChannel::particleEnters(), G4INCL::Particle::setType(), and G4INCL::StandardPropagationModel::shootParticle().

◆ setIncomingAngularMomentum()

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

inline

Set the incoming angular-momentum vector.

Definition at line 257 of file G4INCLNucleus.hh.

                                                          {
      incomingAngularMomentum = j;
    }

References incomingAngularMomentum.

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

◆ setIncomingMomentum()

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

inline

Set the incoming momentum vector.

Definition at line 265 of file G4INCLNucleus.hh.

                                                   {
      incomingMomentum = p;
    }

References incomingMomentum.

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

◆ setInitialEnergy()

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

inline

Set the initial energy.

Definition at line 275 of file G4INCLNucleus.hh.

275{ initialEnergy = e; }

References initialEnergy.

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

◆ setMass()

void G4INCL::Particle::setMass ( G4double mass )

inlineinherited

Set the mass of the particle in MeV/c^2.

Definition at line 778 of file G4INCLParticle.hh.

    {
      this->theMass = mass;
    }

References G4INCL::Particle::theMass.

Referenced by computeRecoilKinematics(), G4INCL::DeltaProductionChannel::fillFinalState(), G4INCL::NDeltaEtaProductionChannel::fillFinalState(), G4INCL::NDeltaOmegaProductionChannel::fillFinalState(), G4INCL::NDeltaToDeltaLKChannel::fillFinalState(), G4INCL::NDeltaToDeltaSKChannel::fillFinalState(), finalizeProjectileRemnant(), G4INCL::PhaseSpaceKopylov::generate(), G4INCL::INCL::makeCompoundNucleus(), G4INCL::Particle::Particle(), G4INCL::Particle::setINCLMass(), G4INCL::Particle::setRealMass(), G4INCL::Particle::setTableMass(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::threeBodyDecay(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::twoBodyDecay(), and useFusionKinematics().

◆ setMomentum()

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

inlinevirtualinherited

Set the momentum vector.

Definition at line 808 of file G4INCLParticle.hh.

    {
      this->theMomentum = momentum;
    };

References G4INCL::Particle::theMomentum.

◆ setNucleusNucleusCollision()

void G4INCL::Nucleus::setNucleusNucleusCollision ( )

inline

Set a nucleus-nucleus collision.

Definition at line 427 of file G4INCLNucleus.hh.

427{ isNucleusNucleus=true; }

References isNucleusNucleus.

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

◆ setNumberOfCollisions()

void G4INCL::Particle::setNumberOfCollisions ( G4int n )

inlineinherited

Set the number of collisions undergone by the particle.

Definition at line 837 of file G4INCLParticle.hh.

837{ nCollisions = n; }

CLHEP::detail::n

n

Definition: Ranlux64Engine.cc:90

References CLHEP::detail::n, and G4INCL::Particle::nCollisions.

◆ setNumberOfDecays()

void G4INCL::Particle::setNumberOfDecays ( G4int n )

inlineinherited

Set the number of decays undergone by the particle.

Definition at line 846 of file G4INCLParticle.hh.

846{ nDecays = n; }

References CLHEP::detail::n, and G4INCL::Particle::nDecays.

Referenced by G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::threeBodyDecay(), and G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::twoBodyDecay().

◆ setNumberOfKaon()

void G4INCL::Particle::setNumberOfKaon ( const G4int NK )

inlineinherited

Definition at line 1054 of file G4INCLParticle.hh.

1054{ theNKaon = NK; }

References G4INCL::Particle::theNKaon.

Referenced by G4INCL::ParticleEntryChannel::fillFinalState(), G4INCL::SurfaceAvatar::getChannel(), and G4INCL::BinaryCollisionAvatar::postInteraction().

◆ setOutOfWell()

void G4INCL::Particle::setOutOfWell ( )

inlineinherited

Mark the particle as out of its potential well.

This flag is used to control pions created outside their potential well in delta decay. The pion potential checks it and returns zero if it is true (necessary in order to correctly enforce energy conservation). The Nucleus::applyFinalState() method uses it to determine whether new avatars should be generated for the particle.

Definition at line 859 of file G4INCLParticle.hh.

859{ outOfWell = true; }

References G4INCL::Particle::outOfWell.

◆ setParticipantType()

void G4INCL::Particle::setParticipantType ( ParticipantType const p )

inlineinherited

Definition at line 314 of file G4INCLParticle.hh.

                                                     {
      theParticipantType = p;
    }

References G4INCL::Particle::theParticipantType.

◆ setParticleBias()

void G4INCL::Particle::setParticleBias ( G4double ParticleBias )

inlineinherited

Set the particle bias.

Definition at line 1035 of file G4INCLParticle.hh.

1035{ this->theParticleBias = ParticleBias; }

References G4INCL::Particle::theParticleBias.

Referenced by G4INCL::Particle::setBiasCollisionVector().

◆ setParticleNucleusCollision()

void G4INCL::Nucleus::setParticleNucleusCollision ( )

inline

Set a particle-nucleus collision.

Definition at line 430 of file G4INCLNucleus.hh.

430{ isNucleusNucleus=false; }

References isNucleusNucleus.

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

◆ setPosition()

void G4INCL::Cluster::setPosition ( const ThreeVector & position )

inlinevirtualinherited

Set the position of the cluster.

This overloads the Particle method to take into account that the positions of the cluster members must be updated as well.

Reimplemented from G4INCL::Particle.

Definition at line 328 of file G4INCLCluster.hh.

                                                  {
      ThreeVector shift(position-thePosition);
      Particle::setPosition(position);
      for(ParticleIter p=particles.begin(), e=particles.end(); p!=e; ++p) {
        (*p)->setPosition((*p)->getPosition()+shift);
      }
    }

References G4INCL::Cluster::particles, G4INCL::Particle::setPosition(), and G4INCL::Particle::thePosition.

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

◆ setPotentialEnergy()

void G4INCL::Particle::setPotentialEnergy ( G4double v )

inlineinherited

Set the particle potential energy.

Definition at line 765 of file G4INCLParticle.hh.

765{ thePotentialEnergy = v; }

References G4INCL::Particle::thePotentialEnergy.

Referenced by emitInsideKaon(), emitInsideLambda(), emitInsidePions(), emitInsideStrangeParticles(), G4INCL::Store::loadParticles(), G4INCL::ParticleEntryChannel::particleEnters(), G4INCL::TransmissionChannel::particleLeaves(), and updatePotentialEnergy().

◆ setProjectileRemnant()

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

inline

Set the projectile remnant.

Definition at line 433 of file G4INCLNucleus.hh.

                                                           {
      delete theProjectileRemnant;
      theProjectileRemnant = c;
    }

References theProjectileRemnant.

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

◆ setRealMass()

void G4INCL::Particle::setRealMass ( )

inlineinherited

Set the mass of the Particle to its real mass.

Definition at line 583 of file G4INCLParticle.hh.

583{ setMass(getRealMass()); }

References G4INCL::Particle::getRealMass(), and G4INCL::Particle::setMass().

Referenced by G4INCL::ClusterDecay::decay(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::phaseSpaceDecay(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::threeBodyDecay(), and G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::twoBodyDecay().

◆ setS()

void G4INCL::Cluster::setS ( const G4int S )

inlineinherited

Set the strangess number of the cluster.

Definition at line 152 of file G4INCLCluster.hh.

152{ theS = S; }

References S(), and G4INCL::Particle::theS.

Referenced by G4INCL::INCL::makeCompoundNucleus(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::phaseSpaceDecay(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::threeBodyDecay(), and G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::twoBodyDecay().

◆ setSpin()

void G4INCL::Cluster::setSpin ( const ThreeVector & j )

inlineinherited

Set the spin of the nucleus.

Definition at line 425 of file G4INCLCluster.hh.

425{ theSpin = j; }

References G4INCL::Cluster::theSpin.

Referenced by finalizeProjectileRemnant(), G4INCL::INCL::makeCompoundNucleus(), and useFusionKinematics().

◆ setStore()

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

inline

Definition at line 348 of file G4INCLNucleus.hh.

                              {
      delete theStore;
      theStore = str;
    };

References theStore.

◆ setTableMass()

void G4INCL::Particle::setTableMass ( )

inlineinherited

Set the mass of the Particle to its table mass.

Definition at line 586 of file G4INCLParticle.hh.

586{ setMass(getTableMass()); }

References G4INCL::Particle::getTableMass(), and G4INCL::Particle::setMass().

Referenced by decayOutgoingPionResonances(), decayOutgoingSigmaZero(), emitInsideKaon(), emitInsideLambda(), emitInsidePions(), emitInsideStrangeParticles(), G4INCL::TransmissionChannel::particleLeaves(), G4INCL::ProjectileRemnant::ProjectileRemnant(), and G4INCL::ProjectileRemnant::reset().

◆ setType()

void G4INCL::Particle::setType ( ParticleType t )

inlineinherited

Definition at line 187 of file G4INCLParticle.hh.

                                 {
      theType = t;
      switch(theType)
      {
        case DeltaPlusPlus:
          theA = 1;
          theZ = 2;
          theS = 0;
          break;
        case Proton:
        case DeltaPlus:
          theA = 1;
          theZ = 1;
          theS = 0;
          break;
        case Neutron:
        case DeltaZero:
          theA = 1;
          theZ = 0;
          theS = 0;
          break;
        case DeltaMinus:
          theA = 1;
          theZ = -1;
          theS = 0;
          break;
        case PiPlus:
          theA = 0;
          theZ = 1;
          theS = 0;
          break;
        case PiZero:
        case Eta:
        case Omega:
        case EtaPrime:
        case Photon:
          theA = 0;
          theZ = 0;
          theS = 0;
          break;
        case PiMinus:
          theA = 0;
          theZ = -1;
          theS = 0;
          break;
        case Lambda:
          theA = 1;
          theZ = 0;
          theS = -1;
          break;
        case SigmaPlus:
          theA = 1;
          theZ = 1;
          theS = -1;
          break;
        case SigmaZero:
          theA = 1;
          theZ = 0;
          theS = -1;
          break;
        case SigmaMinus:
          theA = 1;
          theZ = -1;
          theS = -1;
          break;
        case KPlus:
          theA = 0;
          theZ = 1;
          theS = 1;
          break;
        case KZero:
          theA = 0;
          theZ = 0;
          theS = 1;
          break;
        case KZeroBar:
          theA = 0;
          theZ = 0;
          theS = -1;
          break;
        case KShort:
          theA = 0;
          theZ = 0;
//        theS should not be defined
          break;
        case KLong:
          theA = 0;
          theZ = 0;
//        theS should not be defined
          break;
        case KMinus:
          theA = 0;
          theZ = -1;
          theS = -1;
          break;
        case Composite:
         // INCL_ERROR("Trying to set particle type to Composite! Construct a Cluster object instead" << '\n');
          theA = 0;
          theZ = 0;
          theS = 0;
          break;
        case UnknownParticle:
          theA = 0;
          theZ = 0;
          theS = 0;
          INCL_ERROR("Trying to set particle type to Unknown!" << '\n');
          break;
      }
 
      if( !isResonance() && t!=Composite )
        setINCLMass();
    }

References G4INCL::Composite, G4INCL::DeltaMinus, G4INCL::DeltaPlus, G4INCL::DeltaPlusPlus, G4INCL::DeltaZero, G4INCL::Eta, G4INCL::EtaPrime, INCL_ERROR, G4INCL::Particle::isResonance(), G4INCL::KLong, G4INCL::KMinus, G4INCL::KPlus, G4INCL::KShort, G4INCL::KZero, G4INCL::KZeroBar, G4INCL::Lambda, G4INCL::Neutron, G4INCL::Omega, G4INCL::Photon, G4INCL::PiMinus, G4INCL::PiPlus, G4INCL::PiZero, G4INCL::Proton, G4INCL::Particle::setINCLMass(), G4INCL::SigmaMinus, G4INCL::SigmaPlus, G4INCL::SigmaZero, G4INCL::Particle::theA, G4INCL::Particle::theS, G4INCL::Particle::theType, G4INCL::Particle::theZ, and G4INCL::UnknownParticle.

◆ setUncorrelatedMomentum()

void G4INCL::Particle::setUncorrelatedMomentum ( const G4double p )

inlineinherited

Set the uncorrelated momentum.

Definition at line 1005 of file G4INCLParticle.hh.

1005{ uncorrelatedMomentum = p; }

References G4INCL::Particle::uncorrelatedMomentum.

Referenced by G4INCL::ParticleSampler::sampleOneParticleWithFuzzyRPCorrelation(), and G4INCL::ParticleSampler::sampleOneParticleWithRPCorrelation().

◆ setUniverseRadius()

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

inline

Setter for theUniverseRadius.

Definition at line 421 of file G4INCLNucleus.hh.

421{ theUniverseRadius=universeRadius; }

References theUniverseRadius.

◆ setZ()

void G4INCL::Cluster::setZ ( const G4int Z )

inlineinherited

Set the charge number of the cluster.

Definition at line 146 of file G4INCLCluster.hh.

146{ theZ = Z; }

Z

const G4int Z[17]

Definition: G4WaterStopping.cc:51

References G4INCL::Particle::theZ, and Z.

Referenced by G4INCL::INCL::makeCompoundNucleus(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::phaseSpaceDecay(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::threeBodyDecay(), and G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::twoBodyDecay().

◆ swap() [1/2]

void G4INCL::Cluster::swap ( Cluster & rhs )

inlineinherited

Helper method for the assignment operator.

Definition at line 122 of file G4INCLCluster.hh.

                            {
      Particle::swap(rhs);
      std::swap(theExcitationEnergy, rhs.theExcitationEnergy);
      std::swap(theSpin, rhs.theSpin);
      // std::swap is overloaded by std::list and guaranteed to operate in
      // constant time
      std::swap(particles, rhs.particles);
      std::swap(theParticleSampler, rhs.theParticleSampler);
    }

References G4INCL::Cluster::particles, G4INCL::Particle::swap(), G4INCL::Cluster::theExcitationEnergy, G4INCL::Cluster::theParticleSampler, and G4INCL::Cluster::theSpin.

Referenced by G4INCL::Cluster::operator=().

◆ swap() [2/2]

void G4INCL::Particle::swap ( Particle & rhs )

inlineprotectedinherited

Helper method for the assignment operator.

Definition at line 125 of file G4INCLParticle.hh.

                             {
      std::swap(theZ, rhs.theZ);
      std::swap(theA, rhs.theA);
      std::swap(theS, rhs.theS);
      std::swap(theParticipantType, rhs.theParticipantType);
      std::swap(theType, rhs.theType);
      if(rhs.thePropagationEnergy == &(rhs.theFrozenEnergy))
        thePropagationEnergy = &theFrozenEnergy;
      else
        thePropagationEnergy = &theEnergy;
      std::swap(theEnergy, rhs.theEnergy);
      std::swap(theFrozenEnergy, rhs.theFrozenEnergy);
      if(rhs.thePropagationMomentum == &(rhs.theFrozenMomentum))
        thePropagationMomentum = &theFrozenMomentum;
      else
        thePropagationMomentum = &theMomentum;
      std::swap(theMomentum, rhs.theMomentum);
      std::swap(theFrozenMomentum, rhs.theFrozenMomentum);
      std::swap(thePosition, rhs.thePosition);
      std::swap(nCollisions, rhs.nCollisions);
      std::swap(nDecays, rhs.nDecays);
      std::swap(thePotentialEnergy, rhs.thePotentialEnergy);
      // ID intentionally not swapped
 
      std::swap(theHelicity, rhs.theHelicity);
      std::swap(emissionTime, rhs.emissionTime);
      std::swap(outOfWell, rhs.outOfWell);
 
      std::swap(theMass, rhs.theMass);
      std::swap(rpCorrelated, rhs.rpCorrelated);
      std::swap(uncorrelatedMomentum, rhs.uncorrelatedMomentum);
      
      std::swap(theParticleBias, rhs.theParticleBias);
      std::swap(theBiasCollisionVector, rhs.theBiasCollisionVector);
 
    }

Referenced by G4INCL::Particle::operator=(), and G4INCL::Cluster::swap().

◆ thawPropagation()

void G4INCL::Particle::thawPropagation ( )

inlineinherited

Unfreeze particle propagation.

Make the particle use theMomentum and theEnergy for propagation. Call this method to restore the normal propagation if the freezePropagation() method has been called.

Definition at line 919 of file G4INCLParticle.hh.

                           {
      thePropagationMomentum = &theMomentum;
      thePropagationEnergy = &theEnergy;
    }

References G4INCL::Particle::theEnergy, G4INCL::Particle::theMomentum, G4INCL::Particle::thePropagationEnergy, and G4INCL::Particle::thePropagationMomentum.

Referenced by G4INCL::ReflectionChannel::fillFinalState().

◆ updateClusterParameters()

void G4INCL::Cluster::updateClusterParameters ( )

inlineinherited

Set total cluster mass, energy, size, etc. from the particles.

Definition at line 191 of file G4INCLCluster.hh.

                                   {
      theEnergy = 0.;
      thePotentialEnergy = 0.;
      theMomentum = ThreeVector();
      thePosition = ThreeVector();
      theA = 0;
      theZ = 0;
      theS = 0;
      nCollisions = 0;
      for(ParticleIter p=particles.begin(), e=particles.end(); p!=e; ++p) {
        theEnergy += (*p)->getEnergy();
        thePotentialEnergy += (*p)->getPotentialEnergy();
        theMomentum += (*p)->getMomentum();
        thePosition += (*p)->getPosition();
        theA += (*p)->getA();
        theZ += (*p)->getZ();
        theS += (*p)->getS();
        nCollisions += (*p)->getNumberOfCollisions();
      }
    }

References G4INCL::Particle::nCollisions, G4INCL::Cluster::particles, G4INCL::Particle::theA, G4INCL::Particle::theEnergy, G4INCL::Particle::theMomentum, G4INCL::Particle::thePosition, G4INCL::Particle::thePotentialEnergy, G4INCL::Particle::theS, and G4INCL::Particle::theZ.

Referenced by G4INCL::Cluster::addParticles(), and G4INCL::Cluster::initializeParticles().

◆ updatePotentialEnergy()

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

inline

Update the particle potential energy.

Definition at line 458 of file G4INCLNucleus.hh.

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

References G4INCL::NuclearPotential::INuclearPotential::computePotentialEnergy(), G4INCL::Particle::setPotentialEnergy(), and thePotential.

Referenced by G4INCL::ReflectionChannel::fillFinalState(), initializeParticles(), G4INCL::InteractionAvatar::ViolationEMomentumFunctor::scaleParticleMomenta(), and G4INCL::InteractionAvatar::ViolationEEnergyFunctor::setParticleEnergy().

◆ useFusionKinematics()

void G4INCL::Nucleus::useFusionKinematics ( )

Adjust the kinematics for complete-fusion events.

Definition at line 1203 of file G4INCLNucleus.cc.

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

References G4INCL::Cluster::getTableMass(), incomingAngularMomentum, incomingMomentum, initialEnergy, G4INCL::ThreeVector::mag2(), G4INCL::Particle::setEnergy(), G4INCL::Particle::setMass(), G4INCL::Particle::setMomentum(), G4INCL::Cluster::setSpin(), G4INCL::Particle::theEnergy, G4INCL::Cluster::theExcitationEnergy, and G4INCL::Particle::theMomentum.

Referenced by G4INCL::INCL::postCascade().

Field Documentation

◆ emissionTime

G4double G4INCL::Particle::emissionTime

privateinherited

Definition at line 1104 of file G4INCLParticle.hh.

Referenced by G4INCL::Particle::getEmissionTime(), G4INCL::Particle::setEmissionTime(), and G4INCL::Particle::swap().

◆ ID

long G4INCL::Particle::ID

protectedinherited

Definition at line 1093 of file G4INCLParticle.hh.

Referenced by G4INCL::Particle::dump(), G4INCL::Particle::getID(), G4INCL::Particle::Particle(), G4INCL::Cluster::print(), and G4INCL::Particle::print().

◆ INCLBiasVector

std::vector< G4double > G4INCL::Particle::INCLBiasVector

staticinherited

Time ordered vector of all bias applied.

/!\ Caution /!\ methods Assotiated to G4VectorCache<T> are: Push_back(…), operator[], Begin(), End(), Clear(), Size() and Pop_back()

Definition at line 1071 of file G4INCLParticle.hh.

Referenced by G4INCL::Particle::FillINCLBiasVector(), G4INCL::Particle::getBiasFromVector(), G4INCL::Particle::getTotalBias(), G4INCL::INCL::processEvent(), and G4INCL::Particle::setINCLBiasVector().

◆ incomingAngularMomentum

ThreeVector G4INCL::Nucleus::incomingAngularMomentum

private

Definition at line 497 of file G4INCLNucleus.hh.

Referenced by computeRecoilKinematics(), getIncomingAngularMomentum(), setIncomingAngularMomentum(), and useFusionKinematics().

◆ incomingMomentum

ThreeVector G4INCL::Nucleus::incomingMomentum

private

Definition at line 497 of file G4INCLNucleus.hh.

Referenced by computeOneNucleonRecoilKinematics(), computeRecoilKinematics(), getIncomingMomentum(), setIncomingMomentum(), and useFusionKinematics().

◆ initialCenterOfMass

ThreeVector G4INCL::Nucleus::initialCenterOfMass

private

Definition at line 498 of file G4INCLNucleus.hh.

Referenced by computeRecoilKinematics(), and initializeParticles().

◆ initialEnergy

G4double G4INCL::Nucleus::initialEnergy

private

Definition at line 501 of file G4INCLNucleus.hh.

Referenced by computeOneNucleonRecoilKinematics(), getInitialEnergy(), setInitialEnergy(), and useFusionKinematics().

◆ initialInternalEnergy

G4double G4INCL::Nucleus::initialInternalEnergy

private

Definition at line 496 of file G4INCLNucleus.hh.

Referenced by computeExcitationEnergy(), getInitialInternalEnergy(), and initializeParticles().

◆ 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 519 of file G4INCLNucleus.hh.

Referenced by isNucleusNucleusCollision(), setNucleusNucleusCollision(), and setParticleNucleusCollision().

◆ nCollisions

G4int G4INCL::Particle::nCollisions

protectedinherited

Definition at line 1090 of file G4INCLParticle.hh.

Referenced by G4INCL::Cluster::addParticle(), G4INCL::Particle::getNumberOfCollisions(), G4INCL::Particle::incrementNumberOfCollisions(), G4INCL::ProjectileRemnant::reset(), G4INCL::Particle::setNumberOfCollisions(), G4INCL::Particle::swap(), and G4INCL::Cluster::updateClusterParameters().

◆ nDecays

G4int G4INCL::Particle::nDecays

protectedinherited

Definition at line 1091 of file G4INCLParticle.hh.

Referenced by G4INCL::Particle::getNumberOfDecays(), G4INCL::Particle::incrementNumberOfDecays(), G4INCL::Particle::setNumberOfDecays(), and G4INCL::Particle::swap().

◆ nextBiasedCollisionID

G4ThreadLocal G4int G4INCL::Particle::nextBiasedCollisionID = 0

staticinherited

Definition at line 1077 of file G4INCLParticle.hh.

Referenced by G4INCL::Particle::FillINCLBiasVector(), G4INCL::InteractionAvatar::postInteraction(), and G4INCL::INCL::processEvent().

◆ nextID

G4ThreadLocal long G4INCL::Particle::nextID = 1

staticprivateinherited

Definition at line 1111 of file G4INCLParticle.hh.

Referenced by G4INCL::Particle::Particle().

◆ outOfWell

G4bool G4INCL::Particle::outOfWell

privateinherited

Definition at line 1105 of file G4INCLParticle.hh.

Referenced by G4INCL::Particle::isOutOfWell(), G4INCL::Particle::setOutOfWell(), and G4INCL::Particle::swap().

◆ particles

ParticleList G4INCL::Cluster::particles

protectedinherited

Definition at line 451 of file G4INCLCluster.hh.

◆ projectileA

G4int G4INCL::Nucleus::projectileA

private

The mass number of the projectile.

Definition at line 508 of file G4INCLNucleus.hh.

◆ projectileS

G4int G4INCL::Nucleus::projectileS

private

The strangeness number of the projectile.

Definition at line 510 of file G4INCLNucleus.hh.

◆ projectileZ

G4int G4INCL::Nucleus::projectileZ

private

The charge number of the projectile.

Definition at line 506 of file G4INCLNucleus.hh.

◆ remnant

G4bool G4INCL::Nucleus::remnant

private

Definition at line 499 of file G4INCLNucleus.hh.

Referenced by computeRecoilKinematics(), and hasRemnant().

◆ rpCorrelated

G4bool G4INCL::Particle::rpCorrelated

protectedinherited

Definition at line 1095 of file G4INCLParticle.hh.

Referenced by G4INCL::Particle::getReflectionMomentum(), G4INCL::Particle::rpCorrelate(), G4INCL::Particle::rpDecorrelate(), and G4INCL::Particle::swap().

◆ theA

G4int G4INCL::Particle::theA

protectedinherited

Definition at line 1080 of file G4INCLParticle.hh.

Referenced by G4INCL::ProjectileRemnant::addAllDynamicalSpectators(), G4INCL::ProjectileRemnant::addDynamicalSpectator(), G4INCL::ProjectileRemnant::addMostDynamicalSpectators(), G4INCL::Cluster::addParticle(), applyFinalState(), G4INCL::Cluster::Cluster(), G4INCL::Cluster::computeDynamicalPotential(), computeOneNucleonRecoilKinematics(), computeRecoilKinematics(), decayInsideDeltas(), decayInsideStrangeParticles(), decayMe(), emitInsideKaon(), emitInsideLambda(), emitInsidePions(), emitInsideStrangeParticles(), G4INCL::Particle::getA(), G4INCL::Particle::getEmissionQValueCorrection(), G4INCL::Particle::getINCLMass(), G4INCL::Particle::getRealMass(), G4INCL::Cluster::getSpecies(), G4INCL::Particle::getTableMass(), G4INCL::Particle::getTransferQValueCorrection(), G4INCL::Cluster::initializeParticles(), insertParticle(), G4INCL::Cluster::internalBoostToCM(), Nucleus(), G4INCL::Cluster::print(), G4INCL::ProjectileRemnant::removeParticle(), G4INCL::ProjectileRemnant::reset(), G4INCL::Cluster::setA(), G4INCL::Particle::setType(), G4INCL::Particle::swap(), and G4INCL::Cluster::updateClusterParameters().

◆ theBiasCollisionVector

std::vector<G4int> G4INCL::Particle::theBiasCollisionVector

privateinherited

Time ordered vector of all biased vertices on the particle path.

Definition at line 1108 of file G4INCLParticle.hh.

Referenced by G4INCL::Particle::getBiasCollisionVector(), G4INCL::Particle::Particle(), G4INCL::Particle::setBiasCollisionVector(), and G4INCL::Particle::swap().

◆ theDensity

NuclearDensity const* G4INCL::Nucleus::theDensity

private

Pointer to the NuclearDensity object.

Definition at line 528 of file G4INCLNucleus.hh.

Referenced by getDensity(), getSurfaceRadius(), getTransmissionBarrier(), Nucleus(), and setDensity().

◆ theEnergy

G4double G4INCL::Particle::theEnergy

protectedinherited

Definition at line 1083 of file G4INCLParticle.hh.

◆ theExcitationEnergy

G4double G4INCL::Cluster::theExcitationEnergy

protectedinherited

Definition at line 452 of file G4INCLCluster.hh.

Referenced by computeOneNucleonRecoilKinematics(), computeRecoilKinematics(), G4INCL::Cluster::getExcitationEnergy(), getExcitationEnergy(), G4INCL::Cluster::setExcitationEnergy(), G4INCL::Cluster::swap(), and useFusionKinematics().

◆ theFrozenEnergy

G4double G4INCL::Particle::theFrozenEnergy

protectedinherited

Definition at line 1085 of file G4INCLParticle.hh.

Referenced by G4INCL::Particle::freezePropagation(), G4INCL::Particle::getFrozenEnergy(), G4INCL::Particle::Particle(), G4INCL::Particle::setFrozenEnergy(), and G4INCL::Particle::swap().

◆ theFrozenMomentum

G4INCL::ThreeVector G4INCL::Particle::theFrozenMomentum

protectedinherited

Definition at line 1088 of file G4INCLParticle.hh.

Referenced by G4INCL::Particle::freezePropagation(), G4INCL::Particle::getFrozenMomentum(), G4INCL::Particle::Particle(), G4INCL::Particle::rotateMomentum(), G4INCL::Particle::setFrozenMomentum(), and G4INCL::Particle::swap().

◆ theHelicity

G4double G4INCL::Particle::theHelicity

privateinherited

Definition at line 1103 of file G4INCLParticle.hh.

Referenced by G4INCL::Particle::getHelicity(), G4INCL::Particle::setHelicity(), and G4INCL::Particle::swap().

◆ theInitialA

G4int G4INCL::Nucleus::theInitialA

private

Definition at line 484 of file G4INCLNucleus.hh.

Referenced by getInitialA().

◆ theInitialS

G4int G4INCL::Nucleus::theInitialS

private

Definition at line 484 of file G4INCLNucleus.hh.

Referenced by getInitialS().

◆ theInitialZ

G4int G4INCL::Nucleus::theInitialZ

private

Definition at line 484 of file G4INCLNucleus.hh.

Referenced by getInitialZ().

◆ theMass

G4double G4INCL::Particle::theMass

privateinherited

Definition at line 1110 of file G4INCLParticle.hh.

Referenced by G4INCL::Particle::adjustEnergyFromMomentum(), G4INCL::Particle::adjustMomentumFromEnergy(), G4INCL::Particle::getINCLMass(), G4INCL::Particle::getKineticEnergy(), G4INCL::Particle::getMass(), G4INCL::Particle::getRealMass(), G4INCL::Particle::getTableMass(), G4INCL::Particle::Particle(), G4INCL::Particle::setMass(), and G4INCL::Particle::swap().

◆ theMomentum

G4INCL::ThreeVector G4INCL::Particle::theMomentum

protectedinherited

Definition at line 1086 of file G4INCLParticle.hh.

◆ theNKaon

G4int G4INCL::Particle::theNKaon

protectedinherited

The number of Kaons inside the nucleus (update during the cascade)

Definition at line 1100 of file G4INCLParticle.hh.

Referenced by emitInsideKaon(), G4INCL::Particle::getNumberOfKaon(), and G4INCL::Particle::setNumberOfKaon().

◆ theNkaonminusInitial

G4int G4INCL::Nucleus::theNkaonminusInitial

private

Definition at line 494 of file G4INCLNucleus.hh.

Referenced by computeSeparationEnergyBalance(), getNumberOfEnteringKaons(), and insertParticle().

◆ theNkaonplusInitial

G4int G4INCL::Nucleus::theNkaonplusInitial

private

The number of entering kaons.

Definition at line 493 of file G4INCLNucleus.hh.

Referenced by computeSeparationEnergyBalance(), getNumberOfEnteringKaons(), and insertParticle().

◆ theNnInitial

G4int G4INCL::Nucleus::theNnInitial

private

The number of entering neutrons.

Definition at line 488 of file G4INCLNucleus.hh.

Referenced by computeSeparationEnergyBalance(), getNumberOfEnteringNeutrons(), and insertParticle().

◆ theNpInitial

G4int G4INCL::Nucleus::theNpInitial

private

The number of entering protons.

Definition at line 486 of file G4INCLNucleus.hh.

Referenced by computeSeparationEnergyBalance(), getNumberOfEnteringProtons(), and insertParticle().

◆ theNpionminusInitial

G4int G4INCL::Nucleus::theNpionminusInitial

private

Definition at line 491 of file G4INCLNucleus.hh.

Referenced by computeSeparationEnergyBalance(), getNumberOfEnteringPions(), and insertParticle().

◆ theNpionplusInitial

G4int G4INCL::Nucleus::theNpionplusInitial

private

The number of entering pions.

Definition at line 490 of file G4INCLNucleus.hh.

Referenced by computeSeparationEnergyBalance(), getNumberOfEnteringPions(), and insertParticle().

◆ theParticipantType

ParticipantType G4INCL::Particle::theParticipantType

protectedinherited

Definition at line 1081 of file G4INCLParticle.hh.

Referenced by G4INCL::Particle::getParticipantType(), G4INCL::Particle::isParticipant(), G4INCL::Particle::isProjectileSpectator(), G4INCL::Particle::isTargetSpectator(), G4INCL::Particle::makeParticipant(), G4INCL::Particle::makeProjectileSpectator(), G4INCL::Particle::makeTargetSpectator(), G4INCL::Particle::Particle(), G4INCL::Particle::setParticipantType(), and G4INCL::Particle::swap().

◆ theParticleBias

G4double G4INCL::Particle::theParticleBias

protectedinherited

Definition at line 1098 of file G4INCLParticle.hh.

Referenced by G4INCL::Particle::getParticleBias(), G4INCL::Particle::setParticleBias(), and G4INCL::Particle::swap().

◆ theParticleSampler

ParticleSampler* G4INCL::Cluster::theParticleSampler

protectedinherited

Definition at line 454 of file G4INCLCluster.hh.

Referenced by G4INCL::Cluster::Cluster(), G4INCL::Cluster::initializeParticles(), Nucleus(), setDensity(), G4INCL::Cluster::swap(), and G4INCL::Cluster::~Cluster().

◆ thePosition

G4INCL::ThreeVector G4INCL::Particle::thePosition

protectedinherited

◆ thePotential

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

private

Pointer to the NuclearPotential object.

Definition at line 531 of file G4INCLNucleus.hh.

Referenced by computeSeparationEnergyBalance(), decayInsideDeltas(), getPotential(), getSurfaceRadius(), Nucleus(), and updatePotentialEnergy().

◆ thePotentialEnergy

G4double G4INCL::Particle::thePotentialEnergy

protectedinherited

Definition at line 1092 of file G4INCLParticle.hh.

Referenced by G4INCL::Cluster::addParticle(), G4INCL::Particle::getPotentialEnergy(), G4INCL::ProjectileRemnant::reset(), G4INCL::Particle::setPotentialEnergy(), G4INCL::Particle::swap(), and G4INCL::Cluster::updateClusterParameters().

◆ theProjectileRemnant

ProjectileRemnant* G4INCL::Nucleus::theProjectileRemnant

private

Pointer to the quasi-projectile.

Owned by the Nucleus object.

Definition at line 525 of file G4INCLNucleus.hh.

Referenced by computeRecoilKinematics(), deleteProjectileRemnant(), fillEventInfo(), finalizeProjectileRemnant(), getConservationBalance(), getProjectileRemnant(), initializeParticles(), and setProjectileRemnant().

◆ thePropagationEnergy

G4double* G4INCL::Particle::thePropagationEnergy

protectedinherited

Definition at line 1084 of file G4INCLParticle.hh.

Referenced by G4INCL::Particle::freezePropagation(), G4INCL::Particle::Particle(), G4INCL::Particle::swap(), and G4INCL::Particle::thawPropagation().

◆ thePropagationMomentum

G4INCL::ThreeVector* G4INCL::Particle::thePropagationMomentum

protectedinherited

Definition at line 1087 of file G4INCLParticle.hh.

Referenced by G4INCL::Particle::freezePropagation(), G4INCL::Particle::getCosRPAngle(), G4INCL::Particle::getLongitudinalPosition(), G4INCL::Particle::getPropagationVelocity(), G4INCL::Particle::Particle(), G4INCL::Particle::swap(), and G4INCL::Particle::thawPropagation().

◆ theS

G4int G4INCL::Particle::theS

protectedinherited

Definition at line 1080 of file G4INCLParticle.hh.

Referenced by G4INCL::ProjectileRemnant::addAllDynamicalSpectators(), G4INCL::ProjectileRemnant::addDynamicalSpectator(), G4INCL::ProjectileRemnant::addMostDynamicalSpectators(), G4INCL::Cluster::addParticle(), applyFinalState(), G4INCL::Cluster::Cluster(), computeOneNucleonRecoilKinematics(), computeRecoilKinematics(), decayMe(), emitInsideKaon(), emitInsideLambda(), emitInsidePions(), emitInsideStrangeParticles(), G4INCL::Particle::getEmissionQValueCorrection(), G4INCL::Particle::getINCLMass(), G4INCL::Particle::getRealMass(), G4INCL::Particle::getS(), G4INCL::Cluster::getSpecies(), G4INCL::Particle::getTableMass(), G4INCL::Particle::getTransferQValueCorrection(), insertParticle(), Nucleus(), G4INCL::Cluster::print(), G4INCL::ProjectileRemnant::removeParticle(), G4INCL::Cluster::setS(), G4INCL::Particle::setType(), G4INCL::Particle::swap(), and G4INCL::Cluster::updateClusterParameters().

◆ theSpin

ThreeVector G4INCL::Cluster::theSpin

protectedinherited

Definition at line 453 of file G4INCLCluster.hh.

Referenced by computeRecoilKinematics(), G4INCL::Cluster::getSpin(), G4INCL::Cluster::setSpin(), and G4INCL::Cluster::swap().

◆ theStore

Store* G4INCL::Nucleus::theStore

private

◆ theType

G4INCL::ParticleType G4INCL::Particle::theType

protectedinherited

◆ theUniverseRadius

G4double G4INCL::Nucleus::theUniverseRadius

private

The radius of the universe.

Definition at line 513 of file G4INCLNucleus.hh.

Referenced by getUniverseRadius(), Nucleus(), and setUniverseRadius().

◆ theZ

G4int G4INCL::Particle::theZ

protectedinherited

Definition at line 1080 of file G4INCLParticle.hh.

◆ tryCN

G4bool G4INCL::Nucleus::tryCN

private

Definition at line 503 of file G4INCLNucleus.hh.

Referenced by applyFinalState(), fillEventInfo(), and getTryCompoundNucleus().

◆ uncorrelatedMomentum

G4double G4INCL::Particle::uncorrelatedMomentum

protectedinherited

Definition at line 1096 of file G4INCLParticle.hh.

Referenced by G4INCL::Particle::getReflectionMomentum(), G4INCL::Particle::setUncorrelatedMomentum(), and G4INCL::Particle::swap().

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

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

Data Structures

Public Member Functions

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ Nucleus() [1/2]

◆ ~Nucleus()

◆ Nucleus() [2/2]

Member Function Documentation

◆ addParticle()

◆ addParticles()

◆ adjustEnergyFromMomentum()

◆ adjustMomentumFromEnergy()

◆ applyFinalState()

◆ boost()

◆ boostVector()

◆ clearParticles()

◆ computeCenterOfMass()

◆ computeDynamicalPotential()

◆ computeExcitationEnergy()

◆ computeOneNucleonRecoilKinematics()

◆ computeRecoilKinematics()

◆ computeSeparationEnergyBalance()

◆ computeTotalEnergy()

◆ containsAntiKaon()

◆ containsDeltas()

◆ containsEtas()

◆ containsKaon()

◆ containsLambda()

◆ containsOmegas()

◆ containsSigma()

◆ decayInsideDeltas()

◆ decayInsideStrangeParticles()

◆ decayMe()

◆ decayOutgoingClusters()

◆ decayOutgoingDeltas()

◆ decayOutgoingNeutralKaon()

◆ decayOutgoingPionResonances()

◆ decayOutgoingSigmaZero()

◆ deleteParticles()

◆ deleteProjectileRemnant()

◆ dump()

◆ emitInsideKaon()

◆ emitInsideLambda()

◆ emitInsidePions()

◆ emitInsideStrangeParticles()

◆ fillEventInfo()

◆ FillINCLBiasVector()

◆ finalizeProjectileRemnant()

◆ freezeInternalMotion()

◆ freezePropagation()

◆ getA()

◆ getAngularMomentum()

◆ getBeta()

◆ getBiasCollisionVector()

◆ getBiasFromVector()

◆ getConservationBalance()

◆ getCosRPAngle()

◆ getDensity()

◆ getEmissionQValueCorrection() [1/2]

◆ getEmissionQValueCorrection() [2/2]

◆ getEmissionTime()

◆ getEnergy()

◆ getExcitationEnergy()

◆ getFrozenEnergy()

◆ getFrozenMomentum()

◆ getHelicity()

◆ getID()

◆ getINCLMass()

◆ getIncomingAngularMomentum()

◆ getIncomingMomentum()

◆ getInitialA()

◆ getInitialEnergy()

◆ getInitialInternalEnergy()