33 #define INCLXX_IN_GEANT4_MODE 1
44 #ifndef G4INCLNUCLEUS_HH_
45 #define G4INCLNUCLEUS_HH_
131 for(
ParticleIter i=outgoing.begin(), e=outgoing.end(); i!=e; ++i) {
245 for(
ParticleIter i=inside.begin(), e=inside.end(); i!=e; ++i)
246 if((*i)->isDelta())
return true;
G4int getNumberOfEnteringNeutrons() const
G4int getA() const
Returns the baryon number.
ConservationBalance getConservationBalance(EventInfo const &theEventInfo, const G4bool afterRecoil) const
Compute charge, mass, energy and momentum balance.
void initializeParticles()
Call the Cluster method to generate the initial distribution of particles.
G4bool containsDeltas()
Returns true if the nucleus contains any deltas.
The INCL configuration object.
void updatePotentialEnergy(Particle *p) const
Update the particle potential energy.
ThreeVector incomingMomentum
G4double getReflectionMomentum() const
Return the reflection momentum.
void setIncomingAngularMomentum(const ThreeVector &j)
Set the incoming angular-momentum vector.
void setIncomingMomentum(const ThreeVector &p)
Set the incoming momentum vector.
const G4double eSquared
Coulomb conversion factor [MeV*fm].
ThreeVector incomingAngularMomentum
G4int getInitialZ() const
ParticleList const & getParticles() const
Return the list of "active" particles (i.e.
std::string print()
Print the nucleus info.
G4int getProjectileChargeNumber() const
Return the charge number of the projectile.
const ThreeVector & getIncomingMomentum() const
Get the incoming momentum vector.
G4double getFermiMomentum(const Particle *const p) const
Return the Fermi momentum for a particle.
G4bool decayOutgoingClusters()
Force the decay of unstable outgoing clusters.
void setInitialEnergy(const G4double e)
Set the initial energy.
G4bool isNucleusNucleusCollision() const
Is it a nucleus-nucleus collision?
G4bool isNucleusNucleus
true if running a nucleus-nucleus collision
void computeOneNucleonRecoilKinematics()
Compute the recoil kinematics for a 1-nucleon remnant.
G4int getInitialA() const
void applyFinalState(FinalState *)
Apply reaction final state information to the nucleus.
void setDensity(NuclearDensity const *const d)
Setter for theDensity.
G4bool isTargetSpectator() const
G4int projectileA
The mass number of the projectile.
const ThreeVector & getIncomingAngularMomentum() const
Get the incoming angular-momentum vector.
G4double initialInternalEnergy
Class for constructing a projectile-like remnant.
void setParticleNucleusCollision()
Set a particle-nucleus collision.
G4int getProjectileMassNumber() const
Return the mass number of the projectile.
Struct for conservation laws.
G4double getSurfaceRadius(Particle const *const particle) const
Get the maximum allowed radius for a given particle.
G4int projectileZ
The charge number of the projectile.
G4bool hasRemnant() const
Does the nucleus give a cascade remnant?
Nucleus & operator=(const Nucleus &rhs)
Dummy assignment operator to silence Coverity warning.
Cluster is a particle (inherits from the Particle class) that is actually a collection of elementary ...
G4double getInitialInternalEnergy() const
void useFusionKinematics()
Adjust the kinematics for complete-fusion events.
G4double getInitialEnergy() const
Get the initial energy.
Final state of an interaction.
G4double getSeparationEnergy(const Particle *const p) const
Return the separation energy for a particle.
G4bool decayOutgoingDeltas()
Force the decay of outgoing deltas.
ParticleList const & getOutgoingParticles() const
Return the list of outgoing particles (i.e.
Simple container for output of event results.
void setDensity(NuclearDensity const *const d)
Setter for theDensity.
Book & getBook()
Return the pointer to the Book object which keeps track of various counters.
G4int getZ() const
Returns the charge number.
NuclearDensity const * theDensity
Pointer to the NuclearDensity object.
void setPotentialEnergy(G4double v)
Set the particle potential energy.
G4double computeSeparationEnergyBalance() const
Outgoing - incoming separation energies.
void setUniverseRadius(const G4double universeRadius)
Setter for theUniverseRadius.
G4double getExcitationEnergy() const
Get the excitation energy of the nucleus.
ParticleList const & getUpdatedParticles() const
Get the list of particles that were updated by the last applied final state.
G4double computeExcitationEnergy() const
Compute the current excitation energy.
G4int getNumberOfEnteringProtons() const
void deleteProjectileRemnant()
Delete the projectile remnant.
ParticleList const & getCreatedParticles() const
Get the list of particles that were created by the last applied final state.
G4int theNpInitial
The number of entering protons.
G4bool decayMe()
Force the phase-space decay of the Nucleus.
ParticleSampler * theParticleSampler
G4int theNnInitial
The number of entering neutrons.
ThreeVector computeCenterOfMass() const
Compute the current center-of-mass position.
void incrementCascading()
G4bool getTryCompoundNucleus()
NuclearPotential::INuclearPotential const * getPotential() const
Getter for thePotential.
G4double theExcitationEnergy
void setProjectileChargeNumber(G4int n)
Set the charge number of the projectile.
void propagateParticles(G4double step)
Propagate the particles one time step.
G4bool decayInsideDeltas()
Force the decay of deltas inside the nucleus.
void finalizeProjectileRemnant(const G4double emissionTime)
Finalise the projectile remnant.
NuclearPotential::INuclearPotential const * thePotential
Pointer to the NuclearPotential object.
G4double getTransmissionBarrier(Particle const *const p)
Get the transmission barrier.
G4int getIsospin(const ParticleType t)
Get the isospin of a particle.
void emitInsidePions()
Force emission of all pions inside the nucleus.
G4INCL::ParticleType getType() const
Get the particle type.
NuclearDensity const * getDensity() const
Getter for theDensity.
void setProjectileRemnant(ProjectileRemnant *const c)
Set the projectile remnant.
G4bool isEventTransparent() const
Is the event transparent?
void fillEventInfo(EventInfo *eventInfo)
Fill the event info which contains INCL output data.
G4bool isNucleon() const
Is this a nucleon?
G4double computeTotalEnergy() const
Compute the current total energy.
ProjectileRemnant * getProjectileRemnant() const
Get the projectile remnant.
void setProjectileMassNumber(G4int n)
Set the mass number of the projectile.
G4double getUniverseRadius() const
Getter for theUniverseRadius.
The purpose of the Store object is to act as a "particle manager" that keeps track ofall the particle...
void insertParticle(Particle *p)
Insert a new particle (e.g. a projectile) in the nucleus.
virtual G4double computePotentialEnergy(const Particle *const p) const =0
Nucleus(G4int mass, G4int charge, Config const *const conf, const G4double universeRadius=-1.)
G4double getTransmissionRadius(Particle const *const p) const
The radius used for calculating the transmission coefficient.
void computeRecoilKinematics()
Compute the recoil momentum and spin of the nucleus.
ThreeVector initialCenterOfMass
Particle * getBlockedDelta() const
Get the delta that could not decay.
G4bool isPion() const
Is this a pion?
void setNucleusNucleusCollision()
Set a nucleus-nucleus collision.
ProjectileRemnant * theProjectileRemnant
Pointer to the quasi-projectile.
G4double getMaxRFromP(const ParticleType t, const G4double p) const
Get the maximum allowed radius for a given momentum.
ParticleList::const_iterator ParticleIter
void particleHasEntered(Particle *const particle)
Move a particle from incoming to inside.
G4double theUniverseRadius
The radius of the universe.