#include <G4INCLProjectileRemnant.hh>

Inheritance diagram for G4INCL::ProjectileRemnant:

Public Types
typedef std::map< long, G4double >	EnergyLevelMap

typedef std::vector< G4double >	EnergyLevels

Public Member Functions
ParticleList	addAllDynamicalSpectators (ParticleList const &pL)
	Add back all dynamical spectators to the projectile remnant. More...

ParticleList	addDynamicalSpectators (ParticleList pL)
	Add back dynamical spectators to the projectile remnant. More...

ParticleList	addMostDynamicalSpectators (ParticleList pL)
	Add back dynamical spectators to the projectile remnant. More...

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	boost (const ThreeVector &aBoostVector)
	Boost the cluster with the indicated velocity. More...

ThreeVector	boostVector () const

void	clearEnergyLevels ()
	Clear the stored energy levels. More...

void	clearParticles ()

void	clearStoredComponents ()
	Clear the stored projectile components. More...

G4double	computeExcitationEnergyExcept (const long exceptID) const
	Compute the excitation energy when a nucleon is removed. More...

G4double	computeExcitationEnergyWith (const ParticleList &pL) const
	Compute the excitation energy if some nucleons are put back. More...

void	deleteParticles ()

void	deleteStoredComponents ()
	Clear the stored projectile components and delete the particles. More...

std::string	dump () const

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...

G4double	getCosRPAngle () const
	Get the cosine of the angle between position and momentum. 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 cluster. More...

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

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

EnergyLevels const &	getGroundStateEnergies () const

G4double	getHelicity ()

long	getID () const

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

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	getNumberOfKaon () const
	Number of Kaon inside de nucleus. More...

G4int	getNumberStoredComponents () const
	Get the number of the stored components. 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

G4double	getPotentialEnergy () const
	Get the particle potential energy. 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...

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...

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

G4INCL::ParticleType	getType () const

G4int	getZ () const
	Returns the charge number. 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...

virtual void	initializeParticles ()
	Initialise the NuclearDensity pointer and sample the particles. 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	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	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...

std::string	print () const

	ProjectileRemnant (ParticleSpecies const &species, const G4double kineticEnergy)

void	propagate (G4double step)

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

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

void	removeParticle (Particle *const p, const G4double theProjectileCorrection)
	Remove a nucleon from the projectile remnant. More...

void	reset ()
	Reset the projectile remnant to the state at the beginning of the cascade. 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	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	setMass (G4double mass)

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

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	setPosition (const ThreeVector &position)
	Set the position of the cluster. More...

void	setPotentialEnergy (G4double v)
	Set the particle potential energy. 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	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	setZ (const G4int Z)
	Set the charge number of the cluster. More...

void	storeComponents ()
	Store the projectile components. More...

void	storeEnergyLevels ()
	Store the energy levels. 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...

	~ProjectileRemnant ()

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
G4bool	addDynamicalSpectator (Particle *const p)
	Add back a nucleon to the projectile remnant. More...

G4double	computeDynamicalPotential ()
	Compute the dynamical cluster potential. More...

G4double	computeExcitationEnergy (const EnergyLevels &levels) const
	Compute the excitation energy for a given configuration. More...

EnergyLevels	getPresentEnergyLevelsExcept (const long exceptID) const

EnergyLevels	getPresentEnergyLevelsWith (const ParticleList &pL) const

ParticleList	getStoredComponents () const

ThreeVector const &	getStoredMomentum (Particle const *const p) const
	Return the stored momentum of a given projectile component. More...

ParticleList	shuffleStoredComponents ()
	Shuffle the list of stored projectile components. More...

Private Attributes
G4double	emissionTime

G4bool	outOfWell

std::map< long, Particle * >	storedComponents
	Return the stored energy of a given projectile component. More...

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

EnergyLevels	theGroundStateEnergies
	Ground-state energies for any number of nucleons. More...

G4double	theHelicity

EnergyLevelMap	theInitialEnergyLevels
	Initial energy levels of the projectile. More...

G4double	theMass

Static Private Attributes
static G4ThreadLocal long	nextID = 1

Detailed Description

Definition at line 58 of file G4INCLProjectileRemnant.hh.

Member Typedef Documentation

◆ EnergyLevelMap

typedef std::map<long, G4double> G4INCL::ProjectileRemnant::EnergyLevelMap

Definition at line 63 of file G4INCLProjectileRemnant.hh.

◆ EnergyLevels

typedef std::vector<G4double> G4INCL::ProjectileRemnant::EnergyLevels

Definition at line 62 of file G4INCLProjectileRemnant.hh.

Constructor & Destructor Documentation

◆ ProjectileRemnant()

G4INCL::ProjectileRemnant::ProjectileRemnant	(	ParticleSpecies const &	species,
		const G4double	kineticEnergy
	)

inline

Definition at line 65 of file G4INCLProjectileRemnant.hh.

      : Cluster(species.theZ, species.theA, species.theS) {
 
      // Use the table mass
      setTableMass();
 
      // Set the kinematics
      const G4double projectileMass = getMass();
      const G4double energy = kineticEnergy + projectileMass;
      const G4double momentumZ = std::sqrt(energy*energy - projectileMass*projectileMass);
 
      // Initialise the particles
      initializeParticles();
      internalBoostToCM();
      putParticlesOffShell();
 
      // Store the energy levels of the ProjectileRemnant (used to compute its
      // excitation energy)
      storeEnergyLevels();
 
      // Boost the whole thing
      const ThreeVector aBoostVector = ThreeVector(0.0, 0.0, momentumZ / energy);
      boost(-aBoostVector);
 
      // Freeze the internal motion of the particles
      freezeInternalMotion();
 
      // Set as projectile spectator
      makeProjectileSpectator();
    }

References G4INCL::Cluster::boost(), G4INCL::KinematicsUtils::energy(), G4INCL::Cluster::freezeInternalMotion(), G4INCL::Particle::getMass(), G4INCL::Cluster::initializeParticles(), G4INCL::Cluster::internalBoostToCM(), G4INCL::Cluster::makeProjectileSpectator(), G4INCL::Cluster::putParticlesOffShell(), G4INCL::Particle::setTableMass(), and storeEnergyLevels().

◆ ~ProjectileRemnant()

G4INCL::ProjectileRemnant::~ProjectileRemnant ( )

inline

Definition at line 96 of file G4INCLProjectileRemnant.hh.

                         {
      deleteStoredComponents();
      // The ProjectileRemnant owns its particles
      deleteParticles();
      clearEnergyLevels();
    }

References clearEnergyLevels(), G4INCL::Cluster::deleteParticles(), and deleteStoredComponents().

Member Function Documentation

◆ addAllDynamicalSpectators()

ParticleList G4INCL::ProjectileRemnant::addAllDynamicalSpectators ( ParticleList const & pL )

Add back all dynamical spectators to the projectile remnant.

Return a list of rejected dynamical spectators.

Definition at line 145 of file G4INCLProjectileRemnant.cc.

                                                                                  {
    // Put all the spectators in the projectile
    ThreeVector theNewMomentum = theMomentum;
    G4double theNewEnergy = theEnergy;
    G4int theNewA = theA;
    G4int theNewZ = theZ;
    G4int theNewS = theS;
    for(ParticleIter p=pL.begin(), e=pL.end(); p!=e; ++p) {
// assert((*p)->isNucleonorLambda());
      // Add the initial (off-shell) momentum and energy to the projectile remnant
      theNewMomentum += getStoredMomentum(*p);
      theNewEnergy += (*p)->getEnergy();
      theNewA += (*p)->getA();
      theNewZ += (*p)->getZ();
      theNewS += (*p)->getS();
    }
 
    // Check that the excitation energy of the new projectile remnant is non-negative
    const G4double theNewMass = ParticleTable::getTableMass(theNewA,theNewZ,theNewS);
    const G4double theNewExcitationEnergy = computeExcitationEnergyWith(pL);
    const G4double theNewEffectiveMass = theNewMass + theNewExcitationEnergy;
 
    // If this condition is satisfied, there is no solution. Fall back on the
    // "most" method
    if(theNewEnergy<theNewEffectiveMass) {
      INCL_WARN("Could not add all the dynamical spectators back into the projectile remnant."
           << " Falling back to the \"most\" method." << '\n');
      return addMostDynamicalSpectators(pL);
    }
 
    // Add all the participants to the projectile remnant
    for(ParticleIter p=pL.begin(), e=pL.end(); p!=e; ++p) {
      particles.push_back(*p);
    }
 
    // Rescale the momentum of the projectile remnant so that sqrt(s) has the
    // correct value
    const G4double scalingFactorSquared = (theNewEnergy*theNewEnergy-theNewEffectiveMass*theNewEffectiveMass)/theNewMomentum.mag2();
    const G4double scalingFactor = std::sqrt(scalingFactorSquared);
    INCL_DEBUG("Scaling factor for the projectile-remnant momentum = " << scalingFactor << '\n');
 
    theA = theNewA;
    theZ = theNewZ;
    theS = theNewS;
    theMomentum = theNewMomentum * scalingFactor;
    theEnergy = theNewEnergy;
 
    return ParticleList();
  }

References addMostDynamicalSpectators(), computeExcitationEnergyWith(), getStoredMomentum(), G4INCL::ParticleTable::getTableMass, INCL_DEBUG, INCL_WARN, G4INCL::ThreeVector::mag2(), G4INCL::Cluster::particles, G4INCL::Particle::theA, G4INCL::Particle::theEnergy, G4INCL::Particle::theMomentum, G4INCL::Particle::theS, and G4INCL::Particle::theZ.

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

◆ addDynamicalSpectator()

G4bool G4INCL::ProjectileRemnant::addDynamicalSpectator ( Particle *const p )

private

Add back a nucleon to the projectile remnant.

Try to add a dynamical spectator back to the projectile remnant. Refuse to do so if this leads to a negative projectile excitation energy. Return true on success, false on failure.

Definition at line 296 of file G4INCLProjectileRemnant.cc.

                                                                    {
// assert(p->isNucleon());
 
    // Add the initial (off-shell) momentum and energy to the projectile remnant
    ThreeVector const &oldMomentum = getStoredMomentum(p);
    const ThreeVector theNewMomentum = theMomentum + oldMomentum;
    const G4double oldEnergy = p->getEnergy();
    const G4double theNewEnergy = theEnergy + oldEnergy;
 
    // Check that the excitation energy of the new projectile remnant is non-negative
    const G4double theNewMass = ParticleTable::getTableMass(theA+p->getA(),theZ+p->getZ(),theS+p->getS());
    const G4double theNewInvariantMassSquared = theNewEnergy*theNewEnergy-theNewMomentum.mag2();
 
    if(theNewInvariantMassSquared<0.)
      return false;
 
    const G4double theNewInvariantMass = std::sqrt(theNewInvariantMassSquared);
 
    if(theNewInvariantMass-theNewMass<-1.e-5)
      return false; // negative excitation energy here
 
    // Add the spectator to the projectile remnant
    theA += p->getA();
    theZ += p->getZ();
    theMomentum = theNewMomentum;
    theEnergy = theNewEnergy;
    particles.push_back(p);
    return true;
  }

References G4INCL::Particle::getA(), G4INCL::Particle::getEnergy(), G4INCL::Particle::getS(), getStoredMomentum(), G4INCL::ParticleTable::getTableMass, G4INCL::Particle::getZ(), G4INCL::ThreeVector::mag2(), G4INCL::Cluster::particles, G4INCL::Particle::theA, G4INCL::Particle::theEnergy, G4INCL::Particle::theMomentum, G4INCL::Particle::theS, and G4INCL::Particle::theZ.

Referenced by addDynamicalSpectators().

◆ addDynamicalSpectators()

ParticleList G4INCL::ProjectileRemnant::addDynamicalSpectators ( ParticleList pL )

Add back dynamical spectators to the projectile remnant.

Try to add the dynamical spectators back to the projectile remnant. Refuse to do so if this leads to a negative projectile excitation energy.

Return a list of rejected dynamical spectators.

Definition at line 122 of file G4INCLProjectileRemnant.cc.

                                                                        {
    // Try as hard as possible to add back all the dynamical spectators.
    // Don't add spectators that lead to negative excitation energies, but
    // iterate over the spectators as many times as possible, until
    // absolutely sure that all of them were rejected.
    unsigned int accepted;
    unsigned long loopCounter = 0;
    const unsigned long maxLoopCounter = 10000000;
    do {
      accepted = 0;
      ParticleList toBeAdded = pL;
      for(ParticleIter p=toBeAdded.begin(), e=toBeAdded.end(); p!=e; ++p) {
        G4bool isAccepted = addDynamicalSpectator(*p);
        if(isAccepted) {
          pL.remove(*p);
          accepted++;
        }
      }
      ++loopCounter;
    } while(loopCounter<maxLoopCounter && accepted > 0); /* Loop checking, 10.07.2015, D.Mancusi */
    return pL;
  }

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

◆ addMostDynamicalSpectators()

ParticleList G4INCL::ProjectileRemnant::addMostDynamicalSpectators ( ParticleList pL )

Add back dynamical spectators to the projectile remnant.

Try as hard as possible to add back all the dynamical spectators. Don't add spectators that lead to negative excitation energies. Start by adding all of them, and repeatedly remove the most troublesome one until the excitation energy becomes non-negative.

Return a list of rejected dynamical spectators.

Definition at line 195 of file G4INCLProjectileRemnant.cc.

                                                                            {
    // Try as hard as possible to add back all the dynamical spectators.
    // Don't add spectators that lead to negative excitation energies. Start by
    // adding all of them, and repeatedly remove the most troublesome one until
    // the excitation energy becomes non-negative.
 
    // Put all the spectators in the projectile
    ThreeVector theNewMomentum = theMomentum;
    G4double theNewEnergy = theEnergy;
    G4int theNewA = theA;
    G4int theNewZ = theZ;
    G4int theNewS = theS;
    for(ParticleIter p=pL.begin(), e=pL.end(); p!=e; ++p) {
// assert((*p)->isNucleonorLambda());
      // Add the initial (off-shell) momentum and energy to the projectile remnant
      theNewMomentum += getStoredMomentum(*p);
      theNewEnergy += (*p)->getEnergy();
      theNewA += (*p)->getA();
      theNewZ += (*p)->getZ();
      theNewS += (*p)->getS();
    }
 
    // Check that the excitation energy of the new projectile remnant is non-negative
    const G4double theNewMass = ParticleTable::getTableMass(theNewA,theNewZ,theNewS);
    const G4double theNewInvariantMassSquared = theNewEnergy*theNewEnergy-theNewMomentum.mag2();
 
    G4bool positiveExcitationEnergy = false;
    if(theNewInvariantMassSquared>=0.) {
      const G4double theNewInvariantMass = std::sqrt(theNewInvariantMassSquared);
      positiveExcitationEnergy = (theNewInvariantMass-theNewMass>-1.e-5);
    }
 
    // Keep removing nucleons from the projectile remnant until we achieve a
    // non-negative excitation energy.
    ParticleList rejected;
    while(!positiveExcitationEnergy && !pL.empty()) { /* Loop checking, 10.07.2015, D.Mancusi */
      G4double maxExcitationEnergy = -1.E30;
      ParticleMutableIter best = pL.end();
      ThreeVector bestMomentum;
      G4double bestEnergy = -1.;
      G4int bestA = -1, bestZ = -1, bestS = 0;
      for(ParticleList::iterator p=pL.begin(), e=pL.end(); p!=e; ++p) {
        // Subtract the initial (off-shell) momentum and energy from the new
        // projectile remnant
        const ThreeVector theNewerMomentum = theNewMomentum - getStoredMomentum(*p);
        const G4double theNewerEnergy = theNewEnergy - (*p)->getEnergy();
        const G4int theNewerA = theNewA - (*p)->getA();
        const G4int theNewerZ = theNewZ - (*p)->getZ();
        const G4int theNewerS = theNewS - (*p)->getS();
 
        const G4double theNewerMass = ParticleTable::getTableMass(theNewerA,theNewerZ,theNewerS);
        const G4double theNewerInvariantMassSquared = theNewerEnergy*theNewerEnergy-theNewerMomentum.mag2();
 
        if(theNewerInvariantMassSquared>=-1.e-5) {
          const G4double theNewerInvariantMass = std::sqrt(std::max(0.,theNewerInvariantMassSquared));
          const G4double theNewerExcitationEnergy = ((theNewerA>1) ? theNewerInvariantMass-theNewerMass : 0.);
          // Pick the nucleon that maximises the excitation energy of the
          // ProjectileRemnant
          if(theNewerExcitationEnergy>maxExcitationEnergy) {
            best = p;
            maxExcitationEnergy = theNewerExcitationEnergy;
            bestMomentum = theNewerMomentum;
            bestEnergy = theNewerEnergy;
            bestA = theNewerA;
            bestZ = theNewerZ;
            bestS = theNewerS;
          }
        }
      }
 
      // If we couldn't even calculate the excitation energy, fail miserably
      if(best==pL.end())
        return pL;
 
      rejected.push_back(*best);
      pL.erase(best);
      theNewMomentum = bestMomentum;
      theNewEnergy = bestEnergy;
      theNewA = bestA;
      theNewZ = bestZ;
      theNewS = bestS;
 
      if(maxExcitationEnergy>0.) {
        // Stop here
        positiveExcitationEnergy = true;
      }
    }
 
    // Add the accepted participants to the projectile remnant
    for(ParticleIter p=pL.begin(), e=pL.end(); p!=e; ++p) {
      particles.push_back(*p);
    }
    theA = theNewA;
    theZ = theNewZ;
    theS = theNewS;
    theMomentum = theNewMomentum;
    theEnergy = theNewEnergy;
 
    return rejected;
  }

References getStoredMomentum(), G4INCL::ParticleTable::getTableMass, G4INCL::ThreeVector::mag2(), G4INCL::Math::max(), G4INCL::Cluster::particles, G4INCL::Particle::theA, G4INCL::Particle::theEnergy, G4INCL::Particle::theMomentum, G4INCL::Particle::theS, and G4INCL::Particle::theZ.

Referenced by addAllDynamicalSpectators().

◆ 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 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(), G4INCL::Nucleus::emitInsideKaon(), G4INCL::Nucleus::emitInsideLambda(), G4INCL::Nucleus::emitInsidePions(), G4INCL::Nucleus::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().

◆ 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 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().

◆ clearEnergyLevels()

void G4INCL::ProjectileRemnant::clearEnergyLevels ( )

inline

Clear the stored energy levels.

Definition at line 153 of file G4INCLProjectileRemnant.hh.

                             {
      theInitialEnergyLevels.clear();
      theGroundStateEnergies.clear();
    }

References theGroundStateEnergies, and theInitialEnergyLevels.

Referenced by ~ProjectileRemnant().

◆ 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().

◆ clearStoredComponents()

void G4INCL::ProjectileRemnant::clearStoredComponents ( )

inline

Clear the stored projectile components.

Definition at line 148 of file G4INCLProjectileRemnant.hh.

                                 {
      storedComponents.clear();
    }

References storedComponents.

Referenced by deleteStoredComponents().

◆ 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::ProjectileRemnant::computeExcitationEnergy ( const EnergyLevels & levels ) const

private

Compute the excitation energy for a given configuration.

The function that does the real job of calculating the excitation energy for a given configuration of energy levels.

Parameters

levels a configuration of energy levels

Returns: the excitation energy

Definition at line 336 of file G4INCLProjectileRemnant.cc.

                                                                                      {
    // The ground-state energy is the sum of the A smallest initial projectile
    // energies.
    // For the last nucleon, return 0 so that the algorithm will just put it on
    // shell.
    const unsigned theNewA = levels.size();
// assert(theNewA>0);
    if(theNewA==1)
      return 0.;
 
    const G4double groundState = theGroundStateEnergies.at(theNewA-1);
 
    // Compute the sum of the presently occupied energy levels
    const G4double excitedState = std::accumulate(
        levels.begin(),
        levels.end(),
        0.);
 
    return excitedState-groundState;
  }

References theGroundStateEnergies.

Referenced by computeExcitationEnergyExcept(), and computeExcitationEnergyWith().

◆ computeExcitationEnergyExcept()

G4double G4INCL::ProjectileRemnant::computeExcitationEnergyExcept ( const long exceptID ) const

Compute the excitation energy when a nucleon is removed.

Compute the excitation energy of the projectile-like remnant as the difference between the initial and the present configuration. This follows the algorithm proposed by A. Boudard in INCL4.2-HI, as implemented in Geant4.

Returns: the excitation energy

Definition at line 326 of file G4INCLProjectileRemnant.cc.

                                                                                     {
    const EnergyLevels theEnergyLevels = getPresentEnergyLevelsExcept(exceptID);
    return computeExcitationEnergy(theEnergyLevels);
  }

References computeExcitationEnergy(), and getPresentEnergyLevelsExcept().

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

◆ computeExcitationEnergyWith()

G4double G4INCL::ProjectileRemnant::computeExcitationEnergyWith ( const ParticleList & pL ) const

Compute the excitation energy if some nucleons are put back.

Returns: the excitation energy

Definition at line 331 of file G4INCLProjectileRemnant.cc.

                                                                                      {
    const EnergyLevels theEnergyLevels = getPresentEnergyLevelsWith(pL);
    return computeExcitationEnergy(theEnergyLevels);
  }

References computeExcitationEnergy(), and getPresentEnergyLevelsWith().

Referenced by addAllDynamicalSpectators().

◆ 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(), G4INCL::Nucleus::decayOutgoingClusters(), reset(), and ~ProjectileRemnant().

◆ deleteStoredComponents()

void G4INCL::ProjectileRemnant::deleteStoredComponents ( )

inline

Clear the stored projectile components and delete the particles.

Definition at line 141 of file G4INCLProjectileRemnant.hh.

                                  {
      for(std::map<long,Particle*>::const_iterator p=storedComponents.begin(), e=storedComponents.end(); p!=e; ++p)
        delete p->second;
      clearStoredComponents();
    }

References clearStoredComponents(), and storedComponents.

Referenced by ~ProjectileRemnant().

◆ 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().

◆ 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().

◆ 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 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 addDynamicalSpectator(), G4INCL::Cluster::addParticle(), G4INCL::ClusteringModelIntercomparison::clusterCanEscape(), G4INCL::Nucleus::computeOneNucleonRecoilKinematics(), G4INCL::INCL::continueCascade(), G4INCL::ClusterDecay::decay(), G4INCL::Nucleus::emitInsideLambda(), G4INCL::Nucleus::emitInsideStrangeParticles(), G4INCL::Nucleus::fillEventInfo(), G4INCL::ParticleEntryChannel::fillFinalState(), G4INCL::TransmissionChannel::fillFinalState(), G4INCL::Nucleus::finalizeProjectileRemnant(), G4INCL::ClusteringModelIntercomparison::findClusterStartingFrom(), G4INCL::PauliStandard::getBlockingProbability(), G4INCL::ClusteringModelIntercomparison::getCluster(), G4INCL::Nucleus::getConservationBalance(), G4INCL::SurfaceAvatar::getTransmissionProbability(), G4INCL::NuclearDensity::getTransmissionRadius(), G4INCL::TransmissionChannel::initializeKineticEnergyOutside(), G4INCL::Nucleus::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(), 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 G4INCL::Nucleus::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(), G4INCL::Nucleus::decayMe(), G4INCL::Nucleus::decayOutgoingClusters(), G4INCL::Nucleus::decayOutgoingPionResonances(), G4INCL::Nucleus::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().

◆ 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().

◆ 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 G4INCL::Nucleus::emitInsideKaon(), G4INCL::Nucleus::emitInsideLambda(), G4INCL::Nucleus::emitInsidePions(), G4INCL::Nucleus::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 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(), removeParticle(), G4INCL::StandardPropagationModel::shootComposite(), G4INCL::StandardPropagationModel::shootParticle(), G4INCL::KinematicsUtils::squareTotalEnergyInCM(), G4INCL::KinematicsUtils::transformToLocalEnergyFrame(), and G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::twoBodyDecay().

◆ getExcitationEnergy()

G4double G4INCL::Cluster::getExcitationEnergy ( ) const

inlineinherited

Get the excitation energy of the cluster.

Definition at line 155 of file G4INCLCluster.hh.

155{ return theExcitationEnergy; }

G4INCL::Cluster::theExcitationEnergy

G4double theExcitationEnergy

Definition: G4INCLCluster.hh:452

References G4INCL::Cluster::theExcitationEnergy.

Referenced by G4INCL::Nucleus::fillEventInfo(), G4INCL::Nucleus::getConservationBalance(), and G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::recursiveDecay().

◆ 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.

◆ getGroundStateEnergies()

EnergyLevels const & G4INCL::ProjectileRemnant::getGroundStateEnergies ( ) const

inline

Definition at line 207 of file G4INCLProjectileRemnant.hh.

                                                       {
      return theGroundStateEnergies;
    }

References theGroundStateEnergies.

◆ 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(), getStoredMomentum(), G4INCL::ParticleEntryChannel::particleEnters(), and 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().

◆ 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 G4INCL::Nucleus::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(), G4INCL::Nucleus::emitInsideKaon(), G4INCL::Nucleus::emitInsideLambda(), G4INCL::Nucleus::emitInsidePions(), G4INCL::Nucleus::emitInsideStrangeParticles(), G4INCL::Nucleus::fillEventInfo(), G4INCL::ParticleEntryChannel::fillFinalState(), G4INCL::SurfaceAvatar::getChannel(), G4INCL::Nucleus::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.

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

G4INCL::Particle::theNKaon

G4int theNKaon

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

Definition: G4INCLParticle.hh:1100

References G4INCL::Particle::theNKaon.

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

◆ getNumberStoredComponents()

G4int G4INCL::ProjectileRemnant::getNumberStoredComponents ( ) const

inline

Get the number of the stored components.

Definition at line 184 of file G4INCLProjectileRemnant.hh.

                                            {
      return storedComponents.size();
    }

References storedComponents.

◆ 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 G4INCL::Nucleus::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.

◆ 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(), G4INCL::Nucleus::emitInsideKaon(), G4INCL::Nucleus::emitInsideLambda(), G4INCL::Nucleus::emitInsidePions(), G4INCL::Nucleus::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().

◆ getPresentEnergyLevelsExcept()

ProjectileRemnant::EnergyLevels G4INCL::ProjectileRemnant::getPresentEnergyLevelsExcept ( const long exceptID ) const

private

Definition at line 357 of file G4INCLProjectileRemnant.cc.

                                                                                                         {
    EnergyLevels theEnergyLevels;
    for(ParticleIter p=particles.begin(), e=particles.end(); p!=e; ++p) {
      if((*p)->getID()!=exceptID) {
        EnergyLevelMap::const_iterator i = theInitialEnergyLevels.find((*p)->getID());
// assert(i!=theInitialEnergyLevels.end());
        theEnergyLevels.push_back(i->second);
      }
    }
// assert(theEnergyLevels.size()==particles.size()-1);
    return theEnergyLevels;
  }

References G4INCL::Cluster::particles, and theInitialEnergyLevels.

Referenced by computeExcitationEnergyExcept().

◆ getPresentEnergyLevelsWith()

ProjectileRemnant::EnergyLevels G4INCL::ProjectileRemnant::getPresentEnergyLevelsWith ( const ParticleList & pL ) const

private

Definition at line 370 of file G4INCLProjectileRemnant.cc.

                                                                                                          {
    EnergyLevels theEnergyLevels;
    for(ParticleIter p=particles.begin(), e=particles.end(); p!=e; ++p) {
      EnergyLevelMap::const_iterator i = theInitialEnergyLevels.find((*p)->getID());
// assert(i!=theInitialEnergyLevels.end());
      theEnergyLevels.push_back(i->second);
    }
    for(ParticleIter p=pL.begin(), e=pL.end(); p!=e; ++p) {
      EnergyLevelMap::const_iterator i = theInitialEnergyLevels.find((*p)->getID());
// assert(i!=theInitialEnergyLevels.end());
      theEnergyLevels.push_back(i->second);
    }
 
// assert(theEnergyLevels.size()==particles.size()+pL.size());
    return theEnergyLevels;
  }

References G4INCL::Cluster::particles, and theInitialEnergyLevels.

Referenced by computeExcitationEnergyWith().

◆ 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 G4INCL::Nucleus::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 addDynamicalSpectator(), G4INCL::Cluster::addParticle(), G4INCL::ClusteringModelIntercomparison::clusterCanEscape(), G4INCL::Nucleus::computeOneNucleonRecoilKinematics(), G4INCL::ClusterDecay::decay(), G4INCL::Nucleus::emitInsideKaon(), G4INCL::Nucleus::emitInsideLambda(), G4INCL::Nucleus::emitInsideStrangeParticles(), G4INCL::Nucleus::fillEventInfo(), G4INCL::ParticleEntryChannel::fillFinalState(), G4INCL::TransmissionChannel::fillFinalState(), G4INCL::Nucleus::finalizeProjectileRemnant(), G4INCL::Nucleus::getConservationBalance(), G4INCL::SurfaceAvatar::getTransmissionProbability(), G4INCL::TransmissionChannel::initializeKineticEnergyOutside(), G4INCL::Nucleus::insertParticle(), G4INCL::ClusterDecay::isStable(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::phaseSpaceDecay(), G4INCL::INCL::preCascade(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::recursiveDecay(), 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; }

G4INCL::Cluster::theSpin

ThreeVector theSpin

Definition: G4INCLCluster.hh:453

References G4INCL::Cluster::theSpin.

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

◆ getStoredComponents()

ParticleList G4INCL::ProjectileRemnant::getStoredComponents ( ) const

inlineprivate

Definition at line 234 of file G4INCLProjectileRemnant.hh.

                                             {
      ParticleList pL;
      for(std::map<long,Particle*>::const_iterator p=storedComponents.begin(), e=storedComponents.end(); p!=e; ++p)
        pL.push_back(p->second);
      return pL;
    }

References storedComponents.

Referenced by shuffleStoredComponents().

◆ getStoredMomentum()

ThreeVector const & G4INCL::ProjectileRemnant::getStoredMomentum ( Particle const *const p ) const

inlineprivate

Return the stored momentum of a given projectile component.

Definition at line 242 of file G4INCLProjectileRemnant.hh.

                                                                         {
      std::map<long,Particle*>::const_iterator i = storedComponents.find(p->getID());
      if(i==storedComponents.end()) {
        INCL_ERROR("Couldn't find particle " << p->getID() << " in the list of projectile components" << '\n');
        return p->getMomentum();
      } else {
        return i->second->getMomentum();
      }
    }

References G4INCL::Particle::getID(), G4INCL::Particle::getMomentum(), INCL_ERROR, and storedComponents.

Referenced by addAllDynamicalSpectators(), addDynamicalSpectator(), and addMostDynamicalSpectators().

◆ 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 G4INCL::Nucleus::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 G4INCL::Nucleus::emitInsideKaon(), G4INCL::Nucleus::emitInsideLambda(), G4INCL::Nucleus::emitInsidePions(), G4INCL::Nucleus::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.

◆ 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().

◆ 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.

◆ 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 addDynamicalSpectator(), G4INCL::Cluster::addParticle(), G4INCL::CoulombNonRelativistic::bringToSurface(), G4INCL::Nucleus::computeOneNucleonRecoilKinematics(), G4INCL::ClusterDecay::decay(), G4INCL::CoulombNonRelativistic::distortOut(), G4INCL::Nucleus::emitInsideKaon(), G4INCL::Nucleus::emitInsidePions(), G4INCL::Nucleus::emitInsideStrangeParticles(), G4INCL::Nucleus::fillEventInfo(), G4INCL::ParticleEntryChannel::fillFinalState(), G4INCL::TransmissionChannel::fillFinalState(), G4INCL::Nucleus::finalizeProjectileRemnant(), G4INCL::ClusteringModelIntercomparison::findClusterStartingFrom(), G4INCL::PauliStandard::getBlockingProbability(), G4INCL::ClusteringModelIntercomparison::getCluster(), G4INCL::Nucleus::getConservationBalance(), G4INCL::Nucleus::getTransmissionBarrier(), G4INCL::SurfaceAvatar::getTransmissionProbability(), G4INCL::NuclearDensity::getTransmissionRadius(), G4INCL::TransmissionChannel::initializeKineticEnergyOutside(), G4INCL::Nucleus::insertParticle(), G4INCL::ClusterDecay::isStable(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::phaseSpaceDecay(), G4INCL::INCL::preCascade(), G4INCL::ClusterDecay::anonymous_namespace{G4INCLClusterDecay.cc}::recursiveDecay(), 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().

◆ 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::Cluster::initializeParticles ( )

virtualinherited

Initialise the NuclearDensity pointer and sample the particles.

Reimplemented in G4INCL::Nucleus.

Definition at line 43 of file G4INCLCluster.cc.

                                    {
// assert(theA>=2);
    const ThreeVector oldPosition = thePosition;
    theParticleSampler->sampleParticlesIntoList(thePosition, particles);
#if !defined(NDEBUG) && !defined(INCLXX_IN_GEANT4_MODE)
    const G4int theMassNumber = theA;
    const G4int theChargeNumber = theZ;
#endif
    updateClusterParameters();
    thePosition = oldPosition;
// assert(theMassNumber==theA && theChargeNumber==theZ);
    INCL_DEBUG("Cluster initialized:" << '\n' << print());
  }

References INCL_DEBUG, G4INCL::Cluster::particles, G4INCL::Cluster::print(), G4INCL::ParticleSampler::sampleParticlesIntoList(), G4INCL::Particle::theA, G4INCL::Cluster::theParticleSampler, G4INCL::Particle::thePosition, G4INCL::Particle::theZ, and G4INCL::Cluster::updateClusterParameters().

Referenced by G4INCL::Nucleus::initializeParticles(), and ProjectileRemnant().

◆ 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 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().

◆ 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(), G4INCL::Nucleus::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(), G4INCL::Nucleus::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().

◆ 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(), G4INCL::Nucleus::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(), G4INCL::Nucleus::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 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().

◆ print()

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(), removeParticle(), 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().

◆ 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 ProjectileRemnant().

◆ removeParticle() [1/2]

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 removeParticle().

◆ removeParticle() [2/2]

void G4INCL::ProjectileRemnant::removeParticle	(	Particle *const	p,
		const G4double	theProjectileCorrection
	)

Remove a nucleon from the projectile remnant.

Parameters

p	particle to be removed
theProjectileCorrection	correction to be given to the projectile total energy

Definition at line 76 of file G4INCLProjectileRemnant.cc.

                                                                                                   {
// assert(p->isNucleon() || p->isLambda());
 
    INCL_DEBUG("The following Particle is about to be removed from the ProjectileRemnant:"
        << '\n' << p->print()
        << "theProjectileCorrection=" << theProjectileCorrection << '\n');
    // Update A, Z, S, momentum, and energy of the projectile remnant
    theA -= p->getA();
    theZ -= p->getZ();
    theS -= p->getS();
 
    ThreeVector const &oldMomentum = p->getMomentum();
    const G4double oldEnergy = p->getEnergy();
    Cluster::removeParticle(p);
 
#if !defined(NDEBUG) && !defined(INCLXX_IN_GEANT4_MODE)
    ThreeVector theTotalMomentum;
    G4double theTotalEnergy = 0.;
    const G4double theThreshold = 0.1;
#endif
 
    if(getA()>0) { // if there are any particles left
// assert((unsigned int)getA()==particles.size());
 
      const G4double theProjectileCorrectionPerNucleon = theProjectileCorrection / particles.size();
 
      // Update the kinematics of the components
      for(ParticleIter i=particles.begin(), e=particles.end(); i!=e; ++i) {
        (*i)->setEnergy((*i)->getEnergy() + theProjectileCorrectionPerNucleon);
        (*i)->setMass((*i)->getInvariantMass());
#if !defined(NDEBUG) && !defined(INCLXX_IN_GEANT4_MODE)
        theTotalMomentum += (*i)->getMomentum();
        theTotalEnergy += (*i)->getEnergy();
#endif
      }
    }
 
    theMomentum -= oldMomentum;
    theEnergy -= oldEnergy - theProjectileCorrection;
 
// assert(std::abs((theTotalMomentum-theMomentum).mag())<theThreshold);
// assert(std::abs(theTotalEnergy-theEnergy)<theThreshold);
    INCL_DEBUG("After Particle removal, the ProjectileRemnant looks like this:"
        << '\n' << print());
  }

References G4INCL::Particle::getA(), G4INCL::Particle::getEnergy(), G4INCL::Particle::getMomentum(), G4INCL::Particle::getS(), G4INCL::Particle::getZ(), INCL_DEBUG, G4INCL::Cluster::particles, G4INCL::Cluster::print(), G4INCL::Particle::print(), G4INCL::Cluster::removeParticle(), G4INCL::Particle::theA, G4INCL::Particle::theEnergy, G4INCL::Particle::theMomentum, G4INCL::Particle::theS, and G4INCL::Particle::theZ.

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

◆ reset()

void G4INCL::ProjectileRemnant::reset ( )

Reset the projectile remnant to the state at the beginning of the cascade.

Definition at line 51 of file G4INCLProjectileRemnant.cc.

                                {
    deleteParticles();
    thePosition = ThreeVector();
    theMomentum = ThreeVector();
    theEnergy = 0.0;
    thePotentialEnergy = 0.0;
    theA = 0;
    theZ = 0;
    nCollisions = 0;
 
    for(std::map<long, Particle*>::const_iterator i=storedComponents.begin(); i!=storedComponents.end(); ++i) {
      Particle *p = new Particle(*(i->second));
      EnergyLevelMap::iterator energyIter = theInitialEnergyLevels.find(i->first);
// assert(energyIter!=theInitialEnergyLevels.end());
      const G4double energyLevel = energyIter->second;
      theInitialEnergyLevels.erase(energyIter);
      theInitialEnergyLevels[p->getID()] = energyLevel;
      addParticle(p);
    }
    if(theA>0)
      thePosition /= theA;
    setTableMass();
    INCL_DEBUG("ProjectileRemnant object was reset:" << '\n' << print());
  }

References G4INCL::Cluster::addParticle(), G4INCL::Cluster::deleteParticles(), G4INCL::Particle::getID(), INCL_DEBUG, G4INCL::Particle::nCollisions, G4INCL::Particle::Particle(), G4INCL::Cluster::print(), G4INCL::Particle::setTableMass(), storedComponents, G4INCL::Particle::theA, G4INCL::Particle::theEnergy, theInitialEnergyLevels, G4INCL::Particle::theMomentum, G4INCL::Particle::thePosition, G4INCL::Particle::thePotentialEnergy, and G4INCL::Particle::theZ.

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

◆ 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 G4INCL::Nucleus::decayOutgoingPionResonances(), G4INCL::Nucleus::decayOutgoingSigmaZero(), and G4INCL::SurfaceAvatar::postInteraction().

◆ 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 G4INCL::Nucleus::computeOneNucleonRecoilKinematics(), G4INCL::Nucleus::emitInsideKaon(), G4INCL::Nucleus::emitInsideLambda(), G4INCL::Nucleus::emitInsidePions(), G4INCL::Nucleus::emitInsideStrangeParticles(), G4INCL::Nucleus::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 G4INCL::Nucleus::emitInsideKaon(), G4INCL::Nucleus::emitInsideLambda(), G4INCL::Nucleus::emitInsidePions(), G4INCL::Nucleus::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 G4INCL::Nucleus::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 G4INCL::Nucleus::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()); }

G4INCL::Particle::setMass

void setMass(G4double mass)

Definition: G4INCLParticle.hh:778

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

Referenced by G4INCL::Cluster::Cluster(), G4INCL::ParticleEntryChannel::particleEnters(), G4INCL::Particle::setType(), 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 G4INCL::Nucleus::computeRecoilKinematics(), G4INCL::DeltaProductionChannel::fillFinalState(), G4INCL::NDeltaEtaProductionChannel::fillFinalState(), G4INCL::NDeltaOmegaProductionChannel::fillFinalState(), G4INCL::NDeltaToDeltaLKChannel::fillFinalState(), G4INCL::NDeltaToDeltaSKChannel::fillFinalState(), G4INCL::Nucleus::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 G4INCL::Nucleus::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.

◆ 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().

◆ 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 G4INCL::Nucleus::emitInsideKaon(), G4INCL::Nucleus::emitInsideLambda(), G4INCL::Nucleus::emitInsidePions(), G4INCL::Nucleus::emitInsideStrangeParticles(), G4INCL::Store::loadParticles(), G4INCL::ParticleEntryChannel::particleEnters(), G4INCL::TransmissionChannel::particleLeaves(), and G4INCL::Nucleus::updatePotentialEnergy().

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

S

G4double S(G4double temp)

Definition: G4DNAElectronHoleRecombination.cc:76

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 G4INCL::Nucleus::finalizeProjectileRemnant(), G4INCL::INCL::makeCompoundNucleus(), and G4INCL::Nucleus::useFusionKinematics().

◆ 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 G4INCL::Nucleus::decayOutgoingPionResonances(), G4INCL::Nucleus::decayOutgoingSigmaZero(), G4INCL::Nucleus::emitInsideKaon(), G4INCL::Nucleus::emitInsideLambda(), G4INCL::Nucleus::emitInsidePions(), G4INCL::Nucleus::emitInsideStrangeParticles(), G4INCL::TransmissionChannel::particleLeaves(), ProjectileRemnant(), and 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().

◆ 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().

◆ shuffleStoredComponents()

ParticleList G4INCL::ProjectileRemnant::shuffleStoredComponents ( )

inlineprivate

Shuffle the list of stored projectile components.

Definition at line 228 of file G4INCLProjectileRemnant.hh.

                                           {
      ParticleList pL = getStoredComponents();
      std::shuffle(pL.begin(), pL.end(), Random::getAdapter());
      return pL;
    }

References G4INCL::Random::getAdapter(), and getStoredComponents().

◆ storeComponents()

void G4INCL::ProjectileRemnant::storeComponents ( )

inline

Store the projectile components.

Definition at line 176 of file G4INCLProjectileRemnant.hh.

                           {
      for(ParticleIter p=particles.begin(), e=particles.end(); p!=e; ++p) {
        // Store the particles (needed for forced CN)
        storedComponents[(*p)->getID()]=new Particle(**p);
      }
    }

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

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

◆ storeEnergyLevels()

void G4INCL::ProjectileRemnant::storeEnergyLevels ( )

inline

Store the energy levels.

Definition at line 189 of file G4INCLProjectileRemnant.hh.

                             {
      EnergyLevels energies;
 
      for(ParticleIter p=particles.begin(), e=particles.end(); p!=e; ++p) {
        const G4double theCMEnergy = (*p)->getEnergy();
        // Store the CM energy in the EnergyLevels map
        theInitialEnergyLevels[(*p)->getID()] = theCMEnergy;
        energies.push_back(theCMEnergy);
      }
 
      std::sort(energies.begin(), energies.end());
// assert(energies.size()==(unsigned int)theA);
      theGroundStateEnergies.resize(energies.size());
      // Compute the partial sums of the CM energies -- they are our reference
      // ground-state energies for any number of nucleons
      std::partial_sum(energies.begin(), energies.end(), theGroundStateEnergies.begin());
    }

References G4INCL::Cluster::particles, theGroundStateEnergies, and theInitialEnergyLevels.

Referenced by ProjectileRemnant().

◆ 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().

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().

◆ 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(), 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.

◆ 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().

◆ storedComponents

std::map<long, Particle*> G4INCL::ProjectileRemnant::storedComponents

private

Return the stored energy of a given projectile component.

Stored projectile components

These particles are owned by the ProjectileRemnant.

Definition at line 275 of file G4INCLProjectileRemnant.hh.

Referenced by clearStoredComponents(), deleteStoredComponents(), getNumberStoredComponents(), getStoredComponents(), getStoredMomentum(), reset(), and storeComponents().

◆ theA

G4int G4INCL::Particle::theA

protectedinherited

Definition at line 1080 of file G4INCLParticle.hh.

◆ 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().

◆ theEnergy

G4double G4INCL::Particle::theEnergy

protectedinherited

◆ theExcitationEnergy

G4double G4INCL::Cluster::theExcitationEnergy

protectedinherited

Definition at line 452 of file G4INCLCluster.hh.

Referenced by G4INCL::Nucleus::computeOneNucleonRecoilKinematics(), G4INCL::Nucleus::computeRecoilKinematics(), G4INCL::Cluster::getExcitationEnergy(), G4INCL::Nucleus::getExcitationEnergy(), G4INCL::Cluster::setExcitationEnergy(), G4INCL::Cluster::swap(), and G4INCL::Nucleus::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().

◆ theGroundStateEnergies

EnergyLevels G4INCL::ProjectileRemnant::theGroundStateEnergies

private

Ground-state energies for any number of nucleons.

Definition at line 281 of file G4INCLProjectileRemnant.hh.

Referenced by clearEnergyLevels(), computeExcitationEnergy(), getGroundStateEnergies(), and storeEnergyLevels().

◆ 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().

◆ theInitialEnergyLevels

EnergyLevelMap G4INCL::ProjectileRemnant::theInitialEnergyLevels

private

Initial energy levels of the projectile.

Definition at line 278 of file G4INCLProjectileRemnant.hh.

Referenced by clearEnergyLevels(), getPresentEnergyLevelsExcept(), getPresentEnergyLevelsWith(), reset(), and storeEnergyLevels().

◆ 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 G4INCL::Nucleus::emitInsideKaon(), G4INCL::Particle::getNumberOfKaon(), and G4INCL::Particle::setNumberOfKaon().

◆ 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(), G4INCL::Nucleus::Nucleus(), G4INCL::Nucleus::setDensity(), G4INCL::Cluster::swap(), and G4INCL::Cluster::~Cluster().

◆ thePosition

G4INCL::ThreeVector G4INCL::Particle::thePosition

protectedinherited

Definition at line 1089 of file G4INCLParticle.hh.

◆ thePotentialEnergy

G4double G4INCL::Particle::thePotentialEnergy

protectedinherited

Definition at line 1092 of file G4INCLParticle.hh.

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

◆ 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.

◆ theSpin

ThreeVector G4INCL::Cluster::theSpin

protectedinherited

Definition at line 453 of file G4INCLCluster.hh.

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

◆ theType

G4INCL::ParticleType G4INCL::Particle::theType

protectedinherited

◆ theZ

G4int G4INCL::Particle::theZ

protectedinherited

Definition at line 1080 of file G4INCLParticle.hh.

◆ 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/G4INCLProjectileRemnant.hh
source/processes/hadronic/models/inclxx/incl_physics/src/G4INCLProjectileRemnant.cc

Public Types

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

Member Typedef Documentation

◆ EnergyLevelMap

◆ EnergyLevels

Constructor & Destructor Documentation

◆ ProjectileRemnant()

◆ ~ProjectileRemnant()

Member Function Documentation

◆ addAllDynamicalSpectators()

◆ addDynamicalSpectator()

◆ addDynamicalSpectators()

◆ addMostDynamicalSpectators()

◆ addParticle()

◆ addParticles()

◆ adjustEnergyFromMomentum()

◆ adjustMomentumFromEnergy()

◆ boost()

◆ boostVector()

◆ clearEnergyLevels()

◆ clearParticles()

◆ clearStoredComponents()

◆ computeDynamicalPotential()

◆ computeExcitationEnergy()

◆ computeExcitationEnergyExcept()

◆ computeExcitationEnergyWith()

◆ deleteParticles()

◆ deleteStoredComponents()

◆ dump()

◆ FillINCLBiasVector()

◆ freezeInternalMotion()

◆ freezePropagation()

◆ getA()

◆ getAngularMomentum()

◆ getBeta()

◆ getBiasCollisionVector()

◆ getBiasFromVector()

◆ getCosRPAngle()

◆ getEmissionQValueCorrection() [1/2]

◆ getEmissionQValueCorrection() [2/2]

◆ getEmissionTime()

◆ getEnergy()

◆ getExcitationEnergy()

◆ getFrozenEnergy()

◆ getFrozenMomentum()

◆ getGroundStateEnergies()

◆ getHelicity()

◆ getID()

◆ getINCLMass()

◆ getInvariantMass()

◆ getKineticEnergy()

◆ getLongitudinalPosition()

◆ getMass()

◆ getMomentum()

◆ getNumberOfCollisions()

◆ getNumberOfDecays()

◆ getNumberOfKaon()

◆ getNumberStoredComponents()

◆ getParticipantType()

◆ getParticleBias()

◆ getParticleList()

◆ getParticles()

◆ getPosition()

◆ getPotentialEnergy()

◆ getPresentEnergyLevelsExcept()

◆ getPresentEnergyLevelsWith()

◆ getPropagationVelocity()

◆ getRealMass()

◆ getReflectionMomentum()

◆ getS()

◆ getSpecies()

◆ getSpin()