#include <G4INCLParticle.hh>

Inheritance diagram for G4INCL::Particle:

Public Member Functions
	Particle ()

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

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

virtual	~Particle ()

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

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

G4INCL::ParticleType	getType () const

virtual G4INCL::ParticleSpecies	getSpecies () const
	Get the particle species. More...

void	setType (ParticleType t)

G4bool	isNucleon () const

ParticipantType	getParticipantType () const

void	setParticipantType (ParticipantType const p)

G4bool	isParticipant () const

G4bool	isTargetSpectator () const

G4bool	isProjectileSpectator () const

virtual void	makeParticipant ()

virtual void	makeTargetSpectator ()

virtual void	makeProjectileSpectator ()

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

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

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

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

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

G4double	getBeta () const

ThreeVector	boostVector () const

void	boost (const ThreeVector &aBoostVector)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

G4double	getEnergy () const

void	setMass (G4double mass)

void	setEnergy (G4double energy)

const G4INCL::ThreeVector &	getMomentum () const

virtual G4INCL::ThreeVector	getAngularMomentum () const

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

const G4INCL::ThreeVector &	getPosition () const

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

G4double	getHelicity ()

void	setHelicity (G4double h)

void	propagate (G4double step)

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

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

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

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

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

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

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

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

void	setEmissionTime (G4double t)

G4double	getEmissionTime ()

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

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

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

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

G4bool	isInList (ParticleList const &l) const
	Check if the particle belongs to a given list. More...

G4bool	isCluster () const

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

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

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

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

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

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

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

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

std::string	print () const

std::string	dump () const

long	getID () const

ParticleList const *	getParticles () const

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

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

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

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

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

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

Detailed Description

Definition at line 94 of file G4INCLParticle.hh.

Constructor & Destructor Documentation

G4INCL::Particle::Particle ( )

Definition at line 51 of file G4INCLParticle.cc.

References ID.

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

     : theZ(0), theA(0),
     theParticipantType(TargetSpectator),
     theType(UnknownParticle),
     theEnergy(0.0),
     thePropagationEnergy(&theEnergy),
     theFrozenEnergy(theEnergy),
     theMomentum(ThreeVector(0.,0.,0.)),
     thePropagationMomentum(&theMomentum),
     theFrozenMomentum(theMomentum),
     thePosition(ThreeVector(0.,0.,0.)),
     nCollisions(0),
     nDecays(0),
     thePotentialEnergy(0.0),
     rpCorrelated(false),
     uncorrelatedMomentum(0.),
     theHelicity(0.0),
     emissionTime(0.0),
     outOfWell(false),
     theMass(0.)
   {
     ID = nextID;
     nextID++;
   }

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

Definition at line 76 of file G4INCLParticle.cc.

References getInvariantMass(), ID, INCL_WARN, setMass(), setType(), G4INCL::TargetSpectator, theEnergy, and theParticipantType.

     : theEnergy(energy),
     thePropagationEnergy(&theEnergy),
     theFrozenEnergy(theEnergy),
     theMomentum(momentum),
     thePropagationMomentum(&theMomentum),
     theFrozenMomentum(theMomentum),
     thePosition(position),
     nCollisions(0), nDecays(0),
       thePotentialEnergy(0.),
       rpCorrelated(false),
       uncorrelatedMomentum(theMomentum.mag()),
       theHelicity(0.0),
       emissionTime(0.0), outOfWell(false)
   {
     theParticipantType = TargetSpectator;
     ID = nextID;
     nextID++;
     if(theEnergy <= 0.0) {
       INCL_WARN("Particle with energy " << theEnergy << " created." << std::endl);
     }
     setType(t);
     setMass(getInvariantMass());
   }

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

Definition at line 102 of file G4INCLParticle.cc.

References energy(), ID, INCL_ERROR, isResonance(), G4INCL::ThreeVector::mag2(), setType(), G4INCL::TargetSpectator, theEnergy, theFrozenEnergy, theMomentum, and theParticipantType.

     : thePropagationEnergy(&theEnergy),
     theMomentum(momentum),
     thePropagationMomentum(&theMomentum),
     theFrozenMomentum(theMomentum),
     thePosition(position),
     nCollisions(0), nDecays(0),
       thePotentialEnergy(0.),
       rpCorrelated(false),
       uncorrelatedMomentum(theMomentum.mag()),
       theHelicity(0.0),
       emissionTime(0.0), outOfWell(false)
   {
     theParticipantType = TargetSpectator;
     ID = nextID;
     nextID++;
     setType(t);
     if( isResonance() ) {
       INCL_ERROR("Cannot create resonance without specifying its momentum four-vector." << std::endl);
     }
     G4double energy = std::sqrt(theMomentum.mag2() + theMass*theMass);
     theEnergy = energy;
     theFrozenEnergy = theEnergy;
   }

virtual G4INCL::Particle::~Particle ( )

inlinevirtual

Definition at line 99 of file G4INCLParticle.hh.

99 {}

G4INCL::Particle::Particle ( const Particle & rhs )

inline

Copy constructor.

Does not copy the particle ID.

Definition at line 105 of file G4INCLParticle.hh.

References ID, theEnergy, theFrozenEnergy, theFrozenMomentum, theMomentum, thePropagationEnergy, and thePropagationMomentum.

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

Member Function Documentation

G4double G4INCL::Particle::adjustEnergyFromMomentum ( )

Recompute the energy to match the momentum.

Definition at line 141 of file G4INCLParticle.cc.

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

Referenced by G4INCL::Nucleus::computeRecoilKinematics(), G4INCL::PhaseSpaceDecay::decay(), G4INCL::Nucleus::decayOutgoingDeltas(), G4INCL::DeltaDecayChannel::getFinalState(), and G4INCL::RecombinationChannel::getFinalState().

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

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

Rescale the momentum to match the total energy.

Definition at line 128 of file G4INCLParticle.cc.

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

Referenced by G4INCL::Cluster::Cluster(), G4INCL::Nucleus::emitInsidePions(), G4INCL::CrossSections::interactionDistanceNN(), G4INCL::CrossSections::interactionDistancePiN(), G4INCL::StandardPropagationModel::shootParticle(), and G4INCL::KinematicsUtils::transformToLocalEnergyFrame().

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

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

inline

Boost the particle using a boost vector.

Example (go to the particle rest frame): particle->boost(particle->boostVector());

Definition at line 327 of file G4INCLParticle.hh.

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

Referenced by G4INCL::Cluster::boost(), G4INCL::PhaseSpaceDecay::decay(), G4INCL::Nucleus::decayOutgoingDeltas(), and G4INCL::InteractionAvatar::preInteraction().

                                                 {
       const G4double beta2 = aBoostVector.mag2();
       const G4double gamma = 1.0 / std::sqrt(1.0 - beta2);
       const G4double bp = theMomentum.dot(aBoostVector);
       const G4double alpha = (gamma*gamma)/(1.0 + gamma);
 
       theMomentum = theMomentum + aBoostVector * (alpha * bp - gamma * theEnergy);
       theEnergy = gamma * (theEnergy - bp);
     }

ThreeVector G4INCL::Particle::boostVector ( ) const

inline

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 317 of file G4INCLParticle.hh.

References theEnergy, and theMomentum.

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

                                     {
       return theMomentum / theEnergy;
     }

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

inline

Definition at line 735 of file G4INCLParticle.hh.

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

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

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

void G4INCL::Particle::freezePropagation ( )

inline

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 693 of file G4INCLParticle.hh.

References theFrozenEnergy, theFrozenMomentum, thePropagationEnergy, and thePropagationMomentum.

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

G4int G4INCL::Particle::getA ( ) const

inline

Returns the baryon number.

Definition at line 301 of file G4INCLParticle.hh.

References theA.

Referenced by G4INCL::Cluster::addParticle(), G4INCL::ClusteringModelIntercomparison::clusterCanEscape(), G4INCL::ClusterDecay::decay(), G4INCL::Nucleus::fillEventInfo(), G4INCL::Nucleus::finalizeProjectileRemnant(), G4INCL::PauliStandard::getBlockingProbability(), G4INCL::ClusteringModelIntercomparison::getCluster(), G4INCL::Nucleus::getConservationBalance(), G4INCL::ParticleEntryChannel::getFinalState(), G4INCL::TransmissionChannel::getFinalState(), G4INCL::SurfaceAvatar::getTransmissionProbability(), G4INCL::NuclearDensity::getTransmissionRadius(), G4INCL::Nucleus::insertParticle(), G4INCL::ClusterDecay::isStable(), G4INCL::ProjectileRemnant::removeParticle(), G4INCL::StandardPropagationModel::shootComposite(), and G4INCL::StandardPropagationModel::shootParticle().

301 { return theA; }

G4INCL::Particle::theA

G4int theA

Definition: G4INCLParticle.hh:790

virtual G4INCL::ThreeVector G4INCL::Particle::getAngularMomentum ( ) const

inlinevirtual

Get the angular momentum w.r.t. the origin

Reimplemented in G4INCL::Cluster.

Definition at line 580 of file G4INCLParticle.hh.

References theMomentum, thePosition, and G4INCL::ThreeVector::vector().

Referenced by G4INCL::Cluster::getAngularMomentum(), and G4INCL::StandardPropagationModel::shootParticle().

     {
       return thePosition.vector(theMomentum);
     };

G4double G4INCL::Particle::getBeta ( ) const

inline

Definition at line 306 of file G4INCLParticle.hh.

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

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

G4double G4INCL::Particle::getCosRPAngle ( ) const

inline

Get the cosine of the angle between position and momentum.

Definition at line 781 of file G4INCLParticle.hh.

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

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

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

inline

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 465 of file G4INCLParticle.hh.

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

Referenced by G4INCL::Nucleus::emitInsidePions(), G4INCL::ParticleEntryChannel::getFinalState(), G4INCL::TransmissionChannel::getFinalState(), and G4INCL::SurfaceAvatar::getTransmissionProbability().

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

G4double G4INCL::Particle::getEmissionTime ( )

inline

Definition at line 645 of file G4INCLParticle.hh.

645 { return emissionTime; };

G4double G4INCL::Particle::getEnergy ( ) const

inline

Get the energy of the particle in MeV.

Definition at line 550 of file G4INCLParticle.hh.

References theEnergy.

Referenced by G4INCL::Cluster::addParticle(), G4INCL::DeltaDecayChannel::computeDecayTime(), G4INCL::ClusteringModelIntercomparison::getCluster(), G4INCL::PionNucleonChannel::getFinalState(), G4INCL::ReflectionChannel::getFinalState(), G4INCL::ParticleEntryChannel::getFinalState(), G4INCL::TransmissionChannel::getFinalState(), G4INCL::KinematicsUtils::makeBoostVector(), G4INCL::KinematicsUtils::momentumInCM(), G4INCL::InteractionAvatar::preInteraction(), G4INCL::InteractionAvatar::preInteractionBlocking(), G4INCL::ProjectileRemnant::removeParticle(), G4INCL::StandardPropagationModel::shootComposite(), G4INCL::StandardPropagationModel::shootParticle(), G4INCL::KinematicsUtils::squareTotalEnergyInCM(), and G4INCL::KinematicsUtils::transformToLocalEnergyFrame().

     {
       return theEnergy;
     };

G4double G4INCL::Particle::getFrozenEnergy ( ) const

inline

Get the frozen particle momentum.

Definition at line 682 of file G4INCLParticle.hh.

References theFrozenEnergy.

682 { return theFrozenEnergy; }

G4INCL::Particle::theFrozenEnergy

G4double theFrozenEnergy

Definition: G4INCLParticle.hh:795

ThreeVector G4INCL::Particle::getFrozenMomentum ( ) const

inline

Get the frozen particle momentum.

Definition at line 679 of file G4INCLParticle.hh.

References theFrozenMomentum.

679 { return theFrozenMomentum; }

G4INCL::Particle::theFrozenMomentum

G4INCL::ThreeVector theFrozenMomentum

Definition: G4INCLParticle.hh:798

G4double G4INCL::Particle::getHelicity ( )

inline

Definition at line 606 of file G4INCLParticle.hh.

606 { return theHelicity; };

long G4INCL::Particle::getID ( ) const

inline

Definition at line 748 of file G4INCLParticle.hh.

References ID.

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

748 { return ID; };

G4INCL::Particle::ID

long ID

Definition: G4INCLParticle.hh:803

G4double G4INCL::Particle::getINCLMass ( ) const

inline

Get the INCL particle mass.

Definition at line 359 of file G4INCLParticle.hh.

References G4INCL::Composite, G4INCL::DeltaMinus, G4INCL::DeltaPlus, G4INCL::DeltaPlusPlus, G4INCL::DeltaZero, G4INCL::ParticleTable::getINCLMass(), INCL_ERROR, G4INCL::Neutron, G4INCL::PiMinus, G4INCL::PiPlus, G4INCL::PiZero, G4INCL::Proton, theA, theType, and theZ.

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

                                         {
       switch(theType) {
         case Proton:
         case Neutron:
         case PiPlus:
         case PiMinus:
         case PiZero:
           return ParticleTable::getINCLMass(theType);
           break;
 
         case DeltaPlusPlus:
         case DeltaPlus:
         case DeltaZero:
         case DeltaMinus:
           return theMass;
           break;
 
         case Composite:
           return ParticleTable::getINCLMass(theA,theZ);
           break;
 
         default:
           INCL_ERROR("Particle::getINCLMass: Unknown particle type." << std::endl);
           return 0.0;
           break;
       }
     }

G4double G4INCL::Particle::getInvariantMass ( ) const

inline

Get the the particle invariant mass.

Uses the relativistic invariant

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

Definition at line 528 of file G4INCLParticle.hh.

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

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

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

G4double G4INCL::Particle::getKineticEnergy ( ) const

inline

Get the particle kinetic energy.

Definition at line 539 of file G4INCLParticle.hh.

References theEnergy.

Referenced by G4INCL::NuclearPotential::NuclearPotentialEnergyIsospin::computePotentialEnergy(), G4INCL::NuclearPotential::NuclearPotentialEnergyIsospinSmooth::computePotentialEnergy(), G4INCL::Nucleus::emitInsidePions(), G4INCL::Nucleus::fillEventInfo(), G4INCL::SurfaceAvatar::getChannel(), G4INCL::Nucleus::getConservationBalance(), G4INCL::ParticleEntryChannel::getFinalState(), G4INCL::KinematicsUtils::getLocalEnergy(), G4INCL::SurfaceAvatar::getTransmissionProbability(), G4INCL::CoulombDistortion::maxImpactParameter(), and G4INCL::CDPP::processOneParticle().

539 { return theEnergy - theMass; }

G4INCL::Particle::theEnergy

G4double theEnergy

Definition: G4INCLParticle.hh:793

ThreeVector G4INCL::Particle::getLongitudinalPosition ( ) const

inline

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

Definition at line 653 of file G4INCLParticle.hh.

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

Referenced by getTransversePosition().

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

G4double G4INCL::Particle::getMass ( ) const

inline

Get the cached particle mass.

Definition at line 356 of file G4INCLParticle.hh.

Referenced by G4INCL::DeltaDecayChannel::computeDecayTime(), G4INCL::PhaseSpaceDecay::decay(), G4INCL::Nucleus::emitInsidePions(), G4INCL::InteractionAvatar::enforceEnergyConservation(), G4INCL::Cluster::freezeInternalMotion(), G4INCL::NuclearPotential::INuclearPotential::getFermiMomentum(), G4INCL::DeltaDecayChannel::getFinalState(), G4INCL::TransmissionChannel::getFinalState(), G4INCL::RecombinationChannel::getFinalState(), G4INCL::KinematicsUtils::getLocalEnergy(), G4INCL::SurfaceAvatar::getTransmissionProbability(), G4INCL::CrossSections::interactionDistanceNN(), G4INCL::CrossSections::interactionDistancePiN(), G4INCL::Cluster::internalBoostToCM(), G4INCL::KinematicsUtils::momentumInCM(), G4INCL::KinematicsUtils::momentumInLab(), G4INCL::Cluster::print(), G4INCL::ProjectileRemnant::ProjectileRemnant(), G4INCL::CrossSectionsINCL46::recombination(), and G4INCL::StandardPropagationModel::shootParticle().

356 { return theMass; }

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

inline

Get the momentum vector.

Definition at line 574 of file G4INCLParticle.hh.

References theMomentum.

Referenced by G4INCL::Cluster::addParticle(), G4INCL::ClusteringModelIntercomparison::clusterCanEscape(), G4INCL::Nucleus::computeRecoilKinematics(), G4INCL::Nucleus::decayOutgoingDeltas(), G4INCL::Nucleus::fillEventInfo(), G4INCL::PauliStandard::getBlockingProbability(), G4INCL::ClusteringModelIntercomparison::getCluster(), G4INCL::Nucleus::getConservationBalance(), G4INCL::PionNucleonChannel::getFinalState(), G4INCL::ReflectionChannel::getFinalState(), G4INCL::DeltaProductionChannel::getFinalState(), G4INCL::ParticleEntryChannel::getFinalState(), G4INCL::ElasticChannel::getFinalState(), G4INCL::KinematicsUtils::getLocalEnergy(), G4INCL::SurfaceAvatar::getTransmissionProbability(), G4INCL::KinematicsUtils::makeBoostVector(), G4INCL::KinematicsUtils::momentumInCM(), G4INCL::BinaryCollisionAvatar::postInteraction(), G4INCL::InteractionAvatar::preInteraction(), G4INCL::ProjectileRemnant::removeParticle(), G4INCL::ParticleSampler::sampleParticles(), G4INCL::StandardPropagationModel::shootComposite(), and G4INCL::StandardPropagationModel::shootParticle().

     {
       return theMomentum;
     };

G4int G4INCL::Particle::getNumberOfCollisions ( ) const

inline

Return the number of collisions undergone by the particle.

Definition at line 614 of file G4INCLParticle.hh.

References nCollisions.

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

614 { return nCollisions; }

G4INCL::Particle::nCollisions

G4int nCollisions

Definition: G4INCLParticle.hh:800

G4int G4INCL::Particle::getNumberOfDecays ( ) const

inline

Return the number of decays undergone by the particle.

Definition at line 623 of file G4INCLParticle.hh.

References nDecays.

623 { return nDecays; }

G4INCL::Particle::nDecays

G4int nDecays

Definition: G4INCLParticle.hh:801

ParticipantType G4INCL::Particle::getParticipantType ( ) const

inline

Definition at line 256 of file G4INCLParticle.hh.

References theParticipantType.

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

                                                {
       return theParticipantType;
     }

ParticleList const* G4INCL::Particle::getParticles ( ) const

inline

Return a NULL pointer

Definition at line 753 of file G4INCLParticle.hh.

References INCL_WARN.

                                              {
       INCL_WARN("Particle::getParticles() method was called on a Particle object" << std::endl);
       return 0;
     }

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

inline

Set the position vector.

Definition at line 596 of file G4INCLParticle.hh.

References thePosition.

Referenced by G4INCL::Cluster::addParticle(), G4INCL::InteractionAvatar::bringParticleInside(), G4INCL::CoulombNone::bringToSurface(), G4INCL::ClusteringModelIntercomparison::clusterCanEscape(), G4INCL::StandardPropagationModel::generateBinaryCollisionAvatar(), G4INCL::PauliStandard::getBlockingProbability(), G4INCL::BinaryCollisionAvatar::getChannel(), G4INCL::ClusteringModelIntercomparison::getCluster(), G4INCL::ReflectionChannel::getFinalState(), G4INCL::DeltaDecayChannel::getFinalState(), G4INCL::KinematicsUtils::getLocalEnergy(), G4INCL::StandardPropagationModel::getReflectionTime(), G4INCL::StandardPropagationModel::getTime(), G4INCL::BinaryCollisionAvatar::postInteraction(), and G4INCL::ParticleSampler::sampleParticles().

     {
       return thePosition;
     };

G4double G4INCL::Particle::getPotentialEnergy ( ) const

inline

Get the particle potential energy.

Definition at line 542 of file G4INCLParticle.hh.

References thePotentialEnergy.

Referenced by G4INCL::Cluster::addParticle(), G4INCL::Nucleus::emitInsidePions(), G4INCL::ClusteringModelIntercomparison::getCluster(), G4INCL::ReflectionChannel::getFinalState(), G4INCL::PionNucleonChannel::getFinalState(), G4INCL::TransmissionChannel::getFinalState(), G4INCL::KinematicsUtils::getLocalEnergy(), G4INCL::SurfaceAvatar::getTransmissionProbability(), G4INCL::InteractionAvatar::preInteractionBlocking(), and G4INCL::CDPP::processOneParticle().

542 { return thePotentialEnergy; }

G4INCL::Particle::thePotentialEnergy

G4double thePotentialEnergy

Definition: G4INCLParticle.hh:802

ThreeVector G4INCL::Particle::getPropagationVelocity ( ) const

inline

Get the propagation velocity of the particle.

Definition at line 685 of file G4INCLParticle.hh.

References thePropagationMomentum.

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

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

G4INCL::Particle::thePropagationMomentum

G4INCL::ThreeVector * thePropagationMomentum

Definition: G4INCLParticle.hh:797

G4double G4INCL::Particle::getRealMass ( ) const

inline

Get the real particle mass.

Definition at line 417 of file G4INCLParticle.hh.

References G4INCL::Composite, G4INCL::DeltaMinus, G4INCL::DeltaPlus, G4INCL::DeltaPlusPlus, G4INCL::DeltaZero, G4INCL::ParticleTable::getRealMass(), INCL_ERROR, G4INCL::Neutron, G4INCL::PiMinus, G4INCL::PiPlus, G4INCL::PiZero, G4INCL::Proton, theA, theType, and theZ.

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

                                         {
       switch(theType) {
         case Proton:
         case Neutron:
         case PiPlus:
         case PiMinus:
         case PiZero:
           return ParticleTable::getRealMass(theType);
           break;
 
         case DeltaPlusPlus:
         case DeltaPlus:
         case DeltaZero:
         case DeltaMinus:
           return theMass;
           break;
 
         case Composite:
           return ParticleTable::getRealMass(theA,theZ);
           break;
 
         default:
           INCL_ERROR("Particle::getRealMass: Unknown particle type." << std::endl);
           return 0.0;
           break;
       }
     }

G4double G4INCL::Particle::getReflectionMomentum ( ) const

inline

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 764 of file G4INCLParticle.hh.

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

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

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

virtual G4INCL::ParticleSpecies G4INCL::Particle::getSpecies ( ) const

inlinevirtual

Get the particle species.

Reimplemented in G4INCL::Cluster.

Definition at line 192 of file G4INCLParticle.hh.

References theType.

Referenced by G4INCL::CoulombDistortion::maxImpactParameter(), and G4INCL::StandardPropagationModel::shootParticle().

                                                      {
       return ParticleSpecies(theType);
     };

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

inlinevirtual

Get the tabulated particle mass.

Reimplemented in G4INCL::Cluster.

Definition at line 388 of file G4INCLParticle.hh.

References G4INCL::Composite, G4INCL::DeltaMinus, G4INCL::DeltaPlus, G4INCL::DeltaPlusPlus, G4INCL::DeltaZero, G4INCL::ParticleTable::getTableMass, G4INCL::ParticleTable::getTableParticleMass, INCL_ERROR, G4INCL::Neutron, G4INCL::PiMinus, G4INCL::PiPlus, G4INCL::PiZero, G4INCL::Proton, theA, theType, and theZ.

Referenced by G4INCL::Nucleus::decayOutgoingDeltas(), getEmissionQValueCorrection(), G4INCL::ParticleEntryChannel::getFinalState(), G4INCL::TransmissionChannel::getFinalState(), getTransferQValueCorrection(), and setTableMass().

                                                  {
       switch(theType) {
         case Proton:
         case Neutron:
         case PiPlus:
         case PiMinus:
         case PiZero:
           return ParticleTable::getTableParticleMass(theType);
           break;
 
         case DeltaPlusPlus:
         case DeltaPlus:
         case DeltaZero:
         case DeltaMinus:
           return theMass;
           break;
 
         case Composite:
           return ParticleTable::getTableMass(theA,theZ);
           break;
 
         default:
           INCL_ERROR("Particle::getTableMass: Unknown particle type." << std::endl);
           return 0.0;
           break;
       }
     }

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

inline

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 503 of file G4INCLParticle.hh.

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

                                                                                                                        {
       const G4int AFromDaughter = AFrom - theA;
       const G4int ZFromDaughter = ZFrom - theZ;
       const G4int AToDaughter = ATo + theA;
       const G4int ZToDaughter = ZTo + theZ;
       const G4double theQValue = ParticleTable::getTableQValue(AToDaughter,ZToDaughter,AFromDaughter,ZFromDaughter,AFrom,ZFrom);
 
       const G4double massINCLTo = ParticleTable::getINCLMass(ATo,ZTo);
       const G4double massINCLToDaughter = ParticleTable::getINCLMass(AToDaughter,ZToDaughter);
       /* 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);
     }

ThreeVector G4INCL::Particle::getTransversePosition ( ) const

inline

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

Definition at line 648 of file G4INCLParticle.hh.

References getLongitudinalPosition(), and thePosition.

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

                                               {
       return thePosition - getLongitudinalPosition();
     }

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

inline

Get the particle type.

See also: G4INCL::ParticleType

Definition at line 187 of file G4INCLParticle.hh.

References theType.

Referenced by G4INCL::NuclearPotential::INuclearPotential::computePionPotentialEnergy(), G4INCL::NuclearPotential::NuclearPotentialIsospin::computePotentialEnergy(), G4INCL::NuclearPotential::NuclearPotentialConstant::computePotentialEnergy(), G4INCL::CrossSectionsINCL46::deltaProduction(), G4INCL::CrossSectionsINCL46::elasticNNLegacy(), G4INCL::PauliStandard::getBlockingProbability(), G4INCL::NuclearPotential::INuclearPotential::getFermiEnergy(), G4INCL::NuclearPotential::INuclearPotential::getFermiMomentum(), G4INCL::DeltaProductionChannel::getFinalState(), G4INCL::DeltaDecayChannel::getFinalState(), G4INCL::ElasticChannel::getFinalState(), G4INCL::RecombinationChannel::getFinalState(), G4INCL::KinematicsUtils::getLocalEnergy(), G4INCL::NuclearPotential::INuclearPotential::getSeparationEnergy(), G4INCL::Nucleus::getSurfaceRadius(), G4INCL::NuclearDensity::getTransmissionRadius(), G4INCL::Nucleus::insertParticle(), G4INCL::ParticleConfig::isPair(), G4INCL::CrossSectionsINCL46::pionNucleon(), and G4INCL::CrossSectionsINCL46::recombination().

                                        {
       return theType;
     };

G4int G4INCL::Particle::getZ ( ) const

inline

Returns the charge number.

Definition at line 304 of file G4INCLParticle.hh.

References theZ.

Referenced by G4INCL::Cluster::addParticle(), G4INCL::CoulombNonRelativistic::bringToSurface(), G4INCL::ClusterDecay::decay(), G4INCL::CoulombNonRelativistic::distortOut(), G4INCL::Nucleus::emitInsidePions(), G4INCL::Nucleus::fillEventInfo(), G4INCL::Nucleus::finalizeProjectileRemnant(), G4INCL::PauliStandard::getBlockingProbability(), G4INCL::ClusteringModelIntercomparison::getCluster(), G4INCL::Nucleus::getConservationBalance(), G4INCL::ParticleEntryChannel::getFinalState(), G4INCL::TransmissionChannel::getFinalState(), G4INCL::Nucleus::getTransmissionBarrier(), G4INCL::SurfaceAvatar::getTransmissionProbability(), G4INCL::NuclearDensity::getTransmissionRadius(), G4INCL::Nucleus::insertParticle(), G4INCL::ClusterDecay::isStable(), G4INCL::ProjectileRemnant::removeParticle(), G4INCL::StandardPropagationModel::shootComposite(), and G4INCL::StandardPropagationModel::shootParticle().

304 { return theZ; }

G4INCL::Particle::theZ

G4int theZ

Definition: G4INCLParticle.hh:790

void G4INCL::Particle::incrementNumberOfCollisions ( )

inline

Increment the number of collisions undergone by the particle.

Definition at line 620 of file G4INCLParticle.hh.

References nCollisions.

620 { nCollisions++; }

G4INCL::Particle::nCollisions

G4int nCollisions

Definition: G4INCLParticle.hh:800

void G4INCL::Particle::incrementNumberOfDecays ( )

inline

Increment the number of decays undergone by the particle.

Definition at line 629 of file G4INCLParticle.hh.

References nDecays.

629 { nDecays++; }

G4INCL::Particle::nDecays

G4int nDecays

Definition: G4INCLParticle.hh:801

G4bool G4INCL::Particle::isCluster ( ) const

inline

Definition at line 668 of file G4INCLParticle.hh.

References G4INCL::Composite, and theType.

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

                              {
       return (theType == Composite);
     }

G4bool G4INCL::Particle::isDelta ( ) const

inline

Is it a Delta?

Definition at line 295 of file G4INCLParticle.hh.

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

Referenced by G4INCL::BinaryCollisionAvatar::getChannel(), G4INCL::DecayAvatar::getChannel(), G4INCL::NuclearPotential::INuclearPotential::getFermiMomentum(), G4INCL::DeltaProductionChannel::getFinalState(), isResonance(), G4INCL::CrossSectionsINCL46::recombination(), G4INCL::RecombinationChannel::RecombinationChannel(), and G4INCL::CrossSectionsINCL46::total().

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

G4bool G4INCL::Particle::isInList ( ParticleList const & l ) const

inline

Check if the particle belongs to a given list.

Definition at line 664 of file G4INCLParticle.hh.

                                                  {
       return (std::find(l.begin(), l.end(), this)!=l.end());
     }

G4bool G4INCL::Particle::isNucleon ( ) const

inline

Is this a nucleon?

Definition at line 249 of file G4INCLParticle.hh.

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

Referenced by G4INCL::NuclearPotential::NuclearPotentialEnergyIsospin::computePotentialEnergy(), G4INCL::NuclearPotential::NuclearPotentialEnergyIsospinSmooth::computePotentialEnergy(), G4INCL::CrossSectionsINCL46::elasticNNLegacy(), G4INCL::StandardPropagationModel::generateBinaryCollisionAvatar(), G4INCL::BinaryCollisionAvatar::getChannel(), G4INCL::SurfaceAvatar::getChannel(), G4INCL::PionNucleonChannel::getFinalState(), G4INCL::ParticleEntryChannel::getFinalState(), G4INCL::Nucleus::insertParticle(), G4INCL::CrossSectionsINCL46::pionNucleon(), G4INCL::CDPP::processOneParticle(), G4INCL::StandardPropagationModel::shootParticle(), and G4INCL::CrossSectionsINCL46::total().

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

G4bool G4INCL::Particle::isOutOfWell ( ) const

inline

Check if the particle is out of its potential well.

Definition at line 642 of file G4INCLParticle.hh.

Referenced by G4INCL::NuclearPotential::INuclearPotential::computePionPotentialEnergy().

642 { return outOfWell; }

G4bool G4INCL::Particle::isParticipant ( ) const

inline

Definition at line 264 of file G4INCLParticle.hh.

References G4INCL::Participant, and theParticipantType.

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

                                  {
       return (theParticipantType==Participant);
     }

G4bool G4INCL::Particle::isPion ( ) const

inline

Is this a pion?

Definition at line 289 of file G4INCLParticle.hh.

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

Referenced by G4INCL::CrossSectionsINCL46::elastic(), G4INCL::StandardPropagationModel::generateBinaryCollisionAvatar(), G4INCL::BinaryCollisionAvatar::getChannel(), G4INCL::Nucleus::getSurfaceRadius(), G4INCL::CrossSectionsINCL46::pionNucleon(), G4INCL::InteractionAvatar::preInteractionLocalEnergy(), G4INCL::CDPP::processOneParticle(), and G4INCL::CrossSectionsINCL46::total().

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

G4INCL::PiPlus

Definition: G4INCLParticleType.hh:52

G4INCL::PiZero

Definition: G4INCLParticleType.hh:54

G4INCL::PiMinus

Definition: G4INCLParticleType.hh:53

G4INCL::Particle::theType

G4INCL::ParticleType theType

Definition: G4INCLParticle.hh:792

G4bool G4INCL::Particle::isProjectileSpectator ( ) const

inline

Definition at line 272 of file G4INCLParticle.hh.

References G4INCL::ProjectileSpectator, and theParticipantType.

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

                                          {
       return (theParticipantType==ProjectileSpectator);
     }

G4bool G4INCL::Particle::isResonance ( ) const

inline

Is it a resonance?

Definition at line 292 of file G4INCLParticle.hh.

References isDelta().

Referenced by G4INCL::StandardPropagationModel::generateBinaryCollisionAvatar(), G4INCL::SurfaceAvatar::getChannel(), Particle(), G4INCL::CDPP::processOneParticle(), and setType().

292 { return isDelta(); }

G4INCL::Particle::isDelta

G4bool isDelta() const

Is it a Delta?

Definition: G4INCLParticle.hh:295

G4bool G4INCL::Particle::isTargetSpectator ( ) const

inline

Definition at line 268 of file G4INCLParticle.hh.

References G4INCL::TargetSpectator, and theParticipantType.

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

                                      {
       return (theParticipantType==TargetSpectator);
     }

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

inline

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 345 of file G4INCLParticle.hh.

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

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

virtual void G4INCL::Particle::makeParticipant ( )

inlinevirtual

Reimplemented in G4INCL::Cluster.

Definition at line 276 of file G4INCLParticle.hh.

References G4INCL::Participant, and theParticipantType.

Referenced by G4INCL::Store::loadParticles(), and G4INCL::Cluster::makeParticipant().

                                    {
       theParticipantType = Participant;
     }

virtual void G4INCL::Particle::makeProjectileSpectator ( )

inlinevirtual

Reimplemented in G4INCL::Cluster.

Definition at line 284 of file G4INCLParticle.hh.

References G4INCL::ProjectileSpectator, and theParticipantType.

Referenced by G4INCL::Cluster::makeProjectileSpectator(), and G4INCL::StandardPropagationModel::shootParticle().

                                            {
       theParticipantType = ProjectileSpectator;
     }

virtual void G4INCL::Particle::makeTargetSpectator ( )

inlinevirtual

Reimplemented in G4INCL::Cluster.

Definition at line 280 of file G4INCLParticle.hh.

References G4INCL::TargetSpectator, and theParticipantType.

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

                                        {
       theParticipantType = TargetSpectator;
     }

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

inline

Assignment operator.

Does not copy the particle ID.

Definition at line 177 of file G4INCLParticle.hh.

References swap().

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

                                              {
       Particle temporaryParticle(rhs);
       swap(temporaryParticle);
       return *this;
     }

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

inline

Definition at line 720 of file G4INCLParticle.hh.

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

Referenced by adjustMomentumFromEnergy(), G4INCL::Nucleus::emitInsidePions(), G4INCL::BinaryCollisionAvatar::getChannel(), G4INCL::ParticleEntryChannel::getFinalState(), G4INCL::KinematicsUtils::getLocalEnergy(), G4INCL::StandardPropagationModel::getReflectionTime(), and G4INCL::ProjectileRemnant::removeParticle().

                             {
       std::stringstream ss;
       ss << "Particle (ID = " << ID << ") type = ";
       ss << ParticleTable::getName(theType);
       ss << std::endl
         << "   energy = " << theEnergy << std::endl
         << "   momentum = "
         << theMomentum.print()
         << std::endl
         << "   position = "
         << thePosition.print()
         << std::endl;
       return ss.str();
     };

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

inline

Definition at line 609 of file G4INCLParticle.hh.

References thePosition.

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

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

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

inlinevirtual

Rotate the particle position and momentum.

Parameters

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

Reimplemented in G4INCL::Cluster.

Definition at line 714 of file G4INCLParticle.hh.

References G4INCL::ThreeVector::rotate(), theFrozenMomentum, theMomentum, and thePosition.

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

                                                                        {
       thePosition.rotate(angle, axis);
       theMomentum.rotate(angle, axis);
       theFrozenMomentum.rotate(angle, axis);
     }

void G4INCL::Particle::rpCorrelate ( )

inline

Make the particle follow a strict r-p correlation.

Definition at line 775 of file G4INCLParticle.hh.

References rpCorrelated.

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

775 { rpCorrelated = true; }

G4INCL::Particle::rpCorrelated

G4bool rpCorrelated

Definition: G4INCLParticle.hh:805

void G4INCL::Particle::rpDecorrelate ( )

inline

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

Definition at line 778 of file G4INCLParticle.hh.

References rpCorrelated.

778 { rpCorrelated = false; }

G4INCL::Particle::rpCorrelated

G4bool rpCorrelated

Definition: G4INCLParticle.hh:805

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

inline

Definition at line 644 of file G4INCLParticle.hh.

Referenced by G4INCL::Nucleus::decayOutgoingDeltas(), G4INCL::Nucleus::emitInsidePions(), and G4INCL::Nucleus::finalizeProjectileRemnant().

644 { emissionTime = t; }

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

inline

Set the energy of the particle in MeV.

Definition at line 566 of file G4INCLParticle.hh.

References energy(), and theEnergy.

Referenced by G4INCL::Nucleus::emitInsidePions(), G4INCL::DeltaProductionChannel::getFinalState(), G4INCL::PionNucleonChannel::getFinalState(), G4INCL::CrossSections::interactionDistanceNN(), G4INCL::CrossSections::interactionDistancePiN(), G4INCL::StandardPropagationModel::shootParticle(), G4INCL::KinematicsUtils::transformToLocalEnergyFrame(), and G4INCL::Nucleus::useFusionKinematics().

     {
       this->theEnergy = energy;
     };

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

inline

Set the frozen particle momentum.

Definition at line 676 of file G4INCLParticle.hh.

References energy(), and theFrozenEnergy.

676 { theFrozenEnergy = energy; }

energy

double precision function energy(A, Z)

Definition: dpm25nuc6.f:4106

G4INCL::Particle::theFrozenEnergy

G4double theFrozenEnergy

Definition: G4INCLParticle.hh:795

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

inline

Set the frozen particle momentum.

Definition at line 673 of file G4INCLParticle.hh.

References theFrozenMomentum.

673 { theFrozenMomentum = momentum; }

G4INCL::Particle::theFrozenMomentum

G4INCL::ThreeVector theFrozenMomentum

Definition: G4INCLParticle.hh:798

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

inline

Definition at line 607 of file G4INCLParticle.hh.

Referenced by G4INCL::DeltaProductionChannel::getFinalState(), and G4INCL::DeltaDecayChannel::getFinalState().

607 { theHelicity = h; };

void G4INCL::Particle::setINCLMass ( )

inline

Set the mass of the Particle to its table mass.

Definition at line 452 of file G4INCLParticle.hh.

References getINCLMass(), and setMass().

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

452 { setMass(getINCLMass()); }

G4INCL::Particle::setMass

void setMass(G4double mass)

Definition: G4INCLParticle.hh:558

G4INCL::Particle::getINCLMass

G4double getINCLMass() const

Get the INCL particle mass.

Definition: G4INCLParticle.hh:359

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

inline

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

Definition at line 558 of file G4INCLParticle.hh.

Referenced by G4INCL::Nucleus::computeRecoilKinematics(), G4INCL::PhaseSpaceDecay::decay(), G4INCL::Nucleus::finalizeProjectileRemnant(), G4INCL::DeltaProductionChannel::getFinalState(), G4INCL::PionNucleonChannel::getFinalState(), Particle(), setINCLMass(), setRealMass(), setTableMass(), and G4INCL::Nucleus::useFusionKinematics().

     {
       this->theMass = mass;
     }

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

inlinevirtual

Set the momentum vector.

Definition at line 588 of file G4INCLParticle.hh.

References theMomentum.

Referenced by G4INCL::PhaseSpaceDecay::decay(), G4INCL::Nucleus::decayOutgoingDeltas(), G4INCL::DeltaProductionChannel::getFinalState(), G4INCL::ReflectionChannel::getFinalState(), G4INCL::PionNucleonChannel::getFinalState(), G4INCL::DeltaDecayChannel::getFinalState(), G4INCL::ElasticChannel::getFinalState(), G4INCL::RecombinationChannel::getFinalState(), and G4INCL::Nucleus::useFusionKinematics().

     {
       this->theMomentum = momentum;
     };

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

inline

Set the number of collisions undergone by the particle.

Definition at line 617 of file G4INCLParticle.hh.

References n, and nCollisions.

617 { nCollisions = n; }

G4INCL::Particle::nCollisions

G4int nCollisions

Definition: G4INCLParticle.hh:800

n

const G4int n

Definition: G4UrQMD1_3Interface.hh:144

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

inline

Set the number of decays undergone by the particle.

Definition at line 626 of file G4INCLParticle.hh.

References n, and nDecays.

626 { nDecays = n; }

n

const G4int n

Definition: G4UrQMD1_3Interface.hh:144

G4INCL::Particle::nDecays

G4int nDecays

Definition: G4INCLParticle.hh:801

void G4INCL::Particle::setOutOfWell ( )

inline

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 639 of file G4INCLParticle.hh.

639 { outOfWell = true; }

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

inline

Definition at line 260 of file G4INCLParticle.hh.

References theParticipantType.

                                                      {
       theParticipantType = p;
     }

virtual void G4INCL::Particle::setPosition ( const G4INCL::ThreeVector & position )

inlinevirtual

Reimplemented in G4INCL::Cluster.

Definition at line 601 of file G4INCLParticle.hh.

References position, and thePosition.

Referenced by G4INCL::InteractionAvatar::bringParticleInside(), G4INCL::CoulombNone::bringToSurface(), G4INCL::ClusteringModelIntercomparison::getCluster(), G4INCL::ReflectionChannel::getFinalState(), G4INCL::ParticleSampler::sampleParticles(), G4INCL::Cluster::setPosition(), and G4INCL::StandardPropagationModel::shootParticle().

     {
       this->thePosition = position;
     };

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

inline

Set the particle potential energy.

Definition at line 545 of file G4INCLParticle.hh.

References thePotentialEnergy, and test::v.

Referenced by G4INCL::Nucleus::emitInsidePions(), G4INCL::Store::loadParticles(), and G4INCL::Nucleus::updatePotentialEnergy().

545 { thePotentialEnergy = v; }

test.v

tuple v

Definition: tests/test12/test.py:18

G4INCL::Particle::thePotentialEnergy

G4double thePotentialEnergy

Definition: G4INCLParticle.hh:802

void G4INCL::Particle::setRealMass ( )

inline

Set the mass of the Particle to its real mass.

Definition at line 446 of file G4INCLParticle.hh.

References getRealMass(), and setMass().

Referenced by G4INCL::ClusterDecay::decay().

446 { setMass(getRealMass()); }

G4INCL::Particle::setMass

void setMass(G4double mass)

Definition: G4INCLParticle.hh:558

G4INCL::Particle::getRealMass

G4double getRealMass() const

Get the real particle mass.

Definition: G4INCLParticle.hh:417

void G4INCL::Particle::setTableMass ( )

inline

Set the mass of the Particle to its table mass.

Definition at line 449 of file G4INCLParticle.hh.

References getTableMass(), and setMass().

Referenced by G4INCL::Nucleus::decayOutgoingDeltas(), G4INCL::Nucleus::emitInsidePions(), G4INCL::ProjectileRemnant::ProjectileRemnant(), and G4INCL::ProjectileRemnant::reset().

449 { setMass(getTableMass()); }

G4INCL::Particle::setMass

void setMass(G4double mass)

Definition: G4INCLParticle.hh:558

G4INCL::Particle::getTableMass

virtual G4double getTableMass() const

Get the tabulated particle mass.

Definition: G4INCLParticle.hh:388

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

inline

Definition at line 196 of file G4INCLParticle.hh.

References G4INCL::Composite, G4INCL::DeltaMinus, G4INCL::DeltaPlus, G4INCL::DeltaPlusPlus, G4INCL::DeltaZero, INCL_ERROR, isResonance(), G4INCL::Neutron, G4INCL::PiMinus, G4INCL::PiPlus, G4INCL::PiZero, G4INCL::Proton, setINCLMass(), theA, theType, theZ, and G4INCL::UnknownParticle.

Referenced by G4INCL::Cluster::Cluster(), G4INCL::ClusterDecay::decay(), G4INCL::DeltaProductionChannel::getFinalState(), G4INCL::PionNucleonChannel::getFinalState(), G4INCL::ElasticChannel::getFinalState(), G4INCL::DeltaDecayChannel::getFinalState(), G4INCL::RecombinationChannel::getFinalState(), and Particle().

                                  {
       theType = t;
       switch(theType)
       {
         case DeltaPlusPlus:
           theA = 1;
           theZ = 2;
           break;
         case Proton:
         case DeltaPlus:
           theA = 1;
           theZ = 1;
           break;
         case Neutron:
         case DeltaZero:
           theA = 1;
           theZ = 0;
           break;
         case DeltaMinus:
           theA = 1;
           theZ = -1;
           break;
         case PiPlus:
           theA = 0;
           theZ = 1;
           break;
         case PiZero:
           theA = 0;
           theZ = 0;
           break;
         case PiMinus:
           theA = 0;
           theZ = -1;
           break;
         case Composite:
          // INCL_ERROR("Trying to set particle type to Composite! Construct a Cluster object instead" << std::endl);
           theA = 0;
           theZ = 0;
           break;
         case UnknownParticle:
           theA = 0;
           theZ = 0;
           INCL_ERROR("Trying to set particle type to Unknown!" << std::endl);
           break;
       }
 
       if( !isResonance() && t!=Composite )
         setINCLMass();
     }

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

inline

Set the uncorrelated momentum.

Definition at line 772 of file G4INCLParticle.hh.

References uncorrelatedMomentum.

772 { uncorrelatedMomentum = p; }

p

const char * p

Definition: xmltok.h:285

G4INCL::Particle::uncorrelatedMomentum

G4double uncorrelatedMomentum

Definition: G4INCLParticle.hh:806

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

inlineprotected

Helper method for the assignment operator.

Definition at line 139 of file G4INCLParticle.hh.

References nCollisions, nDecays, rpCorrelated, CLHEP::swap(), theA, theEnergy, theFrozenEnergy, theFrozenMomentum, theMomentum, theParticipantType, thePosition, thePotentialEnergy, thePropagationEnergy, thePropagationMomentum, theType, theZ, and uncorrelatedMomentum.

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

                              {
       std::swap(theZ, rhs.theZ);
       std::swap(theA, rhs.theA);
       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);
     }

void G4INCL::Particle::thawPropagation ( )

inline

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 704 of file G4INCLParticle.hh.

References theEnergy, theMomentum, thePropagationEnergy, and thePropagationMomentum.

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

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

Field Documentation

long G4INCL::Particle::ID

protected

Definition at line 803 of file G4INCLParticle.hh.

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

G4int G4INCL::Particle::nCollisions

protected

Definition at line 800 of file G4INCLParticle.hh.

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

G4int G4INCL::Particle::nDecays

protected

Definition at line 801 of file G4INCLParticle.hh.

Referenced by getNumberOfDecays(), incrementNumberOfDecays(), setNumberOfDecays(), and swap().

G4bool G4INCL::Particle::rpCorrelated

protected

Definition at line 805 of file G4INCLParticle.hh.

Referenced by getReflectionMomentum(), rpCorrelate(), rpDecorrelate(), and swap().

G4int G4INCL::Particle::theA

protected

Definition at line 790 of file G4INCLParticle.hh.

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

G4double G4INCL::Particle::theEnergy

protected

G4double G4INCL::Particle::theFrozenEnergy

protected

Definition at line 795 of file G4INCLParticle.hh.

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

G4INCL::ThreeVector G4INCL::Particle::theFrozenMomentum

protected

Definition at line 798 of file G4INCLParticle.hh.

Referenced by freezePropagation(), getFrozenMomentum(), Particle(), rotate(), setFrozenMomentum(), and swap().

G4INCL::ThreeVector G4INCL::Particle::theMomentum

protected

ParticipantType G4INCL::Particle::theParticipantType

protected

Definition at line 791 of file G4INCLParticle.hh.

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

G4INCL::ThreeVector G4INCL::Particle::thePosition

protected

G4double G4INCL::Particle::thePotentialEnergy

protected

Definition at line 802 of file G4INCLParticle.hh.

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

G4double* G4INCL::Particle::thePropagationEnergy

protected

Definition at line 794 of file G4INCLParticle.hh.

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

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

protected

Definition at line 797 of file G4INCLParticle.hh.

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

G4INCL::ParticleType G4INCL::Particle::theType

protected

Definition at line 792 of file G4INCLParticle.hh.

Referenced by dump(), getINCLMass(), getRealMass(), getSpecies(), getTableMass(), getType(), isCluster(), isDelta(), isNucleon(), isPion(), G4INCL::Cluster::print(), print(), setType(), and swap().

G4int G4INCL::Particle::theZ

protected

Definition at line 790 of file G4INCLParticle.hh.

G4double G4INCL::Particle::uncorrelatedMomentum

protected

Definition at line 806 of file G4INCLParticle.hh.

Referenced by getReflectionMomentum(), setUncorrelatedMomentum(), and swap().

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

Protected Attributes
G4int	theZ

G4int	theA

ParticipantType	theParticipantType

G4INCL::ParticleType	theType

G4double	theEnergy

G4double *	thePropagationEnergy

G4double	theFrozenEnergy

G4INCL::ThreeVector	theMomentum

G4INCL::ThreeVector *	thePropagationMomentum

G4INCL::ThreeVector	theFrozenMomentum

G4INCL::ThreeVector	thePosition

G4int	nCollisions

G4int	nDecays

G4double	thePotentialEnergy

long	ID

G4bool	rpCorrelated

G4double	uncorrelatedMomentum

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Field Documentation